Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection

Resource

Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection Highlights

Authors +

d

Examined chromatin accessibility in endogenous CD8 T cells during viral infection

d

Many changes from naive T cells are stably maintained in differentiated cell types

d

Groups of differentially accessible regions are shared between different subsets

d

Exhaustion-specific accessible regions are enriched for NFAT and Nr4a binding sites

James P. Scott-Browne, Isaac F. Lo´pez-Moyado, Sara Trifari, Victor Wong, Lukas Chavez, Anjana Rao, Renata M. Pereira

Correspondence [email protected]

In Brief Scott-Browne et al. assess genome-wide changes in chromatin accessibility and gene expression in endogenous CD8+ T cells responding to acute and chronic viral infection. They describe shared patterns of differentially accessible regions between naive, effector, memory, and exhausted subsets that are associated with binding sites for known and novel transcription factors.

Accession Numbers GSE88987

Scott-Browne et al., 2016, Immunity 45, 1–14 December 20, 2016 ª 2016 Elsevier Inc. http://dx.doi.org/10.1016/j.immuni.2016.10.028

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

Immunity

Resource Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection James P. Scott-Browne,1 Isaac F. Lo´pez-Moyado,1 Sara Trifari,1 Victor Wong,1 Lukas Chavez,1,4 Anjana Rao,1,2,3,6,* and Renata M. Pereira1,5 1Division

of Signaling and Gene Expression, La Jolla Institute, 9420 Athena Circle, San Diego, CA 92037, USA Consortium for Regenerative Medicine, La Jolla, CA 92037, USA 3Department of Pharmacology and Moores Cancer Center, University of California San Diego, La Jolla, CA 92037, USA 4Computational Oncoepigenomics Group, Division of Pediatric Neurooncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, 69120 Heidelberg, Germany 5Present address: Instituto de Microbiologia Prof. Paulo de Goes, Universidade Federal do Rio de Janeiro, 21941-900 Rio de Janeiro, RJ, Brasil 6Lead Contact *Correspondence: [email protected] http://dx.doi.org/10.1016/j.immuni.2016.10.028 2Sanford

SUMMARY

In response to acute infection, naive CD8+ T cells expand, differentiate into effector cells, and then contract to a long-lived pool of memory cells after pathogen clearance. During chronic infections or in tumors, CD8+ T cells acquire an ‘‘exhausted’’ phenotype. Here we present genome-wide comparisons of chromatin accessibility and gene expression from endogenous CD8+ T cells responding to acute and chronic viral infection using ATAC-seq and RNAseq techniques. Acquisition of effector, memory, or exhausted phenotypes was associated with stable changes in chromatin accessibility away from the naive T cell state. Regions differentially accessible between functional subsets in vivo were enriched for binding sites of transcription factors known to regulate these subsets, including E2A, BATF, IRF4, T-bet, and TCF1. Exhaustion-specific accessible regions were enriched for consensus binding sites for NFAT and Nr4a family members, indicating that chronic stimulation confers a unique accessibility profile on exhausted cells.

INTRODUCTION During acute viral infection, naive antigen-specific CD8+ T cells expand and differentiate to yield effector T cells that enable resolution of the infection (Williams and Bevan, 2007). After viral clearance, most expanded cells die, but a small proportion survive as a long-lived memory population that rapidly produces cytokines and reacquires cytotoxic activity upon secondary exposure to antigen, thus providing protective immunity. These phases of the immune response are typical of acute infections where the virus is eliminated, such as with the Armstrong (Arm5) strain of lymphocytic choriomeningitis virus (LCMV) (Wherry and Ahmed, 2004). When virus persists during chronic or latent infection—as in mice or humans infected with LCMV clone 13, human immunode-

ficiency virus, hepatitis B virus, or hepatitis C virus—CD8+ T cells often enter a state of unresponsiveness to further stimulation that has been termed ‘‘exhaustion’’ (Moskophidis et al., 1993; Wherry, 2011). Exhausted CD8+ T cells are hyporesponsive to stimulation with decreased cytokine production, reduced ability to lyse target cells, and increased expression of several inhibitory cell surface receptors including PD-1 (programmed death 1), LAG3 (lymphocyte-activation gene 3), TIM3 (T cell immunoglobulin and mucin domain-containing protein 3), and CTLA-4 (cytotoxic T lymphocyte-associated protein 4) (Barber et al., 2006; Blackburn et al., 2009; Wherry et al., 2003; Yamamoto et al., 2011). The exhausted phenotype is also common in tumor-infiltrating CD8+ cells (Crespo et al., 2013), and antibody therapy targeting inhibitory receptors (‘‘checkpoint blockade’’) has been remarkably effective in cancer immunotherapy (Pardoll, 2012). Although several transcription factors (TFs), including BATF, IRF4, and the T-box family members T-bet and eomesodermin (Eomes), are known to participate in the formation and function of effector and memory CD8+ T cells (Dominguez et al., 2015; Intlekofer et al., 2005; Kurachi et al., 2014), the molecular mechanisms that determine the exhausted phenotype are not well understood (Speiser et al., 2014; Wherry and Kurachi, 2015). Persistent antigen stimulation appears to be a dominant factor in inducing exhaustion in tumor-infiltrating T cells (Schietinger et al., 2016); consistent with these findings, we previously implicated a gene expression program driven by the antigen-activated TF NFAT in CD8+ T cell exhaustion (Martinez et al., 2015). However, given the limiting cell numbers available, it was technically difficult to confirm these findings by performing ChIP-seq (chromatin immunoprecipitation followed by sequencing) for NFAT1 and other TFs in exhausted cells. Changes in transcriptional programs are controlled not only through the action of TFs near transcription start sites (TSSs), but also through epigenetic changes in a variety of DNA and histone modifications at regulatory elements throughout the genome (Kouzarides, 2007). Active regulatory elements that bind TFs can be defined operationally by their ‘‘accessibility’’ to enzymes such as DNase I and micrococcal nuclease (ENCODE Project Consortium, 2012; Vierstra et al., 2014). Although DNase I hypersensitivity assays are relatively cumbersome and require large numbers of cells (Pipkin et al., 2010), but accessible regions Immunity 45, 1–14, December 20, 2016 ª 2016 Elsevier Inc. 1

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

