key: cord-0767459-fbodme9o
authors: Hafler, David; Sumida, Tomokazu; Dulberg, Shai; Schupp, Jonas; Stillwell, Helen; Axisa, Pierre-Paul; Comi, Michela; Lincoln, Matthew; Unterman, Avraham; Kaminski, Naftali; Madi, Asaf; Kuchroo, Vijay
title: Type I Interferon Transcriptional Network Regulates Expression of Coinhibitory Receptors in Human T cells
date: 2021-06-08
journal: Res Sq
DOI: 10.21203/rs.3.rs-133494/v1
sha: 36bdb87acb5145ab05ad27c11e56d9e1650547ff
doc_id: 767459
cord_uid: fbodme9o

While inhibition of T cell co-inhibitory receptors has revolutionized cancer therapy, the mechanisms governing their expression on human T cells have not been elucidated. Type 1 interferon (IFN-I) modulates T cell immunity in viral infection, autoimmunity, and cancer, and may facilitate induction of T cell exhaustion in chronic viral infection. Here we show that IFN-I regulates co-inhibitory receptor expression on human T cells, inducing PD-1/TIM-3/LAG-3 while surprisingly inhibiting TIGIT expression. High-temporal-resolution mRNA profiling of IFN-I responses enabled the construction of dynamic transcriptional regulatory networks uncovering three temporal transcriptional waves. Perturbation of key transcription factors on human primary T cells revealed unique regulators that control expression of co-inhibitory receptors. We found that the dynamic IFN-I response in vitro closely mirrored T cell features with IFN-I linked acute SARS-CoV-2 infection in human, with high LAG3 and decreased TIGIT expression. Finally, our gene regulatory network identified SP140 as a key regulator for differential LAG3 and TIGIT expression, which were validated at the level of protein expression. The construction of IFN-I regulatory networks with identification of unique transcription factors controlling co-inhibitory receptor expression may provide targets for enhancement of immunotherapy in cancer, infectious diseases, and autoimmunity.

CD4 + and CD8 + T cells (Figure 1a, b, Supplementary Figure 1b) . Unexpectedly, both IL-27 and IFN-b suppressed the expression of TIGIT in CD4 + and CD8 + T cells (Supplementary Figure 1c) . We also observed the increased production of IL-10 by IFN-b; however, IL-10 induction by IL-27 was modest in our in vitro culture settings, which may re ect the difference between mouse and human T cell responses toward IL-27 stimulation (Supplementary Figure 1d) . To determine whether these observations stemmed from the effect of IFN-b on cellular proliferation, we performed a proliferation assay using cell trace violet dye. We found there was no differences in cellular division in naïve CD4 + T cells between control and the IFN-b condition. Additionally, there was even less proliferation of memory CD4 + T cells in the IFN-b condition, indicating that the induction of co-inhibitory receptors by IFN-b is not driven by a state of higher proliferation in T cells (Supplementary Figure 2) , consistent with previous studies 18, 19 . We further determined the impact of IFN-b treatment on gene expression kinetics for co-inhibitory receptors by qPCR. Gene expression dynamics for core co-inhibitory receptors (HAVCR2, LAG3, PDCD1) were upregulated by IFN-b for most time points. In contrast TIGIT was downregulated, which was con rmed by protein expression using ow cytometry (Figure 1c , d). We examined the expression of other co-inhibitory receptors, and found IFN-b induced co-expression of multiple co-inhibitory receptors (e.g. HAVCR2, PDCD1, LAG3), but inhibited expression of others (TIGIT, CD160, BTLA) ( Figure 1e ). Collectively, these data elucidate a role for IFN-b as a cytokine that can directly control multiple co-inhibitory receptors in both human CD4 + and CD8 + T cells in vitro.

