Epigenetic tuning of PD-1 expression improves exhausted T cell function and viral control.

PD-1 is a key negative regulator of CD8+ T cell activation and is highly expressed by exhausted T cells in cancer and chronic viral infection. Although PD-1 blockade can improve viral and tumor control, physiological PD-1 expression prevents immunopathology and improves memory formation. The mechanisms driving high PD-1 expression in exhaustion are not well understood and could be critical to disentangling its beneficial and detrimental effects. Here, we functionally interrogated the epigenetic regulation of PD-1 using a mouse model with deletion of an exhaustion-specific PD-1 enhancer. Enhancer deletion exclusively alters PD-1 expression in CD8+ T cells in chronic infection, creating a 'sweet spot' of intermediate expression where T cell function is optimized compared to wild-type and Pdcd1-knockout cells. This permits improved control of chronic infection without additional immunopathology. Together, these results demonstrate that tuning PD-1 via epigenetic editing can reduce CD8+ T cell dysfunction while avoiding excess immunopathology.

Primary nasal influenza infection rewires tissue-scale memory response dynamics.

The nasal mucosa is often the initial site of respiratory viral infection, replication, and transmission. Understanding how infection shapes tissue-scale primary and memory responses is critical for designing mucosal therapeutics and vaccines. We generated a single-cell RNA-sequencing atlas of the murine nasal mucosa, sampling three regions during primary influenza infection and rechallenge. Compositional analysis revealed restricted infection to the respiratory mucosa with stepwise changes in immune and epithelial cell subsets and states. We identified and characterized a rare subset of Krt13+ nasal immune-interacting floor epithelial (KNIIFE) cells, which concurrently increased with tissue-resident memory T (TRM)-like cells. Proportionality analysis, cell-cell communication inference, and microscopy underscored the CXCL16-CXCR6 axis between KNIIFE and TRM cells. Secondary influenza challenge induced accelerated and coordinated myeloid and lymphoid responses without epithelial proliferation. Together, this atlas serves as a reference for viral infection in the upper respiratory tract and highlights the efficacy of local coordinated memory responses.

Alveolar, Endothelial, and Organ Injury Marker Dynamics in Severe COVID-19.

Rationale: Alveolar and endothelial injury may be differentially associated with coronavirus disease (COVID-19) severity over time. Objectives: To describe alveolar and endothelial injury dynamics and associations with COVID-19 severity, cardiorenovascular injury, and outcomes. Methods: This single-center observational study enrolled patients with COVID-19 requiring respiratory support at emergency department presentation. More than 40 markers of alveolar (including receptor for advanced glycation endproducts [RAGE]), endothelial (including angiopoietin-2), and cardiorenovascular injury (including renin, kidney injury molecule-1, and troponin-I) were serially compared between invasively and spontaneously ventilated patients using mixed-effects repeated-measures models. Ventilatory ratios were calculated for intubated patients. Associations of biomarkers with modified World Health Organization scale at Day 28 were determined with multivariable proportional-odds regression. Measurements and Main Results: Of 225 patients, 74 (33%) received invasive ventilation at Day 0. RAGE was 1.80-fold higher in invasive ventilation patients at Day 0 (95% confidence interval [CI], 1.50-2.17) versus spontaneous ventilation, but decreased over time in all patients. Changes in alveolar markers did not correlate with changes in endothelial, cardiac, or renal injury markers. In contrast, endothelial markers were similar to lower at Day 0 for invasive ventilation versus spontaneous ventilation, but then increased over time only among intubated patients. In intubated patients, angiopoietin-2 was similar (fold difference, 1.02; 95% CI, 0.89-1.17) to nonintubated patients at Day 0 but 1.80-fold higher (95% CI, 1.56-2.06) at Day 3; cardiorenovascular injury markers showed similar patterns. Endothelial markers were not consistently associated with ventilatory ratios. Endothelial markers were more often significantly associated with 28-day outcomes than alveolar markers. Conclusions: Alveolar injury markers increase early. Endothelial injury markers increase later and are associated with cardiorenovascular injury and 28-day outcome. Alveolar and endothelial injury likely contribute at different times to disease progression in severe COVID-19.

