|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
TISSUE-SPECIFIC STEM CELLS |
aExperimental Haematology, Division of Cancer Sciences, University of Glasgow, Glasgow Royal Infirmary, Glasgow, United Kingdom;
bSection for Bioinformatics, Cancer Research U.K., Beatson Laboratories, Garscube Estate, Glasgow, United Kingdom;
cDivision of Immunology, Infection and Inflammation, University of Glasgow, Glasgow Biomedical Research Centre, Glasgow, United Kingdom
Key Words. Stem cells • Hemopoiesis • Leukemia • Chemokines
Correspondence: Gerard J. Graham, Ph.D., Division of Immunology, Infection and Inflammation, University of Glasgow, Glasgow Biomedical Research Centre, 120 University Place, Glasgow G12 8TA, United Kingdom. Telephone: 44-141-330-3982; Fax: 44-141-330-4297; e-mail: g.graham{at}clinmed.gla.ac.uk; Tessa L. Holyoake, Experimental Haematology, Division of Cancer Sciences, University of Glasgow, Glasgow Royal Infirmary, 10 Alexandra Parade, Glasgow G31 2ER, United Kingdom. Telephone: 44-141-211-4676; Fax: 44-141-211-0414; e-mail: Tlh1g{at}clinmed.gla.ac.uk
Received on April 4, 2007;
accepted for publication on August 14, 2007.
First published online in STEM CELLS EXPRESS August 23, 2007.
| ABSTRACT |
|---|
|
|
|---|
Disclosure of potential conflicts of interest is found at the end of this article.
| INTRODUCTION |
|---|
|
|
|---|
In contrast, chronic myeloid leukemic (CML) stem cells are constitutively present in the circulation. We have also recently shown the presence of a significant quiescent (G0) population within the CD34+ population of hemopoietic cells from patients with CML. These G0, CD34+, CML cells are fully functional and are able to repopulate sublethally irradiated NOD/SCID mice [10]. This quiescent pool is likely to be therapeutically relevant as, in contrast to proliferating CD34+ CML cells, it appears to be resistant to the BCR-ABL inhibitor imatinib mesylate [11]. This suggests that these cells may be important mediators of post-therapy relapse in patients with CML [12–15].
To understand the molecular basis for the differences in migratory properties of normal and leukemic and G0, and dividing HSC, we have examined the transcriptome of human CD34+ peripheral blood (PB) cells sorted for cell cycle position. To do this, we have isolated G0 and dividing CD34+ hemopoietic cell populations from normal and patient samples using Hoechst (HST) and pyronin Y (Py) staining [10–12, 15]. This has allowed us to define distinct gene expression profiles for CD34+ cells at specific cell cycle positions and to contrast the transcriptional profile of normal and CML cells in both proliferative states. This study represents the first assessment of the transcriptome of primary quiescent and proliferating CD34+ hemopoietic cells and of its alteration in leukemia, and we believe that the data set we report will be of value in advancing our understanding of the transcriptional profile of primitive hemopoietic cells. Furthermore, our data highlight differences in the expression of genes reflecting the proliferative activity of normal and CML CD34+ cells and identify chemokines as the major class of genes that is differentially expressed between G0 and dividing CD34+ hemopoietic cells, both normal and leukemic.
| MATERIALS AND METHODS |
|---|
|
|
|---|
CD34+ Cell Selection
CD34+ cells were enriched using the Isolex 50 system (NEXELL International, Brussels, Belgium, http://www.nexellinc.com) and cryopreserved in 10% dimethyl sulfoxide (Sigma-Aldrich, Poole, U.K., http://www.sigmaaldrich.com) in ALBA (4.5% human albumin solution; Scottish National Blood Transfusion Service). Flow cytometry analysis demonstrated that the enriched CML cells used for Hoechst/pyronin Y sorting were between 85% and 98% CD34+ and the normal cells between 76% and 85% CD34+. The CD34+ CML cells were also shown to be between 93% and 100% BCR-ABL-positive by fluorescence in situ hybridization.
Cell Culture
Cells were recovered from liquid nitrogen and washed once in phosphate-buffered saline (PBS)/2% fetal calf serum (FCS; Life Technologies, Paisley, U.K., http://www.lifetech.com) (PBS/2%). Cells were cultured overnight in Iscove's modified Dulbecco's medium (Sigma-Aldrich) supplemented with bovine serum albumin, insulin, and transferrin serum substitute (StemCell Technologies, Vancouver, BC, Canada, http://www.stemcell.com), 40 µg/ml low-density lipoprotein, and 10–4 M 2-mercaptoethanol. Serum-free medium (SFM) was further supplemented with a five-growth factor cocktail containing recombinant human, Flt3 ligand, and stem cell factor (100 ng/ml), interleukin (IL) 3, IL6 (StemCell Technologies), and G-CSF (Chugai Pharmaceutical, London, http://www.chugai-pharm.co.jp/hc/chugai_top_en.jsp), all at 20 ng/ml.
Flow Cytometry and Cell Sorting
Cells were recovered following overnight culture, washed in PBS/2%, and stained with 10 µM HST (Cambridge Bioscience, London, http://www.bioscience.co.uk) for 90 minutes at 37°C. Py (5 µg/ml) (Sigma-Aldrich) was added and incubated at 37°C for 45 minutes. Anti-CD34-fluorescein isothiocyanate (Becton, Dickinson and Company, Oxford, U.K., http://www.bd.com) was added for 15 minutes at room temperature before the cells were washed twice in PBS/2% containing HST and Py at equivalent concentrations. Prior to fluorescence-activated cell sorting (FACS) analysis, 1 µg/ml propidium iodide (PI) (Sigma-Aldrich) was added to the final wash. Cells were resuspended in PBS/2% containing HST and Py and kept on ice in the dark until they were sorted using the FACSVantage cell sorter (Becton Dickinson). Cell populations were sorted by gating first on the PI–/CD34+ cells and then on the HSTlo/Pylo (G0) and HST+/Py+ (G1/S/G2/M) cells. For further analysis of cell division, sorted cells were stained with 1 µM carboxyfluorescein diacetate succinimidyl diester (CFSE) (Cambridge Bioscience) in PBS/2% for 10 minutes at 37°C before being washed in ice-cold PBS/20% FCS and then in PBS/2%. Cells were then cultured for up to 3 days in supplemented SFM as described. The cells were recovered and stained with CD34-PE (Becton Dickinson) and PI before FACS analysis.
Isolation of Total RNA
Cells were placed in Trizol (Invitrogen, Paisley, U.K., http://www.invitrogen.com) and stored at –80°C. Total RNA was extracted from all samples according to the manufacturer's protocol. Contaminating genomic DNA was removed using the DNA-free kit (Ambion, Austin, TX, http://www.ambion.com).
Real-Time Polymerase Chain Reaction
TaqMan (Applied Biosystems, Warrington, U.K., http://www.appliedbiosystems.com) real-time polymerase chain reaction (PCR) was used for validation of microarray data. cDNA from samples not used for microarrays (processed using the same protocol) was applied to custom microfluidics cards containing probe/primer sets for the genes of interest. The cards were run on an ABI 7900 instrument (Applied Biosystems) using standard settings. GAPDH was included on the card as a housekeeping gene, and analysis was carried out using the relative quantitation 
–2 cycle threshhold method [16].
Microarrays
RNA for application to Affymetrix gene chips (Affymetrix UK, High Wycombe, U.K., http://www.affymetrix.com) was processed as described. The quality and quantity of RNA was assessed using the BioAnalyzer 2100 (Agilent Technologies, Cheshire, U.K., http://www.agilent.com) and the NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, DE, http://www.nanodrop.com). The specific gene chips used were Affymetrix U133A human gene chips, which are capable of measuring the expression of 14,500 well-characterized human genes.
Thirty-five nanograms of each sample was subject to linear amplification using the Two Cycle Target Labeling Kit (Affymetrix). Processing was carried out at the Sir Henry Wellcome Functional Genomics Facility (SHWFGF) (University of Glasgow, Glasgow, U.K.).
Microarray Analysis
Raw data were normalized using robust multichip average and used to produce lists of differentially expressed genes using rank products (RP) [17, 18]. RP involves the production of lists of up- and downregulated genes and assigning a statistical measure to these differences according to the consistency of the differential expression levels across replicates. Importantly, for single-channel arrays such as the Affymetrix gene chips used in the present study, all possible pairwise comparisons of gene expression changes are calculated, and thus the statistical analysis implicitly reflects the consistency of any reported changes across replicate chips. The gene lists prepared using RP are then processed using iterative group analysis (IGA) to identify statistically significantly changed groups of genes [19]. IGA assigns differentially expressed genes to gene groups defined according to Gene Ontology (http://www.geneontology.org) and measures the statistical strength of these changes in gene class according to both the numbers of differentially expressed genes assigned to each class and the relative position of the genes in the RP lists. Importantly, by focusing on groups of genes instead of single genes, one can use group members as internal replicates to determine statistical significance. An additional advantage of the IGA statistical process is that it does not require all members of a gene grouping to change simultaneously or even in the same direction and is thus less restrictive than other Gene Ontology-based methods. All of the above procedures were done in the SHWFGF at the University of Glasgow and comprise the routine array data analysis pipeline of that facility (http://www.gla.ac.uk/functionalgenomics/rp/affy_analysis.html).
GeneSpring software (Agilent Technologies) was used for the generation of the gene lists reported in supplemental online Tables 2–7. Lists of genes with least variance were produced using the R package (http://www.r-project.org). Other statistics were generated using Mann-Whitney tests or analysis of variance. An arbitrary cut off of threefold differential expression was selected for the presentation of the data generated from the GeneSpring-based analysis. The reason for the choice of threefold was that reducing the fold expression cutoff markedly increased the numbers of genes, in some cases to unmanageable levels. Importantly, the IGA analysis reported in Table 2 is independent of assignment of arbitrary cutoff levels and relates to the statistical robustness of the differences in gene expression as assessed using RP.
| RESULTS |
|---|
|
|
|---|
|
Comparison of the Affymetrix data identified 860 genes as being similarly and consistently expressed across all four groups (normal and CML, G0 and dividing). In addition to genes involved in general metabolic pathways, this cohort contained numerous genes associated with primitive hemopoietic cell function and self-renewal. For example, seven members of the Notch/Wnt pathways were present, as well as a number of other primitive cell-associated genes, including GATA2, RUNX1, HOXC4, ITAG6, ALDH, and E2F4 (data not shown). Many other regulators of transcription, signal transduction, and chromatin modification, all previously associated with a primitive hemopoietic cell phenotype, were also present. Thus, a common primitive hemopoietic cell transcriptional signature was evident in each of the four sorted populations.
Analysis of Differentially Expressed Genes
Table 1 shows the number of genes changed more than threefold between the sample groups; individual gene lists are presented in supplemental online Tables 2–7. The greatest numbers of genes changed were in the normal G0 versus normal dividing (188) and normal G0 versus CML dividing (332) groups, with only low numbers being changed between CML G0 and dividing cells (37), suggesting that CML G0 cells are closer to their cycling counterparts than normal G0 cells are to theirs. Interestingly, comparison of CML versus normal G0 cells revealed a high degree of transcriptional variability, again indicating a fundamental difference between quiescent normal and CML cells. In contrast, normal and CML dividing cells were transcriptionally very similar. Euclidean distance modeling of the relatedness of the sample groups (Fig. 2) confirmed that normal and CML dividing cells were closest in gene expression, with CML G0 cells transcriptionally closer to dividing cells than to normal G0 cells. These general conclusions were also borne out by IGA [19], which identified only limited numbers of gene groupings as being significantly changed between dividing normal and CML cells (Table 2).
|
|
|
|
Classification of Other Differentially Expressed Genes
G0 Versus Dividing Cells. In addition to differences in cell cycle gene expression, normal dividing CD34+ cells displayed higher levels of transcripts for genes involved in DNA repair and oxygen transport than their quiescent counterparts (Table 2). In contrast, dividing CD34+ CML cells displayed very few significant transcriptional differences versus G0 cells other than cell cycle genes (Tables 2, 3; supplemental online Tables 2–7).
Much more similarity was seen in the reciprocal comparison of genes upregulated in G0 normal and CML cells compared with their dividing counterparts. Normal G0 cells displayed increased levels of expression of many histone genes, as well as genes involved in regulation of cell death and antigen presentation, whereas CML G0 cells expressed increased levels of erythroid genes and the aldo-keto reductase genes AKR1C1 and AKR1C2.
Strikingly, the family of genes exhibiting the most significant upregulation, in both normal and CML G0 cells compared with their dividing counterparts, was the chemokine family [20] (Fig. 3). The differentially expressed chemokines were predominantly inflammatory CXC chemokines (CXCR2 and CXCR3 ligands; Fig. 3). Interestingly, among the chemokines, the CXCR2 ligands (CXCLs 1, 2, 3, and 6) were most strongly differentially expressed. These results, which we obtained from GeneSpring analysis, were confirmed by IGA analysis, which identified additional group members and ranked chemokines statistically as the top changed group in the G0 versus dividing cells for both normal and CML samples (Table 2; supplemental online Tables 2–7).
|
CML Versus Normal CD34+ Cells. Aside from cell cycle genes, IGA analysis did not identify any other significantly altered classes of genes as being upregulated in dividing CML compared with normal cells. GeneSpring analysis highlighted only four genes as being more highly expressed by the dividing CML cells than their normal counterparts (Table 2). In examining the genes that were expressed at lower levels in CML, IGA analysis identified genes involved in antigen presentation and processing as being the most significantly altered. In addition, the dividing CML cells expressed lower levels of chemokines than normal cells, although the differences in expression levels were low. GeneSpring analysis identified 12 genes as being expressed at level more than threefold lower in dividing CML compared with normal cells. Again, many of these were genes associated with antigen presentation and processing.
In contrast to dividing cells, the transcriptional profiles of CML and normal G0 cells were more different (Fig. 2). IGA analysis (Table 2) identified significantly higher expression of genes associated with active biological processes in CML compared with normal G0 cells. The top changed groups were hemoglobin complex and oxygen transport. GeneSpring analysis identified 83 genes expressed at a level more than threefold higher in CML compared with normal G0 cells. These genes included many involved in erythroid cell function. IGA analyses again revealed genes involved in antigen processing and presentation, as well as chemokine genes, as being expressed at lower levels in CML G0 compared with normal G0 cells. GeneSpring analysis highlighted 55 genes as being expressed at lower levels by G0 CML compared with G0 normal cells. These again included genes involved in antigen processing and presentation but also included genes such as the HSC marker AC133 [21].
Validation of Microarray Results
Given the difficulty in accessing normal samples and the very poor yields of RNA from such samples, these validation assays were performed on PB cells from three independent CML patients. Initial attempts at measuring chemokines released by these cells using enzyme-linked immunosorbent assay or Luminex approaches were unsuccessful because of the very small numbers of sorted cells. Thus, validation of the array analyses was confirmed by quantitative PCR using microfluidic cards. Expression of genes identified as being differentially expressed in CML G0 cells compared with CML dividing cells was examined using RNA from dividing and G0 sorted CD34+ CML cell populations. RNA from unsorted CD34+ CML cells was used as a calibrator and expression in G0 or dividing cells is reported relative to the expression levels found in these unsorted cells. Assessment of CD34 expression was included in this analysis, and as expected, similar expression levels were detected in the unsorted G0 and dividing cell populations (Fig. 4). Additional results (Fig. 4) confirmed the array analyses, showing differential expression of the chemokines (CCL19, CXCL13, CXCL12, CXCL3, and CXCL5) and the receptor CXCR4 in CML G0 compared with dividing cells. These quantitative PCR results therefore validate the differences in gene expression revealed by the microarray study.
|
To further examine the differences in division kinetics of CML and normal CD34+ cells, they were sorted into G0 and dividing populations, stained with CFSE, and cultured in growth factor-supplemented medium. The cells were then assessed for cell division at 36 hours and again at 3 days. These time points were chosen because we have previously shown substantial entry of quiescent CML cells into cell cycle following 3 days of exposure to cytokines [10, 11]. Thirty-six hours therefore represents a time point intermediate between the initiation and the establishment of this entry into cell cycle. Figure 5 shows that for CD34+ CML cells at 36 hours, the G0 cells had undergone up to three divisions, with only 14% remaining undivided. By 3 days, almost all of the CML G0 cells had gone into division, with only 0.9% remaining in the undivided peak. In contrast, for normal cells, at 36 hours 62% of the G0 cells remained undivided, and 54% remained undivided by 3 days. The dividing CML cells also underwent up to three divisions by 36 hours, with only 1.8% remaining undivided. By 3 days, only 0.2% from the dividing sample remained undivided. In contrast, 46% of the dividing normal cells had not completed a further cell division by 36 hours, and after 3 days, 21% of the dividing cells had not divided further. These observations are consistent with the gene data, suggesting that CML cells demonstrate an increased tendency to enter cell division.
|
| DISCUSSION |
|---|
|
|
|---|
Examination of the transcriptional data that we generated provides a number of insights into the functional differences between normal and CML quiescent and proliferating cells. For example, analysis of the expression of genes involved in cell cycle regulation indicates differences in proliferative responses of normal and CML cells to overnight cytokine stimulation. Specifically, normal CD34+ cells express a cohort of cell cycle genes indicative of early positioning within cell cycle. In contrast, CML cells express later cell cycle genes, suggestive of a more rapid transit through cell cycle following growth factor stimulation. We examined the functional relevance of this, and we show that both G0 and dividing CML cells proliferate more rapidly following cytokine stimulation than their normal counterparts. These transcriptional and biological data are therefore highly suggestive of an intrinsic, BCR-ABL-dependent enhancement of cell cycle transit of CML compared with normal CD34+ cells. Although our findings in this regard are similar to those of other groups [23, 24], they differ significantly from the results from Buckle et al. [22]. An important difference between our experimental approach and that of Buckle et al. [22] is that we have examined the proliferative response of the normal and CML-derived CD34+ cells to growth factor stimulation and not simply in unstimulated cells. Thus, we believe that our data provide a more accurate assessment of the proliferative potential of CML versus normal CD34+ cells. Together with the other previous studies [23, 24], it therefore appears that CD34+ cells from CML sources have a higher intrinsic proliferative activity than their normal counterparts.
In addition to cell cycle genes, other gene differences were highlighted by both IGA and GeneSpring analyses. The implications of the enhanced erythroid gene expression levels in G0 CML cells are not clear at present, although they may relate to the predisposition of HSC to erythroid differentiation (demonstrated in previous studies on the hemopoietic differentiation potential of murine embryonic stem cells [41]), which appears, therefore, to be exaggerated in G0 CML CD34+ cells. G0 CML cells also expressed higher levels of AKR1C1 and AKR1C2 than their proliferating counterparts. It has been reported that the AKR1C enzyme is present in myeloid leukemic cell lines and that its activity inhibits their differentiation [42]. We postulate that expression of these genes may indicate one mechanism that maintains a primitive phenotype in CML G0 CD34+ cells and prevents exhaustion of the malignant clone.
Of the most note were the marked differences in expression of chemokine genes, which were much more strongly expressed by G0 cells than by proliferating cells. The overexpressed chemokines, in general, were ligands for CXCR2 and CXCR3. The CXCR2 ligands were most significantly overexpressed by G0 CD34+ cells and with CXCR2 being the predominant neutrophil-expressed chemokine receptor, this suggests a relationship between G0 CD34+ cells and neutrophils. Notably, all the CD34+ cells used in this study were from PB and in the case of normal cells were the product of mobilization regimens. Recent data on HSC mobilization clearly implicate neutrophils as a key cell type in the mobilization process [43]. The likelihood is that neutrophil-derived proteases are required to disrupt local molecular adhesions anchoring HSC to stromal cells within BM niches [44–46]. One outstanding question relates to the mechanism of targeting neutrophils to the HSC niche during mobilization. As CXCR2 is the major neutrophil chemokine receptor, we postulate that chemokines expressed by quiescent HSC may be involved in neutrophil recruitment to the niche during mobilization. This would account for the preferential mobilization of G0 HSC and the requirement for CXCR2 in mobilization [47] and would ensure localization of neutrophil-protease release to the BM niche. This hypothesis is currently being tested in our laboratories.
| CONCLUSION |
|---|
|
|
|---|
| DISCLOSURE OF POTENTIAL CONFLICTS OF INTEREST |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. S. M. Yong, K. Keyvanfar, N. Hensel, R. Eniafe, B. N. Savani, M. Berg, A. Lundqvist, S. Adams, E. M. Sloand, J. M. Goldman, et al. Primitive quiescent CD34+ cells in chronic myeloid leukemia are targeted by in vitro expanded natural killer cells, which are functionally enhanced by bortezomib Blood, January 22, 2009; 113(4): 875 - 882. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. H. Jamieson Chronic Myeloid Leukemia Stem Cells Hematology, January 1, 2008; 2008(1): 436 - 442. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| STEM CELLS | THE ONCOLOGIST | CME | ALPHAMED PRESS JOURNALS |