To uncover the regulatory mechanisms underlying the IFN-b response in human primary T cells, we generated a transcriptional pro le at high temporal resolution. We used bulk mRNA-seq at ten time points along a 96-hour time course with and without IFN-b treatment (Supplementary Figure 3a) . To avoid interindividual variation, we selected one healthy subject whose T cells exhibited a stable response to IFN-b, and repeated the experiment three times at a two-week interval for each experiment. We identi ed 1,831 (for CD4 + T cells) and 1,571 (for CD8 + T cells) differentially expressed genes (DEGs) across time points with IFN-b treatment, revealing a temporal shift of gene expression patterns in both CD4 + and CD8 + T cells (Figure 2a, Supplementary Figure 3b Interestingly, OSM, which is reported to amplify IFN-b response and suppress Th17 differentiation 20 , was signi cantly induced by IFN-b from the early phase, and maintained induction in all time points (Supplementary Figure 3c) . Among DEGs, we identi ed dynamic expression of 134 TFs for CD4 + T cells and 100 TFs for CD8 + T cells, which were both up-and down-regulated over the course of differentiation (Figure 2c ).

To narrow the list of TFs for perturbation, we prioritized TFs that are differentially expressed in both CD4 + and CD8 + T cells. We then further selected TFs associated with: 1) human tumor in ltrating T cells (TILs) [21] [22] [23] [24] ; 2) HIV speci c T cell signature in progressive patients 25 ; and 3) IL-27 driven co-inhibitory regulators 3 (see Methods). We con rmed that these TFs are identi ed as interferon stimulated genes (ISGs) in human immune cells by the Interferome dataset 26 (Figure 3a ). In total, 31 TFs were listed as candidates based on the overlap between ISG, TIL and IL-27 signatures and we chose 19 of them for perturbation. Since TCF-1 (encoded by TCF7) and Blimp-1 (encoded by PRDM1) are known to express functionally distinct isoforms, we also targeted a unique sequence for the long isoforms (TCF7L and PRDM1L, respectively), resulting in perturbation of 21 different targets in total.

Considering the majority of these regulators are induced at the early (1-2h) and intermediate (4-16h) phases, it was important to perform gene deletion prior to T cell receptor activation. For human primary T cells gene knockdown, we adopted lentiviral delivery of shRNAs with lentiviral gene product X (Vpx) containing virus-like particles (VLPs) system to e ciently transduce lentivirus into unstimulated primary human naïve T cells 27, 28 (Figure 3b ). Spinoculation with Vpx-VLPs signi cantly increased the number of GFP expressing T cells (~30-60%) compared to normal LV particle transduction without Vpx-VLPs (~1-5%), and resulted in successful transduction of lentiviral vectors into non-blasting/quiescent cells (Supplementary Figure 4a) . We achieved e cient knockdown of at least 60% gene expression for 21 target TFs in human naïve CD4 + T cells (Figure 3c ).