Primary nasal viral infection rewires the tissue-scale memory response.

The nasal mucosa is frequently the initial site of respiratory viral infection, replication, and transmission. Recent work has started to clarify the independent responses of epithelial, myeloid, and lymphoid cells to viral infection in the nasal mucosa, but their spatiotemporal coordination and relative contributions remain unclear. Furthermore, understanding whether and how primary infection shapes tissue-scale memory responses to secondary challenge is critical for the rational design of nasal-targeting therapeutics and vaccines. Here, we generated a single-cell RNA-sequencing (scRNA-seq) atlas of the murine nasal mucosa sampling three distinct regions before and during primary and secondary influenza infection. Primary infection was largely restricted to respiratory mucosa and induced stepwise changes in cell type, subset, and state composition over time. Type I Interferon (IFN)-responsive neutrophils appeared 2 days post infection (dpi) and preceded transient IFN-responsive/cycling epithelial cell responses 5 dpi, which coincided with broader antiviral monocyte and NK cell accumulation. By 8 dpi, monocyte-derived macrophages (MDMs) expressing Cxcl9 and Cxcl16 arose alongside effector cytotoxic CD8 and Ifng-expressing CD4 T cells. Following viral clearance (14 dpi), rare, previously undescribed Krt13+ nasal immune-interacting floor epithelial (KNIIFE) cells expressing multiple genes with immune communication potential increased concurrently with tissue-resident memory T (TRM)-like cells and early IgG+/IgA+ plasmablasts. Proportionality analysis coupled with cell-cell communication inference, alongside validation by in situ microscopy, underscored the CXCL16-CXCR6 signaling axis between MDMs and effector CD8 T cells 8dpi and KNIIFE cells and TRM cells 14 dpi. Secondary influenza challenge with a homologous or heterologous strain administered 60 dpi induced an accelerated and coordinated myeloid and lymphoid response without epithelial proliferation, illustrating how tissue-scale memory to natural infection engages both myeloid and lymphoid cells to reduce epithelial regenerative burden. Together, this atlas serves as a reference for viral infection in the upper respiratory tract and highlights the efficacy of local coordinated memory responses upon rechallenge.

Cryoablation and post-progression immune checkpoint inhibition in metastatic melanoma: a phase II trial.

Image-guided percutaneous cryoablation is an established minimally invasive oncologic treatment. We hypothesized that cryoablation may modify the immune microenvironment through direct modulation of the tumor, thereby generating an anti-tumor response in tumors refractory to immune checkpoint inhibition (ICI). In this non-randomized phase II single-center study (NCT03290677), subjects with unresectable melanoma progressing on ICI underwent cryoablation of an enlarging metastasis, and ICI was continued for a minimum of two additional cycles. The primary endpoints were safety, feasibility and tumor response in non-ablated lesions. From May 2018 through July 2020, 17 patients were treated on study. The study met its primary endpoints with the combination strategy found to be safe and feasible with an objective response rate of 23.5% and disease control rate of 41% (4 partial response, 3 stable disease). Our data support further study of this synergistic therapeutic approach.

Characterization of cell states in biliary tract cancers identifies mechanisms of therapeutic resistance in a phase II trial of DKN-01/nivolumab.

Biliary tract cancers demonstrate profound therapeutic resistance, and broadly effective therapies for refractory disease are lacking. We conducted a single-arm, second-line phase II trial combining DKN-01, a humanized monoclonal antibody targeting Dickkopf-1 (DKK-1), and nivolumab to treat patients with advanced biliary tract cancer (NCT04057365). No objective responses were seen. To identify mechanisms of treatment failure, we analyzed paired pre-treatment and on-treatment biopsies using scRNA-seq and constructed a detailed molecular classification of malignant and immune cells. We annotated five biliary tract cancer malignant cell states: classical, basal, mesenchymal, neural-like, and endothelial-like. Neural-like and endothelial-like states, which drive therapeutic resistance in other cancers, have not previously been described in BTC. Malignant cell states co-varied with distinct immune cell states, revealing diverse mechanisms of myeloid and T-cell mediated immune suppression, including M2 myeloid and terminally exhausted T cell programs that were induced by DKN-01/nivolumab. Here, we provide the first systematic classification of functionally annotated cell states in biliary tract cancer and provide new insight into resistance mechanisms to an immunotherapy combination that can inform the next generation of trials.

