An annotated list of bivalent chromatin regions in human ES cells: a new tool for cancer epigenetic research

ChIP-DIP: A multiplexed method for mapping hundreds of proteins to DNA uncovers diverse regulatory elements controlling gene expression

Gene expression is controlled by the dynamic localization of thousands of distinct regulatory proteins to precise regions of DNA. Understanding this cell-type specific process has been a goal of molecular biology for decades yet remains challenging because most current DNA-protein mapping methods study one protein at a time. To overcome this, we developed ChIP-DIP (ChIP Done In Parallel), a split-pool based method that enables simultaneous, genome-wide mapping of hundreds of diverse regulatory proteins in a single experiment. We demonstrate that ChIP-DIP generates highly accurate maps for all classes of DNA-associated proteins, including histone modifications, chromatin regulators, transcription factors, and RNA Polymerases. Using these data, we explore quantitative combinations of protein localization on genomic DNA to define distinct classes of regulatory elements and their functional activity. Our data demonstrate that ChIP-DIP enables the generation of 'consortium level', context-specific protein localization maps within any molecular biology lab.

Electrostatic encoding of genome organization principles within single native nucleosomes

The eukaryotic genome, first packed into nucleosomes of about 150 bp around the histone core, is organized into euchromatin and heterochromatin, corresponding to the A and B compartments, respectively. Here, we asked if individual nucleosomes in vivo know where to go. That is, do mono-nucleosomes by themselves contain A/B compartment information, associated with transcription activity, in their biophysical properties? We purified native mono-nucleosomes to high monodispersity and used physiological concentrations of biological polyamines to determine their condensability. The chromosomal regions known to partition into A compartments have low condensability and vice versa. In silico chromatin polymer simulations using condensability as the only input showed that biophysical information needed to form compartments is all contained in single native nucleosomes and no other factors are needed. Condensability is also strongly anticorrelated with gene expression, and especially so near the promoter region and in a cell type dependent manner. Therefore, individual nucleosomes in the promoter know whether the gene is on or off, and that information is contained in their biophysical properties. Comparison with genetic and epigenetic features suggest that nucleosome condensability is a very meaningful axis onto which to project the high dimensional cellular chromatin state. Analysis of condensability using various condensing agents including those that are protein-based suggests that genome organization principle encoded into individual nucleosomes is electrostatic in nature. Polyamine depletion in mouse T cells, by either knocking out ornithine decarboxylase (ODC) or inhibiting ODC, results in hyperpolarized condensability, suggesting that when cells cannot rely on polyamines to translate biophysical properties of nucleosomes to control gene expression and 3D genome organization, they accentuate condensability contrast, which may explain dysfunction known to occur with polyamine deficiency.

The pediatric leukemia oncoprotein NUP98-KDM5A induces genomic instability that may facilitate malignant transformation

Pediatric Acute Myeloid Leukemia (AML) is a rare and heterogeneous disease characterized by a high prevalence of gene fusions as driver mutations. Despite the improvement of survival in the last years, about 50% of patients still experience a relapse. It is not possible to improve prognosis only with further intensification of chemotherapy, as come with a severe cost to the health of patients, often resulting in treatment-related death or long-term sequels. To design more effective and less toxic therapies we need a better understanding of pediatric AML biology. The NUP98-KDM5A chimeric protein is exclusively found in a particular subgroup of young pediatric AML patients with complex karyotypes and poor prognosis. In this study, we investigated the impact of NUP98-KDM5A expression on cellular processes in human Pluripotent Stem Cell models and a patient-derived cell line. We found that NUP98-KDM5A generates genomic instability through two complementary mechanisms that involve accumulation of DNA damage and direct interference of RAE1 activity during mitosis. Overall, our data support that NUP98-KDM5A promotes genomic instability and likely contributes to malignant transformation.

Pan-cancer analysis revealed H3K4me1 at bivalent promoters premarks DNA hypermethylation during tumor development and identified the regulatory role of DNA methylation in relation to histone modifications

BACKGROUND

DNA hypermethylation at promoter CpG islands (CGIs) is a hallmark of cancers and could lead to dysregulation of gene expression in the development of cancers, however, its dynamics and regulatory mechanisms remain elusive. Bivalent genes, that direct development and differentiation of stem cells, are found to be frequent targets of hypermethylation in cancers.

RESULTS

Here we performed comprehensive analysis across multiple cancer types and identified that the decrease in H3K4me1 levels coincides with DNA hypermethylation at the bivalent promoter CGIs during tumorigenesis. Removal of DNA hypermethylation leads to increment of H3K4me1 at promoter CGIs with preference for bivalent genes. Nevertheless, the alteration of H3K4me1 by overexpressing or knockout LSD1, the demethylase of H3K4, doesn't change the level or pattern of DNA methylation. Moreover, LSD1 was found to regulate the expression of a bivalent gene OVOL2 to promote tumorigenesis. Knockdown of OVOL2 in LSD1 knockout HCT116 cells restored the cancer cell phenotype.

CONCLUSION

In summary, our work identified a universal indicator that can pre-mark DNA hypermethylation in cancer cells, and dissected the interplay between H3K4me1 and DNA hypermethylation in detail. Current study also reveals a novel mechanism underlying the oncogenic role of LSD1, providing clues for cancer therapies.

Bat pluripotent stem cells reveal unusual entanglement between host and viruses

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

The Role of Mitotic Slippage in Creating a "Female Pregnancy-like System" in a Single Polyploid Giant Cancer Cell