of chromatin can be identified reliably even in small cell numbers by ATAC-seq (assay for transposase-accessible chromatin using sequencing), a technique that measures accessibility to an engineered Tn5 transposon containing flanking sequencing adapters (Buenrostro et al., 2013). Here we have mapped accessible regulatory elements by ATAC-seq in naive, effector, memory, and exhausted CD8+ T cells from mice with acute or chronic LCMV infection. We identified dynamic changes in chromatin accessibility in CD8+ T cells, with clusters of regions with shared accessibility profiles between different subsets. By motif enrichment analysis of the regulatory elements and comparison to transcriptional profiles obtained by RNA-sequencing (RNA-seq), we have confirmed the involvement of NFAT in the exhaustion program and implicated new NFAT-induced TFs in CD8+ T cell exhaustion. Our data constitute a comprehensive analysis of the similarities and differences among functionally distinct CD8+ T cell subsets and provide a valuable resource for future comparisons of these cell types in different tissues and disease models or after genetic manipulation or antibody blockade. RESULTS Identification of Accessible Regions in CD8+ T Cells We used ATAC-seq to assess the accessible regions genomewide in CD8+ T cells responding to viral infection. Antigen-specific CD8+ T cells, defined as CD44hi and staining with H-2Db tetramers containing a peptide from LCMV glycoprotein amino acids 33–41 (H-2Db-gp33-41), were collected by fluorescenceactivated cell sorting (FACS) from mice after LMCV Arm5 or clone 13 infection (Figures 1A, S1A, and S1C–S1E). Total effector cells and KLRG1+IL7Ra– (short-lived effector cells [SLECs]) or KLRG1–IL7Ra+ (memory precursor effector cells [MPECs]) subsets (Figures S1A–S1C) were collected 8 days after LCMV Arm5 infection. The frequencies of SLECs and MPECs were around 50% and 8% of total H-2Db-gp33-41 tetramer-positive cells, respectively (Figures S1B and S1C). At 35 days after infection with LMCV Arm5, we sorted cells expressing IL7Ra (memory) or KLRG1 (d35 KLRG1+) (Figure S1D). Cells with an exhausted phenotype (PD-1hiTim3hi) were isolated from mice infected with LCMV clone 13 for 20 days (Figure S1E). We also compared ATAC-seq signals in naive splenic CD8+ T cells and in cultured CD8+ T cells after in vitro stimulation under various conditions or after transduction with retroviral constructs (Figure 1A). As in previous studies (Buenrostro et al., 2013), ATAC-seq coverage in these samples was concentrated with high signal in narrow regions throughout the genome (illustrated for Ifng in Figure 1B). We called peaks on individual replicates (Figure S1F), merged these to 91,099 regions, counted the number of Tn5 insertions from each replicate in all peaks, and normalized these data using signals at regions with low variance with DESeq2 (Supplemental Experiment Procedures; Love et al., 2014). Regions with higher signal had lower inter-replicate error (Figures S1G and S1H), so we restricted our analysis to 70,532 regions that met a minimum threshold of mean ATAC-seq signal in at least one cell type (Table S1 and Supplemental Experimental Procedures). Comparing the signal at all regions between all samples, each biological replicate clustered with the other replicates from the same group, and two technical replicates were as 2 Immunity 45, 1–14, December 20, 2016

similar to each other as they were to independent biological replicates (Figure S1I). We also profiled transcriptional differences by RNA-seq of naive, effector, memory, and exhausted cells; each biological replicate clustered with replicates in the same cell type (Figure S1J and Table S1). In 70 kb of the Ifng locus, we identified 19 peaks that were accessible in at least one cell type (Figure 1B). The conserved region located at
22 kb relative to the transcription start site (TSS) was accessible in all samples, whereas other conserved regions were not accessible in naive T cells but were highly accessible in effector, memory, and exhausted T cells (Figure 1B). Cells that were pre-activated, cultured, and then briefly reactivated for short periods (2 hr) in vitro showed variable accessibility at these regions, as did cells transduced with a constitutively active NFAT1 that cannot cooperate with AP-1 (NFAT-CA-RIT) (Figure 1B; Martinez et al., 2015). The sequences around ATAC-seq peaks were highly conserved (Figure S1K), and the peaks were frequent in promoters, introns, and distal intergenic regions (Figure S1L). Although the Ifng locus illustrates changes in accessibility near the TSS and distal conserved elements, the signal at most TSSs was less variable compared to distal elements (Figure 1C). Considering all peaks, T cells that had responded to LCMV infection were most similar to each other, whereas naive T cells or cells activated and cultured in vitro were part of a separate cluster (Figures 1D and S1I). To identify individual regions that were differentially accessible between CD8+ T cell subsets, we computed differential coverage using DESeq2. Comparing CD8+ T cells activated and cultured in vitro to effector cells at day 8 after LMCV Arm5 infection, at least 9,000 regions were differentially accessible (Figure S1M). Of these, differentially accessible regions that were more accessible in cells cultured in vitro than in ex vivo isolated effector cells were also more accessible in naive CD8+ T cells; conversely, regions that were less accessible in cells cultured in vitro than in ex vivo isolated effector cells also tended to be less accessible in naive cells, compared to effector, memory, and exhausted CD8+ T cells (Figure S1N). Thus, the accessibility profile of CD8+ T cells activated and cultured in vitro under these conditions remained more similar to that of naive CD8+ T cells and did not completely recapitulate the accessibility profile of effector CD8+ T cells that develop over several days of viral infection in vivo. We previously found that high (100 U/mL) or low (10 U/mL) concentrations of IL2 can bias CD8+ T cells cultured in vitro toward effector or memory phenotypes, respectively (Pipkin et al., 2010). However, very few regions were differentially accessible between high and low IL-2 conditions (Figure S1O), indicating that exposure to IL-2 under these conditions did not alter chromatin accessibility substantially, despite clear differences in transcriptional profiles. Stable and Dynamic Changes in Accessible Chromatin of CD8+ T Cells in Acute Viral Infection We assessed chromatin accessibility changes during in vivo responses to acute viral infection by comparing naive, effector, and memory CD8+ T cells. Of the 45,439 regions that were accessible in any of these subsets, more than 12,000 were differentially accessible when comparing naive and effector T cells (Figure 2A); fewer (5,100) were differentially accessible between naive and memory (Figure 2B); and effector and memory

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

A

B

C

D

Figure 1. High ATAC-Seq Signal in CD8+ T Cells at Conserved Regions in Promoters and Distal Regulatory Elements (A) CD8+ T cell populations collected for ATAC-seq comparison. (B) Mean ATAC-seq coverage at the 70 kb Ifng locus with a scale of 0–1,200 for all tracks. (C) k-means clustered heatmap of mean normalized counts or log2 fold-change from global mean at all peaks. (D) Pairwise euclidian distance comparison of asinh-transformed ATAC-seq signal per peak for all populations using all peaks accessible in at least one cell type. Data in (B)–(D) are from mean of at least two independent samples, except for a single d35 KLRG1+ replicate. See also Figure S1.

T cells were the most similar with only 3,000 differentially accessible regions (Figure 2C). Most of these differentially accessible regions were located distal to the TSS, whereas commonly accessible regions were more evenly distributed

between TSS proximal and distal elements (Figures S2A–S2C). Similarly, comparison of RNA-seq data identified more than 1,500 differentially expressed genes between naive and effector cells (Figure S2D) but only 600–800 differentially expressed Immunity 45, 1–14, December 20, 2016 3

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

A

D

G

H

B

E I

J C

F

K

L

M

Figure 2. Dynamic Changes in Chromatin Accessibility Occur in Antigen-Specific Effector and Memory CD8+ T Cells Responding to Acute Viral Infection (A–C) Scatterplots of mean ATAC-seq counts per peak comparing the indicated samples. (D–F) Boxplots of ATAC-seq counts per peak from the indicated samples (labeled at bottom) at common or differentially accessible regions from the comparison labeled above. Box indicates interquartile range with whiskers ±1.5 times this range and outlier points. (G–K) Mean ATAC-seq coverage at Il7r (G), Ccr7 (H), Gzma (I), Gzmk (J), and Dmrta1 (K) loci with a scale of 0–1,200 (left) or RNA-seq gene expression for the indicated genes (right). (L and M) Venn diagrams illustrating intersection of differentially accessible regions from pairwise comparisons of naive, effector, and memory CD8+ T cells characterizing regions ‘‘specific’’ to a subset (L) or ‘‘not’’ in a subset (M) with p values and odds ratios from Fisher’s test comparisons. ATAC-seq data in (A)–(K) are from at least two independent replicates. RNA-seq data in (G)–(K) are mean of two independent replicates for RNA-seq. See also Figure S2.

4 Immunity 45, 1–14, December 20, 2016

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

