Trimethylation at histone H3 lysine 4 (H3K4me3) and lysine 27 (H3K27me3) controls gene activity during development and differentiation. Whether H3K4me3 and H3K27me3 changes dynamically in response to altered microenvironmental conditions, including low-oxygen conditions commonly present in solid tumors, is relatively unknown. Demethylation of H3K4me3 and H3K27me3 is mediated by oxygen and 2-oxoglutarate dioxygenases enzymes, suggesting that oxygen deprivation (hypoxia) may influence histone trimethylation. Using the MCF7 breast epithelial adenocarcinoma cell model, we have determined the relationship between epigenomic and transcriptomic reprogramming as a function of fluctuating oxygen tension. We find that in MCF7, H3K4me3 and H3K27me3 marks rapidly increase at specific locations throughout the genome and are largely reversed upon reoxygenation. Whereas dynamic changes are relatively highest for H3K27me3 marking under hypoxic conditions, H3K4me3 occupation is identified as the defining epigenetic marker of transcriptional control. In agreement with the global increase of H3K27 trimethylation, we provide direct evidence that the histone H3K27me3 demethylase KDM6B/JMJD3 is inactivated by limited oxygen. In situ immunohistochemical analysis confirms a marked rise of histone trimethylation in hypoxic tumor areas. Acquisition of H3K27me3 at H3K4me3-marked loci results in a striking increase in “bivalent” epigenetic marking. Hypoxia-induced bivalency substantially overlaps with embryonal stem cell-associated genic bivalency and is retained at numerous loci upon reoxygenation. Transcriptional activity is selectively and progressively dampened at bivalently marked loci upon repeated exposure to hypoxia, indicating that this subset of genes uniquely maintains the potential for epigenetic regulation by KDM activity. These data suggest that dynamic regulation of the epigenetic state within the tumor environment may have important consequences for tumor plasticity and biology.

Introduction: Relatively little is known about cellular subpopulations in the mature nucleus pulposus (NP). Detailed understanding of the ontogenetic, cellular and molecular characteristics of functional intervertebral disc (IVD) cell populations is pivotal to the successful development of cell replacement therapies and IVD regeneration. In this study, we aimed to investigate whether phenotypically distinct clonal cell lines representing different subpopulations in the human NP could be generated using immortalization strategies. Methods: Nondegenerate healthy disc material (age range, 8 to 15 years) was obtained as surplus surgical material. Early passage NP monolayer cell cultures were initially characterized using a recently established NP marker set. NP cells were immortalized by simian virus 40 large T antigen (SV40LTag) and human telomerase reverse transcriptase expression. Immortalized cells were clonally expanded and characterized based on collagen type I, collagen type II, α1 (COL2A1), and SRY-box 9 (SOX9) protein expression profiles, as well as on expression of a subset of established in vivo NP cell lineage markers. Results: A total of 54 immortal clones were generated. Profiling of a set of novel NP markers (CD24, CA12, PAX1, PTN, FOXF1 and KRT19 mRNA) in a representative set of subclones substantiated successful immortalization of multiple cellular subpopulations from primary isolates and confirmed their NP origin and/or phenotype. We were able to identify two predominant clonal NP subtypes based on their morphological characteristics and their ability to induce SOX9 and COL2A1 under conventional differentiation conditions. In addition, cluster of differentiation 24 (CD24)–negative NP responder clones formed spheroid structures in various culture systems, suggesting the preservation of a more immature phenotype compared to CD24-positive nonresponder clones. Conclusions: Here we report the generation of clonal NP cell lines from nondegenerate human IVD tissue and present a detailed characterization of NP cellular subpopulations. Differential cell surface marker expression and divergent responses to differentiation conditions suggest that the NP subtypes may correspond to distinct maturation stages and represent distinct NP cell subpopulations. Hence, we provide evidence that the immortalization strategy that we applied is capable of detecting cell heterogeneity in the NP. Our cell lines yield novel insights into NP biology and provide promising new tools for studies of IVD development, cell function and disease.

Initiation of and progression through chondrogenesis is driven by changes in the cellular microenvironment. At the onset of chondrogenesis, resting mesenchymal stem cells are mobilized in vivo and a complex, step-wise chondrogenic differentiation program is initiated. Differentiation requires coordinated transcriptomic reprogramming and increased progenitor proliferation; both processes require chromatin remodeling. The nature of early molecular responses that relay differentiation signals to chromatin is poorly understood. We here show that immediate early genes are rapidly and transiently induced in response to differentiation stimuli in vitro. Functional ablation of the immediate early factor EGR1 severely deregulates expression of key chondrogenic control genes at the onset of differentiation. In addition, differentiating cells accumulate DNA damage, activate a DNA damage response and undergo a cell cycle arrest and prevent differentiation associated hyper-proliferation. Failed differentiation in the absence of EGR1 affects global acetylation and terminates in overall histone hypermethylation. We report novel molecular connections between EGR1 and Polycomb Group function: Polycomb associated histone H3 lysine27 trimethylation (H3K27me3) blocks chromatin access of EGR1. In addition, EGR1 ablation results in abnormal Ezh2 and Bmi1 expression. Consistent with this functional interaction, we identify a number of co-regulated targets genes in a chondrogenic gene network. We here describe an important role for EGR1 in early chondrogenic epigenetic programming to accommodate early gene-environment interactions in chondrogenesis.