To identify the effect of perturbation for each regulator, Principal Component Analysis (PCA) was applied to changes in RNA expression associated with each transcription factor knock down (Figure 3d). PC1   divided the impact of perturbation into two modules of regulators; BATF, MAF, ETS2, HOPX, SP140, BCL3,  ID3, and BATF3 constitute 'IFN-I regulator module 1', and IRF1, IRF2, IRF4, STAT1, STAT3, ARID5A, ARID5B,  TCF7, PRDM1, PRDM1L, KLF5, and TCF7L constitute a distinct 'IFN-I regulator module 2' . To visualize the contribution of the selected genes to the PCs and the directionality of the contribution, a PCA biplot analysis was adopted. We found that ISGs are divided into two groups; classical ISGs that are correlated with 'IFN-I regulator module 1' (depicted in green arrows in Figure 3e ), and ISGs that are correlated with 'IFN-I regulator module 2' (depicted in orange arrows in Figure 3e ), which is predominantly in PC1. These results suggest that ISGs are bi-directionally regulated by different modules of TFs (Supplementary To characterize the impact of differentially expressed transcription factors (DETFs) in response to IFN-b, we generated transcriptional regulatory networks describing TFs and their target genes for each of the transcriptional waves identi ed (Figure 4a , Methods). When comparing the three regulatory networks, early and late networks had similar numbers of TFs (46 and 42 TFs, respectively), while the intermediate network contained 73 TFs (Figure 4b, top) . Interestingly, the ratio between up and down-regulated TFs differs between the three regulatory waves. The early and intermediate network contained more upregulated TFs than down-regulated TFs; in contrast, the late network had more down-regulated than upregulated TFs. Thus, IFN-I induced differentiation involves dominance of up-regulated TFs in the rst 16 hours, replaced by the dominance of down-regulated TFs after 48 hours. We next ranked the TFs based on the enrichment of their target genes and their centrality in the networks (Methods), highlighting the signi cance of each TF to the network (Figure 4b , middle). In the early regulatory network, MYC and CDKN1B/KDM5B were among the most dominant up and down-regulated TFs, respectively, in this transcriptional wave. These data indicate that T cell metabolic activation, cell cycle regulation, and transcriptional activation are promoted by IFN-b treatment 29, 30 . Interestingly, FLI1, which is novel as a IFN-I downstream TF and was recently reported to control effector response in T cells 31 , also dominantly regulated the early transcriptional wave. In the intermediate regulatory network, MYC, MAF, IRF1, AFF1, ATF3, and TBX21 were among the most dominant up-regulated TFs for this transcriptional wave. In the late network, effector function related regulators are upregulated (PRDM1, RUNX2, MAF, BCL3); in contrast, the TFs associated with Treg differentiation and maintenance (STAT5A, FOXP1, MYB) were down-regulated, suggesting the skewed differentiation toward effector-like signature.

To further study the relationships between the DETFs, we generated hierarchical backbone networks in order to represent their relationships (Figure 4b, bottom) . Interestingly, the top TFs in all transcriptional time waves were down-regulated in response to IFN-b, while TFs lower in the network hierarchy were more up-regulated. A few examples from the early transcriptional wave hierarchical backbone include CDKN1B and MAZ, which appeared at the top of the hierarchy, whereas KLF5, MYC, and FLI1 were lower in the hierarchy. It was of interest that these TFs were also highly dominant in the regulatory network. These data suggest that loss of suppression of TFs at higher hierarchy triggers the activation of downstream effector TFs under IFN-I response, which was also observed in the intermediate and late regulatory network. The elucidation of this backbone network enables us to shed light on the regulatory interactions within each component of the transcriptional network, providing further depth to the extent of interactions within the network.

While T cell differentiation under IFN-b is characterized by three major transcriptional waves, we hypothesized that there are key TFs that bridge each wave to the next. To this end, we speci cally identi ed TFs that participate in more than one of these transcriptional waves, and termed these 'Bridging TFs' (Figure 4c ). Examples of dominant 'Bridging TFs' between early and intermediate waves include KLF5 and STAT2. Examples of intermediate to late waves include MAF, PRDM1, and MYB. Finally, there are TFs that were upregulated throughout the entire differentiation; such as STAT1, HIF1A, and TBX21. Generally, 'Bridging TFs' tend to be more dominant than other TFs; thus, it is possible that 'Bridging TFs' play an important role in the transition between different transcriptional waves. Indeed, our perturbation experiment demonstrated the critical roles of those 'Bridging TFs' in the regulation of ISGs and coinhibitory receptors (Figure 4c ). Our computational analysis revealed the temporal dynamics of complex regulatory interactions during the IFN-I response and highlighted the usefulness of our approach in discovering this new aspect of IFN-I induced transcriptional regulation.

In vivo validation of regulatory modules controlling IFN-I/co-inhibitory receptors axis in human T cells As it was important to provide direct in vivo evidence for the role of IFN-I on T cell co-inhibitory receptor expression, we sought to validate our regulatory network in the human setting where the IFN-I response of T cells is induced acutely. As acute viral infections are strongly associated with IFN-I responses, we examined a number of clinical models where viral infection is closely linked to IFN-I T cell response. We found that our analysis of single cell RNA seq (scRNA-seq) analysis of T cells in COVID-19 patients revealed an extremely high correlation between viral load and IFN-I score (r=0.8) and time difference between paired samples and the respective change in IFN-I score (r=0.97) 32 , providing a unique opportunity to generate a rich dataset to determine whether the in vitro T cell response to IFN-I can be validated during an acute viral human infection strongly associated with a IFN-I signal.

By using our scRNA-seq data, we subclustered T cell populations into 13 subpopulations and identi ed ve CD4 + T cell and ve CD8 + T cell subsets (Figure 5a, Supplementary Figure 5a , b). We rst focused on total CD4 + T cells and CD8 + T cells and con rmed that the IFN-I response signature is higher in progressive patients who required admission to the ICU and eventually succumbed to the disease ( Figure  5b ). Expression of co-inhibitory receptors differed across disease conditions, but the trend was conserved between CD4 + and CD8 + T cells. We observed a strikingly similar pattern of co-inhibitory receptor expression with IFN-I stimulation in vitro and in vivo. Indeed, we observed the upregulation of 'IFN-I up coinhibitory receptors' (LAG3/HAVCR2) and the downregulation of 'IFN-I down co-inhibitory receptors' (TIGIT/LAIR1/SLAMF6) in T cells from COVID-19 patients (Figure 5c , d). As expected, 'IFN-I up coinhibitory receptors' were positively correlated with canonical ISGs expression, but 'IFN-I down coinhibitory receptors' were not, suggesting that there are different regulatory mechanisms dictating coinhibitory receptor expression patterns ( Figure 5e ). Next, we investigated which subpopulation of CD4 + and CD8 + T cells is more affected by the IFN-I response, and computed the IFN-I score across subpopulations for each T cell subtype in COVID-19 patients (Figure 5f ). Within the subpopulations that exhibited higher IFN-I scores, dividing CD4 + /CD8 + T cells and ISG + CD8 + T cells were uniquely increased in COVID-19 patients, particularly in severe patients 32, 36 . Moreover, these subpopulations and effector T cells expressed higher level of co-inhibitory receptors and 'IFN-I regulator module-1' compared to the other subpopulations ( Figure 5g, Supplementary Figure 5c ).

We then examined which subpopulations were more enriched in the three transcriptional waves of IFN-I response. DEGs speci c for each wave were used to compute the scores for the CD4 + and CD8 + T cell subpopulations in COVID-19 patients (Table 1 , Methods). We found that T cells induced in vitro with IFN-I strongly mirrored the intermediate wave score on dividing CD4 + and CD8 + T cells, and the late wave score on effector CD4 + T cells and ISG + CD8 + /effector CD8 + T cells (Figure 5h Figure 5h ). These ndings strongly suggest that in vitro regulatory network can be utilized as a strong tool to explore human acute viral response in vivo.

Here, our systematic, computational and biological approach identi es IFN-I as a major driver of coinhibitory receptor regulation in human T cells. While classical ISG induction has been extensively studied, those investigations have focused primarily on the canonical JAK-STAT pathway downstream of IFN-I receptor. Given that IFN-I exhibits multiple functions in context-dependent roles, a more complex understanding of the IFN-I response beyond this canonical pathway with a more extensive analysis of ISG transcriptional regulation in T cells is critical for elucidating the mechanism of co-inhibitory receptor regulation.

In these studies, we build a dynamic gene regulatory network that controls IFN-I response, and identi ed key regulatory modules of ISG transcription in T cell responses to IFN-b. Our approach unveiled two mutually antagonistic modules of ISG regulators, which, when acting concordantly, may explain how the harmonized IFN-I induced T cell response is achieved. Within the two modules, we highlighted SP140 as a potential regulator that controls LAG3 and TIGIT in an opposing manner, and STAT3 as a unique positive regulator for TIM-3. These ndings provide novel insight into the landscape of the ISG transcriptional network, and sheds light on the large contribution of the noncanonical IFN-I pathway during IFN-I response in T cells [34] [35] [36] . Although the newly identi ed regulators (e.g. SP140, BCL3) in this study are not necessarily directly downstream of the conventional JAK/STAT pathway and may act differently depending on the context, they are nevertheless attractive targets for manipulation of speci c downstream functional molecules such as co-inhibitory receptors in T cells.

We demonstrate the relevance of our in vitro T cell IFN-I response by integrating scRNA-seq from COVID-19 patients, where a predominant T cell IFN-I response was observed. Intriguingly, the expression pattern of co-inhibitory receptors on T cells in vitro are highly replicated in severe COVID-19 cases, and classical ISGs were well correlated with one module of co-inhibitory receptors (LAG3/PDCD1/HAVCR2), but not with the other modules (TIGIT/CD160/BTLA/LAIR1). While dynamics of IFN-I on T cells from COVID-19 patients should be taken into account with higher temporal pro ling, we con rmed that the IFN-I response was clearly reduced at a later time point; thus, our data based on earlier collection of blood should re ect the active ISG transcriptomics during human acute viral response. Given the IFN-I response has been shown to contribute to chronic viral infection and cancer, the novel regulators we identi ed can be examined in these diseases. Further investigations of factors that characterize acute viral infections is likely to differ from chronic viral infections and will be of interest to explore in the context of long-term human infections not well modeled by COVID-19.

In conclusion, our systems biology approach identi es the cytokine signals and regulatory mechanisms that drive expression of co-inhibitory receptors in humans, and provides a pathway to comprehensively capture the dynamics of their expression in humans. Our results will also advance the understanding of the host immune response to a variety of viral infections, and could serve as a resource for mining of existing datasets. Uncovering novel ISG regulators controlling co-inhibitory receptors will create a foundation for further development of new therapeutics for a multitude of different malignant and infectious diseases.

Peripheral blood was drawn from healthy controls who were recruited as part of an Institutional Review Board (IRB)-approved study at Yale University, and written consent was obtained. All experiments conformed to the principles set out in the WMA Declaration of Helsinki and the Department of Health and Human Services Belmont Report.

Peripheral blood mononuclear cells (PBMCs) were prepared from whole blood by Ficoll gradient centrifugation (Lymphoprep, STEMCELL Technologies) and used directly for total T cell enrichment by negative magnetic selection using Easysep magnetic separation kits (STEMCELL Technologies). Cell suspension was stained with anti-CD4 (RPA-T4), anti-CD8 (RPA-T8), anti-CD25 (clone 2A3), anti-CD45RO (UCHL1), anti-CD45RA (HI100) and anti-CD127 (hIL-7R-M21, all from BD Biosciences) for 30 minutes at 4°C. Naïve CD4 + T cells (CD4 + /CD25 neg /CD127 + /CD45RO neg /CD45RA + ) and naïve CD8 + T cells (CD8 ++ /CD25 neg /CD127 + /CD45RO neg /CD45RA + ) were sorted on a FACSAria (BD Biosciences). Sorted cells were plated in 96-well round-bottom plates (Corning) and cultured in RPMI 1640 medium supplemented with 5 % Human serum, 2 nM L-glutamine, 5 mM HEPES, and 100 U/ml penicillin, 100 mg/ml streptomycin, 0.5 mM sodium pyruvate, 0.05 mM nonessential amino acids, and 5% human AB serum (Gemini Bio-Products). Cells were seeded (30,000-50,000/wells) into wells pre-coated with antihuman CD3 (2 mg/ml, clone UCHT1, BD Biosciences) along with soluble anti-human CD28 (1mg/ml, clone 28.2, BD Biosciences) in the presence or absence of human IFN-b (500 U/ml: Pestka Biomedical Laboratories) or IL-27 (100 ng/ml: BioLegend) without adding IL-2.

Lentiviral and Vpx-VPLs production Lentiviral plasmids encoding shRNA were obtained from Sigma-Aldrich. Each plasmid was transformed into One Shot Stbl3 chemically competent cells (Invitrogen) and puri ed by ZymoPURE plasmid Maxiprep kit (Zymo research). Lentiviral pseudoparticles were obtained after plasmid transfection of 293FT cells using Lipofectamine 2000 (Invitrogen) or TurboFectin 8.0 Transfection Reagent (Origene). To prepare Vpx-VLPs, 293T cells were co-transfected by Lipofectamine 2000 or TurboFectin 8.0 Transfection Reagent with the 5 mg pMDL-X, 2.5 mg pcRSV-Rev, 3.5 mg X4-tropic HIV Env, and 1 mg pcVpx/myc, as described previously with some modi cations 37, 38 . The medium was replaced after 6-12 h with fresh media with 1X Viral boost (Alstem). The lentivirus or Vpx-VLPs containing media was harvested 72 h after transfection and concentrated 80 times using Lenti-X concentrator (Takara Clontech) or Lenti Concentrator (Origene). LV particles were then resuspended in RPMI 1640 media without serum and stored at -80°C before use. Virus titer was determined by using Jurkat T cells and Lenti-X GoStix Plus (Takara Clontech).

Two step Vpx-VLP and LV transduction was performed as described previously with some modiciations 36 . Vpx are pseudotyped with X4-tropic HIV Env to promote e cient entry of Vpx-VLPs into quiescent human T cells 37 RNA-seq library preparation and data analysis cDNAs were generated from isolated RNAs using SMART-Seq v4 Ultra Low Input RNA Kit for sequencing (Takara/Clontech). Barcoded libraries were generated by the Nextera XT DNA Library Preparation kit (Illumina) and sequenced with a 2x100 bp paired-end protocol on the HiSeq 4000 Sequencing System (Illumina).

After sequencing, adapter sequences and poor-quality bases (quality score < 3) were trimmed with Trimmomatic. Remaining bases were trimmed if their average quality score in a 4 bp sliding window fell below 5. FastQC was used to obtain quality control metrics before and after trimming. Remaining reads were aligned to the GRCh38 human genome with STAR 2.5.2 39 . We used Picard to remove optical duplicates and to compile alignment summary statistics and RNA-seq summary statistics. After alignment, reads were quantitated to gene level with RSEM 40 using the Ensembl annotation.

The correlation matrix is created by Pearson correlating 41 the IFN-b expression pro le of each time point with all the other time points, creating a symmetric matrix of Pearson correlation coe cients.

The differential expression (DE) 42 of the sequenced genes from every time point was calculated. The DE was calculated using the DEseq2 43 R package. An in-house decision algorithm was built to determine which genes are DE. The algorithm used three separate testing methods available in DEseq2 for calculating DE genes: Wald 44 , likelihood ratio test (LRT) 45 , and time-course 46 . For each of the three methods, genes with false discovery rate (FDR) adjusted P.value bellow 0.05, are regarded as DE. The algorithm de ned genes as DE differently for CD4 + and CD8 + T cells. For CD4 + T cells, if TFs appeared in two out of the three calculating methods (agree by two), they were regarded as DE. For CD8 + T cells, if TFs appeared in any of the methods above, they were regarded as DE.

The list of TFs for perturbation was selected based on following aspects: 1) overlapped differentially expressed TFs across the time point between CD4 + and CD8 + T cells in our results. Intersection of DETFs in our in vitro data were chosen; 2) differentially expressed TFs in human tumor in ltrated T cells 21-24 that were signi cantly correlated with exhausted T cell cluster where LAG-3/PD-1/TIM-3 were highly upregulated. 'Human TIL score' for each gene was calculated by the number of times it was shared between the four different human cancer TIL datasets [21] [22] [23] [24] ; 3) HIV speci c T cell signature in progressive patients compared to stable patients 25 . 4) TFs that were included by IL-27 and categorized as IL-27 driven co-inhibitory receptor modules 3 . 'Human ISG score' for each gene was calculated by the number of times it was shared between the three different categories (T cells, PBMCs, and all immune cells) of human ISGs identi ed by Interferome database. All perturbed TFs were con rmed as IFN-I responsible genes that showed 'human ISG score' more than 1.

