An IL-4 signalling axis in bone marrow drives pro-tumorigenic myelopoiesis

Myeloid cells are known to suppress antitumour immunity1. However, the molecular drivers of immunosuppressive myeloid cell states are not well defined. Here we used single-cell RNA sequencing of human and mouse non-small cell lung cancer (NSCLC) lesions, and found that in both species the type 2 cytokine interleukin-4 (IL-4) was predicted to be the primary driver of the tumour-infiltrating monocyte-derived macrophage phenotype. Using a panel of conditional knockout mice, we found that only deletion of the IL-4 receptor IL-4Rα in early myeloid progenitors in bone marrow reduced tumour burden, whereas deletion of IL-4Rα in downstream mature myeloid cells had no effect. Mechanistically, IL-4 derived from bone marrow basophils and eosinophils acted on granulocyte-monocyte progenitors to transcriptionally programme the development of immunosuppressive tumour-promoting myeloid cells. Consequentially, depletion of basophils profoundly reduced tumour burden and normalized myelopoiesis. We subsequently initiated a clinical trial of the IL-4Rα blocking antibody dupilumab2-5 given in conjunction with PD-1/PD-L1 checkpoint blockade in patients with relapsed or refractory NSCLC who had progressed on PD-1/PD-L1 blockade alone (ClinicalTrials.gov identifier NCT05013450 ). Dupilumab supplementation reduced circulating monocytes, expanded tumour-infiltrating CD8 T cells, and in one out of six patients, drove a near-complete clinical response two months after treatment. Our study defines a central role for IL-4 in controlling immunosuppressive myelopoiesis in cancer, identifies a novel combination therapy for immune checkpoint blockade in humans, and highlights cancer as a systemic malady that requires therapeutic strategies beyond the primary disease site.

Circulating senescent myeloid cells infiltrate the brain and cause neurodegeneration in histiocytic disorders

Neurodegenerative diseases (ND) are characterized by progressive loss of neuronal function. Mechanisms of ND pathogenesis are incompletely understood, hampering the development of effective therapies. Langerhans cell histiocytosis (LCH) is an inflammatory neoplastic disorder caused by hematopoietic progenitors expressing mitogen-activated protein kinase (MAPK)-activating mutations that differentiate into senescent myeloid cells that drive lesion formation. Some individuals with LCH subsequently develop progressive and incurable neurodegeneration (LCH-ND). Here, we showed that LCH-ND was caused by myeloid cells that were clonal with peripheral LCH cells. Circulating BRAFV600E+ myeloid cells caused the breakdown of the blood-brain barrier (BBB), enhancing migration into the brain parenchyma where they differentiated into senescent, inflammatory CD11a+ macrophages that accumulated in the brainstem and cerebellum. Blocking MAPK activity and senescence programs reduced peripheral inflammation, brain parenchymal infiltration, neuroinflammation, neuronal damage and improved neurological outcome in preclinical LCH-ND. MAPK activation and senescence programs in circulating myeloid cells represent targetable mechanisms of LCH-ND.

Circulating senescent myeloid cells drive blood brain barrier breakdown and neurodegeneration

Neurodegenerative diseases (ND) are characterized by progressive loss of neuronal function. Mechanisms of ND pathogenesis are incompletely understood, hampering the development of effective therapies. Langerhans cell histiocytosis (LCH) is an inflammatory neoplastic disorder caused by hematopoietic progenitors expressing MAPK activating mutations that differentiate into senescent myeloid cells that drive lesion formation. Some patients with LCH subsequently develop progressive and incurable neurodegeneration (LCH-ND). Here, we show that LCH-ND is caused by myeloid cells that are clonal with peripheral LCH cells. We discovered that circulating BRAF V600E + myeloid cells cause the breakdown of the blood-brain barrier (BBB), enhancing migration into the brain parenchyma where they differentiate into senescent, inflammatory CD11a + macrophages that accumulate in the brainstem and cerebellum. Blocking MAPK activity and senescence programs reduced parenchymal infiltration, neuroinflammation, neuronal damage and improved neurological outcome in preclinical LCH-ND. MAPK activation and senescence programs in circulating myeloid cells represent novel and targetable mechanisms of ND.

