CpG island clusters and pro-epigenetic selection for CpGs in protein-coding exons of HOX and other transcription factors

Immune gene expression and epigenetic potential affect the consumption of risky food by female house sparrows

When organisms move into new areas, they are likely to encounter novel food resources. Even if they are nutritious, these foods can also be risky, as they might be contaminated by parasites. The behavioural immune system of animals could help them avoid the negative effects of contaminated resources, but our understanding of behavioural immunity is limited, particularly whether and how behavioural immunity interacts with physiological immunity. Here, we asked about the potential for interplay between these two traits, specifically how the propensity of an individual house sparrow (Passer domesticus) to take foraging risks was related to its ability to regulate a key facet of its immune response to bacterial pathogens. Previously, we found that sparrows at expanding geographic range edges were more exploratory and less risk-averse to novel foods; in those same populations, birds tended to over-express Toll-like receptor 4 (TLR4), a pattern-recognition receptor that distinguishes cell-wall components of Gram-negative bacteria, making it the major sensor of potentially lethal gut microbial infections including salmonellosis. When we investigated how birds would respond to a typical diet (i.e., mixed seeds) spiked with domesticated chicken faeces, birds that expressed more TLR4 or had higher epigenetic potential for TLR4 (more CpG dinucleotides in the putative gene promoter) ate more food, spiked or not. Females expressing abundant TLR4 were also willing to take more foraging risks and ate more spiked food. In males, TLR4 expression was not associated with risk-taking. Altogether, our results indicate that behaviour and immunity covary among individual house sparrows, particularly in females where those birds that maintain more immune surveillance also are more disposed to take foraging risks.

Interconnected Codons: Unravelling the Epigenetic Significance of Flanking Sequences in CpG Dyads

Hypothesizing that CpG codon dyads, formed by consecutive codons containing a cytosine-guanine pair (NNC-GNN), may play a crucial role in gene function, we conducted an extensive analysis to investigate their distribution and conservation within mammalian genes. Our findings reveal that genes characterized by a high density of CpG codon dyads are notably associated with homeobox domains and RNA polymerase II transcription factors. Conversely, genes exhibiting low CpG codon dyad density have links to DNA damage repair and mitosis. Importantly, our study identifies a re-markable increase in expressed genes that harbor CpG during embryonic development, suggesting their potential involvement in gene regulation at these developmental stages. These results under-score the functional significance of CpG codon dyads in DNA methylation and gene expression, fur-ther demonstrating the coevolution of consecutive codons and their contribution to codon usage bias.

Transgenerational Epigenetic DNA Methylation Editing and Human Disease

During gestation, maternal (F0), embryonic (F1), and migrating primordial germ cell (F2) genomes can be simultaneously exposed to environmental influences. Accumulating evidence suggests that operating epi- or above the genetic DNA sequence, covalent DNA methylation (DNAme) can be recorded onto DNA in response to environmental insults, some sites which escape normal germline erasure. These appear to intrinsically regulate future disease propensity, even transgenerationally. Thus, an organism's genome can undergo epigenetic adjustment based on environmental influences experienced by prior generations. During the earliest stages of mammalian development, the three-dimensional presentation of the genome is dramatically changed, and DNAme is removed genome wide. Why, then, do some pathological DNAme patterns appear to be heritable? Are these correctable? In the following sections, I review concepts of transgenerational epigenetics and recent work towards programming transgenerational DNAme. A framework for editing heritable DNAme and challenges are discussed, and ethics in human research is introduced.

Ratiometric Electrochemical Biosensing of Methyltransferase Activity

In this work, a novel ratiometric electrochemical readout platform was proposed and developed for the fast and flexible analysis of M.SssI methyltransferase (MTase) activity. In this platform, two hairpin DNAs (H1 and H2) were designed. H1 contains the palindromic sequence of 5′-CCGG-3′ in its stem which could be methylated and hybridize with H2 labeled by methylene blue (MB) as one of the signal reporters on a gold electrode (GE) in the presence of M.SssI MTase. Additionally, a specific immunoreaction was introduced by conjugating an anti-5-methylcytosine antibody, a DNA CpG methylation recognition unit, with 1,3-ferrocenedicarboxylic acid (Fc) as the second signal reporter. The results showed that when the Fc tag approaches, the MB tag was far from the gold electrode surface, resulting in a decrease in the oxidation peak current of MB (IMB) and an increase in the oxidation peak current of Fc (IFc). The ratiometric electrochemical method above shows the linear range of detection was 0 U/mL 40 U/mL with a detection limit of 0.083 U/mL (the mean signal of blank measures þ3s).

The Mutagenic Consequences of DNA Methylation within and across Generations

DNA methylation is an epigenetic modification with wide-ranging consequences across the life of an organism. This modification can be stable, persisting through development despite changing environmental conditions. However, in other contexts, DNA methylation can also be flexible, underlying organismal phenotypic plasticity. One underappreciated aspect of DNA methylation is that it is a potent mutagen; methylated cytosines mutate at a much faster rate than other genetic motifs. This mutagenic property of DNA methylation has been largely ignored in eco-evolutionary literature, despite its prevalence. Here, we explore how DNA methylation induced by environmental and other factors could promote mutation and lead to evolutionary change at a more rapid rate and in a more directed manner than through stochastic genetic mutations alone. We argue for future research on the evolutionary implications of DNA methylation driven mutations both within the lifetime of organisms, as well as across timescales.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates

Background

Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken).

Results

Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates.

Conclusions

Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01270-x.

Disruption of NIPBL/Scc2 in Cornelia de Lange Syndrome provokes cohesin genome-wide redistribution with an impact in the transcriptome

