# The transcription factors Runx3 and ThPOK cross-regulate acquisition of cytotoxic function by human Th1 lymphocytes

Yasmina Serroukh, Chunyan Gu-Trantien, Baharak Hooshiar Kashani, Matthieu Defrance, Thien-Phong Vu Manh, Abdulkader Azouz, Aurelie Detavernier, Alice Hoyois, Jishnu Das, Martin Bizet, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Yasmina Serroukh, Chunyan Gu-Trantien, Baharak Hooshiar Kashani, Matthieu Defrance, Thien-Phong Vu Manh, et al. The transcription factors Runx3 and ThPOK cross-regulate acquisition of cytotoxic function by human Th1 lymphocytes. eLife, 2018, 7, pp.e30496. 10.7554/eLife.30496. hal-02118048

HAL Id: hal-02118048

https://hal.science/hal-02118048v1

Submitted on 18 Oct 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# The transcription factors Runx3 and ThPOK cross-regulate acquisition of cytotoxic function by human Th1 lymphocytes

Yasmina Serroukh<sup>1†</sup>, Chunyan Gu-Trantien<sup>1†</sup>, Baharak Hooshiar Kashani<sup>1†</sup>, Matthieu Defrance<sup>2</sup>, Thien-Phong Vu Manh<sup>3</sup>, Abdulkader Azouz<sup>1</sup>, Aurélie Detavernier<sup>1</sup>, Alice Hoyois<sup>1</sup>, Jishnu Das<sup>4,5</sup>, Martin Bizet<sup>2</sup>, Emeline Pollet<sup>3</sup>, Tressy Tabbuso<sup>1</sup>, Emilie Calonne<sup>2</sup>, Klaas van Gisbergen<sup>6</sup>, Marc Dalod<sup>3</sup>, François Fuks<sup>2</sup>, Stanislas Goriely<sup>1</sup>, Arnaud Marchant<sup>1*</sup>

<sup>1</sup>Institute for Medical Immunology, Université Libre de Bruxelles, Charleroi, Belgium; <sup>2</sup>Laboratoire d'Epigénétique du Cancer, Université Libre de Bruxelles, Bruxelles, Belgium; <sup>3</sup>Centre d'Immunologie de Marseille-Luminy 13288, Aix Marseille Université UM2, Marseille, France; <sup>4</sup>Ragon Institute of MGH, MIT and Harvard University, Cambridge, United States; <sup>5</sup>Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, United States; <sup>6</sup>Department of Haematopoiesis, Sanquin Research and Landsteiner Laboratory, Amsterdam, Netherlands

**Abstract** Cytotoxic CD4 (CD4<sub>CTX</sub>) T cells are emerging as an important component of antiviral and antitumor immunity, but the molecular basis of their development remains poorly understood. In the context of human cytomegalovirus infection, a significant proportion of CD4 T cells displays cytotoxic functions. We observed that the transcriptional program of these cells was enriched in CD8 T cell lineage genes despite the absence of ThPOK downregulation. We further show that establishment of CD4<sub>CTX</sub>-specific transcriptional and epigenetic programs occurred in a stepwise fashion along the Th1-differentiation pathway. In vitro, prolonged activation of naive CD4 T cells in presence of Th1 polarizing cytokines led to the acquisition of perforin-dependent cytotoxic activity. This process was dependent on the Th1 transcription factor Runx3 and was limited by the sustained expression of ThPOK. This work elucidates the molecular program of human CD4<sub>CTX</sub> T cells and identifies potential targets for immunotherapy against viral infections and cancer.

DOI: https://doi.org/10.7554/eLife.30496.001

## Introduction

The thymic differentiation of helper CD4 and cytotoxic CD8 T lymphocytes results from the opposite activity of key transcription factors (TF) including ThPOK and Runx3 repressing the expression of CD8 and CD4 T cell lineage genes, respectively. After emigration from the thymus, naive CD8 and CD4 T cells maintain the expression of Runx3 or ThPOK, respectively, suggesting that the lineage-defining role of these TF is also active in the periphery (Vacchio and Bosselut, 2016). However, the repression of a cytotoxic program in peripheral CD4 T cells is not absolute as these cells can acquire perforin-dependent cytotoxic activity (Appay et al., 2002a; van de Berg et al., 2008; Cheroutre and Husain, 2013).

*For correspondence: arnaud.marchant@ulb.ac.be

†These authors contributed equally to this work

**Competing interests:** The authors declare that no competing interests exist.

**Funding:** See page 21

**Received:** 18 July 2017

**Accepted:** 20 February 2018

**Published:** 28 February 2018

**Reviewing editor:** Ronald N Germain, National Institute of Allergy and Infectious Diseases, United States

© Copyright Serroukh et al. This article is distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use and redistribution provided that the original author and source are credited.
Initially considered as a phenomenon of peripheral importance, the acquisition of cytotoxic function by CD4 T cells is now recognized as a key component of immunity against viruses and tumors (Swain et al., 2012) and correlates with positive outcome in multiple human and animal models (Brown et al., 2012; Wilkinson et al., 2012; Johnson et al., 2015; Ma et al., 2015; Weiskopf et al., 2015; Verma et al., 2016; Xie et al., 2010; Quezada et al., 2010; Fu et al., 2013). Beside their role in viral infections and cancer, cytotoxic CD4 (CD4<sub>CTX</sub>) T cells may also have a pathogenic role in chronic inflammatory disorders (Dumitriu, 2015).

The development of CD4<sub>CTX</sub> T cells remains incompletely understood. In humans, CD4<sub>CTX</sub> T cell function is a hallmark of terminally differentiated antigen-experienced cells producing large amounts of gamma interferon (IFNγ) and low levels of interleukin-2 (IL-2) (Appay et al., 2002a; Casazza et al., 2006). This suggests that CD4<sub>CTX</sub> T cell differentiation might be induced by the prolonged stimulation of Th1 lymphocytes. On the other hand, recent studies suggest that CD4<sub>CTX</sub> T cells form a distinct lineage emerging from precursors expressing class-I restricted T cell-associated molecule (CRTAM) (Takeuchi et al., 2016). Several TF probably contribute to the differentiation of CD4<sub>CTX</sub> T cells within or outside the Th1 pathway. In the mouse intestine, the acquisition of cytotoxicity by CD4 T cells is associated with the down regulation of ThPOK and the upregulation of Runx3 (Mucida et al., 2013; Reis et al., 2013; Sujino et al., 2016). In CD8 T cells, the cytotoxic program is activated by the cooperation of TF from the T-box family (T-bet and Eomes) and Runx3, (Cruz-Guilloty et al., 2009; Pearce et al., 2003). In mouse models of cancer or neuroinflammation, Eomes was required for the induction of CD4<sub>CTX</sub> T cells (Curran et al., 2013; Raveney et al., 2015). In murine influenza infection, Blimp1 was required for CD4<sub>CTX</sub> T cell differentiation (Marshall et al., 2017). The role of these TF in human CD4<sub>CTX</sub> T cell differentiation remains to be determined.

The transcriptional program and lineage fate of effector T cells are established and maintained at the epigenetic level through DNA and histone modifications (Wilson et al., 2009; Araki et al., 2008; Sellars et al., 2015). Hypomethylation of the *PRF1* promoter, the gene encoding perforin, is associated with increased perforin expression in human CD4 T cells (Kaplan et al., 2004). The epigenetic modifications underlying the differentiation of CD4<sub>CTX</sub> T cells have not been determined.

Here, we studied circulating CD4<sub>CTX</sub> T cells isolated from the peripheral blood of cytomegalovirus-seropositive (CMV<sup>+</sup>) healthy adults. Compared to mouse models of infection (Brown et al., 2012) or cancer (Curran et al., 2013), this situation allows access to a significant number of cells presenting a fully established cytotoxic functional program at steady state. Using transcriptomic and epigenomic approaches, we defined the molecular events that dictate human CD4<sub>CTX</sub> differentiation. We further show that the increased expression of Runx3 and T-bet and key epigenetic modifications at the *PRF1* promoter, without downregulation of ThPOK, underlie the acquisition of cytotoxic function by human Th1 lymphocytes.

## Results

### Phenotype and function of in vivo differentiated perforin<sup>+</sup> human CD4 T cells

In healthy humans, chronic CMV infection is associated with the expansion of perforin<sup>+</sup>/granzyme B<sup>+</sup> CD4 (Figure 1a–b and Figure 1—source data 1) (van Leeuwen et al., 2004). In order to use this model for transcriptomic and epigenomic analyses of CD4<sub>CTX</sub> T lymphocytes, we characterized their phenotype and function in CMV<sup>+</sup> healthy adults and compared them to cytotoxic CD8 T cells. As previously reported, high perforin expression was observed in terminally differentiated CD4 and CD8 T cells that had downregulated the co-stimulatory molecules CD28 and CD27, respectively (Figure 1c) (van de Berg et al., 2008; Appay et al., 2002b). Perforin<sup>+</sup> CD4 T cells were CD8β-negative and a minority expressed low levels of CD8α (Figure 1—figure supplement 1a–b). Further analyses were conducted on sorted naive (CD45RO<sup>-</sup>CD28<sup>+</sup>) and terminally differentiated (CD28<sup>-</sup>) CD4 T cells and on naive (CD45RO<sup>-</sup>CD27<sup>+</sup>) and terminally differentiated (CD27<sup>-</sup>) CD8 T cells (Figure 1d). Increased PRF1 gene expression by CD28<sup>-</sup> CD4 and CD27<sup>-</sup> CD8 T cells was confirmed by mRNA quantification and was associated with potent cytotoxic activity in a polyclonal cell lysis assay (Figure 1e, Figure 1f and Figure 1—source data 1). This activity was abolished by Concanamycin A, supporting a perforin-dependent mechanism (Kataoka et al., 1996). Bisulphite sequencing indicated an inverse correlation between the expression of the PRF1 gene and the DNA methylation
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Phenotype and function of in vivo differentiated perforin⁺human CD4 T cells. (a–b). The expression of perforin (P) and granzyme B (G~B~) was analyzed by flow cytometry in total CD4 T cells of CMV seropositive (CMV⁺) and seronegative (CMV⁻) healthy adults. (a) Representative dot plots of log10 fluorescence. (b) Proportions (median ±interquartile range) of P⁺G~B~⁺ CD4 T cells in 11 CMV⁺ and 5 CMV⁻ subjects. (c) Perforin expression in CD28⁻ CD4 T cells and CD27⁻ CD8 T cells of CMV⁺ subjects. Numbers indicate cell proportions in each quadrant. (d) Sorting strategy of naive and Figure 1 continued on next page</figcaption></figure>
Figure 1 continued

cytotoxic T cells according to expression of CD45RO, CD28 and CD27. Representative dot plot of log10 fluorescence. (e) The expression of PRF1 mRNA was measured by qPCR in purified T cell subsets of 5 CMV⁺ subjects. Results are median ± interquartile range of log2 fold change as compared to naive CD4 T cells. **:p<0,01 and ***:p<0,01. (f) The cytolytic activity of purified T cell subsets against anti-CD3-loaded target cells was assessed with or without pre-incubation with Concanamycin A (ConA). Data are mean ± SEM of three independent experiments on cells from different donors. (f) The methylation status of the PRF1 promoter was assessed in T cell subsets by bisulphite pyrosequencing. Data are median ± interquartile range of five donors for CD28⁻CD4 T cells and HEL-299 and of 9 donors for the other indicated subsets. Grey lines indicate three regions with distinct methylation profiles. See also **Figure 1—figure supplement 1** and Source data file.

DOI: https://doi.org/10.7554/eLife.30496.002

The following source data and figure supplement are available for figure 1:

**Source data 1.** Phenotype and function of in vivo differentiated perforin +human CD4 T cells.

DOI: https://doi.org/10.7554/eLife.30496.004

**Figure supplement 1.** Expression of CD8 subunits by CD4<sub>CTX</sub> T cells.

DOI: https://doi.org/10.7554/eLife.30496.003

status of its promoter region (Figure 1g and Figure 1—source data 1). Whereas PRF1 promoter was hypermethylated in a perforin⁻ fibroblastic cell line (HEL-299), all CpG sites were hypomethylated in CD28⁻ CD4 and CD27⁻ CD8 T cells. Low DNA methylation levels were detected at intermediate (16 to 28; middle grey line on Figure 1g) CpG sites in naive CD8 T cells and at proximal (CpG sites 29 to 34; right grey line) sites in both naive CD4 and CD8 T cells, suggesting that the PRF1 gene is transcriptionally poised in naive T lymphocytes. Together, these results indicate that CD28⁻ CD4 T cells exert a cytotoxic activity comparable to CD27⁻ CD8 T cells and that this subset can therefore be used as a relevant model of in vivo differentiated CD4<sub>CTX</sub> T cells.

## The transcriptional program of CD4<sub>CTX</sub> T cells is enriched in CD8 T cell lineage genes without down regulation of ThPOK