Loss of annulus fibrosus (AF) integrity predisposes to disc herniation and is associated with IVD degeneration. Successful implementation of biomedical intervention therapy requires in-depth knowledge of IVD cell biology. We recently generated unique clonal human nucleus pulposus (NP) cell lines. Recurring functional cellular phenotypes from independent donors provided pivotal evidence for cell heterogeneity in the mature human NP. In this study we aimed to generate and characterize immortal cell lines for the human AF from matched donors. Non-degenerate healthy disc material was obtained as surplus surgical material. AF cells were immortalized by simian virus Large T antigen (SV40LTAg) and human telomerase (hTERT) expression. Early passage cells and immortalized cell clones were characterized based on marker gene expression under standardized culturing and in the presence of Transforming Growth factor β (TGFβ). The AF-specific expression signature included COL1A1, COL5A1, COL12A1, SFRP2 and was largely maintained in immortal AF cell lines. Remarkably, TGFβ induced rapid 3D sheet formation in a subgroup of AF clones. This phenotype was associated with inherent differences in Procollagen type I processing and maturation, and correlated with differential mRNA expression of Prolyl 4-hydroxylase alpha polypeptide 1 and 3 (P4HA1,3) and Lysyl oxidase (LOX) between clones and differential P4HA3 protein expression between AF cells in histological sections. We report for the first time the generation of representative human AF cell lines. Gene expression profile analysis and functional comparison of AF clones revealed variation between immortalized cells and suggests phenotypic heterogeneity in the human AF. Future characterization of AF cellular (sub-)populations aims to combine identification of additional specific AF marker genes and their biological relevance. Ultimately this knowledge will contribute to clinical application of cell-based technology in IVD repair.

A comprehensive assessment of the epigenetic dynamics in cancer cells is the key to understanding the molecular mechanisms underlying cancer and to improving cancer diagnostics, prognostics and treatment. By combining genome-wide ChIP-seq epigenomics and microarray transcriptomics, we studied the effects of oxygen deprivation and subsequent reoxygenation on histone 3 trimethylation of lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in a breast cancer cell line, serving as a model for abnormal oxygenation in solid tumors. A priori, epigenetic markings and gene expression levels not only are expected to vary greatly between hypoxic and normoxic conditions, but also display a large degree of heterogeneity across the cell population. Where traditionally ChIP-seq data are often treated as dichotomous data, the model and experiment here necessitate a quantitative, data-driven analysis of both datasets. We first identified genomic regions with sustained epigenetic markings, which provided a sample-specific reference enabling quantitative ChIP-seq data analysis. Sustained H3K27me3 marking was located around centromeres and intergenic regions, while sustained H3K4me3 marking is associated with genes involved in RNA binding, translation and protein transport and localization. Dynamic marking with both H3K4me3 and H3K27me3 (hypoxia-induced bivalency) was found in CpG-rich regions at loci encoding factors that control developmental processes, congruent with observations in embryonic stem cells. In silico-identified epigenetically sustained and dynamic genomic regions were confirmed through ChIP-PCR in vitro, and obtained results are corroborated by published data and current insights regarding epigenetic regulation.

Background: Gene-environment interactions are mediated by epigenetic mechanisms. Polycomb Group proteins constitute part of an epigenetic cellular transcriptional memory system that is subject to dynamic modulation during differentiation. Molecular insight in processes that control dynamic chromatin association and dissociation of Polycomb repressive complexes during and beyond development is limited. We recently showed that MK3 interacts with Polycomb repressive complex 1 (PRC1). The functional relevance of this interaction, however, remained poorly understood. MK3 is activated downstream of mitogen- and stress-activated protein kinases (M/SAPKs), all of which fulfill crucial roles during development. We here use activation of the immediate-early response gene ATF3, a bona fide PRC1 target gene, as a model to study how MK3 and its effector kinases MAPK/ERK and SAPK/P38 are involved in regulation of PRC1-dependent ATF3 transcription. Results: Our current data show that mitogenic signaling through ERK, P38 and MK3 regulates ATF3 expression by PRC1/chromatin dissociation and epigenetic modulation. Mitogenic stimulation results in transient P38-dependent H3S28 phosphorylation and ERK-driven PRC1/chromatin dissociation at PRC1 targets. H3S28 phosphorylation by itself appears not sufficient to induce PRC1/chromatin dissociation, nor ATF3 transcription, as inhibition of MEK/ERK signaling blocks BMI1/chromatin dissociation and ATF3 expression, despite induced H3S28 phosphorylation. In addition, we establish that concomitant loss of local H3K27me3 promoter marking is not required for ATF3 activation. We identify pERK as a novel signaling-induced binding partner of PRC1, and provide evidence that MK3 controls ATF3 expression in cultured cells via negative regulatory feedback on M/SAPKs. Dramatically increased ectopic wing vein formation in the absence of Drosophila MK in a Drosophila ERK gain-of-function wing vein patterning model, supports the existence of MK-mediated negative feedback regulation on pERK. Conclusion: We here identify and characterize important actors in a PRC1-dependent epigenetic signal/response mechanism, some of which appear to be nonspecific global responses, whereas others provide modular specificity. Our findings provide novel insight into a Polycomb-mediated epigenetic mechanism that dynamically controls gene transcription and support a direct link between PRC1 and cellular responses to changes in the microenvironment.