Protocol for bulk RNA sequencing of enriched human neutrophils from whole blood and estimation of sample purity.

Although neutrophils are the most abundant leukocyte in healthy individuals and impact outcomes of diseases ranging from sepsis to cancer, they remain understudied due to technical constraints of isolation, preservation, and sequencing. We present a modified Smart-Seq2 protocol for bulk RNA sequencing of neutrophils enriched from whole blood. We describe steps for neutrophil isolation, cDNA generation, library preparation, and sample purity estimation via a bioinformatic approach. Our approach permits the collection of large cohorts and enables detection of neutrophil transcriptomic subtypes. For complete details on the use and execution of this protocol, please refer to LaSalle et al. (2022)1 and Boribong et al. (2022).2.

Downregulation of KEAP1 in melanoma promotes resistance to immune checkpoint blockade.

Immune checkpoint blockade (ICB) has demonstrated efficacy in patients with melanoma, but many exhibit poor responses. Using single cell RNA sequencing of melanoma patient-derived circulating tumor cells (CTCs) and functional characterization using mouse melanoma models, we show that the KEAP1/NRF2 pathway modulates sensitivity to ICB, independently of tumorigenesis. The NRF2 negative regulator, KEAP1, shows intrinsic variation in expression, leading to tumor heterogeneity and subclonal resistance.

Neutrophil profiles of pediatric COVID-19 and multisystem inflammatory syndrome in children.

Multisystem inflammatory syndrome in children (MIS-C) is a delayed-onset, COVID-19-related hyperinflammatory illness characterized by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antigenemia, cytokine storm, and immune dysregulation. In severe COVID-19, neutrophil activation is central to hyperinflammatory complications, yet the role of neutrophils in MIS-C is undefined. Here, we collect blood from 152 children: 31 cases of MIS-C, 43 cases of acute pediatric COVID-19, and 78 pediatric controls. We find that MIS-C neutrophils display a granulocytic myeloid-derived suppressor cell (G-MDSC) signature with highly altered metabolism that is distinct from the neutrophil interferon-stimulated gene (ISG) response we observe in pediatric COVID-19. Moreover, we observe extensive spontaneous neutrophil extracellular trap (NET) formation in MIS-C, and we identify neutrophil activation and degranulation signatures. Mechanistically, we determine that SARS-CoV-2 immune complexes are sufficient to trigger NETosis. Our findings suggest that hyperinflammatory presentation during MIS-C could be mechanistically linked to persistent SARS-CoV-2 antigenemia, driven by uncontrolled neutrophil activation and NET release in the vasculature.

Translesion DNA synthesis mediates acquired resistance to olaparib plus temozolomide in small cell lung cancer.

In small cell lung cancer (SCLC), acquired resistance to DNA-damaging therapy is challenging to study because rebiopsy is rarely performed. We used patient-derived xenograft models, established before therapy and after progression, to dissect acquired resistance to olaparib plus temozolomide (OT), a promising experimental therapy for relapsed SCLC. These pairs of serial models reveal alterations in both cell cycle kinetics and DNA replication and demonstrate both inter- and intratumoral heterogeneity in mechanisms of resistance. In one model pair, up-regulation of translesion DNA synthesis (TLS) enabled tolerance of OT-induced damage during DNA replication. TLS inhibitors restored sensitivity to OT both in vitro and in vivo, and similar synergistic effects were seen in additional SCLC cell lines. This represents the first described mechanism of acquired resistance to DNA damage in a patient with SCLC and highlights the potential of the serial model approach to investigate and overcome resistance to therapy in SCLC.

Longitudinal characterization of circulating neutrophils uncovers phenotypes associated with severity in hospitalized COVID-19 patients.