In order to elucidate the molecular basis of CD4<sub>CTX</sub> T cell differentiation, their transcriptome was first compared to that of naive CD4, naive CD8 and CD8<sub>CTX</sub> T cells. Unsupervised analysis of transcriptional programs indicated that naive and cytotoxic T cells formed separate clusters and that CD4<sub>CTX</sub> and CD8<sub>CTX</sub> T cells were more closely related than their naive counterparts (Figure 2a–b). Gene set enrichment analysis (GSEA) was used to quantify the degree of sharing of the transcriptional program of CD4<sub>CTX</sub> T cells with the CD4 and CD8 T cell lineages. Genes expressed at higher levels in CD4<sub>CTX</sub> and CD8<sub>CTX</sub> T cells as compared to their naive counterparts were identified and their enrichment in naive CD8 and CD4 T cell transcriptomes was assessed (Supplementary file 1 and Figure 2c). As expected, genes that were upregulated in CD8<sub>CTX</sub> T cells were significantly enriched in genes of the CD8 T cell lineage (Figure 2c). Strikingly, genes upregulated in CD4<sub>CTX</sub> T cells were also enriched in CD8 rather than CD4 T cell lineage genes. The transcriptional program common to CD4<sub>CTX</sub> and CD8<sub>CTX</sub> T cells included RUNX3, TBX21 (T-bet) and EOMES, TF known to promote effector and memory functions in CD8 T cells (Supplementary file 1, Figure 2d and Figure 2—source data 1) (Cruz-Guillot et al., 2009).

Strikingly, the enrichment in CD8 T cell lineage genes by CD4<sub>CTX</sub> T cells was not associated with the downregulation of ZBTB7B (ThPOK) (Figure 2d and Figure 2—source data 1). As expected, ThPOK expression was higher in naive CD4 as compared to naive CD8 T cells. However, in contrast to mouse intestinal CD4 T cells (Mucida et al., 2013), ThPOK gene expression was not downregulated by human circulating CD4 T cells, CD4<sub>CTX</sub> T cells expressing higher levels of ThPOK mRNA than naive CD4 T cells (Figure 2d and Figure 2—source data 1). Protein expression analysis by flow cytometry confirmed mRNA expression profiles for Runx3, T-bet and Eomes (Figure 2e and Figure 2—source data 1). This analysis showed high and similar expression of ThPOK in CD4<sub>CTX</sub> and naive CD4 T cells and a higher expression of ThPOK in CD8<sub>CTX</sub> as compared to naive CD8 T cells (Figure 2e and Figure 2—source data 1).

We explored the epigenetic basis of the transcriptional program of CD4<sub>CTX</sub> T cells by analysing their DNA methylome and comparing it to that of naive CD4 and CD8<sub>CTX</sub> T cells. The acquisition of cytotoxic function by CD4<sub>CTX</sub> and CD8<sub>CTX</sub> T cells was associated with changes in methylation, primarily hypomethylation, of large numbers of genes (Figure 2—figure supplement 1a–d). Unsupervised analysis of DNA methylomes indicated that naive and cytotoxic T cells formed separate
<figure><img src="image_5.png" /><figcaption>**Figure 2.** The transcriptional program of CD4<sub>CTX</sub> T cells is enriched in CD8 lineage genes without downregulation of ThPOK. The transcriptome of T cell subsets from five CMV-seropositive donors (four for CD8<sub>CTX</sub> T cells) was analyzed by gene expression arrays. Log2 expression values of 14,372 probes with a variance >0.01 corresponding to 11,200 unique genes were submitted to unsupervised clustering (a) and principal component (b) analyses. (c) GSEA was used to test the enrichment of CD8<sub>CTX</sub> and CD4<sub>CTX</sub> T cell GeneSets (Supplementary file 1) in naive CD8 and naive CD4 T cell expression datasets. Genes showing no differential expression in CD8<sub>CTX</sub> and CD4<sub>CTX</sub> T cells were used as negative controls (n = 379). Bar codes show the ranking of the log2 fold change of gene expression values in naive CD8 versus naive CD4 T cells. Green lines represent enrichment profiles. False discovery rates (FDR) below 0.25 were considered significant and are indicated in red. NES: normalized enrichment score. (d) TF mRNA expression by T Figure 2 continued on next page</figcaption></figure>
Figure 2 continued

cell subsets purified from four to six donors, as indicated, was assessed by qPCR (upper panels). Results are expressed as median ±interquartile range of the log2 fold change as compared to naive CD4 or naive CD8 T cells. #:p<0.05 and ##:p<0.01 as compared to CD4 T cell counterparts; *:p<0.05 and **:p<0.01 as compared to naive counterparts. (e) TF protein expression was analyzed in T cell subsets by flow cytometry Upper panels: co-expression with perforin from one representative subject. Gated populations include naive CD4 (red), CD4<sub>CTX</sub> (green), naive CD8 (purple) and CD8<sub>CTX</sub> (blue) T cells. Lower panels: individual median intensity of fluorescence (MFI) of 5 CMV<sup>+</sup> subjects. *:p<0.05. NS: not significant. Naive and cytotoxic T cell subsets were gated using the markers and strategy illustrated in Figure 1d. (f) GSEA was used to determine the correlation between gene expression and DNA methylation for each indicated T cell subset. Graphs show the enrichment of indicated GeneSets (Supplementary file 1) in genes that were either hypo or hypermethylated at the level of the promoter in the corresponding T cell subset. Black rectangles represent saturated bar codes of ranked delta beta values of 205,783 probes between the indicated subsets. (g) Heatmap of methylation beta values of RUNX3, TBX21 and ZBTB7B gene promoter in T cell subsets of five donors.. See also Figure 2—figure supplement 1 and Source data file.

DOI: https://doi.org/10.7554/eLife.30496.005

The following source data and figure supplement are available for figure 2:

**Source data 1.** The transcriptional program of CD4CTX T cells is enriched in CD8 lineage genes without down regulation of ThPOK.

DOI: https://doi.org/10.7554/eLife.30496.007

**Figure supplement 1.** Whole genome methylation profiles of naive and cytotoxic CD4 and CD8 T cells.

DOI: https://doi.org/10.7554/eLife.30496.006

