key: cord-351431-xcqz9d3t
authors: Wang, Jun; Li, Qian; Yin, Yongmei; Zhang, Yingying; Cao, Yingying; Lin, Xiaoming; Huang, Lihua; Hoffmann, Daniel; Lu, Mengji; Qiu, Yuanwang
title: Excessive Neutrophils and Neutrophil Extracellular Traps in COVID-19
date: 2020-08-18
journal: Front Immunol
DOI: 10.3389/fimmu.2020.02063
sha: 
doc_id: 351431
cord_uid: xcqz9d3t

Background: Cases of excessive neutrophil counts in the blood in severe coronavirus disease (COVID-19) patients have drawn significant attention. Neutrophil infiltration was also noted on the pathological findings from autopsies. It is urgent to clarify the pathogenesis of neutrophils leading to severe pneumonia in COVID-19. Methods: A retrospective analysis was performed on 55 COVID-19 patients classified as mild (n = 22), moderate (n = 25), and severe (n = 8) according to the Guidelines released by the National Health Commission of China. Trends relating leukocyte counts and lungs examined by chest CT scan were quantified by Bayesian inference. Transcriptional signatures of host immune cells of four COVID19 patients were analyzed by RNA sequencing of lung specimens and BALF. Results: Neutrophilia occurred in 6 of 8 severe patients at 7–19 days after symptom onset, coinciding with lesion progression. Increasing neutrophil counts paralleled lesion CT values (slope: 0.8 and 0.3–1.2), reflecting neutrophilia-induced lung injury in severe patients. Transcriptome analysis revealed that neutrophil activation was correlated with 17 neutrophil extracellular trap (NET)-associated genes in COVID-19 patients, which was related to innate immunity and interacted with T/NK/B cells, as supported by a protein–protein interaction network analysis. Conclusion: Excessive neutrophils and associated NETs could explain the pathogenesis of lung injury in COVID-19 pneumonia.

As of early May 2020, more than 3 million cases of coronavirus disease 2019 (COVID-19) have been confirmed worldwide, resulting in hundreds of thousands of deaths (1) . According to the Guidelines of the Diagnosis and Treatment of New Coronavirus Pneumonia (version 7) published by the National Health Commission of China, COVID-19 patients can be classified as mild, moderate, and severe cases. Severe patients easily develop acute respiratory distress syndrome (ARDS) or multiple organ failure, with a 4-15% death rate (2, 3) It is not well-understood what drives the exacerbated host response involving a cytokine storm in severe COVID-19 (4) . Specifically, it is unclear what initiates and propagates the cytokine storm. Neutrophil infiltration was noted in three recent reports on the pathological findings from autopsied COVID-19 patients (5) (6) (7) . Neutrophil infiltration in pulmonary capillaries, acute capillaritis with fibrin deposition, extravasation of neutrophils into the alveolar space, and neutrophilic mucositis were observed. Similarly, increased neutrophil counts were reported to occur simultaneously in the peripheral blood of severe and non-surviving COVID-19 patients (3, 8) . Neutrophilia predicts poor outcomes in patients with COVID-19, and our previous research also indicated the neutrophil-to-lymphocyte ratio (NLR) is an independent risk factor for severe disease (8, 9) .

Recently, two serum markers of neutrophil extracellular traps (NETs), myeloperoxidase (MPO)-DNA, and citrullinated histone H3 (Cit-H3) levels were found to be elevated in the serum of COVID-19 patients (10) . This suggested that neutrophilia and excessive NETs may contribute to cytokine release and respiratory failure. As a contributor to pathological inflammation of pneumonia, excessive neutrophils lead to tissue injury by oxidative burst, phagocytosis, and the formation of neutrophil NETs, known as NETosis. NETs are composed of extracellular webs of DNA, histones, microbicidal proteins, and oxidative enzymes that are released by neutrophils to corral infections (11) (12) (13) (14) (15) . The ability of NETs to damage tissues is well-documented in infection and sterile disease. NETs directly kill epithelial and endothelial cells (16, 17) , and excessive NETosis damages the epithelium in pulmonary fungal infection (18) and the endothelium in transfusion-related acute lung injury (19) .