genes between naive and memory (Figure S2E) or memory and effector (Figure S2F) cells. The differentially accessible regions were positively associated with changes in gene expression at genes located within 25 kb of the peaks (Figures S2G–S2I); for example, genes near peaks with higher ATAC-seq signal in effector cells than in naive had higher expression in effector cells compared to naive cells (Figure S2G). Regions that were less accessible in effector than in naive cells also tended to be less accessible in exhausted cells (Figure 2D). In memory cells, some of these regions had a high signal similar to naive (Il7r upstream region, Figure 2G) whereas others had a low signal similar to effector (Ccr7 downstream region, Figure 2H). Similar patterns were apparent at regions more accessible in effector than naive. For example, within the Gzma locus, regions such as the +4 kb and +19 kb peaks had higher signal in effector than naive cells, but the +19 kb peak was high in memory cells and the +4 kb peak was low in memory cells (Figure 2I). The signals at some activation-induced genes, including Fasl and Prf1, were highly similar between naive, effector, and memory CD8+ T cells despite differences in gene expression (Figure S2J), suggesting that changes in expression of these genes would most likely result from differential activation by TFs binding to these commonly accessible regions. We next focused on differentially accessible regions between memory cells and either naive or effector cells. There were more differentially accessible regions between naive and memory (Figure 2B) than between effector and memory (Figure 2C) cells. Regions with lower signal in memory than naive also tended to have lower signal in effector and exhausted cells, while peaks with higher signal in memory than naive tended to have higher signal in effector and exhausted cells (Figure 2E). Taken together, these comparisons identified regions, similar to the –12 kb and +2.5 kb peaks in the Gzmk locus and the +19 kb peak in the Gzma locus, that were stably altered in naive T cells after viral infection (Figures 2I and 2J). Despite the low number of differentially accessible regions between effector and memory cells (Figure 2C), many of the regions were associated with the greatest variation among naive, effector, memory, and exhausted T cell populations (Figure 2F). In an effort to identify accessible regions that were biased toward naive, effector, or memory CD8+ T cells, we assessed the intersection of the differentially accessible regions from each comparison (Figure 2L). Although a few regions were memory specific (Figure 2L), such as the Dmrta1 TSS (Figure 2K), substantial fractions of differentially accessible regions between memory and naive or effector cells were shared with other cell types (Figure 2M) and fewer than expected were absent in memory but shared by effector and naive cells (Figure 2M). These data indicated that at least 25% of the accessible regions changed at some point during CD8+ T cell activation and differentiation during acute infection and that 5%–10% of the CD8+ T cell accessible landscape underwent a stable change from a naive T cell after response to an acute viral infection. The effector population had some unique changes, with increased accessibility at some regions and reduced accessibility at others compared to naive and memory CD8+ T cells. Finally, most of the accessibility profile of memory CD8+ T cells was shared with either naive or effector T cells, with only a small fraction of total accessible regions unique to this population.

Chromatin Accessibility Profiles of SLECs and MPECs Are Similar to Those of Effector and Memory Cells We next compared cells with a terminal effector phenotype (SLECs) or a memory precursor phenotype (MPECs) at day 8 after LMCV Arm5 infection to assess the similarities of these subsets to effector and memory cells. Based on the signal at all accessible regions, SLECs were very similar to total effector cells and MPECs were less similar to both effector and memory cells (Figures 1D and S1I). In the direct comparison of SLECs and MPECs, we identified fewer than 1,000 differentially accessible regions (Figure 3A). Regions that were more accessible in SLECs than MPECs tended to have high accessibility in effector cells with progressively lower signal in exhausted, memory, and naive cells, supporting the overall similarity between effector cells and SLECs (Figure 3B). Regions that were more accessible in MPECs than SLECs tended to have high signal in naive and memory cells and lower signal in effector and exhausted cells (Figure 3B). Together these data suggested that although MPECs were very similar to effector cells, there was a partial bias toward the memory accessibility profile. This trend was also apparent when comparing the change in signal between SLECs and MPECs at differentially accessible regions between memory and effector cells. Peaks with higher signal in memory than effector cells also had higher signal in MPECs than SLECs (Figure 3C), while regions higher in effector than memory tended to have increased signal in SLECs compared to MPECs (Figure 3C), including the region near the Klrg1 TSS (Figure 3D). This bias was incomplete, as exemplified by the region 4–5 kb upstream of Aurkb, which included several peaks with comparable signal between MPECs, SLECs, and effector cells but lower in memory and naive. We also compared the accessibility between SLECs and effector cells as well as MPECs and memory cells (Figures S3A and S3B). Although SLECs were collected from female mice and effector cells were collected from male mice, there were fewer than 50 differentially accessible regions between these two subsets (Figure S3A). In addition, all but two of the regions more accessible in female SLECs than male effector cells were on the X chromosome (data not shown) and there were many female-specific peaks within the Xist/Tsix and Firre loci (Figure S3C). By comparison, there were more than 700 regions differentially accessible between female MPECs and male memory cells (Figure S3B). Those that were higher in MPECs tended to also have high signal in effector cells while those that were higher in memory tended to have high signal in memory and naive cells, indicating that these regions were largely associated with different CD8+ T cell function rather than reflecting the sex differences. Altogether, these data indicate that the accessibility assessed using the total effector cell population is highly similar to the SLEC subset, largely due to the relatively small fraction of MPECs (Figure S1B) and their similarity to SLECs (Figure 3A). MPECs were slightly biased toward the memory cell profile, but the chromatin profile of MPECs was more similar to that of SLECs (Figures 3A and S3B). Exhaustion-Specific Changes in Chromatin Accessibility We next compared effector and memory CD8+ T cells to CD8+ T cells with an exhausted phenotype, isolated from mice chronically infected with LCMV clone 13. There were approximately Immunity 45, 1–14, December 20, 2016 5

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

A

B

D

C

E

Figure 3. Memory Precursor Effector Cells Are Similar to Short-Lived Effector Cells with a Slight Bias toward Memory (A) Scatterplot of mean ATAC-seq counts per peak comparing the SLECs and MPECs. (B) Boxplot of ATAC-seq counts per peak from the indicated samples (labeled at bottom) at common or differentially accessible regions from the comparison labeled above. Box indicates interquartile range with whiskers ±1.5 times this range and outlier points. (C) Histograms of the log2 fold-change between effector and memory cells (top) or SLECs and MPECs (bottom) at regions differentially accessible between effector and memory. (D) Mean ATAC-seq coverage at Klrg1 and Aurkb loci with a scale of 0–1,200. (E) RNA-seq gene expression for Klrg1 and Aurkb. Data in (A)–(D) are from three independent replicates and in (E) is mean of two independent replicates. See also Figure S3.

5,000 and 3,000 differentially accessible regions in cells with an exhausted phenotype compared to memory and effector CD8+ T cells, respectively (Figures 4A and 4B). As observed in comparisons of naive, effector, and memory CD8+ T cells (Figures S2A– S2C), most of these differentially accessible regions were distal to the TSS (Figures S4A and S4B). Similar trends were observed for RNA-seq comparisons between exhausted cells and either effector or memory cells (Figures S4C and S4D), where more genes were differentially expressed between exhausted and memory cells than between exhausted and effector cells. As observed for comparisons of naive, effector, and memory cells, the changes in accessibility were more likely to be positively associated with changes in expression at nearby genes (Figures S4E and S4F). The accessibility profiles of effector and exhausted CD8+ T cells were similar, exemplified by the –7 kb and +19 kb peaks in the Havcr2 (encoding Tim3) locus (Figure 4E). Regions with 6 Immunity 45, 1–14, December 20, 2016

unique accessibility in exhausted cells were more clearly identifiable in comparison with effector cells (Figures 4B and 4D) than when compared with memory cells (Figures 4A and 4C), such as the –19 kb and –61 kb peaks near the Tox2 locus (Figure 4F) associated with high Tox2 gene expression (Figure 4F). Cells with an exhausted phenotype were also characterized by loss of accessibility compared to effector and memory CD8+ T cells (Figures 4C and 4D), such as intronic peaks in the Satb1 locus, associated with low Satb1 expression (Figure 4G). To identify features that were specific to, or absent from, effector, memory, and exhausted CD8+ T cells, we compared by intersections the differentially accessible regions among these groups. More regions were specific to memory or exhausted cells than to effector cells (Figure 4H). Relatively few regions were less accessible in effector than either memory or exhausted cells (Figure 4I), but there were 762 regions with low accessibility in exhausted cells compared to both effector and memory.

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

