Surface density of CD27 and CD11b subdivides mouse natural killer (NK) cells into 4 subsets: CD11blowCD27low, CD11blowCD27high, CD11bhighCD27high, and CD11bhighCD27low. To determine the developmental relationship between these 4 subsets, we used several complementary approaches. First, we took advantage of NDE transgenic mice that express enhanced green fluorescent protein (EGFP) and diphtheria toxin receptor specifically in NK cells. Diphtheria toxin injection leads to a transient depletion of NK cells, allowing the monitoring of the phenotype of developing EGFP+ NK cells after diphtheria toxin injection. Second, we evaluated the overall proximity between NK-cell subsets based on their global gene profile. Third, we compared the proliferative capacity of NK-cell subsets at steady state or during replenishment of the NK-cell pool. Fourth, we performed adoptive transfers of EGFP+ NK cell subsets from NDE mice into unirradiated mice and followed the fate of transferred cells. The results of these various experiments collectively support a 4-stage model of NK-cell maturation CD11blowCD27low → CD11blowCD27high → CD11bhighCD27high → CD11bhighCD27low. This developmental program appears to be associated with a progressive acquisition of NK-cell effector functions.
Introduction
Natural killer (NK) cells are lymphocytes of the innate immunity with effector functions, such as perforin-dependent cytotoxicity and interferon-γ secretion. Several lines of evidence implicate them in the early control of virus infection, in tumor immunosurveillance, and in the regulation of immune responses.1,2 Mature NK cells express a wide variety of cell-surface receptors that enable them to recognize targets expressing low surface amounts of major histocompatibility complex class I molecules or high surface amounts of molecules induced by stress or microbial molecules.3,4 The NK-cell receptors specific for the major histocompatibility complex class I molecules include the human killer-cell Ig-like (KIR) and mouse Ly49 receptors, and CD94/NKG2 complexes in both species.5
NK cells are heterogeneous as to their phenotype, function, and anatomic distribution. In human, 2 subsets of NK cells differing by their level of CD56 expression have been identified and well characterized.6 CD56bright NK cells are found mostly in lymphoid organs.7 They express few surface KIR molecules, are poorly cytolytic, but proliferate vigorously after interaction with activated dendritic cells. By contrast, CD56dim NK cells are predominant in peripheral blood, express high levels of KIR molecules, and display potent cytotoxic function. It has long been proposed that CD56bright NK cells were precursors of CD56dim NK cells.8 In mice, 2 populations of NK cells with low and high expression of CD11b have been identified. In fetal and neonatal mice, CD11blow NK cells composed most of the NK-cell population.9 In the adult, CD11blow NK cells were found mostly in the bone marrow, lymph nodes, and liver and had a high rate of homeostatic proliferation. CD11bhigh NK cells were present rather at peripheral sites, such as spleen, peripheral blood mononuclear cells (PBMCs), and lung. The latter expressed high levels of Ly49 receptors and displayed potent effector functions compared with the CD11blow NK cells.10 The CD11bhigh NK-cell subset could be further dissected in 2 populations based on CD27 expression: CD11bhighCD27high and CD11bhighCD27lowcells. CD11bhighCD27high NK cells were found to be the most potent effector cells and responded better to CXCR3 or CXCR4 ligands in vitro.11 CD11bhighCD27low NK cells expressed Ly49 receptors at the highest frequency and expressed unique receptors for chemoattractants, such as CX3CR112 or S1P5.13
It has been proposed that CD11blowCD27high, CD11bhighCD27high, and CD11bhighCD27low (referred to hereafter as CD11blow, DP [for double-positive], and CD27low NK cells, respectively) were discrete stages of in vivo maturation following the pathway CD11blow → DP → CD27low.11,14 Besides correlative evidence, the most compelling support for this model lies in adoptive transfer of sorted CD11blow NK cells into irradiated recipient mice, which suggested that CD11bhigh NK cells were derived from CD11blow NK cells.9,10,14,15 However, irradiation of recipient mice is a harsh conditioning regimen that can impact on transferred cells by inducing their “homeostatic” proliferation and differentiation.16 Thus, whether subsets defined by CD11b and CD27 expression represent distinct NK-cell development stages has not been formally demonstrated yet.
We sought to address the developmental relationship between NK-cell subsets using a combination of 4 complementary approaches assessing the kinetics of appearance of each subset during development, their gene expression profile at the pan-genomic level, their rate of proliferation, and their potential to give rise to the other subsets in vivo. Our results show that NK-cell maturation is a 4-stage process that starts at a CD11blowCD27low stage and leads them through the following stages: CD11blowCD27high → CD11bhighCD27high → CD11bhighCD27low.
Methods
Mice and NK-cell depletion
The generation of NDE mice has been described previously.17 More information on the characterization of these mice may be found on the CIML website (http://www.ciml.univ-mrs.fr/Lab/Vivier/Protocoles/NKDTR-EGFP_mice.pdf). These mice were 6 generations backcrossed to C57BL/6 mice. NDE mice with at least 80% enhanced green fluorescent protein (EGFP)-positive NK cells were selected for all the depletion experiments in this study. NK-cell depletion was obtained by intravenous injection of 2 μg per mouse of diphtheria toxin (DT; Calbiochem, La Jolla, CA). We verified that DT is short lived (24 hours) and does not induce adverse effects in NDE mice (data not shown). C57BL/6 mice and RAG-2 KO mice on C57BL/6 background were obtained from Charles River Laboratories (Les Oncins, France). All mice were used between 6 and 12 weeks of age unless specified. All experiments have been performed in accordance with protocols approved by the local ethical committee at Université de la Méditerranée.
Flow cytometry and lymphocyte preparation
Lymphocyte suspensions were prepared from different organs as described previously.17 For cell-surface staining, anti-CD16/CD32 (2.4G.2) hybridoma supernatant was added. NKp46 expression was detected using the rat 29A1.4 monoclonal antibody (mAb)17 coupled with Alexa 647 fluorochrome (Invitrogen, Carlsbad, CA). Pacific Blue anti–mouse CD3 complex was purchased from BioLegend (San Diego, CA). Phycoerythrin (PE)-Cy5 anti–mouse CD127, biotin anti–mouse NKG2D, and PE anti–mouse granzyme B were from eBioscience (San Diego, CA). All other antibodies were from BD Biosciences PharMingen (San Diego, CA). Streptavidin-peridinin chlorophyll protein from BD Biosciences PharMingen was used as a secondary reagent for biotinylated antibodies. Before intracellular staining, samples were treated with BD Cytofix/Cytoperm and Perm/Wash solutions (BD Biosciences, San Jose, CA). Samples were run on a fluorescence-activated cell sorter (FACS) Canto II or a FACS LSR II cytometers (BD Biosciences) and analyzed using FlowJo 8.7 software (TreeStar, Ashland, OR).
Cell sorting and transfer
Lymphocytes obtained from spleen and liver of NDE mice were incubated for 20 minutes at 4°C with a cocktail of mAb: rat anti–mouse CD4 (H129.19), CD5 (53-7.3), and CD8 (H59.101.2; all hybridoma supernatants), purified anti–mouse I-A/I-E, and purified rat anti–mouse TER-119 (both from BD Biosciences PharMingen). Cells were then washed and incubated with BioMag goat anti–rat IgG-coated magnetic beads (Polysciences, Warrington, PA) for 45 minutes at 4°C. Cells attached to beads were removed with a magnetic separator. A total of 80% pure NK cells were obtained using this procedure. Cells were then incubated with PE-conjugated anti-CD27 antibody and PE-Cy7–conjugated anti-CD11b antibody. NK cells were identified as EGFP+ and sorted into different subsets using a FACSAria Cell Sorter (BD Biosciences). Purity of sorted cell populations was more than 98% as checked by flow cytometry. For cell transfer, 8 × 103 to 2 × 105 sorted cells were injected intravenously into unirradiated wild-type littermate recipients. The presence of transferred cells was analyzed at days 7 and 14 after transfer by cytofluorimetric analysis of NK cell–enriched splenocytes. The expression of CD27 and CD11b was analyzed on CD19− CD3− NK1.1+ EGFP+ cells.
Measurement of BrdU incorporation
Mice were given 2 intraperitoneal injections of 0.5 mg bromodeoxyuridine (BrdU; BD Biosciences PharMingen) within a 2-hour interval. Six hours after the first BrdU injection, mice were killed and organs harvested. Cells derived from the bone marrow were stained with mAb as described in “Flow cytometry and lymphocyte preparation.” After fixation and permeabilization, cells were treated with DNase and stained with fluorescein isothiocyanate anti-BrdU mAb, according to the manufacturer's instructions. BrdU incorporation for different cell populations was measured by flow cytometry.
Gene profiling and data mining
NK-cell subsets were sorted by flow cytometry from RAG2−/− mice. We determined that the distribution of NK-cell subsets in RAG2−/− mice was similar to that in C57BL/6 mice. RNA was extracted from each sorted NK-cell subset with the QIAGEN (Courtaboeuf, France) micro-RNAeasy kit. Quality and absence of genomic DNA contamination were assessed with a Bioanalyser (Agilent Technologies, Palo Alto, CA). RNA from each sample was used to synthesize probes, using 2 successive rounds of cRNA amplification with appropriate quality control to ensure full-length synthesis according to standard Affymetrix protocols, and hybridized to mouse 430 2.0 chips (Affymetrix, Santa Clara, CA). Raw data were transformed with the Robust Multiarray Analysis (RMA) algorithm to normalize expression value. The software Cluster and Treeview were used to identify differentially expressed genes (1.5 cutoff) and classify cell subsets according to the proximity of their gene expression pattern as assessed by hierarchical clustering with complete linkage. Gene lists were analyzed using the DAVID tools accessible on the National Institute of Allergy and Infectious Diseases website. Individual significant annotations encompassing many common genes or similar biologic processes were regrouped using the “Functional annotation clustering” tool of the DAVID software.18 More information on this type of analysis is available on the DAVID website. Microarray data have been deposited at the National Center for Biotechnology Information GEO repository under accession number GSE13229 (released in March 2009).19
Statistical analyses
Statistical analyses were performed using GraphPad Prism (GraphPad Software, San Diego, CA).
Results
CD11b and CD27 surface expression define 4 subsets of NK cells
As a first step to investigate the developmental relationship between mouse NK-cell subsets defined by CD27 and CD11b surface expression, we took advantage of EGFP expression in NDE transgenic mice (Figure 1A). Indeed, these mice express EGFP and DT receptor under the control of the promoter region of the NK-cell marker NKp46/NCR1/CD335.17 We verified that nearly all EGFP+ cells were also NK1.1+CD122+ (Figure 1A). Surprisingly, we found that the CD11b and CD27 staining dissected 4 subsets, including the aforementioned CD11blow, DP, and CD27low subsets11 and an additional CD11blowCD27low or double-negative (DN) subset (Figure 1A). DN NK cells were clearly visible in previously published flow cytometry plots of CD11b and CD27 expression11 but had not been characterized in detail. We therefore investigated their distribution and surface phenotype. DN NK cells were present at very low frequency in all lymphoid and nonlymphoid organs analyzed of NDE or C57BL/6 mice (Figure 1B and data not shown). They were not NK-like T cells20 as they were also present in RAG−/− mice (data not shown). DN cells expressed NK1.1, NKG2D, NKp46, and granzyme B at the same level as the other subsets (Figure 1C and data not shown). However, they expressed low levels of Ly49 receptors, low levels of CD43 and DX5 (CD49b), but reciprocally high levels of interleukin-7 (IL-7) receptor (CD127), and NKG2A/C/E surface proteins compared with the other NK-cell subsets (Figure 1C). Thus, DN NK cells may represent immature NK cells present in various numbers in bone marrow and peripheral sites, such as spleen, liver, lymph nodes, thymus, and peripheral blood.
DN, CD11blow, DP, and CD27low NK cells develop sequentially in NDE mice after DT injection or during development of C57BL/6 mice
We previously described that DT injection leads to a transient depletion of EGFP+ NK cells in NDE mice.17 Complete depletion of NK cells occurs in 3 days. EGFP+ NK cells start to reappear at day 4 after DT injection, and complete reconstitution takes 3 to 4 weeks (Figure S1A, available on the Blood website; see the Supplemental Materials link at the top of the online article; and data not shown). We took advantage of this unique model to investigate the kinetics of appearance of each NK-cell subset during reconstitution. For this, we monitored the expression of CD27 and CD11b by EGFP+ NK cells from different organs at various times after DT injection. Our results show that the CD11b and CD27 profiles of EGFP+ NK cells changed drastically over time, as illustrated for spleen NK cells in Figure 2A.
When looking at the frequency of each subset among NK cells at different times after DT injection (Figure 2B), several observations can be made. (1) At early time points (days 4 and 5), the DN and CD11blow fractions were overrepresented compared with the later time points. (2) In a second phase (days 7-10), DP NK cells predominated in NK cells from all organs; and (3) in the later phase of the kinetics (days 14-28), CD27low NK cells progressively represented the majority of the cells, especially in the spleen, liver, and PBMCs. When looking at cell numbers in the spleen and bone marrow (Figure S1B), a similar phenomenon was observed, namely, that DN/CD11blow NK cells were overrepresented at day 5, whereas DP and later CD27low NK cells successively accumulated. Next, we wanted to determine whether similar kinetics also occurred during the development of the NK-cell pool in the first weeks of life. To test this, we investigated the CD11b and CD27 phenotypes of splenic NK cells of wild-type mice of different ages. As shown in Figure 2C, the frequency of DN and CD11blow subsets in NK cells was very high at early (3-4 weeks) time points. The frequency of DP NK cells was also high at early time points but decreased more slowly to reach homeostasis at 12 weeks of age. Finally, the frequency of CD27low NK cells rose progressively and started to plateau at 8 to 12 weeks of age. Altogether, these results show that, after DT-induced depletion or during normal development in wild-type mice, DN and CD11blow NK cells arise first, before DP and CD27low NK cells, which accumulate gradually during consecutive phases.
Global transcriptomic analyses support a CD11blow → DP → CD27low pathway of differentiation
To confirm the relevance of CD11b/CD27 cell-surface expression in the identification of NK-cell subsets, we compared the gene expression profile of flow cytometry–sorted NK-cell subsets using pan-genomic microarrays. This analysis was restricted to CD11blow, DP, and CD27low NK cells because DN NK-cell numbers were too limiting for this type of experiment. We performed a hierarchical clustering with complete linkage of differentially expressed genes (Figure 3A). Differentially expressed genes were represented in only 5 clusters, of which 2 (clusters I and IV, Figure 3A) predominated (39% and 43% of differential probe sets, respectively; Figure 3A). Strikingly, clusters I and IV, which included 466 and 341 probe sets, respectively (Table S1), corresponded to genes that followed either the relationship CD11blow > DP > CD27low or the reverse relationship, CD27low > DP > CD11blow. These data thus provided a strong supportive evidence of both the relevance of CD11b/CD27 cell-surface phenotype as markers of NK-cell differentiation and the CD11blow → DP → CD27low model of differentiation.
This pan-genomic analysis also permitted mining for “unbiased” functional differences between NK-cell subsets, focusing on the 2 main clusters of differentially expressed genes (clusters I and IV). For each cluster (CD11blow > DP > CD27low and CD27low > DP > CD11blow), we analyzed the gene ontology terms using DAVID data mining tools that include gene ontology analyses and PubMed searches18 (Table S2). “Cell cycle” was by far the most statistically significant term associated with the CD11blow NK-cell signature (cluster IV), as a large fraction of the corresponding genes were involved in cell cycle, such as cyclins and cyclin-dependent kinases, in chromosome condensation or in microtubule-based movement (Figure 3B). Moreover, when we focused on the transcription factors differentially expressed by CD11blow and CD27low subsets (Tables 1,2), we found in the CD11blow cluster, proto-oncogenes, such as Myc and Rel, as well as genes that have been associated with clonogenic potential in lymphocytes, such as Myb,21 or in pluripotency, such as Klf4.22 Less significantly, the analysis of the gene signature of CD11blow NK cells also revealed their expression of several transcripts associated with organogenesis or T-cell activation (Figure 3B). The “organogenesis” cluster contained genes, such as lymphotoxins α and β, which are also expressed at high levels by lymphoid tissue inducer cells or by gut CD3−NKp46+ cells.23 The “T-cell activation” cluster contains genes shared with the T-cell lineage, such as IL-7 receptors and CD27 and somatic T-cell receptor-γ/δ transcripts, which have been previously shown to be expressed by NK cells.20
Probe set . | Unigene ID . | Gene title . | Gene symbol . |
---|---|---|---|
1430435_at | Mm.336679 | AF4/FMR2 family, member 3 | Aff3 |
1416135_at | Mm.203 | Apurinic/apyrimidinic endonuclease 1 | Apex1 |
1419959_s_at | NA | RIKEN cDNA C330003B14 gene | C330003B14Rik |
1441579_at | Mm.130167 | Doublesex and mab-3 related transcription factor like family A1 | Dmrta1 |
1441107_at | Mm.32825 | Doublesex and mab-3 related transcription factor like family A2 | Dmrta2 |
1440551_at | Mm.246674 | DnaJ (Hsp40) homolog, subfamily C, member 1 | Dnajc1 |
1437187_at | Mm.11747 | E2F transcription factor 7 | E2f7 |
1418635_at | Mm.219460 | ets variant gene 3 | Etv3 |
1417541_at | Mm.57223 | Helicase, lymphoid specific | Hells |
1426114_at | NA | Heterogeneous nuclear ribonucleoprotein A/B | Hnrpab |
1439998_at | NA | Jumonji domain containing 1C | Jmjd1c |
1417394_at | Mm.4325 | Kruppel-like factor 4 (gut) | Klf4 |
1418936_at | Mm.86646 | v-maf musculoaponeurotic fibrosarcoma oncogene family, protein F (avian) | Maff |
1446570_at | Mm.116898 | Mastermind-like 2 (Drosophila) | Maml2 |
1439065_x_at | Mm.270762 | Similar to zinc finger protein 665 | MGC117846 |
1421317_x_at | Mm.52109 | Myeloblastosis oncogene | Myb |
1424942_a_at | Mm.2444 | Myelocytomatosis oncogene | Myc |
1420710_at | Mm.4869 | Reticuloendotheliosis oncogene | Rel |
1436325_at | NA | RAR-related orphan receptor alpha | Rora |
1427677_a_at | Mm.434350 | SRY-box containing gene 6 | Sox6 |
1429951_at | Mm.343095 | Single-stranded DNA binding protein 2 | Ssbp2 |
1433471_at | Mm.31630 | Transcription factor 7, T-cell specific | Tcf7 |
1453076_at | Mm.6922 | RIKEN cDNA 9130211I03 gene | 9130211I03Rik |
Probe set . | Unigene ID . | Gene title . | Gene symbol . |
---|---|---|---|
1430435_at | Mm.336679 | AF4/FMR2 family, member 3 | Aff3 |
1416135_at | Mm.203 | Apurinic/apyrimidinic endonuclease 1 | Apex1 |
1419959_s_at | NA | RIKEN cDNA C330003B14 gene | C330003B14Rik |
1441579_at | Mm.130167 | Doublesex and mab-3 related transcription factor like family A1 | Dmrta1 |
1441107_at | Mm.32825 | Doublesex and mab-3 related transcription factor like family A2 | Dmrta2 |
1440551_at | Mm.246674 | DnaJ (Hsp40) homolog, subfamily C, member 1 | Dnajc1 |
1437187_at | Mm.11747 | E2F transcription factor 7 | E2f7 |
1418635_at | Mm.219460 | ets variant gene 3 | Etv3 |
1417541_at | Mm.57223 | Helicase, lymphoid specific | Hells |
1426114_at | NA | Heterogeneous nuclear ribonucleoprotein A/B | Hnrpab |
1439998_at | NA | Jumonji domain containing 1C | Jmjd1c |
1417394_at | Mm.4325 | Kruppel-like factor 4 (gut) | Klf4 |
1418936_at | Mm.86646 | v-maf musculoaponeurotic fibrosarcoma oncogene family, protein F (avian) | Maff |
1446570_at | Mm.116898 | Mastermind-like 2 (Drosophila) | Maml2 |
1439065_x_at | Mm.270762 | Similar to zinc finger protein 665 | MGC117846 |
1421317_x_at | Mm.52109 | Myeloblastosis oncogene | Myb |
1424942_a_at | Mm.2444 | Myelocytomatosis oncogene | Myc |
1420710_at | Mm.4869 | Reticuloendotheliosis oncogene | Rel |
1436325_at | NA | RAR-related orphan receptor alpha | Rora |
1427677_a_at | Mm.434350 | SRY-box containing gene 6 | Sox6 |
1429951_at | Mm.343095 | Single-stranded DNA binding protein 2 | Ssbp2 |
1433471_at | Mm.31630 | Transcription factor 7, T-cell specific | Tcf7 |
1453076_at | Mm.6922 | RIKEN cDNA 9130211I03 gene | 9130211I03Rik |
From the gene expression analysis described in Figure 3, differentially expressed transcription factors were identified.
Probe set . | Unigene ID . | Gene title . | Gene symbol . |
---|---|---|---|
1455621_at | Mm.26636 | cDNA sequence BC066107 | BC066107 |
1438757_at | Mm.129275 | RIKEN cDNA C130069I09 gene | C130069I09Rik |
1455334_at | Mm.41364 | RIKEN cDNA D330038O06 gene | D330038O06Rik |
1422537_a_at | Mm.34871 | Inhibitor of DNA binding 2 | Id2 |
1416714_at | Mm.334861 | Interferon regulatory factor 8 | Irf8 |
1439846_at | NA | Kruppel-like factor 12 | Klf12 |
1448890_at | Mm.26938 | Kruppel-like factor 2 (lung) | Klf2 |
1452001_at | Mm.255151 | Nuclear factor, erythroid derived 2 | Nfe2 |
1416542_at | Mm.258953 | PHD finger protein 1 | Phf1 |
1447864_s_at | Mm.192094 | pogo transposable element with KRAB domain | Pogk |
1418622_at | Mm.397471 | RAB2, member RAS oncogene family | Rab2 |
1449361_at | Mm.94519 | T-box 21 | Tbx21 |
1432158_a_at | Mm.279752 | Trafficking protein particle complex 2 | Trappc2 |
1443880_at | Mm.27258 | Zinc finger and BTB domain containing 39 | Zbtb39 |
1422748_at | Mm.259595 | Zinc finger homeobox 1b | Zfhx1b |
1453212_at | Mm.333743 | Zinc finger protein 383 | Zfp383 |
1449414_at | Mm.42140 | Zinc finger protein 53 | Zfp53 |
1419239_at | Mm.12916 | Zinc finger protein 54 | Zfp54 |
1424261_at | Mm.72124 | Zinc finger protein 672 | Zfp672 |
1456324_at | Mm.422831 | Zinc finger protein 748 | Zfp748 |
1455945_at | NA | RIKEN cDNA A530094I17 gene | A530094I17Rik |
1457566_at | Mm.33443 | RIKEN cDNA A830058L05 gene | A830058L05Rik |
Probe set . | Unigene ID . | Gene title . | Gene symbol . |
---|---|---|---|
1455621_at | Mm.26636 | cDNA sequence BC066107 | BC066107 |
1438757_at | Mm.129275 | RIKEN cDNA C130069I09 gene | C130069I09Rik |
1455334_at | Mm.41364 | RIKEN cDNA D330038O06 gene | D330038O06Rik |
1422537_a_at | Mm.34871 | Inhibitor of DNA binding 2 | Id2 |
1416714_at | Mm.334861 | Interferon regulatory factor 8 | Irf8 |
1439846_at | NA | Kruppel-like factor 12 | Klf12 |
1448890_at | Mm.26938 | Kruppel-like factor 2 (lung) | Klf2 |
1452001_at | Mm.255151 | Nuclear factor, erythroid derived 2 | Nfe2 |
1416542_at | Mm.258953 | PHD finger protein 1 | Phf1 |
1447864_s_at | Mm.192094 | pogo transposable element with KRAB domain | Pogk |
1418622_at | Mm.397471 | RAB2, member RAS oncogene family | Rab2 |
1449361_at | Mm.94519 | T-box 21 | Tbx21 |
1432158_a_at | Mm.279752 | Trafficking protein particle complex 2 | Trappc2 |
1443880_at | Mm.27258 | Zinc finger and BTB domain containing 39 | Zbtb39 |
1422748_at | Mm.259595 | Zinc finger homeobox 1b | Zfhx1b |
1453212_at | Mm.333743 | Zinc finger protein 383 | Zfp383 |
1449414_at | Mm.42140 | Zinc finger protein 53 | Zfp53 |
1419239_at | Mm.12916 | Zinc finger protein 54 | Zfp54 |
1424261_at | Mm.72124 | Zinc finger protein 672 | Zfp672 |
1456324_at | Mm.422831 | Zinc finger protein 748 | Zfp748 |
1455945_at | NA | RIKEN cDNA A530094I17 gene | A530094I17Rik |
1457566_at | Mm.33443 | RIKEN cDNA A830058L05 gene | A830058L05Rik |
From the gene expression analysis described in Figure 3, differentially expressed transcription factors were identified.
Conversely to CD11blow NK cells, the most statistically significant terms associated with CD27low NK cells were relevant to effector functions (Table S2; Figure 3B). In particular, several of the most important genes in NK cell–mediated cytotoxicity, such as granzyme B, perforin, and various NK-cell receptors, were up-regulated in CD27low NK cells (Figure 3B). A large fraction of genes involved in the secretory pathway, such as Rab proteins, were also up-regulated in this subset and may correspond to genes involved in the trafficking of cytotoxic granules. These results correlated with the high expression by CD27low NK cells of the T-bet (Tbx21) transcription factor (Tables 1,2), which has been shown to be required for the cytotoxic function.24 A further look at the CD27low signature (Table S1) also revealed the presence of various molecules involved in trafficking (CX3CR1), cell-cell adhesion (ITGAM, JAM4, MCAM), cytokine response (IL-12 receptor), or in the induction of inflammation (chemokines CCL3, CCL4), further highlighting the “effector” nature of the DP and CD27low NK-cell subsets. Thus, according to gene ontology analysis, CD11blow and CD27low NK cells largely differ by their specialization in proliferation and effector functions, respectively.
DN NK cells display the highest rate of homeostatic proliferation among NK-cell subsets
Next, we compared the capacity of NK-cell subsets to proliferate at steady state or during development of the NK-cell pool. For this, NDE mice first received DT or phosphate-buffered saline. Seven days later, mice were given BrdU twice at 3-hour intervals, and BrdU incorporation was then measured in NK-cell subsets in the bone marrow. We chose this short-pulse protocol to minimize the transition of BrdU-positive NK cells from one stage of differentiation to the other. At steady state (day 0, Figure 4A), DN NK cells displayed the greater rate of proliferation, followed by CD11blow NK cells, whereas DP and CD27low NK cells displayed low proliferative capacities (Figure 4A). During NK-cell replenishment, the proliferation of all NK-cell subsets was increased (Figure 4B), but NK-cell subsets ranked in the same order as in steady-state conditions comparing their proliferation, although the statistics were less significant. As expected, most CD122+ NK-cell precursors (NKp) were also found to be BrdU positive (data not shown). Thus, the proliferative potential of NK-cell subsets decreased as follows: NKp > DN > CD11blow > DP > CD27low, further supporting that NK cells differentiate according to this developmental path.
Adoptive transfers support a DN → CD11blow → DP → CD27low pathway of differentiation
The capacity of each NK-cell subset to give rise to the other cell subsets was also assessed using adoptive transfer settings. For this, we sorted NK-cell subsets from NDE mice according to CD11b, CD27, and EGFP expression (Figure 5A) and adoptively transferred them individually into unirradiated adult C57BL/6 mice. Two weeks later, transferred cells and their progeny were identified in the spleen of recipient mice on the basis of EGFP, and we measured their levels of CD27 and CD11b expression. As shown in Figure 5B, transferred DN NK cells gave rise to all NK-cell subsets; CD11blow NK cells produced DP and CD27low NK cells, DP NK cells yielded DP and CD27low NK cells, and CD27low NK cells remained CD27low. These results show that DN NK cells are precursors of all other NK-cell subsets and strongly support the DN → CD11blow → DP → CD27low model of NK-cell differentiation.
Discussion
Surface density of CD27 and CD11b subdivides mouse NK cells into 4 subsets: DN, CD11blow, DP, and CD27low.11 Here, we first characterized the DN subset. DN NK cells belong to the NK-cell lineage as shown by their expression of activating receptors NKp46, NKG2D, and NK1.1, their high level of granzyme B, and their development in RAG−/− mice. However, DN NK cells displayed a very immature phenotype, according to a scheme of maturation proposed earlier.25 Indeed, DN cells expressed low levels of Ly49 receptors and DX5 but reciprocally high levels of IL-7R. They composed a fraction of the DX5low population of NK cells previously found to express tumor necrosis factor-related apoptosis-inducing ligand and to be abundant in fetal and neonatal mice.9
The differentiation of CD11blow NK cells into CD11bhigh NK-cells has been previously addressed by several groups using cell transfer into irradiated lymphopenic recipients.9,10,14,15 Irradiation of recipient mice clearly improves grafting of transferred cells. However, it can also induce “homeostatic” proliferation and phenotypic changes of transferred NK or T cells,16,26,–28 which complicate the interpretation of such experiments. Here, we present compelling evidence that not only demonstrate the validity of the CD11blow → DP → CD27low model of NK-cell differentiation but also complete it with the description of DN NK cells. First, we used NDE mice, an experimental model that presents 2 important advantages: (1) it allows the monitoring of NK-cell development in the absence of cell transplantation; and (2) the tracking of NK cells can be performed by EGFP expression, avoiding the use of antibodies that may bind nonspecifically to other cells, an important issue when looking at very rare subsets. Using this unique protocol, we found that DN, CD11blow, DP, and CD27low subsets arise in a sequential manner on replenishment of the NK-cell pool in NDE animals. A similar chronology was observed when looking at developing NK cells after the first weeks of postnatal life in wild-type mice. These results suggest that DP and CD27low NK cells differentiate from DN and/or from CD11blow cells. In these experiments, no clear difference in the timing of appearance of DN and CD11blow NK cells was observed, suggesting that these 2 cell types are closely related. Second, a global analysis of gene expression profiles of CD11blow, DP, and CD27low revealed that CD11blow and CD27low NK cells are the most distant subsets and that DP cells display an intermediate gene profile between those 2 subsets. These transcriptomic analyses also revealed that CD11blow and CD27low NK cells differed mainly by their proliferative potential and by their effector functions. This dissociation is a general feature of cellular differentiation and can also be taken as a supportive evidence of the precursor nature of CD11blow NK cells. These data were confirmed by measurements of BrdU incorporation by NK-cell subsets. Indeed, either at steady-state or during replenishment of the NK-cell pool, DN and CD11blow NK cells displayed a higher proliferation rate than the other subsets, which confirms and extends previous findings.10,29 Importantly, DN NK cells proliferated more than CD11blow NK cells but less than CD122+ NKp, suggesting that they represent an intermediate stage between NKp and CD11blow NK cells. Finally, we sorted NK-cell subsets on the basis of EGFP, CD27, and CD11b expression and transferred them into normal nonirradiated recipients. Two weeks later, the analysis of transferred cells in the spleen of recipient mice demonstrated that DN NK cells could give rise to all other subsets, whereas this potential was gradually lost in other subsets: CD11blow gave rise to themselves, DP and CD27low NK cells, whereas DP gave rise to CD27low and CD27low remained CD27low. Collectively, our data strongly suggest that mouse NK-cell maturation can be modeled in a 4-stage program with the following intermediates: DN → CD11blow → DP → CD27low. CD27low NK cells are probably the end stage of NK-cell differentiation before cell death, as suggested by a previous study that showed that KLRG1+ NK cells, which are mostly CD27low, are very susceptible to spontaneous apoptosis on ex vivo isolation.30 Importantly, CD27 has been recently shown to distinguish between 2 NK-cell subsets in human PBMCs. CD27high and CD27low NK cells corresponded to the previously described CD56bright and CD56dim NK cells,31,32 suggesting that CD27 is regulated similarly during mouse and human NK-cell differentiation. Early stages of NK-cell development have also been identified in human lymph nodes and tonsils that express low levels of the CD11b integrin, unlike mature NK cells.33 Thus, the cell-surface density of CD11b and CD27 appears to be relevant in both mouse and human to dissect NK-cell developmental stages.
Where do NK cells mature and what are the signals involved? As presented here, and in previous reports, NK-cell subsets display some tissue preferences. It is therefore possible that tissue-specific environments impact on their differentiation. Dendritic cells could act as regulators of NK-cell maturation as they have already been shown to control NK-cell proliferation and effector functions.34,35
Is there a functional specialization of the NK-cell subsets? For other leukocyte subsets, maturation stages have been associated with various functions. This has been shown for dendritic cells that have distinct but complementary functions in antigen presentation.36 In an attempt to address this issue, we compared gene expression profiles between CD11blow and CD27low NK cells. CD11blow NK cells were found to preferentially express genes involved in “organogenesis,” such as lymphotoxins or bone morphogenetic proteins. Moreover, we found that DN NK cells expressed RORγt transcripts at high levels (data not shown). This suggests that immature NK cells may share organogenetic properties with lymphoid tissue inducer cells, an issue that starts to be appreciated.23,37,,–40 The most striking difference between CD11blow and CD27low NK cells according to gene profiling analysis was the dissociation between proliferative capacity and effector functions. Indeed, “cell cycle” was the most significant term found associated with the signature of the CD11blow subset, whereas, conversely, genes involved in effector functions, such as cytotoxicity, were up-regulated in CD27low NK cells. This result suggests that the specific function of DN/CD11blow NK cells is probably to act as precursors capable of giving rise to effector cells when needed. Their low number in all organs supports this conclusion. However, DP and CD27low NK cells are also clearly capable of proliferating as shown by BrdU experiments, suggesting that these more differentiated NK cells can also contribute to NK-cell expansion in case of inflammation. Regarding their effector functions, DP and CD27low NK cells, which are very abundant in all organs, have been shown to exhibit comparable ex vivo responses to cytokines or cross-linking of activating receptors13 but to express different chemotactic receptors.11,–13 Further studies will be required to determine whether these 2 subsets have specific roles during immune responses.
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 USC section 1734.
Acknowledgments
E.V.'s laboratory is supported by the European Union (FP6, LSHB-CT-2004-503319-Allostem), Ligue Nationale contre le Cancer (Equipe labellisée La Ligue), Agence Nationale de la Recherche, Institut National du Cancer, Institut National de la Santé et de la Recherche Médicale, Centre National de la Recherche Scientifique, and Ministère de l'Enseignement Supérieur et de la Recherche. E.V. is a scholar of the Institut Universitaire de France. L.C. is sponsored by Institut National du Cancer.
Authorship
Contribution: L.C., T.W., and E.V. designed the research; L.C., N.F., J.C., C.R., and T.W. performed experiments; and L.C., T.W., and E.V. wrote the manuscript.
Conflict-of-interest disclosure: N.F. is an employee of Innate-Pharma. E.V. is a founder and shareholder of Innate-Pharma. The remaining authors declare no competing financial interests.
Correspondence: Thierry Walzer or Eric Vivier, CIML CNRS-Inserm-Université de la Méditerranée, Campus de Luminy, Case 906, 13288 Marseille Cedex 09, France; e-mail: [email protected] or [email protected].