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.

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.

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.

Targeting TBK1 to overcome resistance to cancer immunotherapy.

Despite the success of PD-1 blockade in melanoma and other cancers, effective treatment strategies to overcome resistance to cancer immunotherapy are lacking1,2. Here we identify the innate immune kinase TANK-binding kinase 1 (TBK1)3 as a candidate immune-evasion gene in a pooled genetic screen4. Using a suite of genetic and pharmacological tools across multiple experimental model systems, we confirm a role for TBK1 as an immune-evasion gene. Targeting TBK1 enhances responses to PD-1 blockade by decreasing the cytotoxicity threshold to effector cytokines (TNF and IFNγ). TBK1 inhibition in combination with PD-1 blockade also demonstrated efficacy using patient-derived tumour models, with concordant findings in matched patient-derived organotypic tumour spheroids and matched patient-derived organoids. Tumour cells lacking TBK1 are primed to undergo RIPK- and caspase-dependent cell death in response to TNF and IFNγ in a JAK-STAT-dependent manner. Taken together, our results demonstrate that targeting TBK1 is an effective strategy to overcome resistance to cancer immunotherapy.

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.

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.

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.

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.

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.

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.

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.

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 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.

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.

A unique subset of glycolytic tumour-propagating cells drives squamous cell carcinoma.

Head and neck squamous cell carcinoma (SCC) remains among the most aggressive human cancers. Tumour progression and aggressiveness in SCC are largely driven by tumour-propagating cells (TPCs). Aerobic glycolysis, also known as the Warburg effect, is a characteristic of many cancers; however, whether this adaptation is functionally important in SCC, and at which stage, remains poorly understood. Here, we show that the NAD+-dependent histone deacetylase sirtuin 6 is a robust tumour suppressor in SCC, acting as a modulator of glycolysis in these tumours. Remarkably, rather than a late adaptation, we find enhanced glycolysis specifically in TPCs. More importantly, using single-cell RNA sequencing of TPCs, we identify a subset of TPCs with higher glycolysis and enhanced pentose phosphate pathway and glutathione metabolism, characteristics that are strongly associated with a better antioxidant response. Together, our studies uncover enhanced glycolysis as a main driver in SCC, and, more importantly, identify a subset of TPCs as the cell of origin for the Warburg effect, defining metabolism as a key feature of intra-tumour heterogeneity.

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.

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.

Granzyme B PET imaging of immune-mediated tumor killing as a tool for understanding immunotherapy response.

BACKGROUND: Cancer immunotherapy research is expanding to include a more robust understanding of the mechanisms of treatment response and resistance. Identification of drivers of pro-tumor and anti-tumor immunity during treatment offers new strategies for effective alternative or combination immunotherapies. Currently, tissue or blood samples are collected and analyzed, then dichotomized based on clinical end points that may occur months or years after tissue is collected. While overall survival is ultimately the desired clinical outcome, this dichotomization fails to incorporate the nuances that may occur during an anti-tumor response. By failing to directly measure immune activation at the time of sampling, tumors may be misclassified and potentially obscure important biological information. Non-invasive techniques, such as positron emission tomography (PET), allow for global and quantitative measurements of cancer specific processes and are widely used clinically to help manage disease. METHODS: We have previously developed a novel PET agent that can non-invasively quantify granzyme B release in tumors and have demonstrated its ability to predict response to checkpoint inhibitor therapy in multiple murine models of cancer. Here, we used the quantitative measurement of granzyme B release as a direct and time-matched marker of immune cell activation in order to determine immune cell types and cytokines that correlate with effective checkpoint inhibitor therapy in both tumors and tumor-draining lymph nodes. RESULTS: Through PET imaging, we were able to successfully distinguish distinct microenvironments, based on tumor type, which influenced immune cell subpopulations and cytokine release. Although each tumor was marked by functionally distinct pathways of immune cell activation and inflammation, they also shared commonalities that ultimately resulted in granzyme B release and tumor killing. CONCLUSIONS: These results suggest that discrete tumor immune microenvironments can be identified in both responsive and non-responsive tumors and offers strategic targets for intervention to overcome checkpoint inhibitor resistance.