In our recent work, we observed that triple-negative breast cancer MDA-MB-231 cells respond to doxorubicin (DOX) via "mitotic slippage" (MS), discarding cytosolic damaged DNA during the process that provides their resistance to this genotoxic treatment. We also noted two populations of polyploid giant cells: those budding surviving offspring, versus those reaching huge ploidy by repeated MS and persisting for several weeks. Their separate roles in the recovery from treatment remained unclear. The current study was devoted to characterising the origin and relationship of these two sub-populations in the context of MS. MS was hallmarked by the emergence of nuclear YAP1/OCT4A/MOS/EMI2-positivity featuring a soma-germ transition to the meiotic-metaphase-arrested "maternal germ cell". In silico, the link between modules identified in the inflammatory innate immune response to cytosolic DNA and the reproductive module of female pregnancy (upregulating placenta developmental genes) was observed in polyploid giant cells. Asymmetry of the two subnuclei types, one repairing DNA and releasing buds enriched by CDC42/ACTIN/TUBULIN and the other persisting and degrading DNA in a polyploid giant cell, was revealed. We propose that when arrested in MS, a "maternal cancer germ cell" may be parthenogenetically stimulated by the placental proto-oncogene parathyroid-hormone-like-hormone, increasing calcium, thus creating a "female pregnancy-like" system within a single polyploid giant cancer cell.

Gradistics: An underappreciated dimension in evolutionary space

The growth of complexity is an unsolved and underappreciated problem. We consider possible causes of this growth, hypotheses testing, molecular mechanisms, complexity measures, cases of simplification, and significance for biomedicine. We focus on a general ability of regulation, which is based on the growing information storage and processing capacities, as the main proxy of complexity. Natural selection is indifferent to complexity. However, complexification can be inferred from the same first principle, on which natural selection is founded. Natural selection depends on potentially unlimited reproduction under limited environmental conditions. Because of the demographic pressure, the simple ecological niches become fulfilled and diversified (due to species splitting and divergence). Diversification increases complexity of biocenoses. After the filling and diversification of simple niches, the more complex niches can arise. This is the 'atomic orbitals' (AO) model. Complexity has many shortcomings but it has an advantage. This advantage is ability to regulatory adaptation, including behavioral, formed in the evolution by means of genetic adaptation. Regulatory adaptation is much faster than genetic one because it is based on the information previously accumulated via genetic adaptation and learning. Regulatory adaptation further increases complexity of biocenoses. This is the 'regulatory advantage' (RA) model. The comparison of both models allows testable predictions. We focus on the animal evolution because of the appearance of higher regulatory level (nervous system), which is absent in other lineages, and relevance to humans (including biomedical aspects).

Inheritance of paternal lifestyles and exposures through sperm DNA methylation

Many different lifestyle factors and chemicals present in the environment are a threat to the reproductive tracts of humans. The potential for parental preconception exposure to alter gametes and for these alterations to be passed on to offspring and negatively affect embryo growth and development is of concern. The connection between maternal exposures and offspring health is a frequent focus in epidemiological studies, but paternal preconception exposures are much less frequently considered and are also very important determinants of offspring health. Several environmental and lifestyle factors in men have been found to alter sperm epigenetics, which can regulate gene expression during early embryonic development. Epigenetic information is thought to be a mechanism that evolved for organisms to pass on information about their lived experiences to offspring. DNA methylation is a well-studied epigenetic regulator that is sensitive to environmental exposures in somatic cells and sperm. The continuous production of sperm from spermatogonial stem cells throughout a man's adult life and the presence of spermatogonial stem cells outside of the blood-testis barrier makes them susceptible to environmental insults. Furthermore, altered sperm DNA methylation patterns can be maintained throughout development and ultimately result in impairments, which could predispose offspring to disease. Innovations in human stem cell-based spermatogenic models can be used to elucidate the paternal origins of health and disease.

Cannabis alters DNA methylation at maternally imprinted and autism candidate genes in spermatogenic cells

Cannabis use in the United States is increasing, with highest consumption among men at their peak reproductive years. We previously demonstrated widespread changes in sperm DNA methylation with cannabis exposure in humans and rats, including genes important in neurodevelopment. Here, we use an in vitro human spermatogenesis model to recapitulate chronic cannabis use and assess DNA methylation at imprinted and autism spectrum disorder (ASD) candidate genes in spermatogonial stem cell (SSC)- and spermatid-like cells. Methylation at maternally imprinted genes SGCE and GRB10 was significantly altered in SSC- and spermatid-like cells, respectively, while PEG3 was significantly differentially methylated in spermatid-like cells. Two of ten randomly selected ASD candidate genes, HCN1 and NR4A2, had significantly altered methylation with cannabis exposure in SSC-like cells. These results support our findings in human cohorts and provide a new tool with which to gain mechanistic insights into the association between paternal cannabis use and risk of ASD in offspring.

Polyploidy and Myc Proto-Oncogenes Promote Stress Adaptation via Epigenetic Plasticity and Gene Regulatory Network Rewiring

Polyploid cells demonstrate biological plasticity and stress adaptation in evolution; development; and pathologies, including cardiovascular diseases, neurodegeneration, and cancer. The nature of ploidy-related advantages is still not completely understood. Here, we summarize the literature on molecular mechanisms underlying ploidy-related adaptive features. Polyploidy can regulate gene expression via chromatin opening, reawakening ancient evolutionary programs of embryonality. Chromatin opening switches on genes with bivalent chromatin domains that promote adaptation via rapid induction in response to signals of stress or morphogenesis. Therefore, stress-associated polyploidy can activate Myc proto-oncogenes, which further promote chromatin opening. Moreover, Myc proto-oncogenes can trigger polyploidization de novo and accelerate genome accumulation in already polyploid cells. As a result of these cooperative effects, polyploidy can increase the ability of cells to search for adaptive states of cellular programs through gene regulatory network rewiring. This ability is manifested in epigenetic plasticity associated with traits of stemness, unicellularity, flexible energy metabolism, and a complex system of DNA damage protection, combining primitive error-prone unicellular repair pathways, advanced error-free multicellular repair pathways, and DNA damage-buffering ability. These three features can be considered important components of the increased adaptability of polyploid cells. The evidence presented here contribute to the understanding of the nature of stress resistance associated with ploidy and may be useful in the development of new methods for the prevention and treatment of cardiovascular and oncological diseases.