Cornelia de Lange syndrome (CdLS) is a rare disease affecting multiple organs and systems during development. Mutations in the cohesin loader, NIPBL/Scc2, were first described and are the most frequent in clinically diagnosed CdLS patients. The molecular mechanisms driving CdLS phenotypes are not understood. In addition to its canonical role in sister chromatid cohesion, cohesin is implicated in the spatial organization of the genome. Here, we investigate the transcriptome of CdLS patient-derived primary fibroblasts and observe the downregulation of genes involved in development and system skeletal organization, providing a link to the developmental alterations and limb abnormalities characteristic of CdLS patients. Genome-wide distribution studies demonstrate a global reduction of NIPBL at the NIPBL-associated high GC content regions in CdLS-derived cells. In addition, cohesin accumulates at NIPBL-occupied sites at CpG islands potentially due to reduced cohesin translocation along chromosomes, and fewer cohesin peaks colocalize with CTCF.

Epigenetic potential affects immune gene expression in house sparrows

Epigenetic mechanisms may play a central role in mediating phenotypic plasticity, especially during range expansions, when populations face a suite of novel environmental conditions. Individuals may differ in their epigenetic potential (EP; their capacity for epigenetic modifications of gene expression), which may affect their ability to colonize new areas. One form of EP, the number of CpG sites, is higher in introduced house sparrows (Passer domesticus) than in native birds in the promoter region of a microbial surveillance gene, Toll-like Receptor 4 (TLR4), which may allow invading birds to fine-tune their immune responses to unfamiliar parasites. Here, we compared TLR4 gene expression from whole blood, liver and spleen in house sparrows with different EP, first challenging some birds with lipopolysaccharide (LPS), to increase gene expression by simulating a natural infection. We expected that high EP would predict high inducibility and reversibility of TLR4 expression in the blood of birds treated with LPS, but we did not make directional predictions regarding organs, as we could not repeatedly sample these tissues. We found that EP was predictive of TLR4 expression in all tissues. Birds with high EP expressed more TLR4 in the blood than individuals with low EP, regardless of treatment with LPS. Only females with high EP exhibited reversibility in gene expression. Further, the effect of EP varied between sexes and among tissues. Together, these data support EP as one regulator of TLR4 expression.

Understanding the interplay between CpG island-associated gene promoters and H3K4 methylation

The precise regulation of gene transcription is required to establish and maintain cell type-specific gene expression programs during multicellular development. In addition to transcription factors, chromatin, and its chemical modification, play a central role in regulating gene expression. In vertebrates, DNA is pervasively methylated at CG dinucleotides, a modification that is repressive to transcription. However, approximately 70% of vertebrate gene promoters are associated with DNA elements called CpG islands (CGIs) that are refractory to DNA methylation. CGIs integrate the activity of a range of chromatin-regulating factors that can post-translationally modify histones and modulate gene expression. This is exemplified by the trimethylation of histone H3 at lysine 4 (H3K4me3), which is enriched at CGI-associated gene promoters and correlates with transcriptional activity. Through studying H3K4me3 at CGIs it has become clear that CGIs shape the distribution of H3K4me3 and, in turn, H3K4me3 influences the chromatin landscape at CGIs. Here we will discuss our understanding of the emerging relationship between CGIs, H3K4me3, and gene expression.

Epigenetic Potential in Native and Introduced Populations of House Sparrows (Passer domesticus)

Epigenetic potential, defined as the capacity for epigenetically-mediated phenotypic plasticity, may play an important role during range expansions. During range expansions, populations may encounter relatively novel challenges while experiencing lower genetic diversity. Phenotypic plasticity via epigenetic potential might be selectively advantageous at the time of initial introduction or during spread into new areas, enabling introduced organisms to cope rapidly with novel challenges. Here, we asked whether one form of epigenetic potential (i.e., the abundance of CpG sites) in three microbial surveillance genes: Toll-like receptors (TLRs) 1B (TLR1B), 2A (TLR2A), and 4 (TLR4) varied between native and introduced house sparrows (Passer domesticus). Using an opportunistic approach based on samples collected from sparrow populations around the world, we found that introduced birds had more CpG sites in TLR2A and TLR4, but not TLR1B, than native ones. Introduced birds also lost more CpG sites in TLR1B, gained more CpG sites in TLR2A, and lost fewer CpG sites in TLR4 compared to native birds. These results were not driven by differences in genetic diversity or population genetic structure, and many CpG sites fell within predicted transcription factor binding sites (TFBS), with losses and gains of CpG sites altering predicted TFBS. Although we lacked statistical power to conduct the most rigorous possible analyses, these results suggest that epigenetic potential may play a role in house sparrow range expansions, but additional work will be critical to elucidating how epigenetic potential affects gene expression and hence phenotypic plasticity at the individual, population, and species levels.

Sensitive fluorescent detection of methyltransferase based on thermosensitive poly(N-isopropylacrylamide)

DNA methyltransferase (MTase) has a crucial role in many biological processes, its abnormal expression level has been regarded as a predictive cancer biomarker. Herein, a sensitive fluorescence method based on thermosensitive poly (N-isopr-opylacrylamide) was developed to assay of M.SssI activity. When the M.SssI was introduced, dsDNA was methylated at palindromic sequence 5'-C^mCGG-3' and became resistant to cleavage by the endonuclease HpaII. Therefore, a biotin modified ssDNA and a FAM modified ssDNA were designed including the recognized sites for both methyltransferase M.SssI and endonuclease HpaII. By SA-biotin intereaction, the DNA was conjugated to thermosensitive poly (N-isopropylacrylamide) modified by SA, the methylated substrate fluorescence was increased with the concentration of M.SssI increasing. The proposed method has a low detection limit of 0.18 U/mL. This simple method can be a useful tool to apply in diagnosis and biomedical research, which was successfully investigated in the serum sample.

Insights Into de novo Mutation Variation in Lithuanian Exome