Heatmap of perturbed TFs 21 TFs were perturbed using lentiviral shRNA, together with a scramble shRNA control (SCR). In order to better understand the effect of perturbing said TFs, selected genes of interest (GOIs) were analyzed. The DE analysis was conducted for perturbation against control, genes who yielded n FDR adjusted P.value lower than 0.05 were regarded as signi cant and display a white plus on their tile.

The PCA was conducted on a DEGs de ned as above as variables, and perturbed genes as observations. The data was normalized by the scramble shRNA control identically to the Heatmap of perturbed TFs.

Although PDCD1 was not de ned as DE genes across the time course, it was clearly differentially expressed at later time points in mRNA-seq and qPCR, thus we manually depicted in Figure 5f . Genes of IFN Score B 47 (see Table 1 ) were represented as ISGs in Figure 3e . The PCA and biplot analysis were calculated and visualized using the R package FactoMineR 48 .

Following DE analysis in each transcriptional wave, the DE genes were separated to TFs and their targets (from->to). The targets of all DETFs were determined using ChIP-seq data from the database GTRD 49 .

TFs and targets de ned as DE were added as network nodes, and edges (connections) were added between them. The network gures were created using the software Cytoscape 50 .

We ranked the DE TFs of each transcriptional time wave to identify which are the most dominant in the overall differentiation process. HG stands for hyper-geometric, the value in the heatmap is the of a hypergeometric enrichment test. The targets of each TF are tested for enrichment of DE targets in the network, relatively to targets that aren't DE in the network. The HG calculation was conducted using the python SciPy package 51 . Cent stands for centrality, which is a parameter that is given to each node, based on the shortest path from the node to the other nodes in the network. It represents how central and connected a node is in the rest of the network 52 . The centrality calculation was conducted using the python NetworkX 53 package. The rank column is an average of both HG and Cent values, after normalization.