m6A Topological Transition Coupled to Developmental Regulation of Gene Expression During Mammalian Tissue Development

N6-methyladenosine (m6A) is the most prevalent internal modification and reversible epitranscriptomic mark in messenger RNAs (mRNAs) and plays essential roles in a variety of biological processes. However, the dynamic distribution patterns of m6A and their significance during mammalian tissue development are poorly understood. Here, we found that based on m6A distribution patterns, protein-coding genes were classified into five groups with significantly distinct biological features and functions. Strikingly, comparison of the m6A methylomes of multiple mammalian tissues between fetal and adult stages revealed dynamic m6A topological transition during mammalian tissue development, and identified large numbers of genes with significant m6A loss in 5′UTRs or m6A gain around stop codons. The genes with m6A loss in 5′UTRs were highly enriched in developmental stage-specific genes, and their m6A topological transitions were strongly associated with gene expression regulation during tissue development. The genes with m6A gain around the stop codons were associated with tissue-specific functions. Our findings revealed the existence of different m6A topologies among protein-coding genes that were associated with distinct characteristics. More importantly, these genes with m6A topological transitions were crucial for tissue development via regulation of gene expression, suggesting the importance of dynamic m6A topological transitions during mammalian tissue development.

Polycomb repressive complex 2 shields naïve human pluripotent cells from trophectoderm differentiation

The first lineage choice in human embryo development separates trophectoderm from the inner cell mass. Naïve human embryonic stem cells are derived from the inner cell mass and offer possibilities to explore how lineage integrity is maintained. Here, we discover that polycomb repressive complex 2 (PRC2) maintains naïve pluripotency and restricts differentiation to trophectoderm and mesoderm lineages. Through quantitative epigenome profiling, we found that a broad gain of histone H3 lysine 27 trimethylation (H3K27me3) is a distinct feature of naïve pluripotency. We define shared and naïve-specific bivalent promoters featuring PRC2-mediated H3K27me3 concomitant with H3K4me3. Naïve bivalency maintains key trophectoderm and mesoderm transcription factors in a transcriptionally poised state. Inhibition of PRC2 forces naïve human embryonic stem cells into an ‘activated’ state, characterized by co-expression of pluripotency and lineage-specific transcription factors, followed by differentiation into either trophectoderm or mesoderm lineages. In summary, PRC2-mediated repression provides a highly adaptive mechanism to restrict lineage potential during early human development.

Two side-by-side papers report that H3K27me3 deposited by polycomb repressive complex 2 represents an epigenetic barrier that restricts naïve human pluripotent cell differentiation into alternative lineages including trophoblasts.

Polyploidy as a Fundamental Phenomenon in Evolution, Development, Adaptation and Diseases

DNA replication during cell proliferation is 'vertical' copying, which reproduces an initial amount of genetic information. Polyploidy, which results from whole-genome duplication, is a fundamental complement to vertical copying. Both organismal and cell polyploidy can emerge via premature cell cycle exit or via cell-cell fusion, the latter giving rise to polyploid hybrid organisms and epigenetic hybrids of somatic cells. Polyploidy-related increase in biological plasticity, adaptation, and stress resistance manifests in evolution, development, regeneration, aging, oncogenesis, and cardiovascular diseases. Despite the prevalence in nature and importance for medicine, agri- and aquaculture, biological processes and epigenetic mechanisms underlying these fundamental features largely remain unknown. The evolutionarily conserved features of polyploidy include activation of transcription, response to stress, DNA damage and hypoxia, and induction of programs of morphogenesis, unicellularity, and longevity, suggesting that these common features confer adaptive plasticity, viability, and stress resistance to polyploid cells and organisms. By increasing cell viability, polyploidization can provide survival under stressful conditions where diploid cells cannot survive. However, in somatic cells it occurs at the expense of specific function, thus promoting developmental programming of adult cardiovascular diseases and increasing the risk of cancer. Notably, genes arising via evolutionary polyploidization are heavily involved in cancer and other diseases. Ploidy-related changes of gene expression presumably originate from chromatin modifications and the derepression of bivalent genes. The provided evidence elucidates the role of polyploidy in evolution, development, aging, and carcinogenesis, and may contribute to the development of new strategies for promoting regeneration and preventing cardiovascular diseases and cancer.

Differential brain ADRA2A and ADRA2C gene expression and epigenetic regulation in schizophrenia. Effect of antipsychotic drug treatment

Postsynaptic α_2A-adrenoceptor density is enhanced in the dorsolateral prefrontal cortex (DLPFC) of antipsychotic-treated schizophrenia subjects. This alteration might be due to transcriptional activation, and could be regulated by epigenetic mechanisms such as histone posttranslational modifications (PTMs). The aim of this study was to evaluate ADRA2A and ADRA2C gene expression (codifying for α₂-adrenoceptor subtypes), and permissive and repressive histone PTMs at gene promoter regions in the DLPFC of subjects with schizophrenia and matched controls (n = 24 pairs). We studied the effect of antipsychotic (AP) treatment in AP-free (n = 12) and AP-treated (n = 12) subgroups of schizophrenia subjects and in rats acutely and chronically treated with typical and atypical antipsychotics. ADRA2A mRNA expression was selectively upregulated in AP-treated schizophrenia subjects (+93%) whereas ADRA2C mRNA expression was upregulated in all schizophrenia subjects (+53%) regardless of antipsychotic treatment. Acute and chronic clozapine treatment in rats did not alter brain cortex Adra2a mRNA expression but increased Adra2c mRNA expression. Both ADRA2A and ADRA2C promoter regions showed epigenetic modification by histone methylation and acetylation in human DLPFC. The upregulation of ADRA2A expression in AP-treated schizophrenia subjects might be related to observed bivalent chromatin at ADRA2A promoter region in schizophrenia (depicted by increased permissive H3K4me3 and repressive H3K27me3) and could be triggered by the enhanced H4K16ac at ADRA2A promoter. In conclusion, epigenetic predisposition differentially modulated ADRA2A and ADRA2C mRNA expression in DLPFC of schizophrenia subjects.