Mechanisms of neutrophil involvement in severe coronavirus disease 2019 (COVID-19) remain incompletely understood. Here, we collect longitudinal blood samples from 306 hospitalized COVID-19+ patients and 86 controls and perform bulk RNA sequencing of enriched neutrophils, plasma proteomics, and high-throughput antibody profiling to investigate relationships between neutrophil states and disease severity. We identify dynamic switches between six distinct neutrophil subtypes. At days 3 and 7 post-hospitalization, patients with severe disease display a granulocytic myeloid-derived suppressor cell-like gene expression signature, while patients with resolving disease show a neutrophil progenitor-like signature. Humoral responses are identified as potential drivers of neutrophil effector functions, with elevated severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific immunoglobulin G1 (IgG1)-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirm that while patient-derived IgG antibodies induce phagocytosis in healthy donor neutrophils, IgA antibodies predominantly induce neutrophil cell death. Overall, our study demonstrates a dysregulated myelopoietic response in severe COVID-19 and a potential role for IgA-dominant responses contributing to mortality.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Neutrophil Profiles of Pediatric COVID-19 and Multisystem Inflammatory Syndrome in Children.

Multisystem Inflammatory Syndrome in Children (MIS-C) is a delayed-onset, COVID-19-related hyperinflammatory systemic illness characterized by SARS-CoV-2 antigenemia, cytokine storm and immune dysregulation; however, the role of the neutrophil has yet to be defined. In adults with severe COVID-19, neutrophil activation has been shown to be central to overactive inflammatory responses and complications. Thus, we sought to define neutrophil activation in children with MIS-C and acute COVID-19. We collected samples from 141 children: 31 cases of MIS-C, 43 cases of acute pediatric COVID-19, and 67 pediatric controls. We found that MIS-C neutrophils display a granulocytic myeloid-derived suppressor cell (G-MDSC) signature with highly altered metabolism, which is markedly different than the neutrophil interferon-stimulated gene (ISG) response observed in pediatric patients during acute SARS-CoV-2 infection. Moreover, we identified signatures of neutrophil activation and degranulation with high levels of spontaneous neutrophil extracellular trap (NET) formation in neutrophils isolated from fresh whole blood of MIS-C patients. Mechanistically, we determined that SARS-CoV-2 immune complexes are sufficient to trigger NETosis. Overall, our findings suggest that the hyperinflammatory presentation of MIS-C could be mechanistically linked to persistent SARS-CoV-2 antigenemia through uncontrolled neutrophil activation and NET release in the vasculature. ONE SENTENCE SUMMARY: Circulating SARS-CoV-2 antigen:antibody immune complexes in Multisystem Inflammatory Syndrome in Children (MIS-C) drive hyperinflammatory and coagulopathic neutrophil extracellular trap (NET) formation and neutrophil activation pathways, providing insight into disease pathology and establishing a divergence from neutrophil signaling seen in acute pediatric COVID-19.

SARS-CoV-2 Viremia is Associated with Distinct Proteomic Pathways and Predicts COVID-19 Outcomes.

BACKGROUND: Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) plasma viremia has been associated with severe disease and death in coronavirus disease 2019 (COVID-19) in small-scale cohort studies. The mechanisms behind this association remain elusive. METHODS: We evaluated the relationship between SARS-CoV-2 viremia, disease outcome, inflammatory and proteomic profiles in a cohort of COVID-19 emergency department participants. SARS-CoV-2 viral load was measured using qRT-PCR based platform. Proteomic data were generated with Proximity Extension Assay (PEA) using the Olink platform. RESULTS: Three hundred participants with nucleic acid test-confirmed COVID-19 were included in this study. Levels of plasma SARS-CoV-2 viremia at the time of presentation predicted adverse disease outcomes, with an adjusted odds ratio (aOR) of 10.6 (95% confidence interval [CI] 4.4, 25.5, P<0.001) for severe disease (mechanical ventilation and/or 28-day mortality) and aOR of 3.9 (95%CI 1.5, 10.1, P=0.006) for 28-day mortality. Proteomic analyses revealed prominent proteomic pathways associated with SARS-CoV-2 viremia, including upregulation of SARS-CoV-2 entry factors (ACE2, CTSL, FURIN), heightened markers of tissue damage to the lungs, gastrointestinal tract, endothelium/vasculature and alterations in coagulation pathways. CONCLUSIONS: These results highlight the cascade of vascular and tissue damage associated with SARS-CoV-2 plasma viremia that underlies its ability to predict COVID-19 disease outcomes.