Integration of Perturbation Data to Regulatory Networks DE analysis was conducted for perturbation against control. Genes that were signi cantly affected by a TF perturbation were added as a "validated" edges between the perturbed TF and the target gene. If a gene was up-regulated by a TF perturbation, the interaction between them is registered as downregulation. If a gene was down-regulated by a TF perturbation, the interaction between them is registered as up-regulation.

Using the software Cytoscape, we implemented a hierarchical layout, which takes into account the directionality of the connections between the TFs. A TF which has only outgoing connections will be placed at the top of the hierarchy, while a TF which has only incoming connections will be placed at the bottom.

DETFs and their targets from the three transcriptional waves were combined to create a comprehensive network of the dynamic between transcriptional waves. TFs and their targets were annotated by the transcriptional wave in which they are DE. TFs that appear in more than one transcriptional wave are regarded as bridging TFs.

Reanalysis of COVID-19 single-cell RNA sequencing data A PBMC single cell RNA seq data set of 10 COVID-19 patients and 13 matched controls was reanalyzed which had been previously performed and reported by us 32 . We have described the full cohort and detailed methods elsewhere 32 Libraries had been sequenced on an Illumina Novaseq 6000 platform. Raw reads had been demultiplexed and processed using Cell Ranger (v3.1) mapping to the GRCh38 (Ensembl 93) reference genome.

