Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Ligand-free mitochondria-localized mutant AR-induced cytotoxicity in spinal bulbar muscular atrophy.

Spinal bulbar muscular atrophy (SBMA), the first identified CAG-repeat expansion disorder, is an X-linked neuromuscular disorder involving CAG-repeat-expansion mutations in the androgen receptor (AR) gene. We utilized CRISPR-Cas9 gene editing to engineer novel isogenic human induced pluripotent stem cell (hiPSC) models, consisting of isogenic AR knockout, control and disease lines expressing mutant AR with distinct repeat lengths, as well as control and disease lines expressing FLAG-tagged wild-type and mutant AR, respectively. Adapting a small-molecule cocktail-directed approach, we differentiate the isogenic hiPSC models into motor neuron-like cells with a highly enriched population to uncover cell-type-specific mechanisms underlying SBMA and to distinguish gain- from loss-of-function properties of mutant AR in disease motor neurons. We demonstrate that ligand-free mutant AR causes drastic mitochondrial dysfunction in neurites of differentiated disease motor neurons due to gain-of-function mechanisms and such cytotoxicity can be amplified upon ligand (androgens) treatment. We further show that aberrant interaction between ligand-free, mitochondria-localized mutant AR and F-ATP synthase is associated with compromised mitochondrial respiration and multiple other mitochondrial impairments. These findings counter the established notion that androgens are requisite for mutant AR-induced cytotoxicity in SBMA, reveal a compelling mechanistic link between ligand-free mutant AR, F-ATP synthase and mitochondrial dysfunction, and provide innovative insights into motor neuron-specific therapeutic interventions for SBMA.

Snapshot: a package for clustering and visualizing epigenetic history during cell differentiation.

BACKGROUND: Epigenetic modification of chromatin plays a pivotal role in regulating gene expression during cell differentiation. The scale and complexity of epigenetic data pose significant challenges for biologists to identify the regulatory events controlling cell differentiation. RESULTS: To reduce the complexity, we developed a package, called Snapshot, for clustering and visualizing candidate cis-regulatory elements (cCREs) based on their epigenetic signals during cell differentiation. This package first introduces a binarized indexing strategy for clustering the cCREs. It then provides a series of easily interpretable figures for visualizing the signal and epigenetic state patterns of the cCREs clusters during the cell differentiation. It can also use different hierarchies of cell types to highlight the epigenetic history specific to any particular cell lineage. We demonstrate the utility of Snapshot using data from a consortium project for ValIdated Systematic IntegratiON (VISION) of epigenomic data in hematopoiesis. CONCLUSION: The package Snapshot can identify all distinct clusters of genomic locations with unique epigenetic signal patterns during cell differentiation. It outperforms other methods in terms of interpreting and reproducing the identified cCREs clusters. The package of Snapshot is available at GitHub: https://github.com/guanjue/Snapshot .

Functional and epigenetic studies reveal multistep differentiation and plasticity of in vitro-generated and in vivo-derived follicular T helper cells.

Follicular T helper (Tfh) cells provide critical help to B cells for germinal center (GC) formation. Mutations affecting SLAM-associated protein (SAP) prevent GC formation because of defective T cell-B cell interactions, yet effects on Tfh cell differentiation remain unclear. We describe the in vitro differentiation of functionally competent "Tfh-like" cells that expressed interleukin-21, Tfh cell markers, and Bcl6 and rescued GC formation in SAP-deficient hosts better than other T helper (Th) cells. SAP-deficient Tfh-like cells appeared virtually indistinguishable from wild-type, yet failed to support GCs in vivo. Interestingly, both Tfh-like and in vivo-derived Tfh cells could produce effector cytokines in response to polarizing conditions. Moreover, Th1, Th2, and Th17 cells could be reprogrammed to obtain Tfh cell characteristics. ChIP-Seq analyses revealed positive epigenetic markings on Tbx21, Gata3, and Rorc in Tfh-like and ex vivo Tfh cells and on Bcl6 in non-Tfh cells, supporting the concept of plasticity between Tfh and other Th cell populations.

Molecular convergence in ex vivo models of Diamond-Blackfan anemia.