clusters (Figure 2—figure supplement 1e). GSEA indicated that genes expressed at higher levels in CD4<sub>CTX</sub> or CD8<sub>CTX</sub> T cells as compared to their naive counterparts were significantly enriched in hypomethylated CpG located in their promoter regions (Figure 2f). In contrast, genes that were downregulated in CD4<sub>CTX</sub> or CD8<sub>CTX</sub> T cells were not enriched in hypermethylated CpG, suggesting distinct epigenetic mechanisms in the up- and downregulation of genes upon differentiation of cytotoxic T cells. Analysis of probes located in the promoter region of TF up regulated in cytotoxic as compared to naive T cells suggested significant hypomethylation for RUNX3 and TBX21 in both CD4<sub>CTX</sub> and CD8<sub>CTX</sub> (Figure 2g and Figure 2—source data 1), but not for EOMES (data not shown). As observed during thymopoiesis (Vacchio and Bosselut, 2016), ZBTB7B promoter was significantly hypomethylated in naive CD4 as compared to naive CD8 T cells, but no significant difference was observed between CD4<sub>CTX</sub> and CD8<sub>CTX</sub> cells (Figure 2g and Figure 2—source data 1). Most analyzed RUNX3 probes (77%) were located in the distal promoter (TSS1500), a region suspected to be involved in the control of RUNX3 mRNA translation (Kim et al., 2015), whereas most TBX21 (77%) and ZBTB7B probes (76%) were located in the proximal promoter (5'UTR, 1st exon and TSS200). Together, these results indicate that the transcriptional program of human CD4<sub>CTX</sub> T cells is enriched in CD8 T cell lineage genes. Acquisition of this program involves extensive hypomethylation of the promoter regions of a large number of genes, including TF, and is not accompanied by ThPOK downregulation.

## Stepwise differentiation of CD4<sub>CTX</sub> T cells within the Th1 lymphocyte lineage

In order to decipher the molecular pathways involved in the differentiation of CD4<sub>CTX</sub> T cells, we measured the expression of cytotoxicity-related genes in subsets of memory CD4 T cells. The production of perforin was associated with the expression of the Th1 chemokine receptors CCR5 and, to a lower extend, CXCR3 but not with Th17 or Th2 receptors CCR6, CCR4 or CRTh2 (Figure 3a). (Sallusto et al., 1998; Cosmi et al., 2000; Rivino et al., 2004; Cohen et al., 2011; Couturier et al., 2014). To determine at which stage of their differentiation Th1 cells initiate the production of perforin, central memory (CM) and CD28<sup>+</sup> effector memory (EM) T cells expressing Th1 chemokine receptors were compared to naive and CD4<sub>CTX</sub> T cells (Figure 3b). No perforin<sup>+</sup> cells were detected by flow cytometry among CM CD4 T cells. A small proportion of perforin<sup>+</sup> cells were detected among CD28<sup>+</sup> EM CD4 T cells in some donors and these cells were CCR5<sup>+</sup>. Mean fluorescence intensity analysis indicated high perforin expression by CD4<sub>CTX</sub> cells and low and comparable expression in naive, CCR5<sup>+</sup> CM (CM<sub>Th1</sub>) and CCR5<sup>+</sup> CD28<sup>+</sup> EM (EM28<sup>+</sup><sub>Th1</sub>) CD4 T cells (Figure 3c and Figure 3—source data 1). In contrast, PRF1 mRNA expression analysis of sorted cells (sorting strategy shown in Figure 3—figure supplement 1a) indicated that the PRF1 gene was already expressed at higher levels in CM<sub>Th1</sub> as compared to naive cells and was further upregulated in EM28<sup>+</sup><sub>Th1</sub> cells and in CD4<sub>CTX</sub> cells (Figure 3d and Figure 3—source data 1). Notably, this pattern
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Differentiation of CD4<sub>CTX</sub> T cells within the Th1 lineage. (a–b). Expression of perforin and chemokine receptors was assessed by flow cytometry in total CD4 T cells (a) central memory (CM), effector memory (EM) CD28<sup>+</sup> and CD4<sub>CTX</sub> T cells (b). Log10 fluorescence of one representative donor out of 9. (c) Median intensity of fluorescence (MFI) of perforin expression in T cell subsets from 7 CMV<sup>+</sup> subjects. (d) mRNA expression of indicated genes was assessed by qPCR in T cell subsets purified from two to five donors, as indicated. CM<sub>Th1</sub> and EM28<sup>+</sup><sub>Th1</sub> were CCR5<sup>+</sup>CCR6<sup>-</sup> as Figure 3 continued on next page</figcaption></figure>
Figure 3 continued

illustrated in Supplementary Figure 3. Results are median ±interquartile range of the log2 fold change as compared to naive CD4 T cells. \*:p<0.05 and \*:p<0.01. (e) Methylation status of *PRF1* promoter was assessed by bisulphite pyrosequencing in indicated purified T cell subsets. Results are presented as median ±interquartile range of two to five donors depending on the CpG site, as detailed in *Supplementary file 2a*. Grey lines indicate three regions with distinct methylation profiles. (e) Histone modifications at indicated regions of the *PRF1* locus were studied by ChIP-qPCR in purified CD4 T cell subsets of three donors. Results are % of input after normalization for the enrichment in pan H3. Letters refer to indicated distances from the transcription start site (TSS) (left panel). See also *Figure 3—figure supplements 1* and *2* and Source data file.

DOI: https://doi.org/10.7554/eLife.30496.008

The following source data and figure supplements are available for figure 3:

**Source data 1.** Differentiation of CD4CTX T cells within the Th1 lineage.

DOI: https://doi.org/10.7554/eLife.30496.011

**Figure supplement 1.** Sorting strategy of memory CD4 T cell subsets and expression of candidate genes.

DOI: https://doi.org/10.7554/eLife.30496.009

**Figure supplement 2.** Transcriptional program underlying the expression of perforin in Th1 cells.

DOI: https://doi.org/10.7554/eLife.30496.010

of progressive up regulation of gene expression was similar for *RUNX3* and *TBX21* and contrasted with the high expression of *IFNG* and *EOMES* mRNA already detected at the *CM<sub>Th1</sub>* stage of differentiation (*Figure 3d*, and *Figure 3—source data 1* and *Supplementary file 1*). Analysis of sorted CM Th1, Th2 and Th17 cells indicated that the increased expression of perforin and associated TF was specific to the Th1 lineage (*Figure 3—figure supplement 1b-c*). Epigenetic analyses of the *PRF1* promoter further supported the stepwise acquisition of perforin expression within the Th1 lineage. DNA methylation levels at intermediate CpGs (sites 16 to 28) progressively decreased from naive to *CM<sub>Th1</sub>*<sup>+</sup>, *EM28<sup>+</sup><sub>Th1</sub>* and *CD4<sub>CTX</sub>* T cells (*Figure 3e* and *Supplementary file 2a*). As observed in naive and *CD4<sub>CTX</sub>* T cells, proximal CpGs (sites 29 to 34) were hypomethylated in *CM<sub>Th1</sub>* and *CD28<sup>+</sup>*<sup>EM<sub>Th1</sub></sup> T cells. H3K4me3 (active promoter) and H3K27ac (active promoter and enhancer [*Shlyueva et al., 2014*]) enrichment was analyzed in previously identified putative regulatory regions (*Figure 3f*, left panel) (*Pipkin et al., 2010*; *Adams et al., 2012*). In agreement with *PRF1* gene expression and promoter methylation, high enrichment of H3K4me3 and of H3K27ac was detected at the *PRF1* proximal promoter region in *CD4<sub>CTX</sub>* T cells (*Figure 3f*, middle and right panels and *Figure 3—source data 1*). Interestingly, an H3K27ac enrichment was observed in a region located 5,5 to 7,7 kb upstream of the TSS in Th1 cell subsets, suggesting the presence of an active enhancer.

Together, these results indicate the progressive acquisition of *PRF1* gene expression from *CM<sub>Th1</sub>* to *EM<sub>Th1</sub>* lymphocytes. Single-cell PCR analysis revealed that this process is related to a progressive increase in the proportion of *PRF1* mRNA<sup>+</sup> Th1 cells (*Figure 4a* and *Figure 4—source data 1*). Gene co-expression analysis at the single-cell level indicated that *PRF1* mRNA was co-expressed with distinct sets of TF in Th1 cell subsets (*Figure 4b* and *Figure 4—source data 1*). In *CM<sub>Th1</sub>* cells, *PRF1* expression was co-expressed with a relatively restricted set of TF, including *PRDM1*, *RUNX3* and *EOMES*. In *CD4<sub>CTX</sub>* T cells, a larger set of co-expressed TF was identified, including *TBX21*, *HOPX*, *ZNF683* (Hobit), *PRDM1* and *RUNX3* and *EOMES*. Notably, a lower proportion of perforin<sup>+</sup> *CD4<sub>CTX</sub>* T cells co-expressed *EOMES* as compared to the other co-expressed TF. In conclusion, the analysis of *PRF1* gene expression in vivo suggests a model in which Th1 lymphocytes acquire permis- sive modifications of the local chromatin environment and a network of TF factors that could underlie the acquisition of cytotoxic functions.

## Transcriptional program underlying the expression of perforin in Th1 cells

In order to identify the key steps that drive the acquisition of cytotoxic functions along the Th1 pathway, we compared the gene expression profile of *CD4<sub>CTX</sub>* T cells to that of *CM<sub>Th1</sub>* cells. Unsupervised gene expression analysis indicated that *CM<sub>Th1</sub>* cells formed a cluster separated from naive and *CD4<sub>CTX</sub>* T cells and were more closely related to *CD4<sub>CTX</sub>* T cells than to naive cells (*Figure 5a*). Including CD8 T cell subsets in this analysis indicated that the transcriptome of *CM<sub>Th1</sub>* cells was more distant to *CD8<sub>CTX</sub>* T cells than *CD4<sub>CTX</sub>* T cells (*Figure 3—figure supplement 2a-b*). 693 genes were significantly upregulated in *CM<sub>Th1</sub>* as compared to naive cells and 322 genes were up regu- lated in *CD4<sub>CTX</sub>* T cells as compared to *CM<sub>Th1</sub>* cells (*Figure 5b* and *Supplementary file 1*). Similar
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Proportion of *PRF1* mRNA⁺ cells and co-expression with transcription factors in the Th1 lineage. The expression of *PRF1* and transcription factors (TF) mRNA was analyzed in 43 single naive, central memory (CM) Th1, effector memory (EM) CD28⁺ Th1 and CD4<sub>CTX</sub> T cells from one donor. (a) *PRF1* mRNA expression in single cells and proportions of *PRF1*⁺ cells in Th1 cell subsets. (b) Heat map of *PRF1*⁺ cells co-expressing individual TF in Th1 cell subsets. (calculated proportion: double-positive cells/perforin-positive cells). See also **Figure 4—source data 1**.</figcaption></figure>

DOI: https://doi.org/10.7554/eLife.30496.012

The following source data is available for figure 4:

**Source data 1.** Proportion of PRF1 mRNA⁺ cells and co-expression with transcription factors in the Th1 lineage.

DOI: https://doi.org/10.7554/eLife.30496.013

gene numbers were downregulated in the two subsets. As expected, genes upregulated in CD4<sub>CTX</sub> T cells as compared to CM<sub>Th1</sub> cells included cytotoxicity-related molecules, among which GNLY (granulysin), granzymes, CD107a (LAMP1), CX3CR1 and CD8A as well as TF RUNX3 and EOMES (Figure 5c–d). T-bet was not up regulated in the transcriptome dataset, in line with reported lack of sensitivity of the Illumina array for this gene (Dimova et al., 2015). Expression of IFNG, TNF and
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Transcriptional program underlying the expression of perforin in Th1 cells. The transcriptome of CD4 T cell subsets from five donors was analyzed by gene expression arrays. (a) Log2 expression values of 13,551 probes with a variance >0.01 corresponding to 10,669 unique genes were submitted to unsupervised cluster analysis. See also Supplementary Figure 4. (b) Genes differentially expressed by naive CD4 (CD4<sub>N</sub>), central memory Figure 5 continued on next page</figcaption></figure>
Figure 5 continued

(CM) Th1 and CD4<sub>CTX</sub> T cells were identified. (c-d) Heatmaps of mean log2 expression values of all probes for each selected effector molecules and receptors (c) and transcription factors (d) upregulated in CD4<sub>CTX</sub> T cells as compared to CM<sub>Th1</sub> cells. See also Figure 5—figure supplement 1. (e-f) Heatmaps of methylation beta values of all the probes located in the promoter region of effector molecules and receptors (e) and transcription factors (f) significantly hypomethylated in CD4<sub>CTX</sub> as compared to CM<sub>Th1</sub> as assessed by a two-way ANOVA with Tukey's multiple comparisons (p<0.01). See also Figure 5—figure supplement 2 and Source data file.

DOI: https://doi.org/10.7554/eLife.30496.014

The following source data and figure supplements are available for figure 5:

**Source data 1.** Transcriptional program underlying the expression of perforin in Th1 cells.

DOI: https://doi.org/10.7554/eLife.30496.017

**Figure supplement 1.** Transcriptional program of CM<sub>Th1</sub> and CD4<sub>CTX</sub> T cells.

DOI: https://doi.org/10.7554/eLife.30496.015

**Figure supplement 2.** Epigenetic and transcriptional program underlying the expression of perforin in Th1 cells.

DOI: https://doi.org/10.7554/eLife.30496.016

GZMK was increased in CM<sub>Th1</sub> as compared to naive cells but was not further up regulated in CD4<sub>CTX</sub> T cells (Supplementary file 1).

Additional TF were upregulated in CD4<sub>CTX</sub> T cells (Figure 5d). Most of them were selectively upregulated in CD4<sub>CTX</sub> T cells as compared to CM<sub>Th1</sub> cells or were up regulated from naive T cells to CM<sub>Th1</sub> cells and from CM<sub>Th1</sub> cells to CD4<sub>CTX</sub> T cells. Some TF were first downregulated in CM<sub>Th1</sub> as compared to naive cells and were up regulated in CD4<sub>CTX</sub> T cells as compared to CM<sub>Th1</sub> cells, indicating a more complex regulation pathway (Figure 5—figure supplement 1 and Supplementary file 1). Methylome analysis indicated that gene up regulation in CM<sub>Th1</sub> and CD4<sub>CTX</sub> T cells was primarily associated with DNA hypomethylation, as observed in the naive to CD4<sub>CTX</sub> T cell differentiation (Figure 5—figure supplement 2a). Analysis of probes located in the promoter region of effector molecules and receptors and of TF upregulated in CD4<sub>CTX</sub> cells were also hypo-methylated in these cells (Figure 5e–f and Figure 5—source data 1).

Among the upregulated TF identified through transcriptomic analysis, ASCL2, HOPX, ZEB2 and ZNF683 (Hobit) were selected for further analyses because of the high magnitude of their expression in CD4<sub>CTX</sub> T cells (Figure 5—figure supplement 2b and Supplementary file 1) and because they have been previously studied for their effector functions in T cells (Liu et al., 2014; Albrecht et al., 2010; Dominguez et al., 2015; van Gisbergen et al., 2012; Vieira Braga et al., 2015a2015a).

## TF controlling the expression of perforin in CD4 T cells

The role of selected TF in the acquisition of cytotoxic function by CD4 T cells was analyzed in a novel in vitro model of CD4<sub>CTX</sub> T cell differentiation involving the stimulation of naive CD4 T cells in the presence of Th1, Th2 or Th17 polarizing cytokines and the analysis of perforin and granzyme B expression at multiple time points (Figure 6a). In line with our in vivo observations, cytotoxicity was acquired under Th1 and not under Th2 or Th17 culture conditions (Figure 6b to e and Figure 6—source data 1). In vitro differentiated Th1 cells included a diverse repertoire of clonotypes with similar proportions of large and intermediate expansions as compared to Th1 and Th2 cells (Figure 6—figure supplement 1a). Furthermore, Th1 cells expressed low levels of the PLZF TF (Figure 6—figure supplement 1b) indicating that they were conventional and not innate-type effector T lymphocytes.

Frequencies of perforin<sup>+</sup> granzyme B<sup>+</sup> cells and levels of PRF1 mRNA increased from day 7 to day 21 in Th1 polarized cells (Figure 6b–c–e and Figure 6—source data 1). GNLY and TNF mRNA levels increased following a similar kinetics, whereas GZMB, IFNG and GZMK already reached maximum levels at day 7 (Figure 6e, Figure 6—figure supplement 2, and Figure 6—source data 1). CX3CR1 mRNA presented a unique expression profile with an initial upregulation followed by a marked downregulation (Figure 6—figure supplement 2a). At day 14, Th1 polarized cells acquired potent Concanamycin A-sensitive cytotoxic activity (Figure 6d and Figure 6—source data 1).

Most selected TF were induced by cell activation but their pattern of expression was differently associated with polarizing conditions (Figure 6e and Figure 6—source data 1). TBX21 and HOPX expression was specifically induced under Th1 conditions. RUNX3 and ZBTB7B mRNA were upregulated in all conditions but reached higher levels in Th1 cells. EOMES was induced in both Th1 and
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Naive CD4 T cells differentiate into CD4<sub>CTX</sub> in Th1 culture conditions in vitro. (a) Naive CD4 T cells were stimulated polyclonally in the presence of Th1, Th2 or Th17 polarizing cytokines. Flow cytometry (FACS), gene expression (RNA) and cytotoxicity analyses were performed at the indicated time points. (b-c) Perforin (P) and GranzymeB (G<sub>B</sub>) expression was assessed by flow cytometry. (d) Dot plots (Log<sub>10</sub> fluorescence of one Perforin<sup>+</sup> and Perforin<sup>-</sup> cells) are shown for each time point.

*Figure 6 continued on next page*</figcaption></figure>
Figure 6 continued

representative out of six different donors. Numbers indicate % of cells in individual quadrants. (c) Mean ±SEM of % of perforin (P) and granzyme B (G<sub>B</sub>) positive cells from six independent experiments on different donors (only 4 for day 21 and 5 for day 3). (d) Cytotoxic activity of in vitro differentiated effector T cells against anti-CD3-loaded target cells was assessed at indicated effector/target ratios with or without pre-incubation with Concanamycin A. Figure shows one representative out of four experiments on different donors. (e) mRNA expression of indicated genes was quantified by qPCR. Results are mean ±SEM fold change as compared to naive CD4 T cells from six independent experiments on different donors (only 3 for day 3). *: p<0.05 and **:p<0.01 as compared to Th1 condition at the corresponding time point. (f) Expression of TF was measured by flow cytometry in perforin<sup>high</sup> and perforin<sup>low</sup> Th1 cells after 14 days of in vitro stimulation. Histograms of one representative out of six experiments. Numbers are median fluorescence intensity (MFI) of six experiments (only four for ThPOK). See also **Figure 6—figure supplements 1 and 2** and Source data file.

DOI: https://doi.org/10.7554/eLife.30496.018

The following source data and figure supplements are available for figure 6:

**Source data 1.** Naive CD4 T cells differentiate into CD4<sub>CTX</sub> in Th1 culture conditions in vitro.

DOI: https://doi.org/10.7554/eLife.30496.021

**Figure supplement 1.** TCR repertoire of and PLZF expression by in vitro differentiated CD4 T cells.

DOI: https://doi.org/10.7554/eLife.30496.019

**Figure supplement 2.** Expression of TF and effector molecules by in vitro differentiated CD4 T cells.

DOI: https://doi.org/10.7554/eLife.30496.020

Th2 cells, suggesting that expression of this TF is not sufficient to promote cytotoxic function in CD4 T cells. The pattern of expression of ASCL2 and ZEB2 was inconsistent and ZNF683 mRNA was not detected in in vitro differentiated effector T cells (Figure 6—figure supplement 2a). Flow cytometry analysis of Th1 cells indicated that the expression of perforin was correlated with those of T-bet, Eomes, Runx3 and ThPOK (Figure 6f and Figure 6—figure supplement 2b).

Based on these results, we further evaluated the role of Runx3, T-bet, Eomes, Hopx and ThPOK in this model using shRNA silencing (Figure 7a). Significant knockdown was achieved for each target TF at mRNA and protein levels (Figure 7—figure supplement 1a), whereas the non-silencing (N-S) shRNA and the empty vector (EV) had no significant effect on any of the studied genes (Figure 7—figure supplement 1b–c). Knockdown of Runx3 and, to a lesser extend, T-bet resulted in the decreased expression of perforin and granzyme B. Knockdown of Hopx affected GZMB but not PRF1 mRNA expression, whereas knockdown of Eomes had no significant effect on the expression of either perforin or granzyme B (Figure 7b–d and Figure 7—source data 1). Of note, IFNG expression was principally dependent on T-bet with no significant effect of Runx3 knockdown (Figure 7b). T-bet also controlled TNF, CX3CR1 and GZMK expression whereas GNLY was controlled by both Runx3 and T-bet (Figure 7—figure supplement 1b). Significant interactions were observed between TF, with Runx3 and T-bet controlling the expression of HOPX and T-bet and Hopx controlling the expression of EOMES (Figure 7—figure supplement 1c). In line with its central influence on the expression of cytotoxic molecules, Runx3 also significantly modulated the acquisition of cytotoxic activity by CD4 T cells, whereas T-bet knockdown reduced cytotoxicity of CD4 T cells only at low E/T ratios (Figure 7e and Figure 7—source data 1). ThPOK globally acted as a negative regulator of the cytotoxic program (Figure 7f–g, Figure 7—figure supplement 1b–c and Figure 7—source data 1). ThPOK knockdown markedly increased PRF1 and GZMB expression and cytotoxic activity of CD4 T cells and also upregulated the expression of all the other TF studied and of CD8A.

Together, these results indicate that the acquisition of a cytotoxic program by naive CD4 T cells is dependent on Runx3 and, to a lesser extend, T-bet and is limited by the sustained expression of the CD4 lineage TF ThPOK. Because ThPOK was upregulated in CD8<sub>CTX</sub> T cells in vivo (Figure 2d–e and Figure 2—source data 1), we studied its expression and role in the in vitro differentiation of CD8 T cells. In vitro activation of naive CD8 T cells induced the differentiation of perforin<sup>+</sup> cells that co-expressed ThPOK (Figure 7—figure supplement 2a). Knockdown of ThPOK in differentiated CD8 T cells did not significantly influence the expression of perforin but significantly increased the expression of granzyme B as compared to N-S shRNA (Figure 7—figure supplement 2b), suggesting that ThPOK may limit the cytotoxic function of human CD8 T cells.
Discussion

This study demonstrates that the acquisition of cytotoxic function by human CD4 T cells is an integral part of the Th1 linear differentiation pathway. Several concordant observations support this conclusion. First, *PRF1* gene expression was detected in all subsets of memory Th1 cells. The proportion of *PRF1* expressing cells increased from CM to terminally differentiated EM Th1 cells and this process was associated with the diversification of co-expressed TF networks. The epigenetic modifications detected at the *PRF1* gene promoter reflect this increase in perforin expression and suggest that the local chromatin environment becomes progressively more favorable from naive T cells to CM Th1 cells to terminally differentiated EM Th1 cells. Progressive acquisition of specific transcriptional and epigenetic marks is a hallmark of the linear differentiation model of CD4 T cell memory development (*DEEP Consortium et al., 2016*). Our work suggests that acquisition of cytotoxic functions by CD4 T cells follows a similar stepwise program.

The accessibility and expression of PRF1 gene in all subsets of Th1 cells also provide a basis for the classical observation that memory Th1 cells of diverse antigen-specificities that do not express the perforin protein ex vivo become perforin positive and cytotoxic upon in vitro expansion and cloning (Parronchi et al., 1992; Riaz et al., 2016). This in vitro acquisition of perforin expression likely reflects pre-established PRF1 chromatin modifications in precursors of CD4<sub>CTX</sub> T cells.

The cytotoxic potential of Th1 cells is also supported by the in vitro model of CD4<sub>CTX</sub> T cell differentiation. In this model, human naive CD4 T cells stimulated in the presence of Th1, and not Th2 or Th17, polarizing cytokines differentiated in perforin<sup>+</sup>granzyme B<sup>+</sup> cells with potent cytotoxic activity. This could be due to the promotion of the cytotoxic phenotype by Th1 cytokines or to its repression by the Th2 or Th17 transcriptional programs, or both (Parronchi et al., 1992; Xiong et al., 2013; Ciucci et al., 2017). In contrast to the rapid acquisition of IFNG expression, naive CD4 T cells expressed high levels of PRF1<sup>+</sup> to 2 weeks after stimulation. This observation also contrasts with the more rapid acquisition of cytotoxic function by CD8 T lymphocytes (Araki et al., 2008) and suggests that the initial steps of Th1 cell differentiation provide the required epigenetic and transcriptional signals promoting the expression of PRF1. The relatively delayed acquisition of cytotoxicity also suggests that CD4<sub>CTX</sub> T cells are induced by prolonged antigen stimulation in vivo and may intervene when Th1 and CD8 T cells do not adequately control pathogens.

The acquisition of cytotoxic function within the Th1 lineage was promoted by Runx3 and T-bet. Runx3 knockdown reduced the expression of perforin, granzyme B and granulysin by CD4<sub>CTX</sub> and decreased their cytotoxic activity. In contrast to mouse Th1 cells, Runx3 did not influence IFN-γ mRNA expression by human Th1 cells (Djuretic et al., 2007; Wang et al., 2014). This discrepancy may be related to inter-species differences or to incompleteness of the knockdown in our experimental conditions. T-bet knockdown reduced the expression of IFN-γ and cytotoxic molecules but had a more moderate impact than Runx3 knockdown on the acquisition of cytotoxic activity. This result is in line with a recent report indicating a role of T-bet in the induction of cytotoxic molecules by TCR-engineered tumor-specific CD4 T cell lines (Jha et al., 2015). Together these results indicate that human Th1 cells acquire cytotoxic functions under the control of the master regulator Runx3 cooperating with T-bet (Cruz-Guilloty et al., 2009; Jha et al., 2015). Runx3 and T-bet also controlled the expression of other TF that were not directly involved in the acquisition of cytotoxic function in these experimental conditions, including Eomes and Hopx.

Eomes was upregulated in CD4<sub>CTX</sub> T cells as compared to naive cells in vivo but its expression was less correlated with perforin than Runx3 and T-bet. In vitro, Eomes was neither sufficient nor necessary to induce cytotoxicity as it was upregulated in Th2 cells that did not express perforin and its knockdown did not impact the expression of perforin in Th1 cells. Together, these results suggest that the role of Eomes in the acquisition of cytotoxic function by human CD4 T cells may be limited. This contrasts with a report indicating that overexpression of Eomes induces cytotoxic function in CD4 T cell lines, suggesting that this TF may promote cytotoxicity when expressed at high levels in CD4 T cells, as observed in CD8 T cells (Pearce et al., 2003; Eshima et al., 2012; Intlekofer et al., 2008). The role of Eomes could also be restricted to specific conditions of co-stimulation as recently reported (Mittal et al., 2018). The expression of Hopx by CD4<sub>CTX</sub> T cells showed a similar pattern as T-bet and Runx3 both in vivo and in vitro. Its knockdown reduced the expression of granzyme B but did not significantly impact perforin. Therefore, the primary role of Hopx in CD4<sub>CTX</sub> T cells may not be the induction of cytotoxicity but may include other functions, including survival (Albrecht et al., 2010).
Two other TF, Hobit (ZNF683) and ZEB2, were specifically expressed by CD4<sub>CTX</sub> T cells in vivo but not in vitro. Hobit is upregulated in human effector CD8 T cells and murine NKT cells and therefore appears to be part of a signature common to cytotoxic lymphoid cells (van Gisbergen et al., 2012; Vieira Braga et al., 2015a, 2015b; Vieira Braga et al., 2015b). A recent study revealed that Hobit induces a transcriptional program promoting tissue residency of memory T cells and suggests that it could regulate cytotoxic function in murine NKT1 cells (Mackay et al., 2016). Its role in the promotion of cytotoxic function by CD4 T cells in vivo therefore remains to be established. Similarly, further studies should establish the role of ZEB2 in terminally differentiated CD4<sup>+</sup> T lymphocytes (Dominguez et al., 2015; Omilusik et al., 2015).

In parallel with the positive regulation operated by T-bet and Runx3, the cytotoxic function of human CD4 T cells was negatively regulated by ThPOK. In CD4<sup>+</sup>CD8<sup>lo</sup> thymocytes, ThPOK decreases the expression of Runx3 and perforin and promotes the development of the CD4 T cell lineage (Liu et al., 2005). In the periphery, ThPOK also inhibits the expression of Runx3, T-bet and Eomes by mouse CD4 T cells and restricts their cytotoxic differentiation (Wang et al., 2008). We observed that the expression of Runx3 and ThPOK is not mutually exclusive in human CD4<sub>CTX</sub> and CD8<sub>CTX</sub> T cells. Co-expression of Runx3 and ThPOK has been observed in murine Th1 cells (Djuretic et al., 2007) and in simian MHCII-restricted CD8<sub>αα</sub> T cells after CD4 downregulation (Vinton et al., 2017). Also, ThPOK is up regulated by mouse effector CD8 T cells during acute viral infection and promotes their expansion and effector function upon rechallenge (Setoguchi et al., 2009). On the other hand, ThPOK is required for the development of other murine lymphoid subsets with cytotoxic potential including CD4<sup>+</sup> NKT cells and γδ T cells (Wang et al., 2010; Park et al., 2010). Together, these observations indicate that the expression of high levels of ThPOK is compatible with the expression of Runx3 and with cytotoxic function. Yet, the in vitro model of T cell differentiation revealed that ThPOK is an important regulator of the cytotoxic activity of human CD4 and possibly CD8 T lymphocytes. This observation suggests that ThPOK may contribute to the regulation of antiviral, anticancer and immunopathological properties of Th1 cells in vivo.

In conclusion, this study shows that Runx3 and ThPOK cross-regulate the acquisition of cytotoxic function by Th1 lymphocytes and therefore represent targets for interventions against viral infections, cancer and autoimmune disorders.

## Materials and methods

### Key resources table

<table><thead><tr><th>Reagent type (species) or resource</th><th>Designation</th><th>Source or reference</th><th>Identifiers</th></tr></thead><tbody><tr><td>Cell line (mouse leukemia)</td><td>RAW 264.7</td><td>ATCC Cat# TIB-71</td><td>RRID:CVCL_0493</td></tr><tr><td>Cell line (human kidney cell line)</td><td>HEK-293 human kidney cell line</td><td>ATCC Cat# CRL-1573,</td><td>RRID:CVCL_0045</td></tr><tr><td>Cell line (human)</td><td>HEL 299</td><td>ATCC Cat# CCL-137,</td><td>RRID:CVCL_2480</td></tr><tr><td>Transfected constructs</td><td colspan="3">References of all shRNAs are listed in <i>Supplementary file 2f</i></td></tr><tr><td>Antibody-ChIP</td><td colspan="3">References of used antibodies are indicated in the method section</td></tr><tr><td>Antibody-cytometry</td><td colspan="3">References of all used cytometry antibodies, including company and clone are listed in <i>Supplementary file 2b</i>.</td></tr><tr><td>Sequence-based reagent</td><td colspan="3">All sequences are listed in <i>Supplementary files 2c, d and e</i></td></tr><tr><td>Chemical compound, drug</td><td>Concanamycin A</td><td>Sigma-Aldrich (Merck, Germany)</td><td>C 9705</td></tr><tr><td>Other</td><td>PKH26 staining</td><td>Sigma-Aldrich (Merck, Germany)</td><td>Catalog numbers MINI26 and PKH26GL</td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** TF controlling the expression of perforin in CD4 T cells. (a) Naive CD4 T cells were stimulated polyclonally in the presence of Th1 polarizing cytokines and were transduced on day 2 with shRNA-GFP. GFP+ cells were purified on day 14 for flow cytometry, mRNA expression and cytotoxicity analyses. (b) Expression of PRF1, GZMB and IFNG mRNA was quantified by qPCR in cells transduced with indicated targeting shRNA, non-silencing (N- Figure 7 continued on next page</figcaption></figure>
Figure 7 continued

S) shRNA or empty vector (EV). Data are mean ± SD log2 fold change as compared to controls (mean of N-S and EV), from seven biological replicates generated in four independent experiments on different donors. (c-d) Expression of perforin (P) and granzyme B (G<sub>B</sub>) was measured by flow cytometry in cells transduced with N-S or indicated gene targeting shRNA. (c) Representative dot plot (log10 fluorescence) of the proportions of double positive cells from seven independent experiments on different donors. Numbers indicate % of cells in individual quadrants. (d) Perforin expression following knockdown of indicated TF. Data are median fluorescence intensity (MFI) normalized for perforin MFI in cells transduced with N-S shRNA (dotted line, 100%). Experiments in which gene knockdown was below 10% were excluded from the analysis. (e) Cytotoxic activity of transduced cells against anti-CD3-loaded target cells. Data are mean ± SEM of five independent experiments on different donors. *:p<0.05 as compared to N-S shRNA. (f-g) The effect of ThPOK knockdown on the expression of indicated TF, CD8A and CD4 was studied by qPCR and flow cytometry. (f) Data are mean ± SD log2 fold change as compared to mean RNA expression in control shRNA (N-S and EV) from seven biological replicates generated in four independent experiments on different donors. *:p<0.05; **:p<0.01 as compared to control shRNA. (g) Histograms (Log10 fluorescence) from one representative donor. Data are MFI normalized for TF MFI in cells transduced with N-S shRNA (dotted line, 100%) of seven independent experiments on different donors. See also **Figure 7—figure supplements 1 and 2** and Source data file.