In the last decade, one of the biggest challenges in genomics research has been to distinguish definitive pathogenic variants from all likely pathogenic variants identified by next-generation sequencing. This task is particularly complex because of our lack of knowledge regarding overall genome variation and pathogenicity of the variants. Therefore, obtaining sufficient information about genome variants in the general population is necessary as such data could be used for the interpretation of de novo mutations (DNMs) in the context of patient's phenotype in cases of sporadic genetic disease. In this study, data from whole-exome sequencing of the general population in Lithuania were directly examined. In total, 84 (VarScan) and 95 (VarSeqTM) DNMs were identified and validated using different algorithms. Thirty-nine of these mutations were considered likely to be pathogenic based on gene function, evolutionary conservation, and mutation impact. The mutation rate estimated per position pair per generation was 2.74 × 10-8 [95% CI: 2.24 × 10-8-3.35 × 10-8] (VarScan) and 2.4 × 10-8 [95% CI: 1.96 × 10-8-2.99 × 10-8] (VarSeqTM), with 1.77 × 10-8 [95% CI: 6.03 × 10-9-5.2 × 10-8] de novo indels per position per generation. The rate of germline DNMs in the Lithuanian population and the effects of the genomic and epigenetic context on DNM formation were calculated for the first time in this study, providing a basis for further analysis of DNMs in individuals with genetic diseases. Considering these findings, additional studies in patient groups with genetic diseases with unclear etiology may facilitate our ability to distinguish certain pathogenic or adaptive DNMs from tolerated background DNMs and to reliably identify disease-causing DNMs by their properties through direct observation.

Caenorhabditis elegans as an emerging model system in environmental epigenetics

The roundworm Caenorhabitis elegans has been an established model organism for the study of genetics and developmental biology, including studies of transcriptional regulation, since the 1970s. This model organism has continued to be used as a classical model system as the field of transcriptional regulation has expanded to include scientific advances in epigenetics and chromatin biology. In the last several decades, C. elegans has emerged as a powerful model for environmental toxicology, particularly for the study of chemical genotoxicity. Here, we outline the utility and applicability of C. elegans as a powerful model organism for mechanistic studies of environmental influences on the epigenome. Our goal in this article is to inform the field of environmental epigenetics of the strengths and limitations of the well-established C. elegans model organism as an emerging model for medium-throughput, in vivo exploration of the role of exogenous chemical stimuli in transcriptional regulation, developmental epigenetic reprogramming, and epigenetic memory and inheritance. As the field of environmental epigenetics matures, and research begins to map mechanisms underlying observed associations, new toolkits and model systems, particularly manipulable, scalable in vivo systems that accurately model human transcriptional regulatory circuits, will provide an essential experimental bridge between in vitro biochemical experiments and mammalian model systems. Environ. Mol. Mutagen. 59:560-575, 2018. © 2018 Wiley Periodicals, Inc.

Genome-wide analyses reveal a role of Polycomb in promoting hypomethylation of DNA methylation valleys

BACKGROUND

Previous studies showed that the majority of developmental genes are devoid of DNA methylation at promoters even when they are repressed. Such hypomethylated regions at developmental genes are unusually large and extend well beyond proximal promoters, forming DNA methylation valleys (DMVs) or DNA methylation canyons. However, it remains elusive how most developmental genes can evade DNA methylation regardless of their transcriptional states.

RESULTS

We show that DMVs are hypomethylated in development and are highly conserved across vertebrates. Importantly, DMVs are hotspots of regulatory regions for key developmental genes and show low levels of deamination mutation rates. By analyzing a panel of DNA methylomes from mouse tissues, we identify a subset of DMVs that are dynamically methylated. These DMVs are strongly enriched for Polycomb-deposited H3K27me3 when the associated genes are silenced, and surprisingly show elevated DNA methylation upon gene activation. 4C-seq analyses indicates that Polycomb-bound DMVs form insulated and self-interacting chromatin domains. Further investigations show that DNA hypomethylation is better correlated with the binding of Polycomb than with H3K27me3. In support of a role of Polycomb in DMV hypomethylation, we observe aberrant methylation in DMVs in mouse embryonic stem cells deficient in the EED protein. Finally, we show that Polycomb regulates hypomethylation of DMVs likely through ten-eleven translocation (TET) proteins.

CONCLUSIONS

We show that Polycomb promotes the hypomethylation of DMVs near key developmental genes. These data reveal a delicate interplay between histone modifiers and DNA methylation, which contributes to their division at distinct gene targets, allowing lineage-specifying genes to largely maintain DNA methylation-free at regulatory elements.

Single quantum dot-based nanosensor for sensitive detection of 5-methylcytosine at both CpG and non-CpG sites

DNA methylation is an important epigenetic modification in human genomes. Herein, we develop a single quantum dot (QD)-based nanosensor for sensitive detection of DNA methylation at both CpG and non-CpG sites using tricyclic ligation chain reaction (LCR)-mediated QD-based fluorescence resonance energy transfer (FRET). We design two sets of DNA probes (X and Y, X' and Y') for methylated DNA assay. In the presence of thermostable DNA ligase, probes X and Y may adjacently hybridize with the methylated DNA to obtain the ligated XY products which may function as the templates for probes X' and Y' to generate the X'Y' products. The resultant X'Y' products may in turn act as the templates to ligate probes X and Y for the generation of XY products, consequently inducing tricyclic LCR amplification under thermal denaturation conditions to generate a large number of XY products. The subsequent hybridization of XY products with the capture and reporter probes results in the formation of sandwich hybrids which may assemble on the 605QD surface to obtain 605QD-oligonucleotide-Cy5 nanostructures, inducing efficient FRET from the 605QD to Cy5 and the emission of Cy5. This nanosensor can detect DNA methylation at single 5-methylcytosine (5-mC) resolution with a detection limit of as low as 1.0 aM and a large dynamic range of 7 orders of magnitude. Moreover, this nanosensor can distinguish as low as a 0.01% methylation level, and it can detect DNA methylation in human lung cancer cells as well, holding great potential for accurate epigenetic evaluation and early cancer diagnosis.

