SUMMARY While many genetic variants have been associated with risk for human diseases, how these variants affect gene expression in various cell types remains largely unknown. To address this gap, the DICE (Database of Immune Cell Expression, Expression quantitative trait loci (eQTLs) and Epigenomics) project was established. Considering all human immune cell types and conditions studied, we identified cis-eQTLs for a total of 12,254 unique genes, which represent 61% of all protein-coding genes expressed in these cell types. Strikingly, a large fraction (41%) of these genes showed a strong cis-association with genotype only in a single cell type. We also found that biological sex is associated with major differences in immune cell gene expression in a highly cell-specific manner. These datasets will help reveal the effects of disease risk-associated genetic polymorphisms on specific immune cell types, providing mechanistic insights into how they might influence pathogenesis (http://dice-database.org).

The contribution of CD4 + T cells to protective or pathogenic immune responses to SARS-CoV-2 infection remains unknown. Here, we present single-cell transcriptomic analysis of >100,000 viral antigen-reactive CD4 + T cells from 40 COVID-19 patients. In hospitalized patients compared to non-hospitalized patients, we found increased proportions of cytotoxic follicular helper (T FH ) cells and cytotoxic T helper cells (CD4-CTLs) responding to SARS-CoV-2, and reduced proportion of SARS-CoV-2-reactive regulatory T cells (T REG ). Importantly, in hospitalized COVID-19 patients, a strong cytotoxic T FH response was observed early in the illness which correlated negatively with antibody levels to SARS-CoV-2 spike protein. Polyfunctional T helper (T H )1 and T H 17 cell subsets were underrepresented in the repertoire of SARS-CoV-2-reactive CD4 + T cells compared to influenza-reactive CD4 + T cells. Together, our analyses provide insights into the gene expression patterns of SARS-CoV-2-reactive CD4 + T cells in distinct disease severities.

Therapies that boost the anti-tumor responses of cytotoxic T lymphocytes (CTLs) have shown promise; however, clinical responses to the immunotherapeutic agents currently available vary considerably, and the molecular basis of this is unclear. We performed transcriptomic profiling of tumor-infiltrating CTLs from treatment-naive patients with lung cancer to define the molecular features associated with the robustness of anti-tumor immune responses. We observed considerable heterogeneity in the expression of molecules associated with activation of the T cell antigen receptor (TCR) and of immunological-checkpoint molecules such as 4-1BB, PD-1 and TIM-3. Tumors with a high density of CTLs showed enrichment for transcripts linked to tissue-resident memory cells (TRM cells), such as CD103, and CTLs from CD103hi tumors displayed features of enhanced cytotoxicity. A greater density of TRM cells in tumors was predictive of a better survival outcome in lung cancer, and this effect was independent of that conferred by CTL density. Here we define the ‘molecular fingerprint’ of tumor-infiltrating CTLs and identify potentially new targets for immunotherapy.

TET2 and TET3 redundantly regulate Foxp3 stability, and their activity can be modulated by vitamin C.

Natural killer T cells (NKT cells) have stimulatory or inhibitory effects on the immune response that can be attributed in part to the existence of functional subsets of NKT cells. These subsets have been characterized only on the basis of the differential expression of a few transcription factors and cell-surface molecules. Here we have analyzed purified populations of thymic NKT cell subsets at both the transcriptomic level and epigenomic level and by single-cell RNA sequencing. Our data indicated that despite their similar antigen specificity, the functional NKT cell subsets were highly divergent populations with many gene-expression and epigenetic differences. Therefore, the thymus ‘imprints’ distinct gene programs on subsets of innate-like NKT cells that probably impart differences in proliferative capacity, homing, and effector functions.

CD4+ cytotoxic T lymphocytes (CD4-CTLs) have been reported to play a protective role in several viral infections. However, little is known in humans about the biology of CD4-CTL generation, their functional properties, and heterogeneity, especially in relation to other well-described CD4+ memory T cell subsets. We performed single-cell RNA sequencing in more than 9000 cells to unravel CD4-CTL heterogeneity, transcriptional profile, and clonality in humans. Single-cell differential gene expression analysis revealed a spectrum of known transcripts, including several linked to cytotoxic and costimulatory function that are expressed at higher levels in the TEMRA (effector memory T cells expressing CD45RA) subset, which is highly enriched for CD4-CTLs, compared with CD4+ T cells in the central memory (TCM) and effector memory (TEM) subsets. Simultaneous T cell antigen receptor (TCR) analysis in single cells and bulk subsets revealed that CD4-TEMRA cells show marked clonal expansion compared with TCM and TEM cells and that most of CD4-TEMRA were dengue virus (DENV)–specific in donors with previous DENV infection. The profile of CD4-TEMRA was highly heterogeneous across donors, with four distinct clusters identified by the single-cell analysis. We identified distinct clusters of CD4-CTL effector and precursor cells in the TEMRA subset; the precursor cells shared TCR clonotypes with CD4-CTL effectors and were distinguished by high expression of the interleukin-7 receptor. Our identification of a CD4-CTL precursor population may allow further investigation of how CD4-CTLs arise in humans and, thus, could provide insights into the mechanisms that may be used to generate durable and effective CD4-CTL immunity.

The expression of CD45RA is generally associated with naive T cells. However, a subset of effector memory T cells re-expresses CD45RA (termed TEMRA) after antigenic stimulation with unknown molecular characteristics and functions. CD4 TEMRA cells have been implicated in protective immunity against pathogens such as dengue virus (DENV). Here we show that not only the frequency but also the phenotype of CD4 TEMRA cells are heterogeneous between individuals. These cells can be subdivided into two major subsets based on the expression of the adhesion G protein-coupled receptor GPR56, and GPR56 + TEMRA cells display a transcriptional and proteomic program with cytotoxic features that is distinct from effector memory T cells. Moreover, GPR56 + TEMRA cells have higher levels of clonal expansion and contain the majority of virus-specific TEMRA cells. Overall, this study reveals the heterogeneity of CD4 TEMRA cells and provides insights into T-cell responses against DENV and other viral pathogens.

The molecular properties of CD8+ T cells that respond to SARS-CoV-2 infection are not fully known. Here, we report on the single-cell transcriptomes of >80,000 virus-reactive CD8+ T cells, obtained using a modified Antigen-Reactive T cell Enrichment (ARTE) assay, from 39 COVID-19 patients and 10 healthy subjects. COVID-19 patients segregated into two groups based on whether the dominant CD8+ T cell response to SARS-CoV-2 was ‘exhausted’ or not. SARS-CoV-2-reactive cells in the exhausted subset were increased in frequency and displayed lesser cytotoxicity and inflammatory features in COVID-19 patients with mild compared to severe illness. In contrast, SARS-CoV-2-reactive cells in the dominant non-exhausted subset from patients with severe disease showed enrichment of transcripts linked to co-stimulation, pro-survival NF-κB signaling, and anti-apoptotic pathways, suggesting the generation of robust CD8+ T cell memory responses in patients with severe COVID-19 illness. CD8+ T cells reactive to influenza and respiratory syncytial virus from healthy subjects displayed polyfunctional features and enhanced glycolysis. Cells with such features were largely absent in SARS-CoV-2-reactive cells from both COVID-19 patients and healthy controls non-exposed to SARS-CoV-2. Overall, our single-cell analysis revealed substantial diversity in the nature of CD8+ T cells responding to SARS-CoV-2.