A

C

E

F

B

D G

H

I

Figure 4. Chronic Activation Profile Identified by Comparison of Viral Antigen-Specific Effector, Memory, and Exhausted CD8+ T Cells (A and B) Scatterplots of mean ATAC-seq counts per peak comparing the indicated samples. (C and D) Boxplots of ATAC-seq counts per peak from the indicated samples (labeled at bottom) at common or differentially accessible regions from the comparison labeled above. Box indicates interquartile range with whiskers ±1.5 times this range and outlier points. (E–G) Mean ATAC-seq coverage at Havcr2 (E), Tox2 (F), and Satb1 (G) loci with a scale of 0–1,200 (left) or RNA-seq gene expression for the indicated genes (right). (H and I) Venn diagrams illustrating intersection of differentially accessible regions from pairwise comparisons of effector, memory, and exhausted CD8+ T cells characterizing regions ‘‘specific’’ to a subset (H) or ‘‘not’’ in a subset (I) with p values and odds ratios from Fisher’s test comparisons. ATAC-seq data in (A)–(G) are from at least two independent replicates. RNA-seq data in (E)–(G) are mean of two independent replicates. See also Figure S4.

Despite exhaustion-specific changes, these data indicated that CD8+ T cells with an exhausted phenotype responding to a chronic infection shared a substantial fraction of chromatin accessibility with effector CD8+ T cells responding to acute infection. In particular, the accessibility at key effector-related genes including Ifng, Gzma, Gzmk, Fasl, and Prf1 (Figures 1B, 2I, 2J, and S2J) and inhibitory receptor loci including Havcr2/ Tim3, Lag3, and Ctla4 (Figures 4E and S4G) was similar between effector and exhausted cells. Differentially Accessible Regions Are Associated with Regulators of CD8+ T Cell Differentiation We identified 18,043 regions that were differentially accessible in at least one pairwise comparison of naive cells, effector cells, SLECs, MPECs, memory cells, and exhausted cells and compared the relative location of individual replicates on multidimensional scaling plots computed using ATAC-seq signal at these regions. SLECs, effector cells, and exhausted cells were

most distant from naive cells, whereas MPECs and memory cells were closer to naive (Figure 5A). The substantial overlap between groups of differentially accessible regions (Figures 2L, 2M, 4H, and 4I) suggested that there are functionally related groups of accessible regions that may explain the dissimilarity (Figure 5A). Using k-means clustering, we partitioned all 18,043 differentially accessible regions into 12 groups with shared profiles between naive cells, effector cells, SLECs, MPECs, memory cells, and exhausted CD8+ T cells (Figure 5B). In this comparison, relatively few effector- or memory-specific regions were apparent, with slight biases in clusters 11 and 8 in favor of effector and memory, respectively. In contrast, cluster 10 and clusters 3 and 6 identified groups that were largely biased toward naive or exhausted, respectively (Figure 5B). Regions in clusters 1 and 4 that were shared by naive and memory cells were enriched for HMG motifs and a high frequency of TCF1 ChIP-seq peaks from naive CD8+ T cells overlapped these regions (Figures 5C and 5D; Steinke et al., 2014). Immunity 45, 1–14, December 20, 2016 7

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

A

B

C

D

F

E

Figure 5. Differentially Accessible Regions in CD8+ T Cells Are Associated with bhLH, bZIP, HMG, T-box, NR, and RHD Family TFs (A) Two-dimensional multidimensional scaling plot of ATAC-seq signal for all replicates of naive cells, effector cells, SLECs, MPECs, memory cells, and exhausted cells at 18,043 regions differentially accessible regions identified from comparisons of naive, effector, memory, and exhausted cells. (B) k-means clustered log2 fold-change from mean ATAC-seq signal for all differentially accessible regions identified from comparisons between naive cells, effector cells, SLECs, MPECs, memory cells, and exhausted CD8+ T cells. (C) Enrichment of all known motifs within each cluster of differentially accessible regions compared to all differentially accessible regions in naive, effector, memory, and exhausted CD8+ T cells. All motifs with an enrichment log p value less than –35 and found in 10% or more regions in at least one cluster are shown. (D) Percent of ChIP-seq peaks overlapping differentially accessible regions in each cluster or percent of ATAC-seq peaks in each cluster (left column). The total number of ChIP-seq peaks for each TF and the fraction of these that overlap any of these differentially accessible regions are shown below the plot. (E) log2 fold-change from mean RNA-seq counts per transcript are shown for all expressed TFs from families associated with each enriched motif. (F) MeDIP-seq coverage compared to input for naive and effector CD8+ T cells 8 days after LCMV Arm5 infection. The top graph is for all accessible regions in CD8+ T cells, where each graph below is associated with the clusters indicated at left in (B). ATAC-seq data in (A) and (B) are mean of at least two independent replicates and RNA-seq data in (E) are mean of two independent replicates. See also Figure S5.

Gene expression data indicated that among the HMG-containing TFs, the expression of Lef1 and Tcf7 was biased toward naive and memory cells (Figure 5E). Altogether, these data agree with the well-characterized contribution of these two TFs to the function of naive and memory CD8+ T cells (Steinke et al., 2014; Zhou et al., 2010) and suggested that their contribution is related to binding at these differentially accessible regions. The same clusters, along with clusters 9 and 10 (which also shared naive- and memory-biased accessibility), also overlapped E2A ChIP-seq peaks from mouse thymocytes (Leong et al., 2016). Another major group including clusters 2, 5, 6, 7, and 11 had higher signal in effector, memory, and exhausted cells than naive (Figure 5B). These groups of accessible regions were enriched for bZIP and T-box motifs compared to all differentially accessible regions in CD8+ T cells (Figure 5C) and a substantial fraction of ChIP-seq peaks identified for BATF, IRF4, and T-bet from 8 Immunity 45, 1–14, December 20, 2016

effector CD8+ T cells overlapped these regions (Figure 5D), supporting the established role for these TFs in effector T cell differentiation (Dominguez et al., 2015; Intlekofer et al., 2005; Kurachi et al., 2014). The regions with accessibility biased most toward exhausted cells (clusters 3 and 6) were enriched for consensus binding motifs characteristic of nuclear hormone receptor (NR)- and rel homology domain (RHD)-containing proteins (Figure 5C). Among NR-family TFs, we found that the expression of Nr4a2 was biased toward the exhausted cells compared to naive, effector, and memory cells (Figures 5E and S5A). Among RHD-family TFs, we did not identify any genes that were primarily expressed in exhausted cells, including those encoding NFAT family members (Figure S5B), although Rel was reduced in exhausted cells compared to the other subsets (Figure 5E).

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