Intratumoral dendritic cell-CD4+ T helper cell niches enable CD8+ T cell differentiation following PD-1 blockade in hepatocellular carcinoma

Despite no apparent defects in T cell priming and recruitment to tumors, a large subset of T cell rich tumors fail to respond to immune checkpoint blockade (ICB). We leveraged a neoadjuvant anti-PD-1 trial in patients with hepatocellular carcinoma (HCC), as well as additional samples collected from patients treated off-label, to explore correlates of response to ICB within T cell-rich tumors. We show that ICB response correlated with the clonal expansion of intratumoral CXCL13+CH25H+IL-21+PD-1+CD4+ T helper cells ("CXCL13+ TH") and Granzyme K+ PD-1+ effector-like CD8+ T cells, whereas terminally exhausted CD39hiTOXhiPD-1hiCD8+ T cells dominated in nonresponders. CD4+ and CD8+ T cell clones that expanded post-treatment were found in pretreatment biopsies. Notably, PD-1+TCF-1+ (Progenitor-exhausted) CD8+ T cells shared clones mainly with effector-like cells in responders or terminally exhausted cells in nonresponders, suggesting that local CD8+ T cell differentiation occurs upon ICB. We found that these Progenitor CD8+ T cells interact with CXCL13+ TH within cellular triads around dendritic cells enriched in maturation and regulatory molecules, or "mregDC". These results suggest that discrete intratumoral niches that include mregDC and CXCL13+ TH control the differentiation of tumor-specific Progenitor exhasuted CD8+ T cells following ICB.

Molecular states during acute COVID-19 reveal distinct etiologies of long-term sequelae

Post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are debilitating, clinically heterogeneous and of unknown molecular etiology. A transcriptome-wide investigation was performed in 165 acutely infected hospitalized individuals who were followed clinically into the post-acute period. Distinct gene expression signatures of post-acute sequelae were already present in whole blood during acute infection, with innate and adaptive immune cells implicated in different symptoms. Two clusters of sequelae exhibited divergent plasma-cell-associated gene expression patterns. In one cluster, sequelae associated with higher expression of immunoglobulin-related genes in an anti-spike antibody titer-dependent manner. In the other, sequelae associated independently of these titers with lower expression of immunoglobulin-related genes, indicating lower non-specific antibody production in individuals with these sequelae. This relationship between lower total immunoglobulins and sequelae was validated in an external cohort. Altogether, multiple etiologies of post-acute sequelae were already detectable during SARS-CoV-2 infection, directly linking these sequelae with the acute host response to the virus and providing early insights into their development.

Abstract 3611: Determining the regulation and role of mreg DC in tumor immunity

Therapeutic immunity against tumors is mainly determined by effector T cells. However, T cells are educated by antigen presenting cells, particularly dendritic cells (DC), that control their differentiation into effector T cells able to eliminate target tumor cells. While most studies focused on understanding T cell features leading to antitumor immunity, little is known about DC determinants required to promote effective antitumor T cell function. Yet several studies are revealing that high DC content in tumor lesions correlate with improved tumor outcome. To identify molecular program that modulate DC functionality in tissues we have profiled tumor-associated DC in mice and tumor lesions. We identified a DC molecular program that was shared between mice and human tumors that included maturation genes associated with T cell simulation and migration to the LN (Cd40, Cd80, Fscn, Ccr7), but also genes associated with immunoregulation (Socs1, Socs2, Socs3, Cd200, Pd-l1, Pd-l2, Fasl), we named mregDC (Maier et al., Nature 2020). We also found that mregDC molecular state is induced in DC1 and DC2 upon capture of tumor debris leading to hypothesize that the mregDC state identifies antigen-charged DC that migrate to lymphoid structures, including tertiary lymphoid structures and tumor draining lymph nodes to modulate tumor specific T cell immune responses. While the role of DC in the draining LN has been established the exact role played by DC in the tumor microenvironment (TME) remains unclear. To explore this question, we used multiplex imaging and spatial transcriptomics platforms to identify the distribution and spatial interactions of DC1, DC2 and mregDC within the TME of human NSCLC and HCC lesions during treatment with PD-1 blockade. We found that mregDC accumulate within the tertiary lymphoid structures together with activated T cell and B cell programs. We are continuing to generate additional data on mregDC interactions in human tumors tissues that respond or resist PD-1 blockade and should be able to provide additional information at the time of the meeting.

