Abstract
Helper versus cytotoxic T lineage decision in the thymus has been studied as a model for silencing of alternative lineage genes. Although the transcription factor RUNX3 is required for the initiation of Cd4 silencing in developing CD8 T cells, it is unknown how silencing of Cd4 and other helper T lineage genes is maintained. We show that the histone methyltransferase G9a is necessary for silencing helper T lineage genes in proliferating mouse CD8 T cells. Despite normal initial Cd4 downregulation, G9a-deficient CD8 T cells derepress Cd4 and other helper lineage genes during repeated division in lymphopenia or in response to tumor Ag. However, G9a was dispensable for continued silencing of those genes in CD8 T cells that respond to infection by Listeria monocytogenes. These results demonstrate that G9a facilitates maintenance of cellular identity of CD8 T cells during cell division, which is further reinforced by inflammatory signals.
Introduction
During a binary fate decision, genes related to the opposing lineage are heritably silenced (1, 2). This silencing is achieved through the constitutive activity of transcription factors involved in the lineage-determining process or by recruitment of epigenetic machinery in a locus-specific manner, presumably by those lineage-specific transcription factors. The differentiation of the common thymocyte precursor to the helper or the cytotoxic T cell lineage in the thymus has been studied to understand the requirements for transcription factors and epigenetic gene regulation for stable lineage decisions (1–9). CD4+CD8+ double-positive (DP) thymocytes are subjected to positive selection of rearranged TCRαβ by self-peptides presented on MHC class II or MHC class I and differentiate into cells in the helper or cytotoxic T lineages, respectively. MHC class I–selected thymocytes express the transcription factor RUNX3 that establishes the silencing of helper T lineage genes, including Cd4 and Zbtb7b (10–14).
However, it is poorly understood how helper T lineage–associated genes are heritably silenced in mature cytotoxic T cells. During thymocyte development, Cd4 is transiently repressed by RUNX1, an ortholog of RUNX3, in CD4−CD8− double-negative thymocytes via direct binding to the silencer cis-element in the locus (11). This repression is subsequently reversed upon selection of a successfully rearranged Tcrb locus (15). Cd4 is expressed uniformly in DP thymocytes that give rise to helper and cytotoxic T lineage cells after positive selection. Whereas Cd4 continues to be expressed in helper lineage T cells, CD8+ cytotoxic T cells terminate Cd4 transcription by upregulating RUNX3, which binds the identical cis-element as RUNX1 (11). Deletion of the silencer element or disrupting RUNX binding sites in the silencer results in continued Cd4 expression in CD8 T cells (16, 17). However, deletion of the silencer or Runx3 in differentiated CD8 T cells does not reactivate Cd4, indicating that the initial repression but not maintenance of Cd4 silencing requires RUNX3 (16, 18). RUNX3 may therefore recruit epigenetic modifications to the Cd4 locus and loci encoding helper lineage genes, which are maintained independently of RUNX3. Because the Cd4 locus is only reversibly repressed in double-negative thymocytes but irreversibly silenced in CD8 T cells, we hypothesized that the irreversible silencing is mediated by epigenetic modifiers that specifically interact with RUNX3 but not RUNX1.
In this study, we identified the histone methyltransferase (HMT) G9a as an epigenetic modifier that preferentially interacts with RUNX3 to RUNX1 and is necessary for continued silencing of helper lineage genes in dividing CD8 T cells under noninflammatory conditions. G9a deficiency resulted in derepression of several genes, which are otherwise expressed only in CD4 T cells, although it was compensated for by the inflammatory cytokine IL-12. These results suggest that G9a and inflammatory cues cooperatively maintain the identity of CD8 T cells during their division.
Materials and Methods
Mice
C57BL/6N and B6-CD45.1 mice were purchased from Charles River Laboratories. OT-I mice (19) were purchased from Taconic. Cd4-cre (10) and CD8-cre (E8I-cre) (20) mice were obtained from D. Littman (New York University). Ehmt2-flox mice were previously described (21). All mice were generated in or have been backcrossed more than eight times to C57BL/6 mice. Unless otherwise specified, littermate cre+ Ehmt2+/+ or cre− Ehmt2F/F were used as control. All mice were maintained in the specific pathogen-free facility at Washington University School of Medicine. All experiments were conducted following a protocol approved by the Washington University Animal Studies Committee.
Coimmunoprecipitation
1200M and AKR1 cell lines were transduced with MSCV-based retrovirus as described (22). For interactome analyses, 1200M cells in which endogenous Runx1 expression had been knocked down (22) were transduced with FLAG-hemagglutinin–tagged RUNX1 or RUNX3. RUNX1- and RUNX3-interacting proteins in nuclear extract were immunoprecipitated with anti-FLAG beads (M2; Sigma-Aldrich), eluted with 3× FLAG peptide (GenScript), and analyzed by mass spectrometry at the Taplin Mass Spectrometry Facility at Harvard University. For analytical immunoprecipitation, nuclear proteins were extracted from AKR1 cells that were transduced with RUNX1, RUNX3, RUNX1, and FLAG-tagged G9a retrovirus (FLAG-G9a), or RUNX3 and FLAG-G9a. Immune complexes containing FLAG-tagged protein were precipitated with anti-FLAG, followed by immunoblotting using anti-FLAG and anti–pan-RUNX Abs (12).
Flow cytometry
T cell transfer, tumor inoculation, and Listeria monocytogenes infection
Naive CD8 cells from OT-I TCR transgenic mice were sorted by flow cytometry as Vα2+Vβ5+CD62L+CD44lo/−CD4−CD8+ cells. Cells (2 × 105) were transferred i.v. to Tcrb−/−Tcrd−/− mice. PBMCs and splenocytes were isolated from the recipient mice 4 wk after transfer and analyzed for surface marker and gene expression. For experiments with transplanted tumors, 1 × 106 E.G7-OVA (no. CRL-2113; American Type Culture Collection) cells were injected s.c. in the flank of B6-CD45.1 mice. Five days later, 1 × 106 OT-I cells were transferred i.v. to the tumor-bearing mice followed by analysis of T cells collected from the draining lymph node 7 d later. For L. monocytogenes infection, 5 × 103 OT-I cells were transferred i.v. into B6-CD45.1 mice, which were infected i.v. with 2 × 104 CFU of L. monocytogenes expressing OVA (Lm-OVA) on the next day as described (23).
Quantitative RT-PCR
Total RNA was purified using TRIzol (Life Technologies) and reverse transcribed using qScript (QuantaBio). Gene expression was quantitated using a Luminaris SYBR Green reagent (Thermo Fisher) and a Roche LightCycler 480. Primer sequences are listed in Supplemental Table I.
RNA sequencing
RNA sequencing (RNA-seq) was done essentially as described previously (24) using total RNA extracted from 5 × 104 Ehmt2−/− and Ehmt2+/+ Vα2+CD8+ T cells purified from Tcrb−/−Tcrd−/− recipient mice 4 wk after transfer. Sequenced tags were mapped to the mouse genome mm9 using TopHat (25) with default setting, followed by transcript assembly and estimation of expression levels using Cufflinks (26–29) on Galaxy (https://usegalaxy.org/).
Chromatin immunoprecipitation
Mononucleosomes were prepared from cultured CD8 T cells by micrococcal nuclease digestion as described (30). Trimethylation of histone 3 lysine 9 (H3K9me3)–modified nucleosomes were immunoprecipitated using anti-H3K9me3 (8898; Abcam) conjugated with protein G magnetic beads (Life Technologies). For genome-wide analysis, purified DNA from precipitated nucleosomes was sequenced with a HiSeq 2500 sequencer (Illumina) with a 50-bp single end read option as described (23).
Statistical analysis
All statistics were performed using Graphpad Prism (version 7.0) using nonparametric two-tailed Student t tests for comparing two groups. Multiple groups were analyzed using two-tailed ANOVA. All other statistics were performed as described in the manuscript. Statistical analyses are shown with the mean ± SD. A p value < 0.05 was considered significant.
Results and Discussion
G9a forms a complex with RUNX3
To identify candidate epigenetic modifiers that are recruited by RUNX3, we analyzed RUNX1 and RUNX3 interactomes in 1200M thymoma cells, which have active Cd4 silencing machinery (15). Among 71 DNA-binding proteins and epigenetic modifiers identified in RUNX1 or RUNX3 interactomes, all three components of the G9a HMT complex, G9a, GLP, and WIZ, as well as a known G9a-interacting transcription factor CUX1 (31) were found predominantly in the RUNX3 interactome (Fig. 1A, 1B). Preferential interaction between G9a and RUNX3 compared with RUNX1 was confirmed by coimmunoprecipitation against FLAG-G9a and immunoblotting with anti–pan-RUNX Ab that detects both RUNX1 and RUNX3 (Fig. 1C). These data show that G9a and RUNX3 form a complex, which may deliver the HMT activity to Cd4 and other helper lineage–related genes repressed by RUNX3 in developing CD8 T cells.
G9a preferentially forms a complex with RUNX3 to RUNX1. (A) Heat map showing the 71 DNA-interacting proteins that were coimmunoprecipitated with only RUNX3 (cluster 1), RUNX1 and RUNX3 (cluster 2), or only RUNX1 (cluster 3) from transduced 1200M cells. (B) List of the transcription factors, DNA methyltransferases, and histone-modifying enzymes in (A). (C) Immunoblotting (IB) for RUNX proteins coimmunoprecipitated with FLAG-G9a from lysates of AKR1 cells transduced with RUNX1, RUNX3, or FLAG-G9a. Blots are representative of two experiments.
G9a is required for silencing of helper lineage–associated genes in proliferating CD8 T cells in vivo under noninflammatory conditions
To define the role of G9a in CD8 T cells, we conditionally inactivated Ehmt2, encoding G9a, at the DP stage of thymocyte development using Cd4-cre. Ehmt2 mRNA was barely detectable in CD8+ mature thymocytes from Ehmt2F/F Cd4-cre mice (data not shown) (referred to as Ehmt2−/− CD8 T cells hereafter). As previously reported using pLck-cre (32), numbers and frequencies of total thymocytes and CD4+ and CD8+ splenic T cells were comparable between Ehmt2−/− and control Ehmt2+/+ mice (data not shown). Different from the reported phenotype in Runx3−/− thymocytes (11, 14), CD4 was normally downregulated in mature CD8 thymocytes and was not expressed in splenic naive or memory CD8 T cells under steady-state conditions (data not shown).
To determine whether G9a is required for maintaining Cd4 silencing during cell division, we adoptively transferred Ehmt2−/− or control Ehmt2+/+ naive CD4−CD8+ T cells expressing the OT-I TCR transgene into Tcrb−/−Tcrd−/− mice, in which donor-derived CD8 T cells divide under noninflammatory conditions. By 4 wk after transfer, both Ehmt2−/− CD8 T cells expanded at similar rates if not faster than control Ehmt2+/+ cells as determined by CFSE, and they repopulated in the recipients’ peripheral blood (Fig. 2A, data not shown). Although Ehmt2+/+ CD8 T cells remained CD4−, ∼30% of Ehmt2−/− CD8 T cells that had diluted CFSE beyond the limit of detection upregulated CD4 (Fig. 2A, 2B). Because a similar result was observed with transferred Vα2+ CD8 T cells from Ehmt2F/F CD8-cre OT-I mice, in which Ehmt2 was deleted after positive selection (Fig. 2B), it is unlikely that Cd4 derepression is secondary to deregulated thymocyte selection in the absence of G9a.
G9a is required to maintain silencing of helper lineage genes in CD8 T cells during lymphopenia- or tumor Ag–driven proliferation. (A and B) CD4 expression and CFSE dilution of CD8 T cells in PBMC of Tcrb−/−Tcrd−/− mice that received Ehmt2−/− or Ehmt2+/+ OT-I T cells 4 wk prior to the analysis. Data are pooled from three experiments in which one donor of each genotype was transferred into two to three recipients. (C and D) RNA-seq analysis of CD4+CD8+ Ehmt2−/−, CD4−CD8+ Ehmt2−/−, and CD4−CD8+ Ehmt2+/+ OT-I T cells harvested from Tcrb−/−Tcrd−/− mice 4 wk after transfer. Quantification of genes with ≥1 fragments per kb per million in Ehmt2−/− or Ehmt2+/+ samples and >2-fold difference in expression is indicated for each genotype. Dashed red lines indicate 2-fold change between genotypes. (E) Heat maps showing genes differentially expressed between CD4+CD8+ or CD4−CD8+ Ehmt2−/− and control Ehmt2+/+ CD8 T cells. Values represent the log2 fold change of the mean of two to four mice compared with Ehmt2+/+ CD8 T cells. Highlighted genes represent genes differentially expressed between CD4 and CD8 memory T cells from Immunological Genome Project datasets. (F and G) Expression of CD4 of OT-I T cells in the lymph node draining transplanted E.G7-OVA tumors. n = 6–8 in two experiments.
To determine whether Ehmt2−/− CD8 T cells derepressed additional helper lineage genes, global gene expression in CD4+CD8+ and CD4−CD8+ Ehmt2−/− T cells as well as control Ehmt2+/+ CD8+ T cells 4 wk after transfer was profiled by RNA-seq. Approximately 1100 genes were differentially expressed by >2-fold with most (637 genes) being upregulated in Ehmt2−/− CD8 T cells compared with Ehmt2+/+ CD8 T cells (Fig. 2C, 2D). Among the genes that were differentially expressed between CD4 and CD8 memory T cells in the Immunological Genome Project datasets (>1.8-fold difference), 92 genes that are more highly expressed in CD4 T cells, including Cd4, Foxp3, Cd40lg, Rorc, Rora, Zbtb7b, and Il21, were derepressed in the Ehmt2−/− CD4+CD8+ T cells (Fig. 2E). We also observed downregulation of 18 genes that are more highly expressed in CD8 T cells in the absence of G9a (Fig. 2E), suggesting that G9a also contributes to turning on genes in CD8 T cells directly or indirectly. A similar change in gene expression, albeit to lesser extent, was also observed in Ehmt2−/− CD4−CD8+ T cells in which we have confirmed Ehmt2 deletion was also nearly complete (Fig. 2E, data not shown). In Ehmt2−/− CD8 T cells, the level of total dimethylation of histone 3 lysine 9 (H3K9me2) was substantially reduced whereas that of H3K9me3 was unchanged (Supplemental Fig. 1A). Additionally, we did not observe a reduction in H3K9me3 deposition near transcriptional start sites of the genes that were upregulated in CD4+CD8+ Ehmt2−/− T cells compared with Ehmt2+/+ CD8 T cells (Supplemental Fig. 1B). Cd4 upregulation was also observed in Ehmt2−/− OT-I T cells in response to E.G7-OVA tumor cells transplanted to wild-type mice (Fig. 2F, 2G). These results suggest that derepression of helper lineage genes in proliferating Ehmt2−/− CD8 T cells occurs also in lymphocyte-replete mice although it may not continue once they slow down or stop division, such as memory CD8 T cells under steady-state conditions. Collectively, these data indicate that G9a is required for continued silencing of a subset of helper lineage–associated genes in dividing CD8 T cells, which appears independent of H3K9me3. Because the G9a/GLP complex recruits PRC2 to its repressive target loci in embryonic stem cells (33), these helper lineage–associated genes may be kept silenced through the G9a-mediated recruitment of PRC2 activity.
G9a is dispensable for silencing helper lineage genes in the presence of strong TCR or IL-12R signals
To determine whether G9a is required for continued Cd4 silencing in dividing CD8 T cells in response to infection, Ehmt2−/− or Ehmt2+/+ OT-I T cells (Thy1.2/CD45.2) mixed with internal control congenic OT-I T cells (Thy1.1/CD45.2) were transferred into CD45.1 wild-type mice, which were subsequently infected with Lm-OVA. In contrast to lymphopenia- and tumor-driven proliferation, Ehmt2−/− OT-I T cells remained CD4− (Fig. 3A). Additionally, we observed comparable expansion and the ability to produce IFN-γ of Ehmt2−/− and Ehmt2+/+ OT-I T cells relative to control OT-I T cells (Fig. 3B, data not shown). These results indicate that G9a is dispensable for CD8 T cells under inflammatory conditions to maintain Cd4 silencing and express IFN-γ, which is also RUNX3-dependent (34).
Increased TCR signaling compensates for G9a deficiency in silencing helper lineage genes in CD8 T cells. (A and B) CD4 and CD8 expression of Ehmt2−/− and Ehmt2+/+ OT-I T cells (Thy1.2/CD45.2) that were cotransferred as a 1:1 mixture with Thy1.1/CD45.2 OT-I T cells into CD45.1 mice 4 d after Lm-OVA infection. (B) Ratios of Thy1.2+ to Thy1.1+ donor cells 4 d postinfection. Data points represent individual recipient in three experiments, in which cells from one donor per genotype were transferred into three recipients. (C and D) Expression of CD4 and CD8 on Ehmt2+/+ (C) or Ehmt2−/− (D) T cells cultured in the presence of indicated concentrations of anti-CD3 and anti-CD28 Abs. Data are representative of three experiments (n = 2 per genotype per experiment). (E) Percentage of CD4+ cells in cultured Ehmt2−/− or Ehmt2+/+ CD8 T cells shown by mean ± SD. (F) Quantitative PCR analysis of gene expression in Ehmt2−/− and Ehmt2+/+ CD8 T cells cultured in either the CD28hi or CD28lo condition. Ehmt2+/+ Th1 CD4 T cells were used as control (Th1).
Distinct dependency of silencing of helper lineage genes on G9a expressed in CD8 T cells between inflammatory and noninflammatory conditions suggests that cell-extrinsic signals through TCR, costimulatory molecules, or cytokine receptors engage compensatory pathways that reinforce gene silencing. To define such cell-extrinsic determinants, we first cultured naive polyclonal CD8 T cells with varying concentrations of anti-CD3 and anti-CD28 Abs and determined whether distinct intensities of signaling through TCR or CD28 alter Cd4 derepression in Ehmt2−/− CD8 T cells. As seen in the lymphopenic condition, a substantial fraction of Ehmt2−/− CD8 T cells derepressed CD4 when they were cultured with low anti-CD3 and high anti-CD28 Ab concentrations (CD28hi condition), whereas CD4 was barely expressed in Ehmt2+/+ CD8 T cells (Fig. 3C–E). Ehmt2−/− CD8 T cells also upregulated additional helper lineage–related genes, such as Il21 and Rorc (Fig. 3F). In contrast, when Ehmt2−/− CD8 T cells were stimulated with high anti-CD3 and low anti-CD28 Ab concentrations (CD28lo condition), derepression of the helper lineage genes was markedly reduced (Fig. 3D–F).
Furthermore, Cd4 derepression in Ehmt2−/− CD8 T cells cultured in the CD28hi condition was significantly inhibited by the cytokine IL-12 (Fig. 4A, 4B). We detected elevated H3K9me2 in Ehmt2−/− CD8 T cells cultured in the presence of IL-12 compared with those without IL-12 (Fig. 4C, 4D). The IL-12 treatment upregulated the H3K9me3 demethylase Kdm4c as well as GLP/Ehmt1 by 3-fold (Fig. 4E), thus possibly maintaining H3K9me2-dependent gene regulation by increasing demethylation of H3K9me3 by KDM4C or by elevating residual HMT activity of GLP. These results suggest that the inflammation-dependent compensation may reinforce stable lineage-specific gene expression signature in CD8 T cells that proliferate in response to infection.
Signals through IL-12 receptor compensate for G9a deficiency in repression of helper lineage genes in CD8 T cells. (A and B) Percentages of CD4+ cells in Ehmt2−/− or control Ehmt2+/+ CD8 T cells cultured in the CD28lo condition with or without murine IL-12 (mIL-12; 10 ng/ml). Plots are representative of two independent experiments (two to three mice per experiment). (C) Immunoblotting for total H3 and H3K9me2 of lysates from Ehmt2−/− or Ehmt2+/+ CD8 T cells cultured in the CD28hi or CD28lo condition. Data are representative of two experiments. (D) Relative expression of Ehmt2 in CD8 T cells from Ehmt2−/− and Ehmt2+/+ mice cultured with the CD28hi condition with or without mIL-12. (E) Quantitative PCR analysis of expression of H3K9 HMTs and demethylases in Ehmt2−/− CD8 T cells cultured in the CD28hi condition with or without mIL-12. n = 3.
Our study has demonstrated that G9a is required for maintaining silencing of multiple helper lineage–associated genes, such as Cd4, in dividing CD8 T cells in response to lymphopenia or tumor Ag. In developing CD8 T cells, Cd4 is shut off by RUNX3 in a G9a-independent manner. However, the continued silenced state, which is independent of RUNX3, is not maintained in the absence of G9a, suggesting that transient RUNX3-dependent recruitment of G9a establishes the heritably silenced states of the locus in cooperation with other factors, such as additional methyltransferases. Alternatively, G9a is constitutively recruited to the Cd4 locus initially by RUNX3 and subsequently by a RUNX3-independent mechanism. All the G9a-dependent repression targets are not RUNX3 targets (18, 34), and Ehmt2−/− CD8 T cells are also able to proliferate and express IFN-γ, which is dependent on RUNX3. Therefore, there are multiple distinct RUNX3- or G9a-containing complexes that regulate gene activation or repression in CD8 T cells, whereas Cd4 silencing and repression of some of TFH-signature genes (18), including Icos, Cxcr5, and Il21, appear to be dependent on a complex containing both.
The absence of Cd4 derepression in CD8 T cells responding to Lm-OVA infection may be explained by high levels of Ag and IL-12, both of which compensate for G9a deficiency in CD8 T cells in vitro. Although the exact mechanism is unknown, our data suggest that IL-12R signaling alters the balance between histone methylation and demethylation to increase the levels of H3K9me2 independent of G9a, possibly facilitating heritable gene silencing in dividing CD8 T cells through the compensatory pathways. CD8 T cells may thus engage multiple epigenetic pathways in a context-dependent manner to shape their gene expression signature.
Disclosures
The authors have no financial conflicts of interest.
Acknowledgments
We thank D. Littman for cell lines, Cd4-cre mice, and CD8-cre mice, C.-S. Hsieh and E.M. Oltz for discussions, and J. Chaudhuri, P. Collins, M. Holmgren, S. Hsiung, S. Raju, and E. Russler-Germain for technical assistance.
Footnotes
This work was supported by National Institutes of Health Grants R01AI097244, R56AI114593 (to T.E.), and T32GM007200 (to D.J.V.), by the Leukemia and Lymphoma Society Scholar award (to T.E.), and by the Shawn Hu and Angela Zeng Fellowship (to D.J.V.).
The mass spectrometry data presented in this article have been submitted to PRIDE (https://www.ebi.ac.uk/pride/) under accession number PXD007511. The RNA sequencing data and chromatin immunoprecipitation sequencing data presented in this article have been submitted to the National Center for Biotechnology Information’s Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE101730) under accession number GSE101730.
The online version of this article contains supplemental material.
Abbreviations used in this article:
- DP
- double-positive
- FLAG-G9a
- FLAG-tagged G9a retrovirus
- H3K9me2
- dimethylation of histone 3 lysine 9
- H3K9me3
- trimethylation of histone 3 lysine 9
- HMT
- histone methyltransferase
- Lm-OVA
- L. monocytogenes expressing OVA
- RNA-seq
- RNA sequencing.
- Received December 11, 2017.
- Accepted April 11, 2018.
- Copyright © 2018 by The American Association of Immunologists, Inc.