Electrochemical signal-amplified detection of 5-methylcytosine and 5-hydroxymethylcytosine in DNA using glucose modification coupled with restriction endonucleases

A sensitive electrochemical assay based on glucose modification coupled with restriction endonucleases was fabricated for the detection of 5-methylcytosine and 5-hydroxymethylcytosine in DNA.

Genomic landscape of CpG rich elements in human

Background

The studies on CpG islands (CGI) and Alu elements functions, evolution, and distribution in the genome started since the discovery in nineteen eighties (1981, 1986, correspondingly). Their highly skewed genome wide distribution implies the non-random retrotransposition pattern. Besides CGIs in gene promoters, CGIs clusters were observed in the homeobox gene regions and in the macrosatellites, but the whole picture of their distribution specifics was not grasped. Attempts to identify any causative features upon their (genome wide) distribution, such as the DNA context mediated preferred insertion sites of Alu repeats, have been made to ascribe their clusters location.

Methods

Recent emergence of high resolution 3D map of human genome allowed segregating the genome into the large scale chromatin domains of naturally observable nuclear subcompartments, or Topologically Associated Domains (TADs), designated by spatial chromatin distribution. We utilized the chromatin map to elucidate relations between large scale chromatin state and CpG rich elements landscape.

In the course of analysis it was confirmed that genes, Alu and CGI clusters maintain obvious, albeit different in strength, preference for open chromatin. For the first time it was clearly shown that the clusters density of the Alu and CGIs monotonically depend on the chromatin accessibility rate. In particular, the highest density of these elements is found in A1 euchromatin regions characterized by a high density of small length genes replicating in the early S-phase. It implies that these elements mediate (CGIs) or are a side element (Alus) of chromatin accessibility.

Results

We elucidated that both methylated and non-methylated CGIs display the affinity to chromatin accessibility. As a part of comparative genomics section, we elucidated that the dog’s genome non-canonical structure, outstanding in mammals for its high CGIs abundance compared to gene number, is explained by the presence of dense tandem CGI extended hotspots (500 kb on average) in subtelomeric and pericentromeric regions with highly skewed CG content, and not by CGIs global distribution pattern shift.

Conclusions

The study underlines the close association of CG-rich elements distribution with the newly introduced large scale chromatin state map, proposing a refined standpoint on interrelation of aforementioned genome elements and the chromatin state. To our expertise, the TAD-associated partition model employed in the study is likely the most substantial one regarding CpG rich clusters distribution among the whole genome chromatin/isochores maps available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0864-0) contains supplementary material, which is available to authorized users.

CTCF-mediated Chromatin Loop for the Posterior Hoxc Gene Expression in MEF Cells

Modulation of chromatin structure has been proposed as a molecular mechanism underlying the spatiotemporal collinear expression of Hox genes during development. CCCTC-binding factor (CTCF)-mediated chromatin organization is now recognized as a crucial epigenetic mechanism for transcriptional regulation. Thus, we examined whether CTCF-mediated chromosomal conformation is involved in Hoxc gene expression by comparing wild-type mouse embryonic fibroblast (MEF) cells expressing anterior Hoxc genes with Akt1 null MEFs expressing anterior as well as posterior Hoxc genes. We found that CTCF binding between Hoxc11 and -c12 is important for CTCF-mediated chromosomal loop formation and concomitant posterior Hoxc gene expression. Hypomethylation at this site increased CTCF binding and recapitulated the chromosomal conformation and posterior Hoxc gene expression patterns observed in Akt1 null MEFs. From this work we found that CTCF at the C12|11 does not function as a barrier/boundary, instead let the posterior Hoxc genes switch their interaction from inactive centromeric to active telomeric genomic niche, and concomitant posterior Hoxc gene expression. Although it is not clear whether CTCF affects Hoxc gene expression solely through its looping activity, CTCF-mediated chromatin structural modulation could be an another tier of Hox gene regulation during development. © 2016 IUBMB Life, 68(6):436-444, 2016.

Epigenetics: A primer for clinicians

With recent advances in cellular biology, we now appreciate that modifications to DNA and histones can have a profound impact on transcription and function, even in the absence of changes to DNA sequence. These modifications, now commonly referred to as "epigenetic" alterations, have changed how we understand cell behavior, reprogramming and differentiation and have provided significant insight into the mechanisms underlying carcinogenesis. Epigenetic alterations, to this point, are largely identified by changes in DNA methylation and hydroxymethylation as well as methylation, acetylation, and phosphorylation of histone tails. These modifications enable significant flexibility in gene expression, rather than just turning genes "ON" or "OFF." Herein we describe the epigenetic landscape in the regulation of gene expression with a particular focus on interrogating DNA methylation in myeloid malignancy.

Mapping Human Pluripotent-to-Cardiomyocyte Differentiation: Methylomes, Transcriptomes, and Exon DNA Methylation "Memories"

The directed differentiation of human cardiomyocytes (CMs) from pluripotent cells provides an invaluable model for understanding mechanisms of cell fate determination and offers considerable promise in cardiac regenerative medicine. Here, we utilize a human embryonic stem cell suspension bank, produced according to a good manufacturing practice, to generate CMs using a fully defined and small molecule-based differentiation strategy. Primitive and cardiac mesoderm purification was used to remove non-committing and multi-lineage populations and this significantly aided the identification of key transcription factors, lncRNAs, and essential signaling pathways that define cardiomyogenesis. Global methylation profiles reflect CM development and we report on CM exon DNA methylation "memories" persisting beyond transcription repression and marking the expression history of numerous developmentally regulated genes, especially transcription factors.

Regulation of Hox orthologues in the oyster Crassostrea gigas evidences a functional role for promoter DNA methylation in an invertebrate