Citation Format: Raphael Mattiuz, Nima Assad, Pauline Hamon, Assaf Magen, Jessica Le Berichel, Miriam Merad. Determining the regulation and role of mreg DC in tumor immunity [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3611.

A phase 1b/2 trial of dupilumab given in conjunction with PD-(L)1 blockade in the treatment of relapsed/refractory metastatic NSCLC

TPS9139

Background: Although tumor microenvironments may contain chronic inflammatory cells, this milieu can lead to immune evasion and may contribute to resistance to immunotherapy. Dendritic cells are one of the most potent cross presenters of tumor antigen, and their function is crucial to tumor-directed adaptive immune responses. In the KP mouse model of lung adenocarcinoma (KrasG12D; Tp53-/-), we recently identified interleukin-4 (IL-4)-driven suppression of tumor-antigen charged dendritic cells, which was similarly seen in human lung cancer lesions (Maier, B., et al., A conserved dendritic-cell regulatory program limits antitumour immunity. Nature, 2020.580[7802]:257-262). In the mouse model IL-4 blockade resulted in an increase in IL-12, IFNg and TNF in CD8+ T cells and decreased tumor burden, and further antitumor activity was seen when combined with PD-L1 blockade. Dupilumab is a fully human monoclonal antibody that blocks the IL-4 receptor alpha subunit which disrupts signaling through receptors for both IL-4 and IL-13; it is FDA approved for patients with multiple atopic conditions, in whom series adverse effects are rare. Based on this pre-clinical data, we hypothesize that the addition of dupilumab to anti PD-(L)1 therapy will be well tolerated, and will rescue the anti-tumor effect of immune checkpoint blockade. Methods: This is a phase Ib/II trial. Patients with relapsed/refractory NSCLC who have received prior anti PD-(L)1 treatment are eligible for enrollment. In a phase Ib safety run-in, six patients will be enrolled in a modified set-dose 3+3 design. If no more than 1 DLT is observed during this phase, the trial will proceed to phase II. Phase II will enroll a further 15 patients in a minimax design with early stopping for futility to a maximum total of 21 patients. Patients will undergo pre-treatment biopsies and peripheral blood sampling prior to receiving 3 doses of dupilumab, administered every three weeks, in conjunction with continuing standard-of-care anti-PD-(L)1 therapy. Patients will undergo repeat biopsies 4 weeks after starting therapy. After completion of dupilumab (at 9 weeks) patients will undergo repeat staging. The primary endpoint of phase 1b is safety as measured by frequency and severity of adverse effects. The primary endpoint of phase II is overall response rate as assessed using RECIST 1.1 criteria at the time of post-dupilumab imaging. Exploratory endpoints include analysis of peripheral blood by CyTOF and O-link, and tissue biopsies will be analyzed by multiplex-IHC and bulk-RNA-sequencing. T cell receptor sequencing will be performed on tumor and matched peripheral blood samples, and circulating tumor DNA will be assessed at multiple time points. Phase 1b is currently enrolling as planned. Clinical trial information: NCT05013450.

Spatial CRISPR genomics identifies regulators of the tumor microenvironment

While CRISPR screens are helping uncover genes regulating many cell-intrinsic processes, existing approaches are suboptimal for identifying extracellular gene functions, particularly in the tissue context. Here, we developed an approach for spatial functional genomics called Perturb-map. We applied Perturb-map to knock out dozens of genes in parallel in a mouse model of lung cancer and simultaneously assessed how each knockout influenced tumor growth, histopathology, and immune composition. Moreover, we paired Perturb-map and spatial transcriptomics for unbiased analysis of CRISPR-edited tumors. We found that in Tgfbr2 knockout tumors, the tumor microenvironment (TME) was converted to a fibro-mucinous state, and T cells excluded, concomitant with upregulated TGFβ and TGFβ-mediated fibroblast activation, indicating that TGFβ-receptor loss on cancer cells increased TGFβ bioavailability and its immunosuppressive effects on the TME. These studies establish Perturb-map for functional genomics within the tissue at single-cell resolution with spatial architecture preserved and provide insight into how TGFβ responsiveness of cancer cells can affect the TME.

Single-cell analysis of human non-small cell lung cancer lesions refines tumor classification and patient stratification

Immunotherapy is a mainstay of non-small cell lung cancer (NSCLC) management. While tumor mutational burden (TMB) correlates with response to immunotherapy, little is known about the relationship between the baseline immune response and tumor genotype. Using single-cell RNA sequencing, we profiled 361,929 cells from 35 early-stage NSCLC lesions. We identified a cellular module consisting of PDCD1⁺CXCL13⁺ activated T cells, IgG⁺ plasma cells, and SPP1⁺ macrophages, referred to as the lung cancer activation module (LCAM^hi). We confirmed LCAM^hi enrichment in multiple NSCLC cohorts, and paired CITE-seq established an antibody panel to identify LCAM^hi lesions. LCAM presence was found to be independent of overall immune cell content and correlated with TMB, cancer testis antigens, and TP53 mutations. High baseline LCAM scores correlated with enhanced NSCLC response to immunotherapy even in patients with above median TMB, suggesting that immune cell composition, while correlated with TMB, may be a nonredundant biomarker of response to immunotherapy.

685 Characterization of molecular and spatial diversity of macrophages in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) has a dismal prognosis, and though checkpoint blocking antibodies have significantly improved patient outcome, many patients remain left out, highlighting the need to identify additional immune target to enhance therapeutic immunity. Macrophages (MF) are an abundant and heterogeneous population in the tumor microenvironment (TME), and are associated with a poor prognosis in multiple tumor types, including HCC, however, their molecular and functional diversity is still poorly understood.