Growth of Biological Complexity from Prokaryotes to Hominids Reflected in the Human Genome

The growth of complexity in evolution is a most intriguing phenomenon. Using gene phylostratigraphy, we showed this growth (as reflected in regulatory mechanisms) in the human genome, tracing the path from prokaryotes to hominids. Generally, the different regulatory gene families expanded at different times, yet only up to the Euteleostomi (bony vertebrates). The only exception was the expansion of transcription factors (TF) in placentals; however, we argue that this was not related to increase in general complexity. Surprisingly, although TF originated in the Prokaryota while chromatin appeared only in the Eukaryota, the expansion of epigenetic factors predated the expansion of TF. Signaling receptors, tumor suppressors, oncogenes, and aging- and disease-associated genes (indicating vulnerabilities in terms of complex organization and strongly enrichment in regulatory genes) also expanded only up to the Euteleostomi. The complexity-related gene properties (protein size, number of alternative splicing mRNA, length of untranslated mRNA, number of biological processes per gene, number of disordered regions in a protein, and density of TF–TF interactions) rose in multicellular organisms and declined after the Euteleostomi, and possibly earlier. At the same time, the speed of protein sequence evolution sharply increased in the genes that originated after the Euteleostomi. Thus, several lines of evidence indicate that molecular mechanisms of complexity growth were changing with time, and in the phyletic lineage leading to humans, the most salient shift occurred after the basic vertebrate body plan was fixed with bony skeleton. The obtained results can be useful for evolutionary medicine.

The regulatory function of tandem repeat VNTR2-1 in hTERT gene involves basic Helix-loop-helix family transcription factors

Telomerase is a ribonucleoprotein enzyme that synthesizes telomere end sequence. The expression of hTERT gene, encoding the catalytic subunit of human telomerase, is restricted to highly proliferative tissues and is undetectable in most somatic cells. Abnormal activation of hTERT gene is found in 90% of human tumors. Previously, we identified tandem repeat of 42-bp/unit, VNTR2-1, in intron 2 of the hTERT gene, as a novel regulatory element important for hTERT transcription in cancer cells. In the current study, we found that multiple 42-bp repeats of VNTR2-1 activated luciferase gene in reporter plasmids. Mutation of the predicted cis-regulatory elements within the 42-bp repeats, including a E-box motif, resulted in a partial or complete loss of its enhancer activity. Moreover, MYC family proteins, c-MYC, MAX, and MNT, regulated hTERT gene transcription through both VNTR2-1 and E-boxes at the proximal hTERT promoter. Chromatin segmentation analysis of published ChIP-sequencing data from K562 cells indicated that VNTR2-1 was a bivalent enhancer. In telomerase-expressing human melanoma cell line MelJuSo, deletion of VNTR2-1 caused the hTERT promoter chromatin status to change from an active state to a repressed state, accompanied by increases of H3K27me3 and H3K9me3 marks. Therefore, we provided additional evidence for VNTR2-1 as a functional regulatory element that regulated hTERT expression by MYC family transcription factors. These results have improved our knowledge on the functions of repetitive genomic DNAs and the regulatory mechanisms of human telomerase gene.

Widespread overexpression from the four DNA hypermethylated HOX clusters in aggressive (IDHwt) glioma is associated with H3K27me3 depletion and alternative promoter usage

In human, the 39 coding HOX genes and 18 referenced noncoding antisense transcripts are arranged in four genomic clusters named HOXA, B, C, and D. This highly conserved family belongs to the homeobox class of genes that encode transcription factors required for normal development. Therefore, HOX gene deregulation might contribute to the development of many cancer types. Here, we study HOX gene deregulation in adult glioma, a common type of primary brain tumor. We performed extensive molecular analysis of tumor samples, classified according to their isocitrate dehydrogenase (IDH1) gene mutation status, and of glioma stem cells. We found widespread expression of sense and antisense HOX transcripts only in aggressive (IDHwt) glioma samples, although the four HOX clusters displayed DNA hypermethylation. Integrative analysis of expression, DNA methylation, and histone modification signatures along the clusters revealed that HOX gene upregulation relies on canonical and alternative bivalent CpG island promoters that escape hypermethylation. H3K27me3 loss at these promoters emerges as the main cause of widespread HOX gene upregulation in IDHwt glioma cell lines and tumors. Our study provides the first comprehensive description of the epigenetic changes at HOX clusters and their contribution to the transcriptional changes observed in adult glioma. It also identified putative 'master' HOX proteins that might contribute to the tumorigenic potential of glioma stem cells.

Stable DNMT3L overexpression in SH-SY5Y neurons recreates a facet of the genome-wide Down syndrome DNA methylation signature

Background

Down syndrome (DS) is characterized by a genome-wide profile of differential DNA methylation that is skewed towards hypermethylation in most tissues, including brain, and includes pan-tissue differential methylation. The molecular mechanisms involve the overexpression of genes related to DNA methylation on chromosome 21. Here, we stably overexpressed the chromosome 21 gene DNA methyltransferase 3L (DNMT3L) in the human SH-SY5Y neuroblastoma cell line and assayed DNA methylation at over 26 million CpGs by whole genome bisulfite sequencing (WGBS) at three different developmental phases (undifferentiated, differentiating, and differentiated).

Results

DNMT3L overexpression resulted in global CpG and CpG island hypermethylation as well as thousands of differentially methylated regions (DMRs). The DNMT3L DMRs were skewed towards hypermethylation and mapped to genes involved in neurodevelopment, cellular signaling, and gene regulation. Consensus DNMT3L DMRs showed that cell lines clustered by genotype and then differentiation phase, demonstrating sets of common genes affected across neuronal differentiation. The hypermethylated DNMT3L DMRs from all pairwise comparisons were enriched for regions of bivalent chromatin marked by H3K4me3 as well as differentially methylated sites from previous DS studies of diverse tissues. In contrast, the hypomethylated DNMT3L DMRs from all pairwise comparisons displayed a tissue-specific profile enriched for regions of heterochromatin marked by H3K9me3 during embryonic development.