In the present study, first, the dynamics of neutrophil counts in COVID-19 patients (n = 23) during hospitalization were examined, together with the corresponding lung injury, to clinically define the relationship between lung injury and leukocyte counts. Second, transcriptional signatures of host immune cells from COVID-19 patients (n = 4) were analyzed by RNA sequencing of lung specimens or bronchoalveolar lavage fluids (BALF). Immune cell frequency was analyzed by MCPcouter. We used average expression of genes enriched in neutrophil degranulation and activation to screen highly correlated genes and further identified NET associated genes in the correlated gene list to construct an interactive network from the STRING database.

The study was approved by the Ethics Committee of the Fifth People's Hospital, Wuxi (No. 2020-006-1). The 55 confirmed COVID-19 patients were enrolled in this retrospective study from January 23 to March 15, 2020. Written informed consent was obtained from all patients from the Fifth People's Hospital, Wuxi, China.

The clinical handling of COVID-19 patients was performed according to the Guidelines of the Diagnosis and Treatment of New Coronavirus Pneumonia (version 7) published by the National Health Commission of China. Mild, moderate, and severe cases were defined by the following conditions: (1) epidemiological history, (2) fever or other respiratory symptoms, (3) frequency of typical CT image abnormalities of viral pneumonia, and (4) positive RT-PCR result for SARS-CoV-2 RNA. In addition, mild cases were diagnosed if no typical CT image abnormality of viral pneumonia (#3 above) was seen and severe patients also met at least one of the following conditions: (1) shortness of breath, respiratory rate (RR) ≥30 times/min, (2) oxygen saturation (resting state) ≤93%, or (3) PaO 2 /FiO 2 ≤300 mm Hg.

All medical records including epidemiological, demographic, clinical manifestation, laboratory data, radiological characteristics, treatment, and outcome data were reviewed and collected. Laboratory confirmation of SARS-CoV-2 infection was performed by real-time RT-PCR (Bojie Ltd, 119 Shanghai, China) according to Chinese CDC approval. Five sets of RNAseq data from BALF of two COVID-19 patients were acquired from BIG Data Center (accession number CRA002390), and corresponding data of three healthy controls were from the NCBI SRA database (accession numbers SRR10571724, SRR10571730, and SRR10571732). Four RNA-seq data from lung specimens of two COVID-19 patients and two healthy controls were acquired from the GEO database (accession numbers GSM4462416, GSM4462415, GSM4462414, and GSM4462413).

All images were obtained on the CT system (Somatom Definition AS+, Siemens Healthineers, Germany) with patients in supine position. The main scanning parameters were as follows: tube voltage = 120 kV, automatic tube current modulation (about 95 mAs), pitch = 1.2 mm, slice thickness = 7 mm, field of view = 350 mm × 350 mm. All images were then reconstructed with a slice thickness of 0.6 mm with the same increment.

Two professional radiologists (Y.M.Y. and X.M.L.), who were blinded to the laboratory test data, reported chest CT features and assessed the CT features by consensus. The lesion CT values were assessed using the Skyview pacs system. The region-of-interest was selected manually marking the area of highest intensity (most restricted area) of the lesion in CT images.

Kallisto was used to pseudoalign the RNA-seq reads and perform bootstrap analysis using an index based on the ENSEMBL GRCh38 Homo sapiens release 99 transcriptomes (20) . Gene expression levels were then calculated as transcripts per million (TPM). Sleuth (version 0.30.0) (21) was used to perform differential gene expression (DEGs) analysis with the Wald test. Benjamini-Hochberg-adjusted false discovery rate (q < 0.1) was used to correct for multiple comparisons.

To compare lung and BALF samples of COVID-19 patients with healthy controls, differentially expressed genes were exhibited in a scaled heatmap using pheatmap (22) . MCP-counter was used to characterize immune cell subpopulations (23) . The MCP-counter scores obtained from the three underlying transcriptome platforms (Affymetrix Human Genome U133 Plus 2.0, Affymetrix 133A, and Illumina HiSeq) were used to estimate the expression of each cell population. Functional enrichment analysis of the 29 upregulated marker genes of neutrophils was conducted with Metascape (http://metascape.org/) (24) . Gene set enrichment analysis (GSEA) was performed in pre-ranked list mode with 1,000 permutations and weighted enrichment statistic (25) . The gene interaction was analyzed by STRING (26). Gene interaction networks were visualized with eXamine (27) .

Quantitative parameters are described as the median value followed by the inter-quartile range (IQR) in parentheses. Principal component analysis was performed with R package "FactoMineR" to identify those clinical parameters that contribute most to distinguishing severe, moderate, and mild cases of COVID-19 (28) . Figures were produced with R package "ggplot2" (29) . Logistic regression was conducted with R package "rstanarm" (30) to identify associations of laboratory parameters with severity of cases.

Severe cases were typed as severe and others (moderate and mild cases) as non-severe. The generalized linear model was then used to calculate coefficients (mean value with 5%, 95% confidence interval) of all parameters for severe. Finally, we used the function of exp [exp(x) = ex] for coefficients. The results were an odd's ratio (mean, 5-95% credible interval). Receiver operating characteristic curves (ROC) were calculated by R package "pROC." The area under the ROC curve (AUC) and cut-off values of selected parameters were used to distinguish mild and severe cases (31) . Numerical Bayesian linear regression was carried out with Stan using Hamiltonian Monte Carlo (Supplemental Materials; Supplementary Figure 1 ) (32) .

Fifty-five confirmed COVID-19 patients were hospitalized in The Fifth People's Hospital of Wuxi from Jan 23 to Mar 15, 2020. The median age of patients was 45 years (IQR 25-61), and 27 (49%) were male. Based on the previously described guidelines, 22 (40%), 25 (45%), and 8 (15%) of the 55 COVID-19 patients were classified as mild, moderate, and severe cases, respectively. There were five patients with diabetes (9%), 13 with hypertension (24%), eight with surgical history (15%), and two with co-infections (4%). The most common symptoms at onset were fever in 28 cases (51%), sputum production in 13 cases (24%), cough in 22 cases (40%), and fatigue in 17 cases (31%) ( Table 1) .

The clinical handling and relevant time-points of 33 patients including eight severe and 25 moderate cases are shown in Figure 1 . The median time from the date of onset of symptoms to hospital admission, lymphopenia, ARDS, and neutrophilia was 3, 7, 8, and 9 d, respectively. Lymphopenia occurred in seven of eight severe patients and 11 of 25 moderate cases within 7 d, ARDS occurred in all eight severe patients within 8 d, and neutrophilia occurred in six of eight severe patients and one of 25 moderate cases within 9 d (Figure 1) .

The laboratory test of each patient on the day of hospital admission showed that the median neutrophil count in severe COVID-19 patients (3.4, IQR: 1.8-6.7) was higher than in the moderate (3.0, 2.4-3.6) and mild (2.9, 2.3-3.5) groups. In contrast, lymphocyte and monocyte counts in severe COVID-19 patients were lower than in the other two groups ( 

Principal component analysis was performed to visualize the contribution of all mentioned clinical parameters on disease severity (Figure 2A) . Nine variables contributed most strongly. Among them, higher CRP, FIB, neutrophil count, and NLR, and lower lymphocyte count were associated with increased disease severity. These parameters may therefore be used for prognosis.

To assess the diagnostic value of the top two contributors, CRP and lymphocytes, the AUC and cut-off values from the ROC curves were calculated for the severe and mild cases, respectively (Supplementary Figure 2B) . The cut-off values for severe patients were CRP (26.1) and lymphocytes (1.0), and for mild patients the values were CRP (2.2) and lymphocytes (1.4) (see dashed lines in Figure 2B ). Next, dynamic changes of neutrophil, lymphocyte, and monocyte counts in the peripheral blood of COVID-19 patients were monitored (Figure 2C ). Dramatically increased neutrophil counts were found in severe COVID-19 patients in comparison to the other two groups. In contrast, lymphocyte counts persisted at lower values in severe COVID-19 patients. Monocyte counts were lower in severe cases, although the monocyte count fluctuated over a wide range. Timing of the occurrence of maximum neutrophil, minimum lymphocyte, and minimum monocyte counts, and the corresponding counts in COVID-19 patients, during hospitalization are shown in Figure 2D . From day 7 to day 9 after symptom onset, neutrophil counts erupted (>7.7 × 10 9 /L) and peaked in six of eight severe COVID-19 patients. In contrast, only one moderate (1/26) COVID-19 patient was found with neutrophilia. Lymphopenia occurred in seven of eight severe patients but only in four mild (4/22) COVID-19 patients. Monopenia (<1 × 10 8 /L) was found in three moderate (3/25) and four severe (4/8) COVID-19 patients. Overall, monitoring blood cell parameters revealed neutrophilia as a characteristic of severe COVID-19 patients.

Neutrophilia and lymphopenia obviously occurred in severe COVID-19 patients during hospitalization. Here was a case of severe patient. The CRP level remained low when neutrophilia occurred, and the D-dimer levels increased after neutrophilia. Series of chest CT images exhibited enlarged patches and groundglass nodules in the sub-pleura area of both lungs during neutrophilia. Interestingly, all observed lesions were reduced or gradually absorbed along with the return of neutrophils to normal levels after neutrophilia (Figures 3A,B) . The CT value of lesions, reflecting lung lesions, was further demonstrated to have the same trend with neutrophils but the opposite trend with lymphocytes ( Figure 3C) .

To estimate the overall correlation of CT value with neutrophil and lymphocyte counts across patients with a visual inspection of possible trends, linear models were fitted to summarize the dependency of z-values of CT value (CTz, see Figure 3D . Overall, the results showed that the CTz value has no average trend with changing neutrophil and lymphocyte counts for moderate cases (green). However, for the severe cases (red), there are clear trends for CTz value with changing cell counts; specifically, CTz value increased for increasing neutrophil counts, whereas CTz value decreased for increasing lymphocyte counts (Figure 3D ).

Immune cell transcriptional signatures were established from RNA-seq data of BALF and lung specimens of COVID-19 patients and healthy controls. Marker genes of neutrophils, T cells, monocytes, and B cells were identified from Microenvironment Cell Populations-counter (MCP-counter). Their representation in the RNA-seq data were exhibited using a scaled heatmap by comparing both lung and BALF samples of COVID-19 patients to healthy controls ( Figure 4A) .

The results revealed that 112 marker genes represented four immune populations: neutrophils (46 genes), T cells (13 genes), monocytes (10 genes), and B cells (43 genes). For lung 

White blood cell (×10 9 /L) 3.5-9. tissue, the most up-regulated marker genes were enriched in neutrophils, second in monocytes, and only a small proportion were enriched in B cells. Marker genes of T cells were almost all lowly expressed. For BALF, the most upregulated marker genes were similarly enriched in neutrophils, but more up-regulated genes in monocytes and B cells were observed in COVID-19 patients compared to healthy controls, which is different from the lung samples.

Functional enrichment analysis of the 27 upregulated marker genes of neutrophils were further conducted with Metascape. The enrichment analysis revealed that five gene sets with lowest qvalue were related to neutrophil degranulation and activation ( Figure 4B ) and there were 15 marker genes involved. Then, we calculated the average expression of these genes as an evaluating score for neutrophil activation (NAS).

To further assess the abundance of infiltrating immune cells of the lung and BALF in COVID-19 patients, the MCP-counter score was used to quantify the absolute abundance of immune cell subpopulations. Notably, the neutrophil scores were higher and T cell scores were lower in lung samples of COVID-19 patients. The higher abundance of cytotoxic T lymphocytes contributed for cell injury, not for anti-virus. Due to the marker genes for cytotoxic T lymphocytes was KLRC1 (Killer Cell Lectin Like Receptor C1). For the BALF samples, the score of neutrophils, cytotoxic lymphocytes, B cells, monocytes, and dendritic cells were found to be higher in one of the COVID-19 patients compared to the three healthy controls (Figure 4C ).

To explore the outcome of neutrophil activation in COVID-19, we further analyzed the correlation of NAS with 1,363 DEGs that overlapped in both the lung and BALF samples. The spearman correlation was used separately for COVID-19 patients and healthy controls. Then, the R value for every single gene was acquired for COVID-19 patients (R1) and healthy cases (R2). All DEGs were ranked based on R (R1-R2). The "R value" of the top 84 genes (R1 > 0) in the two groups are displayed in Figure 5A . Of these 84 genes, 16 genes were NETs associated genes ( Figure 5B ; Table 3 ) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) Of the 16 genes, LGALS9, HCK, LCP1, CEACAM1 were involved in the cytokine-mediated signaling pathway. S100A8, LGALS9, and CTSC were involved in regulation of apoptotic signal by enrichment annotation from the Metascape tool ( Figure 5B ; Table 3 ) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) .

To further investigate the role of NETs in COVID-19, we generated a gene set termed "NET-associated genes" based on genes coding for proteins enriched in NETs released from human neutrophils with mass spectrometry (Supplementary Table 1) . Pre-ranked GSEA by R resulted in significant enriched gene sets of "NET-associated genes" (Enrichment Score = 0.80) and "Regulation of inflammatory response" (Enrichment Score = 0.72) (Figure 5C ).

As known, the formation of NETs could induce direct lung injury (17). There were 16 NETs associated genes related with neutrophils activation in COVID-19 patients. To further illustrate the interaction between these NETs associated genes with other neutrophils activation related genes, we constructed a protein-to-protein interaction network from the STRING database (Figure 6) . We found that the NETs interacted with STAT1 induced Interferon stimulated genes by IL2RG, implying that NETs associated genes may be triggered by IFN signaling. Besides, NETs in turn may activate B cells via TNFSF13B and inhibit the function of T and NK cells via LGAS9 and CEACAM1, which are negative regulators for T and NK cells.

LGAS9 is a possible promoter of protein-arginine deiminase type 4 (PAD4). PAD4, a key NETs associated gene, lies downstream of ROS and promotes chromatin decondensation (47, 48) . Of note, we also observed ROS related genes including HCK, RAC2, and NCF2 among NETs associated genes (Figure 6) .

To annotate the function of NETs associated genes, they were categorized as metabolic enzymes (RAC2, NCF2), structural proteins (LCP1), anti-microbial related (TREM1), peroxisomal (SH3BGRL3), and others (C1QC, LGALS9, SERPINA1, C1QB, CCL7, CCL8, CEACAM1, HCK, and CXCL16) ( Table 3) . Thus, we speculate that NETs may be activated by innate immunity such as IFN signaling, in COVID-19 patients. NETs may negatively regulate the immune function of T cells and NK cells via LGAS9 and CEACAM1, respectively, leading to insufficient anti-viral immunity and injuring the lung tissue directly.

In this study, a set of laboratory test parameters and the corresponding chest CT images of 55 COVID-19 patients were collected during hospitalization. Among these variables, excessive neutrophils were associated with disease severity, as shown by principal component analysis. Bayesian inference across patients quantified that the increased trend of pneumonia lung injury, as represented by CT values, was in accord with the increased trend in neutrophil counts. Transcriptome analysis of lung specimens and BALF from COVID-19 patients also indicated the most up-regulated marker genes were neutrophil related. Importantly, many neutrophil activation genes were categorized as NET-associated genes. These genes were further assessed to interact with T and NK cells via negative regulatory Functional enrichment analysis of these 84 genes, of which 16 genes were NETs associated genes. (C) Nets associated genes set (Enrichment Score, 0.80) and the GO term of regulation of inflammatory (Enrichment Score, 0.72) by GSEA with DEGs from pre-ranked by R. molecules in COVID-19 patients leading to insufficient anti-viral response and lung injury (Figure 6) .

Our previous study also found an increased neutrophil-tolymphocyte ratio in the most severe disease cases (9) . Recently, neutrophil infiltration was also noted in the lung tissue of autopsied COVID-19 patients (5-7). Since neutrophilia predicts poor outcomes in patients with , we propose that the change in neutrophil counts in peripheral blood or tissues may be closely associated with pathological injury in COVID-19 patients. We demonstrated here that the dynamics of neutrophil counts in COVID-19 patients during hospitalization exhibited the same trend as the corresponding lung injury. NETs, as confirmed contributors to pathological inflammation of pneumonia, can damage tissues by killing epithelial and endothelial cells (16, 17) of pulmonary tissue in infection and sterile disease. Recently, two elevated NETs markers have been observed in serum from COVID-19 patients, which suggests that neutrophilia and excessive NETs may contribute to cytokine release and respiratory failure in COVID19 patients (10) . However, evidence is still lacking regarding NETosis in lungs. We analyzed the differentially expressed genes in lung tissue and BALF samples from COVID-19 patient in comparison to healthy controls. Among all up-regulated genes in neutrophil modules in COVID-19 patients, we found 17 genes derived from the neutrophil activation pathway were NETs associated genes. Thus, NETs may be activated in the lung of COVID-19 patients. It is also poorly understood how NETosis induces the cytokine storm or modulates the host immune response. Our STRING analysis suggests that NETs associated genes could interact with T, NK, and B cells through regulation of LGALS9, CEACAM1, and TNFSF13B expressions, respectively. We suspect that the progression of lesions in COVID-19 patients may be induced by NETs as well as NETs-T/NK/B cell interactions.

In conclusion, the clear trend of lung injury in accord with the trend of increasing neutrophils was quantified by Bayesian inference analysis in COVID-19 patients. The transcriptome signature of immune cells also indicated elevated neutrophil markers in the lung and BALF samples of COVID-19 patients. Frontiers in Immunology | www.frontiersin.org Importantly, among the excessive neutrophil activated genes, 17 were NETs associated genes and these genes interacted with T cells and NK cells through negative regulation. Therefore, we posit that NETosis in lung tissue leads to an insufficient anti-viral response in COVID-19 patients. We hope that future studies will investigate the predictive power of circulating NETs in well-phenotyped longitudinal cohorts.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

The studies involving human participants were reviewed and approved by the Ethics Committee of the Fifth People's Hospital, Wuxi (No. 2020-006-1). The patients provided their written informed consent to participate in this study.

JW, YQ, and QL conceived and designed the experiments. QL, JW, DH, and ML drafted and revised the manuscript. YY, YZ, and XL carried out the data collection. JW, DH, and YC carried out the data analysis and interpretation. DH, YQ, ML, and LH contributed reagents, materials, and analysis tools. All authors contributed to the article and approved the submitted version.

This work was supported by the foundation of Wuxi Medical Development Discipline for Infectious Disease (FZXK006) and Wuxi Young Medical Talents (QNRC072), Health and Science Bureau of Wuxi (MS201731, CSE31N1712, Q201743). The funding source was not involved in the study design; in the collection, analysis, and interpretation of data, in the writing of the report, and in the decision to submit the paper for publication. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.

WHO. Coronavirus Disease (COVID-2019) Situation Reports. Available online at

Clinical Management of Severe Acute Respiratory Infection When Novel Coronavirus (nCoV) Infection Is Suspected. Geneva: World Health Organization (2020)

Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan

COVID-19: consider cytokine storm syndromes and immunosuppression

Vander Heide RS. Pulmonary and cardiac pathology in Covid-19: the first autopsy series from New Orleans

A pathological report of three COVID-19 cases by minimally invasive autopsies

Targeting potential drivers of COVID-19: neutrophil extracellular traps

Predictors of mortality for patients with COVID-19 pneumonia caused by SARS-CoV-2: a prospective cohort study

Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients

Neutrophil extracellular traps (NETs) as markers of disease severity in COVID-19

Neutrophil extracellular traps are indirectly triggered by lipopolysaccharide and contribute to acute lung injury

NET balancing: a problem in inflammatory lung diseases

Neutrophil extracellular traps in immunity and disease

Neutrophil extracellular traps (NETs) are increased in the alveolar spaces of patients with ventilator-associated pneumonia

Plasma myeloperoxidase-conjugated DNA level predicts outcomes and organ dysfunction in patients with septic shock

Neutrophil extracellular traps directly induce epithelial and endothelial cell death: a predominant role of histones

Netting neutrophils induce endothelial damage, infiltrate tissues, and expose immunostimulatory molecules in systemic lupus erythematosus

Neutrophils sense microbe size and selectively release neutrophil extracellular traps in response to large pathogens

Extracellular DNA traps are associated with the pathogenesis of TRALI in humans and mice

Near-optimal probabilistic RNAseq quantification

Differential analysis of RNA-seq incorporating quantification uncertainty

Pretty Heatmaps. R package v. 16. R Foundation for Statistical Computing

Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression

Metascape provides a biologist-oriented resource for the analysis of systemslevel datasets

Gene set enrichment analysis: a knowledgebased approach for interpreting genome-wide expression profiles

STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

eXamine: exploring annotated modules in networks

FactoMineR: an R package for multivariate analysis

Elegant Graphics for Data Analysis

Bayesian Applied Regression Modeling via Stan

pROC: an open-source package for R and S+ to analyze and compare ROC curves

Stan: a probabilistic programming language

Rac2 is required for the formation of neutrophil extracellular traps

Haploinsufficiency of NADPH oxidase subunit neutrophil cytosolic factor 2 is sufficient to accelerate full-blown lupus in NZM 2328 mice

TREM-1 promotes survival during Klebsiella pneumoniae liver abscess in mice

Neutrophil extracellular traps that are not degraded in systemic lupus erythematosus activate complement exacerbating the disease

Neutrophil extracellular traps protein composition is specific for patients with Lupus nephritis and includes methyl-oxidized αenolase (methionine sulfoxide 93)

Galectin-9 is a possible promoter of immunopathology in rheumatoid arthritis by activation of peptidyl arginine deiminase 4 (PAD-4) in granulocytes

Disruption of neutrophil extracellular traps (NETs) links mechanical strain to post-traumatic inflammation

Cystic fibrosis sputum DNA has NETosis characteristics and neutrophil extracellular trap release is regulated by macrophage migration-inhibitory factor

Characterization of neutrophil function in Papillon-Lefèvre syndrome

Src family kinases and Syk are required for neutrophil extracellular trap formation in response to β-glucan particles

Protective role of mincle in bacterial pneumonia by regulation of neutrophil mediated phagocytosis and extracellular trap formation

Neutrophil extracellular trap-associated CEACAM1 as a putative therapeutic target to prevent metastatic progression of colon carcinoma

Siglec-5 and Siglec-14 are polymorphic paired receptors that modulate neutrophil and amnion signaling responses to group B Streptococcus

Histone hypercitrullination mediates chromatin decondensation and neutrophil extracellular trap formation

We are grateful to the doctors, nurses, disease control workers, and researchers for their fight against COVID-19 under extreme conditions. Some of them have lost their lives in this fight.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu. 2020.02063/full#supplementary-material