Methods

We analyzed the molecular and spatial organization patterns of immune cells within the TME and adjacent tissue of 26 resected HCC lesions using single-cell RNA sequencing and multiplex immunohistochemistry (IHC).

Results

We found that Kupffer cells, the self-renewing tissue-resident macrophages in liver tissue, are lacking from the TME, which is dominated by monocyte-derived macrophages. ScRNAseq followed by high-resolution clustering identified distinct MF molecular programs within the monocyte-derived macrophage compartment. One MF subset expressed a shared signature with monocytes including FCN1, S100A8 and T cell activation genes like CXCL9 and IL32. Conversely, one subset of MF expressed FOLR2, SEPP1 and genes of the complement (C1Qs), a program shared with KC in the adjacent tissue, and include another intratumoral subset enriched for the expression of TREM2 and GPNMB. Guided by these results, we are developing an IHC antibody panel that allows to visualize distinct MF localization in the TME. Intratumoral MF interface with, and potentially regulate, the T cell compartment within the TME. We are analyzing HCC tumor lesions in our treatment-naïve cohort and in patients treated with neoadjuvant anti-PD-1 therapy (NCT03916627) to study co-localization and direct interaction of MF and T cells using physically-interacting cell sequencing (PICseq). This analysis enables us to identify MF with direct cell-cell contact with T cells, and our preliminary analysis demonstrates an enrichment in MF with an immunosuppressive phenotype. We are also using spatial transcriptomic to map molecular programs of MF in the TME, that we will corroborated with additional patients.

Conclusions

Taken together, our data provide a new understanding of intratumoral MF diversity and highlight the presence of specific immunoregulatory MF programs unique to tumor lesions, with subsets of these MF found to be directly interacting with T cells, potentially modulating anti-tumor responsiveness. Our analysis of resected tumor from anti-PD-1 treated patients, will allow us to correlate MF programs, and direct T cell interaction, with clinical response, and will inform therapeutic trials targeting specific MF populations so as to improve clinical efficacy of cancer immunotherapy.

Trial Registration

NCT03916627

Ethics Approval

Samples of tumor and non-involved liver were obtained from surgical specimens of patients undergoing resection at Mount Sinai Hospital (New York, NY) after obtaining informed consent in accordance with a protocol reviewed and approved by the Institutional Review Board at the Icahn School of Medicine at Mount Sinai (RUTH Human Subjects Electronic Submission System 18–00407 and 20–04150) and in collaboration with the Biorepository and Department of Pathology.

Acute COVID-19 gene-expression profiles show multiple etiologies of long-term sequelae