Diamond-Blackfan anemia (DBA) is a congenital bone marrow failure syndrome characterized by erythroid hypoplasia, usually without perturbation of other hematopoietic lineages. Approximately 65% of DBA patients with autosomal dominant inheritance have heterozygous mutations or deletions in ribosomal protein (RP) genes while <1% of patients with X-linked inheritance have been identified with mutations in the transcription factor GATA1 Erythroid cells from patients with DBA have not been well characterized, and the mechanisms underlying the erythroid specific effects of either RP or GATA1 associated DBA remain unclear. We have developed an ex vivo culture system to expand peripheral blood CD34+ progenitor cells from patients with DBA and differentiate them into erythroid cells. Cells from patients with RP or GATA1 mutations showed decreased proliferation and delayed erythroid differentiation in comparison with controls. RNA transcript analyses of erythroid cells from controls and patients with RP or GATA1 mutations showed distinctive differences, with upregulation of heme biosynthesis genes prominently in RP-mediated DBA and failure to upregulate components of the translational apparatus in GATA1-mediated DBA. Our data show that dysregulation of translation is a common feature of DBA caused by both RP and GATA1 mutations. This trial was registered at www.clinicaltrials.gov as #NCT00106015.

Differential expression of interleukin-17A and -17F is coupled to T cell receptor signaling via inducible T cell kinase.

T helper 17 (Th17) cells play major roles in autoimmunity and bacterial infections, yet how T cell receptor (TCR) signaling affects Th17 cell differentiation is relatively unknown. We demonstrate that CD4(+) T cells lacking Itk, a tyrosine kinase required for full TCR-induced phospholipase C-gamma (PLC-gamma1) activation, exhibit decreased interleukin-17A (IL-17A) expression in vitro and in vivo, despite relatively normal expression of retinoic acid receptor-related orphan receptor-gammaT (ROR-gammaT) and IL-17F. IL-17A expression was rescued by pharmacologically induced Ca(2+) influx or constitutively activated nuclear factor of activated T cells (NFAT). Conversely, decreased TCR stimulation or calcineurin inhibition preferentially reduced IL-17A expression. We further found that the promoter of Il17a but not Il17f has a conserved NFAT binding site that bound NFATc1 in wild-type but not Itk-deficient cells, even though both exhibited open chromatin conformations. Finally, Itk(-/-) mice also showed differential regulation of IL-17A and IL-17F in vivo. Our results suggest that Itk specifically couples TCR signaling to Il17a expression and the differential regulation of Th17 cell cytokines through NFATc1.

Cystic cerebellar dysplasia and biallelic LAMA1 mutations: a lamininopathy associated with tics, obsessive compulsive traits and myopia due to cell adhesion and migration defects.

BACKGROUND: Laminins are heterotrimeric complexes, consisting of α, ß and γ subunits that form a major component of basement membranes and extracellular matrix. Laminin complexes have different, but often overlapping, distributions and functions. METHODS: Under our clinical protocol, NCT00068224, we have performed extensive clinical and neuropsychiatric phenotyping, neuroimaging and molecular analysis in patients with laminin α1 (LAMA1)-associated lamininopathy. We investigated the consequence of mutations in LAMA1 using patient-derived fibroblasts and neuronal cells derived from neuronal stem cells. RESULTS: In this paper we describe individuals with biallelic mutations in LAMA1, all of whom had the cerebellar dysplasia, myopia and retinal dystrophy, in addition to obsessive compulsive traits, tics and anxiety. Patient-derived fibroblasts have impaired adhesion, reduced migration, abnormal morphology and increased apoptosis due to impaired activation of Cdc42, a member of the Rho family of GTPases that is involved in cytoskeletal dynamics. LAMA1 knockdown in human neuronal cells also showed abnormal morphology and filopodia formation, supporting the importance of LAMA1 in neuronal migration, and marking these cells potentially useful tools for disease modelling and therapeutic target discovery. CONCLUSION: This paper broadens the phenotypes associated with LAMA1 mutations. We demonstrate that LAMA1 deficiency can lead to alteration in cytoskeletal dynamics, which may invariably lead to alteration in dendrite growth and axonal formation. Estimation of disease prevalence based on population studies in LAMA1 reveals a prevalence of 1-20 in 1 000 000. TRIAL REGISTRATION NUMBER: NCT00068224.

Lineage and species-specific long noncoding RNAs during erythro-megakaryocytic development.