Conclusions

Taken together, these results support a mechanism whereby regions of bivalent chromatin that lose H3K4me3 during neuronal differentiation are targeted by excess DNMT3L and become hypermethylated. Overall, these findings demonstrate that DNMT3L overexpression during neurodevelopment recreates a facet of the genome-wide DS DNA methylation signature by targeting known genes and gene clusters that display pan-tissue differential methylation in DS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00387-7.

The Role of H3K4 Trimethylation in CpG Islands Hypermethylation in Cancer

CpG methylation in transposons, exons, introns and intergenic regions is important for long-term silencing, silencing of parasitic sequences and alternative promoters, regulating imprinted gene expression and determining X chromosome inactivation. Promoter CpG islands, although rich in CpG dinucleotides, are unmethylated and remain so during all phases of mammalian embryogenesis and development, except in specific cases. The biological mechanisms that contribute to the maintenance of the unmethylated state of CpG islands remain elusive, but the modification of established DNA methylation patterns is a common feature in all types of tumors and is considered as an event that intrinsically, or in association with genetic lesions, feeds carcinogenesis. In this review, we focus on the latest results describing the role that the levels of H3K4 trimethylation may have in determining the aberrant hypermethylation of CpG islands in tumors.

Differentiating cancer cells reveal early large-scale genome regulation by pericentric domains

Finding out how cells prepare for fate change during differentiation commitment was our task. To address whether the constitutive pericentromere-associated domains (PADs) may be involved, we used a model system with known transcriptome data, MCF-7 breast cancer cells treated with the ErbB3 ligand heregulin (HRG), which induces differentiation and is used in the therapy of cancer. PAD-repressive heterochromatin (H3K9me3), centromere-associated-protein-specific, and active euchromatin (H3K4me3) antibodies, real-time PCR, acridine orange DNA structural test (AOT), and microscopic image analysis were applied. We found a two-step DNA unfolding after 15–20 and 60 min of HRG treatment, respectively. This behavior was consistent with biphasic activation of the early response genes (c-fos - fosL1/myc) and the timing of two transcriptome avalanches reported in the literature. In control, the average number of PADs negatively correlated with their size by scale-free distribution, and centromere clustering in turn correlated with PAD size, both indicating that PADs may create and modulate a suprachromosomal network by fusing and splitting a constant proportion of the constitutive heterochromatin. By 15 min of HRG treatment, the bursting unraveling of PADs from the nucleolus boundary occurred, coinciding with the first step of H3K4me3 chromatin unfolding, confirmed by AOT. The second step after 60 min of HRG treatment was associated with transcription of long noncoding RNA from PADs and peaking of fosL1/c-myc response. We hypothesize that the bursting of PAD clusters under a critical silencing threshold pushes the first transcription avalanche, whereas the destruction of the PAD network enables genome rewiring needed for differentiation repatterning, mediated by early response bivalent genes.

DNA Methylation Differences Between Zona Pellucida-Bound and Manually Selected Spermatozoa Are Associated With Autism Susceptibility

Children conceived through intracytoplasmic sperm injection (ICSI) have been reported to have a higher risk of many abnormalities and disorders, including autism and intellectual disability, which may be due to bypassing of the natural sperm selection process during ICSI. Zona pellucida (ZP)-bound spermatozoa (ZPBS) have normal morphology and nuclear DNA. Using these spermatozoa for ICSI results in better outcomes compared with conventional ICSI. However, differences besides morphology that exist between sperm selected by ZP and by an embryologist and whether these differences affect the risk of autism in offspring after ICSI are unclear. To explore these questions, we compared genome-wide DNA methylation profiles between ZPBS and manually selected spermatozoa (MSS)using single-cell bisulfite sequencing. Global DNA methylation levels were significantly lower in ZPBS than in MSS. Using gene ontology (GO) analysis, genes overlapping differentially methylated regions (DMRs) were enriched in biological processes involving neurogenesis. Furthermore, we found that 47.8% of autism candidate genes were associated with DMRs, compared with 37.1% of matched background genes (P<0.001). This was mainly because of the high proportion of autism candidate genes with bivalent chromatin structure. In conclusion, bivalent chromatin structure results in large differences in the methylation of autism genes between MSS and ZPBS. ICSI using MSS, which increases the risk of methylation mutations compared with ZPBS, may lead to a higher risk of autism in offspring.

Paradoxes of cancer: Survival at the brink

The fundamental understanding of how Cancer initiates, persists and then progresses is evolving. High-resolution technologies, including single-cell mutation and gene expression measurements, are now attainable, providing an ever-increasing insight into the molecular details. However, this higher resolution has shown that somatic mutation theory itself cannot explain the extraordinary resistance of cancer to extinction. There is a need for a more Systems-based framework of understanding cancer complexity, which in particular explains the regulation of gene expression during cell-fate decisions. Cancer displays a series of paradoxes. Here we attempt to approach them from the view-point of adaptive exploration of gene regulatory networks at the edge of order and chaos, where cell-fate is changed by oscillations between alternative regulators of cellular senescence and reprogramming operating through self-organisation. On this background, the role of polyploidy in accessing the phylogenetically pre-programmed "oncofetal attractor" state, related to unicellularity, and the de-selection of unsuitable variants at the brink of cell survival is highlighted. The concepts of the embryological and atavistic theory of cancer, cancer cell "life-cycle", and cancer aneuploidy paradox are dissected under this lense. Finally, we challenge researchers to consider that cancer "defects" are mostly the adaptation tools of survival programs that have arisen during evolution and are intrinsic of cancer. Recognition of these features should help in the development of more successful anti-cancer treatments.

Suppression of poised oncogenes by ZMYND8 promotes chemo-sensitization