Two years into the SARS-CoV-2 pandemic, the post-acute sequelae of infection are compounding the global health crisis. Often debilitating, these sequelae are clinically heterogeneous and of unknown molecular etiology. Here, a transcriptome-wide investigation of this new condition was performed in a large cohort of acutely infected patients followed clinically into the post-acute period. Gene expression signatures of post-acute sequelae were already present in whole blood during the acute phase of infection, with both innate and adaptive immune cells involved. Plasma cells stood out as driving at least two distinct clusters of sequelae, one largely dependent on circulating antibodies against the SARS-CoV-2 spike protein and the other antibody-independent. Altogether, multiple etiologies of post-acute sequelae were found concomitant with SARS-CoV-2 infection, directly linking the emergence of these sequelae with the host response to the virus.

Downregulation of exhausted cytotoxic T cells in gene expression networks of multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and pathology of multiple organs in individuals under 21 years of age in the weeks following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although an autoimmune pathogenesis has been proposed, the genes, pathways and cell types causal to this new disease remain unknown. Here we perform RNA sequencing of blood from patients with MIS-C and controls to find disease-associated genes clustered in a co-expression module annotated to CD56dimCD57+ natural killer (NK) cells and exhausted CD8+ T cells. A similar transcriptome signature is replicated in an independent cohort of Kawasaki disease (KD), the related condition after which MIS-C was initially named. Probing a probabilistic causal network previously constructed from over 1,000 blood transcriptomes both validates the structure of this module and reveals nine key regulators, including TBX21, a central coordinator of exhausted CD8+ T cell differentiation. Together, this unbiased, transcriptome-wide survey implicates downregulation of NK cells and cytotoxic T cell exhaustion in the pathogenesis of MIS-C.

Abstract 64: Characterization of molecular and spatial diversity of macrophages in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) has a dismal prognosis, and though checkpoint blocking antibodies have significantly improved patient outcome, many patients remain left out highlighting the need to identify additional immune target to enhance therapeutic immunity. Macrophages (MΦ) are an abundant and heterogeneous population in the tumor microenvironment (TME), and are associated with a poor prognosis in multiple tumor types, including HCC, however, their molecular and functional diversity is still poorly understood. We analyzed the molecular and spatial organization patterns of immune cells within the TME and adjacent tissue of 23 resected HCC lesions using single-cell RNA sequencing and multiplex immunohistochemistry (IHC). Strikingly, we found that Kupffer cells (KC), the self-renewing tissue-resident macrophages in liver tissue, are lacking from the TME, which is dominated by monocyte-derived macrophages. ScRNAseq followed by high-resolution clustering identified distinct MΦ molecular programs within the monocyte-derived macrophage compartment. One MΦ subset expressed a shared signature with monocytes including FCN1 and S100A8 in addition to T cell activation genes like CXCL9 and IL32. Conversely, one subset of MΦ expressed FOLR2, SEPP1 and genes of the complement (C1Qs), a program shared with KC in the adjacent uninvolved tissue, and include another intratumoral subset enriched for the expression of TREM2 and GPNMB. Guided by these results, we are developing an IHC antibody panel that allows to visualize distinct MΦ localization in the TME. Intratumoral MΦ interface with, and potentially regulate, the T cell compartment within the TME. We are analyzing HCC tumor lesions in our treatment-naïve cohort and in patients treated with neoadjuvant anti-PD-1 therapy (NCT03916627) to study direct interaction of MΦ and T cells using an innovative technology of physically-interacting cell sequencing (PICseq). This analysis enables us to identify MΦ with direct cell-cell contact with T cells, and our preliminary analysis demonstrates an enrichment in MΦ with a highly immunosuppressive phenotype. We are also using PICseq to map the molecular program of MΦ that are interacting with T cells and that we will corroborated with additional patients. Taken together, our data provide a new understanding of intratumoral MΦ diversity and highlight the presence of specific immunoregulatory MΦ programs unique to tumor lesions, with subsets of these MΦ found to be directly interacting with T cells, potentially modulating anti-tumor responsiveness. Our analysis of resected tumor from anti-PD-1 treated patients, will allow us to correlate MΦ programs, and direct T cell interaction, with clinical response, and will inform therapeutic trials targeting specific MΦ populations so as to improve clinical efficacy of cancer immunotherapy.