Resulting gene-cell matrices had been analyzed using the package Seurat 54 in the software R 55 including integration of data, clustering, multiplet identi cation and cell type annotation. The nal annotated R object was used and re-analyzed in Seurat with default settings -unless otherwise speci ed -as follows:

The three cell populations "Dividing T & NK", "Effector T" and "Memory CD4 & MAIT" were each subsetted and reclustered to obtain a ner cell type granularity as they included a mix of CD4, CD8, MAIT and gamma delta T cells. Per subset, the top 500 variable genes were determined by the "FindVariableFeatures" function using the "vst" method. Data was scaled using the "ScaleData" function regressing out the total number of UMI and the percentage of UMIs arising from the mitochondrial genome. After Principal Component (PC) Analysis, the rst 10 Principal Components (PCs) were utilized to detect the nearest neighbors using the "FindNeighbors" function and clustered by Seurat's Louvain algorithm implementation "FindClusters" using a resolution of 0.2 for "Dividing T & NK", of 0.3 for "Effector T" and of 0.1 for "Memory CD4 & MAIT" subsets. Cluster-speci c gene expression pro les were established using the "FindAllMarkers" per cluster and per subset to annotate the clusters. New cell type annotations were then transferred back to the full dataset.

A new Uniform Approximation and Projection (UMAP) embedding was created by integrating the datasets on a subject level as follows: A subset containing all T cells was generated, which was then split by subject. For each subject, the top 2000 variable genes were selected, then integration anchors determined by "FindIntegrationAnchors" (with k. lter = 150). These anchors were used to integrate the data using the "IntegrateData" function with top 30 dimensions. The integrated data was scaled, subjected to a PC analysis and the top 13 PCs used as input for the "RunUMAP" function on 75 nearest neighbors 56 .