DOI: https://doi.org/10.7554/eLife.30496.022

The following source data and figure supplements are available for figure 7:

**Source data 1.** TF controlling the expression of perforin in CD4 T cells.

DOI: https://doi.org/10.7554/eLife.30496.025

**Figure supplement 1.** Knockdown of transcription factors (TF) expression by shRNA.

DOI: https://doi.org/10.7554/eLife.30496.023

**Figure supplement 2.** Influence of ThPOK on the differentiation of cytotoxic CD8 T cells.

DOI: https://doi.org/10.7554/eLife.30496.024

## Cell collection and purification

Blood samples were collected from CMV-seropositive or seronegative healthy adult volunteers, aged 26 to 60 (median (interquartile range): 45 (39.5–51.25)) years. Volunteers were recruited by the research centre ImmuneHealth, CHU Tivoli, La Louvière. Clinical staff informed the volunteers about the objectives of the study and obtained their written consent to use the human material for research purposes. The study and the informed consent form were approved by the Ethics committee of the CHU Tivoli, La Louvière, Belgium (Reference B09620097253). The study followed the Good Clinical Practice (ICH/GCP) guidelines, the Belgian Law and the declaration of Helsinki ('World Medical Association Declaration of Helsinki; Ethical Principles for Medical Research Involving Human Subjects'). Transcriptomic and methylomic analyses were conducted on 3 CMV-seropositive women, aged 38, 48 and 60 years, and 2 CMV-seropositive men, aged 40 and 52 years.