We observed similar patterns among all 2,466 differentially expressed genes, with prominent clusters of naive- or exhaustedbiased genes, as well as groups of genes with shared expression between naive and memory cells, effector and exhausted cells, or effector, memory, and exhausted cells (Figure S5C). In an effort to associate the groups of differentially accessible regions with global changes in gene expression, we compared the expression of genes within 25 kb of peaks in each cluster from Figure 5B. The changes in accessibility were largely positively correlated with differences in gene expression, although there were examples of some peak-gene pairs with a negatively correlated association (Figure S5D). For example, some differentially accessible peaks near the Tcf7 locus tended to have high signal in naive and memory cells, which was associated with high Tcf7 expression in these cell types (Figure S5D, clusters 4 and 9), but there were other peaks near Tcf7 with high signal in effector cells despite lower expression in these cells (Figure S5D, clusters 11 and 2). Despite these examples, there were many positive associations between changes in accessibility and activation-associated differentially expressed genes; for instance, Ccl5, Tbx21, Gzma, and Batf all showed accessibility and gene expression biased toward effector and exhausted cells (Figure S5D, clusters 2, 5, 6, 7). To associate accessibility with local epigenetic modifications, we compared the MeDIP-seq signal from naive cells and effector cells after LCMV Arm5 infection (Scharer et al., 2013) around these groups of differentially accessible peaks; this technique assesses DNA methylation by immunoprecipitation with a methylcytosine-specific antibody. Most enhancers have been found to have moderate to low DNA methylation (Thurman et al., 2012) and around the center of all ATAC-seq peaks in both naive and effector cells, we also found low average MeDIPseq coverage compared to input (Figure 5F). Examining the MeDIP-seq signal around the clusters of differentially accessible peaks, we identified several patterns. First, clusters 1, 4, 9, and 10, which had higher ATAC-seq signal in naive and memory T cells than in effector cells, were associated with low MeDIPseq signal in both naive and effector cells, indicating that methylation is not strongly linked with reduced accessibility at these regions in effector cells. Similarly, cluster 2, which had higher accessibility in effector than naive cells, was also associated with low MeDIP-seq signal, suggesting that these regions had low methylation already in naive cells. Clusters 7 and 11 were the only clusters in which increased accessibility from naive to effector was associated with reduced MeDIP-seq signal in effector compared to naive. Finally, clusters 3 and 6, whose accessibility is biased toward exhausted, did not have substantial depletion of MeDIP-seq coverage compared to input; these regions may be demethylated only in exhausted cells, for which MeDIP-seq data are not yet available. NFAT Activity Contributes Directly and Indirectly to the Exhaustion Profile We previously showed that constitutive NFAT activity induces an exhaustion-like profile in T cells (Martinez et al., 2015). To assess how NFAT activity influences chromatin accessibility and the exhaustion phenotype, we expressed a constitutively active form of NFAT (NFAT-CA-RIT) with three amino acid substitutions (R468, I469, and T535) that disrupt the interaction between NFAT

and AP-1 (Figure 1A). For comparison, cells were transduced with a version of NFAT-CA-RIT unable to bind DNA (NFAT-CARIT-DBDmut) or left untransduced (Mock). More than 4,000 regions were differentially accessible between mock-transduced cells and cells expressing NFAT-CA-RIT (Figure 6A); in contrast, expression of NFAT-CA-RIT-DBDmut did not change accessibility substantially (Figure S6A). Although regions with increased accessibility in NFAT-CA-RIT-expressing cells were also more accessible in effector, memory, and exhausted cells than in naive cells, the greatest increase in accessibility was observed in exhausted cells. Conversely, regions with reduced accessibility in NFAT-CA-RIT-expressing cells were least accessible in exhausted cells (Figure 6B). Moreover, regions with higher accessibility in exhausted than in effector or memory cells had higher ATAC-seq signal in cells expressing NFAT-CA-RIT compared to mock-transduced cells (Figure S6B). Focusing on regions with altered accessibility in exhausted cells in vivo, we asked whether they showed higher enrichment for NFAT-CA-RIT binding based on ChIP-seq analysis. Of the 993 regions highly accessible in exhausted cells (compared to both effector and memory cells, Figure 4H), most showed increased accessibility in NFAT-CA-RIT-expressing compared to mock-transduced cells (Figure 6C, bottom, data points shifted to right of the vertical line), many of which had high NFAT-CA-RIT binding based on ChIP-seq (Figure 6C, bottom; data points shifted above the horizontal line). Conversely, of the 762 regions with low accessibility in exhausted cells (compared to both effector and memory, Figure 4I), most showed reduced accessibility in NFAT-CA-RIT-expressing compared to mock-transduced cells (Figure 6C, top, data points shifted to left of the vertical line) and fewer of these sites showed high ChIP-seq signal for NFAT-CA-RIT (Figure 6C, top, data points shifted above the horizontal line). Together these results extend our prior finding that NFAT activates the program of CD8+ T cell exhaustion and indicate that continuous NFAT activity contributes to the exhaustion profile during chronic infection in vivo (Martinez et al., 2015). In contrast, brief restimulation of CD8+ T cells activated and cultured in vitro only partially recapitulated the effect of NFATCA-RIT in inducing increased accessibility at exhaustion-related regions. The combination of phorbol myristate acetate (PMA) and ionomycin (Iono), or Iono alone, induced rapid (within 2 hr) changes in accessibility at 3,600 and 700 regions, respectively (Figures S6C and S6D). These changes were partially blocked by addition of the calcineurin inhibitor cyclosporin (CsA), suggesting NFAT involvement. Unlike changes induced by chronic expression of NFAT-CA-RIT (Figure 6B, right), these accessible regions induced by acute stimulation had similar signals in effector, memory, and exhausted cells (Figures S6C and S6D, right). The changes in accessibility observed at regions with high or low accessibility in exhausted cells (Figures 4H and 4I) were more apparent in NFAT-CA-RIT-expressing cells than in PMA-Iono-stimulated cells (Figure S6E). The Pdcd1 locus illustrates several aspects of the contribution of NFAT to changes in accessibility in exhausted-phenotype cells (Figure 6D). Three accessible regions in the Pdcd1 locus had higher signal in exhausted cells compared to naive, memory, or effector cells: two of these were located 1.1 and 2.7 kb 50 of the Pdcd1 TSS whereas the third was located 22 kb 50 of the Immunity 45, 1–14, December 20, 2016 9

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

A

B

C

D

E

Figure 6. Constitutively Active NFAT Partially Recapitulates the Chronic Activation Profile In Vitro (A) Scatterplot of ATAC-seq counts per peak comparing in vitro cultured CD8+ T cells after transduction with retroviruses expressing the NFAT-CA-RIT mutant or left untransduced (Mock). (B) Boxplots of ATAC-seq counts per peak in naive, effector, memory, and exhausted CD8+ T cells at differentially accessible regions between mock and NFATCA-RIT mutant expressing cells. Box indicates interquartile range with whiskers ±1.5 times this range and outlier points. (C) Scatterplot of NFAT-CA-RIT ChIP-seq coverage with log2 fold-change ATAC-seq signal between mock and NFAT-CA-RIT mutant expressing cells at regions with lower (top) or higher (bottom) ATAC-seq signal in exhausted compared to effector and memory CD8+ T cells. (D) Mean ATAC-seq and NFAT ChIP-seq coverage at the Pdcd1 locus with a scale of 0–1,200 for ATAC-seq tracks. (E) Nr4 family member gene expression in CD8+ T cells overexpressing the NFAT-CA-RIT mutant or left untransduced (Mock) showing mean plus range. ATAC-seq data in (A)–(D) are from at least two independent replicates. See also Figure S6.

TSS (Figure 6D). All three peaks were observed in NFAT-CA-RITexpressing cells, but only the two peaks proximal to the TSS were induced by stimulation with PMA and ionomycin and bound both NFAT and NFAT-CA-RIT. In contrast, the accessible region at +22 kb showed no NFAT binding and minimal change in accessibility after brief stimulation in vitro. Taken together, these data suggest that NFAT activity can directly and indirectly influ10 Immunity 45, 1–14, December 20, 2016

ence the chromatin accessibility profile associated with T cell exhaustion. We compared the consensus TF binding motifs enriched in differentially accessible regions after NFAT-CA-RIT expression or after restimulation of in vitro-cultured cells. Accessible regions associated with NFAT-CA-RIT expression and stimulation with either Iono alone or PMA and Iono together were enriched, as

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

