Genome-wide map of quantified epigenetic changes during in vitrochondrogenic differentiation of primary human mesenchymal stem cells
© Herlofsen et al.; licensee BioMed Central Ltd. 2013
Received: 3 December 2012
Accepted: 12 February 2013
Published: 15 February 2013
For safe clinical application of engineered cartilage made from mesenchymal stem cells (MSCs), molecular mechanisms for chondrogenic differentiation must be known in detail. Changes in gene expression and extracellular matrix synthesis have been extensively studied, but the epigenomic modifications underlying these changes have not been described. To this end we performed whole-genome chromatin immunoprecipitation and deep sequencing to quantify six histone modifications, reduced representation bisulphite sequencing to quantify DNA methylation and mRNA microarrays to quantify gene expression before and after 7 days of chondrogenic differentiation of MSCs in an alginate scaffold. To add to the clinical relevance of our observations, the study is based on primary bone marrow-derived MSCs from four donors, allowing us to investigate inter-individual variations.
We see two levels of relationship between epigenetic marking and gene expression. First, a large number of genes ontogenetically linked to MSC properties and the musculoskeletal system are epigenetically prepatterned by moderate changes in H3K4me3 and H3K9ac near transcription start sites. Most of these genes remain transcriptionally unaltered. Second, transcriptionally upregulated genes, more closely associated with chondrogenesis, are marked by H3K36me3 in gene bodies, highly increased H3K4me3 and H3K9ac on promoters and 5' end of genes, and increased H3K27ac and H3K4me1 marking in at least one enhancer region per upregulated gene. Within the 7-day time frame, changes in promoter DNA methylation do not correlate significantly with changes in gene expression. Inter-donor variability analysis shows high level of similarity between the donors for this data set.
Histone modifications, rather than DNA methylation, provide the primary epigenetic control of early differentiation of MSCs towards the chondrogenic lineage.
The best treatment of lesions of hyaline cartilage known to date is transplantation of in vitro expanded autologous chondrocytes . However, as articular chondrocytes dedifferentiate during in vitro culture, the resulting tissue frequently consists of a mixture of hyaline and fibrous cartilage. A better strategy may be to produce an implant of hyaline cartilage by in vitro tissue engineering using mesenchymal stem cells (MSCs) and a biomaterial. Several approaches including a variety of scaffold systems have been used to induce chondrogenic differentiation of MSCs, but so far all have failed to produce perfect hyaline cartilage for clinical use [1, 2]. To achieve this, it is important to understand the processes that regulate chondrogenic differentiation of MSCs.
Changes in epigenetic marks are known to be important regulatory factors in stem cell differentiation. Gene expression is correlated with the level of post-translational modifications of histone tails and DNA methylation of promoter regions. For embryonic stem cells (ESCs), lineage commitment is associated with repression of genes important for differentiation along other lineages [3–5]. This process frequently involves the polycomb-group mediated trimethylation of lysine 27 of histone 3 (H3K27me3) and DNA methylation of the promoter regions [6, 7]. Transcriptionally active genes are marked by monoubiquitination of H2B followed by trimethylation of H3K4 (H3K4me3) and H3K9 acetylation (H3K9ac) on nucleosomes close to the transcription start site (TSS). Other histone modifications, including H3K36me3, are associated with transcriptional elongation and occur throughout the body of transcribed genes. Additional transcriptional control is provided by cis-acting regulatory regions marked by H3K4me1 and H3K27ac [4, 5, 8, 9].
Less is known about the epigenetic control of gene expression during differentiation of multipotent stem cells. Analysis of genes regulating adipogenic differentiation of human MSCs revealed dynamic changes in histone marks reminiscent of those seen in ESCs  and modest changes in promoter DNA methylation . Recent studies of differentiation of murine myogenic precursor lines and in vivo differentiation of murine hair follicle stem cells have focused on the importance of H3K27me3 [12, 13]. Likewise, studies of human hematopoietic stem cells have focused on changes in genes concomitantly marked by H3K4me3 and H3K27me3 and changes in DNA methylation [6, 14–16]. Except for a few studies of the DNA methylation status of some genes involved during in vitro chondrogenesis [17–19], nothing is known to date about epigenetic changes during chondrogenic differentiation of MSCs.
Embryological chondrogenesis is controlled by a complex interplay of growth- and transcription factors. Among the essential secreted factors are Tgfβ and Bmp proteins . For in vitro differentiation of human MSCs a combination of TGFβ and BMP2 yields the most hyaline-like cartilage . In vivo, the Sox9 transcription factor (TF) is essential through the entire process of chondrogenesis  but the co-activators L-Sox5 and Sox6 are also important for cartilage formation [23, 24]. In vitro, transfection of SOX9, 5 and 6 together is sufficient to induce permanent cartilage in mesenchymal precursors . Using TGFβ1 and BMP2 to induce chondrogenic differentiation in hMSCs in a self-gelling alginate scaffold system, we have recently shown that the differentiated cells upregulate the extracellular matrix (ECM) molecules typical of hyaline cartilage, and downregulate molecules typical of MSC functionality . In this system, the SOX TFs are upregulated sufficiently to allow transcription of all essential ECM molecules. We reasoned that a system where availability of TFs is not rate limiting should be suitable to study the relationship between changes in gene expression and concomitant epigenetic changes. To this end we used whole-genome chromatin immunoprecipitation and deep sequencing (ChIP-seq) to quantify H3K4me3, H3K9ac, H327me3, H3K36me3, H3K4me1 and H3K27ac, and reduced representation bisulphite sequencing (RRBS) to quantify DNA methylation in hMSCs before and after 7 days of chondrogenic differentiation. To add to the clinical relevance of our observations, the study is based on primary bone marrow-derived MSCs (BM-MSC) from four donors, which allowed us to also investigate the inter-individual variations in the dynamics of epigenetic marks. We find evidence of two levels of relationship between epigenetic marking and gene expression during chondrogenic differentiation of MSCs: a prepatterning level marked by moderate changes in H3K4me3 and H3K9ac near TSSs and no or little change in gene expression, and a level associated with transcriptional upregulation marked by increased H3K36me3 along the gene body and highly increased H3K4me3 and H3K9ac. The first level is generally associated with the musculoskeletal system and MSC functionality, while the second level is more specifically associated with chondrogenesis. Within the 7-day time frame, changes in promoter DNA methylation do not correlate significantly with changes in gene expression.
Chondrogenic differentiation of hBM-MSCs in 3D scaffold
To first focus our epigenetic analyses on one set of genes typical of differentiated chondrocytes and another set typical of undifferentiated MSCs, we chose 10 canonical chondrogenic (COL2A1, COL9A2, COL9A3, COL11A1, COL11A2, ACAN, COMP, SOX5, SOX6 and SOX9) and mesenchymal (CXCL12, IL6, VCAM1, CCL2, PLAU, CLCF1, ADAMTS1, EFEMP1, ACTA2 and TPM1) signature genes (Figure 1D). The chondrogenic signature genes were selected because collagens II, IX and XI are the molecules required for synthesis of the hyaline cartilage collagen heterofibrils, which make up approximately 40% of the dry weight of hyaline cartilage. Aggrecan and COMP are two other essential components of hyaline ECM, while SOX9, 6 and 5 are essential transcription factors. The MSC signature genes were selected because they encode chemokines, cytokines and ECM degradation molecules important for the role played by MSCs in inflammation (CXCL12, IL6, VCAM1, CCL2, PLAU, CLCF1, ADAMTS1), or represent molecules involved in differentiation along other lineages (ACTA2, TPM1). All these genes are widely distributed throughout the genome, suggesting that if genes respond similarly to the differentiation process, this impact occurs genome-wide. We also included ten control genes representing highly (NUCB1, RPL3 and HSP90B1) and moderately (E2F4, TBP, SDHA) expressed housekeeping genes as well as non-expressed endodermal (GATA5, PDX1) and ectodermal (SOX2, GFAP) genes. The change in mRNA levels for these genes between day 0 and day 7 is shown in Figure 1D, and shows upregulation of the chondrogenic signature genes, downregulation of the MSC signature genes and no change for the control genes. Lists of all up- and downregulated genes in this differentiation system are found in Additional file 1 and Additional file 2.
Histone modifications in promoters and gene bodies
To extend the comparison to include both increases and decreases in histone marks, and to relate these changes to changes in gene expression, heat maps of genome-wide changes in facilitating promoter marks H3K4me3 and H3K9ac and restrictive promoter mark H4K27me3 are shown in Figure 2B, and a heat map of changes in H3K36me3 is shown in Figure 2C, in both cases directly comparable to changes in gene expression. As expected the maps representing the different promoter marks cluster together, with no apparent order between the donors. Upregulated H3K4me3 marks are seen to be closely related to promoters of genes with increased expression, while upregulated H3K9ac marks genes with increased expression, but also unchanged and perhaps even some downregulated genes. Upregulated H3K27me3 clusters with downregulated genes (Figure 2B). For H3K36me3 the relationship between enhanced histone mark and increased gene expression seems even stronger (Figure 2C).
Genome wide correlation between changes in H3K4me3 and changes in gene expression
Based on the association observed between increased H3K4me3 and increased expression for the chondrogenic signature genes, we next investigated if increased H3K4me3 could be found for all genes that were transcriptionally upregulated. Of the 488 genes with increased mRNA levels, 90% were associated with increased H3K4me3 around the TSS, with approximately 65% showing a >2-fold increase (data not shown). We next determined if changes in H3K4me3 were restricted to genes with concomitant changes in mRNA levels. While approximately 1000 genes were significantly either up- or downregulated upon chondrogenic differentiation, >2600 genes showed significant changes in H3K4me3 near the TSS (Additional file 4A). Gene ontology (GO)-term analysis of genes with increased H3K4me3 showed enrichment of functions linked to skeletal system development, ECM and chondrocyte differentiation. The functional terms for the genes with decreasing H3K4me3 were more general and associated with regulation of cytoskeleton, apoptosis, cell development or metabolic processes (Additional file 4B). Genome-wide, only 25% of the genes associated with a >2-fold change in H3K4me3 showed significant change in expression on day 7 (p<0.05) (Additional file 5A). Hardly any of these genes would change in expression over the next two weeks of differentiation, as the vast majority of genes that show change in expression during the first three weeks do so during the first week in this differentiation system. These results indicate that H3K4me3 enrichment occurs during chondrogenic differentiation on genes that are ontologically close to cartilage and musculoskeletal systems, but that these genes are not transcriptionally activated at the ≥2-fold level defined as significant in this study.
Correlation between combination of histone modifications and gene expression
Most of the genes that were transcriptionally unchanged but with highly increased H3K4me3 had no increase in either H3K9ac or H3K36me3 or both. This suggests that lack of H3K9ac and/or H3K36me3 occupancy could explain why transcription of these genes did not change despite their increased H3K4me3, or that changes in H3K9ac and H3K36me3 did not occur because transcription did not change despite changes in H3K4me3 marking. Based on this observation, we determined the relationship between increased gene expression and the presence of all the gene-associated histone modifications, alone or in combination (Additional file 5). Not surprisingly, H3K4me3 and H3K9ac were co-regulated around a large number of genes. However, the histone mark that most closely correlated with increased gene expression was H3K36me3. Of the 488 genes with increased mRNA levels, 94% had the H3K36me3 modification, 85% of these were numerically increased from day 0 to day 7 of differentiation, while 26% were >2-fold increased (Additional file 5 and data not shown). Taken together, this shows that highly increased levels of H3K4me3 and H3K9ac and upregulated levels of H3K36me3 are found near practically every transcriptionally upregulated gene in this differentiation system.
The impact of cis-acting regulatory elements
In order to identify potential chondrogenic enhancer regions, we investigated changes in these regulatory marks within 50 kb of the TSS of genes that changed transcriptionally during chondrogenic differentiation. For our signature genes, we found at least one site of >2 fold enrichment in both histone modifications near every chondrogenic signature gene, and at least one site with a significant >2 fold decrease of these marks near every MSC signature gene (Figure 6B). None of the unchanged control genes showed changes in H3K4me1 or H3K27ac (data not shown). The changes in H3K4me1 and H3K27ac near the chondrogenic marker gene COL2A1 are shown in Figure 6C. We detected six areas with enrichment in both H3K4me1 and H3K27ac in the chondrogenically differentiated cells. Two of these could be identified as described and validated enhancer regions with binding sites for SOX9 [30–33]. The other four regions (Figure 6C, marked in red) have not yet been described, and are likely to be new regulatory elements of importance for the regulation of COL2A1 expression.
Transcription factor binding sites in cis-acting regulatory regions - Motif search
The impact of changes in cis-acting regulatory elements on the regulation of gene expression most likely occurs through the binding of TFs. In order to identify TFs which may be involved in the regulation of chondrogenesis we performed a motif search in all areas with increase of H3K27ac during the differentiation process, using JASPAR, Transfac and Uniprobe databases. We found enrichment in motifs which could bind TFs belonging to several families (Figure 6D). The SRY family is especially important in the context of chondrogenesis, as it contains SOX5, 6 and 9, which are essential and sufficient to start the differentiation process [23, 34]. The literature also contains evidence of involvement in chondrogenesis of members of all the other families with enriched binding motifs in these regions [35–38]. Interestingly, for most of these TF families binding motifs could be found also within the new regulatory areas close to the COL2A1 gene, further supporting their potential role in chondrogenic differentiation of MSCs (Figure 6D). Most significantly, new regions 1, 2 and 4 could be shown to contain 3, 8 and 2 sequences identical to the SOX9 binding motif (AACAAT), respectively, while new region 3 contains no such sequences.
To determine if there is a relationship between DNA methylation and changes in gene expression before and after 7 days of chondrogenic differentiation, genome-wide DNA methylation in CpG-rich regions was mapped using RRBS. More than 2 × 106 unique CpGs were observed with at least 10X coverage in most samples. The overall methylation patterns were highly similar between samples and followed expected patterns. CpGs in CpG-rich regions (including CpG islands) tended to be hypomethylated (<10% median methylation across covered CpGs) while CpGs in CpG-poor regions tended to be hypermethylated (>50% median methylation; Additional file 10A). This relationship was preserved when the promoter regions were examined specifically. A highly significant correlation was found between the average methylation values for the MSC and the average methylation values for the chondrocytes in promoter regions (r2 = 0.78), with the methylation states of the vast majority of promoters remaining unchanged . For the MSC and chondrocyte signature genes, promoters were hypomethylated both before and after chondrogenesis, which suggests that DNA methylation is unlikely to play a major role in regulating these key genes during chondrogenesis. Genome-wide, only 37 promoters (0.14% of those examined), including those of 9 microRNAs, showed a mean increase of at least 50% in methylation across the donors, and only 284 (1.1%) showed more than a 25% increase. However, except for GAS6 (growth-arrest specific 6, five-fold downregulated, increased promoter methylation), these genes were generally expressed at low levels both before and after differentiation and no correlation to changes in the investigated histone marks could be detected. Moreover, at least 15 of the 37 genes that gained >50% methylation across the donors were already hypermethylated in one of the donors prior to differentiation, which suggests that the observed promoter methylation events are not likely to be directly linked to regulation of chondrogenesis.
Chondrogenic differentiation opens up a new epigenetic landscape in MSCs
In vitro tissue engineering of hyaline cartilage implants using MSCs and a biomaterial may well represent the best future treatment option for patients with lesions of articular cartilage. To be applicable to a wide range of patients the tissue engineering strategy should be robust, with little between-donor variability. Preferably, the differentiation process should be induced synchronously in all the cells used. A priori, this may be difficult to achieve: MSCs from different donors behave differently in cell expansion and differentiation cultures, and MSCs are known to be heterogenous both phenotypically and functionally [39, 40]. In contrast, using the alginate disc differentiation procedure described here, hyaline ECM molecules are produced by the vast majority of cells in all donors tested so far  (Küchler AM and Brinchmann JE, unpublished results). However, there is a delay between upregulation of mRNA and detection of ECM molecules, and ECM molecule synthesis does not seem to be synchronous. Against this background, we observe a remarkable similarity between the histone modifications induced after seven days of chondrogenic differentiation using bone marrow-derived MSCs from four donors. At this time, a new epigenetic landscape is opened up in these cells. Changes in the H3K4me3 and H3K9ac modifications occur near the TSS in thousands of genes, many of which are ontogenetically linked to skeletal system development, ECM, chondrocyte differentiation and MSC functionality. Although these modifications are known to be associated with initiation of transcription [4, 41], no significant change in mRNA levels is actually observed in the majority of these genes. With the level of change in gene expression set as significant in this study (≥2-fold) it is possible that a small change in gene expression may accompany these changes in histone marks, or perhaps a greater change in a very small subset of cells. However, these changes are most reminiscent of the prepatterning of developmental gene expression by H3K4me3 recently observed in zebrafish before zygotic genome activation . It is possible that every chondrogenic differentiation protocol would induce an epigenetic state that maintains activation potential of genes required both for the desired and also for related differentiation pathways . It is also possible that certain components of the present differentiation cocktail, e.g. BMP2, may improve chondrogenic differentiation (Jakobsen RB, Østrup E, Mikkelsen TS and Brinchmann JE, unpublished observations) while at the same time setting an epigenetic stage also for bone differentiation.
Close correlation between histone marks and gene expression
Although many genes within the general domain of the musculoskeletal system show increased H3K4me3 following 7 days of chondrogenic differentiation, only the genes with the highest increase in H3K4me3 could be correlated with increased gene expression. The distribution of H3K4me3 and H3K9ac was, in many instances, very similar, suggesting co-occurrence on the same histone tails. H3K4me3 and H3K9ac may poise genes for transcription without transcription actually occurring. However, increased H3K4me3 and H3K9ac in promoter regions is, on its own, unlikely to explain the close correlation between highly increased H3K4me3 and H3K9ac and increased transcriptional activity. Increase in the distribution of these modifications within the gene body may conceivably explain the increased rate of transcription, as H3K4me3 has also been suggested to stimulate transcriptional elongation . However, the modification most closely correlated with increased transcription is H3K36me3. The combination of H3K4me3 with H3K36me3 is observed at practically every gene with increased mRNA level, and at few other genes. The following model may reconcile our results with current knowledge about the relationship between histone modifications and transcription. Trimethylation of H3K4 is regulated by several factors, including trithorax group proteins and cyclin-dependent kinase-9 (CDK9) of the positive transcription elongation factor-b complex (P-TEFb) [43, 44]. The regulation of H3K9ac is not well known, but its co-occurrence with H3K4me3 has been described [7, 45], and the two modifications may have common chromatin readers [46, 47]. These modifications frequently co-occur with RNA polymerase II . However, transcriptional elongation does not occur in the absence of trimethylated H3K36 along the gene body. Methylation of H3K36 is performed by the methyltransferase Set2 and is, like H3K4 trimethylation, dependent on CDK9 [44, 48]. CDK9 and Set2 are part of a large complex of proteins which is involved in the release of paused polymerase II . With the factors responsible for trimethylation of H3K36 being part of the transcription elongation complex, it is uncertain at this time whether H3K36me3 is a prerequisite for, or a consequence of the passage of this complex along the gene body. In either case transcriptional elongation, heralded by H3K36me3, is the most specific correlate of actively transcribed genes in the present data set. Further research will, hopefully, explain how transcriptional elongation is restricted to a set of genes with particular relevance for hyaline chondrogenesis in this differentiation model.
Apparent co-occurrence of active and repressive marks
Genes regulated both by H3K4me3 and H3K27me3 in the promoter region are called bivalent genes. ESCs are known to have a high proportion of bivalent genes, particularly genes responsible for specification along lineages [4, 5, 8]. These genes are not transcribed as long as differentiation is inhibited. On a cell population basis, we also observed what looks like bivalent genes, both among the upregulated chondrogenic signature genes and among the downregulated MSC signature genes. On day 7 of differentiation, transcription of all these genes was observed. We cannot rule out the possibility that there are genes among these that are actively transcribed in the presence of both H3K4me3 and H3K27me3 near the TSS. More likely, however, is an explanation based on non-synchronous changes in histone modifications and gene transcription between the different cells in a heterogeneous population. This would imply, for COL2A1, for instance, that one subset of MSCs lose H3K27me3 and gain H3K4me3 (and H3K36me3) at day 7. This subset is responsible for the observed increased level of COL2A1 mRNA, while another subset with retained H3K27me3 near the COL2A1 promoter does not support transcription. The level of COL2A1 mRNA increases further through days 14 and 21, which is consistent with subsequent loss of H3K27me3 and gain of active marks in cells that are still repressed at day 7.
Relationship between distal cis-regulatory elements and gene expression
The genome-wide analysis of all 50 kb areas around TSS revealed a correlation between gene expression and changes in the regulatory histone marks H3K4me1 or H3K27ac. For the signature genes, every upregulated chondrogenic gene is associated with at least one upregulated enhancer region and every downregulated MSC-specific gene is associated with at least one downregulated enhancer region, whereas none of the control genes show any such change within the 50 kb region. For the COL2A1 marker gene, 4 new cis-regulatory regions were discovered, all containing binding motifs for TFs essential for chondrogenesis. This strongly suggests that distal cis-regulatory regions impact importantly on the regulation of gene expression in this differentiation model. However, most of the histone modifications marking cis-regulatory regions occurred >50 kb away from genes with changed expression, suggesting that the enhancers may be far removed from the target gene.
We have performed genome-wide quantitative epigenetic analysis on biological replicates of primary human MSCs undergoing chondrogenic differentiation in alginate disc scaffolds relevant for tissue engineering of hyaline cartilage. The differentiation cocktail and 3D culture induce changes in H3K4me3 and H3K9ac in genes associated with the musculoskeletal system, but mostly in the absence of changes in gene expression. Increased gene expression may be predicted by particularly high increase in H3K4me3 and by increase in H3K36me3. This implicates transcriptional elongation as the specificity conferring, rate limiting step in the differentiation process. Cis-regulatory activity marked by increasing H3K27ac and H3K4me1 is involved in transcriptional regulation in the vast majority of genes with altered expression. Within the 7-day time frame, changes in promoter DNA methylation do not correlate significantly with changes in gene expression.
Culture and validation of MSCs
The study was approved by the Regional Committee for Medical Research Ethics, Southern Norway. Bone marrow was obtained from four healthy voluntary donors following informed, written consent. The MSCs were isolated and cultured as described in . The cells were validated as MSCs by expression of CD105, CD73, CD90 and CD44, and failure to express CD34, HLA- DR, CD45, CD14 and CD19. The multipotency of the MSCs was validated by differentiation to osteogenic and adipogenic lineages as previously described in .
In vitro chondrogenic differentiation
The chondrogenic differentiation system was described previously . Briefly, MSCs were established in a self-gelling system containing PRONOVA-LVG alginate and calcium-alginate-particles to a final alginate concentration of 1% and a final cell concentration of 5 × 106 cells/ml. Chondrogenic differentiation was induced by 500 ng/ml BMP2, 10 ng/ml TGFB1 and 0.1 μM dexamethasone in serum free medium containing supplements as described .
Validation of extracellular matrix synthesis
Fluorescence immunohistochemistry was performed as described  with reagents as specified in Additional file 11. The supernatants from disc cultures were analyzed for sulfated glycosaminoglycan content in technical quadruplicates for each donor using the Blyscan™ kit (Biocolor Ltd).
RNA preparation and expression analysis
Total RNA was prepared using TRIzol (Invitrogen). Depolymerisation of the alginate discs to obtain single cells for RNA isolation was performed by enzymatic digestion using G-Lyase (kindly provided by professor Gudmund Skjåk-Bræk) as described in . RT-qPCR was performed using TaqMan arrays (Applied Biosystems) and normalized to GAPDH expression. Results are presented as mean plus standard deviation of technical triplicates. Microarray assays were performed in uniplicates per donor and differentiation state using the Human-6 v3 Expression Beadchips (Illumina), and the analysis was performed using J-express 2009 .
Tracks were uploaded in a UCSC Genome Browser, using unique reads extended to 200 bp. Experiments were normalized by enforcing equal read counts. 500 bp bins were generated, and the coverage in each bin was calculated for both chip and control. For each bin, a fold change value and a poisson p-value was generated. P-values were FDR corrected. The read counts from all experiments were quantile normalized. Bins with a log2 fold change >0.2 or <−0.2 in at least three of the four replicates were defined as having an increased or decreased signal, respectively. For quantification of H3K27me3, H3K4me3, and H3K9ac the sum of aligned reads within 2 kb on both sides of the TSS for each gene was calculated. This was performed for all biological replicates for both MSCs and differentiated cells, and a mean for each cell population was calculated. For H3K36me3, with known functionality throughout the entire gene , the longest transcript for each Ensembl gene was selected. The total number of aligned reads within this transcript was then calculated for all biological replicates in both cell populations, and a mean for each population was calculated. For the profiling of H3K4me1 and H3K27ac signals the highest ratio of bins associated to a single gene was calculated.
For promoter marks, Ensembl transcripts were used to identify transcripts with a 5′ end within 2 kb of bins. The transcript with the closest 5′ end was selected. For enhancer marks the association was made to the most differentially expressed gene within 50 kb. Microarray expression data were integrated using refseq IDs.
Motif enrichment analysis
For each histone mark, the top 5000 increasing and decreasing bins were selected. Bins next to each other were merged. The regions were scanned with motifs from Jaspar, Transfac and Uniprobe with a score threshold of 80%. Fisher test was used to compare the number of sequences with at least one hit and the number of sequences with no hits between the two sets. Wilcox test was used to compare the rank of total number of hits between the two sets. Motifs with a p-value < 0.05 from both tests were defined as significant.
RRBS was performed as previously described . To analyze potential changes in promoter DNA methylation during chondrogenesis we examined 26,003 distinct TSSs annotated on the UCSC Genome Browser (derived from the gh19 refFlat database). For each sample, we assigned a methylation value to each TSS for which we had at least 5 distinct CpGs with at least 5X coverage each within −1,000 bp to +1,500 bp relative to the TSS. The methylation value for each such TSS was taken to be the median methylation value of all of the CpGs within this region with at least 5X coverage, using at least three of the four donors in each state.
Sequences and data access
All ChIP-Seq and RRBS data have been deposited in public databases and can be accessed from the NIH Roadmap Project on Epigenomics website: http://www.roadmapepigenomics.org/ and in the NCBI GEO database under accession GSE19465.
Reprint author: Jan E. Brinchmann MD, PhD
Institute of Immunology, Oslo University Hospital Rikshospitalet, PO Box 4950 Nydalen, 0424 Oslo, Norway
Telephone: +47 22 84 04 89
Fax: +47 22 85 10 58
Mesenchymal stem cell
Human bone marrow-derived mesenchymal stem cells
Embryonal stem cell
Transcriptional start site
Chromatin immunoprecipitation and deep sequencing
Reduced representation bisulphite sequencing
Positive transcription elongation factor-b complex.
This work was supported by a grant from South-Eastern Norway Regional Health Authority, Storforsk and Stamceller grants from the Research Council of Norway and Gidske og Peter Jacob Sørensens Foundation for the Promotion of Science. Data generation was supported by the NIH Roadmap Epigenomics Program (grant U01ES017155 to Bradley E. Bernstein) and performed by the NIH Reference Epigenome Mapping Center at the Broad Institute. We thank Leif Lindeman for introduction to the ChIP assay and Axel Küchler for performing the immunohistochemical stainings shown in Figure 1B. We acknowledge Professor Gudmund Skjåk-Bræk, Norwegian University of Science and Technology, Trondheim, Norway for providing us with G-Lyase. We thank Charles B. Epstein, Tim Fennell and Noam Shoresh, Broad Institute, Cambridge, MA, for epigenomics pipeline support.
- Ahmed TA, Hincke MT: Strategies for articular cartilage lesion repair and functional restoration. Tissue Eng. 2010, 16 (3): 305-329. 10.1089/ten.teb.2009.0590.View ArticleGoogle Scholar
- Vinatier C, Bouffi C, Merceron C, Gordeladze J, Brondello JM, Jorgensen C, Weiss P, Guicheux J, Noel D: Cartilage tissue engineering: towards a biomaterial-assisted mesenchymal stem cell therapy. Curr Stem Cell Res Ther. 2009, 4 (4): 318-329. 10.2174/157488809789649205.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernstein BE, Kamal M, Lindblad-Toh K, Bekiranov S, Bailey DK, Huebert DJ, McMahon S, Karlsson EK, Kulbokas EJ, Gingeras TR, Schreiber SL, Lander ES: Genomic maps and comparative analysis of histone modifications in human and mouse. Cell. 2005, 120 (2): 169-181. 10.1016/j.cell.2005.01.001.View ArticlePubMedGoogle Scholar
- Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, Lee W, Mendenhall E, O’Donovan A, Presser A, Russ C, Xie X, Meissner A, Wernig M, Jaenisch R, Nusbaum C, Lander ES, Bernstein BE: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007, 448 (7153): 553-560. 10.1038/nature06008.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawkins RD, Hon GC, Lee LK, Ngo Q, Lister R, Pelizzola M, Edsall LE, Kuan S, Luu Y, Klugman S, Antosiewicz-Bourget J, Ye Z, Espinoza C, Agarwahl S, Shen L, Ruotti V, Wang W, Stewart R, Thomson JA, Ecker JR, Ren B: Distinct epigenomic landscapes of pluripotent and lineage-committed human cells. Cell Stem Cell. 2010, 6 (5): 479-491. 10.1016/j.stem.2010.03.018.PubMed CentralView ArticlePubMedGoogle Scholar
- Challen GA, Sun D, Jeong M, Luo M, Jelinek J, Berg JS, Bock C, Vasanthakumar A, Gu H, Xi Y, Liang S, Lu Y, Darlington GJ, Meissner A, Issa JP, Godley LA, Li W, Goodell MA: Dnmt3a is essential for hematopoietic stem cell differentiation. Nat Genet. 2011, 44 (1): 23-31. 10.1038/ng.1009.PubMed CentralView ArticlePubMedGoogle Scholar
- Guenther MG, Levine SS, Boyer LA, Jaenisch R, Young RA: A chromatin landmark and transcription initiation at most promoters in human cells. Cell. 2007, 130 (1): 77-88. 10.1016/j.cell.2007.05.042.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, Jaenisch R, Wagschal A, Feil R, Schreiber SL, Lander ES: A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006, 125 (2): 315-326. 10.1016/j.cell.2006.02.041.View ArticlePubMedGoogle Scholar
- Suganuma T, Workman JL: Signals and combinatorial functions of histone modifications. Annu Rev Biochem. 2011, 80: 473-499. 10.1146/annurev-biochem-061809-175347.View ArticlePubMedGoogle Scholar
- Mikkelsen TS, Xu Z, Zhang X, Wang L, Gimble JM, Lander ES, Rosen ED: Comparative epigenomic analysis of murine and human adipogenesis. Cell. 2010, 143 (1): 156-169. 10.1016/j.cell.2010.09.006.PubMed CentralView ArticlePubMedGoogle Scholar
- Sorensen AL, Jacobsen BM, Reiner AH, Andersen IS, Collas P: Promoter DNA methylation patterns of differentiated cells are largely programmed at the progenitor stage. Mol Biol Cell. 2010, 21 (12): 2066-2077. 10.1091/mbc.E10-01-0018.PubMed CentralView ArticlePubMedGoogle Scholar
- Asp P, Blum R, Vethantham V, Parisi F, Micsinai M, Cheng J, Bowman C, Kluger Y, Dynlacht BD: Genome-wide remodeling of the epigenetic landscape during myogenic differentiation. Proc Natl Acad Sci USA. 2011, 108 (22): E149-158. 10.1073/pnas.1102223108.PubMed CentralView ArticlePubMedGoogle Scholar
- Lien WH, Guo X, Polak L, Lawton LN, Young RA, Zheng D, Fuchs E: Genome-wide maps of histone modifications unwind in vivo chromatin states of the hair follicle lineage. Cell Stem Cell. 2011, 9 (3): 219-232. 10.1016/j.stem.2011.07.015.PubMed CentralView ArticlePubMedGoogle Scholar
- Cedar H, Bergman Y: Epigenetics of haematopoietic cell development. Nat Rev Immunol. 2011, 11 (7): 478-488. 10.1038/nri2991.View ArticlePubMedGoogle Scholar
- Cui K, Zang C, Roh TY, Schones DE, Childs RW, Peng W, Zhao K: Chromatin signatures in multipotent human hematopoietic stem cells indicate the fate of bivalent genes during differentiation. Cell Stem Cell. 2009, 4 (1): 80-93. 10.1016/j.stem.2008.11.011.PubMed CentralView ArticlePubMedGoogle Scholar
- Hodges E, Molaro A, Dos Santos CO, Thekkat P, Song Q, Uren PJ, Park J, Butler J, Rafii S, McCombie WR, Smith AD, Hannon GJ: Directional DNA methylation changes and complex intermediate states accompany lineage specificity in the adult hematopoietic compartment. Mol Cell. 2011, 44 (1): 17-28. 10.1016/j.molcel.2011.08.026.PubMed CentralView ArticlePubMedGoogle Scholar
- Ezura Y, Sekiya I, Koga H, Muneta T, Noda M: Methylation status of CpG islands in the promoter regions of signature genes during chondrogenesis of human synovium-derived mesenchymal stem cells. Arthritis Rheum. 2009, 60 (5): 1416-1426. 10.1002/art.24472.View ArticlePubMedGoogle Scholar
- Furumatsu T, Asahara H: Histone acetylation influences the activity of Sox9-related transcriptional complex. Acta Med Okayama. 2011, 64 (6): 351-357.Google Scholar
- Zimmermann P, Boeuf S, Dickhut A, Boehmer S, Olek S, Richter W: Correlation of COL10A1 induction during chondrogenesis of mesenchymal stem cells with demethylation of two CpG sites in the COL10A1 promoter. Arthritis Rheum. 2008, 58 (9): 2743-2753. 10.1002/art.23736.View ArticlePubMedGoogle Scholar
- Goldring MB, Tsuchimochi K, Ijiri K: The control of chondrogenesis. J Cell Biochem. 2006, 97 (1): 33-44. 10.1002/jcb.20652.View ArticlePubMedGoogle Scholar
- Ichinose S, Yamagata K, Sekiya I, Muneta T, Tagami M: Detailed examination of cartilage formation and endochondral ossification using human mesenchymal stem cells. Clin Exp Pharmacol Physiol. 2005, 32 (7): 561-570. 10.1111/j.1440-1681.2005.04231.x.View ArticlePubMedGoogle Scholar
- Lefebvre V, Smits P: Transcriptional control of chondrocyte fate and differentiation. Birth Defects Res C Embryo Today. 2005, 75 (3): 200-212. 10.1002/bdrc.20048.View ArticlePubMedGoogle Scholar
- Lefebvre V, Behringer RR, de Crombrugghe B: L-Sox5, Sox6 and Sox9 control essential steps of the chondrocyte differentiation pathway. Osteoarthritis and cartilage / OARS, Osteoarthritis Research Society. 2001, 9 (Suppl A): S69-75.View ArticlePubMedGoogle Scholar
- Smits P, Li P, Mandel J, Zhang Z, Deng JM, Behringer RR, de Crombrugghe B, Lefebvre V: The transcription factors L-Sox5 and Sox6 are essential for cartilage formation. Dev Cell. 2001, 1 (2): 277-290. 10.1016/S1534-5807(01)00003-X.View ArticlePubMedGoogle Scholar
- Ikeda T, Zhang J, Chano T, Mabuchi A, Fukuda A, Kawaguchi H, Nakamura K, Ikegawa S: Identification and characterization of the human long form of Sox5 (L-SOX5) gene. Gene. 2002, 298 (1): 59-68. 10.1016/S0378-1119(02)00927-7.View ArticlePubMedGoogle Scholar
- Herlofsen SR, Kuchler AM, Melvik JE, Brinchmann JE: Chondrogenic differentiation of human bone marrow-derived mesenchymal stem cells in self-gelling alginate discs reveals novel chondrogenic signature gene clusters. Tissue Eng Part A. 2011, 17 (7–8): 1003-1013.PubMed CentralView ArticlePubMedGoogle Scholar
- Heintzman ND, Hon GC, Hawkins RD, Kheradpour P, Stark A, Harp LF, Ye Z, Lee LK, Stuart RK, Ching CW, Ching KA, Antosiewicz-Bourget JE, Liu H, Zhang X, Green RD, Lobanenkov VV, Stewart R, Thomson JA, Crawford GE, Kellis M, Ren B: Histone modifications at human enhancers reflect global cell-type-specific gene expression. Nature. 2009, 459 (7243): 108-112. 10.1038/nature07829.PubMed CentralView ArticlePubMedGoogle Scholar
- Heintzman ND, Stuart RK, Hon G, Fu Y, Ching CW, Hawkins RD, Barrera LO, Van Calcar S, Qu C, Ching KA, Wang W, Weng Z, Green RD, Crawford GE, Ren B: Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome. Nat Genet. 2007, 39 (3): 311-318. 10.1038/ng1966.View ArticlePubMedGoogle Scholar
- Tew SR, Clegg PD, Brew CJ, Redmond CM, Hardingham TE: SOX9 transduction of a human chondrocytic cell line identifies novel genes regulated in primary human chondrocytes and in osteoarthritis. Arthritis Res Ther. 2007, 9 (5): R107-10.1186/ar2311.PubMed CentralView ArticlePubMedGoogle Scholar
- Kerney R, Hall BK, Hanken J: Regulatory elements of Xenopus col2a1 drive cartilaginous gene expression in transgenic frogs. Int J Dev Biol. 2010, 54 (1): 141-150. 10.1387/ijdb.092848rk.View ArticlePubMedGoogle Scholar
- Lefebvre V, Huang W, Harley VR, Goodfellow PN, de Crombrugghe B: SOX9 is a potent activator of the chondrocyte-specific enhancer of the pro alpha1(II) collagen gene. Mol Cell Biol. 1997, 17 (4): 2336-2346.PubMed CentralView ArticlePubMedGoogle Scholar
- Oh CD, Maity SN, Lu JF, Zhang J, Liang S, Coustry F, de Crombrugghe B, Yasuda H: Identification of SOX9 interaction sites in the genome of chondrocytes. PLoS One. 2010, 5 (4): e10113-10.1371/journal.pone.0010113.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou G, Lefebvre V, Zhang Z, Eberspaecher H, de Crombrugghe B: Three high mobility group-like sequences within a 48-base pair enhancer of the Col2a1 gene are required for cartilage-specific expression in vivo. J Biol Chem. 1998, 273 (24): 14989-14997. 10.1074/jbc.273.24.14989.View ArticlePubMedGoogle Scholar
- Park JS, Yang HN, Woo DG, Jeon SY, Do HJ, Lim HY, Kim JH, Park KH: Chondrogenesis of human mesenchymal stem cells mediated by the combination of SOX trio SOX5, 6, and 9 genes complexed with PEI-modified PLGA nanoparticles. Biomaterials. 2011, 32 (14): 3679-3688. 10.1016/j.biomaterials.2011.01.063.View ArticlePubMedGoogle Scholar
- Amarilio R, Viukov SV, Sharir A, Eshkar-Oren I, Johnson RS, Zelzer E: HIF1alpha regulation of Sox9 is necessary to maintain differentiation of hypoxic prechondrogenic cells during early skeletogenesis. Development. 2007, 134 (21): 3917-3928. 10.1242/dev.008441.View ArticlePubMedGoogle Scholar
- Greenwel P, Inagaki Y, Hu W, Walsh M, Ramirez F: Sp1 is required for the early response of alpha2(I) collagen to transforming growth factor-beta1. J Biol Chem. 1997, 272 (32): 19738-19745. 10.1074/jbc.272.32.19738.View ArticlePubMedGoogle Scholar
- Iwamoto M, Higuchi Y, Enomoto-Iwamoto M, Kurisu K, Koyama E, Yeh H, Rosenbloom J, Pacifici M: The role of ERG (ets related gene) in cartilage development. Osteoarthritis and cartilage / OARS, Osteoarthritis Research Society. 2001, 9 (Suppl A): S41-47.View ArticlePubMedGoogle Scholar
- Nakada C, Iida A, Tabata Y, Watanabe S: Forkhead transcription factor foxe1 regulates chondrogenesis in zebrafish. J Exp Zool. 2009, 312 (8): 827-840.View ArticleGoogle Scholar
- Boquest AC, Shahdadfar A, Fronsdal K, Sigurjonsson O, Tunheim SH, Collas P, Brinchmann JE: Isolation and transcription profiling of purified uncultured human stromal stem cells: alteration of gene expression after in vitro cell culture. Mol Biol Cell. 2005, 16 (3): 1131-1141. 10.1091/mbc.E04-10-0949.PubMed CentralView ArticlePubMedGoogle Scholar
- Nombela-Arrieta C, Ritz J, Silberstein LE: The elusive nature and function of mesenchymal stem cells. Nature reviews. 2011, 12 (2): 126-131. 10.1038/nrm3049.PubMed CentralView ArticlePubMedGoogle Scholar
- Weake VM, Workman JL: Inducible gene expression: diverse regulatory mechanisms. Nat Rev Genet. 2010, 11 (6): 426-437. 10.1038/nrg2781.View ArticlePubMedGoogle Scholar
- Lindeman LC, Andersen IS, Reiner AH, Li N, Aanes H, Ostrup O, Winata C, Mathavan S, Muller F, Alestrom P, Collas P: Prepatterning of developmental gene expression by modified histones before zygotic genome activation. Dev Cell. 2011, 21 (6): 993-1004. 10.1016/j.devcel.2011.10.008.View ArticlePubMedGoogle Scholar
- Schuettengruber B, Chourrout D, Vervoort M, Leblanc B, Cavalli G: Genome regulation by polycomb and trithorax proteins. Cell. 2007, 128 (4): 735-745. 10.1016/j.cell.2007.02.009.View ArticlePubMedGoogle Scholar
- Pirngruber J, Shchebet A, Johnsen SA: Insights into the function of the human P-TEFb component CDK9 in the regulation of chromatin modifications and co-transcriptional mRNA processing. Cell cycle (Georgetown, Tex). 2009, 8 (22): 3636-3642. 10.4161/cc.8.22.9890.View ArticleGoogle Scholar
- Vermeulen M, Mulder KW, Denissov S, Pijnappel WW, van Schaik FM, Varier RA, Baltissen MP, Stunnenberg HG, Mann M, Timmers HT: Selective anchoring of TFIID to nucleosomes by trimethylation of histone H3 lysine 4. Cell. 2007, 131 (1): 58-69. 10.1016/j.cell.2007.08.016.View ArticlePubMedGoogle Scholar
- Vermeulen M, Eberl HC, Matarese F, Marks H, Denissov S, Butter F, Lee KK, Olsen JV, Hyman AA, Stunnenberg HG, Mann M: Quantitative interaction proteomics and genome-wide profiling of epigenetic histone marks and their readers. Cell. 2010, 142 (6): 967-980. 10.1016/j.cell.2010.08.020.View ArticlePubMedGoogle Scholar
- Yun M, Wu J, Workman JL, Li B: Readers of histone modifications. Cell Res. 2011, 21 (4): 564-578. 10.1038/cr.2011.42.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoh SM, Lucas JS, Jones KA: The Iws1:Spt6:CTD complex controls cotranscriptional mRNA biosynthesis and HYPB/Setd2-mediated histone H3K36 methylation. Genes Dev. 2008, 22 (24): 3422-3434. 10.1101/gad.1720008.PubMed CentralView ArticlePubMedGoogle Scholar
- Dysvik B, Jonassen I: J-Express: exploring gene expression data using Java. Bioinformatics (Oxford, England). 2001, 17 (4): 369-370. 10.1093/bioinformatics/17.4.369.View ArticleGoogle Scholar
- Kolasinska-Zwierz P, Down T, Latorre I, Liu T, Liu XS, Ahringer J: Differential chromatin marking of introns and expressed exons by H3K36me3. Nat Genet. 2009, 41 (3): 376-381. 10.1038/ng.322.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu H, Bock C, Mikkelsen TS, Jager N, Smith ZD, Tomazou E, Gnirke A, Lander ES, Meissner A: Genome-scale DNA methylation mapping of clinical samples at single-nucleotide resolution. Nat Methods. 2010, 7 (2): 133-136. 10.1038/nmeth.1414.PubMed CentralView ArticlePubMedGoogle Scholar
- Murdoch AD, Grady LM, Ablett MP, Katopodi T, Meadows RS, Hardingham TE: Chondrogenic differentiation of human bone marrow stem cells in transwell cultures: generation of scaffold-free cartilage. Stem Cells. 2007, 25 (11): 2786-2796. 10.1634/stemcells.2007-0374.View ArticlePubMedGoogle Scholar