Uncovering the neutrophil inflammatory response in severe COVID-19: implications of type I interferon in inflammasome activation

Patient population

Adult clinical patients with confirmed COVID-19 (RT-PCR positive for SARS-CoV-2) at Helsinki University Hospital (HUH) (hospitalized: n = 34; outpatients: n = 8) were enrolled in the present study. Blood samples were collected during hospitalization for the severe COVID-19 group, and after confirmation of diagnosis for the mild COVID-19 outpatient group. Samples for RNA sequencing were collected in 2020 and representing infections by the original and early SARS-CoV-2 variants, whereas samples for ex vivo culture experiments were collected in 2021-2022 likely representing infections by omicron subvariants of SARS-CoV-2. As controls, healthy blood donors were included for RNA sequencing (n = 7, age 57 ± 7, male/female 3/4) and ex vivo culturing experiments (n = 9, age 38 ± 14, male/female 4/5). The study was approved by the Ethics Committee of the Hospital District of Helsinki and Uusimaa (HUS/853/2020, HUS/1238/2020). All volunteers gave a written informed consent, in accordance with the Declaration of Helsinki. For clinical correlation analysis, severe COVID-19 patients were further categorized by severity based on their need for hospitalization and oxygen supplementation, as described previously ⁷. For each patient, medical history and clinical data were collected through retrospective patient record review and are presented for the severe COVID-19, hospitalized patients in Table 1 and as previously described ⁷. Calprotectin was measured from serum (diluted 1:1000) by ELISA, according to the manufacturer’s protocol (calprotectin/S100A8 DuoSet kit, R&D systems).

Isolation of granulocytes from human blood

Blood samples from COVID-19 patients and healthy controls (HC) were collected in EDTA vacutainer tubes and transported to the laboratory. Peripheral blood mononuclear cells (PBMCs) or polymorphonuclear cells (PMNs) were isolated from whole blood by density gradient centrifugation using either Ficoll-Paque Plus (GE Healthcare) or Polymorphprep (Axis-Shield) respectively, following standard procedures including the use of 2 mM EDTA in PBS and red blood cell lysis with ACK lysis buffer (Lonza by Thermo Fisher). Subsequently, isolation of CD66⁺ granulocytes (low-density granulocytes, LDGs) from the PBMC fraction was performed using the CD66abce MicroBead Kit (Miltenyi Biotec, Germany) with an MS column, according to the manufacturer’s instructions. Both the positively selected CD66⁺ LDGs and the isolated PMNs were then washed and counted, using a TC20^TM Automated Cell Counter (Bio-Rad Laboratories, Inc.) with trypan blue staining for dead cell exclusion. All described procedures in this section were done at room temperature. An aliquot of cells was lysed in Trizol reagent (Thermo Fisher Scientific, USA) and stored at –80°C for later extraction of total RNA and subsequent RNA sequencing (RNA-seq) analysis.

Caspase1 activity

Caspase1 activity was assessed in isolated cells after 2 h of culture (1 million cells/ml) using the caspase-Glo^®1 inflammasome assay (Promega) according to the manufacturer’s protocol, with 2.5 µM nigericin (Invivogen) treatment as the activator. The resulting luminescence was measured by a Hidex Sense microplate reader (Hidex).

Soluble factor stimulation assays