expected, in NFAT and NFAT-AP1 motifs (Figure S6H). In contrast, the NR-family binding motif that was found in exhaustion-biased regions, including the peak located 22 kb 50 of the TSS in the Pdcd1 locus (Figure 6D), was enriched only in accessible regions induced by NFAT-CA-RIT (Figure S6H). Activation of the NR family under these conditions may be a direct effect of NFAT activity, because NFAT-CA-RIT-expressing cells showed increased expression of Nr4a2 and Nr4a3 (Figure 6E) and NFAT-CA-RIT-occupied accessible regions within the Nr4a2 locus (Figure S6I). Altogether, these data indicated that prolonged stimulation during chronic infection was associated with changes in accessibility that were directly and indirectly induced by NFAT activity. DISCUSSION In this study, we compared chromatin accessibility in CD8+ T cells under different conditions of stimulation and differentiation: in effector, memory, and exhausted cells arising in response to in vivo viral infection; in cells pre-activated in culture and assessed before or after brief restimulation in vitro; and in cells subjected to retroviral transduction with constitutively active versions of NFAT. We identified more than 70,000 regions that were accessible in at least one subset of activated or differentiated cells. At least half of these accessible regions had similar signals in most cell types examined, demonstrating that there is a core pattern of chromatin accessibility that is relatively stable in CD8+ T cells regardless of their activation or differentiation state. Many commonly accessible regions were located near the TSS, but most were located within introns or at intergenic regions and could potentially function as enhancers (ENCODE Project Consortium, 2012). Although these regions remain similarly accessible among CD8+ T cell subsets, their function may be regulated by subsequent DNA or histone modifications on flanking nucleosomes (Kouzarides, 2007). Indeed, variable enrichment of the trimethylation of H3K4 and H3K27 at promoters has been associated with differential transcription in naive, effector, and memory CD8+ T cells after influenza infection (Russ et al., 2014). At least 30% of all accessible regions were differentially accessible between any pair of cell types. Most of these regions were located far from TSSs but were linked to differentially expressed genes and may represent enhancers that control CD8+ T cell differentiation and function. Among the populations we have compared here, the largest number of differentially accessible regions were observed between naive and effector T cells: approximately 25% of all accessible regions were differentially accessible between these two cell types. Regions that were differentially accessible between naive, effector, and memory CD8+ T cells were enriched for binding sites of several TFs known to contribute to the state of the relevant cell types. Binding of TCF family members, in particular TCF1, was substantially enriched at regions accessible in naive and memory CD8+ T cells. Conversely, binding sites for bZIP, IRF, and T-box TFs were enriched within accessible regions specific to effector and memory CD8+ T cells. Our analyses did not identify a substantial number of accessible regions specific to memory CD8+ T cells, which were characterized by a mixture of the accessibility profiles observed in

effector and naive cells. Many peaks induced by viral responses at day 8 remained accessible at day 35 in memory cells, providing a potential mechanism for their more rapid responses to antigen compared to naive cells (Williams and Bevan, 2007). However, when compared to effector T cells, memory CD8+ T cells also shared certain aspects of their accessibility profile with naive T cells, accompanied by a corresponding similarity in gene expression, as described previously (Best et al., 2013). Comparison of MPEC and SLEC subsets from day 8 after acute infection revealed that MPECs are substantially similar to SLECs and total effector populations, with a slight bias toward the memory accessibility profile. Our data do not address whether the shift from MPECs to memory cells reflects the selective survival of cells with a memory-like chromatin accessibility profile after viral infection, or whether cells that originally had an effectorlike pattern of chromatin accessibility revert to accessibility profiles that more closely resemble those of memory cells arising later in the infection. Based on comparisons to DNA methylation data from naive and effector T cells after LCMV infection (Scharer et al., 2013), we found low enrichment of DNA methylation around accessible regions in any given cell type, as generally observed in other cells (Thurman et al., 2012). Among regions that were accessible in effector cells after activation, some had a corresponding reduction in DNA methylation from naive to effector, whereas other regions were already demethylated in naive cells. Together, these data support the hypothesis that concurrent regulation of DNA methylation and chromatin accessibility controls the activity of regulatory elements in activated T cells to alter gene expression and influence differentiation. In a chronic infection, virus is not completely cleared; antigenspecific cells are continually subjected to chronic stimulation at the site of infection and they acquire a hyporesponsive or ‘‘exhausted’’ state (Speiser et al., 2014; Wherry and Kurachi, 2015). Based on the chromatin accessibility described here, exhausted cells were most similar to effector cells: for instance, both functional subsets had similar accessibility profiles at the Ifng locus, which is strongly activated in effector cells but poorly induced upon antigen stimulation of exhausted CD8+ T cells (Barber et al., 2006; Kao et al., 2011). These observations may explain the effectiveness of antibody therapies, targeting PD-1, PD-L1, and CTLA4 where blockade of the inhibitory signals transmitted by these receptors may enable exhausted T cells to reactivate enhancers that remained as accessible as in effector cells (Barber et al., 2006; Pardoll, 2012). Regions with similar accessibility between effector and exhausted cells were enriched for T-box motifs and often bound by T-bet in effector cells, consistent with the suggestion that variable expression of T-bet and EOMES among exhausted CD8+ T cells could alter cell function in cells that are similar in their overall accessibility profiles (Doering et al., 2012; Intlekofer et al., 2005; Wherry and Kurachi, 2015). Regions that were uniquely accessible in exhausted cells were enriched for two motifs (NFAT and Nr4 family) that were not substantially enriched in any other viral-specific T cell subset. NFAT can directly induce both Nr4a2 and Nr4a3, whereas only Nr4a2 was induced in exhausted T cells (this study and Martinez et al., 2015). Nr4a family members have been associated with the development and function of regulatory T (Treg) cells (Sekiya Immunity 45, 1–14, December 20, 2016 11

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

et al., 2013); Treg cells also fail to express cytokines characteristic of activated cells after antigen receptor stimulation, and thus Nr4a proteins may mediate similar cell-intrinsic inhibitory programs in Treg cells and exhausted CD8+ cells. We hypothesize that the uniquely exhaustion-specific regions may be induced by continued antigen stimulation during chronic infection, as their accessibility profiles overlap substantially with the profiles of cells activated in vitro and cells expressing constitutively active NFAT, even though only some of the accessible regions bound NFAT in ChIP-seq experiments. Some regions (for instance in the Satb1 and Il7r loci) were less accessible in exhausted cells than in effector and memory cell types, suggesting that at least some aspect of T cell-specific chromatin accessibility is lost in exhausted cells. Loss of chromatin accessibility relative to effector and memory T cells may explain the reported incomplete restoration of function in exhausted cells after checkpoint blockade or after transfer to mice in the absence of antigen (Barber et al., 2006; Schietinger et al., 2012; Utzschneider et al., 2013). Additionally, exhaustion-biased accessible regions showed higher enrichment for DNA methylation in both naive and effector cells, compared to sites accessible in naive and effector cells. These sites may be demethylated specifically in exhausted cells, as shown for two conserved regions in the Pdcd1 locus during chronic infection in mice and humans (Youngblood et al., 2011). Thus DNA methylation is an additional factor that may govern the difference between effector and exhausted phenotypes. The changes in chromatin accessibility observed after activation of naive CD8+ T cells, whether in vitro or in vivo, form a continuum reflecting the diverse ways in which TFs expressed in different cellular subsets can cooperate to influence gene expression and cellular function. Based on their patterns of chromatin accessibility, in vitro activated cells, which have been exposed to stimulatory antibodies with limited cytokine exposure and without antigen-presenting cells, are most similar to naive cells that have never encountered antigen, whereas memory, effector, and exhausted cells, which have been exposed to diverse stimulation in a native environment, are progressively less similar. Future experiments will determine whether continuing antigen stimulation or inflammation induces stable changes in chromatin accessibility across an entire population, or whether stochastic changes in accessibility and TF binding and expression give rise to a heterogeneous cell population from which cells with particular phenotypes survive. EXPERIMENTAL PROCEDURES Mice CD8+ T cells were purified from splenocytes of 6- to 10-week-old C57BL/6J mice after infection with LCMV or left uninfected. Effector and memory CD8+ T cell subsets were collected after infection with LCMV Arm5, based on staining with tetramer and antibodies to CD44, IL7Ra, and KLRG1 (Figures S1A, S1C, and S1D). Exhausted cells were collected after infection LCMV clone 13 based on staining with tetramer and antibodies to PD-1 and Tim3 (Figure S1E). For in vitro stimulation, CD8+ T cells were purified from P14 T cell receptor transgenic mice on a TCRa constant region-deficient background, before stimulation with CD3ε- and CD28-specific antibodies for 2 days and then expanded 10 U/mL or 100 U/mL IL-2 (Pipkin et al., 2010). On day 6, cells were restimulated for 2 hr under various conditions (Supplemental Experimental Procedures). For retroviral transduction experiments, CD8+ T cells