Module scores were calculated using the "AddModuleScore" function using a) all genes within the GO list "RESPONSE TO TYPE I INTERFERON" (GO:0034340) 57 and b) all genes signi cantly associated with either of the three waves in our in vitro perturbation experiments (see Table 1 ). Differential gene expression was established using Seurat's implementation of the Wilcoxon Rank Sum test within the "FindMarkers" function with a Bonferroni correction for multiple testing.

Detailed information about statistical analysis, including tests and values used, is provided in the gure legends. P-values of 0.05 or less were considered signi cant.

Data and software availability

The sequence data generated in this study will be deposited in the Gene Expression Omnibus (GEO) and the accession code will be provided prior to publication. Transcriptional regulatory network under IFN-I response a, Overview of regulatory network generation. A scheme of the pipeline, in order to generate preliminary regulatory network is generated from integrating the gene expression kinetics data coupled with TF-target gene datasets. The key regulators' perturbation data was further integrated to re ne the preliminary network. b, In depth view of the transcriptional regulation at each wave. Top row; the representation of regulatory networks highlighting TFs interaction. Integration of IFN-I regulatory network with T cells signature in COVID-19 a, UMAP representation of T cells from healthy control samples (n = 13) and COVID-19 samples (n = 18). 13 subcluster were identi ed. b, IFN-I score for CD4+ and CD8+ T cells across the three disease conditions. c, Heatmaps for coinhibitory receptors expression in CD4+ and CD8+ T cells across the three disease conditions. d, Expression of key co-inhibitory receptors between control vs COVID-19 for CD4+ and CD8+ T cells.

Timing and magnitude of type I interferon responses by distinct sensors impact CD8 T cell exhaustion and chronic viral infection

Cell Exhaustion During Chronic Viral Infection and Cancer

Induction and transcriptional regulation of the co-inhibitory gene module in T cells

IL-27 and TCR Stimulation Promote T Cell Expression of Multiple Inhibitory Receptors

Immunomodulatory functions of type I interferons

Type I interferons in host defence and in ammatory diseases

Type-I-interferons in infection and cancer: Unanticipated dynamics with therapeutic implications

Type 1 interferons and antiviral CD8 T-cell responses

Type I interferons: bene cial in Th1 and detrimental in Th17 autoimmunity

Persistent LCMV infection is controlled by blockade of type I interferon signaling

Blockade of chronic type I interferon signaling to control persistent LCMV infection

Coinhibitory Pathways in Immunotherapy for Cancer

Lag-3, Tim-3, and TIGIT: Co-inhibitory Receptors with Specialized Functions in Immune Regulation

Tim-3 and its role in regulating anti-tumor immunity

Molecular and transcriptional basis of CD4(+) T cell dysfunction during chronic infection

Targeting type I interferon-mediated activation restores immune function in chronic HIV infection

IL-27 synthesis induced by TLR ligation critically depends on IFN regulatory factor

Antiproliferative action of interferon-alpha requires components of T-cellreceptor signalling

Virus-induced transient immune suppression and the inhibition of T cell proliferation by type I interferon

Oncostatin M Suppresses Activation of IL-17/Th17 via SOCS3 Regulation in CD4+ T Cells

Molecular Pathways of Colon In ammation Induced by Cancer Immunotherapy

De ning T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma

Intratumoral CD4(+) T Cells Mediate Anti-tumor Cytotoxicity in Human Bladder Cancer

Dysfunctional CD8 T Cells Form a Proliferative, Dynamically Regulated Compartment within Human Melanoma

Transcriptional analysis of HIV-speci c CD8+ T cells shows that PD-1 inhibits T cell function by upregulating BATF

Interferome v2.0: an updated database of annotated interferon-regulated genes

SAMHD1 is the dendritic-and myeloid-cell-speci c HIV-1 restriction factor counteracted by Vpx

SAMHD1 restricts HIV-1 infection in resting CD4(+) T cells

c-Myc represses FOXO3a-mediated transcription of the gene encoding the p27(Kip1) cyclin dependent kinase inhibitor

