The molecular mechanism responsible for a decline of stem cell functioning after replicative stress remains unknown. We used mouse embryonic fibroblasts (MEFs) and hematopoietic stem cells (HSCs) to identify genes involved in the process of cellular aging. In proliferating and senescent MEFs one of the most differentially expressed transcripts was Enhancer of zeste homolog 2 (Ezh2), a Polycomb group protein (PcG) involved in histone methylation and deacetylation. Retroviral overexpression of Ezh2 in MEFs resulted in bypassing of the senescence program. More importantly, whereas normal HSCs were rapidly exhausted after serial transplantations, overexpression of Ezh2 completely conserved long-term repopulating potential. Animals that were reconstituted with 3 times serially transplanted control bone marrow cells all died due to hematopoietic failure. In contrast, similarly transplanted Ezh2-overexpressing stem cells restored stem cell quality to normal levels. In a “genetic genomics” screen, we identified novel putative Ezh2 target or partner stem cell genes that are associated with chromatin modification. Our data suggest that stabilization of the chromatin structure preserves HSC potential after replicative stress.
Introduction
Hematopoietic stem cell (HSC) self-renewal is driven by both intrinsic and extrinsic factors, but the molecular mechanism specifying whether developmental potential is lost or retained during asymmetric cell divisions is unknown. Serial transplantation studies have clearly indicated that self-renewal potential of HSCs is impaired after replicative stress.1,2 HSC activity may be irreversibly lost in a single cell division,3,4 indicating that the epigenetic regulation of gene expression, largely dictated by the histone code, may play an important role. Recently, substantial attention has focused on the role of Polycomb group (PcG) proteins in stem cell self-renewal. To date, 2 distinct Polycomb repressive complexes, (PRCs) 1 and 2, have been identified. Mammalian PRC1 contains Cbx, Mph, Ring, Bmi1, and Mel185 and is thought to be important in the maintenance of gene repression. The second complex, PRC2, contains Ezh2, Eed, and Su(z)125 and is thought to be involved in initiation of gene repression. PcG complexes are targeted to cis-regulatory Polycomb response elements (PREs) by DNA-binding transcription factors.
The role of only few PcG proteins in murine hematopoiesis has been established. Melo18 (Pcqf2) negatively regulates self-renewal of HSCs because its loss leads to an increase of HSCs in G0 and to enhanced HSC self-renewal.6 Mph1/Rae28 (Phc1) mutant mice are embryonic lethal because HSC activity in these animals is not sufficient to maintain hematopoiesis during embryonic development.7 Bmi1 (Pcgf4)-/- HSCs derived from fetal liver were not able to contribute to long-term hematopoiesis in competitive repopulation experiments, demonstrating cell autonomous impairment of their self-renewal potential.8 In contrast, overexpression of Bmi-1 extends the replicative lifespan of mouse and human fibroblasts and causes lymphomas in transgenic mice.9 Heterozygosity for a null allele of Eed caused marked myeloproliferative and lymphoproliferative defects, indicating a negative regulation of cell cycle activity of both lymphoid and myeloid progenitor cells.10 Eed, which is contained in the PRC2 complex, seems to have an opposite effect on hematopoiesis than proteins that are present in PRC1.
Deficiency of Ezh2, also present in the PRC2 complex, resulted in lethality at early stages of mouse development. Mutant blastocysts either ceased developing after implantation or failed to complete gastrulation.11 Conditional deletion of Ezh2 in HSCs resulted in impaired T-cell generation after transplantation.12 Expression levels of EZH2 increased with differentiation of bone marrow cells.13 In contrast, EZH2 was down-regulated on differentiation of HL-60 cells to mature granulocytes.14 In humans EZH2 is commonly up-regulated in different cancer types.15-20 Inhibition of Ezh2 levels by antisense oligonucleotides or siRNA experiments led to cell cycle arrest.14,15,21 Recently, it has been shown that EZH2 is down-regulated in senescent human fibroblasts,22 implying a role for EZH2 during cellular senescence.
During development, PcG protein complexes are thought to maintain long-term and heritable gene silencing through local alterations of the chromatin structure.23 Methylation of the N-terminal tail of histone H3 causes stable conformational changes in chromatin that define the status of gene expression.24,25 The Drosophila Enhancer of zeste (E(z)) protein contains histone methyltransferase (HMT) activity in its SET domain, which is specific for lysines 9 and 27 of H3.26 It has been suggested that E(z) plays an important role not only in the establishment of PcG-specific methylation marks but also in their maintenance during replication, as well as after recruiting PRC1.27 E(z)-null Drosophila showed general decondensation of chromatin,28 reflected by increased chromosome breakage and a low mitotic index.29 In mammalian cells, histone deacetylases (HDACs) 1 and 2 are also present in the Eed-Ezh2 complex,30 suggesting that transcriptional repression of this complex is also mediated via histone deacetylation. In addition, PcG genes are known to play an important role in maintaining the silent state of Hox gene expression.31 Overexpression of specific Hox genes resulted in increased HSC pool size.32
Although an important function for several PcG genes in hematopoiesis has been established, the role of Ezh2 has not yet been elucidated. In the current study, Ezh2 was identified as a senescence-preventing gene in a genetic screen using young and senescent mouse embryonic fibroblasts. Moreover, Ezh2 was abundantly expressed in highly purified HSCs and was rapidly down-regulated on differentiation. We established that overexpression of Ezh2 in HSCs fully preserves stem cell potential and prevents HSC exhaustion even after multiple rounds of serial transplantation. Finally, we identified putative novel partner or target genes of Ezh2 in stem cells.
Materials and methods
Animals
Timed pregnant female C57BL/6 mice were used to obtain day 14 postcoitus embryos from which mouse embryonic fibroblasts (MEFs) were isolated. Eight- to 12-week-old female B6.SJL-PtprcaPep3b/BoyJ (CD45.1) mice were used as donors for transplantation and were bred at the Central Animal Facility of the University of Groningen. Eight- to 12-week-old female C57BL/6 (CD45.2) mice were purchased from Harlan (Horst, The Netherlands) and were used as recipients.
MEFs
Day 14 postcoitus embryos were mechanically dissociated into single cells and cultured in DMEM with 10% FCS and penicillin and streptomycin (Invitrogen, Breda, The Netherlands). MEFs were passaged every 3 to 4 days and viable cells were counted by trypan blue exclusion. Population doublings were calculated as per Todaro and Green.33
Gene array analysis
Total RNA was isolated from passage 1 (young p1) and passage 5 (aged p5) MEFs with a RNeasy mini kit (Qiagen, Venlo, The Netherlands). During cDNA synthesis, samples were labeled with [33P]-deoxycytidine triphosphate (dCTP; MP Biochemicals, Irvine, CA). Mouse filter gene arrays (GF400a; Research Genetics, Invitrogen) were hybridized and analysis was performed as described previously.34 Expression patterns were verified by reverse transcription-polymerase chain reaction (RT-PCR).
Retroviral vectors
The MIEV vector (a kind gift of Prof Dr C. Jordan, University of Rochester, Rochester, NY) contains an internal ribosomal entry site (IRES) sequence and the gene for enhanced green fluorescent protein (eGFP). This MIEV vector served as a control and backbone for cloning of an Ezh2 cDNA upstream of the IRES to create MIEV-Ezh2 vector (Figure 1D).
Ecotropic Phoenix packaging cells (Nolan Lab, Stanford, CA) were transfected with 2 μg plasmid DNA and Fugene 6 (Roche, Basel, Switzerland), 24 hours after plating 3 × 105 cells/well in a 6-well plate. Virus-containing supernatants from transfected ecotropic Phoenix packaging cells were used to infect MEFs and HSCs.
Retroviral overexpression of Ezh2 in MEFs
MEFs were plated at a density of 105 cells/well in a 6-well plate 24 hours prior to transduction. Viral supernatant was harvested 48 and 72 hours after transfection and added to the MEFs with addition of 8 μg Polybrene (Sigma, St Louis, MO). Three days later GFP+ cells were selected using a MoFlo flow cytometer (DakoCytomation, Fort Collins, CO). The experiment was performed 3 times in triplicate.
Retroviral overexpression of Ezh2 in HSCs
Primary bone marrow (BM) cells were transduced as previously described with some adjustments.35 Briefly, BM cells were extracted from mice (CD45.1) given intraperitoneal injections of 150 mg/kg 5-fluorouracil (5-FU; Pharmachemie Haarlem, Haarlem, The Netherlands) 4 days earlier. Cells were cultured for 48 hours in StemSpan (StemCell Technologies, Vancouver, BC, Canada) supplemented with 10% FCS, 300 ng/mL PEG-rrSCF (Amgen, Thousand Oaks, CA), 20 ng/mL rmIL-11 (R&D Systems, Minneapolis, MN), 10 ng/mL Flt3 ligand (Amgen), penicillin, and streptomycin. Viral supernatant was harvested 24 and 48 hours after transfection of ecotropic Phoenix cells and inoculated in 6-well plates that were coated with 50 μg retronectin (Takara, Kyoto, Japan). Plates containing the viral supernatant were spun for 1 hour at 800 g at room temperature. Four hours later viral supernatant was removed, and 7.5 × 105 cultured BM cells were inoculated per well together with 4 μg Polybrene (Sigma). At the second transduction 2 μg Polybrene was added. Four days later transduction efficiency was determined by flow cytometry (FACSCalibur; Becton Dickinson, Palo Alto, CA).
Transplantation
After transduction, 3 × 106 CD45.1 BM cells were transplanted into lethally irradiated (9.5Gy, IBL 637 137Cs γ source; CIS Biointernational, Gif-sur-Yvette, France) recipients (CD45.2). GFP+ cells were not selected after the transduction protocol, so transplants consisted of a mixture of transduced and nontransduced cells. At several points after transplantation, blood was drawn from the retro-orbital plexus and white blood cell (WBC) numbers were counted using a Coulter Counter (Beckman Coulter, Fullerton, CA). The percentage GFP+CD45.1+ leukocytes was assessed by flow cytometry (FACSCalibur; Becton Dickinson). Two independent experiments were carried out. The total number of primary recipients was 11 controls and 11 for Ezh2.
Serial transplantation
Four months after primary transplantation, recipients of the control and Ezh2 group were killed and BM cells were isolated by flushing 2 hind legs. Unfractionated BM cells were transplanted into lethally irradiated (9.5 Gy) secondary recipients (CD45.2) in limiting dilution with 106 freshly isolated C57BL/6 competitor BM cells (CD45.2). Reconstitution by the transduced HSCs was assessed for each recipient 3 months after serial transplantation by analysis of peripheral blood cell chimerism. Animals showing more than 1% GFP+ engraftment in myeloid and lymphoid lineages were considered to be reconstituted by the transduced HSCs. The competitive repopulating index (CRI) and the frequency of long-term repopulating cells present in primary recipients (LTRA) were calculated. The formula to calculate the CRI was: CRI = [(% CD45.1+GFP+ cells/% CD45.1+GFP- cells)/(ratio of CD45.1+GFP+/CD45.1+GFP- cells transplanted)]. The entire experiment was repeated twice. In total, 35 secondary recipients were given transplants with Ezh2-overexpressing cells and 35 recipients with control cells. In addition, a secondary transplantation without cotransplantation of competitors was performed. A total of 2 × 106 BM cells were transplanted into 10 lethally irradiated (9.5 Gy) recipients.
Four months after secondary transplantation, a tertiary transplantation was performed, identical to the secondary transplantation. For both control and Ezh2-overexpressing cells, a total of 28 recipients was used for competitive repopulation assays and 12 animals underwent transplantation without additional support of competitor BM cells.
Sorting hematopoietic cell populations
Different populations of hematopoietic cells were sorted using flow cytometry as described before.34 The 4 different populations (Lin-Sca-1- c-kit-, Lin-Sca-1+c-kit-, Lin-Sca-1-c-kit+, Lin-Sca-1+c-kit+) were directly sorted in RNA lysis buffer of the RNeasy mini kit (Qiagen), using a MoFlo flow cytometer (DakoCytomation) for gene expression analysis. In addition, Lin-Sca-1+c-kit+ cells were sorted for further culturing.
To assess clonogenic activity in purified stem cell fractions after primary transplantation, the cobblestone area-forming cell (CAFC) assay was performed as described previously.36 Early-appearing CAFCs day 7 correspond to relatively committed progenitor cells, whereas late-appearing CAFCs day 35 reflect more primitive cell subsets.
In vitro differentiation of purified HSCs
Purified Lin-Sca-1+c-kit+ (LSK) cells were cultured in IMDM (Invitrogen) containing 100 ng/mL PEG-rrSCF (Amgen) and 10 ng/mL GM-CSF (Behringwerke, Marburg, Germany) to stimulate rapid in vitro differentiation. Total cell number was assessed using trypan blue exclusion. At different time points during culturing RNA was isolated for gene expression analysis.
Western blotting
Cells were resuspended in PBS, lysed by addition of sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) sample buffer, and sonicated. Proteins were separated by 10% SDS-PAGE and transferred to nitrocellulose. After blocking in 3% skim milk, membranes were probed with a 1:100 dilution of mouse anti-Ezh2 mAb (M18EZH2, kindly provided by Prof Dr A. P. Otte, Swammerdam Institute for Life Science, University of Amsterdam, The Netherlands37,38 ), washed, and incubated with anti-mouse horseradish peroxidase-conjugated secondary antibody (Amersham Biosciences, Buckinghamshire, United Kingdom). Membranes were developed using enhanced chemiluminescence (ECL) reagents (Amersham Biosciences). Equal loading of membranes was verified with rabbit anti-γ-tubulin (Sigma).
Quantitative PCR
Total RNA was isolated using the RNeasy Mini kit (Qiagen) and standard cDNA synthesis was performed. Quantitative PCRs (qPCRs) were performed in triplicate. PCR amplification, using SYBR Green, was performed in 96-well microtiter plates in an iCycler thermal cycler (Bio-Rad, Hercules, CA). Sample cDNAs were compared with expression of house-keeping genes Gapdh and actin using the relative quantification ΔΔCT technique.39 Relative expression levels in the different LSK populations were estimated by first calculating the number of molecules formed at reaching the CT. By correcting this value for the initial number of cells used, relative expression was determined.
Search for putative targets or partners of Ezh2
To find putative targets or partners of Ezh2, we used The GeneNetwork, an online stem cell database.40 We have recently described the establishment of this database.41 The GNF Hematopoietic Cells U74Av2 (Mar04) RMA database was chosen and information on Ezh2 expression levels in different BXD strains was retrieved. Subsequently, linkage analysis as implemented in GeneNetwork40 was performed to determine by which locus differences in expression levels of Ezh2 were regulated. Furthermore, it is feasible to identify which stem cell transcripts show similar or opposite expression profiles as Ezh2 across the 30 BXD strains. These transcripts are potential partners or targets of Ezh2. The top 100 best correlating transcripts were selected and imported in WebGestalt (http://genereg.ornl.gov/webgestalt/, a web-based Gene Set Analysis Toolkit implemented in GeneNetwork40 ).
Statistical analysis
Statistical differences between means were assessed using the 2-tailed t test assuming unequal variances and ANOVA with 2 factors. Kaplan-Meier survival analysis was performed for mice that received a tertiary transplant of control or Ezh2-transduced HSCs. Differences in survival were tested for significance using a log-rank test. Statistical analysis was performed using the SPSS (Chicago, IL) statistical package.
Results
Ezh2 and cellular senescence in MEFs
An established in vitro model to study cellular senescence involves the serial passaging of MEFs.42 MEFs were cultured from 14-day-old C57BL/6 embryos. As expected, MEFs were able to undergo a limited number (8) of population doublings before they senesced (Figure 1A). To screen for genes that were differently expressed in young versus senescent MEFs, expression-profiling experiments were conducted. One of the most differentially expressed transcripts was Enhancer of zeste homolog 2 (Ezh2) (see Table S1 for all differentially expressed genes; see the Supplemental Table link at the top of the online article, at the Blood website). Differences in expression levels of Ezh2 in young and senescent MEFs were confirmed by RT-PCR (Figure 1B). Protein levels of Ezh2 in MEFs correlated with RNA levels; Ezh2 was more abundant at early passage (p1) compared to senescent MEFs (p9; Figure 1C). We next established the growth kinetics of MEFs in which Ezh2 was ectopically overexpressed using retroviral transductions. The retroviral vectors used to induce overexpression of Ezh2 are shown schematically in Figure 1D. As control we used a vector that only expresses eGFP (MIEV). The MIEV vector (Ezh2) contains Ezh2 and eGFP sequences, separated by an IRES. Western blotting failed to detect expression of Ezh2 in MEFs transduced with the control vector, whereas protein levels were clearly increased in cells transduced with the Ezh2 vector (Figure 1E). At different passages after retroviral transduction expression levels of Ezh2 were assessed by qPCR (Figure 1F). We confirmed that Ezh2 was readily overexpressed in MEFs transduced with the Ezh2 vector relative to cells transduced with the control vector. Levels of overexpression gradually increased from 40-fold to 70-fold at different time points after transduction. MEFs transduced with the control vector senesced rapidly after transduction (Figure 1G). In contrast, MEFs in which Ezh2 was overexpressed were able to escape senescence and showed continuous growth several passages after transduction (Figure 1G).
Ezh2 expression in HSCs
Different members of the PcG gene family have recently been implicated as having a role in HSC self-renewal.6-8 In contrast, the potential function of Ezh2 in long-term repopulating HSCs has not been investigated. To quantify Ezh2 expression in different hematopoietic cell populations, BM cells were isolated and stained with a cocktail of lineage (Lin)-specific antibodies and for the stem cell markers Sca-1 and c-kit. The 5% most lineage-negative fraction (Figure 2A) was further subfractionated in 4 different populations based on Sca-1 and c-kit expression (Figure 2B), and levels of Ezh2 were determined for all 4 fractions by qPCR (Figure 2C). In 2 separate experiments, highest expression levels were found in the Lin-Sca-1-c-kit+ population, which is known to be enriched for progenitors.43 The Lin-Sca-1+c-kit+ (LSK) population is highly enriched for the most primitive HSCs.44-47
Next, LSK cells were purified and differentiation was induced by culture in the presence of GM-CSF and SCF. As expected, LSK cells showed rapid growth during the first few days in culture and stopped proliferating at day 7 (Figure 2D) with accompanying changes in morphology (data not shown). At several time points after initiation of differentiation the relative expression of Ezh2 was assessed by quantitative PCR. Ezh2 expression was rapidly down-regulated on differentiation (Figure 2E).
Overexpression of Ezh2 in HSCs
Primary transplantation. Because Ezh2 was expressed in freshly isolated HSCs, but down-regulated during differentiation-inducing cultures, we wished to assess the functional consequences of enforced overexpression of Ezh2 in HSCs. The 5-FU-treated BM cells (CD45.1) were transduced with the control or Ezh2 vector and 3 × 106 cells were transplanted without sorting for GFP+ cells, in lethally irradiated CD45.2 recipients. In this way, the relative contribution of Ezh2-overexpressing (CD45.1+GFP+) stem cells compared to identically treated, but nontransduced stem cells (CD45.1+GFP-) could be determined. Shortly after transplantation, peripheral WBC counts were somewhat higher (P = .15) in mice given Ezh2 BM cells, but returned to normal levels and remained similar to control values as from 100 days after transplantation (Figure 3A). After transplantation, chimerism was assessed by quantifying the percentage of GFP+ cells within the donor CD45.1+ cell fraction (Figure 3B). No significant skewing in contribution of transduced (CD45.1+GFP+) cells was observed.
Secondary transplantation. Primary recipients were killed 4 months after transplantation (at a time point when no skewing of CD45.1+GFP+ contribution was observed). At this time point protein levels of Ezh2 were analyzed in the spleen. As expected, Ezh2 was present at a higher level in mice given transplants with Ezh2-overexpressing cells compared to controls (Figure 3C). LSK GFP+ cells were isolated from primary recipients and assayed for CAFC activity. Overexpression of Ezh2 resulted in a marked increase in CAFC day 7 and 14 frequencies and a more subtle effect for later appearing subsets (Figure 3D).
Secondary recipients were given transplants with harvested donor BM cells alone or in competitive repopulation with different ratios of freshly isolated CD45.2+ BM cells. At multiple time points after transplantation the recipients were analyzed for presence of CD45.1+GFP+ cells. In the 2 independent experiments serially passaged retrovirally transduced control cells (CD45.1+GFP+) showed little to no contribution to the peripheral blood in these secondary recipients, whereas Ezh2-overexpressing stem cells contributed to significantly higher (P < .001) levels (Figure 3E). The CRI, comparing stem cell activity of transduced cells (CD45.1+GFP+) with nontransduced cells (CD45.1+GFP-), was calculated 3 months after secondary transplantation. The CRI is a measurement of the quality of reconstitution of the transplanted cells, where a CRI of 1 by definition indicates equal stem cell potential of the 2 populations. The CRI was about 12-fold higher (P < .02) in Ezh2-overexpressing stem cells, when these were compared with serially passaged, but nontransduced (CD45.1+GFP-) cells (Figure 3F). The CRI was also calculated for transduced cells (CD45.1+GFP+) compared to freshly isolated CD45.2+ BM cells (Figure 3G). A 10-fold decrease of stem cell quality (ie, a CRI of 0.1) can usually be expected after a single serial transplantation1 and was indeed observed for control cells. Strikingly, overexpression of Ezh2 in HSCs completely prevented loss of stem cell quality after serial transplantation, as the CRI levels were maintained at a value of 1 (Figure 3G). CRI measurements do not necessarily correlate with actual stem cell frequency. To quantify stem cell numbers, we used limiting dilution analysis of chimerism data of all secondary recipients. Animals that showed more than 1% donor cell contribution in both myeloid and lymphoid compartments were considered engrafted. The BM of secondary recipients given transplants with Ezh2-overexpressing stem cells contained an 8- to 14-fold higher number of stem cells with long-term repopulating ability, showing that overexpression of Ezh2 results in an increase of the stem cell pool (Figure 3H).
Tertiary transplantation. Four months after secondary transplantation, recipients were killed and BM cells were once more serially transplanted in a third cohort of lethally irradiated recipients with addition of different ratios of freshly isolated CD45.2+ competitor BM cells. Figure 4A shows that virtually all myeloid cells in the peripheral blood were derived from Ezh2-overexpressing HSCs, whereas approximately half of the lymphoid lineage was donor-derived. Similar to primary and secondary recipients, peripheral blood cell values were normal in tertiary recipients of Ezh2-overexpressing stem cells (data not shown). The chimerism levels of all individual recipients of control and Ezh2-overexpressing stem cells 3 months after transplantation are shown in Figure 4B. Virtually no donor contribution was detected in recipients receiving control cells, whereas very high chimerism levels originated from Ezh2-overexpressing stem cells. In agreement, recipients that received serially transplanted control cells without competitor cells all died within 40 days after transplantation, whereas Ezh2-overexpressing stem cells provided radioprotection and long-term repopulation to tertiary recipients (Figure 4C). Similarly as after the secondary transplantation, the CRI was calculated for CD45.1+GFP+ cells over CD45.1+GFP- cells. Figure 4D shows these values for both secondary and tertiary recipients and documents that stem cells transduced with Ezh2 show a 10- and 25-fold increase in stem cell quality, respectively. Remarkably, when the quality of stem cells that were transduced with Ezh2 was compared with freshly isolated BM cells, no loss of stem cell quality could be demonstrated (Figure 4E).
Stem cell gene networks involving Ezh2
We have previously shown that different inbred strains of mice vary substantially with respect to numerous stem cell traits, including variation in stem cell frequency.48 In a panel of 30 BXD recombinant inbred strains, offspring of DBA/2J mice (containing high numbers of stem cells) and C57BL/6 mice (containing low numbers of stem cells), genomic loci that contribute to these traits were mapped. We identified a quantitative trait locus (QTL) on chromosome 18 that was associated with variation in stem cell pool size (Figure 5A).48 Recently, we have also measured genome-wide variation in gene expression levels in HSCs isolated from these BXD strains.41 This genetic genomic analysis has enabled us to identify clusters of genes that are transcriptionally controlled by a common locus. These data have been made accessible at GeneNetwork,40 an online genetic database that allows the query of coregulated transcripts for stem cell genes of interest. In the present study we determined that Ezh2 was highly, but differentially, expressed in HSCs of the various BXD strains (Figure 5B). Strikingly, variation in Ezh2 levels was largely accounted for by the exact chromosome 18 locus that we had previously associated with variation in stem cell frequency (Figure 5C). The gene encoding Ezh2 is located on chromosome 6 in mice, so variation in Ezh2 expression is regulated in trans by an as yet unidentified gene on chromosome 18. Further analysis showed that mice that had inherited a C57BL/6 allele at D18Mit83 (the marker most strongly associated with variation in Ezh2 levels) had significantly higher expression of Ezh2 (P < .003) than mice carrying a DBA/2J allele (Figure 5D). The limited genetic resolution precludes direct identification of the chromosome 18 gene that is responsible for the variation in Ezh2 gene expression. However, Ezh2 coregulated transcripts in stem cells can easily be retrieved from GeneNetwork.40 The 100 genes of which the transcript levels correlated best with variation in Ezh2 expression were identified. These genes were classified according to their functional annotation using WebGestalt (http://genereg.ornl.gov/webgestalt). We identified 3 functional categories for which transcripts correlating with variation in Ezh2 levels were significantly enriched (Figure 6A). A high number of coregulated stem cell genes was involved in DNA modification, DNA packaging, and DNA replication (Table 1). Two other categories of transcripts that were significantly more often observed than expected were associated with cell cycle/cytokinesis and, surprisingly, nucleocytoplasmic transport (Table 1).
Symbol . | Name . | QTL chromosome . | Mb, if on chromosome 18 . | Correlation . |
---|---|---|---|---|
DNA modification, packaging, or replication | ||||
Mcm4 | Minichromosome maintenance deficient 4 homolog (Saccharomyces cerevisiae) | 14, 18 | 20-23 | 0.7545 |
Mcm3 | Minichromosome maintenance deficient 3 (S cerevisiae) | 9, 18 | 20-23 | 0.7284 |
Pcna | Proliferating cell nuclear antigen | 18 | 20-23 | 0.8141 |
Prim1 | DNA primase, p49 subunit | 18 | 12-14, 20-25 | 0.7986 |
Hat1 | Histone aminotransferase 1 | 11 | — | 0.7973 |
Recc1 | Replication factor C 1 | 2 | — | 0.7911 |
Eed | Embryonic ectoderm development | 18 | 47-50 | 0.7555 |
Pole2 | Polymerase (DNA directed), ϵ 2 (p59 subunit) | 18 | 20-23 | 0.751 |
Mcm2 | Minichromosome maintenance deficient 2 mitotin (S cerevisiae) | 14, 18 | 20-23 | 0.739 |
Rrm1 | Ribonucleotide reductase M1 | 4 | — | 0.7359 |
Nap1/1 | Nucleosome assembly protein 1-like 1 | 18 | 47-50 | 0.7294 |
Rpa3 | Replication protein A3 | 18 | 17-20, 35-38 | 0.7239 |
H2afz | H2A histone family, member Z | 18 | 12-14 | 0.7115 |
Setdb1 | SET domain, bifurcated 1 | 6, 17 | — | 0.7102 |
Ahcy | S-adenosylhomocysteine hydrolase | 11 | — | 0.6855 |
Dntt | Deoxynucleotidyltransferase, terminal | 15 | — | 0.6851 |
Aurkb | Aurora kinase B | — | — | 0.718 |
Nuclear organization and biogenesis, cell cycle, and cytokinesis | ||||
Top2a | Topoisomerase (DNA) II α | 18 | 20-23 | 0.8733 |
Smc2/1 | SMC2 structural maintenance of chromosomes 2-like 1 (yeast) | — | — | 0.8507 |
Topbp1 | Topoisomerase (DNA) II β-binding protein | 18 | 20-25 | 0.8497 |
Mki67 | Antigen identified by monoclonal antibody Ki 67 | 2, 18 | 20-23 | 0.837 |
Tfdp1 | Transcription factor Dp 1 | 7 | — | 0.8299 |
Pole2 | Polymerase (DNA directed), ϵ 2 (p59 subunit) | 18 | 20-23 | 0.7986 |
Ccnb2 | Cyclin B2 | 13 | — | 0.7949 |
Recc1 | Replication factor C 1 | 2 | — | 0.7911 |
Ect2 | Ect2 oncogene | 4 | — | 0.7879 |
Gmnn | Geminin | 18 | 14-23 | 0.7865 |
Cdca5 | Cell division cycle associated 5 | 2 | — | 0.764 |
Prc1 | Protein regulator of cytokinesis 1 | 5, 7 | — | 0.7584 |
Nusap1 | Nucleolar and spindle-associated protein 1 | 14, 18 | 35-38, 71-73 | 0.7584 |
Mcm4 | Minichromosome maintenance deficient 4 homolog (S cerevisiae) | 14, 18 | 20-23 | 0.7545 |
Ppp1cc | Protein phosphatase 1, catalytic subunit, γ isoform | 18 | 50-53 | 0.751 |
Nek2 | NIMA (never in mitosis gene a)-related expressed kinase 2 | 7 | — | 0.751 |
Mad2/1 | MAD2 (mitotic arrest deficient, homolog)-like 1 (yeast) | 18 | 12-13, 35-38 | 0.7428 |
Mcm2 | Minichromosome maintenance deficient 2 mitotin (S cerevisiae) | 14, 18 | 20-23 | 0.739 |
Rrm1 | Ribonucleotide reductase M1 | 4 | — | 0.7359 |
Tacc3 | Transforming, acidic coiled-coil-containing protein 3 | 18 | 12-14 | 0.7335 |
Prim1 | DNA primase, p49 subunit | 18 | 12-14, 20-25 | 0.7321 |
Pcna | Proliferating cell nuclear antigen | 18 | 20-23 | 0.7321 |
Nap1/1 | Nucleosome assembly protein 1-like 1 | 18 | 47-50 | 0.7294 |
Mcm3 | Minichromosome maintenance deficient 3 (S cerevisiae) | 9, 18 | 20-23 | 0.7284 |
Rbl1 | Retinoblastoma-like 1 (p107) | 14 | — | 0.7268 |
Rpa3 | Replication protein A3 | 18 | 17-20, 35-38 | 0.7239 |
Aurkb | Aurora kinase B | — | — | 0.718 |
H2afz | H2A histone family, member Z | 18 | 12-14 | 0.7115 |
Setdb1 | SET domain, bifurcated 1 | 6, 17 | — | 0.7102 |
Stk6 | Serine/threonine kinase 6 (Aurka) | 9 | — | 0.7062 |
Dntt | Deoxynucleotidyltransferase, terminal | 15 | — | 0.6851 |
Cdca7 | Cell division cycle-associated 7 | 9, 18 | 14-23 | 0.6843 |
Cul2 | Cullin 2 | 18 | 50-54 | 0.682 |
Smc4/1 | SMC4 structural maintenance of chromosomes 4-like 1 (yeast) | 3, 18 | 12-14, 25-27 | — |
AA545217 | Expressed sequence AA545217 | 16 | — | — |
Protein import, nuclear transport, and nucleocytoplasmic transport | ||||
Kpna2 | Karyopherin (importin) α 2 | 18 | 51-53 | 0.78 |
Tacc3 | Transforming, acidic coiled-coil-containing protein 3 | 18 | 12-14 | 0.7335 |
Kpnb3 | Karyopherin (importin) β 3 | 5, 8, 13 | — | 0.7277 |
Xpo1 | Exportin 1, CRM1 homolog (yeast) | 14 | — | 0.7581 |
Symbol . | Name . | QTL chromosome . | Mb, if on chromosome 18 . | Correlation . |
---|---|---|---|---|
DNA modification, packaging, or replication | ||||
Mcm4 | Minichromosome maintenance deficient 4 homolog (Saccharomyces cerevisiae) | 14, 18 | 20-23 | 0.7545 |
Mcm3 | Minichromosome maintenance deficient 3 (S cerevisiae) | 9, 18 | 20-23 | 0.7284 |
Pcna | Proliferating cell nuclear antigen | 18 | 20-23 | 0.8141 |
Prim1 | DNA primase, p49 subunit | 18 | 12-14, 20-25 | 0.7986 |
Hat1 | Histone aminotransferase 1 | 11 | — | 0.7973 |
Recc1 | Replication factor C 1 | 2 | — | 0.7911 |
Eed | Embryonic ectoderm development | 18 | 47-50 | 0.7555 |
Pole2 | Polymerase (DNA directed), ϵ 2 (p59 subunit) | 18 | 20-23 | 0.751 |
Mcm2 | Minichromosome maintenance deficient 2 mitotin (S cerevisiae) | 14, 18 | 20-23 | 0.739 |
Rrm1 | Ribonucleotide reductase M1 | 4 | — | 0.7359 |
Nap1/1 | Nucleosome assembly protein 1-like 1 | 18 | 47-50 | 0.7294 |
Rpa3 | Replication protein A3 | 18 | 17-20, 35-38 | 0.7239 |
H2afz | H2A histone family, member Z | 18 | 12-14 | 0.7115 |
Setdb1 | SET domain, bifurcated 1 | 6, 17 | — | 0.7102 |
Ahcy | S-adenosylhomocysteine hydrolase | 11 | — | 0.6855 |
Dntt | Deoxynucleotidyltransferase, terminal | 15 | — | 0.6851 |
Aurkb | Aurora kinase B | — | — | 0.718 |
Nuclear organization and biogenesis, cell cycle, and cytokinesis | ||||
Top2a | Topoisomerase (DNA) II α | 18 | 20-23 | 0.8733 |
Smc2/1 | SMC2 structural maintenance of chromosomes 2-like 1 (yeast) | — | — | 0.8507 |
Topbp1 | Topoisomerase (DNA) II β-binding protein | 18 | 20-25 | 0.8497 |
Mki67 | Antigen identified by monoclonal antibody Ki 67 | 2, 18 | 20-23 | 0.837 |
Tfdp1 | Transcription factor Dp 1 | 7 | — | 0.8299 |
Pole2 | Polymerase (DNA directed), ϵ 2 (p59 subunit) | 18 | 20-23 | 0.7986 |
Ccnb2 | Cyclin B2 | 13 | — | 0.7949 |
Recc1 | Replication factor C 1 | 2 | — | 0.7911 |
Ect2 | Ect2 oncogene | 4 | — | 0.7879 |
Gmnn | Geminin | 18 | 14-23 | 0.7865 |
Cdca5 | Cell division cycle associated 5 | 2 | — | 0.764 |
Prc1 | Protein regulator of cytokinesis 1 | 5, 7 | — | 0.7584 |
Nusap1 | Nucleolar and spindle-associated protein 1 | 14, 18 | 35-38, 71-73 | 0.7584 |
Mcm4 | Minichromosome maintenance deficient 4 homolog (S cerevisiae) | 14, 18 | 20-23 | 0.7545 |
Ppp1cc | Protein phosphatase 1, catalytic subunit, γ isoform | 18 | 50-53 | 0.751 |
Nek2 | NIMA (never in mitosis gene a)-related expressed kinase 2 | 7 | — | 0.751 |
Mad2/1 | MAD2 (mitotic arrest deficient, homolog)-like 1 (yeast) | 18 | 12-13, 35-38 | 0.7428 |
Mcm2 | Minichromosome maintenance deficient 2 mitotin (S cerevisiae) | 14, 18 | 20-23 | 0.739 |
Rrm1 | Ribonucleotide reductase M1 | 4 | — | 0.7359 |
Tacc3 | Transforming, acidic coiled-coil-containing protein 3 | 18 | 12-14 | 0.7335 |
Prim1 | DNA primase, p49 subunit | 18 | 12-14, 20-25 | 0.7321 |
Pcna | Proliferating cell nuclear antigen | 18 | 20-23 | 0.7321 |
Nap1/1 | Nucleosome assembly protein 1-like 1 | 18 | 47-50 | 0.7294 |
Mcm3 | Minichromosome maintenance deficient 3 (S cerevisiae) | 9, 18 | 20-23 | 0.7284 |
Rbl1 | Retinoblastoma-like 1 (p107) | 14 | — | 0.7268 |
Rpa3 | Replication protein A3 | 18 | 17-20, 35-38 | 0.7239 |
Aurkb | Aurora kinase B | — | — | 0.718 |
H2afz | H2A histone family, member Z | 18 | 12-14 | 0.7115 |
Setdb1 | SET domain, bifurcated 1 | 6, 17 | — | 0.7102 |
Stk6 | Serine/threonine kinase 6 (Aurka) | 9 | — | 0.7062 |
Dntt | Deoxynucleotidyltransferase, terminal | 15 | — | 0.6851 |
Cdca7 | Cell division cycle-associated 7 | 9, 18 | 14-23 | 0.6843 |
Cul2 | Cullin 2 | 18 | 50-54 | 0.682 |
Smc4/1 | SMC4 structural maintenance of chromosomes 4-like 1 (yeast) | 3, 18 | 12-14, 25-27 | — |
AA545217 | Expressed sequence AA545217 | 16 | — | — |
Protein import, nuclear transport, and nucleocytoplasmic transport | ||||
Kpna2 | Karyopherin (importin) α 2 | 18 | 51-53 | 0.78 |
Tacc3 | Transforming, acidic coiled-coil-containing protein 3 | 18 | 12-14 | 0.7335 |
Kpnb3 | Karyopherin (importin) β 3 | 5, 8, 13 | — | 0.7277 |
Xpo1 | Exportin 1, CRM1 homolog (yeast) | 14 | — | 0.7581 |
Ezh2 coregulated transcripts in stem cells were retrieved from www.genenetwork.org. The 100 genes of which the transcript levels correlated best with variation in Ezh2 expression were identified. These genes were classified according to their functional annotation using Webgestalt (http://genereg.oml.gov/webgestalt). Three functional categories were identified. — indicates not applicable. Correlation values refer to Pearson correlation coefficients between variation in the gene of interest and Ezh2 levels (Figure 6B-C).
Interestingly, one of the stem cell transcripts whose levels correlated with variable Ezh2 expression was Eed, its well-known partner (Figure 6B). Other interesting novel putative Ezh2 partners or targets included Mcm2, Mcm3, and Mcm4; Aurka and Aurkb; Gmnn (Figure 6C); Prim1; topoisomerases; and polymerases. A substantial number of these transcripts was also controlled by the chromosome 18 locus (Table 1).
Discussion
In this study we show that Ezh2 is down-regulated on serial passaging of MEFs and during differentiation of primitive, highly proliferative HSCs to committed myeloid cells. These 2 complementary models suggest a general suppressive role for Ezh2 during cellular differentiation and senescence. Indeed, we show that overexpression of Ezh2 leads to bypassing of the senescence program in MEFs. More strikingly, enforced expression of Ezh2 in HSCs completely prevents stem cell exhaustion and stabilizes normal stem cell functioning, even after 3 serial transplantations.
It has been previously shown by others and us that exhaustion of the HSC pool is apparent following serial transplantation,2,49,50 chemotherapy,51,52 and also during normal aging.36,49 The aging of HSCs is genetically regulated, because different inbred mouse strains show distinct patterns,36,53 which are maintained when HSCs of different mouse strains are present in the same microenvironment.54,55 However, the underlying molecular mechanism responsible for a decline of HSC functioning during aging and serial transplantation remains to be elucidated. Several groups have identified genes that directly influence HSC self-renewal potential and their differentiation program, such as Hoxb4,35 Notch1,56 and Bmi-1.8 Our present data document that Ezh2 is able to completely preserve stem cell functioning after repeated serial transplantation, notably without resulting in malignant transformation. Interestingly, a recent gene profiling study in which expression patterns in young and aged HSCs were compared, documented that Ezh2 levels decrease during normal aging.57
The natural variation in Ezh2 transcript levels that can be observed in stem cells of recombinant inbred mouse strains (BXDs), allows 2 approaches.41 First, the genomic locus that regulates Ezh2 expression can be mapped and, second, coregulated transcripts can be identified. Using the first approach we show that variation in expression of Ezh2 is regulated by a locus on chromosome 18. Strikingly, exactly the same locus is also associated with HSC frequency.48 Because many transcripts are located in the critical interval, the Ezh2-regulating gene on chromosome 18 remains unknown at present. Using the second approach, genetic networks associated with variation in Ezh2 expression can be elucidated. We found that Eed expression levels were positively correlated with Ezh2 gene expression. Because Eed is a well-known partner of Ezh2,31 this finding proves the potential power of this approach. Furthermore, expression levels of Gmnn, a DNA replication inhibitor with which PRC1 members are associated,58 were also positively correlated with Ezh2 expression. We identified many Ezh2 coregulated stem cell genes that qualify as chromatin remodeling factors and DNA binding proteins. Interestingly, the interval on chromosome 18 that was associated with variation in Ezh2 expression levels also regulated the expression levels of 12 of 42 coregulated genes. Because this locus is involved in regulating HSC pool size,48 this strongly suggests that chromatin remodeling gene networks affecting epigenetic transcription programs are critically involved in stem cell self-renewal.
Ezh2 can act as a transcriptional repressor as a result of its HMT activity, which ultimately results in DNA condensation and a more compact chromatin structure. Overexpression of Ezh2 is therefore expected to change global chromatin structure and consequently alter gene expression profiles. Moreover, because methylation is considered to be relatively irreversible, these epigenetic changes of gene expression are stable. In stem cell homeostasis a delicate balance exists between self-renewal and terminal differentiation. It is likely that during replicative stress (resulting from serial transplantation, but also from normal aging) this balance weighs in favor of terminal differentiation, resulting in exhaustion of the HSC pool. Emerging evidence suggests that this balance is partly regulated by the histone code, which specifies gene transcription and provides cellular memory. It has been suggested that during normal aging methylation patterns are less strictly maintained, resulting in loss of heterochromatin and deregulated gene transcription.59 It is becoming evident that PcG genes play an important role here because aberrant expression of different PcG genes has been shown to change the program of stem cell self-renewal.
In this study we show that stem cell exhaustion is readily induced but can be prevented in HSCs by overexpressing Ezh2. These data suggest that after replicative stress of HSCs chromatin alterations occur, which impede stem cell functioning. In contrast, in HSCs that overexpress Ezh2 the chromatin structure is stabilized, which provides maintenance of HSC quality. Chemical modifiers of higher order chromatin structure may have a functional role in clinical ex vivo stem cell expansion protocols.
Prepublished online as Blood First Edition Paper, November 17, 2005; DOI 10.1182/blood-2005-09-3585.
Supported by a grant from the Dutch Organization of Scientific Research (NWO; grant 901-08-339) and the National Institutes of Health (NIH; grant RO1 HL073710).
L.M.K. designed and performed research, analyzed data, and wrote the paper; L.V.B. performed research, contributed vital new analytical tools, analyzed data, and wrote the paper; A.d.B., S.H., J.D., and E.W. performed research and analyzed data; B.D. performed research; and G.d.H. designed research, contributed vital new analytical tools, and wrote the paper.
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 U.S.C. section 1734.
The authors would like to thank G. Harms and A. van den Berg for their assistance with the mouse filter gene arrays; C. Jordan for providing us with the MIEV vector; A. Otte for providing us with the Ezh2 antibody; G. Mesander and H. Moes for flow cytometry support; and R. van Os for critically reading the manuscript and for advice on CRI calculations.
This feature is available to Subscribers Only
Sign In or Create an Account Close Modal