The major challenge in chemotherapy lies in the gain of therapeutic resistance properties of cancer cells. The relatively small fraction of chemo-resistant cancer cells outgrows and are responsible for tumor relapse, with acquired invasiveness and stemness. We demonstrate that zinc-finger MYND type-8 (ZMYND8), a putative chromatin reader, suppresses stemness, drug resistance, and tumor-promoting genes, which are hallmarks of cancer. Reinstating ZMYND8 suppresses chemotherapeutic drug doxorubicin-induced tumorigenic potential (at a sublethal dose) and drug resistance, thereby resetting the transcriptional program of cells to the epithelial state. The ability of ZMYND8 to chemo-sensitize doxorubicin-treated metastatic breast cancer cells by downregulating tumor-associated genes was further confirmed by transcriptome analysis. Interestingly, we observed that ZMYND8 overexpression in doxorubicin-treated cells stimulated those involved in a good prognosis in breast cancer. Consistently, sensitizing the cancer cells with ZMYND8 followed by doxorubicin treatment led to tumor regression in vivo and revert back the phenotypes associated with drug resistance and stemness. Intriguingly, ZMYND8 modulates the bivalent or poised oncogenes through its association with KDM5C and EZH2, thereby chemo-sensitizing the cells to chemotherapy for better disease-free survival. Collectively, our findings indicate that poised chromatin is instrumental for the acquisition of chemo-resistance by cancer cells and propose ZMYND8 as a suitable epigenetic tool that can re-sensitize the chemo-refractory breast carcinoma.

Phylostratic Shift of Whole-Genome Duplications in Normal Mammalian Tissues towards Unicellularity Is Driven by Developmental Bivalent Genes and Reveals a Link to Cancer

Tumours were recently revealed to undergo a phylostratic and phenotypic shift to unicellularity. As well, aggressive tumours are characterized by an increased proportion of polyploid cells. In order to investigate a possible shared causation of these two features, we performed a comparative phylostratigraphic analysis of ploidy-related genes, obtained from transcriptomic data for polyploid and diploid human and mouse tissues using pairwise cross-species transcriptome comparison and principal component analysis. Our results indicate that polyploidy shifts the evolutionary age balance of the expressed genes from the late metazoan phylostrata towards the upregulation of unicellular and early metazoan phylostrata. The up-regulation of unicellular metabolic and drug-resistance pathways and the downregulation of pathways related to circadian clock were identified. This evolutionary shift was associated with the enrichment of ploidy with bivalent genes (p < 10⁻¹⁶). The protein interactome of activated bivalent genes revealed the increase of the connectivity of unicellulars and (early) multicellulars, while circadian regulators were depressed. The mutual polyploidy-c-MYC-bivalent genes-associated protein network was organized by gene-hubs engaged in both embryonic development and metastatic cancer including driver (proto)-oncogenes of viral origin. Our data suggest that, in cancer, the atavistic shift goes hand-in-hand with polyploidy and is driven by epigenetic mechanisms impinging on development-related bivalent genes.

Bivalent genes that undergo transcriptional switching identify networks of key regulators of embryonic stem cell differentiation

Background

Bivalent promoters marked with both H3K27me3 and H3K4me3 histone modifications are characteristic of poised promoters in embryonic stem (ES) cells. The model of poised promoters postulates that bivalent chromatin in ES cells is resolved to monovalency upon differntiation. With the availability of single-cell RNA sequencing (scRNA-seq) data, subsequent switches in transcriptional state at bivalent promoters can be studied more closely.

Results

We develop an approach for capturing genes undergoing transcriptional switching by detecting ‘bimodal’ gene expression patterns from scRNA-seq data. We integrate the identification of bimodal genes in ES cell differentiation with analysis of chromatin state, and identify clear cell-state dependent patterns of bimodal, bivalent genes. We show that binarization of bimodal genes can be used to identify differentially expressed genes from fractional ON/OFF proportions. In time series data from differentiating cells, we build a pseudotime approximation and use a hidden Markov model to infer gene activity switching pseudotimes, which we use to infer a regulatory network. We identify pathways of switching during differentiation, novel details of those pathway, and transcription factor coordination with downstream targets.

Conclusions

Genes with expression levels too low to be informative in conventional scRNA analysis can be used to infer transcriptional switching networks that connect transcriptional activity to chromatin state. Since chromatin bivalency is a hallmark of gene promoters poised for activity, this approach provides an alternative that complements conventional scRNA-seq analysis while focusing on genes near the ON/OFF boundary of activity. This offers a novel and productive means of inferring regulatory networks from scRNA-seq data.

Globally altered epigenetic landscape and delayed osteogenic differentiation in H3.3-G34W-mutant giant cell tumor of bone

The neoplastic stromal cells of giant cell tumor of bone (GCTB) carry a mutation in H3F3A, leading to a mutant histone variant, H3.3-G34W, as a sole recurrent genetic alteration. We show that in patient-derived stromal cells H3.3-G34W is incorporated into the chromatin and associates with massive epigenetic alterations on the DNA methylation, chromatin accessibility and histone modification level, that can be partially recapitulated in an orthogonal cell line system by the introduction of H3.3-G34W. These epigenetic alterations affect mainly heterochromatic and bivalent regions and provide possible explanations for the genomic instability, as well as the osteolytic phenotype of GCTB. The mutation occurs in differentiating mesenchymal stem cells and associates with an impaired osteogenic differentiation. We propose that the observed epigenetic alterations reflect distinct differentiation stages of H3.3 WT and H3.3 MUT stromal cells and add to H3.3-G34W-associated changes.

The histone variant mutation H3.3-G34W occurs in the majority of giant cell tumor of bone (GCTB). By profiling patient-derived GCTB tumor cells, the authors show that this mutation associates with epigenetic alterations in heterochromatic and bivalent regions that contribute to an impaired osteogenic differentiation and the osteolytic phenotype of GCTB.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