Isolated granulocytes from HC and COVID-19 patients were cultured at 2 million cells/ml in RPMI 1640 supplemented with 10% fetal bovine serum (R10) at 37°C. Cells were primed (1^st signal) with either LPS (20 ng/ml, Sigma Aldrich) or IFN-I (combination of 2.7*10⁴ IU/ml IFN-α and IFN-β, Immunotools) for 4 h, followed by activation (2^nd signal) by 2.5 µM nigericin or monosodium urate crystals (MSU, 100 µg/ml, Invivogen) for an additional 4 h. For the 24 h stimulation experiments, nigericin was added to the cultured cells, in the presence or absence of inflammasome inhibitors MCC995 (2 µg/ml) and Ac-YVAD-FMK (20 µg/ml, both from Invivogen).Cells were pelleted by centrifugation at 400 G for 5 min and stored in Trizol at – 80°C for later RNA extraction whereas supernatants were used to measure IL-1β, IL-18, myeloperoxidase (MPO) and IL-8 by ELISAs according to the manufacturer’s protocols (DuoSet kits from R&D systems). LDH was measured in supernatants using Cyquant LDH cytotoxicity assay (ThermoFisher). HL-60 cells (ATCC #CCL-240) were activated similarly to neutrophils after a 5-day differentiation period induced by 1 % DMSO.

Isolated mouse neutrophils were stimulated in the same fashion as human neutrophils. However, 1 million cells/ml were used, and only mouse IFN-α at 1 µg/ml served as the IFN-I stimulus. IL-1β and IL-18 levels were measured by DuoSet kits (R&D systems).

Virus propagation

The SARS-CoV-2 hCoV-19/Finland/THL-202117309/2021 (delta strain B.1.617.2) and the mouse-adapted strain MaVie ¹⁸ were propagated in VeroE6-TMPRSS2 cells (kidney epithelial cells expressing the transmembrane protease serine 2) ¹⁹ grown in DMEM supplemented with 10% inactivated FCS, 100 IU/mL Penicillin, 100 μg/mL Streptomycin and 2 mM L-glutamine at 37°C. The virus was purified from supernatants by ultracentrifugation (SW28 rotor, 27,000 rpm, 90 min, +4 °C) through a 0.22 μm-filtered 30% ultra-pure sucrose cushion (in PBS), to obtain virus preparations free of cell culture contaminants. Virus titers were calculated by the median tissue culture infectious dose (TCID50) after assessing cytopathic effects by crystal violet staining of cell cultures infected for 5 days with serially diluted virus.

RNA sequencing

Neutrophils isolated from different cohorts comprised three PMN groups (severe COVID-19, mild COVID-19, and healthy controls), and one LDG group (given that these cells were rare in mild COVID-19 patients and HC, only LDGs from patients with severe COVID-19 were included).

cDNA synthesis from total RNA was performed according to Takara SMARTseq v4 Ultra-low input RNA kit for Sequencing user manual (Takara Bio, Mountain View, CA, USA) followed by Illumina Nextera XT Library preparation according to Illumina Nextera XT Reference Guide (Illumina, San Diego, CA, USA). UDI index setup was used for the Nextera XT libraries. Library quality check was performed using LabChip GX Touch HT High Sensitivity assay (PerkinElmer, USA) and libraries were pooled based on the concentrations acquired from the assay. The pooled libraries were quantified for sequencing using KAPA Library Quantification Kit (KAPA Biosystems, Wilmington, MA, USA) and sequenced on the Illumina NovaSeq6000 system for 200 cycles using S1 flow cell (Illumina, San Diego, CA, USA). Read length for the paired-end run was 2×101 bp.

RNA data analysis

Principal Component Analysis (PCA) was performed on the normalized gene expression matrix using the Chipster software ²⁰ to identify patterns in the data and reduce the dimensionality of the dataset. The top principal components were selected based on the percentage of variance explained. For enrichment analyses, Gene Set Enrichment Analysis (GSEA) and Over-Representation Analysis (ORA) were performed on the top 5000 DE genes identified by DESeq2 (adjusted P value < 0.05, log2FC >1) using ExpressAnalyst ²¹. GSEA was used to identify enriched signaling pathways using the Reactome database, while ORA was used to identify enriched pathways using the KEGG database. The resulting p-values were corrected for multiple testing using the Benjamini-Hochberg method, and pathways with a corrected p-value <0.05 were considered significant.

To visualize the expression patterns of the differentially expressed (DE) genes, the data was analyzed using the AltAnalyze software ²², which selected the top 95 genes based on correlation and determined the heatmap clustering, using the Euclidean distance metric and the complete linkage method. Then, the obtained heatmap was re-generated using heatmapper.ca ²³ for better visualization.

CIBERSORTx, a machine learning algorithm that infers cell type proportions using a reference gene expression matrix of known cell types (25) was used to perform RNA-seq deconvolution on the gene expression data to estimate the abundance of immune cell types in the samples ²⁴. The signature matrix used was taken from Lasalle et al. ⁸. This reference matrix made use of a published whole-blood single-cell dataset ⁹, and included the main immune cell types: monocytes, NK cells, T lymphocytes, B lymphocytes, plasmablasts and neutrophils, the latter subclassified into mature and immature. The smaller subsets of granulocytes (eosinophils and basophils) are not considered separately and are most likely categorized as neutrophils in the bulk data deconvolution. Nonetheless, the resulting cell type proportions were used to compare the immune cell composition between groups.

Additionally, the determination of sample purity (>80% of “mature neutrophils” for PMNs or “immature neutrophils” for LDGs) served as a limiting parameter for the visualization of differentially expressed inflammasome related genes from the RNA sequencing results, which were selected and graphed in a heatmap using heatmapper.ca ²³, clustered by complete linkage and ordered by Spearman’s rank.

Volcano Plots

To visualize differentially expressed (DE) genes between groups from human and mice RNA-seq results previously identified by DESeq2, a volcano plot was generated using GraphPad Prism. Genes with a P-adjusted value (padj or FDR) <0.05 were considered significant. Similarly, RNA sequencing data from GSE93996 ²⁵ was reanalyzed, and all DE genes in ATRA-differentiated HL-60 cells were visualized in a volcano plot.

Single cell transcriptomics data analysis

This study made use of the “COVID-19 Immune Atlas: integration of 5 public COVID-19 PBMC single-cell datasets” available online ²⁶. This standardized data collection contains cells from different assays (10x 3’ v2, 10x 3’ v3, 10x technology and Seq-Well) and consists of a total of 239,696 cells from the peripheral blood, 3,693 of which are neutrophils. These neutrophils were further subclassified as mature (59%) and immature (41%), based on the immune atlas predetermined cell classes. This was confirmed by a CD16b expression in mature neutrophils, and a higher CD66b expression in the immature population. This data was obtained from and analyzed in the Chan Zuckerberg CELLxGENE platform ²⁶.

Reverse transcription and quantitative PCR (RT-qPCR) for human selected human genes

Total RNA was extracted from unstimulated or ex vivo stimulated PMNs using the Trizol reagent (Invitrogen, USA) according to the manufacturer’s protocol. Subsequently, cDNA synthesis was performed using the RevertAid RT Reverse Transcription Kit (Thermo Scientific, USA) as per the manufacturer’s instructions. Quantitative PCR (qPCR) was performed using the Stratagene model (Agilent Technologies) and SYBR Green/ROX master mix (Thermo Scientific, USA). The primer sequences for qPCR are presented in Supplementary Table S1.

Primer specificity was confirmed using melting curve analysis and dissociation curves. The relative expression levels of the genes of interest were calculated using the 2-ΔΔCT method and normalized to the expression of the housekeeping gene GAPDH. Baseline gene expressions of unstimulated samples were statistically assessed using the Mann-Whitney test, while the two-way ANOVA Tukey’s multiple comparisons test was performed for the ex vivo stimulated samples.

Mouse infections

Experimental procedures were approved by the Animal Experimental Board of Finland (license number ESAVI/28687/2020). Female BALB/c mice (Envigo, Indianapolis, IN, USA; 7 to 8 weeks, n = 36 in total) were transferred to the University of Helsinki biosafety level-3 (BSL-3) facility and acclimatized to individually ventilated biocontainment cages (ISOcage; Scanbur, Karl Sloanestran, Denmark) for 7 days with ad libitum water and food (rodent pellets). For subsequent experimental infection, the mice were placed under isoflurane anesthesia and inoculated intranasally with 50 µL of SARS-CoV-2 MaVie strain (5*10⁵ TCID50/animal) or PBS (mock-infected control). Daily weighting of all mice was performed, and their well-being was carefully monitored for signs of illness (e.g., changes in posture or behavior, rough coat, apathy, ataxia). Euthanasia was performed by cervical dislocation under terminal isoflurane anesthesia. All animals were dissected immediately after euthanasia, and the lungs were sampled for multiple downstream analyses. The infections were performed as 3 separate experiments (exp) with 12 mice each: 1) Exp 1 included 8 mice infected with MaVie and 4 mock infected mice. At 2 days post infection (dpi), 4 infected and the mock infected mice were euthanized; the remaining infected mice were euthanized at 4 dpi. The right lung was sampled for virus-specific RT-qPCR (1/5) and neutrophil isolation (4/5), the left lung was fixed for histological and immunohistochemical examination. 2) Exp 2 included 8 infected and 4 mock infected mice of which half were euthanized at 2 dpi and 4 dpi, respectively. From these mice, both lung lobes were subjected to neutrophil isolation. 3) Exp 3 included 8 mice that were infected and immediately inoculated intraperitoneally with 250 µg of anti-mouse IFNAR-1 (n = 4) or IgG1 isotype control (n = 4) (Bio-X-Cell, New Hampshire, USA), and 4 mock-infected animals. All mice were euthanized at 2 dpi. Each 1 1/5 of the left lobe was processed for virus-specific RT-qPCR and histology and immunohistochemistry respectively. The remaining approx. 80% of the lungs served for neutrophil isolation.

