Summary
Tissue generation and repair requires a stepwise process of cell fate restriction to ensure that adult stem cells differentiate in a timely and appropriate manner. A crucial role has been implicated for Polycomb-group (PcG) proteins and the H3K27me3 repressive histone mark in coordinating the transcriptional programmes necessary for this process, but the targets and developmental timing for this repression remain unclear. To address these questions, we generated novel genome-wide maps of H3K27me3 and H3K4me3 in freshly isolated muscle stem cells. These data, together with the analysis of two conditional Ezh2-null mouse strains, identified a critical proliferation phase in which Ezh2 activity is essential. Mice lacking Ezh2 in satellite cells exhibited decreased muscle growth, severely impaired regeneration and reduced stem cell number, due to a profound failure of the proliferative progenitor population to expand. Surprisingly, deletion of Ezh2 after the onset of terminal differentiation did not impede muscle repair or homeostasis. Using these knockout models and the RNA-Seq and ChIP-Seq datasets, we show that Ezh2 does not regulate the muscle differentiation process in vivo. These results emphasise the lineage and cell-type-specific functions of Ezh2 and Polycomb repressive complex 2.
Introduction
Maintenance of a correct differentiation programme is essential for tissue homeostasis and repair. Accumulating evidence suggests that the Polycomb group (PcG) proteins, which mediate transcriptionally repressive epigenetic information, represent likely candidates for the regulation of adult stem cell fate and differentiation. Polycomb repressive complex 2 (PRC2) is composed of four core subunits, Ezh2, Suz12, EED and RbAp48, with Ezh2 responsible for the trimethylation of lysine 27 of histone H3 (H3K27me3) (Margueron et al., 2008; Shen et al., 2008). Numerous studies investigating embryonic stem cells (ESCs) have demonstrated that PRC2 and H3K27me3 transcriptionally repress key factors involved in lineage specification (which are de-repressed in PRC2-deficient ESCs) leading to aberrant differentiation programs (Boyer et al., 2006; Chamberlain et al., 2008; Leeb et al., 2010; Mikkelsen et al., 2007; Pan et al., 2007; Pasini et al., 2007; Pasini et al., 2010; Shen et al., 2008; Wang et al., 2008). In ESCs, genes encoding key developmental regulators have been suggested to contain multiple histone modifications, i.e. the repressive H3K27me3 mark and a histone mark associated with gene transcription, H3K4me3. These ‘bivalent’ domains are thought to poise these factors for rapid expression or repression (Bernstein et al., 2006; Mikkelsen et al., 2007; Vastenhouw and Schier, 2012).
The ability to study the function of PRC2 within an adult stem cell population has been hampered by the prenatal lethality of Eed, Suz12 and Ezh2 null embryos (Faust et al., 1995; O'Carroll et al., 2001; Pasini et al., 2004; Schumacher et al., 1996). Conditional deletion of Ezh2, for example, in B cells (Su et al., 2003), neuronal stem cells (Hirabayashi et al., 2009; Pereira et al., 2010), epidermal progenitors (Ezhkova et al., 2011; Ezhkova et al., 2009; Lien et al., 2011) and pancreatic islet β cells (Chen et al., 2009) has provided insights into the function of PRC2 within lineage restricted somatic cells. However, the differing phenotypes within these studies indicate that Ezh2 action may be highly dependent on the differentiation stage of the stem cell and may even regulate the differentiation programme (Benoit et al., 2012). Interestingly, loss of Ezh2 within cardiac progenitors leads to an upregulation of Six1 dependent skeletal muscle genes (Delgado-Olguín et al., 2012), suggesting that Ezh2 may not only be required for stem cell function but may also act to ensure the correct lineage choice of a stem cell.
The process of skeletal muscle development and repair serves as an excellent paradigm with which to investigate adult stem cell self-renewal, lineage specification and terminal differentiation. Skeletal muscle has a high regenerative capacity, with satellite cells (the principal muscle stem cell) being responsible for this capability (Lepper et al., 2011; Murphy et al., 2011; Sacco et al., 2008; Sambasivan et al., 2011). The mechanisms that govern the lineage restricted nature of these stem cells, and the fate choice to either differentiate or self-renew, have until recently been studied in the C2C12 myoblast cell line, which is thought to recapitulate the differentiation programmes that occur in vivo. Indeed, using the C2C12 cell line, Ezh2 has been suggested to be a master regulator of the muscle terminal differentiation process by suppressing key muscle differentiation genes (Juan et al., 2009; Palacios et al., 2010; Seenundun et al., 2010). Recently, Juan et al. (Juan et al., 2011) suggested that loss of Ezh2 in Pax7+ satellite cells causes defects in stem cell self-renewal, although the role of Ezh2 within the terminal differentiation process was not addressed in this study. Surprisingly, expression profiling from the Ezh2 null skeletal muscle tissue (a mixed population of differentiated muscle and non-myogenic cells) showed only small changes in gene expression patterns (123 upregulated transcripts, with the majority only two-fold higher than controls) (Juan et al., 2011). However, studies utilising the C2C12 cell line have shown that over 4,000 gene promoters are marked by H3K27me3 (Asp et al., 2011). Indeed, siRNA-mediated depletion of Ezh2 or Suz12 in C2C12 cells leads to enhanced muscle differentiation and upregulation of muscle differentiation genes (Asp et al., 2011; Juan et al., 2009).
In this study, we use knockout models and genome wide datasets to define the myogenic cell population in which Ezh2 functions, and to determine categorically if Ezh2 regulates the terminal differentiation process in vivo. By coupling these observations with genome-wide mapping of H3K27me3 and H3K4me3 in isolated primary satellite cells, we identify a fundamental role for H3K27me3 in suppressing factors associated with lineage specification within a proliferative progenitor population. This approach has identified that Ezh2 does not function during the terminal differentiation process or regulate terminal differentiation genes.
Results
Deletion of Ezh2 within satellite cells impairs muscle growth and survival
We characterised the functional consequence of deleting Ezh2 within the satellite cell pool using a Pax7-Cre (Ezh2P strain). Skeletal muscle derived from the Ezh2P mouse line exhibited an absence of Ezh2 protein and a global reduction in H3K27me3 at 7 weeks of age (Fig. 1A). In agreement with Juan et al. (Juan et al., 2011), the Ezh2P mouse line exhibited significant growth retardation (Fig. 1B), which was in large part due to a decrease in muscle fibre area (Fig. 1D,E) and subsequent reduction in muscle mass (supplementary material Fig. S1A). However, in contrast to previous work, we observed a high degree of postnatal lethality – only 54% of Ezh2P mice survived until P50 (Fig. 1C). The Ezh2P line is on a C57/B6 background, which may account for the lethality differences observed between this study and that of Juan et al. (Juan et al., 2011). Loss of Ezh2 within the satellite cell lineage has been suggested to lead to a decline in satellite cell number (Juan et al., 2011). We also observed these defects, with a 65% reduction in the number of Pax7+ or M-Cadherin+ (two well-characterised satellite cell markers) satellite cells (Fig. 1E,F). This was confirmed by FACS analysis, with a 63% reduction in Vcam1+ CD34+ CD31− CD45− satellite cells in limb skeletal muscle (Fig. 1G). During embryonic development and early postnatal growth Pax7+ progenitor cells expand rapidly and are required to establish the quiescent satellite cell pool that is maintained throughout adult life (White et al., 2010). We therefore reasoned that the reduction in satellite cell number observed could occur during the establishment of the quiescent satellite cell population. Indeed, at P0 we observed a 15% reduction in mouse weight and a 40% reduction in Pax7+ progenitors in skeletal muscle (Fig. 1H–J), demonstrating that some of the lesions in the Ezh2P line are established prenatally. However, Ezh2P mice are born at Mendelian ratios (25 expected, 22 obtained), which is in contrast to the lethality that occurs postnatally (Fig. 1C). Additionally, the growth defect is exacerbated after birth, with Ezh2P mice being 50% smaller than littermates by P10 (Fig. 1B).
Conditional knockout of Ezh2 impairs skeletal muscle growth. (A) Western blot for Ezh2 (β-actin control) and H3K27me3 (total H3 control). Ezh2fl-P, Ezh2 fl/fl without Pax7-Cre; Ezh2P, Ezh2 fl/fl with Pax7-Cre. (B) Whole-body mouse weights (male); means±s.e.m.; n = 11 per time point. ***P<0.001, significantly different from Ezh2fl-P, assessed by ANOVA. (C) Kaplan–Meier survival curves for the knockout mouse strains; Ezh2P, n = 35; Ezh2fl-P, n = 56. (D) Proportion of myofibres per cross-sectional area (CSA) for the Ezh2P strain (700 myofibres assessed per n; n = 3); mean CSA is shown. (E) Pax7 and laminin staining of cross-sections of TA muscle; nuclei are counterstained with DAPI (blue); arrowheads indicate Pax7+ cells. Scale bars: 100 µm. (F) Quantification of the average number of Pax7- or M-cadherin-positive cells Scale bars: 100 µm. per field of view (40×); means+s.e.m.; n = 3. (G) Flow cytometry analysis of the proportion of Vcam1+ CD34+ CD31− CD45− satellite cells within a satellite cell preparation derived from limb skeletal muscle; percentage calculated from a total of 300,000 events per sample; independent observations are shown. (H) Pax7 staining (red) in cross-sections of newborn (P0) hind leg. DAPI (blue) is used to mark nuclei. Scale bars: 100 µm. (I) The average number of Pax7+ cells per field of view (40×) is significantly reduced in the newborn (P0) Ezh2P mice. Means+s.e.m.; n = 3. **P<0.01 assessed by Student's t-test. (J) Ezh2P mice are significantly smaller than littermate controls at birth (P0); means+s.e.m.; Ezh2P n = 8, Ezh2fl-P n = 17; ***P<0.001 assessed by Student's t-test.
Ezh2 is required for muscle regeneration but not satellite cell self-renewal
The Ezh2P defects (Fig. 1H–J) suggest that the phenotype observed in the adult could be due to a failure to expand the progenitor population during development. To directly test whether or not Ezh2 has a role in satellite cell self-renewal, we employed a muscle regeneration protocol, which measures the ability of stem cells to repair skeletal muscle and maintain the stem cell population. Muscle regeneration was assessed 5 (early phase of repair), 10 (late phase of repair) and 25 (fully regenerated) days following cardiotoxin (CTX)-induced muscle damage. Muscle regeneration was nearly completely absent from the Ezh2P mice 25 days post CTX injection – only a small number of centrally nucleated fibres (a hallmark of newly formed myofibres) were observed (Fig. 2A,B). This failure to regenerate was accompanied by an accumulation of calcium deposits (Alizarin Red), fatty deposits (Oil Red O) and fibrous tissue (Van Gieson) (Fig. 2A). By analysing regeneration across a time course we were able to show that Ezh2 is required for effective muscle regeneration and does not simply delay muscle repair, as had been suggested by analysis of the tissue at 7 days after CTX injection (Juan et al., 2011).
Severely impaired muscle regeneration in Ezh2 null muscle. (A) Cross-sections of TA muscle taken from the Ezh2P strain 25 d after CTX injection and stained for haematoxylin and eosin (H and E), Alizarin Red (red/brown marks calcium deposits), Oil Red O (fat accumulation is marked in red) and Van Gieson (highlights fibrous tissue in red/pink), as indicated. Scale bars: 100 µm. (B) Proportion of muscle fibres within the damaged area with centralised nuclei. Means+s.e.m.; n = 3. ***P<0.001, assessed by Student's t-test. (C) Cross-section of a damaged area 5 days after cardiotoxin injection stained for Pax7 (progenitor cells) or haematoxylin and eosin (H and E), as indicated; nuclei are counterstained with DAPI (blue). Scale bars: 100 µm. (D) The average number of Pax7+ cells was determined from muscle sections; undamaged, 5, 10 and 25 days after CTX injection; means±s.e.m.; n = 3. **P<0.01, ***P<0.001, significantly different from Ezh2fl-P, assessed by ANOVA. (E) Ki67 was used to quantify proliferating cells in the Ezh2P strain; the average number of Ki67+, Pax7+ and double-positive cells was quantified within a field of view (40×) 5 days after CTX injection; means+s.e.m.; n = 3. **P<0.01, ***P<0.001, assessed by Student's t-test. (F) FACS isolated satellite cells cultured for 120 h, stained for Ezh2, Pax7 and MHC; nuclei were counterstained with DAPI (blue). Scale bars: 100 µm.
During regeneration satellite cells expand in response to injury, repair the damaged area and return to a quiescent state to ensure the maintenance of the stem cell pool (Shea et al., 2010). Using the CTX regeneration model it is therefore possible to examine stem cell self-renewal in vivo. The number of Pax7+ cells was quantified during the regeneration time course (Fig. 2C,D). As expected, the number of Ezh2fl-P satellite cells increased ninefold in response to the regenerative stimulus, returning to baseline by 25 days after injury. However, Ezh2 null Pax7+ cells were unable to expand during regeneration and their number remained at pre-regeneration levels (Fig. 2D). This failure to expand in response to regenerative signals is due to significantly reduced proliferation, as indicated by Ki67+ nuclei (Fig. 2E). No increase in apoptosis was observed, as determined by cleaved caspase 3 staining (supplementary material Fig. S1B). Crucially, the numbers of Ezh2 null satellite cells had not declined following the regeneration period (25 days).
These conclusions were also confirmed ex vivo. Ezh2 null satellite cells were isolated by FACS and cultured to induce activation, proliferation and differentiation. Following culture for 120 h after isolation, Ezh2P satellite cells showed a dramatic reduction in myoblast (Pax7+) number (Fig. 2F). However, Ezh2 null satellite cells were still capable of undergoing terminal differentiation (MHC+).
Ezh2 is not required for terminal differentiation or myofibre homeostasis
Without the actions of Ezh2 within Pax7+ progenitors there is a severe defect in muscle growth and regeneration. This appears to be due to the failure of the progenitor cells to expand and thus supply the progeny for terminal differentiation, although it could be the result of a defect during the terminal differentiation process. Indeed, Ezh2 is expressed during muscle terminal differentiation, albeit at reduced levels (supplementary material Fig. S2A–C). Consequently, we set out to further characterise the developmental window in which Ezh2 and H3K27me3 exert their influence. To achieve this, a second conditional mouse line was generated by crossing Myogenin-Cre mice (Li et al., 2005) with the Ezh2 fl/fl strain, termed Ezh2M. This specifically deleted Ezh2 at the onset of terminal differentiation, leaving the quiescent satellite cell and proliferating MyoD+ populations intact. Skeletal muscle derived from the Ezh2M mouse line exhibited an absence of Ezh2 protein, and interestingly, a global reduction in H3K27me3 levels (Fig. 3A). The results observed in the Ezh2P mice were in marked contrast to the phenotype of the Ezh2M mice, which showed no defect in total body weight, fibre size or survival (Fig. 3B,C; data not shown). The number of satellite cells, the ability of Ezh2M skeletal muscle to regenerate (5–25 days) and the self-renewal of the satellite cell population was comparable to the Ezh2fl-M controls (Fig. 3D–F). Thus, importantly, it can be concluded that the phenotype observed in the Ezh2P line was due to intrinsic defects within myoblast populations prior to the onset of terminal differentiation. Consequently, the Ezh2P phenotype was not due to a deficiency in terminal differentiation, myoblast fusion or myofibre hypertrophy.
Conditional knockout of Ezh2 during terminal differentiation. (A) Western blot for Ezh2 (with a β-actin control) and H3K27me3 (with a total H3 control). Ezh2fl-M, Ezh2 fl/fl without Myogenin-Cre; Ezh2M, Ezh2 fl/fl with Myogenin-Cre. (B) Whole-body mouse weights (male); means±s.e.m.; n = 5 per time point. (C) Proportion of myofibres per cross-sectional area (CSA) for the Ezh2M strain (700 myofibres assessed per n, n = 3); mean CSA is shown. (D) Quantification of the average number of Pax7- or M-cadherin-positive cells per field of view (40×); means+s.e.m.; n = 3. (E) Haematoxylin and eosin staining of cross-sections of the TA muscle from the Ezh2M strain after injection of cardiotoxin (CTX). Muscle was harvested at 5, 10 and 25 days after injection. Scale bars: 100 µm. (F) The average number of Pax7+ cells was determined from muscle sections; undamaged, 5, 10 and 25 days after CTX injection; means±s.e.m.; n = 3.
Genome-wide mapping of H3K27me3 and H3K4me3 in muscle satellite cells
The known targets of Ezh2 in stem cells, is limited to hair follicle stem cells (Lien et al., 2011) and haematopoietic stem cells (HSCs) (Adli et al., 2010; Cui et al., 2009), with little information regarding the skeletal muscle lineage. We therefore set out to determine the genomic targets of Ezh2 and the interplay between repressive (H3K27me3) and active (H3K4me3) histone modifications in satellite cells. We optimised a flow-cytometry purification protocol (Fukada et al., 2007) to isolate a near pure population of satellite cells (Fig. 4A). The majority of isolated cells expressed the muscle stem cell markers Pax7 (92±3%; n = 3; Fig. 4B) and M-cadherin (93%+5%; n = 3). We developed a chromatin immunoprecipitation protocol for FACS isolated muscle satellite cells and performed next generation sequencing (ChIP-Seq), which was used to characterise the localisation of H3K4me3 and H3K27me3 within the satellite cell genome (supplementary material Tables S1 and S2, respectively). The majority (96%) of H3K4me3 regions were in close proximity (<2 kb) to annotated genes. H3K27me3 peaks were distributed broadly throughout the genome, with 57% of enriched regions in close proximity (<2 kb) to annotated genes. We validated the ChIP-Seq data sets by examining H3K4me3 and H3K27me3 regions with different levels of enrichment and performed independent ChIP analysis coupled to quantitative PCR (ChIP-qPCR; Fig. 4C,D).
ChIP-Seq on satellite cells. (A) Representative FACS plots for the isolation of satellite cells. The final cell population is DAPI− CD31− CD45− Vcam1+ CD34+. A total of 30,000 events are shown; gates are set according to single-stained controls. (B) FACS-isolated cells were cytospun, and immunofluorescence for Pax7 (red) was performed; nuclei are marked by DAPI staining (blue). The arrowhead indicates a Pax7− cell. (C,D) ChIP-qPCR for indicated gene promoters, for H3K4me3 (C) or H3K27me3 (D); chromatin was obtained from freshly isolated satellite cells; data are normalised to input and presented as fold enrichment relative to an intergenic region on chromosome 1 (intergenic); means±s.e.m; H3K4me3, n = 3; H3K27me3, n = 4. ***P<0.001, significantly different from Oct4 (H3K4me3) or Gapdh (H3K27me3), assessed by ANOVA. (E) Read-density plot for a7-integrin, Cxcr4, Ncam1, Notch1, Sprouty1 and Sox8. Peaks indicate read density. (F) Plot for the average read density±4 kb of the transcriptional start site (TSS) of all annotated genes; H3K4me3, blue; H3K27me3, red; input, green.
Known satellite cell factors Itga7, Ncam1, Cxcr4, Pax7 and Myf5 (Fig. 4E; supplementary material Table S1) were enriched for the active histone modification H3K4me3. Furthermore, factors that have been proposed to regulate satellite cell function, Sprouty1 (Shea et al., 2010), Notch1 (Brack et al., 2008) and Sox8 (Schmidt et al., 2003), were also marked by H3K4me3 (Fig. 4E), suggesting that our data set may provide a useful resource of novel regulators of satellite cell function that could be examined in future functional studies (supplementary material Table S1).
Within the satellite cell genome, H3K4me3 was found within 1 kb of the transcriptional start site (TSS) (Fig. 4F). The H3K27me3 repressive mark was distributed broadly around promoters and deep into gene bodies (Fig. 4F). While this is in contrast to the distribution seen in ESC cells, where there is a peak of H3K27me3 at TSS of genes (Bernstein et al., 2006; Boyer et al., 2006; Guttman et al., 2009; Mikkelsen et al., 2007; Pan et al., 2007; Wang et al., 2008), it mirrors the localisation observed in more lineage restricted cell types, including fibroblasts and CD4+ T cells (Barski et al., 2007; Hawkins et al., 2010). The reduction in sequencing reads around the TSS observed for both H3K27me3 and input samples is a typical feature of ChIP-Seq studies (Asp et al., 2011; Barski et al., 2007; Cui et al., 2009; Pan et al., 2007; Young et al., 2011; Zhao et al., 2007), and most likely occurs due to the difficulties of sequencing GC rich TSSs.
In pluripotent ESCs, many genes encoding key developmental regulators are bivalent (Bernstein et al., 2006). However, we identified only 83 putative ‘bivalent’ domains within the satellite cell genome (supplementary material Table S3), supporting the notion that although bivalent domains may be a key feature of the ESC cell genome (containing >1600 bivalent domains), they may not be prevalent in the genome of the lineage-restricted satellite cell. GO ontology analysis on the bivalent domains (supplementary material Fig. S3B) indicated only a small enrichment for GO terms, suggesting that the bivalent domains that we identified do not share a common function.
H3K27me3 marks regulators of cell fate
The H3K27me3 mark was classically described as regulating the Hox clusters (Duncan, 1982; Soshnikova and Duboule, 2008) and this is also the case for muscle satellite cells, with broad H3K27me3 domains spread across Hox clusters (Fig. 5A). In the satellite cell population, H3K27me3 was also found at known non-myogenic cell fate regulators and genes associated with early development (examples are shown in Fig. 5B; GO analysis, Fig. 5C). However, factors involved in muscle development were not enriched in the H3K27me3 Gene Ontology (GO) analysis. By contrast, GO analysis indicated that the H3K4me3 was enriched at genes involved in muscle organ development and cellular homeostasis (Fig. 5D). Network analysis on the factors involved in muscle organ development suggests that there is crosstalk between these genes (supplementary material Fig. S3B). H3K27me3 was highly enriched at key transcription factors that regulate lineage decisions, such as members of the Pax, Sox, Gata and Foxo families (Fig. 5B; supplementary material Table S2); indeed all Pax transcription factors (Pax1–9), apart from Pax7, were marked by H3K27me3. Thus, in muscle stem cells, H3K27me3 is enriched at key developmental regulators associated with non-muscle lineages.
H3K27me3 marks regulators of cell fate. (A) Read-density plot for Hoxb and Hoxd clusters. Peaks indicate read density. (B) Schematic indicating a representative selection of lineage commitment genes marked by H3K27me3; genes are grouped according to the germ layers in which they function. (C,D) Gene ontology (GO) analysis for the 500 most enriched H3K27me3 (C) H3K4me3 (D) genes, performed using DAVID annotation clustering; representative GO categories are shown and are plotted against −log (p-value).
H3K27me3 and satellite cell identity
Our ChIP-Seq analysis identified that H3K27me3 marks a considerable number of key regulators of non-muscle development (Fig. 5B,C). Upon loss of Ezh2 there is a reduction of H3K27me3 at the promoters of Neurog1 (neuronal differentiation and cell type specification), Hoxb13 (foetal skin development), Gata4 (myocardial differentiation) and p16Ink4a (Fig. 6A). No statistically significant enrichment of H3K4me3 was observed for these loci (supplementary material Fig. S5).
Transcriptome changes in Ezh2-null satellite cells. (A) H3K27me3 ChIP-qPCR for indicated gene promoters; chromatin was obtained from freshly isolated satellite cells from Ezh2fl-P and Ezh2P mice; Ezh2P data are represented as a percentage of the Ezh2fl-P control; means+s.e.m.; n = 4. **P<0.01, *P<0.05 using one-sample t-test. (B) Venn diagram indicating overlap of transcripts, identified by RNA-Seq, upregulated greater than twofold in Ezh2P satellite cells when compared with Ezh2fl-P satellite cells (purple circle); genes (identified by Chip-Seq on Ezh2fl-P satellite cells) containing greater than 4-fold enrichment of H3K27me3 when compared with input (light blue circle). Overlap between two data sets is highlighted in darker blue. Bar charts indicate representative GO terms on analysis performed for each group. (C) Network analysis on upregulated (P<0.05) transcripts and their interactors (defined by STRING analysis). Factors involved in cell cycle progression are marked in blue.
To determine transcriptome changes that occur upon the loss of Ezh2 and subsequent reduction of H3K27me3, we performed RNA-Seq analysis on Ezh2 null FACS isolated satellite cells. We observed significant upregulation of a number of transcripts within the Ezh2P strain compared to the Ezh2Fl-P strain (supplementary material Table S4). We also noted a considerable number of transcripts that had low sequence reads (corresponding to low/no expression) in Ezh2fl-P satellite cells, but which showed greater than twofold increase in the Ezh2P satellite cells (Fig. 6B). However, only 11.4% (71/624) of the genes marked by H3K27me3 in satellite cells were upregulated greater than twofold in Ezh2 null satellite cells (Fig. 6B). Interestingly, these genes corresponded to factors involved in the mesenchymal lineage, cell fate commitment and transcription factor activity. Genes marked by H3K27me3 and not upregulated upon Ezh2 loss, were also highly enriched for factors involved in cell fate commitment. However, a significant proportion was associated with factors involved in embryonic morphogenesis, epithelium development and kidney development – i.e. lineages that are very distinct from the skeletal muscle lineage (Fig. 6B). Thus, the RNA-Seq analysis suggests that only a proportion of the genes marked by H3K27me3 in satellite cells are suppressed by this modification. Instead, it implies that a large proportion of H3K27me3 marked genes may be suppressed by further layers of epigenetic repression or may require transcriptional machinery that is not present in the satellite cell.
To further examine the global consequences of Ezh2 loss we performed network and GO analysis on the transcripts upregulated in Ezh2P satellite cells. The upregulated transcripts have diverse functions (Fig. 6B), although a significant proportion (P = 8.3E-7) were involved in cell cycle regulation (Fig. 6C). The phenotype of the Ezh2P strain suggested that Ezh2 may suppress factors involved in cell cycle progression. One such locus, the Cdkn2a locus, has been identified as a key target of Ezh2 within pancreatic islet β cells (Chen et al., 2009), epidermal stem cells (Ezhkova et al., 2011; Ezhkova et al., 2009) and numerous cancers (Agger et al., 2009; Barradas et al., 2009). Indeed, the only cyclin dependent kinase (CDK) inhibitor marked by H3K27me3 in wild-type satellite cells is p16Ink4a (Fig. 4D; Fig. 7A), which is encoded by the Cdkn2a locus. Interestingly, in Ezh2 null satellite cells (Ezh2P) the Cdkn2a locus exhibited a 70% reduction in H3K27me3 enrichment when compared with wild-type (Ezh2fl-P) (Fig. 6A). This loss of H3K27me3 was accompanied by a considerable increase in p16Ink4a and p19Arf mRNA levels (Fig. 7B).
Ezh2 does not regulate muscle terminal differentiation genes. (A) Read-density plots for Cdkn2a; peaks indicate read density. (B) RT-qPCR for indicated transcripts. mRNA obtained from freshly isolated satellite cells from Ezh2fl-P and Ezh2P mice, normalised to Rpl32 and 18S; Ezh2P data are represented as a fold change relative to the Ezh2fl-P control; means+s.e.m.; n = 4. ***P<0.001, **P<0.01, *P<0.05 assessed using a one-sample t-test. (C) Read-density plots for Ckm, Myogenin and Myh4; peaks indicate read density. (D) FACS-isolated satellite cells cultured for 72 h, stained for Myogenin and Ki67; nuclei were counterstained with DAPI (blue). Scale bars: 100 µm. (E) Cross-section of a resting (undamaged) and damaged (5 days after cardiotoxin injection) TA muscle stained for Pax7 (progenitor cells) and Myogenin (marks terminally differentiating cells), as indicated; nuclei are counterstained with DAPI (blue); higher-magnification images indicate Myogenin+/− nuclei. Scale bars: 100 µm. (F) Quantification of the number of terminally differentiating myoblasts (Myogenin+), average number of Myogenin+ nuclei per field of view (40×); means+s.e.m.; n = 3. ***P<0.001, assessed by Student's t-test.
Ezh2 and H3K27me3 do not suppress muscle terminal differentiation genes
By combining the ChIP-Seq data set with RNA-Seq and RT-qPCR (Fig. 6C) analysis of the Ezh2 null satellite cells, we were able to determine whether Ezh2 suppresses myoblast terminal differentiation. Surprisingly, our ChIP-Seq study shows that H3K27me3 does not in fact mark key muscle differentiation genes in satellite cells. We determined that MhcIIb, Myogenin and Ckm were marked in C2C12 cells (supplementary material Fig. S2B), but no H3K27me3 enrichment was detected by ChIP-Seq (Fig. 7C) or ChIP-qPCR (Fig. 4D) on freshly isolated satellite cells. Importantly, loss of Ezh2 in satellite cells does not lead to an upregulation of Myogenin, Ckm or MhcIIb transcripts (Fig. 7B). GO analysis on the RNA-seq data set indicates that genes associated with muscle development/differentiation are not significantly (P = 0.5) upregulated in Ezh2 null satellite cells (Fig. 6B).
To confirm these results we examined the capability of Ezh2 null satellite cells to commit to terminal differentiation both ex vivo and in vivo. FACS purified Ezh2P satellite cells committed to the terminal differentiation programme (as marked by Myogenin expression) at an equivalent rate to the Ezh2fl-P controls (Fig. 7D). These results were mirrored in vivo – resting muscle of either Ezh2fl-P or Ezh2P mice showed no expression of Myogenin (Fig. 7E). Additionally, the defects within the progenitor population led to a 70% reduction in Myogenin+ cells in the Ezh2P line during muscle regeneration (Fig. 7F). Thus, Ezh2 does not suppress the muscle differentiation programme in satellite cells.
Discussion
During muscle development, progenitor cells undergo multiple rounds of cell fate restriction, ultimately leading to the creation of lineage-committed satellite cells. By mapping the genome-wide pattern of H3K27me3 in freshly isolated satellite cells, we have identified that a considerable number of key developmental factors involved in fate choice and lineage commitment are marked by this repressive histone modification. These data, combined with the characterisation of two conditional knockout mouse models, outline the critical roles for Ezh2 within satellite cells: to suppress a subset of regulators of non-muscle cell fate, and to ensure appropriate progenitor cell expansion prior to the onset of terminal differentiation.
Utilising the above strategy we have conclusively shown that H3K27me3 marks genes associated with alternative lineage choice, but that it does not suppress skeletal muscle terminal differentiation. Myogenin, Ckm and MhcIIb have been identified as key effectors of Ezh2 function in the C2C12 cell line (Asp et al., 2011; Caretti et al., 2004; Juan et al., 2009; Seenundun et al., 2010). However, we have determined that muscle terminal differentiation genes are not marked by H3K27me3 in satellite cells. Furthermore, loss of Ezh2 did not lead to an upregulation of the expression of these genes, and Ezh2 null satellite cells did not spontaneously differentiate [as would have been predicted if master regulators of terminal differentiation (e.g. Myogenin) were expressed]. This disparity between C2C12 cells and satellite cells could be due to several factors. C2C12 cells are immortalised myoblasts that are highly proliferative [they do not contain an intact Cdkn2a locus (Pajcini et al., 2010)]. They are also a heterogeneous population (Olguin and Olwin, 2004). For example, the master regulator of satellite cells, Pax7 (Lepper et al., 2009; Seale et al., 2000) is only expressed in 50% of the cell population (supplementary material Fig. S4D) and (Olguin and Olwin, 2004). In addition, regulators of adipogenesis and osteogenesis, C/EBP, PPAR-y and Runx2, are all marked by the active H3K4me3 mark in C2C12 cells (Asp et al., 2011) but were not marked in satellite cells (supplementary material Table S1). Thus, the H3K27me3 enriched promoters that have been identified within the C2C12 cell line may only exist in a subpopulation of C2C12 cells, and these may not be of true myogenic lineage.
When Ezh2 was ablated via the Pax7 promoter, we observed a significant failure in progenitor myoblast expansion, which resulted in regeneration defects and accumulation of non-muscle tissue. This was not observed when Ezh2 was deleted during terminal differentiation (Ezh2M). Utilisation of the Pax7-Cre can lead to Cre-mediated recombination within the central nervous system (Jostes et al., 1990). However, the satellite cell expansion phenotype was also observed ex vivo, and deletion of Ezh2 with a MyoD-Cre, which is specific to activated satellite cells, has shown a similar phenotype to that of the Ezh2P strain. Therefore, the muscle defects observed in the Ezh2P strain are likely to result from the loss of Ezh2 within the muscle lineage. However, we cannot exclude the possibility that the growth defects and postnatal lethality observed in the Ezh2P line may occur, in part, due to central nervous system defects.
Even after severe regenerative stimulus, no decline in the number of quiescent satellite cells could be detected in the Ezh2P muscle. Thus, the magnitude of proliferation and differentiation of progenitor populations during muscle regeneration may not influence the replenishment of the satellite cell pool. This notion agrees with the identification of a population of asymmetrically dividing satellite cells that maintains the stem cell pool but does not undergo expansion during regeneration (Kuang et al., 2007; Le Grand et al., 2009). A limited regenerative response was observed in the Ezh2P line and a small number of differentiating myoblasts were present during regeneration. Ezh2 null satellite cells are therefore still capable of differentiation, but it appears that the expansion of progenitor cells is impeded, and the supply of differentiating myoblasts is consequently reduced, leading to a very limited regenerative response.
Within pancreatic islet β cells (Chen et al., 2009) and epidermal basal cells (Ezhkova et al., 2009), the regulation of p16Ink4a by Ezh2 had a powerful influence on proliferation and was shown to directly regulate the expansion of progenitor populations. This common function of Ezh2 also occurs within satellite cells, with a 60-fold upregulation of p16Ink4a in Ezh2 null satellite cells. Indeed, knockdown of p16Ink4a can rescue some of the proliferation defects observed in Ezh2 null satellite cells (Juan et al., 2011). Interestingly, p16Ink4a was the only cyclin dependent kinase inhibitor marked by H3K27me3 in satellite cells. However, RNA-Seq analysis indicates that p16Ink4a may be part of a network of factors that regulate cell cycle progression in Ezh2 null satellite cells.
RNA-Seq analysis of the Ezh2 null satellite cells indicated that even though a large number of developmental regulators are marked by H3K27me3, only a subset of these factors become expressed upon the loss of Ezh2. H3K27me3 is therefore fundamental for regulating a subset of target loci, but at a number of loci, Ezh2 is not required for gene repression. It is tantalising that the factors not upregulated upon loss of Ezh2, are associated with very distant lineages (i.e. epithelium) from that of skeletal muscle. It is possible that these factors are ‘permanently’ suppressed by layers of epigenetic modifications, which may include DNA methylation and H3K9me3 (Henikoff and Shilatifard, 2011). However, some factors such as Pax3 and Cdkn2a, are upregulated upon the loss of Ezh2, and thus, H3K27me3 may be a dominant mark regulating their expression. Importantly, we have demonstrated that specific developmental factors targeted and suppressed by PRC2 appear to be highly dependent on stem cell lineage.
The maintenance of transcriptional repression is not only important for stem cell function but for tissue integrity in general. For example, mis-expression of the transcription factors Msx1 and Twist1 (marked by H3K27me3 in satellite cells) in adult myofibres, results in their de-differentiation and loss of tissue integrity (Hjiantoniou et al., 2008; Kumar et al., 2004; Odelberg et al., 2000). Evidence from the C2C12 cell line and from cultured satellite cells has suggested that Ezh2 may regulate gene repression at the onset of, and during, terminal differentiation (Asp et al., 2011; Palacios et al., 2010). However, we did not observe a skeletal muscle phenotype upon deletion of Ezh2 at the onset of myoblast terminal differentiation in vivo, despite a global reduction in H3K27me3 levels. Our work supports the growing body of evidence suggesting that Ezh2 functions within highly proliferative progenitor cell populations (Chen et al., 2009; Ezhkova et al., 2011; Ezhkova et al., 2009; Margueron et al., 2009; Su et al., 2003), as opposed to terminally differentiated tissues or quiescent stem cells. The dependency on Ezh2 during the proliferative phase may reflect the necessity to pass on the information contained within the histone code through multiple cell divisions, which requires the faithful copying of the histone code to progeny during mitosis (Hansen et al., 2008; Petruk et al., 2012). In post-mitotic tissues and quiescent cells it is not necessary for the histone code to be copied and therefore these cell types may not require the actions of Ezh2.
We have identified the developmental window in which Ezh2 acts in the muscle lineage and established that it does not regulate the terminal differentiation process. Instead, Ezh2 promotes satellite cell expansion, ensuring that a critical mass of progenitor cells is achieved prior to the commitment to terminal differentiation. By characterising target loci marked by H3K27me3 and H3K4me3, the study presented here enhances our understanding of the mechanisms that define adult stem cell identity. Determining how tissue specific stem cells are restricted to a particular fate choice may enable greater understanding of the mechanisms underlying stem cell self-renewal, and would aid the use of satellite cells for tissue regeneration, such as in muscular dystrophies.
Materials and Methods
Mouse strains
To conditionally delete Ezh2, mice containing floxed alleles of the SET domain of Ezh2 (Su et al., 2003), provided by Alexander Tarakhovsky (Rockefeller University), were crossed with mice expressing either Pax7-Cre (Keller et al., 2004) [provided by Mario Capecchi (University of Utah)] or Myogenin-Cre [provided by Eric Olson (University of Texas)]. The Pax7-Cre transgene was used to delete Ezh2 within the satellite cell lineage, while the Myogenin-Cre ensured deletion of Ezh2 upon the initiation of differentiation. All mice were homozygous for the Ezh2 flox alleles with (Ezh2P or Ezh2M) or without (Ezh2fl-P or Ezh2fl-M) one copy of the Pax7-Cre or Myogenin-Cre transgene. All experiments were performed on mice 7–9 weeks of age [satellite cells are quiescent at this time point (Fukada et al., 2007)] and compared with littermate controls. For CTX-induced muscle regeneration, mice were killed at 5, 10 or 25 days after injection of 100 µl CTX (10 µM; cardiotoxin) or 100 µl PBS (vehicle control) into the tibialis anterior (TA) muscle. All mouse husbandry and experiments were carried out according to the Babraham Institute animal ethics committee under the terms of the UK Home Office.
FACS isolation of satellite cells
To isolate satellite cells, dissected hind limb muscles were digested in 0.25% collagenase type II (Lorn Labs) and 0.2% Dispase II (Roche) in DMEM for 90 min at 37°C with agitation. Satellite cells were mechanically dissociated from myofibres, by passing the tissue suspension through a 19-gauge needle. These cells were purified by fluorescent activated cell sorting on a FACSAria II (BD Biosciences). Antibodies used were CD31-PE (1∶100, Abcam), CD45-PE (1∶100, B.D. Biosciences), CD34-AF647 (1∶33, eBiosciences), Vcam1-biotin (1∶15, B.D. Biosciences) and Streptavidin-AF488 (1∶15, Invitrogen). Dead cells and debris were excluded by DAPI staining and by gating on forward and side scatter profiles. Cells were gated for the Vcam1+ CD31− and CD45− population; the satellite cell population was further enriched by gating on CD34+ cells. The final population was Vcam1+, CD34+, CD31− and CD45−. Satellite cells were cultured in high serum medium (20% FBS, 1% l-glutamine and 1× antibiotic antimycotic solution in DMEM) supplemented with basic FGF (2.5 ng/ml, Promega).
RT-qPCR
Total RNA was isolated from 3×105 FACS purified satellite cells with TRIzol reagent (Invitrogen) and cDNA synthesis was performed using the Quantitect Reverse Transcriptase kit (Qiagen). Quantitative real-time PCR was performed using SYBR Green Jumpstart qPCR Mix (Sigma). Expression levels were normalised to Rpl32 (60 s ribosomal protein L32) and 18 s (18 rRNA). RT-qPCR primers are detailed in supplementary material Table S5.
Chromatin immunoprecipitation
ChIP was carried out using ChIP grade antibodies, H3K27me3 (07-499, Millipore) and H3K4me3 (pAB-003-050, Diagenode). ChIP was performed essentially as previously described (Adli et al., 2010). Briefly, freshly isolated satellite cells (5×105) or C2C12 cells (8×105) were fixed in 4% formaldehyde for 10 min at room temperature and fixation quenched with 150 mM glycine. Chromatin was sheared by sonication to an average fragment size of between 200 and 300 bp using a Bioruptor sonication instrument (Diagenode). Immunoprecipitation was performed using Protein A DynaBeads (Invitrogen). The beads were washed three times with 900 µl of each of the following buffers at 4°C: RIPA buffer, Radioimmunoprecipitation assay buffer, (10 mM Tris-HCl pH 7.5; 140 mM NaCl; 1 mM EDTA; 0.5 mM EGTA; 1% [v/v] Triton X-100; 0.1% [w/v] SDS; 0.1% [w/v] Na-deoxycholate; 10 mM NaF; 1 mM NaVO3; 5 µg/ml aprotinin; 10 µg/ml leupeptin; 1 mM PMSF), high-salt RIPA buffer (RIPA buffer with 500 mM NaCl) and TE (10 mM Tris-HCl pH 8.0; 1 mM EDTA). After reverse crosslinking, enriched DNA fragments were extracted with phenol/chloroform followed by isopropanol precipitation. Immunoprecipitated DNA was either analysed by real-time qPCR or prepared for sequencing. For ChIP-qPCR, enrichment was quantified by real-time qPCR as a percentage of input DNA and represented as a fold change relative to a control region. Primers employed in this study are detailed in supplementary material Table S5.
ChIP-Seq
For ChIP-seq analysis, 4×106 freshly sorted satellite cells were collected and pooled from ten independent FACS experiments (160 male mice in total); 2×106 satellite cells were used per ChIP. Subsequently, the DNA was adaptor-ligated, amplified using standard protocols (Quail et al., 2008) and sequenced using an Illumina Genome Analyzer GAIIx. H3K4me3-ChIP, H3K27me3-ChIP and input DNA sequence reads were aligned to the mouse genome (assembly NCBIM37/mm9), excluding non-uniquely mapping, and duplicate reads. All analyses were performed using SeqMonk (http://www.bioinformatics.bbsrc.ac.uk/projects/seqmonk/). For identification of peaks, probes (360–2000 bp in length) were designed with a 5-fold enrichment over a random data set and only peaks with a read density that was 4-fold higher than the input were counted. Bivalent domains were defined as containing a H3K27me3 peak (as described above) within 2 kb of a H3K4me3 enriched (4-fold higher than input) promoter. GO analysis was performed using DAVID (Huang et al., 2008).
RNA-Seq
Total RNA was extracted from FACS sorted satellite cells using TRIzol reagent, the quality and the quantity of RNA was checked using Agilent Bioanalyser. Libraries for RNA-Seq were made using Illumina Truseq kit v2 (as per the manufacturer's instructions) and these were then sequenced using Illumina HiSeq. The sequences were then aligned to the mouse genome (NCBIM37/mm9). Three independent replicates, using eight mice per replicate, were used for both Ezh2fl-P and Ezh2P lines. Further analysis and normalization was performed using Seqmonk. Lists of differentially regulated genes were generated using either a statistical significance test (P<0.05) or an intensity difference test (log2 fold change) for the average of three replicates. Network analysis of differentially regulated genes was performed using STRING (Jensen et al., 2009) and Cytoscape (Cline et al., 2007).
Immunofluorescence and histology
Immunofluorescence on cultured C2C12 cells and satellite cells was performed as previously described (Paxton et al., 2011). TA muscles were harvested and snap frozen in liquid-nitrogen-cooled isopentane, mounted in OCT and 10 µm cryosections collected. Histological analysis was performed using haematoxylin and eosin (morphological analysis), Van Gieson (fibrous tissue), Oil Red O (fatty deposits) or Alizarin Red (calcium deposits) staining using standard histological protocols. For immunofluorescence, cryosections were fixed using 4% paraformaldehyde (PFA) for 10 min, and permeabilised using 0.2% Triton X/PBS. Endogenous IgG was blocked using the MOM kit (Vector Labs). The cryosections were further blocked with 5% goat serum and 3% BSA in 0.2% Triton X/PBS. Both tissue sections and cultured myoblasts were incubated with primary antibodies overnight at 4°C.
Brightfield images were acquired using a Qimaging micro publisher 3.3RTV camera mounted on an Olympus BX41 microscope. For postnatal fibre size analysis, immunofluorescence for laminin (α-2 chain) was used to mark myofibre boundary, and cross sectional area was calculated using ImageJ software (NIH). Fibre types were identified by immunofluorescence with antibodies specific for MyHC subtypes. Total fibre number was quantified, with the aid of Image J (NIH), from tiled images of H and E stained cross-sections of the TA muscle.
The primary antibodies were incubated overnight at 4°C. The antibodies used were: Pax7 (DSHB), MyoD (Santa Cruz Biotechnology), Laminin α-2 chain (Enzo Life Sciences), Myogenin (Santa Cruz Biotechnology), Skeletal Fast Myosin (Sigma), Skeletal Slow Myosin (Sigma), pH 3 Ser10 (Cell Signaling Technology), Ki67 (Leica), Ezh2 (Cell Signaling Technology), M-cadherin (Santa Cruz Biotechnology) and Cleaved Caspase 3 (Cell Signaling Technology).
Secondary antibodies used were donkey anti-rabbit AF488, goat anti-mouse (IgG1 specific) AF568 and goat anti-rat AF633 (Invitrogen). These were added at a dilution of 1∶300 for 1 h at room temperature together with DAPI to label nuclei. Slides were mounted with Vectashield hardset (Vector Labs).
Statistical analysis
For all quantitative analyses presented, a minimum of three replicates were performed. Data are presented as means±s.e.m. ANOVA with post-hoc tests and Student's t-tests were performed as appropriate (indicated in the figure legends). Statistical analysis was carried out using GraphPad Prism.
Acknowledgments
We thank Alexander Tarakhovsky, Mario Capecchi and Eric Olson for the transgenic mice strains. We are grateful to Peter Rugg-Gunn for critical reading of the manuscript and helpful discussions. The GAIIX sequencing was performed by Kristina Tabbada. Sequencing analysis was performed with the help and advice of Simon Andrews and Felix Krueger. The authors declare that they have no conflict of interest.
Footnotes
↵*Present address: Centre for Stem Cells and Regenerative Medicine, Kings College London, London SE1 9RT, UK
↵‡Present address: Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
Funding
We thank the Biotechnology and Biological Sciences Research Council (BBSRC) for competitive strategic funding to J.M.P. via the Babraham Institute, studentship funding for S.W. and P.B., and Responsive Mode funding [grant number BB-H019243-1 to D.P. and J.M.P.].
Supplementary material available online at http://jcs.biologists.org/lookup/suppl/doi:10.1242/jcs.114843/-/DC1
- Accepted November 2, 2012.
- © 2013. Published by The Company of Biologists Ltd