DNA methylation within promoter regions (PRDM) controls vertebrate early gene transcription and thereby development, but is neglected outside this group. However, epigenetic features in the oyster Crassostrea gigas suggest functional significance of PDRM in invertebrates. To investigate this, reporter constructs containing in vitro methylated oyster Hox gene promoters were transfected into oyster embryos. The influence of in vivo methylation was studied using bisulfite sequencing and DNA methyltransferase inhibition during development. Our results demonstrate that methylation controls the transcriptional activity of the promoters investigated, unraveling a functional role for PRDM in a lophotrochozoan, an important finding regarding the evolution of epigenetic regulation.

Altered histone mark deposition and DNA methylation at homeobox genes in human oral squamous cell carcinoma

We recently reported a role of polycomb repressive complex 2 (PRC2) and PRC2 trimethylation of histone 3 lysine 27 (H3K27me3) in the regulation of homeobox (HOX) (Marcinkiewicz and Gudas, 2013, Exp Cell Res) gene transcript levels in human oral keratinocytes (OKF6-TERT1R) and tongue squamous cell carcinoma (SCC) cells. Here, we assessed both the levels of various histone modifications at a subset of homeobox genes and genome wide DNA methylation patterns in OKF6-TERT1R and SCC-9 cells by using ERRBS (enhanced reduced representation bisulfite sequencing). We detected the H3K9me3 mark at HOXB7, HOXC10, HOXC13, and HOXD8 at levels higher in OKF6-TERT1R than in SCC-9 cells; at IRX1 and SIX2 the H3K9me3 levels were conversely higher in SCC-9 than in OKF6-TERT1R. The H3K79me3 mark was detectable only at IRX1 in OKF6-TERT1R and at IRX4 in SCC-9 cells. The levels of H3K4me3 and H3K36me3 marks correlate with the transcript levels of the assessed homeobox genes in both OKF6-TERT1R and SCC-9. We detected generally lower CpG methylation levels on DNA in SCC-9 cells at annotated genomic regions which were differentially methylated between OKF6-TERT1R and SCC-9 cells; however, some genomic regions, including the HOX gene clusters, showed DNA methylation at higher levels in SCC-9 than OKF6-TERT1R. Thus, both altered histone modification patterns and changes in DNA methylation are associated with dysregulation of homeobox gene expression in human oral cavity SCC cells, and this dysregulation potentially plays a role in the neoplastic phenotype of oral keratinocytes.

Anti-viral gene induction is absent upon secondary challenge with double-stranded RNA in the Pacific oyster, Crassostrea gigas

Oyster farming is one of the most important aquaculture industries in the world. However, its productivity is increasingly limited by viral disease and we do not yet have management practices, such as protective vaccination, that can control these disease outbreaks. Hence, in the current study we investigated the expression of known anti-viral genes in oysters (Crassostrea gigas) in response to primary and secondary encounter with a virus associated molecular pattern (dsRNA), and tested whether a common form of epigenetic gene regulation (DNA methylation) was associated with the expression of these anti-viral genes. Injection of dsRNA into the adductor muscle resulted in the rapid and transient expression of virus recognition receptors (TLR & MDA5), whereas several anti-viral signalling (IRF & SOC-1) and effector (PKR & viperin) genes were still up-regulated at one week post primary challenge (p < 0.05). This primary encounter with dsRNA appeared to deplete the immune system because anti-viral gene induction was absent in the gills when oysters were given a second injection of dsRNA at 1-week post-primary injection. The expression of DNA methylation genes (DNMT1, DNMT3b, TDG, TET2) and DNA methylation profiles up-stream of specific anti-viral genes (STING, SOC-1 & Viperin) did not change in response to dsRNA injection (p > 0.05). These results collectively suggest that C. gigas does not have an enhanced anti-viral gene response (immune-priming) to secondary dsRNA challenge and that the sustained up-regulation of anti-viral signalling and effector genes following primary challenge is unlikely to be associated with upstream DNA methylation levels.

Akt1 mediates the posterior Hoxc gene expression through epigenetic modifications in mouse embryonic fibroblasts

The evolutionarily conserved Hox genes are organized in clusters and expressed colinearly to specify body patterning during embryonic development. Previously, Akt1 has been identified as a putative Hox gene regulator through in silico analysis. Substantial upregulation of consecutive 5' Hoxc genes has been observed when Akt1 is absent in mouse embryonic fibroblast (MEF) cells. In this study, we provide evidence that Akt1 regulates the 5' Hoxc gene expression by epigenetic modifications. Enrichment of histone H3K9 acetylation and a low level of the H3K27me3 mark were detected at the posterior 5' Hoxc loci when Akt1 is absent. A histone deacetylase (HDAC) inhibitor de-repressed 5' Hoxc gene expression when Akt1 is present, and a DNA demethylating reagent synergistically upregulated HDAC-induced 5' Hoxc gene expression. A knockdown study revealed that Hdac6 is mediated in the Hoxc12 repression through direct binding to the transcription start site (TSS) in the presence of Akt1. Co-immunoprecipitation analysis revealed that endogenous Akt1 directly interacted with Hdac6. Furthermore, exogenous Akt1 was enriched at the promoter region of the posterior Hoxc genes such as Hoxc11 and Hoxc12, not the Akt1-independent Hoxc5 and Hoxd10 loci. The regulation of the H3K27me3 mark by Ezh2 and Kdm6b at the 5' Hoxc gene promoter turned out to be Akt1 dependent. Taken together, these results suggest that Akt1 mediates the posterior 5' Hoxc gene expression through epigenetic modification such as histone methylation and acetylation, and partly through a direct binding to the promoter region of the 5' Hoxc genes and/or Hdac6 in mouse embryonic fibroblast cells.

Enhanced evolution by stochastically variable modification of epigenetic marks in the early embryo