12 Immunity 45, 1–14, December 20, 2016

were isolated by negative selection from C57BL/6J mice. These cells were stimulated for 1 day with antibodies specific for CD3ε and CD28 before infection with retroviruses followed by expansion for three days in 10 U/mL IL-2. All experiments were performed according to protocols approved by the La Jolla Institute animal care and use committee. ATAC-Seq ATAC-seq library preparations were performed as described (Buenrostro et al., 2013). Cells were washed with PBS before treatment with lysis buffer followed by labeling with Nextera enzyme (Illumina), before PCR amplification with 10–12 cycles with barcoded primers and 23 50 cycle paired-end sequencing (Illumina). Reads were mapped to mouse genome (mm9) using bowtie (v.1.0.0) (Langmead et al., 2009). Unmapped reads were processed with trim_galore, re-mapped with bowtie, and merged with previous mapping output. Duplicate reads identified by picard MarkDuplicates and reads mapping to chrM were excluded. Peak summits were identified with MACS2 (Zhang et al., 2008) from individual replicates using short DNA fragments (<100 bp) and expanded to 500 bp regions. We merged overlapping regions from all replicates and excluded regions on chrY or those that intersected ENCODE blacklisted regions. The number of transposase insertions within each region was computed for each replicate and differentially accessible regions were identified with DESeq2 (v.1.6) based on an fdr adjusted p value of less than 1 3 10
3 and an estimated fold change of at least 3. Full details are available in Supplemental Experimental Procedures. E2A (Leong et al., 2016) (GEO: GSE84974), TCF1 (Steinke et al., 2014) (GEO: GSE52070), T-bet (Dominguez et al., 2015) (SRA: SRR2075567-8, SRR2075585), BATF, and IRF4 (Kurachi et al., 2014) (GEO: GSE54191) ChIP-seq and naive and effector MeDIP-seq (Scharer et al., 2013) (GEO: GSE44638), including input controls where available, were downloaded from NCBI SRA database and mapped to the mm9 genome (bowtie v.1.0.0). ChIP-seq peaks were called with MACS2. RNA-Seq SMARTseq2 RNA-seq libraries were prepared from naive, effector, memory, and exhausted cells (Picelli et al., 2014) sequenced with an Illumina HiSeq2500 and mapped to the mm9 genome (TopHat 1.4.1) (Trapnell et al., 2009). Lowquality reads were excluded and feature counts were computed with htseqcount (0.6.0 using the ‘‘union’’ option). Differentially expressed genes were identified with DESeq2 with an fdr adjusted p value of less than 1 3 10
3 and an estimated fold-change of at least 2 and a mean of least 10 transcripts per million transcripts. Full details are available in Supplemental Experimental Procedures. ACCESSION NUMBERS The accession number for the ATAC-seq and RNA-seq data reported in this paper is GEO: GSE88987. SUPPLEMENTAL INFORMATION Supplemental Information includes six figures, one table, and Supplemental Experimental Procedures and can be found with this article online at http:// dx.doi.org/10.1016/j.immuni.2016.10.028. AUTHOR CONTRIBUTIONS J.P.S.-B., R.M.P., S.T., V.W., and A.R. designed experiments; J.P.S.-B., R.M.P., V.W., and S.T. performed experiments; J.P.S.-B., I.F.L.-M., and L.C. analyzed data; and J.P.S.-B., R.M.P., and A.R. wrote the manuscript. ACKNOWLEDGMENTS We thank C. Kim, K. Gunst, and L. Nosworthy for cell sorting, J. Greenbaum, Z. Fu, and J. Day for next-generation sequencing, S. Crotty for LCMV virus (all at La Jolla Institute), the NIH Tetramer Facility for tetramers, and L. Gapin (University of Colorado), P. Marrack (National Jewish Health), and M. Prlic (FHCRC) for providing critical comments. This work was funded by NIH grants R01

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

AI109842 and AI40127 (to A.R.). J.P.S.-B. was the Fraternal Order of Eagles Fellow of the Damon Runyon Cancer Research Foundation, DRG-2069-11. R.M.P. was supported by a fellowship from the Pew Latin American Fellows Program in the Biomedical Sciences. I.F.L.-M. was supported by a UC MEXUS-CONACYT fellowship. Received: May 6, 2016 Revised: September 8, 2016 Accepted: September 23, 2016 Published: December 6, 2016 REFERENCES Barber, D.L., Wherry, E.J., Masopust, D., Zhu, B., Allison, J.P., Sharpe, A.H., Freeman, G.J., and Ahmed, R. (2006). Restoring function in exhausted CD8 T cells during chronic viral infection. Nature 439, 682–687. Best, J.A., Blair, D.A., Knell, J., Yang, E., Mayya, V., Doedens, A., Dustin, M.L., and Goldrath, A.W.; Immunological Genome Project Consortium (2013). Transcriptional insights into the CD8(+) T cell response to infection and memory T cell formation. Nat. Immunol. 14, 404–412. Blackburn, S.D., Shin, H., Haining, W.N., Zou, T., Workman, C.J., Polley, A., Betts, M.R., Freeman, G.J., Vignali, D.A., and Wherry, E.J. (2009). Coregulation of CD8+ T cell exhaustion by multiple inhibitory receptors during chronic viral infection. Nat. Immunol. 10, 29–37. Buenrostro, J.D., Giresi, P.G., Zaba, L.C., Chang, H.Y., and Greenleaf, W.J. (2013). Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat. Methods 10, 1213–1218. ENCODE Project Consortium (2012). An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74. Crespo, J., Sun, H., Welling, T.H., Tian, Z., and Zou, W. (2013). T cell anergy, exhaustion, senescence, and stemness in the tumor microenvironment. Curr. Opin. Immunol. 25, 214–221. Doering, T.A., Crawford, A., Angelosanto, J.M., Paley, M.A., Ziegler, C.G., and Wherry, E.J. (2012). Network analysis reveals centrally connected genes and pathways involved in CD8+ T cell exhaustion versus memory. Immunity 37, 1130–1144. Dominguez, C.X., Amezquita, R.A., Guan, T., Marshall, H.D., Joshi, N.S., Kleinstein, S.H., and Kaech, S.M. (2015). The transcription factors ZEB2 and T-bet cooperate to program cytotoxic T cell terminal differentiation in response to LCMV viral infection. J. Exp. Med. 212, 2041–2056. Intlekofer, A.M., Takemoto, N., Wherry, E.J., Longworth, S.A., Northrup, J.T., Palanivel, V.R., Mullen, A.C., Gasink, C.R., Kaech, S.M., Miller, J.D., et al. (2005). Effector and memory CD8+ T cell fate coupled by T-bet and eomesodermin. Nat. Immunol. 6, 1236–1244. Kao, C., Oestreich, K.J., Paley, M.A., Crawford, A., Angelosanto, J.M., Ali, M.A., Intlekofer, A.M., Boss, J.M., Reiner, S.L., Weinmann, A.S., and Wherry, E.J. (2011). Transcription factor T-bet represses expression of the inhibitory receptor PD-1 and sustains virus-specific CD8+ T cell responses during chronic infection. Nat. Immunol. 12, 663–671. Kouzarides, T. (2007). Chromatin modifications and their function. Cell 128, 693–705. Kurachi, M., Barnitz, R.A., Yosef, N., Odorizzi, P.M., DiIorio, M.A., Lemieux, M.E., Yates, K., Godec, J., Klatt, M.G., Regev, A., et al. (2014). The transcription factor BATF operates as an essential differentiation checkpoint in early effector CD8+ T cells. Nat. Immunol. 15, 373–383. Langmead, B., Trapnell, C., Pop, M., and Salzberg, S.L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25. Leong, Y.A., Chen, Y., Ong, H.S., Wu, D., Man, K., Deleage, C., Minnich, M., Meckiff, B.J., Wei, Y., Hou, Z., et al. (2016). CXCR5(+) follicular cytotoxic T cells control viral infection in B cell follicles. Nat. Immunol. 17, 1187–1196. Love, M.I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550.