Mammals express thousands of long noncoding (lnc) RNAs, a few of which are known to function in tissue development. However, the entire repertoire of lncRNAs in most tissues and species is not defined. Indeed, most lncRNAs are not conserved, raising questions about function. We used RNA sequencing to identify 1109 polyadenylated lncRNAs expressed in erythroblasts, megakaryocytes, and megakaryocyte-erythroid precursors of mice, and 594 in erythroblasts of humans. More than half of these lncRNAs were unannotated, emphasizing the opportunity for new discovery through studies of specialized cell types. Analysis of the mouse erythro-megakaryocytic polyadenylated lncRNA transcriptome indicates that ~75% arise from promoters and 25% from enhancers, many of which are regulated by key transcription factors including GATA1 and TAL1. Erythroid lncRNA expression is largely conserved among 8 different mouse strains, yet only 15% of mouse lncRNAs are expressed in humans and vice versa, reflecting dramatic species-specificity. RNA interference assays of 21 abundant erythroid-specific murine lncRNAs in primary mouse erythroid precursors identified 7 whose knockdown inhibited terminal erythroid maturation. At least 6 of these 7 functional lncRNAs have no detectable expression in human erythroblasts, suggesting that lack of conservation between mammalian species does not predict lack of function.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis -regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our V al i dated S ystematic I ntegrati on (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

A somatic cell defect is associated with the onset of neurological symptoms in a lysosomal storage disease.

Mutations in individuals with the lysosomal storage disorder Niemann-Pick disease, type C1 (NPC1) are heterogeneous, not localized to specific protein domains, and not correlated to time of onset or disease severity. We demonstrate direct correlation of the time of neurological symptom onset with the severity of lysosomal defects in NPC1 patient-derived fibroblasts. This is a novel assay for NPC1 individuals that may be predictive of NPC1 disease progression and broadly applicable to other lysosomal disorders.

Defective inhibition of B-cell proliferation by Wiskott-Aldrich syndrome protein-deficient regulatory T cells.

Wiskott-Aldrich syndrome (WAS) is an inherited immunodeficiency characterized by high incidence of autoantibody-mediated autoimmune complications. Such a feature has been associated with defective suppressor activity of WAS protein-deficient, naturally occurring CD4(+)CD25(+)Foxp3(+) regulatory T cells on responder T cells. However, it remains to be established whether the altered B-cell tolerance reported in WAS patients and Was knockout (WKO) mice is secondary to abnormalities in the direct suppression of B-cell function by nTreg cells or to impaired regulation of T-helper function. Because activated nTreg cells are known to induce granzyme B-mediated B-cell killing, we decided to evaluate the regulatory capabilities of WKO nTregs on B lymphocytes. We found that preactivated WKO nTreg cells failed to effectively suppress B-cell proliferation and that such a defect was associated with reduced killing of B cells and significantly decreased degranulation of granzyme B. Altogether, these results provide additional mechanistic insights into the loss of immune tolerance in WAS.

Genome-wide ChIP-Seq reveals a dramatic shift in the binding of the transcription factor erythroid Kruppel-like factor during erythrocyte differentiation.

Erythropoiesis is dependent on the activity of transcription factors, including the erythroid-specific erythroid Kruppel-like factor (EKLF). ChIP followed by massively parallel sequencing (ChIP-Seq) is a powerful, unbiased method to map trans-factor occupancy. We used ChIP-Seq to study the interactome of EKLF in mouse erythroid progenitor cells and more differentiated erythroblasts. We correlated these results with the nuclear distribution of EKLF, RNA-Seq analysis of the transcriptome, and the occupancy of other erythroid transcription factors. In progenitor cells, EKLF is found predominantly at the periphery of the nucleus, where EKLF primarily occupies the promoter regions of genes and acts as a transcriptional activator. In erythroblasts, EKLF is distributed throughout the nucleus, and erythroblast-specific EKLF occupancy is predominantly in intragenic regions. In progenitor cells, EKLF modulates general cell growth and cell cycle regulatory pathways, whereas in erythroblasts EKLF is associated with repression of these pathways. The EKLF interactome shows very little overlap with the interactomes of GATA1, GATA2, or TAL1, leading to a model in which EKLF directs programs that are independent of those regulated by the GATA factors or TAL1.

Autoimmune lymphoproliferative syndrome due to FAS mutations outside the signal-transducing death domain: molecular mechanisms and clinical penetrance.

PURPOSE: Autoimmune lymphoproliferative syndrome is a disorder of lymphocyte apoptosis. Although FAS molecules bearing mutations in the signal-transducing intracellular death domain exhibit dominant-negative interference with FAS-mediated apoptosis, mechanisms for pathology of non-death domain FAS mutations causing autoimmune lymphoproliferative syndrome are poorly defined. METHODS: RNA stability, protein expression, ligand binding, and ability to transmit apoptosis signals by anti-FAS antibody or FAS ligand were determined for a cohort of 39 patients with non-death domain autoimmune lymphoproliferative syndrome. Correlations between mutation type and disease penetrance were established in mutation-positive family members. RESULTS: Frameshifts or transcriptional stop mutations before exon 7 resulted in messenger RNA haploinsufficiency, whereas an amino-terminal signal sequence mutation and certain intracellular truncations prevented cell surface localization of FAS. All resulted in decreased FAS localization, inability to bind FAS ligand, and reduced FAS ligand-induced apoptosis. Extracellular missense mutations and in-frame deletions expressed defective FAS protein, failed to bind FAS ligand, and exhibited dominant-negative interference with FAS-mediated apoptosis. Mutation-positive relatives with haploinsufficient or extracellular mutations had lower penetrance of autoimmune lymphoproliferative syndrome clinical phenotypes than did relatives with death domain mutations. CONCLUSION: We have defined molecular mechanisms by which non-death domain FAS mutations result in reduced lymphocyte apoptosis, established a hierarchy of genotype-phenotype correlation among mutation-positive relatives of patients with autoimmune lymphoproliferative syndrome, and demonstrated that FAS haploinsufficiency can lead to autoimmune lymphoproliferative syndrome.

Effects of sheared chromatin length on ChIP-seq quality and sensitivity.

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genome-wide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low-quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types.

Establishment of regulatory elements during erythro-megakaryopoiesis identifies hematopoietic lineage-commitment points.

BACKGROUND: Enhancers and promoters are cis-acting regulatory elements associated with lineage-specific gene expression. Previous studies showed that different categories of active regulatory elements are in regions of open chromatin, and each category is associated with a specific subset of post-translationally marked histones. These regulatory elements are systematically activated and repressed to promote commitment of hematopoietic stem cells along separate differentiation paths, including the closely related erythrocyte (ERY) and megakaryocyte (MK) lineages. However, the order in which these decisions are made remains unclear. RESULTS: To characterize the order of cell fate decisions during hematopoiesis, we collected primary cells from mouse bone marrow and isolated 10 hematopoietic populations to generate transcriptomes and genome-wide maps of chromatin accessibility and histone H3 acetylated at lysine 27 binding (H3K27ac). Principle component analysis of transcriptional and open chromatin profiles demonstrated that cells of the megakaryocyte lineage group closely with multipotent progenitor populations, whereas erythroid cells form a separate group distinct from other populations. Using H3K27ac and open chromatin profiles, we showed that 89% of immature MK (iMK)-specific active regulatory regions are present in the most primitive hematopoietic cells, 46% of which contain active enhancer marks. These candidate active enhancers are enriched for transcription factor binding site motifs for megakaryopoiesis-essential proteins, including ERG and ETS1. In comparison, only 64% of ERY-specific active regulatory regions are present in the most primitive hematopoietic cells, 20% of which containing active enhancer marks. These regions were not enriched for any transcription factor consensus sequences. Incorporation of genome-wide DNA methylation identified significant levels of de novo methylation in iMK, but not ERY. CONCLUSIONS: Our results demonstrate that megakaryopoietic profiles are established early in hematopoiesis and are present in the majority of the hematopoietic progenitor population. However, megakaryopoiesis does not constitute a "default" differentiation pathway, as extensive de novo DNA methylation accompanies megakaryopoietic commitment. In contrast, erythropoietic profiles are not established until a later stage of hematopoiesis, and require more dramatic changes to the transcriptional and epigenetic programs. These data provide important insights into lineage commitment and can contribute to ongoing studies related to diseases associated with differentiation defects.

Second-site mutation in the Wiskott-Aldrich syndrome (WAS) protein gene causes somatic mosaicism in two WAS siblings.

Revertant mosaicism due to true back mutations or second-site mutations has been identified in several inherited disorders. The occurrence of revertants is considered rare, and the underlying genetic mechanisms remain mostly unknown. Here we describe somatic mosaicism in two brothers affected with Wiskott-Aldrich syndrome (WAS). The original mutation causing disease in this family is a single base insertion (1305insG) in the WAS protein (WASP) gene, which results in frameshift and abrogates protein expression. Both patients, however, showed expression of WASP in a fraction of their T cells that were demonstrated to carry a second-site mutation causing the deletion of 19 nucleotides from nucleotide 1299 to 1316. This deletion abrogated the effects of the original mutation and restored the WASP reading frame. In vitro expression studies indicated that mutant protein encoded by the second-site mutation was expressed and functional, since it was able to bind to cellular partners and mediate T cell receptor/CD3 downregulation. These observations were consistent with evidence of in vivo selective advantage of WASP-expressing lymphocytes. Molecular analysis revealed that the sequence surrounding the deletion contained two 4-bp direct repeats and that a hairpin structure could be formed by five GC pairs within the deleted fragment. These findings strongly suggest that slipped mispairing was the cause of this second-site mutation and that selective accumulation of WASP-expressing T lymphocytes led to revertant mosaicism in these patients.

Variegated expression from the murine band 3 (AE1) promoter in transgenic mice is associated with mRNA transcript initiation at upstream start sites and can be suppressed by the addition of the chicken beta-globin 5' HS4 insulator element.

The anion exchanger protein 1 (AE1; band 3) is an abundant erythrocyte transmembrane protein that regulates chloride-bicarbonate exchange and provides an attachment site for the erythrocyte membrane skeleton on the cytoplasmic domain. We analyzed the function of the erythroid AE1 gene promoter by using run-on transcription, RNase protection, transient transfection, and transgenic mouse assays. AE1 mRNA was transcribed at a higher level and maintained at a higher steady-state level than either ankyrin or beta-spectrin in mouse fetal liver cells. When linked to a human gamma-globin gene, two different AE1 promoters directed erythroid-specific expression of gamma-globin mRNA in 18 of 18 lines of transgenic mice. However, variegated expression of gamma-globin was observed in 14 of 18 lines. While there was a significant correlation between transgene copy number and the amount of gamma-globin mRNA in all 18 lines, the transgene mRNAs initiated upstream of the start site of the endogenous AE1 mRNA. Addition of the insulator element from 5'HS4 of the chicken beta-globin cluster to the AE1/gamma-globin transgene allowed position-independent, copy-number-dependent expression at levels similar to the AE1 transcription rate in six of six lines of transgenic mice. The mRNA from the insulated AE1/gamma-globin transgene mapped to the start site of the endogenous AE1 mRNA, and gamma-globin protein was expressed in 100% of erythrocytes in all lines. We conclude that the chicken beta-globin 5'HS4 element is necessary for full function of the AE1 promoter and that position effect variegation is associated with RNA transcription from the upstream start sites.

Itk is required for Th9 differentiation via TCR-mediated induction of IL-2 and IRF4.

Th9 cells produce interleukin (IL)-9, a cytokine implicated in allergic asthma and autoimmunity. Here we show that Itk, a mediator of T cell receptor signalling required for Th2 immune responses and the development of asthma, is a positive regulator of Th9 differentiation. In a model of allergic lung disease, Itk-deficient mice show reduced pulmonary inflammation and IL-9 production by T cells and innate lymphoid type 2 cells (ILC2), despite normal early induction of ILC2s. In vitro, Itk(-/-) CD4(+) T cells do not produce IL-9 and have reduced levels of IRF4 (Interferon Regulator Factor 4), a critical transcription factor for effector T cell function. Both IL-9 and IRF4 expression are rescued by either IL-2 or constitutively active STAT5, but not NFATc1. STAT5 binds the Irf4 promoter, demonstrating one mechanism by which IL-2 rescues weakly activated T cells. Itk inhibition also reduces IL-9 expression by human T cells, implicating ITK as a key regulator of Th9 induction.

The TCF1-Bcl6 axis counteracts type I interferon to repress exhaustion and maintain T cell stemness.

During chronic viral infections and in cancer, T cells become dysfunctional, a state known as T cell exhaustion. Although it is well recognized that memory CD8 T cells account for the persistence of CD8 T cell immunity after acute infection, how exhausted T cells persist remains less clear. Using chronic infection with lymphocytic choriomeningitis virus clone 13 and tumor samples, we demonstrate that CD8 T cells differentiate into a less exhausted TCF1high and a more exhausted TCF1low population. Virus-specific TCF1high CD8 T cells, which resemble T follicular helper (TFH) cells, persist and recall better than do TCF1low cells and act as progenitor cells to replenish TCF1low cells. We show that TCF1 is both necessary and sufficient to support this progenitor-like CD8 subset, whereas cell-intrinsic type I interferon signaling suppresses their differentiation. Accordingly, cell-intrinsic TCF1 deficiency led to a loss of these progenitor CD8 T cells, sharp contraction of virus-specific T cells, and uncontrolled viremia. Mechanistically, TCF1 repressed several pro-exhaustion factors and induced Bcl6 in CD8 T cells, which promoted the progenitor fate. We propose that the TCF1-Bcl6 axis counteracts type I interferon to repress T cell exhaustion and maintain T cell stemness, which is critical for persistent antiviral CD8 T cell responses in chronic infection. These findings provide insight into the requirements for persistence of T cell immune responses in the face of exhaustion and suggest mechanisms by which effective T cell-mediated immunity may be enhanced during chronic infections and cancer.