Evolution by gene duplication is generally accepted as one of the crucial driving forces for the gain of new complexity and functions, but the formation of pseudogenes remains a problem for this mechanism. Here we expand on earlier ideas that epigenetic modifications can drive neo- and subfunctionalization in evolution by gene duplication. We explore the effects of stochastic epigenetic modifications on the evolution (and thus development) of complex organisms in a constant environment. Modeling is done both using a modified genetic drift analytical treatment and computer simulations, which were found to agree. A transposon silencing model is also explored. Some key assumptions made include (i) stochastic, incomplete removal (or addition) of repressive epigenetic marks takes place during a window(s) of opportunity in the zygote and early embryo; (ii) there is no statistical variation of the marks after the window closes; and (iii) the genes affected are sensitive to dosage. Our genetic drift treatment takes into account that after gene duplication the prevailing case upon which selection operates is a duplicate/singlet heterozygote; to the best of our knowledge, this has not been considered in previous treatments. We conclude from our modeling that stochastic epigenetic modifications, with rates consistent with experimental observation, can both increase the rate of gene fixation and decrease pseudogenization, thus dramatically improving the efficacy of evolution by gene duplication. We also find that a transposon silencing model is advantageous for fixation of recessive genes in diploid organisms, especially with large effective population sizes.

DNA methylation is crucial for the early development in the Oyster C. gigas

In vertebrates, epigenetic modifications influence gene transcription, and an appropriate DNA methylation is critical in development. Indeed, a precise temporal and spatial pattern of early gene expression is mandatory for a normal embryogenesis. However, such a regulation and its underlying mechanisms remain poorly understood in more distant organisms such as Lophotrochozoa. Thus, despite DNA in the oyster genome being methylated, the role of DNA methylation in development is unknown. To clarify this point, oyster genomic DNA was examined during early embryogenesis and found differentially methylated. Reverse transcriptase quantitative polymerase chain reaction indicated stage-specific levels of transcripts encoding DNA-methyltransferase (DNMT) and methyl-binding domain proteins. In addition, as highlighted by electronic microscopy and immunohistochemistry, the DNMT inhibitor 5-aza-cytidine induced alterations in the quantity and the localisation of methylated DNA and severe dose-dependent development alterations and was lethal after zygotic genome reinitiation. Furthermore, methyl-DNA-immunoprecipitation-quantitative polymerase chain reaction revealed that the transcription level of most of the homeobox gene orthologues examined, but not of the other early genes investigated, was inversely correlated with their specific DNA methylation. Altogether, our results demonstrate that DNA methylation influences gene expression in Crassostrea gigas and is critical for oyster development, possibly by specifically controlling the transcription level of homeobox orthologues. These findings provide evidence for the importance of epigenetic regulation of development in Lophotrochozoans and bring new insights into the early life of C. gigas, one of the most important aquaculture resources worldwide.

A distinct group of CpG islands shows differential DNA methylation between replicas of the same cell line in vitro

BACKGROUND

CpG dinucleotide-rich genomic DNA regions, known as CpG islands (CGIs), can be methylated at their cytosine residues as an epigenetic mark that is stably inherited during cell mitosis. Differentially methylated regions (DMRs) are genomic regions showing different degrees of DNA methylation in multiple samples. In this study, we focused our attention on CGIs showing different DNA methylation between two culture replicas of the same cell line.

RESULTS

We used methylation data of 35 cell lines from the Encyclopedia of DNA Elements (ENCODE) consortium to identify CpG islands that were differentially methylated between replicas of the same cell line and denoted them Inter Replicas Differentially Methylated CpG islands (IRDM-CGIs). We identified a group of IRDM-CGIs that was consistently shared by different cell lines, and denoted it common IRDM-CGIs. X chromosome CGIs were overrepresented among common IRDM-CGIs. Autosomal IRDM-CGIs were preferentially located in gene bodies and intergenic regions had a lower G + C content, a smaller mean length, and a reduced CpG percentage. Functional analysis of the genes associated with autosomal IRDM-CGIs showed that many of them are involved in DNA binding and development.

CONCLUSIONS

Our results show that several specific functional and structural features characterize common IRDM-CGIs. They may represent a specific subset of CGIs that are more prone to being differentially methylated for their intrinsic characteristics.

Identification of a functionally distinct truncated BDNF mRNA splice variant and protein in Trachemys scripta elegans

Brain-derived neurotrophic factor (BDNF) has a diverse functional role and complex pattern of gene expression. Alternative splicing of mRNA transcripts leads to further diversity of mRNAs and protein isoforms. Here, we describe the regulation of BDNF mRNA transcripts in an in vitro model of eyeblink classical conditioning and a unique transcript that forms a functionally distinct truncated BDNF protein isoform. Nine different mRNA transcripts from the BDNF gene of the pond turtle Trachemys scripta elegans (tBDNF) are selectively regulated during classical conditioning: exon I mRNA transcripts show no change, exon II transcripts are downregulated, while exon III transcripts are upregulated. One unique transcript that codes from exon II, tBDNF2a, contains a 40 base pair deletion in the protein coding exon that generates a truncated tBDNF protein. The truncated transcript and protein are expressed in the naïve untrained state and are fully repressed during conditioning when full-length mature tBDNF is expressed, thereby having an alternate pattern of expression in conditioning. Truncated BDNF is not restricted to turtles as a truncated mRNA splice variant has been described for the human BDNF gene. Further studies are required to determine the ubiquity of truncated BDNF alternative splice variants across species and the mechanisms of regulation and function of this newly recognized BDNF protein.

Chromatin-linked determinants of zygotic genome activation

The merging of the maternal and paternal genomes into a single pronucleus after fertilization is accompanied by a remarkable reconfiguration of chromatin in the newly formed zygote. The first stages of embryonic chromatin remodeling take place in the absence of ongoing transcription, during a species-specific developmental time-frame. Once post-fertilization chromatin states are organized, zygotic genome activation (ZGA) is initiated, and embryonic transcripts gradually take control of development. We review here transitions in chromatin modifications associated with the onset of ZGA, and the role of transcription factors and DNA motifs in the regulation of ZGA. We propose a model of sequential chromatin remodeling events preceding ZGA, leading to the onset of embryonic transcription.