Martinez, G.J., Pereira, R.M., Aijo, T., Kim, E.Y., Marangoni, F., Pipkin, M.E., Togher, S., Heissmeyer, V., Zhang, Y.C., Crotty, S., et al. (2015). The transcription factor NFAT promotes exhaustion of activated CD8(+) T cells. Immunity 42, 265–278. Moskophidis, D., Lechner, F., Pircher, H., and Zinkernagel, R.M. (1993). Virus persistence in acutely infected immunocompetent mice by exhaustion of antiviral cytotoxic effector T cells. Nature 362, 758–761. Pardoll, D.M. (2012). The blockade of immune checkpoints in cancer immunotherapy. Nat. Rev. Cancer 12, 252–264. Picelli, S., Faridani, O.R., Bjo¨rklund, A.K., Winberg, G., Sagasser, S., and Sandberg, R. (2014). Full-length RNA-seq from single cells using Smartseq2. Nat. Protoc. 9, 171–181. Pipkin, M.E., Sacks, J.A., Cruz-Guilloty, F., Lichtenheld, M.G., Bevan, M.J., and Rao, A. (2010). Interleukin-2 and inflammation induce distinct transcriptional programs that promote the differentiation of effector cytolytic T cells. Immunity 32, 79–90. Russ, B.E., Olshanksy, M., Smallwood, H.S., Li, J., Denton, A.E., Prier, J.E., Stock, A.T., Croom, H.A., Cullen, J.G., Nguyen, M.L., et al. (2014). Distinct epigenetic signatures delineate transcriptional programs during virus-specific CD8(+) T cell differentiation. Immunity 41, 853–865. Scharer, C.D., Barwick, B.G., Youngblood, B.A., Ahmed, R., and Boss, J.M. (2013). Global DNA methylation remodeling accompanies CD8 T cell effector function. J. Immunol. 191, 3419–3429. Schietinger, A., Delrow, J.J., Basom, R.S., Blattman, J.N., and Greenberg, P.D. (2012). Rescued tolerant CD8 T cells are preprogrammed to reestablish the tolerant state. Science 335, 723–727. Schietinger, A., Philip, M., Krisnawan, V.E., Chiu, E.Y., Delrow, J.J., Basom, €mmerling, G.J., et al. R.S., Lauer, P., Brockstedt, D.G., Knoblaugh, S.E., Ha (2016). Tumor-specific T cell dysfunction is a dynamic antigen-driven differentiation program initiated early during tumorigenesis. Immunity 45, 389–401. Sekiya, T., Kashiwagi, I., Yoshida, R., Fukaya, T., Morita, R., Kimura, A., Ichinose, H., Metzger, D., Chambon, P., and Yoshimura, A. (2013). Nr4a receptors are essential for thymic regulatory T cell development and immune homeostasis. Nat. Immunol. 14, 230–237. €nz, C., Romero, P., and Speiser, D.E., Utzschneider, D.T., Oberle, S.G., Mu Zehn, D. (2014). T cell differentiation in chronic infection and cancer: functional adaptation or exhaustion? Nat. Rev. Immunol. 14, 768–774. Steinke, F.C., Yu, S., Zhou, X., He, B., Yang, W., Zhou, B., Kawamoto, H., Zhu, J., Tan, K., and Xue, H.H. (2014). TCF-1 and LEF-1 act upstream of Th-POK to promote the CD4(+) T cell fate and interact with Runx3 to silence Cd4 in CD8(+) T cells. Nat. Immunol. 15, 646–656. Thurman, R.E., Rynes, E., Humbert, R., Vierstra, J., Maurano, M.T., Haugen, E., Sheffield, N.C., Stergachis, A.B., Wang, H., Vernot, B., et al. (2012). The accessible chromatin landscape of the human genome. Nature 489, 75–82. Trapnell, C., Pachter, L., and Salzberg, S.L. (2009). TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105–1111. Utzschneider, D.T., Legat, A., Fuertes Marraco, S.A., Carrie´, L., Luescher, I., Speiser, D.E., and Zehn, D. (2013). T cells maintain an exhausted phenotype after antigen withdrawal and population reexpansion. Nat. Immunol. 14, 603–610. Vierstra, J., Rynes, E., Sandstrom, R., Zhang, M., Canfield, T., Hansen, R.S., Stehling-Sun, S., Sabo, P.J., Byron, R., Humbert, R., et al. (2014). Mouse regulatory DNA landscapes reveal global principles of cis-regulatory evolution. Science 346, 1007–1012. Wherry, E.J. (2011). T cell exhaustion. Nat. Immunol. 12, 492–499. Wherry, E.J., and Ahmed, R. (2004). Memory CD8 T-cell differentiation during viral infection. J. Virol. 78, 5535–5545. Wherry, E.J., and Kurachi, M. (2015). Molecular and cellular insights into T cell exhaustion. Nat. Rev. Immunol. 15, 486–499. €ber, V., Becker, T.C., Masopust, D., Kaech, S.M., Antia, Wherry, E.J., Teichgra R., von Andrian, U.H., and Ahmed, R. (2003). Lineage relationship and protective immunity of memory CD8 T cell subsets. Nat. Immunol. 4, 225–234.

Immunity 45, 1–14, December 20, 2016 13

Please cite this article in press as: Scott-Browne et al., Dynamic Changes in Chromatin Accessibility Occur in CD8+ T Cells Responding to Viral Infection, Immunity (2016), http://dx.doi.org/10.1016/j.immuni.2016.10.028

Williams, M.A., and Bevan, M.J. (2007). Effector and memory CTL differentiation. Annu. Rev. Immunol. 25, 171–192. Yamamoto, T., Price, D.A., Casazza, J.P., Ferrari, G., Nason, M., Chattopadhyay, P.K., Roederer, M., Gostick, E., Katsikis, P.D., Douek, D.C., et al. (2011). Surface expression patterns of negative regulatory molecules identify determinants of virus-specific CD8+ T-cell exhaustion in HIV infection. Blood 117, 4805–4815. Youngblood, B., Oestreich, K.J., Ha, S.J., Duraiswamy, J., Akondy, R.S., West, E.E., Wei, Z., Lu, P., Austin, J.W., Riley, J.L., et al. (2011). Chronic virus

14 Immunity 45, 1–14, December 20, 2016

infection enforces demethylation of the locus that encodes PD-1 in antigenspecific CD8(+) T cells. Immunity 35, 400–412. Zhang, Y., Liu, T., Meyer, C.A., Eeckhoute, J., Johnson, D.S., Bernstein, B.E., Nusbaum, C., Myers, R.M., Brown, M., Li, W., and Liu, X.S. (2008). Modelbased analysis of ChIP-Seq (MACS). Genome Biol. 9, R137. Zhou, X., Yu, S., Zhao, D.M., Harty, J.T., Badovinac, V.P., and Xue, H.H. (2010). Differentiation and persistence of memory CD8(+) T cells depend on T cell factor 1. Immunity 33, 229–240.