SARS-CoV-2 viremia is associated with distinct proteomic pathways and predicts COVID-19 outcomes.

BACKGROUNDSARS-CoV-2 plasma viremia has been associated with severe disease and death in COVID-19 in small-scale cohort studies. The mechanisms behind this association remain elusive.METHODSWe evaluated the relationship between SARS-CoV-2 viremia, disease outcome, and inflammatory and proteomic profiles in a cohort of COVID-19 emergency department participants. SARS-CoV-2 viral load was measured using a quantitative reverse transcription PCR-based platform. Proteomic data were generated with Proximity Extension Assay using the Olink platform.RESULTSThis study included 300 participants with nucleic acid test-confirmed COVID-19. Plasma SARS-CoV-2 viremia levels at the time of presentation predicted adverse disease outcomes, with an adjusted OR of 10.6 (95% CI 4.4-25.5, P < 0.001) for severe disease (mechanical ventilation and/or 28-day mortality) and 3.9 (95% CI 1.5-10.1, P = 0.006) for 28-day mortality. Proteomic analyses revealed prominent proteomic pathways associated with SARS-CoV-2 viremia, including upregulation of SARS-CoV-2 entry factors (ACE2, CTSL, FURIN), heightened markers of tissue damage to the lungs, gastrointestinal tract, and endothelium/vasculature, and alterations in coagulation pathways.CONCLUSIONThese results highlight the cascade of vascular and tissue damage associated with SARS-CoV-2 plasma viremia that underlies its ability to predict COVID-19 disease outcomes.FUNDINGMark and Lisa Schwartz; the National Institutes of Health (U19AI082630); the American Lung Association; the Executive Committee on Research at Massachusetts General Hospital; the Chan Zuckerberg Initiative; Arthur, Sandra, and Sarah Irving for the David P. Ryan, MD, Endowed Chair in Cancer Research; an EMBO Long-Term Fellowship (ALTF 486-2018); a Cancer Research Institute/Bristol Myers Squibb Fellowship (CRI2993); the Harvard Catalyst/Harvard Clinical and Translational Science Center (National Center for Advancing Translational Sciences, NIH awards UL1TR001102 and UL1TR002541-01); and by the Harvard University Center for AIDS Research (National Institute of Allergy and Infectious Diseases, 5P30AI060354).

Longitudinal characterization of circulating neutrophils uncovers distinct phenotypes associated with disease severity in hospitalized COVID-19 patients.

Multiple studies have identified an association between neutrophils and COVID-19 disease severity; however, the mechanistic basis of this association remains incompletely understood. Here we collected 781 longitudinal blood samples from 306 hospitalized COVID-19 + patients, 78 COVID-19 âˆ' acute respiratory distress syndrome patients, and 8 healthy controls, and performed bulk RNA-sequencing of enriched neutrophils, plasma proteomics, cfDNA measurements and high throughput antibody profiling assays to investigate the relationship between neutrophil states and disease severity or death. We identified dynamic switches between six distinct neutrophil subtypes using non-negative matrix factorization (NMF) clustering. At days 3 and 7 post-hospitalization, patients with severe disease had an enrichment of a granulocytic myeloid derived suppressor cell-like state gene expression signature, while non-severe patients with resolved disease were enriched for a progenitor-like immature neutrophil state signature. Severe disease was associated with gene sets related to neutrophil degranulation, neutrophil extracellular trap (NET) signatures, distinct metabolic signatures, and enhanced neutrophil activation and generation of reactive oxygen species (ROS). We found that the majority of patients had a transient interferon-stimulated gene signature upon presentation to the emergency department (ED) defined here as Day 0, regardless of disease severity, which persisted only in patients who subsequently died. Humoral responses were identified as potential drivers of neutrophil effector functions, as enhanced antibody-dependent neutrophil phagocytosis and reduced NETosis was associated with elevated SARS-CoV-2-specific IgG1-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirmed that while patient-derived IgG antibodies mostly drove neutrophil phagocytosis and ROS production in healthy donor neutrophils, patient-derived IgA antibodies induced a predominant NETosis response. Overall, our study demonstrates neutrophil dysregulation in severe COVID-19 and a potential role for IgA-dominant responses in driving neutrophil effector functions in severe disease and mortality.