Sperm DNA methylation altered by THC and nicotine: Vulnerability of neurodevelopmental genes with bivalent chromatin

Men consume the most nicotine and cannabis products but impacts on sperm epigenetics are poorly characterized. Evidence suggests that preconception exposure to these drugs alters offspring neurodevelopment. Epigenetics may in part facilitate heritability. We therefore compared effects of exposure to tetrahydrocannabinol (THC) and nicotine on DNA methylation in rat sperm at genes involved in neurodevelopment. Reduced representation bisulfite sequencing data from sperm of rats exposed to THC via oral gavage showed that seven neurodevelopmentally active genes were significantly differentially methylated versus controls. Pyrosequencing data revealed majority overlap in differential methylation in sperm from rats exposed to THC via injection as well as those exposed to nicotine. Neurodevelopmental genes including autism candidates are vulnerable to environmental exposures and common features may mediate this vulnerability. We discovered that autism candidate genes are significantly enriched for bivalent chromatin structure, suggesting this configuration may increase vulnerability of genes in sperm to disrupted methylation.

Dynamic Signatures of the Epigenome: Friend or Foe?

Highly dynamic epigenetic signaling is influenced mainly by (micro)environmental stimuli and genetic factors. The exact mechanisms affecting particular epigenomic patterns differ dependently on the context. In the current review, we focus on the causes and effects of the dynamic signatures of the human epigenome as evaluated with the high-throughput profiling data and single-gene approaches. We will discuss three different aspects of phenotypic outcomes occurring as a consequence of epigenetics interplaying with genotype and environment. The first issue is related to the cases of environmental impacts on epigenetic profile, and its adverse and advantageous effects related to human health and evolutionary adaptation. The next topic will present a model of the interwoven co-evolution of genetic and epigenetic patterns exemplified with transposable elements (TEs) and their epigenetic repressors Krüppel-associated box zinc finger proteins (KRAB–ZNFs). The third aspect concentrates on the mitosis-based microevolution that takes place during carcinogenesis, leading to clonal diversity and expansion of tumor cells. The whole picture of epigenome plasticity and its role in distinct biological processes is still incomplete. However, accumulating data define epigenomic dynamics as an essential co-factor driving adaptation at the cellular and inter-species levels with a benefit or disadvantage to the host.

Modifications of H3K4 methylation levels are associated with DNA hypermethylation in acute myeloid leukemia

The 'instructive model' of aberrant DNA methylation in human tumors is based on the observation that CpG islands prone to hypermethylation in cancers are embedded in chromatin enriched in H3K27me3 in human embryonic stem cells (hESC). Recent studies also link methylation of CpG islands to the methylation status of H3K4, where H3K4me3 is inversely correlated with DNA methylation. To provide insight into these conflicting findings, we generated DNA methylation profiles for acute myeloid leukemia samples from patients and leukemic cell lines and integrated them with publicly available ChIp-seq data, containing H3K4me3 and H3K27me3 CpG island occupation in hESC, or hematopoietic stem or progenitor cells (hHSC/MPP). Hypermethylated CpG islands in AML samples displayed H3K27me3 enrichments in hESC and hHSC/MPP; however, ChIp analysis of specific hypermethylated CpG islands revealed a significant reduction in H3K4me3 signal with a concomitant increase in H3K4me0 levels as opposed to a nonsignificant increase in H3K27me3 marks. The integration of AML DNA methylation profiles with the ChIp-seq data in hESC and hHSC/MPP also led to the identification of Iroquois homeobox 2 (IRX2) as a previously unknown factor promoting differentiation of leukemic cells. Our results indicate that in contrast to the 'instructive model', H3K4me3 levels are strongly associated with DNA methylation patterns in AML and have a role in the regulation of critical genes, such as the putative tumor suppressor IRX2.

Transcriptional alterations in glioma result primarily from DNA methylation-independent mechanisms

In cancer cells, aberrant DNA methylation is commonly associated with transcriptional alterations, including silencing of tumor suppressor genes. However, multiple epigenetic mechanisms, including polycomb repressive marks, contribute to gene deregulation in cancer. To dissect the relative contribution of DNA methylation–dependent and –independent mechanisms to transcriptional alterations at CpG island/promoter-associated genes in cancer, we studied 70 samples of adult glioma, a widespread type of brain tumor, classified according to their isocitrate dehydrogenase (IDH1) mutation status. We found that most transcriptional alterations in tumor samples were DNA methylation–independent. Instead, altered histone H3 trimethylation at lysine 27 (H3K27me3) was the predominant molecular defect at deregulated genes. Our results also suggest that the presence of a bivalent chromatin signature at CpG island promoters in stem cells predisposes not only to hypermethylation, as widely documented, but more generally to all types of transcriptional alterations in transformed cells. In addition, the gene expression strength in healthy brain cells influences the choice between DNA methylation- and H3K27me3-associated silencing in glioma. Highly expressed genes were more likely to be repressed by H3K27me3 than by DNA methylation. Our findings support a model in which altered H3K27me3 dynamics, more specifically defects in the interplay between polycomb protein complexes and the brain-specific transcriptional machinery, is the main cause of transcriptional alteration in glioma cells. Our study provides the first comprehensive description of epigenetic changes in glioma and their relative contribution to transcriptional changes. It may be useful for the design of drugs targeting cancer-related epigenetic defects.

DNA hypermethylation in disease: mechanisms and clinical relevance

Increasing numbers of studies implicate abnormal DNA methylation in cancer and many non-malignant diseases. This is consistent with numerous findings about differentiation-associated changes in DNA methylation at promoters, enhancers, gene bodies, and sites that control higher-order chromatin structure. Abnormal increases or decreases in DNA methylation contribute to or are markers for cancer formation and tumour progression. Aberrant DNA methylation is also associated with neurological diseases, immunological diseases, atherosclerosis, and osteoporosis. In this review, I discuss DNA hypermethylation in disease and its interrelationships with normal development as well as proposed mechanisms for the origin of and pathogenic consequences of disease-associated hypermethylation. Disease-linked DNA hypermethylation can help drive oncogenesis partly by its effects on cancer stem cells and by the CpG island methylator phenotype (CIMP); atherosclerosis by disease-related cell transdifferentiation; autoimmune and neurological diseases through abnormal perturbations of cell memory; and diverse age-associated diseases by age-related accumulation of epigenetic alterations.