The number of samples analyzed in each experiment was defined on the basis of previous experience of the investigators or on published literature. No sample size was computed. Peripheral blood mononuclear cells (PBMC) were purified by density gradient centrifugation and stained with titrated conjugated antibodies. Cells were sorted on a BD FACS Aria III or acquired with a BD LSR Fortessa cytometer and data were analyzed using the FlowJo software (v9.2). The RAW 264.7 murine macrophage cell line (RRID:CVCL_0493; obtained from ATCC) and the HEK-293 human kidney cell line (RRID:CVCL_0045, obtained from ATCC) were cultured in DMEM (Lonza) supplemented with 10% fetal calf serum, 1% AAG, 1% Na Pyruvate and 1% Pen/strep. The HEL-299 fibroblastic cell line (RRID:CVCL_2480; obtained from ATCC) was cultured in EMEM (Lonza) supplemented with 10% fetal calf serum, 1% NEAA, 1% Hepes, 1% Glutamine, 1% Na Pyruvate and 1% Pen/strep. All cell lines were tested negative for mycoplasma infection (MycoAlert, Lonza). Because they were used only as tools to produce lentivirus particles, as targets of cytotoxic cells, and as a negative control in one methylation analysis, they were not re-authenticated after purchase.

## FACS-staining

For membrane staining, cells were washed with PBS containing 0.1% bovine serum albumin (BSA). Antibodies were incubated in PBS + 0.1% BSA for 10 min at 37°C or 15 min at room temperature. Cells were then washed with PBS + 0.1% BSA before addition of Cellfix (BD) or intracellular staining. Cells were permeabilized for intranuclear and intravesicular staining using the FoxP3 staining kit (eBiosciences). Active caspase3 staining was performed using the cytofix-cytoperm and Permwash buffers (BD). References of used antibodies are presented as **Supplementary file 2b**.
## FACS-sorting

Naive CD4 T cells were isolated by negative selection for in vitro stimulation. Before cell sorting, fresh PBMC were enriched in CD4 T cells with the Miltenyi human CD4 +T Cell Isolation Kit. Membrane staining was then performed as mentioned above with a dump channel in PE including CD14, CD19, CD16, CD56, TCRgd and CD8 mAbs. CD25-, CD45RO- and CXCR3-negative cells were further selected in order to exclude regulatory T cells, memory and stem cell memory CD4 T cells, respectively. A small fraction of these untouched naive CD4 T cells (CD3⁺CD4⁺CD45RO⁻CCR7⁺-CD28⁺) cells were then stained to verify their naive phenotype. Cell purity was 96 [94-97]% (Median [IQ]). CD4 and CD8 T cell subsets were isolated by positive selection. Before cell sorting, fresh PBMC were depleted of glycophorin A-, CD19- and CD14-positive cells as well as CD8- or CD4-positive cells using an Automacs instrument (Miltenyi). Membrane staining was then performed as mentioned above. Cells were resuspended in complete antibiotics-containing medium, sorted and collected in the same complete medium, centrifuged and lysed in RLT Plus buffer +10 µl betamercaptoethanol for later nucleic acid extraction or immediately tested for cytotoxic activity. Cell purity was 98 [95-99]% (Median [IQ]). Single cells were sorted on a FACS Aria III cell sorter (BD) following staining and suspension in an EDTA-containing sorting buffer. Quality of sorting was assessed using the staining index from the DIVA software version 8.0.

## Cytotoxicity assay

Effector cells were pre-incubated for 1 hr with or without 100 nM Concanamycin A (Sigma-Aldrich-Merck, Germany). RAW target cells were labelled with PKH-26 (Sigma-Aldrich-Merck, Germany) as previously described (Sheehy et al., 2001). Effector cells were added to 5000 RAW cells at appropriate effector/target ratios and incubated for 4 hr in the presence of 2 µg/ml mouse anti-human CD3 antibody (clone OKT3). RAW cells incubated with the anti-CD3 antibody but without effector cells were used as controls. Percentage of lysis was calculated as the percentage of caspase3-positive RAW cells after subtraction of the % of active caspase3 in the control wells (He et al., 2005).

## Quantitative PCR

qPCR was performed using the Taqman RNA Amplification kit or the LightCycler Multiplex RNA Virus Master and a LightCycler 480 instrument (Roche). Raw data were analyzed using the fit points method and fold change was calculated with the Delta-Delta Cp method using the housekeeping gene *EEF1A1* (*EF1*) as a reference. Primers and fluorescent probes were designed using Primer3 and purchased from Eurogentec. A Taqman Gene Expression assay was used for *RUNX3* (Hs00231709_m1, Thermo Fisher) and *PLZF* (ZBTB16, Hs00232313_m1, Thermo Fisher) analyses. Oligonucleotide sequences are presented as *Supplementary file 2c*.

## Bisulphite pyrosequencing

The *PRF1* promoter was divided into 11 amplicons covering 34 CpG sites as previously described (Narasimhan et al., 2009). Genomic DNA was bisulphite-converted using the Qiagen FAST Epitect bisulfite kit and sequenced using a Pyromark Q96 device after isolation of single strand biotinylated DNA from the PCR product using streptavidin and a pyromark Q96 vacuum prep station. Quality assessment and methylation level calculation were performed using the software Pyro Q-CpG and the CpG assay 1.0.9 (Biotage).

## Nucleic acid material

DNA and RNA were extracted using the Qiagen AllPrep DNA/RNA kit. Concentration and purity were assessed by spectrophotometry (nanodrop - Thermoscientific). Median [interquartile range] of A260/A280 ratios were 1.93 [1.77-2.10] and 1.82 (1.71-1.91) for RNA and DNA samples, respectively. For microarray analyses, RNA integrity number (RIN) was measured using the Eukaryote Total RNA Nano assay and a Bioanalyzer (Agilent). One sample out of 25 had a RIN below seven and was excluded from the analyses.

## Gene expression and methylation microarrays