Citation Format: Pauline Hamon, Assaf Magen, Merav Cohen, Alexandra Tabachnikova, Christie Chang, Mark Buckup, Leanna Troncoso, Zhen Zhao, Joel Kim, Amir Giladi, Nausicaa Malissen, Fiona Desland, Jessica Le Berichel, Ido Amit, Ephraim Kenigsberg, Myron Schwartz, Thomas Marron, Miriam Merad. Characterization of molecular and spatial diversity of macrophages in hepatocellular carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 64.

Abstract 515: Heterogeneity of PD-1hi CD8 T cells associates with response to PD-1 blockade in hepatocellular carcinoma

Blockade of the PD-1 pathway is a therapeutic strategy to reinvigorate T cell responses against tumors, and when combined with other biologic therapies in the first line setting this achieves significant clinical response in about 25% of hepatocellular carcinoma (HCC) patients. We hypothesize that phenotypic diversity of tumor infiltrating T cells can explain, at least partially, the disparate clinical responses to immunotherapy. Here, we analyze the molecular diversity of T cells in tumor, adjacent tissue and tumor-draining lymph node (dLN) by single-cell RNA sequencing of tissue from 23 patients with early stage HCC treated by neoadjuvant PD-1 blockade (NCT03916627; final data to be submitted to AACR 2021 as late breaking abstract). We identify distinct subsets of PD-1hi CD8 T cells with varying degrees of exhaustion and effector gene programs. Compared to parallel analysis of untreated HCC tumors, we observed that PD-1 blockade resulted in expansion of PD-1hi CD8 T cells in the tumor, regardless of clinical response. PD-1hi CD8 T cell subsets were highly clonal and enriched in the tumor compared to adjacent tissue, suggesting specificity to tumor antigens. Remarkably, we find an association between a specific transcriptomic phenotype within PD-1hi CD8 T cells and response to PD-1 blockade. Using T cell receptor (TCR) sequencing to study the differentiation patterns between T cell states, we found that expanded clonotypes present among PD-1hi CD8 T cells were also found in CD8 effector cells; these data identify characteristic clonally related T cell populations that are enriched in clinical responders. Furthermore, we find the dLNs harbor clonotypes of PD-1hi CD8 T cells expanded in the tumor. These clonally expanded CD8 T cells in dLN display features of activation and exhaustion, suggesting a continuous role of dLN in anti-tumor responses. These results will be corroborated with 8 additional patient samples in which we will further analyze the molecular characteristics of tumor-specific CD8 T cells in the TME and dLN. Furthermore, our ongoing sequencing of pre- and post-treatment lesions will enable measurement of CD8 T cell clonal expansion and molecular programs during response or resistance to PD-1 blockade. These results should help unravel the contribution of dLN T cells to response to immunotherapy as well as identify T cell molecular programs at baseline that correlate with response.

Citation Format: Assaf Magen, Pauline Hamon, Etienne Humblin, Alessandra Schanoski, Joel Kim, Jessica Le Berichel, Ephraim Kenigsberg, Myron Schwartz, Thomas Marron, Alice Kamphorst, Miriam Merad. Heterogeneity of PD-1hi CD8 T cells associates with response to PD-1 blockade in hepatocellular carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 515.

Cytotoxic lymphocytes are dysregulated in multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and multiple organ involvement in individuals under 21 years following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To identify genes, pathways and cell types driving MIS-C, we sequenced the blood transcriptomes of MIS-C cases, pediatric cases of coronavirus disease 2019, and healthy controls. We define a MIS-C transcriptional signature partially shared with the transcriptional response to SARS-CoV-2 infection and with the signature of Kawasaki disease, a clinically similar condition. By projecting the MIS-C signature onto a co-expression network, we identified disease gene modules and found genes downregulated in MIS-C clustered in a module enriched for the transcriptional signatures of exhausted CD8 + T-cells and CD56 dim CD57 + NK cells. Bayesian network analyses revealed nine key regulators of this module, including TBX21 , a central coordinator of exhausted CD8 + T-cell differentiation. Together, these findings suggest dysregulated cytotoxic lymphocyte response to SARS-Cov-2 infection in MIS-C.