Neutrophil isolation from mouse lungs

Neutrophil isolation was performed from the lungs of all mice. The dissected lung tissue was chopped into small pieces using scissors and enzymatically digested with a cocktail of Liberase (50 ug/ml; Roche #05401020001 from Merck) and DnaseI (100 ug/ml; Roche #11284932001 from Merck) in RPMI-1640 for 30 min at 37 °C. The resulting homogenate was diluted 10-fold in R10 and passed through a 70 µm Cell strainer (Pluriselect) to obtain a single-cell suspension. Neutrophils were isolated by positive selection using Ly6G-binding magnetic beads and MS columns according to the manufacturer’s recommendations (Miltenyi Biotec). Neutrophils were isolated with a purity exceeding 95% based on flow cytometry analysis of Ly6G expression. Subsequently, the isolated neutrophils were utilized for RNA isolation or inflammasome activation, as described below.

RNA sequencing of mouse neutrophils

Mouse neutrophils were isolated, lysed in Trizol (Thermo Scientific) and the RNA extracted in the liquid phase using chloroform. RNA isolation was carried out using the Rneasy micro kit (Qiagen). Isolated RNA (1 ng) underwent whole transcriptome sequencing with ribodepletion. Briefly, RNA sequencing was performed using the Illumina Stranded with RiboZero library preparation method. Sample quality and integrity were assessed using TapeStation RNA analysis. Sequencing was conducted on the Illumina NextSeq platform, followed by standard bioinformatics analysis for gene expression quantification.

The service was provided by the Biomedicum Functional Genomics Unit at the Helsinki Institute of Life Science and Biocenter Finland at the University of Helsinki.

RT-qPCR of mouse samples

RNA was extracted from dissected lung samples (1/10 of the whole lung) of mice in Exp 1 and Exp 3 as well as isolated neutrophils using Trizol (Thermo Scientific) following the manufacturers’ instructions. The isolated RNA was directly subjected to one-step RT-qPCR analysis based on a previously described protocol using primer-probe sets detecting the viral genome encoding for the RNA-dependent RNA polymerase (RdRp) ²⁷, subgenomic E ²⁸ as well as mouse caspase1, IL1b and GAPDH (Applied biosystems #Mm00438023_m1, #Mm00434228_m1 and #Mm99999915_g1 respectively, Thermo scientific). The PCRs were performed with TaqPath 1-step master mix (ThermoFisher Scientific) using AriaMx instrumentation (Agilent, Santa Clara, CA, USA).

Histology and immunohistochemistry

From animals in Exp 1 and Exp 3 the whole left lung (Exp 1) or 1/5 of the left lung (Exp 3) were trimmed for histological examination and routinely paraffin wax embedded. Consecutive sections (3 µm) were prepared and routinely stained with hematoxylin-eosin (HE) or subjected to immunohistochemistry (IHC) for the detection of SARS-CoV-2 nucleoprotein (NP) ²⁹ and Ly6G (neutrophil marker); for Exp 3, a further section of the infected lungs was stained for histone H3 (NET marker) ³⁰. All stains followed previously published protocols ³¹.

Morphometric analyses

For quantification of SARS-CoV-2 antigen expression and the extent of neutrophil influx into the lungs, a morphometric analysis was undertaken on the slides stained for SARS-CoV-2 NP and Ly6G, respectively. The stained slides were scanned using NanoZoomer 2.0-HT (Hamamatsu, Hamamatsu City, Japan), and several sections of the lung of each animal were quantitatively analysed using the Visiopharm 2022.01.3.12053 software (Visiopharm, Hoersholm, Denmark). The morphometric analysis served to quantify the area, in all lung sections of an animal, that showed immunostaining for viral NP and Ly6G, respectively. In Visiopharm, for each section, the lung was manually outlined and annotated as a Region Of Interest (ROI), manually excluding artifactually altered areas. The manual tissue selection was further refined with an Analysis Protocol Package (APP) based on a Decision Forest classifier, with the pixels from the ROI being ultimately classified as either “Tissue” or “Background”. This new “Tissue” ROI, regrouping the different lung samples analysed for each animal, was further quantified by executing two APPs successively. The first APP was based on a Threshold classifier and served to detect and outline areas with immunostaining. The second APP then measured both the surface of the immunostained area (µm²) and the surface of the “Tissue” ROI (µm²). The percentage of immunostained area (%), expressed as the ratio between the immunostained area and the total area, was obtained for each animal in Excel (Microsoft Office 2019; Microsoft, Redmond, Washington, United States), according to the following formula: ([positive area (µm2)]/ [total area (µm2)]) x 100.

Statistical analyses

Statistical analysis was performed using GraphPad Prism 8.3 software (GraphPad Software, San Diego, CA, USA) and R software v3.6.3 (R core team). Statistically significant correlations between parameters were assessed by calculating Spearman’s correlation coefficients, and differences between groups were assessed with Mann-Whitney, Kruskall-Wallis or ordinary one-way or 2-way ANOVA tests, depending on sample distribution and the number of groups analyzed.

Unsupervised RNA-seq analysis reveals a gene expression signature of circulating neutrophils in COVID-19 that is strongly influenced by maturity

With our recent findings on increased frequencies of low-density granulocytes (LDGs, isolated from the PBMC fraction) during COVID-19 and their likely relevant role in disease progression ⁷, we sought to understand in more detail how the transcriptomic profile of LDGs differs from their higher “normal” density counterpart, the circulating polymorphonuclear cells (PMNs) ³², typically consisting mainly of mature neutrophils. Neutrophils isolated from different cohorts comprised three PMN groups (severe COVID-19, mild COVID-19, and healthy controls), and one LDG group. We explored gene expression through RNA sequencing (RNA-seq) and confirmed the distinctive gene signature of LDGs by principal component analysis (PCA) (Figure 1A), which revealed that the gene expression patterns of COVID-19 LDGs differed from those of all PMNs regardless of the patients’ disease state.

• Download figure
• Open in new tab

Figure 1. Comparison of gene expression in granulocyte populations of COVID-19 patients using RNA-seq analysis.

(A) Principal component analysis (PCA) of the RNA-seq samples (n = 7 for HC PMNs, n = 19 for COVID-19 PMNs (n = 11 and 8 for severe and mild, respectively), n = 6 for severe COVID-19 LDGs). (B) Ridgeline diagrams depicting the top 10 enriched signal pathways from the genes differentially expressed by LDGs versus PMNs during severe COVID-19, including an overrepresentation analysis (ORA) using KEGG database and gene-set enrichment analysis (GSEA) according to Reactome database. Both enrichment analyses were made using ExpressAnalyst and are sorted by P-value, obtained from Welch’s t-test. (C) Heatmap of the top 95 differentially expressed genes between PMNs from healthy controls, mild and severe COVID-19, as well as LDGs from severe disease, identified by unsupervised ICGS analysis based on correlation, using AltAnalyze software.

Functional enrichment analyses through gene overrepresentation and gene-set enrichment analyses (ORA and GSEA, respectively) compared LDGs with PMNs from severe COVID-19 patients (Figure 1B). When considering the top 5000 differentially expressed genes (DEG), DNA replication and cell cycle were the most overrepresented pathways identified by ORA (Figure 1B), which supported our previous findings suggesting LDGs to be predominantly immature cells ⁷. However, the next statistically significant result was an underrepresentation of the NLR signaling pathway genes in LDGs in contrast with PMNs, highlighting that LDGs have an inflammatory profile that differs from that of PMNs, also depicted in the GSEA (Figure 1B). Briefly, type I and II interferon signaling were among the most significantly enriched gene sets in PMNs compared to LDGs. Further analysis showed that PMNs from severe COVID-19 patients had significantly enriched DEGs in inflammatory pathways such as NLR, tumor necrosis factor (TNF), retinoic acid-inducible gene-I (RIG-I) like, autophagy, MAPK and JAK/STAT signaling pathways compared to either HC PMNs or severe COVID-19 LDGs (Supplementary Fig. 1A-D). Interestingly, many of the observed DEGs were linked to inflammasome formation (Supplementary Fig. 1E-F). However, mild COVID-19 PMNs did not significantly differ from HC PMNs in their inflammatory profile (Supplementary Fig. 1G-H).

Additionally, we performed unsupervised Iterative Clustering and Guide Gene Selection (ICGS), which identified the top DEGs and classified the samples into two major clusters: a first one containing all isolated LDG samples, and a second one comprising all isolated PMN samples (Figure 1C). The former cluster consisted of many neutrophil marker genes (e.g. CEACAM8, MMP8, DEFA), while in the latter cluster the most significantly upregulated genes in the severe COVID-19 subgroup were mainly interferon inducible (e.g. OAS1, OAS2, IFI44L, IFITM3, IFIT2) as well as some inflammasome related (CASP5, CARD16) genes.

Mature neutrophils from severe COVID-19 display a strong IFN-related gene expression signature

To gain a comprehensive understanding of the cellular composition within PMN and LDG fractions, we performed a deconvolution of the RNA-seq data, which verified that most cells present in the samples were neutrophils. This analysis also demonstrated that cells in the LDG fraction were predominantly immature neutrophils, meanwhile PMNs were composed of mainly mature neutrophils (Figure 2A). Subsequently, gene expression analysis of PMNs and LDGs containing high purity of mature neutrophils and immature neutrophils (over 80%), respectively, confirmed that mature neutrophils from severe COVID-19 displayed an increased expression of type I IFN related genes, while immature neutrophils from severe COVID-19 distinctively lacked this signature (Figure 2B).

• Download figure
• Open in new tab

Figure 2. Increased IFN-I related gene expression in mature COVID-19 neutrophils.

(A) Deconvoluted RNA-seq data. The cellular composition in isolated PMN and LDG fractions was estimated using CIBERSORTx through the identification of cell populations based on RNA-seq. The bar plots in the figure represent the cell composition of each RNA-seq sample, offering insights on sample purity.

(B) Heatmap of differentially expressed IFN-related genes in COVID-19 PMNs and LDGs as compared to HC PMNs. RNA sequencing was performed on purified PMNs from healthy controls, mild COVID-19 and severe COVID-19, as well as severe COVID-19 LDGs. The heatmap was reduced to include only the samples with the highest purity for their cell type, determined by a cell fraction over 0.8 of mature or immature neutrophils for PMNs and LDGs, respectively, as identified by CIBERSORTx. The heatmap was clustered by complete linkage and ordered by Spearman’s rank. n = 13 COVID-19 PMNs (7 severe, 6 mild), n = 4 severe COVID-19 LDGs and n = 5 HC PMNs.

Inflammasomes are activated in severe COVID-19 PMNs, but not directly by SARS-CoV-2

Looking more closely into PMN fractions, ORA identified the inflammasome related NOD-like and RIG-like receptor signaling pathways among the most significantly overrepresented pathways, differentially expressed in severe COVID-19 PMNs versus HC PMNs (Figure 3A) or mild COVID-19 PMNs (Figure 3B). The increased expression of selected IFN-I (OAS1, OAS2, and IFIT1) and inflammasome related genes (CASP1, CASP5, NLRC5 and NAIP) was confirmed by RT-qPCR. However, some inflammasome related genes (IL-1β, NLRP3 and NLRC4) were seemingly downregulated, although not statistically significant (Supplementary Fig. 2).

• Download figure
• Open in new tab

Figure 3. Inflammasome related gene expression and cytokine secretion in PMNs during severe COVID-19.

(A-B) Ridgeline diagrams of overrepresentation analyses (ORA) according to Reactome and KEGG databases, depicting the top 10 enriched signaling pathways in PMNs during severe COVID-19 compared to (A) healthy controls and (B) mild COVID-19. (C) IL-1β and (D) IL-18 levels in 24-h cell culture supernatants from COVID-19 (n = 11 for IL-1β and 9 for IL-18) and HC PMNs (n =6 for both). (E) IL-1β and (F) IL-18 levels in 24-h cell culture supernatant from PMNs exposed or non-exposed to purified SARS-CoV-2 viral particles (10 virus particles / PMN) (n = 3). (G) Caspase1 activity in PMNs following a 2-h stimulation with nigericin or purified SARS-CoV-2 viral particles (10 virus particles / PMN). For HC PMNs, n = 9 for mock and nigericin and n = 6 for SARS-CoV-2 exposure. For COVID-19 PMNs, n = 12 for mock and nigericin and n = 9 for SARS-CoV-2 exposure. *p < 0.05 and **p < 0.01. Data presented as mean ± SD. Tukey’s multiple comparisons test for mixed-effect analysis was applied for (G), meanwhile P values for (C-F) were calculated with the Mann-Whitney U-test.

Given the strong upregulation of many inflammasome related genes during severe COVID-19, we assessed whether PMNs exhibit active inflammasome formation in vivo. To evaluate spontaneous inflammasome mediated cytokine secretion, fresh PMNs isolated from severe COVID-19 patients and HC were cultured ex vivo overnight. We measured the levels of IL-1β and IL-18 in the supernatant and found that IL-1β secretion was significantly increased in the supernatant of severe COVID-19 PMNs compared to HC PMNs (Figure 3C), whereas the IL-18 levels did not differ significantly (Figure 3D). Additionally, since SARS-CoV-2 viral particles were previously implicated to induce inflammasome formation in macrophages ¹⁷, the IL-1β and IL-18 levels after HC PMNs exposure to SARS-CoV-2 were also assessed but no significant effects in the secretion of these cytokines were observed (Figure 3E, F).

The spontaneous secretion of IL-1β by COVID-19 PMNs suggests that these cells are actively producing and releasing IL-1β through inflammasome formation which is dependent on caspase1 activity ³³. We assessed caspase1 activity in response to the second signal required for inflammasome activation, induced by nigericin, and observed increased caspase1 activity in severe COVID-19 PMNs compared to HC PMNs (Figure 3G). These findings suggest that severe COVID-19 PMNs have an increased capacity for inflammasome activation, potentially due to an existing priming signal during acute disease in vivo. However, no significant difference in caspase1 activity between non-exposed and virus-exposed PMNs were observed (Figure 3G), indicating that caspase1 activation in COVID-19 PMNs is not directly triggered by the virus.

Activation of neutrophil inflammasome related pathways during respiratory distress is not specific to COVID-19

We also reanalyzed the RNA-seq data previously generated by LaSalle et al. ⁸, focusing on neutrophil transcriptomics in patients with COVID-19 as compared to non-COVID-19 patients exhibiting similar respiratory symptoms, as well as to healthy controls. Our analysis included IFN-α response, IL-1β production, TLR signaling, NLRP3 inflammasome, and pyroptosis pathways, using the Gene Ontology (GO) database; the NLR signaling pathway using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; and inflammasome pathway using the REACTOME database (Figure 4). These pathways were significantly enriched in COVID-19 patients, supporting our findings. Importantly, the genes from the above-mentioned pathways were also induced in non-COVID-19 patients, suggesting that these pathways represent a general neutrophil response to inflammatory stimuli rather than a COVID-19 specific response.

• Download figure
• Open in new tab

Figure 4. Comparative neutrophil transcriptomics of COVID-19 and non-COVID-19 patients.

Bar graphs represent the activation levels of selected pathways and processes as identified by neutrophil transcriptomics. The analysis includes IFN-α responses, IL-1β production, TLR signaling, NLRP3 inflammasomes, and pyroptosis as determined through the Gene Ontology (GO) database. The NOD-like receptor signaling pathway was investigated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and the inflammasome pathway was explored via the REACTOME database. The graphs compare the activation levels of these pathways in healthy controls (HC), non-COVID patients with similar symptoms (COVID-19 negative), and COVID-19 positive individuals. Statistical significance is denoted as follows: *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. P values were calculated with Kruskall-Wallis test.

Type I IFNs prime PMNs for inflammasome activation

Since PMNs from COVID-19 patients concomitantly display a strong IFN-I signature (Figure 2B) and an increased propensity for inflammasome activation, we hypothesized that IFN-I could act as the priming signal for PMN inflammasomes during COVID-19. Isolated HC PMNs were stimulated ex vivo with exogenous IFN-I and the well-described inflammasome priming (1^st signal) and activator (2^nd signal) agents LPS and nigericin, respectively ^34,35. After stimulation, both priming signals induced pro-IL-1β (31 kDa) in the cell lysates, followed by the release of active IL-1β (17 kDa) into the supernatant in response to nigericin (Figure 5A), confirming the ability of IFN-I to prime PMNs for inflammasome activation.

• Download figure
• Open in new tab

Figure 5. IFN-I primes inflammasome activation while COVID-19 PMNs show defective inflammasome responses ex vivo.

Isolated HC or COVID-19 PMNs were non-stimulated or stimulated 4h with IFN-I (combination of 2.7*10⁴ IU/ml IFN-α and IFN-β) or 20 ng/ml LPS (1^st signal), followed by 4h with 2.5 µM nigericin or purified SARS-CoV-2 (10:1 virus/PMNs) (2^nd signal). Then, (A) western blot of pro-IL-1β (31 kD) and active IL-1β (17 kD) was performed from HC PMNs supernatant and cell lysates, (B) IL-1β (n = 5 HC PMN and 9 COVID-19 PMN) and (C) MPO (n= 5 HC PMN and 9 COVID-19 PMN) were measured from supernatants by ELISA. (D-E) Effect of inflammasome inhibitor MCC950 (2 µg/ml, added simultaneously with nigericin) on IL-1β secretion in HC (D) and severe COVID-19 PMN (E) supernatant (n = 3). (F) LDH and (G) IL-8 in HC and severe COVID-19 PMN supernatants (n = 3). (H-K) RT-qPCR of selected mRNAs in IFN-I or LPS-primed HC and COVID-19 PMNs (n = 6-8 HC PMN and 7-10 COVID-19 PMN). *p < 0.05, **p < 0.01, ***p < 0.001, **** p < 0.0001; n = 2-4. Data were presented as mean ± SD. P values calculated with Kruskall-Wallis test for the comparison between treatments by group (HC or COVID-19 PMNs), and Mann-Whitney test for the comparison between HC and COVID-19 PMNs by individual treatment for (B-G), and Two-way ANOVA Tukey’s multiple comparisons test for (B, H-K). Data presented as mean ± SD.

To assess inflammasome formation in circulating neutrophils during COVID-19, PMNs from HC and COVID-19 patients underwent similar stimulation assays as above, followed by IL-1β measurement from supernatants by ELISA. In addition, to further assess the role of SARS-CoV-2 virus particles in neutrophil inflammasome activation, HC PMNs were cultured in the presence of purified viruses (10 infectious units/PMN). HC PMNs responded to both LPS and IFN-I by increasing their IL-1β secretion, which was exponentiated after exposure to nigericin (Figure 5A), confirming the ability of IFN-Is to prime for inflammasome assembly in PMNs, albeit less efficiently than LPS. Interestingly, COVID-19 PMNs produced less IL-1β than HC PMNs upon exogenous inflammasome activation primed by either LPS or IFN-I, while SARS-CoV-2 particles did not have any effect on PMN inflammasome activation (Figure 5B). As with 24-h cultures (Figure 3D), we did not detect any significant changes in IL-18 expression in either HC or COVID-19 PMNs (data not shown). However, myeloperoxidase (MPO), a marker of neutrophil degranulation and NETosis ³⁶, responded to nigericin indistinctively between COVID-19 PMNs and HC PMNs, and therefore the observed diminished IL-1β release by COVID-19 PMNs is not due to general cellular inertia but may be specific to the ex vivo induced inflammasome pathway. Furthermore, additional stimulation assays in the presence of the NLRP3 inhibitor MCC950 (Figures 5D-E) and caspase1 inhibitor YVAD (Supplementary Fig. 3A) confirmed that induced IL-1β secretion is dependent on canonical NLRP3 inflammasome activation. However, unlike IL-1β (Supplementary Fig. 3B), increased IL-18 secretion was not detectable even after 24-h stimulation (Supplementary Fig. 3C). Furthermore, the observed residual IL-18 was not affected by inflammasome inhibitors, suggesting its secretion to be unrelated to inflammasome activity in PMNs.

We further assessed the specificity of inflammasome activation by measuring LDH and IL-8 levels in the supernatants from the same cells and under the same experimental conditions as shown in Figures 5D-E. The measurements of the former were done to assess inflammasome mediated cell death by pyroptosis in response to nigericin, while the latter was assessed to demonstrate the responsiveness of PMNs to an inflammasome unrelated inflammatory cascade. As with IL-1β secretion, COVID-19 PMNs were less responsive than HC PMNs to nigericin- and LPS-mediated LDH (Figure 5F) and IL-8 (Figure 5G) release, respectively. This suggests that COVID-19 PMNs are generally poorly responsive to inflammatory stimuli while retaining the ability for degranulation and NETosis (see above).

To examine this reduced responsiveness to external inflammatory priming, we evaluated the inflammasome related gene expression following ex vivo stimulation with IFN-I or LPS (Figure 5H-I and Supplementary Fig. 3D-G). OAS1 gene, an interferon stimulated gene (ISG), showed significant upregulation by IFN-I in COVID-19 PMNs as compared to HC PMNs (Figure 5H), while the inflammasome related genes IL-1β (Figure 5I), CASP1 (Figure 5J) and NLRC5 (Figure 5K) were more efficiently induced in HC PMNs than COVID-19 PMNs. This suggests that the inflammasome defect in COVID-19 PMNs is at the transcriptional level when using IFN-I as the priming factor, while high OAS1 gene expression indicates transcriptional defect is restricted to individual genes.

Association between ex vivo inflammasome activation and disease severity

Our analysis of the association between ex vivo inflammasome activation (caspase1 activity and IL-1β release) and clinical markers of disease severity, including neutrophil responses, revealed intriguing links. Calprotectin is marker of neutrophil activation or death ³⁷ but also potentially activates the inflammasome (³⁸. A significant positive correlation between calprotectin plasma levels and PMN caspase1 activity (Figure 6A-B) underscores this latter possibility and highlights the interplay between inflammation and inflammasome activation in PMNs of COVID-19 patients. Furthermore, the negative association of PMN IL-1β levels (after ex vivo stimulation) with disease severity (WHO ordinal scale, Figure 6A) and patient neutrophil counts (Figure 6A, C) supports the exhaustion hypothesis, wherein PMNs from severe COVID-19 patients may be less responsive to stimuli due to prior in vivo activation.. While these findings provide intriguing insights into the complex interplay between calprotectin release, caspase1 activity, and inflammasome activation in COVID-19, additional research is required to further elucidate these connections.

• Download figure
• Open in new tab

Figure 6. Correlation analysis between clinical parameters and ex vivo PMN inflammasome activation.

(A) Spearman’s correlation matrix depicting the relationships among clinical parameters and results of ex vivo experimentation. For the WHO ordinal scale, the baseline parameters were used. (B-C) Linear regression analysis demonstrating the associations between: (B) Positive association between PMN Caspase1 activity, measured after ex vivo nigericin stimulation, and the levels of Calprotectin in the matched patient’s peripheral blood; (C) Negative association between ex vivo stimulated PMN IL-1β levels (LPS+Nig) and the blood neutrophil count in matched patients at the time of sampling (n=12).

LOS = length of stay. WHO = World Health Organization. Min = minimum. Casp1 = caspase1. LPS or IFN + nig = lipopolysaccharide or type I interferon + nigericin ex vivo stimulation.

LDGs differ from PMNs in gene expression and release of inflammasome related interleukins

To assess the inflammasome related inflammatory profile of COVID-19 LDGs in comparison to PMNs, we analyzed the differential expression of inflammasome related genes using RNA-seq (Figure 7A). LDGs differed significantly from PMNs, with the most striking difference being their increased expression of IL-18 and NLRC4, whereas PMNs displayed higher levels of IL-1β, NLRP3 and caspases 1, 4 and 5. Single cell sequencing data from the COVID-19 immune atlas confirmed our transcriptomic results (Figure 7B), from which a detailed gene by gene analysis of the most relevant inflammasome related genes is shown (Figure 7C and Supplementary Fig. 4). Briefly, PYCARD gene coding for the ASC protein was expressed similarly in mature and immature neutrophils (Supplementary Fig. 4), confirming that both cell types have inflammasome forming capacity. However, most of the inflammasome gene expressions differed significantly and in the same manner as in our transcriptomic analysis.

• Download figure
• Open in new tab

Figure 7. Immature neutrophils express IL-18 in response to inflammasome activation.

(A) Heatmap depicting selected differentially expressed inflammasome related genes from RNA sequencing performed in PMNs from HC, mild and severe COVID-19, as well as severe COVID-19 LDGs. Only the samples with the highest purity, determined by a cell fraction over 0.8 of mature or immature neutrophils for PMN and LDG, respectively, (identified by CIBERSORTx) are included. The heatmap was clustered by complete linkage and ordered by Spearman’s rank.

(B) UMAP analysis of the COVID-19 Immune Atlas, which integrates 5 public COVID-19 PBMC single-cell transcriptomics datasets, created using CELLxGENE. (Top) UMAP showing the clustering of CD16⁺ cells (mature, FCGR3B expressing cells) and CD66b⁺ cells (immature, CEACAM8 expressing cells). Each dot represents a single cell colored according to the expression level of a selected gene. The color scale ranges from green (low expression) to purple (high expression). (Bottom) Pie chart summarizing the percentage of mature (black) and immature (blue) cells in the dataset.

(C) The fraction of mature and immature neutrophils cells expressing inflammasome related genes identified in Figure 7B are shown in a bar graph. For each gene, the proportion of expressing cells is shown in light blue, while the proportion of negative or not-expressing cells is shown in gray. Zoomed-in bar graph depicts the proportion of mature and immature cells expressing each gene.

(D-G) Isolated COVID-19 LDGs or HL-60 cells (differentiated for 5 days with 1% DMSO) were non-stimulated or stimulated 4h with IFN-I or LPS (1^st signal), followed by 4h with nigericin (2^nd signal) in the presence or absence of inflammasome inhibitors MCC950 or YVAD as previously. Secretion of (D, F) IL-1β and (E, G) IL-18 were measured from the supernatants by ELISA (n = 2 for LDGs and 3-5 for HL-60).

*p < 0.05 and **p < 0.01. P values calculated with Kruskal-Wallis test. Data presented as mean ± SD.

(H) Volcano plot of differentiated vs undifferentiated HL-60 cells gene expression from GSE93996, with inflammasome related genes marked in blue. Only significant DE genes are shown (adjusted p value < 0.05).

We conducted ex vivo stimulation assays using LDGs isolated from COVID-19 patients, similar to the approach used for PMNs described earlier. Like PMNs, IL-1β secretion by LDGs was elevated in the presence of a priming signal (IFN-I or LPS), which exponentially increased when the inflammasome activation signaling molecule nigericin was added (Figure 7D). Contrary to PMNs and in line with the transcriptomics data, an increased IL-18 secretion was detected (Figure 7E). Additionally, the secretion of both ILs by LDGs was inhibited in the presence of inflammasome specific inhibitors MCC950 and YVAD (Figure 7F-G).

These findings suggested that the outcome of neutrophil inflammasome activation varies based on cellular maturation state. To explore this further, we conducted in vitro stimulation studies using differentiated HL-60 cells, an immature neutrophil-like model ³⁹. Similar to LDGs from COVID-19 patients, HL-60 displayed comparable IL-18 secretion pattern upon LPS or IFN-I stimulation and nigericin-induced activation. Notably, their IL-1β release was only detected with LPS priming (Figures 7F-G). Furthermore, consistent with transcriptomic analysis revealing an upregulation of inflammasome related genes upon differentiation (Figure 7H), the capacity of HL-60 cells to secrete inflammasome related cytokines was differentiation-dependent (data not shown). Overall, these findings suggest that neutrophils may lose the ability to secrete IL-18 during maturation, and release of neutrophil-derived IL-18 occurs primarily in disease states associated with extensive granulopoiesis and increased immature granulocyte counts in the blood, like COVID-19 ⁴⁰.

Neutrophils are recruited to the lungs in SARS-CoV-2 infected mice

Hamsters and human ACE2 expressing mice infected with SARS-CoV-2 develop pulmonary inflammation including neutrophil recruitment ^41–43. To further assess the role of neutrophils in COVID-19, we utilized a recently developed SARS-CoV-2 mouse model ¹⁸. This model employs the MaVie strain, serially passaged in mouse lungs and causing pneumonia similar to human COVID-19 in wild-type BALB-C mice ¹⁸. Infected mice started losing weight by day 2 post-infection, with some mice reaching the clinical endpoint of 20% weight loss by day 4 (Supplementary Fig. 5A). Viral loads were significantly higher at 2 dpi than 4 dpi (Figure 8A), and viral antigen expression, widespread at 2dpi in bronchioles and alveoli, matched this pattern (Figure 8B and D). The extensive viral replication at 2 dpi was associated with degeneration of infected epithelial cells, most prominent in the respiratory epithelium, accompanied by neutrophil (Ly6G+) infiltration (Figure 8D) and a significant increase in the number of neutrophils in the lungs of infected mice compared to PBS-inoculated mice (Figure 8C). Neutrophil numbers significantly decreased by 4dpi but remained higher than controls (Figure 8C). Detailed information on the histological and immunohistochemical features of these mice is provided in Supplementary Table S2. Together, these findings suggest a pivotal role of neutrophils in clearing the virus in SARS-CoV-2 infected mice.

• Download figure
• Open in new tab

Figure 8. Neutrophil accumulation in the lungs correlates with viral loads in SARS-CoV-2 infected mice.

Female BALB/c mice were intranasally inoculated with 5*10⁵ TCID50 SARS-CoV-2 MaVie strain or PBS as control and euthanized at 2 dpi or 4 dpi.

(A) RNA was isolated from lungs and subjected to RT-qPCR targeting viral subE and GAPDH as housekeeping gene. The relative expression of subE was measured using the comparative Ct method as compared to mock-infected control (in which subE was undetectable but set to 40 Ct) **** p < 0.0001. P values calculated with Welch’s t-test.

(B) Quantification based on morphometric analysis that determines the area of immunolabelling for SARS-CoV-2 nucleoprotein in relation to total tissue area. ***p < 0.001. P values calculated with Welch’s t-test.

(C) Quantification of Ly-6G based on morphometric analysis that determines the area of immunolabelling for Ly6G in relation to total tissue area in mock-infected controls. **** p < 0.0001. P values calculated with ordinary one-way ANOVA using Tukey’s multiple comparisons test. Black line represents the mean.

(D) Left column: immunohistochemistry for SARS-CoV-2 nucleoprotein; right column: immunohistochemistry for Ly6G (neutrophil marker), hematoxylin counterstain. Bars = 500 µm (large images) and 50 µm (insets).

At 2 dpi (top), the arrow points at a bronchus with viral antigen expression in epithelial cells. A close–up of the bronchus (bottom; B: bronchial lumen) shows degenerated and slough off antigen positive epithelial cells. Adjacent alveoli exhibit viral antigen expression in typeI (arrowhead) and typeII (arrow) pneumocytes. The overview (top) shows neutrophils between the infected bronchial (arrow) epithelial cells, in parenchymal areas (arrowhead; right inset) and in capillaries (arrowheads). A close-up of the bronchus (bottom; B: bronchial lumen) highlights numerous neutrophils between degenerate (arrowheads) epithelial cells.

At 4 dpi (middle), there are focal areas with antigen expression in alveolar epithelial cells and infiltrating macrophages. Neutrophils are present among the infiltrating cells (arrow) as individual cells (inset: arrows) or in aggregates (inset: arrowhead).

The bottom shows the lung of a mock-infected control animal. There is no viral antigen expression. Staining for Ly6G depicts individual neutrophils in larger vessels (inset: arrow) or in capillaries (inset: arrowheads).

Neutrophils from SARS-CoV-2 infected mice display IFN-I dependent caspase1 activation

We isolated neutrophils from the lungs of infected mice at 2 and 4 dpi, as well as from non-infected mice (inoculated with PBS) for transcriptomic analysis by RNA-seq (Exp 1). PCA showed distinct differences between neutrophils from infected and non-infected mice, with slight variation between the 2 and 4 dpi time points (Figure 9A). These differences were reflected in many DEGs, including several IFN-I responsive and inflammasome related genes, which showed strong upregulation at 2 dpi with slightly lower but still significantly elevated levels at 4 dpi, compared to non-infected mice (highlighted in the DEG heatmap; Figure 9B). The volcano plot (Figure 9C) provided a comprehensive view of the DEG pattern between neutrophils from SARS-CoV-2 infected and mock-infected mice. In addition to confirming the upregulation of IFN-I responsive and inflammasome related genes observed in the heatmap, the plot revealed a broader transcriptional response to viral infection with several additional DEG.

• Download figure
• Open in new tab

Figure 9. Neutrophils of SARS-CoV-2 infected mice display increased caspase1 activation ability that is dependent on IFN-I.

Female BALB/c mice were intranasally inoculated with 5* 10⁵ TCID50 SARS-CoV-2 MaVie strain or PBS as control. Lungs were harvested at 2 and 4 dpi and Ly-6G+ neutrophils isolated based on positive selection with magnetic beads. RNA was isolated and subjected to transcriptomic analysis by RNA-seq.

(A) Principal component analysis (PCA) of the PBS-inoculated control and SARS-CoV-2 infected mice lung neutrophil RNA-seq samples. (B) Heatmap of the top DEGs. (C) Volcano plots of DEGs between neutrophils isolated from SARS-CoV-2 infected mice versus uninfected PBS-inoculated mice. Blue points represent significant terms (adjusted p-value < 0.05), while smaller gray points represent non-significant terms. Relevant inflammasome and interferon related genes are shown with larger and darker blue points.

(D) Caspase1 activity in isolated mice neutrophils following a 2-h stimulation with nigericin was assessed by a bioluminescence method (Caspase-Glo® 1 Inflammasome Assay).

(E-I) Mice were intraperitoneally inoculated with 250 µg anti-IFNAR or IgG1 isotype control directly after infection with SARS-CoV-2 and lung neutrophils isolated at 2 dpi (including also intranasally PBS-inoculated control mice without intraperitoneal injection).

(E) Caspase1 activity was assessed following a 2-h stimulation with nigericin by bioluminescence method.

(F-G) RNA was isolated from isolated neutrophils and fold change mRNA expressions of (C) Caspase1 (Casp1) and (D) IL-1β (IL1b) in isotype control and anti-IFNAR treated infected mice as compared to mock-infected control mice assessed by RT-qPCR.

(H-I) The levels of IL-1β and were measured by ELISA in the supernatants of (H) LPS primed or (I) IFN-α primed (4-h) and nigericin-activated (4-h) neutrophils.

DEG = differentially expressed genes.

*p < 0.05, and **p < 0.01. P values for D, E and H panels were calculated with ordinary one-way ANOVA using Tukey’s multiple comparisons. Welch’s t-test was used for panels F, G and I. Data presented as mean ± SD.

Next, we investigated if neutrophils from SARS-CoV-2 infected mice showed increased inflammasome formation, similar to PMNs from COVID-19 patients. Neutrophils from infected mice, harvested at 2 and 4 dpi, displayed increased caspase1 activity upon nigericin stimulation, compared to neutrophils from non-infected mice (Exp 2; Figure 9D). To assess the role of IFN-I, we inoculated mice with an IFN-I blocking anti-IFNAR monoclonal antibody or an isotype control antibody post-infection (Exp 3). Remarkably, neutrophils from anti-IFNAR treated mice showed diminished nigericin-induced caspase1 activity (Figure 9E). Furthermore, caspase1 and IL-1β gene expressions were lower in anti-IFNAR treated than isotype treated mice (Figure 9F-G). Taken together, IFN-I appears to be responsible for the increased caspase1 activity in neutrophils of infected mice.

Similar to human PMNs and LDGs, we assessed the release of IL-1β and IL-18 from murine neutrophils following LPS priming and subsequent nigericin activation. Despite the increased caspase1 activity and elevated IL-1β gene expression (Fig. 9E-G), neutrophils from isotype treated infected mice released IL-1β levels upon activation comparable to those of non-infected control mice, which is in line with previous experiments using human PMNs. However, surprisingly, neutrophils from anti-IFNAR treated infected mice displayed significantly higher IL-1β levels (Figure 9H). As with human PMNs, IL-18 was not secreted in a fashion similar to that of IL-1β from murine neutrophils. Neutrophils from isotype-treated infected mice displayed a slight increase in IL-18 levels compared to uninfected or anti-IFNAR treated infected mice (Supplementary Fig. 5B), in line with their increased caspase1 activity (Figure 9E).

Neutrophils isolated from isotype and anti-IFNAR treated SARS-CoV-2 infected mice, released IL-1β when primed with IFN-α (Figure 9I). However, as with human neutrophils, IL-1β levels after IFN-α priming were significantly lower than after LPS priming. IL-18 secretion, on the other hand, displayed analogous low levels after both IFN-α and LPS priming, which were slightly higher in the isotype-treated than anti-IFNAR treated mice (Supplementary Fig. 5B). Due to a lower neutrophil yield, these additional IFN-α priming experiments were not conducted with mock-infected mice.

Histology of isotype-treated mice and anti-IFNAR treated mice lungs showed comparable features (Supplementary Fig. 5D) as infections without antibody (Exp 1; Figure 8B-D). Regardless of treatment, some neutrophils in infected mice displayed degeneration and NETosis evidenced by histone H3 staining (Supplementary Fig. 5D; Supplementary Table S2), and viral loads remained consistent between anti-IFNAR or IgG1 isotype control treated mice (Supplementary Fig. 5C). Taken together, blocking IFN-I signaling did not alter virus replication, virus-induced pathological changes, or early neutrophil recruitment following infection.