Development of an Interactive Web Application "Shiny App for Frequency Analysis on Homo sapiens Genome (SAFA-HsG)"

The web application "Shiny App for Frequency Analysis on Homo sapiens Genome (SAFA-HsG)" was developed using R programming-based bioconductor packages and shiny framework. Through the app, preliminary descriptive data analysis on nucleotide frequency, and CpG island, CpG non-island, and CpG island shores and shelves (downstream and upstream) of human reference genome can be carried out, which will help biologists to work on human epigenomics. Table view of these analyses of all chromosomes can be visualized and downloaded by the end users. Similarly, the respective comparative plots can be used for CpG sites comparison. In addition, to introduce the personal genome project, the present study has done a preliminary work on few raw data and are included in the app, which will create interest on personal genome information. The app is hosted on https://SAFA-HsG.shinyapps.io/home/. It is a multi-platform application and can be initiated locally from any computer that has or has not installed R. It is a user-friendly interface, which will allow a biologist, even who has little computer knowledge to access and analyze further.

Targets and genomic constraints of ectopic Dnmt3b expression

DNA methylation plays an essential role in mammalian genomes and expression of the responsible enzymes is tightly controlled. Deregulation of the de novo DNA methyltransferase DNMT3B is frequently observed across cancer types, yet little is known about its ectopic genomic targets. Here, we used an inducible transgenic mouse model to delineate rules for abnormal DNMT3B targeting, as well as the constraints of its activity across different cell types. Our results explain the preferential susceptibility of certain CpG islands to aberrant methylation and point to transcriptional state and the associated chromatin landscape as the strongest predictors. Although DNA methylation and H3K27me3 are usually non-overlapping at CpG islands, H3K27me3 can transiently co-occur with DNMT3B-induced DNA methylation. Our genome-wide data combined with ultra-deep locus-specific bisulfite sequencing suggest a distributive activity of ectopically expressed Dnmt3b that leads to discordant CpG island hypermethylation and provides new insights for interpreting the cancer methylome.

DNA methylation profiling reveals a pathological signature that contributes to transcriptional defects of CD34+ CD15- cells in early chronic-phase chronic myeloid leukemia

Despite the high efficiency of tyrosine kinase inhibitors (TKI), some patients with chronic myeloid leukemia (CML) will display residual disease that can become resistant to treatment, indicating intraclonal heterogeneity in chronic‐phase CML (CP‐CML). To determine the basis of this heterogeneity, we conducted the first exhaustive characterization of the DNA methylation pattern of sorted CP‐CML CD34⁺CD15⁻ (immature) and CD34⁻CD15⁺ (mature) cells at diagnosis (prior to any treatment) and compared it to that of CD34⁺CD15⁻ and CD34⁻CD15⁺ cells isolated from healthy donors (HD). In both cell types, we identified several hundreds of differentially methylated regions (DMRs) showing DNA methylation changes between CP‐CML and HD samples, with only a subset of them in common between CD34⁺CD15⁻ and CD34⁻CD15⁺ cells. This suggested DNA methylation variability within the same CML clone. We also identified 70 genes that could be aberrantly repressed upon hypermethylation and 171 genes that could be aberrantly expressed upon hypomethylation of some of these DMRs in CP‐CML cells, among which 18 and 81, respectively, were in CP‐CML CD34⁺CD15⁻ cells only. We then validated the DNA methylation and expression defects of selected candidate genes. Specifically, we identified GAS2, a candidate oncogene, as a new example of gene the hypomethylation of which is associated with robust overexpression in CP‐CML cells. Altogether, we demonstrated that DNA methylation abnormalities exist at early stages of CML and can affect the transcriptional landscape of malignant cells. These observations could lead to the development of combination treatments with epigenetic drugs and TKI for CP‐CML.

Bivalent complexes of PRC1 with orthologs of BRD4 and MOZ/MORF target developmental genes in Drosophila

Kang et al. confirm PRC1–Br140 and PRC1–Fs(1)h interactions and identify their genomic binding sites. PRC1–Br140 bind developmental genes in fly embryos, with analogous co-occupancy of PRC1 and BRD1 at bivalent loci in human ES cells.

Regulatory decisions in Drosophila require Polycomb group (PcG) proteins to maintain the silent state and Trithorax group (TrxG) proteins to oppose silencing. Since PcG and TrxG are ubiquitous and lack apparent sequence specificity, a long-standing model is that targeting occurs via protein interactions; for instance, between repressors and PcG proteins. Instead, we found that Pc-repressive complex 1 (PRC1) purifies with coactivators Fs(1)h [female sterile (1) homeotic] and Enok/Br140 during embryogenesis. Fs(1)h is a TrxG member and the ortholog of BRD4, a bromodomain protein that binds to acetylated histones and is a key transcriptional coactivator in mammals. Enok and Br140, another bromodomain protein, are orthologous to subunits of a mammalian MOZ/MORF acetyltransferase complex. Here we confirm PRC1–Br140 and PRC1–Fs(1)h interactions and identify their genomic binding sites. PRC1–Br140 bind developmental genes in fly embryos, with analogous co-occupancy of PRC1 and a Br140 ortholog, BRD1, at bivalent loci in human embryonic stem (ES) cells. We propose that identification of PRC1–Br140 “bivalent complexes” in fly embryos supports and extends the bivalency model posited in mammalian cells, in which the coexistence of H3K4me3 and H3K27me3 at developmental promoters represents a poised transcriptional state. We further speculate that local competition between acetylation and deacetylation may play a critical role in the resolution of bivalent protein complexes during development.