The human placenta methylome

Tissue-specific DNA methylation is found at promoters, enhancers, and CpG islands but also over larger genomic regions. In most human tissues, the vast majority of the genome is highly methylated (>70%). Recently, sequencing of bisulfite-treated DNA (MethylC-seq) has revealed large partially methylated domains (PMDs) in some human cell lines. PMDs cover up to 40% of the genome and are associated with gene repression and inactive chromatin marks. However, to date, only cultured cells and cancers have shown evidence for PMDs. Here, we performed MethylC-seq in full-term human placenta and demonstrate it is the first known normal tissue showing clear evidence of PMDs. We found that PMDs cover 37% of the placental genome, are stable throughout gestation and between individuals, and can be observed with lower sensitivity in Illumina 450K Infinium data. RNA-seq analysis confirmed that genes in PMDs are repressed in placenta. Using a hidden Markov model to map placental PMDs genome-wide and compare them to PMDs in other cell lines, we found that genes within placental PMDs have tissue-specific functions. For regulatory regions, methylation levels in promoter CpG islands are actually higher for genes within placental PMDs, despite the lower overall methylation of surrounding regions. Similar to PMDs, polycomb-regulated regions are hypomethylated but smaller and distinct from PMDs, with some being hypermethylated in placenta compared with other tissues. These results suggest that PMDs are a developmentally dynamic feature of the methylome that are relevant for understanding both normal development and cancer and may be of use as epigenetic biomarkers.

Non-canonical functions of the DNA methylome in gene regulation

Methylation of the cytosine base in DNA, DNA methylation, is an essential epigenetic mark in mammals that contributes to the regulation of transcription. Several advances have been made in this area in recent years, leading to a leap forward in our understanding of how this pathway contributes to gene regulation during embryonic development, and the functional consequences of its perturbation in human disease. Critical to these advances is a comprehension of the genomic distribution of modified cytosine bases in unprecedented detail, drawing attention to genomic regions beyond gene promoters. In addition, we have a more complete understanding of the multifactorial manner by which DNA methylation influences gene regulation at the molecular level, and which genes rely directly on the DNA methylome for their normal transcriptional regulation. It is becoming apparent that a major role of DNA modification is to act as a relatively stable, and mitotically heritable, template that contributes to the establishment and maintenance of chromatin states. In this regard, interplay is emerging between DNA methylation and the PcG (Polycomb group) proteins, which act as evolutionarily conserved mediators of cell identity. In the present paper we review these aspects of DNA methylation, and discuss how a multifunctional view of DNA modification as an integral part of chromatin organization is influencing our understanding of this epigenetic mark's contribution to transcriptional regulation.

Early de novo DNA methylation and prolonged demethylation in the muscle lineage

Myogenic cell cultures derived from muscle biopsies are excellent models for human cell differentiation. We report the first comprehensive analysis of myogenesis-specific DNA hyper- and hypo-methylation throughout the genome for human muscle progenitor cells (both myoblasts and myotubes) and skeletal muscle tissue vs. 30 non-muscle samples using reduced representation bisulfite sequencing. We also focused on four genes with extensive hyper- or hypo-methylation in the muscle lineage (PAX3, TBX1, MYH7B/MIR499 and OBSCN) to compare DNA methylation, DNaseI hypersensitivity, histone modification, and CTCF binding profiles. We found that myogenic hypermethylation was strongly associated with homeobox or T-box genes and muscle hypomethylation with contractile fiber genes. Nonetheless, there was no simple relationship between differential gene expression and myogenic differential methylation, rather only for subsets of these genes, such as contractile fiber genes. Skeletal muscle retained ~30% of the hypomethylated sites but only ~3% of hypermethylated sites seen in myogenic progenitor cells. By enzymatic assays, skeletal muscle was 2-fold enriched globally in genomic 5-hydroxymethylcytosine (5-hmC) vs. myoblasts or myotubes and was the only sample type enriched in 5-hmC at tested myogenic hypermethylated sites in PAX3/CCDC140 andTBX1. TET1 and TET2 RNAs, which are involved in generation of 5-hmC and DNA demethylation, were strongly upregulated in myoblasts and myotubes. Our findings implicate de novo methylation predominantly before the myoblast stage and demethylation before and after the myotube stage in control of transcription and co-transcriptional RNA processing. They also suggest that, in muscle, TET1 or TET2 are involved in active demethylation and in formation of stable 5-hmC residues.

Methylation alterations of WT1 and homeobox genes in inflamed muscle biopsy samples from patients with untreated juvenile dermatomyositis suggest self-renewal capacity

OBJECTIVE

To determine the effect of methylation alteration in inflamed muscles from children with juvenile dermatomyositis (DM) and other idiopathic inflammatory myopathies (IIMs).

METHODS

Magnetic resonance imaging-directed diagnostic muscle biopsies yielded samples from 20 children with juvenile DM, which were used for genome-wide DNA methylation profiling, as were muscle biopsy samples from 4 healthy controls. Bisulfite treatment followed by pyrosequencing confirmed methylation status in juvenile DM and other IIMs. Immunohistochemistry defined localization and expression levels of WT1.

RESULTS

Comparison of genome-wide DNA methylation profiling between juvenile DM muscle and normal control muscle revealed 27 genes with a significant methylation difference between the groups. These genes were enriched with transcription factors and/or cell cycle regulators and were unrelated to duration of untreated disease. Six homeobox genes were among them; ALX4, HOXC11, HOXD3, and HOXD4 were hypomethylated, while EMX2 and HOXB1 were hypermethylated. WT1 was significantly hypomethylated in juvenile DM (Δβ = -0.41, P < 0.001). Bisulfite pyrosequencing verification in samples from 56 patients with juvenile DM confirmed the methylation alterations of these genes. Similar methylation alterations were observed in juvenile polymyositis (n = 5) and other IIMs (n = 9). Concordant with the other findings, WT1 protein was increased in juvenile DM muscle, with average positive staining of 11.6%, but was undetectable in normal muscle (P < 0.001).