RSV-Induced H3K4 Demethylase KDM5B Leads to Regulation of Dendritic Cell-Derived Innate Cytokines and Exacerbates Pathogenesis In Vivo

<em>In vivo</em> CRISPR screening identi es Fli1 as a transcriptional safeguard that restrains effector CD8 T cell differentiation during infection and cancer. bioRxiv

Single-Cell Omics Reveals Dyssynchrony of the Innate and Adaptive Immune System in Progressive COVID-19. medRxiv

LAG-3 regulates CD8+ T cell accumulation and effector function in murine selfand tumor-tolerance systems

Biologic consequences of Stat1-independent IFN signaling

Transcriptional Regulation of Antiviral Interferon-Stimulated Genes

Parsing the Interferon Transcriptional Network and Its Disease Associations

E cient delivery of lentiviral vectors into resting human CD4 T cells

Vpx-containing dendritic cell vaccine induces CTLs and reactivates latent HIV-1 in vitro

STAR: ultrafast universal RNA-seq aligner

RSEM: accurate transcript quanti cation from RNA-Seq data with or without a reference genome

Analytic posteriors for Pearson's correlation coe cient

Differentiation of Spermatogonia Stem Cells into Functional Mature Neurons Characterized with Differential Gene Expression

Moderated estimation of fold change and dispersion for RNAseq data with DESeq2

Wald tests for variance-adjusted equivalence assessment with normal endpoints

Determining qualitative effect size ratings using a likelihood ratio scatter matrix in diagnostic test accuracy systematic reviews

A global non-coding RNA system modulates ssion yeast protein levels in response to stress

A novel two-score system for interferon status segregates autoimmune diseases and correlates with clinical features

FactoMineR: An R Package for Multivariate Analysis

GTRD: a database on gene transcription regulation-2019 update

Cytoscape: a software environment for integrated models of biomolecular interaction networks

SciPy 1.0: fundamental algorithms for scienti c computing in Python

Percolation centrality: quantifying graph-theoretic impact of nodes during percolation in networks

Exploring Network Structure, Dynamics, and Function Using NetworkX

Comprehensive Integration of Single-Cell Data

UMAP: Uniform Manifold Approximation and Projection for Dimension Reduction

Gene ontology: tool for the uni cation of biology. The Gene Ontology Consortium

Figures Average expression per subject for each gene is shown. *p < 0.05, **p < 0.01, ***p < 0.001. Kruskal-Wallis test. e, Correlation matrix of ISGs (dark gray) and co-inhibitory receptors (light gray

Computed three transcriptional waves (early, intermediate, and late) score for the subsets of CD4+ and CD8+ T cells in COVID-19 patients. Scores were calculated based on upregulated DEGs of CD4+ and CD8+ T cells for each transcriptional wave. i, Regulatory relationship between regulators in intermediate phase network for LAG3 and TIGIT are shown. Positive regulation (TF to target) is highlighted in red and negative regulations in blue. j, Box plots showing expression of key regulators between control vs COVID-19 for CD4+ T cells. Average expression per subject for each gene is shown

We thank L. Devine and C. Wang for assistance with FACS based cell sorting; G. Wang 

Perturbation of key transcription factors in quiescent human T cells a, Characterization of candidate TFs for perturbation. Perturbed TFs are listed based on overlap between differentially expressed TFs of CD4+ T cells and CD8+ T cells. Human ISG score (top; blue), human TIL co-inhibitory receptors score (green),HIV speci c T cell signature genes in progressive patients (yellow), and IL-27 driven co-inhibitory receptor regulators (orange) are shown for each TFs. b, Experimental work ow of Vpx-VLP supported lentiviral shRNA perturbation. Ex vivo isolated naïve CD4 T cells were transduced with Vpx-VLPs, followed by two times of lentiviral particle transduction before starting T cell activation. T cells were stimulated with anti-CD3/CD28 in the absence or presence of IFN-b (500 U/ml) for 96h and GFP positive cells were sorted by The thicker and darker an edge is the more TF-target connections it represents. Target genes are represented by up and down hexagons, according to their regulatory response to IFN-b. Middle row;heatmaps representing a ranking of the TFs based on 'Cent' stands for centrality and 'HG' stands for hypergeometric test. Bottom row; hierarchical backbone networks. Red circles represent up-regulated TFs, blue circles represent down-regulated TFs. c, Dynamics of TFs regulation across the transcriptional waves. Each hexagon represents targets from each transcriptional wave. Green circles represent regulatory TFs which are differentially expressed only in one transcriptional wave they are connected to, while purple circles represent bridging TFs, which are DE in all transcriptional waves they are connected to. The thicker and darker an edge is, the more TF-target connections it represents.