Total RNA was amplified with the Illumina TotalPrep RNA Amplification Kit (Ambion) and hybridized with the HumanHT-12 v4 array containing 47,323 probes for 44,053 annotated genes, according to
the instructions of the Whole-Genome Gene Expression Direct Hybridization Assay (Illumina). Chips were scanned with the HiScan Reader (Illumina). For methylation analyses, genomic DNA was bisulphite-converted using an EZ DNA methylation Kit (ZYMO). DNA methylation level was measured using the Illumina Infinium HD Methylation Assay. Bisulphite converted DNA was hybridized with the Illumina HumanMethylation450 BeadChip 450K array (12 samples/chip), as described previously (Dedeurwaerder et al., 2011). Data from both arrays are available on GEO (https://www.ncbi.nlm.nih.gov/geo/) under the accession number GSE75406.

## Chromatin immunoprecipitation (ChIP) and ChIP-qPCR

MACS-purified CD4 T cells were stained and fixed with 1% formaldehyde. Glycine was added at a final concentration of 0.125 M to quench the crosslinking reaction. Cells were washed twice with ice-cold PBS and resuspended in complete medium for FACS sorting. Dry pellets of sorted cells were frozen at -80°C. Thawed pellets were lysed in 1% SDS-containing buffer and sonicated to obtain DNA fragments of 300-800 bp using a Bioruptor device (Diagenode). Chromatin of 200,000 cells was immunoprecipitated with an anti-histone antibody and protein G magnetic-activated beads. Chromatin was incubated overnight at 4°C with the following antibodies: 1 µg anti-H3K4me3 (Millipore 17-614 rabbit monoclonal), 0.5 µg anti-H3K27ac (abcam ab4729 rabbit polyclonal) or anti-H3 (diagenode C15310135 rabbit polyclonal). 1% of the IP reaction was collected before adding the antibody and the beads and was used as total chromatin input. Beads were washed five times: once with low-salt buffer, once with high-salt buffer, once with lithium chloride containing buffer and twice with TE buffer. After washing and reverse crosslinking (incubation with NaCl 200 mM for 4 hr at 65°C), chromatin was eluted with a buffer containing 1% SDS and 100 mM NaHCO₃ and treated with RNAse A and Proteinase K for 1 hr at 45°C. IP-DNA was purified using the MinElute PCR purification kit (Qiagen) and then analyzed by qPCR using the Probe Master 480 kit with primers encompassing different regulatory regions of the perforin locus (sequences are presented as *Supplementary file 2d*). The DeltaCp method was used to calculate the % of input for each IP. Results were normalized for the DeltaCp of H3-IP DNA.

## Singe-cell qPCR assay

Single cells were collected in 5 µL lysis buffer (CelluLyser micro lysis buffer from Tataa biocenter), immediately frozen on dry ice and stored at -80°C until used (Svec et al., 2013). Reverse transcription (RT) was performed using the CelluLyser Micro Lysis and cDNA Synthesis Kit following manufacturer's instructions (Tataa biocenter). RT step was validated using the Universal RNA Spike (TATAA Universal RNA Spike I from Tataa biocenter) in each well to ensure the absence of inhibitions. Wells most likely to contain single cells were selected on the basis of housekeeping gene expression and exclusion of outliers. cDNA was then pre-amplified for 20 cycles using the TATAA PreAmp Grand-Master Mix kit from the same company. Single-cell qPCR was performed on 43 cells per subset in 96.96 Dynamic Array IFC plates for Gene Expression (BMK-M-96.96) using the fluidigm technology (Biomark HD). Primer sequences are presented in *Supplementary file 2e*.

## In vitro activation and polarisation of naive CD4 T cells

Ex vivo isolated naive CD4 T cells were activated with 0.5 µg/ml phytohemagglutinin (PHA) in the presence of allogeneic CD4-depleted irradiated PBMC used as feeders at the ratio of 1/1. Cells were then cultured in the presence of IL-2 (R and D, 6 ng/ml), IL-15 (R and D, 5 ng/ml) and Th1, Th2 or Th17 polarising cytokines for 3 to 21 days. Culture medium was RPMI 1640 (Lonza) supplemented with 10% Fetal calf serum, 1% amino acids and glutamine, 1% Penicillin/Streptomycin. All incubation steps were performed at 37°C with 5% CO₂. IL-2 and IL-15 were added on day 2 and medium was replenished when required on the basis of cell proliferation. Polarizing cytokines and neutralizing antibodies were purchased from eBiosciences and used at final concentration of 10 ng/ml and 10 µg/ml, respectively. Th1 polarizing medium contained IFNγ, IL-12, anti-IL-4 (clone MP4-25D2), and anti-IL-17 (clone eBio64CAP17). Th2 polarizing medium contained IL-4, anti-IL-17, anti-IFNγ (clone NIB42), and anti-IL-12 (clone 20C2). Th17 polarizing medium contained IL-1β, IL-6, IL-23, TGFβ, anti-IL-4, anti-IFNγ, and anti-IL-12. On day 3, CD4 T cells were isolated by MACS-positive selection before nucleic acid extraction. On days 7, 14 and 21, cells were collected for downstream analyses.
## Transcription factor knockdown

Lentiviral particles were produced by transient transfection of packaging HEK 293 T cells. pMD2.G and psPAX2 were used as envelop and core packaging plasmids, respectively, together with the gene transfer plasmid (Supplementary file 2f). Before transduction, 50,000 freshly isolated naive CD4 T cells were stimulated during 46 hr in the presence of feeders and PHA (5 µg/ml) in Th1-cytokines containing X-Vivo15 medium (Lonza). For transduction, viral particles were added at a MOI of 10, together with IL-2 and IL-15 in 50 µl of fresh medium. Cells were amplified during 10 days before sorting of GFP+ transduced cells for down-stream analysis.

## TCR CDR3 sequencing

Purified cDNA (AMPure XP Beads (Agencourt)) was obtained from total RNA and used in template-switch anchored RT-PCR experiment with specific alpha and beta chain primers. PCR products were then submitted to high-throughput sequencing as previously described (Van Caeneghem et al., 2017). Briefly, V2 300 kit with 200 bp at the 3' end (read 2) and 100 bp at the 5' end (read 1) were used on the Illumina MiSeq platform.

## Statistical analyses

Data were analyzed with the GraphPad Prism software unless otherwise specified. After one-way ANOVA analysis of variance, a Mann-Whitney test was performed for selected two-by-two comparisons and a Dunnet's test for multiple comparisons when appropriate. For grouped analysis, we used two-way ANOVA with multiple Tukey's tests. Differences were considered statistically significant at p-values<0.05.

## Illumina Expression HT12 Arrays

Raw data were quantile normalized using the normalization method from the lumi package (Du et al., 2008). Unsupervised clustering (Uc) analysis of gene expression datasets was performed using the pvclust package of the R software (R)(Suzuki and Shimodaira, 2006). The robustness of the Uc tree was tested by multiscale bootstrap resampling using Pearson's correlation as distance and Ward.D2 as clustering method, with 1000 iterations. An AU (approximately unbiased) p-value (percentage) was calculated and placed on the nodes of the cluster dendrogram. Principal component analysis (PCA) on the expression dataset was performed using the MultiExperiment Viewer (MeV) tool and the scatter plot function in R. The GeneSign module of the BubbleGUM software (Spinelli et al., 2015) was used with the Min/max method to identify lists of genes specifically expressed in cell subsets. A probe was considered specific of a given cell subset if its minimal expression value across the replicates of the cell subset of interest (test population) was higher than its maximal expression value across the replicates of the cell population used as reference (reference population). To obtain a limited number of genes, probes for which the ratio between the maximal and minimal expression values across all samples was below 1.2 were considered not regulated in any cell subsets and thus excluded from the analysis. Finally, probe lists were transformed into gene lists (or GeneSets). Heatmaps were generated using the heatmap.2 function of gplots package of R.

## Illumina HumanMethylation450 BeadChip arrays

Raw data were filtered using a detection p-value<0.05. Cross-reactive probes were filtered out while probes containing SNPs, which do not introduce an important confounder in intra-individual studies, were kept in the analysis as previously detailed. (Dedeurwaerder et al., 2014) β-values were computed using the following formula: $β-value = M/[U + M]$ where M and U are the raw 'methylated' and 'unmethylated' signals, respectively. The β-values were corrected for type I and type II bias using the peak-based correction (Dedeurwaerder et al., 2011; Dedeurwaerder et al., 2014). The differential analyses were applied according to published recommendations (Dedeurwaerder et al., 2014): first the methylation values were converted to M-values using the following formula: $M-value = log2 (β-value/(1-β-value))$. The statistical significance of the differential methylation was then assessed using a paired t-test applied on these M-values. In parallel, a Delta-β was computed as the absolute difference between the median β-values. Cytosine showing a p-value<1e-4 together with an absolute delta-β>0.1 were reported as differentially methylated. Heatmaps were generated based on the scaled beta values of all the probes located in the promoter region of each
represented gene, using the heatmap.2 function of gplots package of R. Promoter regions included 5'prime, TSS1500, TSS200 and 1st exon.

## Gene Set Enrichment Analysis (GSEA)

GSEA was used to analyze the enrichment of GeneSets obtained by GeneSign on the pairwise comparisons of our gene expression microarray data (Subramanian et al., 2005). Catalog c3 from MsigDB was added to our GeneSets for robust statistical analysis. GSEA was used with 1000 Gene-Set-based permutations and 'difference of classes' as a metric for ranking the genes since the expression values were in Log2 scale. The software quantifies enrichments by computing the Normalized Enrichment Score (NES) and the False Discovery Rate (FDR). FDR below 0.25 was considered significant (Subramanian et al., 2005). Genes showing no differential expression between CD8<sub>CTX</sub> and CD4<sub>CTX</sub> were used as negative control. GSEA was also used to quantify the correlation between gene expression and methylation. The identifiers of the genes included in the GeneSets were transformed to match the DNA methylation probe identifiers (Perl scripts included in Supplementary file 3). GSEA Pre-ranked analysis was then used to assess the enrichment of our expression-based Gene-Sets on the pairwise comparisons of the pre-ranked methylation delta β values. Genesets of catalog c3 from MsigDB were similarly converted to probe identifiers and added to our GeneSet file for robust statistical evaluation.

## Single-cell qPCR

Data mining was performed using the Fluidigm Real-time PCR analysis (V1.4.3), and data analysis was performed using the Genex6 MultiD software as previously described (Ståhlberg et al., 2013). Heat Map was created using the standard function in R.

## TCR repertoire

Raw sequencing reads from fastq files (both reads) were aligned to reference V, D and J genes specifically for 'TRA' or 'TRB' to build CDR3 sequences using the MiXCR software version 1.7. (Bolotin et al., 2015). CDR3 sequences were then assembled and clonotypes were exported and analyzed using the tcR package (Nazarov et al., 2015).

## Acknowledgements

We thank the blood donors and ImmuneHealth for the clinical sample collection. We are grateful to Véronique Olislagers, Sandra Lecomte, Margreet Brouwer, Muriel Nguyen, Séverine Thomas, Inès Vu Duc, Charlotte Givord and lain Welsby for technical assistance, Ariane Huygens for precious expertise in flow cytometry and Charles Dehout, Nicolas Dauby, Muriel Moser, Alain Le Moine and Fabienne Willems for helpful discussions, Sophie Lucas (Institut de Duve – Université Catholique de Louvain) for providing single cell PCR reagents as well as Thierry Voet and Koen Theunis (Katholieke Universiteit Leuven) for single cell PCR analyses. YS is research assistant, SG senior research associate and AM research director at the Fonds de la Recherche Scientifique (FRS – FNRS), Belgium. A.H. was supported by the ASCO grant of the Fonds Erasme. This work was supported by the FRS. – FNRS, by the Belgian Federal Public Planning Service Science Policy and by the European Regional Development Fund (ERDF) and the Walloon Region (Wallonia-Biomed portfolio, 411132–957270).

## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Fonds De La Recherche Scientifique - FNRS</td><td>PhD Student Fellowship</td><td>Yasmina Serroukh</td></tr><tr><td>Belgian Federal Public Planning Service Science Policy</td><td>Research Project Grant</td><td>Stanislas Goriely<br>Arnaud Marchant</td></tr></tbody></table>
<table><tr><td>European Regional Development Fund and Walloon Region</td><td>Research Project Grant (411132-957270)</td><td>Stanislas Goriely Arnaud Marchant</td></tr><tr><td>Fonds De La Recherche Scientifique - FNRS</td><td>Research Project Grant (PDR)</td><td>François Fuks Stanislas Goriely Arnaud Marchant</td></tr><tr><td>Fonds Erasme</td><td>PhD Student Fellowship</td><td>Alice Hoyois</td></tr></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## Author contributions

Yasmina Serroukh, Chunyan Gu-Trantien, Baharak Hooshiar Kashani, Conceptualization, Formal analysis, Investigation, Writing—original draft, Writing—review and editing; Matthieu Defrance, Thien-Phong Vu Manh, Jishnu Das, Martin Bizet, Emeline Pollet, Formal analysis, Writing—review and editing; Abdulkader Azouz, Aurélie Detavernier, Alice Hoyois, Tressy Tabbuso, Emilie Calonne, Investigation, Writing—review and editing; Klaas van Gisbergen, Marc Dalod, François Fuks, Conceptualization, Resources, Writing—review and editing; Stanislas Goriely, Arnaud Marchant, Conceptualization, Supervision, Funding acquisition, Writing—review and editing

## Author ORCIDs

Yasmina Serroukh http://orcid.org/0000-0001-8751-5725

Thien-Phong Vu Manh http://orcid.org/0000-0002-0294-342X

Jishnu Das http://orcid.org/0000-0002-5747-064X

Stanislas Goriely http://orcid.org/0000-0002-7005-6195

Arnaud Marchant http://orcid.org/0000-0003-0578-0467

## Ethics

Human subjects: Volunteers were recruited by the research centre ImmuneHealth, CHU Tivoli, La Louvière. Clinical staff informed the volunteers about the objectives of the study and obtained their written consent to use the human material for research purposes. The study and the informed consent form were approved, following approval by the Ethics committee of the CHU Tivoli, La Louvière, Belgium (Reference B09620097253). The study followed the Good Clinical Practice (ICH/GCP) guidelines, the Belgian Law and the declaration of Helsinki ("World Medical Association Declaration of Helsinki; Ethical Principles for Medical Research Involving Human Subjects").

## Decision letter and Author response

Decision letter https://doi.org/10.7554/eLife.30496.035

Author response https://doi.org/10.7554/eLife.30496.036

## Additional files

### Supplementary files

* Supplementary file 1. Lists of genes included in the 12 GeneSets obtained from CD4 versus CD8 T cell comparison (Figure 2) and naive CD4 T cell versus CMTh1 cell versus CD4CTX T cell comparison (Figure 4). Genes expressed at higher level by the first listed subset as compared to the subset indicated between brackets were identified using the min/max method. Genes are ranked according to mean log2 fold change calculated using the Limma package in R.

DOI: https://doi.org/10.7554/eLife.30496.026

* Supplementary file 2. (a) Level of methylation at individual CpG sites in memory Th1 cell subsets in vivo. The level of methylation was measured by pyrosequencing in memory Th1 cell subsets of two to nine CMV-seropositive healthy adults. Data are median percentage ± interquartile range and were compared with the Mann-Whitney non-parametric test. † CM<sub>Th1</sub> versus naive CD4 T cells. ‡ EM28<sup>+</sup><sub>Th1</sub> versus naive CD4 T cells. § CD4<sub>CTX</sub> versus naive CD4 T cells. ¶ EM28<sup>+</sup><sub>Th1</sub> versus CM<sub>Th1</sub>. || CD4<sub>CTX</sub> versus EM28<sup>+</sup><sub>Th1</sub>. nd : not done, no statistical analysis was performed because of insufficient
number of subjects (n = 2 or 3) ; ns : non significant ;*=p<0.05 ; **=p<0.01 ; ***=p<0.001 ; ****=p<0.0001. (b) Monoclonal antibodies used in flow cytometry experiments. (c) QPCR Oligonucleotide sequences (d) ChIP-PCR oligonucleotide sequences. (e) Single-cell qPCR oligonucleotide sequences. (f) References for the gene transfer plasmids used in transcription factor knockdown experiments.

DOI: https://doi.org/10.7554/eLife.30496.027

* Supplementary file 3. Method to analyze the correlation between gene expression and DNA methylation. Gene Set Enrichment Analysis (GSEA) was used to analyze the enrichment of DNA methylation profiles in the GeneSets expressed by T cell subsets. Five Perl scripts were created to prepare the file required for GSEA.

DOI: https://doi.org/10.7554/eLife.30496.028

* Transparent reporting form

DOI: https://doi.org/10.7554/eLife.30496.029

## Major datasets

The following dataset was generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database, license, and accessibility information</td></tr></thead><tbody><tr><td>Serroukh Y, Kashani BH, Gu-Trantien C, Defrance M, Vu Manh TP, Azouz A, Detavernier A, Das J, Bizet M, Pollet E, Tabbuso T, Ca-lonne E, van Gis-bergen K, Dalod M, Fuks F, Goriely S, Marchant A</td><td>2017</td><td>Genome-wide study of transcriptional and methylation profile of human cytotoxic CD4 T cells</td><td>https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?&acc=GSE75406</td><td>Publicly available at the NCBI Gene Expression Omnibus (accession no: GSE75406)</td></tr></tbody></table>

## References

* Adams D, Altucci L, Antonarakis SE, Ballesteros J, Beck S, Bird A, Bock C, Boehm B. Campo E, Caricasole A, Dahl F, Dermitzakis ET, Enver T, Esteller M, Estivill X, Ferguson-Smith A, Fitzgibbon J, Flicek P, Giehl C, Graf T, et al. 2012. BLUEPRINT to decode the epigenetic signature written in blood. *Nature Biotechnology* **30**:224-226. DOI: https://doi.org/10.1038/nbt.2153, PMID: 22398613

* Albrecht I, Niesner U, Janke M, Menning A, Loddenkemper C, Kühl AA, Lepenies I, Lexberg MH, Westendorf K, Hradilkova K, Grün J, Hamann A, Epstein JA, Chang HD, Tokoyoda K, Radbruch A. 2010. Persistence of effector memory Th1 cells is regulated by Hopx. *European Journal of Immunology* **40**:2993-3006. DOI: https://doi.org/10.1002/eji.201040936, PMID: 21061432

* Appay V, Dunbar PR, Callan M, Klenerman P, Gillespie GM, Papagno L, Ogg GS, King A, Lechner F, Spina CA, Little S, Havlir DV, Richman DD, Gruener N, Pape G, Waters A, Easterbrook P, Salio M, Cerundolo V, McMichael AJ, et al. 2002b. Memory CD8+ T cells vary in differentiation phenotype in different persistent virus infections. *Nature Medicine* **8**:379-385. DOI: https://doi.org/10.1038/nm0402-379, PMID: 11927944

* Appay V, Zaunders JJ, Papagno L, Sutton J, Jaramillo A, Waters A, Easterbrook P, Grey P, Smith D, McMichael AJ, Cooper DA, Rowland-Jones SL, Kelleher AD. 2002a. Characterization of CD4(+) CTLs ex vivo. *The Journal of Immunology* **168**:5954-5958. DOI: https://doi.org/10.4049/jimmunol.168.11.5954, PMID: 12023402

* Araki Y, Fann M, Wersto R, Weng NP. 2008. Histone acetylation facilitates rapid and robust memory CD8 T cell response through differential expression of effector molecules (eomesodermin and its targets: perforin and granzyme B). *The Journal of Immunology* **180**:8102-8108. DOI: https://doi.org/10.4049/jimmunol.180.12.8102, PMID: 18523274

* Bolotin DA, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IV, Putintseva EV, Chudakov DM. 2015. MiXCR: software for comprehensive adaptive immunity profiling. *Nature Methods* **12**:380-381. DOI: https://doi.org/10.1038/nmeth.3364, PMID: 25924071

* Brown DM, Lee S, Garcia-Hernandez ML, Swain SL. 2012. Multifunctional CD4 cells expressing gamma interferon and perforin mediate protection against lethal influenza virus infection. *Journal of Virology* **86**:6792-6803. DOI: https://doi.org/10.1128/JVI.07172-11, PMID: 22491469

* Casazza JP, Betts MR, Price DA, Precopio ML, Ruff LE, Brenchley JM, Hill BJ, Roederer M, Douek DC, Koup RA. 2006. Acquisition of direct antiviral effector functions by CMV-specific CD4+ T lymphocytes with cellular maturation. *The Journal of Experimental Medicine* **203**:2865-2877. DOI: https://doi.org/10.1084/jem.20052246, PMID: 17158960
Cheroutre H, Husain MM. 2013. CD4 CTL: living up to the challenge. *Seminars in Immunology* **25**:273-281.
DOI: https://doi.org/10.1016/j.smim.2013.10.022, PMID: 24246226

Ciucci T, Vacchio MS, Bosselut R. 2017. A STAT3-dependent transcriptional circuitry inhibits cytotoxic gene expression in T cells. *PNAS* **114**:13236-13241. DOI: https://doi.org/10.1073/pnas.1711160114, PMID: 29180433

Cohen CJ, Crome SQ, MacDonald KG, Dai EL, Mager DL, Levings MK. 2011. Human Th1 and Th17 cells exhibit epigenetic stability at signature cytokine and transcription factor Loci. *The Journal of Immunology* **187**:5615-5626. DOI: https://doi.org/10.4049/jimmunol.1101058

Cosmi L, Annunziato F, Galli MIG, Maggi RME, Nagata K, Romagnani S. 2000. CRTH2 is the most reliable marker for the detection of circulating human type 2 Th and type 2 T cytotoxic cells in health and disease. *European Journal of Immunology* **30**:2972-2979. DOI: https://doi.org/10.1002/1521-4141(200010)30:10<2972::AID-IMMU2972>3.0.CO;2-, PMID: 11069080

Couturier J, Hutchison AT, Medina MA, Gingaras C, Urvil P, Yu X, Nguyen C, Mahale P, Lin L, Kozinetz CA, Schmitz JE, Kimata JT, Savidge TC, Lewis DE. 2014. HIV replication in conjunction with granzyme B production by CCR5+ memory CD4 T cells: Implications for bystander cell and tissue pathologies. *Virology* **462-463**:175-188. DOI: https://doi.org/10.1016/j.virol.2014.06.008, PMID: 24999042

Cruz-Guilloty F, Pipkin ME, Djuretic IM, Levanon D, Lotem J, Lichtenheld MG, Groner Y, Rao A. 2009. Runx3 and T-box proteins cooperate to establish the transcriptional program of effector CTLs. *The Journal of Experimental Medicine* **206**:51-59. DOI: https://doi.org/10.1084/jem.20081242, PMID: 19139168

Curran MA, Geiger TL, Montalvo W, Kim M, Reiner SL, Al-Shamkhani A, Sun JC, Allison JP. 2013. Systemic 4-1BB activation induces a novel T cell phenotype driven by high expression of Eomesodermin. *The Journal of Experimental Medicine* **210**:743-755. DOI: https://doi.org/10.1084/jem.20121190, PMID: 23547098

Dedeurwaerder S, Defrance M, Bizet M, Calonne E, Bontempi G, Fuks F. 2014. A comprehensive overview of Infinium HumanMethylation450 data processing. *Briefings in Bioinformatics* **15**:929-941. DOI: https://doi.org/10.1093/bib/bbt054, PMID: 23990268

Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F. 2011. Evaluation of the Infinium Methylation 450K technology. *Epigenomics* **3**:771-784. DOI: https://doi.org/10.2217/epi.11.105, PMID: 22126295

DEEP Consortium, Durek P, Nordström K, Gasparoni G, Salhab A, Kressler C, de Almeida M, Bassler K, Ulas T, Schmidt F, Xiong J, Glažar P, Klironomos F, Sinha A, Kinkley S, Yang X, Arrigoni L, Amirabad AD, Ardakani FB, Feuerbach L, Gorka O, et al. 2016. Epigenomic profiling of human CD4⁺T cells supports a linear differentiation model and highlights molecular regulators of memory development. *Immunity* **45**:1148-1161. DOI: https://doi.org/10.1016/j.immuni.2016.10.022, PMID: 27851915

Dimova T, Brouwer M, Gosselin F, Tassignon J, Leo O, Donner C, Marchant A, Vermijlen D. 2015. Effector Vγ9Vδ2 T cells dominate the human fetal γδ T-cell repertoire. *PNAS* **112**:E556-E565. DOI: https://doi.org/10.1073/pnas.1412058112, PMID: 25617367

Djuretic IM, Levanon D, Negreanu V, Groner Y, Rao A, Ansel KM. 2007. Transcription factors T-bet and Runx3 cooperate to activate Ifng and silence Il4 in T helper type 1 cells. *Nature Immunology* **8**:145-153. DOI: https://doi.org/10.1038/ni1424, PMID: 17195845

Dominguez CX, Amezquita RA, Guan T, Marshall HD, Joshi NS, Kleinstein SH, Kaech SM. 2015. The transcription factors ZEB2 and T-bet cooperate to program cytotoxic T cell terminal differentiation in response to LCMV viral infection. *The Journal of Experimental Medicine* **212**:2041-2056. DOI: https://doi.org/10.1084/jem.20150186, PMID: 26503446

Du P, Kibbe WA, Lin SM. 2008. lumi: a pipeline for processing Illumina microarray. *Bioinformatics* **24**:1547-1548. DOI: https://doi.org/10.1093/bioinformatics/btn224, PMID: 18467348

Dumitriu IE. 2015. The life (and death) of CD4+ CD28(null) T cells in inflammatory diseases. *Immunology* **146**:185-193. DOI: https://doi.org/10.1111/imm.12506, PMID: 26190355

Eshima K, Chiba S, Suzuki H, Kokubo K, Kobayashi H, Iizuka M, Iwabuchi K, Shinohara N. 2012. Ectopic expression of a T-box transcription factor, eomesodermin, renders CD4(+) Th cells cytotoxic by activating both perforin- and FasL-pathways. *Immunology Letters* **144**:7-15. DOI: https://doi.org/10.1016/j.imlet.2012.02.013, PMID: 22425747

Fu J, Zhang Z, Zhou L, Qi Z, Xing S, Lv J, Shi J, Fu B, Liu Z, Zhang JY, Jin L, Zhao Y, Lau GK, Zhao J, Wang FS. 2013. Impairment of CD4+ cytotoxic T cells predicts poor survival and high recurrence rates in patients with hepatocellular carcinoma. *Hepatology* **58**:139-149. DOI: https://doi.org/10.1002/hep.26054, PMID: 22961630

He L, Hakimi J, Salha D, Miron I, Dunn P, Radvanyi L. 2005. A sensitive flow cytometry-based cytotoxic T-lymphocyte assay through detection of cleaved caspase 3 in target cells. *Journal of Immunological Methods* **304**:43-59. DOI: https://doi.org/10.1016/j.jim.2005.06.005, PMID: 16076473

Intlekofer AM, Banerjee A, Takemoto N, Gordon SM, Dejong CS, Shin H, Hunter CA, Wherry EJ, Lindsten T, Reiner SL. 2008. Anomalous type 17 response to viral infection by CD8+ T cells lacking T-bet and eomesodermin. *Science* **321**:408-411. DOI: https://doi.org/10.1126/science.1159806, PMID: 18635804

Jha SS, Chakraborty NG, Singh P, Mukherji B, Dorsky DI. 2015. Knockdown of T-bet expression in Mart-127-35 - specific T-cell-receptor-engineered human CD4(+) CD25(+) and CD8(+) T cells attenuates effector function. *Immunology* **145**:124-135. DOI: https://doi.org/10.1111/imm.12431, PMID: 25495780

Johnson S, Eller M, Teigler JE, Maloveste SM, Schultz BT, Soghoian DZ, Lu R, Oster AF, Chenine AL, Alter G, Dittmer U, Marovich M, Robb ML, Michael NL, Bolton D, Streeck H. 2015. Cooperativity of HIV-specific cytolytic CD4 T cells and CD8 T cells in control of HIV Viremia. *Journal of Virology* **89**:7494-7505. DOI: https://doi.org/10.1128/JVI.00438-15, PMID: 25972560
Kaplan MJ, Lu Q, Wu A, Attwood J, Richardson B. 2004. Demethylation of promoter regulatory elements contributes to perforin overexpression in CD4+ lupus T cells. *Journal of Immunology* **172**:3652-3661.
DOI: https://doi.org/10.4049/jimmunol.172.6.3652

Kataoka T, Shinohara N, Takayama H, Takaku K, Kondo S, Yonehara S, Nagai K. 1996. Concanamycin A, a powerful tool for characterization and estimation of contribution of perforin- and Fas-based lytic pathways in cell-mediated cytotoxicity. *Journal of Immunology* **156**:3678-3686.

Kim B, Sasaki Y, Egawa T. 2015. Restriction of nonpermissive RUNX3 protein expression in T lymphocytes by the Kozak Sequence. *The Journal of Immunology* **195**:1517-1523. DOI: https://doi.org/10.4049/jimmunol.1501039

Liu X, Chen X, Zhong B, Wang A, Wang X, Chu F, Nurieva RI, Yan X, Chen P, van der Flier LG, Nakatsukasa H, Neelapu SS, Chen W, Clevers H, Tian Q, Qi H, Wei L, Dong C. 2014. Transcription factor achaete-scute homologue 2 initiates follicular T-helper-cell development. *Nature* **507**:513-518. DOI: https://doi.org/10.1038/nature12910, PMID: 24463518

Liu X, Taylor BJ, Sun G, Bosselut R. 2005. Analyzing expression of perforin, Runx3, and Thpok genes during positive selection reveals Activation of CD8-differentiation programs by MHC II-signaled thymocytes. *The Journal of Immunology* **175**:4465-4474. DOI: https://doi.org/10.4049/jimmunol.175.7.4465

Ma Y, Yuan B, Zhuang R, Zhang Y, Liu B, Zhang C, Zhang Y, Yu H, Yi J, Yang A, Jin B. 2015. Hantaan virus infection induces both Th1 and ThGranzyme B+ cell immune responses that associated with viral control and clinical outcome in humans. *PLoS pathogens* **11**:e1004788. DOI: https://doi.org/10.1371/journal.ppat.1004788, PMID: 25836633

Mackay LK, Minnich M, Kragten NA, Liao Y, Nota B, Seillet C, Zaid A, Man K, Preston S, Freestone D, Braun A, Wynne-Jones E, Behr FM, Stark R, Pellicci DG, Godfrey DI, Belz GT, Pellegrini M, Gebhardt T, Busslinger M, et al. 2016. Hobit and Blimp1 instruct a universal transcriptional program of tissue residency in lymphocytes. *Science* **352**:459-463. DOI: https://doi.org/10.1126/science.aad2035, PMID: 27102484

Marshall NB, Vong AM, Devarajan P, Brauner MD, Kuang Y, Nayar R, Schutten EA, Castonguay CH, Berg LJ, Nutt SL, Swain SL. 2017. NKG2C/E marks the unique cytotoxic CD4 T Cell Subset, ThCTL, generated by influenza infection. *The Journal of Immunology* **198**:1142-1155. DOI: https://doi.org/10.4049/jimmunol.1601297

Mittal P, Abblett R, Ryan JM, Hagymasi AT, Agyekum-Yamoah A, Svedova J, Reiner SL, St. Rose M-C, Hanley MP, Vella AT, Adler AJ. 2018. An immunotherapeutic CD137 agonist releases eomesodermin from ThPOK repression in CD4 T Cells. *The Journal of Immunology* **200**:1513-1526. DOI: https://doi.org/10.4049/jimmunol.1701039

Mucida D, Husain MM, Muroi S, van Wijk F, Shinnakasu R, Naoe Y, Reis BS, Huang Y, Lambolez F, Docherty M, Attinger A, Shui JW, Kim G, Lena CJ, Sakaguchi S, Miyamoto C, Wang P, Atarashi K, Park Y, Nakayama T, et al. 2013. Transcriptional reprogramming of mature CD4⁺ helper T cells generates distinct MHC class II-restricted cytotoxic T lymphocytes. *Nature Immunology* **14**:281-289. DOI: https://doi.org/10.1038/ni.2523, PMID: 2333478

Narasimhan S, Falkenberg VR, Khin MM, Rajeevan MS. 2009. Determination of quantitative and site-specific DNA methylation of perforin by pyrosequencing. *BMC Research Notes* **2**:104. DOI: https://doi.org/10.1186/1756-0500-2-104, PMID: 19523225

Nazarov VI, Pogorelyy MV, Komech EA, Zvyagin IV, Bolotin DA, Shugay M, Chudakov DM, Lebedev YB, Mamedov IZ. 2015. tcR: an R package for T cell receptor repertoire advanced data analysis. *BMC Bioinformatics* **16**:175. DOI: https://doi.org/10.1186/s12859-015-0613-1, PMID: 26017500

Omilusik KD, Best JA, Yu B, Goossens S, Weidemann A, Nguyen JV, Seuntjens E, Stryjewska A, Zweier C, Roychoudhuri R, Gattinoni L, Bird LM, Higashi Y, Kondoh H, Huylebroeck D, Haigh J, Goldrath AW. 2015. Transcriptional repressor ZEB2 promotes terminal differentiation of CD8+ effector and memory T cell populations during infection. *The Journal of Experimental Medicine* **212**:2027-2039. DOI: https://doi.org/10.1084/jem.20150194, PMID: 26503445

Park K, He X, Lee HO, Hua X, Li Y, Wiest D, Kappes DJ. 2010. TCR-mediated ThPOK induction promotes development of mature (CD24-) gammadelta thymocytes. *The EMBO Journal* **29**:2329-2341. DOI: https://doi.org/10.1038/emboj.2010.113, PMID: 20551904

Parronchi P, De Carli M, Manetti R, Simonelli C, Sampognaro S, Piccinni MP, Macchia D, Maggi E, Del Prete G, Romagnani S. 1992. IL-4 and IFN (alpha and gamma) exert opposite regulatory effects on the development of cytolytic potential by Th1 or Th2 human T cell clones. *Journal of Immunology* **149**:2977-2983.

Pearce EL, Mullen AC, Martins GA, Krawczyk CM, Hutchins AS, Zediak VP, Banica M, DiCioccio CB, Gross DA, Mao CA, Shen H, Cereb N, Yang SY, Lindsten T, Rossant J, Hunter CA, Reiner SL. 2003. Control of effector CD8+ T cell function by the transcription factor Eomesodermin. *Science* **302**:1041-1043. DOI: https://doi.org/10.1126/science.1090148, PMID: 14605368

Pipkin ME, Rao A, Lichtenheld MG. 2010. The transcriptional control of the perforin locus. *Immunological Reviews* **235**:55-72. DOI: https://doi.org/10.1111/j.0105-2896.2010.00905.x, PMID: 20536555

Quezada SA, Simpson TR, Peggs KS, Merghoub T, Vider J, Fan X, Blasberg R, Yagita H, Muranski P, Antony PA, Restifo NP, Allison JP. 2010. Tumor-reactive CD4(+) T cells develop cytotoxic activity and eradicate large established melanoma after transfer into lymphopenic hosts. *The Journal of Experimental Medicine* **207**:637-650. DOI: https://doi.org/10.1084/jem.20091918, PMID: 20156971

Raveney BJ, Oki S, Hohjoh H, Nakamura M, Sato W, Murata M, Yamamura T. 2015. Eomesodermin-expressing T-helper cells are essential for chronic neuroinflammation. *Nature Communications* **6**:8437. DOI: https://doi.org/10.1038/ncomms9437, PMID: 26436530
Reis BS, Rogoz A, Costa-Pinto FA, Taniuchi I, Mucida D. 2013. Mutual expression of the transcription factors Runx3 and ThPOK regulates intestinal CD4⁺ T cell immunity. *Nature Immunology* **14**:271–280. DOI: https://doi.org/10.1038/ni.2518, PMID: 23334789

Riaz T, Sollid LM, Olsen I, de Souza GA. 2016. Quantitative Proteomics of Gut-Derived Th1 and Th1/Th17 Clones Reveal the Presence of CD28+ NKG2D- Th1 Cytotoxic CD4+ T cells. *Molecular & Cellular Proteomics* **15**:1007–1016. DOI: https://doi.org/10.1074/mcp.M115.050138, PMID: 26637539

Rivino L, Messi M, Jarrossay D, Lanzavecchia A, Sallusto F, Geginat J. 2004. Chemokine receptor expression identifies Pre-T helper (Th)1, Pre-Th2, and nonpolarized cells among human CD4+ central memory T cells. *The Journal of Experimental Medicine* **200**:725–735. DOI: https://doi.org/10.1084/jem.20040774, PMID: 15381728

Sallusto F, Lenig D, Mackay CR, Lanzavecchia A. 1998. Flexible programs of chemokine receptor expression on human polarized T helper 1 and 2 lymphocytes. *The Journal of Experimental Medicine* **187**:875–883. DOI: https://doi.org/10.1084/jem.187.6.875, PMID: 9500790

Sellars M, Huh JR, Day K, Issuree PD, Galan C, Gobeil S, Absher D, Green MR, Littman DR. 2015. Regulation of DNA methylation dictates Cd4 expression during the development of helper and cytotoxic T cell lineages. *Nature Immunology* **16**:746–754. DOI: https://doi.org/10.1038/ni.3198, PMID: 26030024

Setoguchi R, Taniuchi I, Bevan MJ. 2009. ThPOK derepression is required for robust CD8 T cell responses to viral infection. *The Journal of Immunology* **183**:4467–4474. DOI: https://doi.org/10.4049/jimmunol.0901428

Sheehy ME, McDermott AB, Furlan SN, Klenerman P, Nixon DF. 2001. A novel technique for the fluorometric assessment of T lymphocyte antigen specific lysis. *Journal of Immunological Methods* **249**:99–110. DOI: https://doi.org/10.1016/S0022-1759(00)00329-X, PMID: 11226468

Shlyueva D, Stampfel G, Stark A. 2014. Transcriptional enhancers: from properties to genome-wide predictions. *Nature Reviews Genetics* **15**:272–286. DOI: https://doi.org/10.1038/nrg3682, PMID: 24614317

Spinelli L, Carpentier S, Montañana Sanchis F, Dalod M, Vu Manh TP. 2015. BubbleGUM: automatic extraction of phenotype molecular signatures and comprehensive visualization of multiple Gene Set Enrichment Analyses. *BMC Genomics* **16**:814. DOI: https://doi.org/10.1186/s12864-015-2012-4, PMID: 26481321

Ståhlberg A, Rusnakova V, Forootan A, Anderova M, Kubista M. 2013. RT-qPCR work-flow for single-cell data analysis. *Methods* **59**:80–88. DOI: https://doi.org/10.1016/j.ymeth.2012.09.007, PMID: 23021995

Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. 2005. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. *PNAS* **102**:15545–15550. DOI: https://doi.org/10.1073/pnas.0506580102, PMID: 16199517

Sujino T, London M, Hoytema van Konijnenburg DP, Rendon T, Buch T, Silva HM, Lafaille JJ, Reis BS, Mucida D. 2016. Tissue adaptation of regulatory and intraepithelial CD4⁺ T cells controls gut inflammation. *Science* **352**: 1581–1586. DOI: https://doi.org/10.1126/science.aaf3892, PMID: 27256884

Suzuki R, Shimodaira H. 2006. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. *Bioinformatics* **22**:1540–1542. DOI: https://doi.org/10.1093/bioinformatics/btl117, PMID: 16595560

Svec D, Andersson D, Pekny M, Sjöback R, Kubista M, Ståhlberg A. 2013. Direct cell lysis for single-cell gene expression profiling. *Frontiers in Oncology* **3**:274. DOI: https://doi.org/10.3389/fonc.2013.00274, PMID: 24224157

Swain SL, McKinstry KK, Strutt TM. 2012. Expanding roles for CD4⁺ T cells in immunity to viruses. *Nature reviews. Immunology* **12**:136–148. DOI: https://doi.org/10.1038/nri3152, PMID: 22266691

Takeuchi A, Badr MS, Miyauchi K, Ishihara C, Onishi R, Guo Z, Sasaki Y, Ike H, Takumi A, Tsuji NM, Murakami Y, Katakai T, Kubo M, Saito T. 2016. CRTAM determines the CD4+ cytotoxic T lymphocyte lineage. *The Journal of Experimental Medicine* **213**:123–138. DOI: https://doi.org/10.1084/jem.20150519, PMID: 26694968

Vacchio MS, Bosselut R. 2016. What happens in the thymus does not stay in the thymus: how T cells recycle the CD4+–CD8+ lineage commitment transcriptional circuitry to control their function. *Journal of Immunology* **196**: 4848–4856. DOI: https://doi.org/10.4049/jimmunol.1600415, PMID: 27260768

Van Caeneghem Y, De Munter S, Tieppo P, Goetgeluk G, Weening K, Verstichel G, Bonte S, Taghon T, Leclercq G, Kerre T, Debets R, Vermijlen D, Abken H, Vandekerckhove B. 2017. Antigen receptor-redirected T cells derived from hematopoietic precursor cells lack expression of the endogenous TCR/CD3 receptor and exhibit specific antitumor capacities. *Oncoimmunology* **6**:e1283460. DOI: https://doi.org/10.1080/2162402X.2017.1283460, PMID: 28405508

van de Berg PJ, van Leeuwen EM, ten Berge IJ, van Lier R. 2008. Cytotoxic human CD4(+) T cells. *Current Opinion in Immunology* **20**:339–343. DOI: https://doi.org/10.1016/j.coi.2008.03.007, PMID: 18440213

van Gisbergen KP, Kragten NA, Hertoghs KM, Wensveen FM, Jonjic S, Hamann J, Nolte MA, van Lier RA. 2012. Mouse Hobit is a homolog of the transcriptional repressor Blimp-1 that regulates NKT cell effector differentiation. *Nature Immunology* **13**:864–871. DOI: https://doi.org/10.1038/ni.2393, PMID: 22885984

van Leeuwen EMM, Remmerswaal EBM, Vossen MTM, Rowshani AT, Wertheim-van Dillen PME, van Lier RAW, ten Berge IJM. 2004. Emergence of a CD4+CD28- Granzyme B+, Cytomegalovirus-Specific T Cell Subset after Recovery of Primary Cytomegalovirus Infection. *The Journal of Immunology* **173**:1834–1841. DOI: https://doi.org/10.4049/jimmunol.173.3.1834

Verma S, Weiskopf D, Gupta A, McDonald B, Peters B, Sette A, Benedict CA. 2016. Cytomegalovirus-Specific CD4 T cells are cytolytic and mediate vaccine protection. *Journal of Virology* **90**:650–658. DOI: https://doi.org/10.1128/JVI.02123-15

Vieira Braga FA, Hertoghs KM, Kragten NA, Doody GM, Barnes NA, Remmerswaal EB, Hsiao CC, Moerland PD, Wouters D, Derks IA, van Stijn A, Demkes M, Hamann J, Eldering E, Nolte MA, Tooze RM, ten Berge IJ, van
Gisbergen KP, van Lier RA. 2015a. Blimp-1 homolog Hobit identifies effector-type lymphocytes in humans. *European Journal of Immunology* **45**:2945-2958. DOI: https://doi.org/10.1002/eji.201545650, PMID: 26179882

Vieira Braga FA, Hertoghs KM, van Lier RA, van Gisbergen KP. 2015b. Molecular characterization of HCMV-specific immune responses: Parallels between CD8(+) T cells, CD4(+) T cells, and NK cells. *European Journal of Immunology* **45**:2433-2445. DOI: https://doi.org/10.1002/eji.201545495, PMID: 26228786

Vinton CL, Ortiz AM, Calantone N, Mudd JC, Deleage C, Morcock DR, Whitted S, Estes JD, Hirsch VM, Brenchley JM. 2017. Cytotoxic T cell functions accumulate When CD4⁺ T cells in african green monkeys. *The Journal of Immunology* **198**:4403-4412. DOI: https://doi.org/10.4049/jimmunol.1700136

Wang L, Carr T, Xiong Y, Wildt KF, Zhu J, Feigenbaum L, Bendelac A, Bosselut R. 2010. The sequential activity of Gata3 and Thpok is required for the differentiation of CD1d-restricted CD4+ NKT cells. *European Journal of Immunology* **40**:2385-2390. DOI: https://doi.org/10.1002/eji.201040534, PMID: 20706986

Wang L, Wildt KF, Castro E, Xiong Y, Feigenbaum L, Tessarollo L, Bosselut R. 2008. The zinc finger transcription factor Zbtb7b represses CD8-lineage gene expression in peripheral CD4+ T cells. *Immunity* **29**:876-887. DOI: https://doi.org/10.1016/j.immuni.2008.09.019, PMID: 19062319

Wang Y, Godec J, Ben-Aissa K, Cui K, Zhao K, Pucsek AB, Lee YK, Weaver CT, Yagi R. Lazarevic V. 2014. The transcription factors T-bet and Runx are required for the ontogeny of pathogenic interferon-γ-producing T helper 17 cells. *Immunity* **40**:355-366. DOI: https://doi.org/10.1016/j.immuni.2014.01.002, PMID: 24530058

Weiskopf D, Bangs DJ, Sidney J, Kolla RV, De Silva AD, de Silva AM, Crotty S, Peters B, Sette A. 2015. Dengue virus infection elicits highly polarized CX3CR1+ cytotoxic CD4+ T cells associated with protective immunity. *PNAS* **112**:E4256-4263. DOI: https://doi.org/10.1073/pnas.1505956112, PMID: 26195744

Wilkinson TM, Li CKF, Chui CSC, Huang AKY, Perkins M, Liebner JC, Lambkin-Williams R, Gilbert A, Oxford J, Nicholas B, Staples KJ, Dong T, Douek DC, McMichael AJ, Xu X-N. 2012. Preexisting influenza-specific CD4+ T cells correlate with disease protection against influenza challenge in humans. *Nature Medicine* **18**:274-280. DOI: https://doi.org/10.1038/nm.2612

Wilson CB, Rowell E, Sekimata M. 2009. Epigenetic control of T-helper-cell differentiation. *Nature reviews. Immunology* **9**:91-105. DOI: https://doi.org/10.1038/nri2487, PMID: 19151746

Xie Y, Akpinarli A, Maris C, Hipkiss EL, Lane M, Kwon EK, Muranski P, Restifo NP, Antony PA. 2010. Naive tumor-specific CD4(+) T cells differentiated in vivo eradicate established melanoma. *The Journal of Experimental Medicine* **207**:651-667. DOI: https://doi.org/10.1084/jem.20091921, PMID: 20156973

Xiong Y, Castro E, Yagi R, Zhu J, Lesourne R, Love PE, Feigenbaum L, Bosselut R. 2013. Thpok-independent repression of Runx3 by Gata3 during CD4+ T-cell differentiation in the thymus. *European Journal of Immunology* **43**:918-928. DOI: https://doi.org/10.1002/eji.201242944, PMID: 23310955