CONCLUSION

These results suggest that affected muscles of children with juvenile DM and IIMs have the capacity to be repaired, and that homeobox and WT1 genes are epigenetically marked to facilitate this repair process, potentially suggesting new avenues of therapeutic intervention.

Human-specific CpG "beacons" identify loci associated with human-specific traits and disease

Regulatory change has long been hypothesized to drive the delineation of the human phenotype from other closely related primates. Here we provide evidence that CpG dinucleotides play a special role in this process. CpGs enable epigenome variability via DNA methylation, and this epigenetic mark functions as a regulatory mechanism. Therefore, species-specific CpGs may influence species-specific regulation. We report non-polymorphic species-specific CpG dinucleotides (termed “CpG beacons”) as a distinct genomic feature associated with CpG island (CGI) evolution, human traits and disease. Using an inter-primate comparison, we identified 21 extreme CpG beacon clusters (≥ 20/kb peaks, empirical p < 1.0 × 10⁻³) in humans, which include associations with four monogenic developmental and neurological disease related genes (Benjamini-Hochberg corrected p = 6.03 × 10⁻³). We also demonstrate that beacon-mediated CpG density gain in CGIs correlates with reduced methylation in these species in orthologous CGIs over time, via human, chimpanzee and macaque MeDIP-seq. Therefore mapping into both the genomic and epigenomic space the identified CpG beacon clusters define points of intersection where a substantial two-way interaction between genetic sequence and epigenetic state has occurred. Taken together, our data support a model for CpG beacons to contribute to CGI evolution from genesis to tissue-specific to constitutively active CGIs.

Developmental features of DNA methylation during activation of the embryonic zebrafish genome

BACKGROUND

Zygotic genome activation (ZGA) occurs at the mid-blastula transition (MBT) in zebrafish and is a period of extensive chromatin remodeling. Genome-scale gametic demethylation and remethylation occurs after fertilization, during blastula stages, but how ZGA relates to promoter DNA methylation states is unknown. Using methylated DNA immunoprecipitation coupled to high-density microarray hybridization, we characterize genome-wide promoter DNA methylation dynamics before, during and after ZGA onset, in relation to changes in post-translational histone modifications and gene expression.

RESULTS

We show methylation of thousands of promoters before ZGA and additional methylation after ZGA, finding more dynamic methylation -1 to 0 kb upstream of the transcription start site than downstream. The MBT is marked by differential methylation of high and low CpG promoters, and we identify hypomethylated promoters that are mostly CG-rich and remain hypomethylated through the MBT. Hypomethylated regions constitute a platform for H3K4me3, whereas H3K9me3 preferentially associates with methylated regions. H3K27me3 associates with either methylation state depending on its coincidence with H3K4me3 or H3K9me3. Cohorts of genes differentially expressed through the MBT period display distinct promoter methylation patterns related to CG content rather than transcriptional fate. Lastly, although a significant proportion of genes methylated in sperm are unmethylated in embryos, over 90% of genes methylated in embryos are also methylated in sperm.

CONCLUSIONS

Our results suggest a pre-patterning of developmental gene expression potential by a combination of DNA hypomethylation and H3K4 trimethylation on CG-rich promoters, and are consistent with a transmission of DNA methylation states from gametes to early embryos.

Complex selection associated with Hox genes in a natural population of lizards

Hox genes are recognized for their explanatory power of bilateral development. However, relatively little is known about natural variation in, and the evolutionary dynamics of, Hox genes within wild populations. Utilizing a natural population of sand lizards (Lacerta agilis), we screened HoxA13 for genetic variation and an association with incidence of offspring malformations. We found significant effects of parental genetic similarity and offspring sex, and their interaction, on risk of hatching malformed as an offspring. We also found within population genetic variation in HoxA13, and identified a significant effect of a three-way interaction among Hox genotype, parental genetic similarity, and offspring sex on the risk of hatching malformation. Since malformed offspring in this population do not survive to maturity, this study reveals complex and ongoing selection associated with Hox genes in a wild reptile population. Importantly, this demonstrates the utility of natural populations in unveiling microevolutionary processes shaping variation in highly conserved genes.

Ribozymes: Flexible molecular devices at work

The discovery of ribozymes, RNAs with catalytic activity, revealed the extraordinary characteristic of this molecule, and corroborated the idea that RNA was the first informative polymer. The "RNA world" hypothesis asserts that the DNA/RNA/PROTEIN world arose from an earlier RNA world in which were present only RNA molecules able to perform both of the two functions performed separately by DNA and proteins in the present-day cells: the ability to transfer genetic information and to carry out catalytic activity. The catalytic properties of ribozymes are exclusively due to the capacity of RNA molecules to assume particular structures. Moreover, the structural versatility of RNA can allow to a single RNA sequence to fold in more than one structure, able to perform more than one function. In the first part of this work we will discuss the RNA plasticity, focusing on "bifunctional" ribozymes isolated by in vitro selection experiments, and on the consequences of this plasticity in the prospective of the emergence of new specific functions. The possibility that one sequence could have more than one structure/function, greatly increase the evolutionary potential of RNA, and the capacity of RNA to switch from a structure/function to another is probably one of the reasons of the evolutionary success also in modern-day cells. Naturally occurring ribozymes discovered in contemporary cells, demonstrate the crucial role that ribozymes still have in the modern protein world. In the second part of this paper we will discuss the capacity of natural ribozymes to modulate gene expression making use of their exclusive catalytic properties. Moreover, we will consider the possibility of their ancient origin.