Longitudinal proteomic analysis of severe COVID-19 reveals survival-associated signatures, tissue-specific cell death, and cell-cell interactions.

Mechanisms underlying severe coronavirus disease 2019 (COVID-19) disease remain poorly understood. We analyze several thousand plasma proteins longitudinally in 306 COVID-19 patients and 78 symptomatic controls, uncovering immune and non-immune proteins linked to COVID-19. Deconvolution of our plasma proteome data using published scRNA-seq datasets reveals contributions from circulating immune and tissue cells. Sixteen percent of patients display reduced inflammation yet comparably poor outcomes. Comparison of patients who died to severely ill survivors identifies dynamic immune-cell-derived and tissue-associated proteins associated with survival, including exocrine pancreatic proteases. Using derived tissue-specific and cell-type-specific intracellular death signatures, cellular angiotensin-converting enzyme 2 (ACE2) expression, and our data, we infer whether organ damage resulted from direct or indirect effects of infection. We propose a model in which interactions among myeloid, epithelial, and T cells drive tissue damage. These datasets provide important insights and a rich resource for analysis of mechanisms of severe COVID-19 disease.

Early cross-coronavirus reactive signatures of humoral immunity against COVID-19.

The introduction of vaccines has inspired hope in the battle against SARS-CoV-2. However, the emergence of viral variants, in the absence of potent antivirals, has left the world struggling with the uncertain nature of this disease. Antibodies currently represent the strongest correlate of immunity against SARS-CoV-2, thus we profiled the earliest humoral signatures in a large cohort of acutely ill (survivors and nonsurvivors) and mild or asymptomatic individuals with COVID-19. Although a SARS-CoV-2­specific immune response evolved rapidly in survivors of COVID-19, nonsurvivors exhibited blunted and delayed humoral immune evolution, particularly with respect to S2-specific antibodies. Given the conservation of S2 across ß-coronaviruses, we found that the early development of SARS-CoV-2­specific immunity occurred in tandem with preexisting common ß-coronavirus OC43 humoral immunity in survivors, which was also selectively expanded in individuals that develop a paucisymptomatic infection. These data point to the importance of cross-coronavirus immunity as a correlate of protection against COVID-19.

Plasma proteomics reveals tissue-specific cell death and mediators of cell-cell interactions in severe COVID-19 patients.

COVID-19 has caused over 1 million deaths globally, yet the cellular mechanisms underlying severe disease remain poorly understood. By analyzing several thousand plasma proteins in 306 COVID-19 patients and 78 symptomatic controls over serial timepoints using two complementary approaches, we uncover COVID-19 host immune and non-immune proteins not previously linked to this disease. Integration of plasma proteomics with nine published scRNAseq datasets shows that SARS-CoV-2 infection upregulates monocyte/macrophage, plasmablast, and T cell effector proteins. By comparing patients who died to severely ill patients who survived, we identify dynamic immunomodulatory and tissue-associated proteins associated with survival, providing insights into which host responses are beneficial and which are detrimental to survival. We identify intracellular death signatures from specific tissues and cell types, and by associating these with angiotensin converting enzyme 2 (ACE2) expression, we map tissue damage associated with severe disease and propose which damage results from direct viral infection rather than from indirect effects of illness. We find that disease severity in lung tissue is driven by myeloid cell phenotypes and cell-cell interactions with lung epithelial cells and T cells. Based on these results, we propose a model of immune and epithelial cell interactions that drive cell-type specific and tissue-specific damage in severe COVID-19.