Front. Immunol. Frontiers in Immunology Front. Immunol. 1664-3224 Frontiers Media S.A. 10.3389/fimmu.2022.1034420 Immunology Original Research Burn-injured skin is marked by a prolonged local acute inflammatory response of innate immune cells and pro-inflammatory cytokines Mulder Patrick P.G. 1 2 * Vlig Marcel 1 Fasse Esther 2 Stoop Matthea M. 3 Pijpe Anouk 1 3 4 5 van Zuijlen Paul P.M. 3 4 5 6 Joosten Irma 2 Boekema Bouke K.H.L. 1 4 * Koenen Hans J.P.M. 2 1 Preclinical & Clinical Research, Association of Dutch Burn Centres (ADBC), Beverwijk, Netherlands 2 Laboratory of Medical Immunology, Department of Laboratory Medicine, Radboud University Medical Center, Nijmegen, Netherlands 3 Burn Center & Department of Plastic and Reconstructive Surgery, Red Cross Hospital, Beverwijk, Netherlands 4 Department of Plastic Reconstructive and Hand Surgery, Amsterdam UMC Vrije Universiteit Amsterdam, Amsterdam, Netherlands 5 Amsterdam Movement Sciences (AMS) Institute, Amsterdam UMC, Amsterdam, Netherlands 6 Paediatric Surgical Centre, Emma Children’s Hospital, Amsterdam UMC University of Amsterdam, Amsterdam, Netherlands

Edited by: Sandra Sacre, Brighton and Sussex Medical School, United Kingdom

Reviewed by: Marc G. Jeschke, Hamilton Health Sciences, Canada; Lars-Peter Kamolz, Medical University of Graz, Austria

*Correspondence: Patrick P.G. Mulder, pmulder@burns.nl; Bouke K.H.L. Boekema, bboekema@burns.nl

This article was submitted to Inflammation, a section of the journal Frontiers in Immunology

 14 11 2022 2022 13 1034420 01 09 2022 28 10 2022 Copyright © 2022 Mulder, Vlig, Fasse, Stoop, Pijpe, van Zuijlen, Joosten, Boekema and Koenen 2022 Mulder, Vlig, Fasse, Stoop, Pijpe, van Zuijlen, Joosten, Boekema and Koenen

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The systemic and local immune response in burn patients is often extreme and derailed. As excessive inflammation can damage healthy tissues and slow down the healing process, modulation of inflammatory responses could limit complications and improve recovery. Due to its complexity, more detailed information on the immune effects of thermal injury is needed to improve patient outcomes. We therefore characterized and quantified subsets of immune cells and mediators present in human burn wound tissue (eschar), sampled at various time points. This study shows that after burn injury, the number of immune cells were persistently increased, unlike the normal wound healing process. There was an immediate, strong increase in neutrophils and a moderate increase in monocytes/macrophages and lymphocytes, especially in the second and third week post burn. The percentage of classical (CD14highCD16-) monocytes/macrophages demonstrated a steady decrease over time, whereas the proportion of intermediate (CD14highCD16+) monocytes/macrophages slowly increased. The absolute numbers of T cells, NK cells and B cells increased up to week 3, while the fraction of γδ T cells was increased only in week 1. Secretome profiling revealed high levels of chemokines and an overall pro-inflammatory cytokine milieu in burn tissue. The local burn immune response shows similarities to the systemic immune reaction, but differs in neutrophil maturity and lymphocyte composition. Altogether, the neutrophil surges, high levels of pro-inflammatory cytokines and limited immunosuppression might be key factors that prolong the inflammation phase and delay the wound healing process in burns.

 Graphical Abstract

Graphical overview of the study. Top part shows the investigation of burn wound tissue using flow cytometry, immunohistochemistry and cytokine analysis. Bottom part shows the relative response of immune cells during normal wound healing (as described in literature) and after burn injury (as shown by the study results) (1).

 cell isolation immune cells flow cytometry neutrophils lymphocytes macrophages inflammation burn wound tissue WO/17.108 Nederlandse Brandwonden Stichting10.13039/501100003503 Introduction

Burn injury is often accompanied by an extensive, derailed immune response in both burn wound tissue and peripheral blood (2, 3). Regardless of infection, burn patients generally show signs of systemic inflammation caused by high levels of cytokines and danger signals that originate from damaged tissue (4, 5). Necrotized and inflamed tissue stimulates the immune system to recruit acute phase immune cells to the affected site (3, 6, 7). Fibroblasts and keratinocytes surrounding the wound site and infiltrating leukocytes release a storm of cytokines, chemokines and growth factors that initiate the inflammation phase (8).

Typically during wound healing, neutrophils and macrophages with a pro-inflammatory (i.e. M1) phenotype will migrate into the wounded skin to remove debris and prevent bacterial colonization (9). Within days, wound neutrophils will disappear through apoptosis and macrophages will differentiate into a state that supports wound healing (i.e. M2 phenotype) (10). Generally within week after injury, lymphocytes will infiltrate the wound site to orchestrate tailored pathogen-eliminating and immune cell regulating responses (11). The reduction, transition and control of immune cells are crucial for dampening of the inflammatory response and for the establishment of a healthy wound healing process. After burn injury however, the immune system can be overactive and is then likely to cause damage to surrounding tissues, delay wound healing and contribute to the severity of scarring (3, 7).

Burn patients who experience persistent inflammation might benefit from immune suppressive treatment, however at the same time they are at risk of contracting infections such as pneumonia or cellulitis, caused by opportunistic bacteria (12). Therefore, innovative and precise interventions that modulate the immune response could be crucial in the relief of secondary illnesses while improving wound healing and preventing infection. Still, there is only little information on the immune response after burn injury and how exactly it differs from normal wound healing, mainly due to its complexity and variation among cases (e.g. burn size, depth and cause) and burn patients (e.g. age, sex and co-morbidities) (13). Moreover, present evidence on the processes that underlie burn injury originates mostly from animal research (14), which is only partially translatable to the human situation (15). We previously showed that in blood from severely burned patients, there was an extreme increase in innate immune cells and pro-inflammatory cytokines (8). In this longitudinal study, we investigated immune cells and soluble factors present in burn wound tissue (eschar) that was surgically debrided as part of standard treatment (16). A better understanding of the immune reactions to burn injury will facilitate the design of improved and more targeted treatment approaches for trauma-induced immune dysfunction.

 Materials and methods Sample collection

Burn wound tissue (eschar) from patients of all ages and thermal burn causes who underwent eschar debridement as part of their treatment at the Burn Center of the Red Cross Hospital in Beverwijk, the Netherlands. Healthy skin samples were obtained from adult patients who underwent cosmetic surgery (abdominoplasty or elective) at the Department of Plastic and Reconstructive Surgery of the Red Cross Hospital. Tissue samples were collected in the period between February 2019 and December 2021. Consent for the use of residual samples was received through the opt-out protocol of the Red Cross Hospital, in accordance with the national guidelines (https://www.coreon.org/). Subjects were informed of this procedure and were able to withdraw at any point. After surgical removal, tissue samples were stored in RPMI 1640 (Gibco, Paisley, UK) containing 1% penicillin and streptomycin (Gibco) as soon as possible to increase cell survival (17, 18). Samples were stored overnight at 4 °C and processed the following morning. Subject and sample characteristics are shown in Supplementary Table 1 .

 Single cell isolation

This protocol was based on the immune cell isolation procedure from He et al. (19). Biopsies were taken from viable areas of the burn tissues, i.e. white or red areas with bleeding spots and not blackened or leathery areas. Approximately 600 mg of tissue was used per cell isolation for flow cytometry (FCM). Tissue samples were cut into smaller pieces and subsequently divided over 2 C-tubes (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany) containing 5 mL of RPMI 1640 containing 1% penicillin and streptomycin. C-tubes were placed on a tissue dissociator (gentleMACS, Miltenyi Biotec) and program “B” was run. Hundred-fifty µL of 80 mg/mL collagenase I (Merck, St. Louis, MO, USA) in PBS (Gibco) was added and the sample was incubated for 1 h in a shaking water bath at 37 °C. After incubation, the C-tube was placed on the tissue dissociator to run program “B”. Samples were passed through a 500 µm and 40 µm cell strainer (pluriSelect, Leipzig, Germany) to obtain a single cell suspension. Suspensions were centrifuged for 10 min at 450 × g, and supernatant was discarded. The cell pellet was resuspended in erythrocyte lysis buffer (1.5 mM NH4Cl, 0.1 mM NaHCO3 and 0.01 mM EDTA in demineralized water) for 10 min at room temperature. Twenty mL of FCM buffer (PBS containing 1% BSA, 0.05% natrium-azide and 1 mM EDTA) was added and the suspension was centrifuged for 10 min at 450 × g. The pellet was resuspended in 5 mL of FACS buffer and cells were counted on the flow cytometer (MACS Quant Analyzer 10, Miltenyi Biotec GmbH, Bergisch Gladbach, Germany).

 Supervised flow cytometry

From the single cells suspensions approximately 2.5 × 105 cells were used per staining panel. Antibodies used for FCM are displayed in Supplementary Table 2 . A solution of 7-AAD (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany) or propidium iodide (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany) were used to calculate viability of cells. Stained cell samples were acquired on the MACS Quant Analyzer 10 and analysis was performed using FlowLogic (Inivai Technologies, Victoria, Australia). Gating strategy is shown in Supplementary Figure 1 . Data was visualized using Graphpad version 5.01 (PRISM, La Jolla, USA) and R (ggplot package).

 Unsupervised flow cytometry analysis

Lymphocyte (panel 1), T cell (panel 2), neutrophil (panel 3) or monocytic (panel 4) populations were gated based on FSC/SSC, CD45, CD3, CD15 in MACSQuantify 2.13.3 software (Miltenyi Biotec). Data of these sole populations were uploaded to Cytobank (20) to create Flow Self-Organizing Map (FlowSOM) clusters.

 Immunohistochemistry

Kryofix (50% ethanol, 3% PEG300)-fixed paraffin-embedded 5 µm thick sections were used after deparaffinization and rehydration. Endogenous peroxidase was blocked in 1% hydrogen peroxide for 15 min at room temperature. Next, antigen retrieval for different antigens was performed. The blocking step was performed using 5% normal goat serum (Merck) diluted in PBS + 1% bovine serum albumin (BSA). Tissue sections were then incubated with the primary antibodies ( Supplementary Table 3 ) for 1 h at RT followed by incubation with a poly-HRP-goat-anti-mouse or rabbit secondary antibody (BrightVision, VWR) for 30 min at RT. Detection of the target protein was established using 3,3′-Diaminobenzidine (BrightDAB, VWR). After immunohistochemical (IHC) DAB staining was successful, sections were counterstained with hematoxylin, dehydrated and mounted with Eukitt Mounting Medium (Merck). Percentage of MPO, CD3 or CD68 positive area was calculated using NIS Elements (Nikon Instruments Europe B.V.) and based on 3 images from representative tissue sections.

 Multiplex imaging and analysis

Formalin-fixed and paraffin-embedded 5 µm sections were deparaffinized using xylene and rehydrated with ethanol and distilled water. Antigen retrieval was performed by boiling in TRIS-borate-EDTA buffer for 10 min. A multiplex staining for the detection of neutrophils and lymphocytes was performed performed using the protocol described by Rodriguez et al. (21).

 Immunoassay of tissue homogenates

Frozen tissue samples of approximately 60 mg were thawed, minced into smaller pieces and further dissociated in M-tubes (Miltenyi Biotec) by adding lysis buffer (PBS containing 0.01 mM EDTA and protease inhibitor (1 tablet per 10 mL; Pierce, Thermo Scientific)) and running program “Protein_01” on a gentleMACS (Miltenyi Biotec). Debris was removed from the samples using a filter plate (Multiscreen, Merck) and samples were diluted to concentration of 12 mg tissue/mL. Luminex assay was performed according to the manufacturer’s instructions (Merck KGaA). The following assay kits were used: HCYTA-60K, TGFBMAG-64K, HCYP2MAG-62K and HTH17MAG-14K. In short, 25 µL of tissue homogenate was used to determine the concentrations of 37 soluble factors, namely MCP-1 (CCL2), MIP-1α (CCL3), MIP-1β (CCL4), MIP-3α (CCL20), GROα (CXCL1), IP-10 (CXCL10), IFN-α2, IFN-γ, TNF-α, TGF-β1, TGF-β2, TGF-β3, CTACK (CCL27), RANTES (CCL5), IL-1α, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-8 (CXCL8), IL-9, IL-10, IL-12p40, IL-12p70, IL-13, IL-17A (CTLA-8), IL-17F, IL-18, IL-21, IL-22, IL-23, IL-33 (NF-HEV), GM-CSF, PDGF-AA, PDGF-AB/BB and VEGF-A. For TGF-β1,2,3 samples were acid-treated prior to the assay, according to the manufacturer’s instructions. Mean fluorescence intensity of samples was measured with a Flexmap 3D System (Luminex Corp, Austin, USA) and concentrations were calculated using Bio-Plex Manager Software (Bio-Rad Laboratories, Veenendaal, The Netherlands). Values below the minimum of the standard were based on extrapolation of the standard curve by the software.

 Statistical analyses

Distribution of the data was checked for normality using the Shapiro Wilk test. For the FCM and IHC data, differences between burn tissue and healthy skin, and between burn tissues of different time intervals after injury (PBW 1, 2, 3 and 4) were explored using the Mann Whitney U test in Graphpad version 5.01 (PRISM, La Jolla, USA). Only statistically significant differences are shown and are indicated by black asterisks (* = p < 0.05; ** = p < 0.01; *** = p < 0.001). The data was visualized using Graphpad version 5.01 (PRISM, La Jolla, USA). Levels of the soluble immune factors in burn tissue were transformed to fold change differences compared to the levels in healthy skin. P values of differences between burn tissue and healthy skin were determined using Mann Whitney U tests. We considered a p value of < 0.05 to be statistically significant. Volcano plots were created using “EnhancedVolcano” version 1.6.0 package in R version 3.6.2.

 Results Burn injury is followed by a strong local increase in granulocytes and moderate increase in monocytes and lymphocytes

Local immune effects of burn trauma were investigated in burn tissue that was debrided during routine surgical procedures (subject and sample characteristics shown in Supplementary Table 1 ). We selected viable sections of tissue biopsies and neglected necrotized or blackened segments to ensure the isolation of viable cells. CD45 immunohistochemical (IHC) staining showed an extreme infiltration of leukocytes (CD45+ cells) in burn tissue ( Figure 1A ). The majority of leukocytes were viable after isolation from healthy skin (90.3% ± 6.6) and burn tissue (89.1% ± 10.5) ( Supplementary Figure 2A ). Flow cytometry (FCM)-based quantification revealed that the increase in leukocyte numbers was most abundant at post burn week (PBW) 2-3 ( Figure 1B ). As a result, the percentage of CD45- cells, which include fibroblasts, keratinocytes and endothelial cells, was lower in burn tissue from PBW 2-3 than in healthy skin ( Supplementary Figure 2B ). The leukocytes isolated from healthy skin consisted of approximately 25% granulocytes, 55% monocytic cells (monocytes and macrophages) and 20% lymphocytes ( Figure 1C ). In burn tissue from PBW 1, there were 52% granulocytes, while for the proportion of monocytic cells was 29%. The lymphocytes fraction in burn tissue was similar to healthy skin (19%). In burn tissue from PBW 2-4, the portion of granulocytes was still enlarged (55-62%), while the fraction of monocytic cells decreased further to 13-16% and the lymphocyte fraction increased (24-31%). During PBW 1-3, absolute number of granulocytes, monocytic cells and lymphocytes rose and declined only at PBW 4 ( Figure 1D ). Multiplex spatial phenotyping of healthy skin and burn tissue sections using CD3 and CD15 revealed dense areas populated with granulocytes and T cells in burn tissue ( Figures 1E, F ).

High number of immune cells infiltrate the skin as response to burn injury. (A) CD45 immunohistochemical DAB staining of a representative section of healthy skin and burn tissue (from 15 days post burn) (black scale bar = 100 µm). Flow cytometry-based quantification of: (B) Absolute number of leukocytes per mg tissue (based on side scatter and CD45); (C) Proportion of granulocytes (Gran), monocytic cells (Mon) and lymphocytes (Lym) in tissue (based on side scatter and CD45); (D) Absolute numbers of granulocytes, monocytic cells and lymphocytes per mg tissue (based on side scatter and CD45). Microscopic image of multiplex DAPI, CD15 (granulocytes) and CD3 (T cells) staining of a representative: (E) Healthy skin sample; (F) Burn tissue sample (from 25 days post burn), shown separately and as composite (black scale bar = 100 µm). P values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01; ***p < 0.001.

 Granulocytes in burn tissue consist mainly of activated mature neutrophils

IHC analysis of myeloperoxidase (MPO) expression, an enzyme abundantly present in azurophilic granules of neutrophils (22), showed an immediate increase in neutrophil numbers in burn tissue at PBW 1, and an even larger increase from PBW 2 onward ( Figures 2A, B ). This was confirmed by FCM analysis of neutrophils (CD15+CD16+ granulocytes) ( Figure 2C ). Eosinophils (CD9+CD15+CD16- granulocytes) were increased at PBW 2-3, but to a lesser extent ( Supplementary Figure 2C ). In both healthy skin and burn tissue neutrophils were almost exclusively CD10+, a marker that is associated with neutrophil maturation (23) ( Figure 2D ). Only in burn tissue from PBW 1 there was a slight increase in immature (CD10-) neutrophils. Activation markers CD11b and CD66b were upregulated in neutrophils at PBW 2-3 ( Figures 2E, F ). Self-organizing map clustering of flow data (FlowSOM) using Cytobank displayed cell populations (nodes) and clusters based on marker expression in an unsupervised manner ( Figure 2G ). This analysis highlights some of the burn-specific changes that occur in wound neutrophils. Burn injury caused significant differences in the percentage of neutrophils per cluster ( Figure 2H ). CD11blowCD14+CD66b- neutrophils (cluster 1) were decreased early after burn injury, while CD11b+CD66blow neutrophils (cluster 2) were increased. Although CD11bhighCD66b+ neutrophils (cluster 3) seemed more represented in burn tissue than in healthy skin, no significant difference was found. A small population of CD16low neutrophils (cluster 4) was significantly increased at PBW 1 and the percentage of CD16lowCD14+ neutrophils (cluster 5) was significantly increased at PBW 4.

Local neutrophil response to burn injury. (A) Myeloperoxidase (MPO) immunohistochemical DAB staining of a representative section of healthy skin and burn tissue (from 15 days post burn) (black scale bar = 100 µm). (B) MPO+ area of tissue sections. Flow cytometry-based quantification of: (C) Absolute number of neutrophils (CD15+CD16+ granulocytes) per mg tissue; (D) Percentage of CD10+ (mature) neutrophils (CD15+CD16+ granulocytes); (E) MFI of CD11b on neutrophils (CD15+CD16+ granulocytes) in tissue; (F) MFI of CD66b on neutrophils (CD15+CD16+ granulocytes) in tissue. (G) Unsupervised clustering of neutrophil (CD15+CD16+ granulocytes) flow data from healthy skin and burn tissue, 5 clusters are highlighted. Node size represents relative size of population and node diagram shows expression level of markers. (H) Percentage of neutrophils within each cluster. Error bars in H show boxplot, p values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01; ***p < 0.001.

 Burn injury increases macrophage numbers and affects differentiation

Macrophage differentiation was assessed by analyzing the CD14 and CD16 expression of the monocytic cell population (among which are monocytes and macrophages) using flow cytometry. In both healthy skin and burn tissue the majority of these cells expressed a classical phenotype (CD14highCD16-) ( Figures 3A, B ). In burn tissue from PBW 3, the proportion of classical monocytic cells decreased while the proportion of intermediate (CD14highCD16+) monocytic cells increased. Next, we analyzed the macrophages (CD68+ cells) within the monocytic cell population and found a steady increase in macrophages over time after burn injury ( Figures 3C–E ). By both IHC and FCM we could detect a significant increase in macrophages at PBW 3. Macrophage phenotype was further investigated by analyzing CD40 and CD80 expression (indicative for pro-inflammatory phenotype) and CD163 and CD206 expression (hallmarks for pro-healing) ( Supplementary Figure 2D ). The only significant difference we observed was a reduction of CD40+ macrophages at PBW 3 ( Supplementary Figure 2D ). Using FlowSOM analysis of the FCM data, we identified macrophage subtypes with different expression patterns: CD163- macrophages (cluster 1), CD163+ macrophages with a low or moderate expression of CD40, CD80 and CD206 (cluster 2) and CD163+ macrophages with a moderate to high expression of CD40, CD80 and CD206 (cluster 3) ( Figure 3F ). A significant increase in macrophages in cluster 2 was observed in burn tissue at PBW 1 and 3 ( Figure 3G ). Overall, this analysis demonstrated that burn injury increased the number of macrophages and changed their composition.

Local macrophage response to burn injury. (A) Flow cytometry gating strategy for detection of differentiation stages of monocytic cells (classical, intermediate or non-classical, as based on CD14 and CD16). (B) Flow cytometry-based quantification of percentage of monocytic cells within classical, intermediate, non-classical gates. (C) CD68 immunohistochemical DAB staining of a representative section of healthy skin and burn tissue (from 15 days post burn) (black scale bar = 100 µm). (D) CD68+ area of tissue sections. (E) Flow cytometry-based quantification of absolute number of macrophages (CD68+ monocytic cells) per mg tissue. (F) Unsupervised clustering of macrophages (CD68+ monocytic cells) in healthy skin and burn tissue, 3 clusters are highlighted. Node size represents relative size of population and node diagram shows expression level of markers. (G) Percentage of macrophages (CD68+ monocytic cells) within each cluster. Error bars in G show boxplot, p values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01.

 Burn injury causes shifts in the lymphocyte composition and increases total T cells at PBW 2-3

T cell (CD3+ lymphocyte) numbers rose significantly at PBW 2-3 ( Figure 4A ), in line with the total lymphocyte increase ( Figure 1D ). A shift towards more CD4+ T cells was detected in burn tissue compared to healthy skin and were highest in burn tissue from PBW 3 as the CD4/CD8 T cell ratio (CD3+CD4+/CD3+CD4- ratio) was higher in burn tissue than in healthy skin ( Figure 4B and Supplementary Figure 2E ). An increase in the proportion of γδ T cells (CD3+CD4-γδTCR+ lymphocytes) was found only at PBW 1 ( Figure 4C ), indicating a fast response of γδ T cells after burn injury. The absolute number of γδ T cells steadily increased over time after burn injury ( Supplementary Figure 2F ). The shift towards a higher abundance of γδ T cells at PBW 1 was confirmed by mapping flow cytometry data of T cells using FlowSOM (clusters 3 and 4; Figures 4D, E ) and shows that the majority of the γδ T cells was CD25+, which is a prominent marker for cellular activation (24). At PBW 1 there was a relative decrease of T cells with a regulatory phenotype (CD25+CD127-; cluster 1). We did not observe considerable alterations in the cluster containing CD3+CD4- T cells (cluster 5).

Local T cell response to burn injury. Flow cytometry-based quantification of: (A) Absolute number of T cells (CD3+ lymphocytes) per mg tissue; (B) CD4+/CD4- T cell (CD3+ lymphocytes) ratio in tissue; (C) Percentage of T cells (CD3+ lymphocytes) that are γδ T cells (γδTCR+CD4- T cells). (D) Unsupervised clustering of T cells (CD3+ lymphocytes) in healthy skin and burn tissue, 5 clusters are highlighted. Node size represents relative size of population and node diagram shows expression level of markers. (E) Percentage of T cells within each cluster. Error bars in E show boxplot, p values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01; ***p < 0.001.

 Absolute number of NK and B cells increase after burn injury

FCM analysis showed that the absolute number of NK cells (CD56+ lymphocytes) was higher in burn tissue from PBW 2-3 ( Figure 5A ). Relative to total leukocyte numbers, NK cells were significantly reduced in burn tissue from PBW 1 and normalized afterwards ( Supplementary Figure 2F ). In both healthy skin and burn tissue, the majority of the NK cells was CD16- ( Figure 5B ), which is opposed to the NK cell composition in peripheral blood where approximately 90% of the NK cells are CD16+ (25). Differences in CD16 expression of the NK cells were not observed between healthy skin and burn tissue or between time points. The absolute number of B cells (CD19+ lymphocytes) were higher in burn tissue from PBW 3 ( Figure 5C ), while the proportion of B cells within the leukocyte population in burn tissue was similar to that of healthy skin ( Supplementary Figure 2G ). We identified 4 clusters using FlowSOM analysis of the FCM data: CD56+CD16+ NK cells, CD56+CD16- NK cells, CD9+CD56+ B cells and CD9lowCD56- B cells ( Figure 5D ). Clustering analysis showed a clear shift towards more CD9lowCD56- B cells in burn tissue but no significant differences were detected over time ( Figure 5E ).

Local NK and B cell response to burn injury. Flow cytometry-based quantification of: (A) Absolute number of NK cells (CD56+ lymphocytes) per mg tissue; (B) Percentage of NK cells that are CD16- and CD16+; (C) Absolute number of B cells (CD19+ lymphocytes) per mg tissue. (D) Unsupervised clustering of NK and B cells in healthy skin and burn tissue, 4 clusters are highlighted. Node size represents relative size of population and node diagram shows expression level of markers. (E) Percentage of NK or B cells within each cluster. Error bars in E show boxplot, p values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01; ***p < 0.001.

 Immune cell infiltration coincides with high levels of cytokines, chemokines and growth factors

The concentrations of 37 soluble immune factors were determined in homogenates of burn tissue using Luminex immunoassay (raw data is presented in Supplementary Figure 3 ). Figure 6 shows an overview of these results using volcano plots and heatmaps at 4 time intervals after burn injury. In burn tissue there was an extremely high expression of IL-6, IL-1β, IFN-γ and TNF-α compared to healthy skin. The levels of these factors were persistently high, but for IL-6 and IFN- γ the levels declined at the later time intervals. Interestingly, increased levels of IL-12p40 and IL-5 were found only late after burn injury (PBW 3-4). As compared to healthy skin, a decrease was found for IL-1 family members IL-1α, IL-33 and IL-18. This is opposed to the level of IL-1β, which is also an IL-1 cytokine. The levels of IL-1α, IL-33 and IL-18 somewhat normalized at PBW 4 to the levels found in healthy skin. Chemokines MCP-1, IL-8, GROα, MIP-1α, MIP-1β, RANTES, IP-10 in burn tissue were increased at all analyzed time intervals, while the levels of T cell attracting chemokines CTACK and MIP-3α were decreased at PBW 1-4 and PBW 4, respectively. Among the growth factors, an increase in VEGF-A and TGF-β1 levels was found at PBW 1-4. From the growth factors, the level of GM-CSF was increased found at PBW 1-3, PDGF-AA at PBW 2-4 and PDGF-AB/BB and TGF-β2 only at PBW 3.

Expression of cytokines, chemokines and growth factors in burn tissue. (A) Volcano plot visualization of the expression of soluble factors in burn tissue from 1, 2, 3 and 4 weeks post burn injury. Dots represent soluble factors in burn tissue as Log2 fold change as found in healthy skin controls. Factors with a statistically significantly different expression (p < 0.05) are labeled (values above the black striped line). (B) Heatmap visualization of fold increase/decrease of soluble factors in burn tissue compared to healthy skin, categorized by cytokines, chemokines and growth factors. P values were calculated using Mann-Whitney U statistical tests.

 Discussion

Next to being a protective, physical barrier, the skin carries out immune surveillance to ensure early and effective defense mechanisms against external threats. Besides fibroblasts and keratinocytes, healthy skin is inhabited mainly by lymphocytes and antigen presenting cells that survey the skin and react to foreign structures and danger signals (26). Here, we aimed to provide detailed insight in the cellular and soluble immune response in burn injured skin during the first four weeks after injury. In this study, we showed that after burn injury, there is a fast, extensive and long-lasting increase in innate immune cells that is present even in burn tissue debrided 3 to 4 weeks after injury. Lymphocytes also rise in numbers, but mainly at PBW 2-4. In addition, the cytokine composition in these burn tissue samples is highly pro-inflammatory and likely continues the attraction and activation of immune cells. Excessive pro-inflammatory immune responses and a lack or delay of anti-inflammatory responses could complicate wound healing and patient recovery. Limitations of this study that should be addressed are minor differences in treatment between patients such as medication and timing of surgery that could have influenced the inflammatory response. In addition, the broad range in subject age, burn cause and TBSA could have increased variation in the responses.

In tissue samples from PBW 1, the proportion of γδ T cells was increased, indicating that γδ T cells could play a role during the early phase of burn-induced response. γδ T cells possess a unique T cell receptor and can, unlike αβ T cells, interact with antigens directly (27). They execute immune surveilling functions and react to damaged cell structures by producing cytokines and chemokines to recruit immune cells (28). Mouse studies have shown that γδ T cells regulate the infiltration of innate immune cells shortly after trauma (29, 30). Our data suggests that next to keratinocytes, fibroblasts, mast cells and platelets (31), γδ T cells could be important inducers of the inflammatory response in humans as well. Within the same timeframe (PBW 1), IL-1β, IL-6, IL-8 (CXCL8), MCP-1 (CCL2), and GRO-α (CXCL1) levels were highly augmented. Others have demonstrated that these cytokines are also elevated in burn wound exudate (32). These factors are known enhancers of the inflammatory response and attract neutrophils and monocytes/macrophages to wound site (33). On the contrary, levels of IL-1α, IL-18 and IL-33 in burn tissue were reduced, especially during at PBW 1. These IL-1 cytokine family members are constitutively produced by keratinocytes to maintain the immune surveillance aspect of the skin (34). Reduction of these factors is presumably caused by extreme loss of keratinocytes due to destruction of the epidermal layer by thermal injury. In burn tissue from PBW 2-4, the levels of IL-1α, IL-18 and IL-33 were returning to the levels in healthy skin, which may be related to the presence of keratinocytes closing the defect. Levels of cytokines, as well as microRNAs (35), could be potential biomarkers to predict disease progression or recovery (36).

The rapid neutrophil response to burn injury is presumably caused by the persisting levels of neutrophil attractants, such as IL-8, MCP-1 and GRO-α. This can also be observed in the circulation of burn patients, where high levels of neutrophils were accompanied by high levels of IL-8 and MCP1, especially early after injury (8). Other studies have also shown that burn tissue contains large numbers of neutrophils in both human (37) and animals (14, 38). The vast majority of neutrophils that infiltrated the wound area were mature, whereas, in peripheral blood from severely burned patients high numbers of immature neutrophils were detected (8). This release of immature neutrophils may well be a compensatory response by the bone marrow (39). Nucleus flexibility and chemotactic activity increases with neutrophil age and could explain the presence of mainly mature neutrophils in burn tissue (40). If only mature neutrophils are able to enter the wound site, immature neutrophils would be trapped in the circulation until they reach maturity. As immature neutrophils are proposedly more active and less predictable in reacting to danger signals (41), they are likely to enhance systemic inflammation, thereby delaying recovery. In burn tissue, we found only a small number of immature neutrophils and only at PBW 1, which could have been released from the blood circulation by capillary leakage caused by the burn injury. Expression of CD11b and CD66b was increased on neutrophils isolated from burn tissue. This highlights the inflammatory state of the infiltrating neutrophils as CD11b and CD66b are important for neutrophil activation, adhesion and migration to inflamed tissue (42, 43). The surges of active neutrophils in the wound could lead to an overproduction of products such as elastase, MPO and ROS which can (further) damage surrounding tissues and organs (44, 45).

Blood monocytes are progenitors of both pro-inflammatory macrophages and wound healing macrophages. Although there is little evidence in this respect, it has been suggested that classical monocytes could be predisposed progenitors to pro-inflammatory macrophages (46), while intermediate and non-classical monocytes are biased progenitors to wound healing macrophages (47, 48). The initial monocyte population in burn tissue consisted mainly of classical monocytes. The relative decrease in classical monocytes in PBW 3 could indicate a relevant shift towards more wound healing macrophages, which is assumed to happen during wound healing (9). In burn tissue, the number of macrophages was increased this population showed a different composition of M1 and M2 markers. CD163+ macrophages with low to moderate expression of CD40, CD80 and CD206 were more abundant in burn tissue. M1 macrophage differentiation factor GM-CSF was increased in burn tissue from PBW 1-3 and mediators that are known to be actively produced by M1 macrophages such as TNF-α, IFN-γ, IL-1β, IL-6, IL-8 and MCP-1 (CCL2) (9, 49), were all increased in these burn tissues. While typical M2 macrophage factors like IL-4, IL-10, IL-13 were unaffected, the levels of TGF-β1 and VEGF-α, which are also described as M2 mediators (48), were increased in burn tissue. Altogether, the monocyte/macrophage composition and cytokine environment possibly supports the generation of macrophages with a pro-inflammatory phenotype. Timely transition towards more suppressive, regenerative macrophages is however essential for a healthy healing process, as these cells support fibroblasts in the formation of collagen and enhance re-epithelization (9, 48). Due to the active, continuing inflammation and the presence of danger signals from tissue damage, macrophage transition might be delayed or insufficient, although more research is required to elucidate this.

Immunosuppression from the adaptive arm of the immune system is essential to create an environment in which fibroblasts and keratinocytes can repair the damaged skin (50). Here, we revealed that lymphocyte numbers, including T cells, NK cells and B cells, were increased at PBW 2-3, which is relatively late after injury (11). This coincided with a high levels of chemokines MIP-1α (CCL3), MIP-1 β (CCL4) and RANTES (CCL5), which are known to attract lymphocytes to injured skin (11). Particularly CCL3, CCL4 and CCL5 are involved in the activation of NK cells (51) and could lead to increased cytokine release by NK cells in burn tissue. Information on the response of NK cells and B cells in burn tissue is very limited at this moment. We here showed that after burn injury there is an increased number of NK cells and B cells in burn tissue, however, functional assays are needed in order to speculate about their behavior and involvement in the burn immune response. The levels of CCL3, CCL4, CCL5, IFN-γ, TNF-α and IP-10 (CXCL10) were associated with the number of T cells at PBW 3. IP-10 promotes chemotaxis and inflammation and is likely induced by IFN-γ. Peters et al. previously described an interplay of keratinocytes and T cells and showed that IP-10 is actively produced by keratinocytes in co-culture, even with relatively low numbers of keratinocytes (52). This interplay is presumably also active during burn wound healing by residual, surrounding or re-epithelializing keratinocytes. Cytokines with anti-inflammatory properties such as IL-4, IL-10 and IL-13 were not detected in these burn tissue samples. Altogether, the soluble factors in burn tissue are likely to support Th1 response, resulting in more attraction of leukocytes to the wound site, while control or suppression of inflammation appears to be limited.

In this study, we showed that after burn injury, the numbers of immune cells were persistently elevated, while during normal wound healing neutrophils disappear within days and lymphocytes counts start to increase in the first week (10, 53, 54). Burn injury often leads to a prolonged hyperinflammatory state (3, 55) and treatment of burn wounds is therefore a difficult and time-consuming process. Damage to the skin is a trigger for the immune system to recruit immune cells en masse and replenish these immune cells in the blood from the bone marrow. Ancillary damage and chemokine production by immune cells and stressed skin cells will trigger the immune system to react, thereby creating a vicious circle of prolonged inflammation in both the skin and in the blood. Therapy is often empiric due to the large diversity among patients and their injuries, e.g. burn type, size, depth and location. In the present study, there was no indication that burn size or burn cause (water versus flame) affected cellular or soluble inflammatory markers (data not shown). Excessive and persistent inflammation is also among the causes of long-term complications such as the formation of hypertrophic scars (7). On top of that, there is a risk of contracting an infection and the activity of the immune system is unpredictable. In clinical practice, patients with burns larger than 15% TBSA are hypermetabolic and often develop SIRS or organ insufficiency. Hence debridement of burn tissue is important to reduce inflammation and promote wound healing while also preventing further tissue necrosis and cellulitis. Possibly, early debridement of burn tissue (noted as post-burn days 2 through 12) or impediment of pro-inflammatory cytokines such as IL-6 might remove inflammatory triggers at an early stage and avoid secondary damage (56, 57). Resolution of excessive inflammation using immune suppressants could increase the patients’ recovery rate, but might increase the risk for infection. Moreover, it can be very difficult to discriminate burn-induced SIRS from sepsis. Our analysis of the local immune reactions to burn injury aids in improving our understanding of burn-induced inflammation. This knowledge is needed to design more sophisticated and effective ways to diagnose and treat immune dysfunction and hyperactive inflammation. Immune modulating treatment targeting the disturbed immune processes will improve patients’ overall health recovery time and scar quality.

In conclusion, through the characterization of immune cell subsets isolated from human burn tissue we demonstrated that burn injury induced a local persistent surge of pro-inflammatory immune cells and cytokines, while immunosuppression appeared to be limited. These burn-induced immune reactions might be key factors that extent the inflammation phase and thereby obstruct the wound healing process in burn injury.

 Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

 Ethics statement

Tissue samples were obtained from planned surgeries that were part of routine patient therapy. Data and material was only used after consent from subjects (or legal representatives) using the opt-out procedure of the Red Cross Hospital in Beverwijk. This procedure is in accordance with the national guidelines (https://www.coreon.org/) and institutional guidelines of the Red Cross Hospital in Beverwijk in the Netherlands. No animals were used in these experiments.

 Author contributions

Conceptualization: PM, MV, IJ, BB, and HK. Methodology: PM, MV, IJ, BB, HK, and AP. Investigation: PM, MV, and EF. Resources: MS, AP, and PVZ. Formal analysis: PM, MV, BB, and HK. Visualization: PM and MV. Supervision: IJ, BB, and HK. Writing – original draft: PM. Writing – review & editing: MV, EF, MS, AP, PVZ, IJ, BB, and HK. All authors have read and approved the submission of this manuscript.

 Funding

Dutch Burns Foundation awarded grant WO/17.108 to BB to execute this project.

 Acknowledgments

We would like to thank Xuehui He for helping with the cell isolation protocol, Bram van Cranenbroek for his assistance on the immunoassay, Mark Gorris for helping with the multiplex imaging and Evi Warmerdam, Rosa Rentenaar and Myrthe Witbaard for their technical assistance during experiments. We are grateful for the work and participation of all involved physicians, surgeons, nurses and patients of the Burn Center and the Department of Plastic and Reconstructive Surgery of the Red Cross Hospital.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

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

Gating strategy of flow cytometric analysis. At the top of the plots the panel numbers are shown for which the gating was performed (P1, 2, 3 and 4).

Cell counts in burn tissue. Flow cytometry-based quantification of: (A) Percentage of leukocytes that stained negative for 7-AAD or propidium iodide (viable cells); (B) Percentage of isolated cells that is CD45- (fibroblasts, keratinocytes, endothelial cells, others); (C) Number of eosinophils (CD9+CD16- granulocytes) per mg tissue; (D) Percentage of macrophages that is positive for CD40, CD80, CD163 or CD206; (E) Percentage of T cells that are CD4- and CD4+. (F) Number of γδ T cells (γδTCR+) per mg tissue; (G) Percentage of NK cells (CD56+ lymphocytes) within leukocyte population; (H) Percentage of B cells (CD19+ lymphocytes) within leukocyte population. P values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01; ***p < 0.001.

Concentrations of soluble factors in burn tissue. Healthy skin was used as controls. Black lines show mean values and the black striped line represents the lowest limit of detection. P values were calculated using Mann-Whitney U statistical tests, significant differences are indicated by black asterisks: *p < 0.05; **p < 0.01; ***p < 0.001.

 References Witte MB Barbul A . General principles of wound healing. Surg Clin North Am (1997) 77:509–28. doi: 10.1016/S0039-6109(05)70566-1 Dahiya P . Burns as a model of SIRS. Front Biosci (2009) 14:4962–7. doi: 10.2741/3580 Bergquist M Hästbacka J Glaumann C Freden F Huss F Lipcsey M . The time-course of the inflammatory response to major burn injury and its relation to organ failure and outcome. Burns (2019) 45:354–63. doi: 10.1016/j.burns.2018.09.001 Wu G Zhuang M Jiang Y Fan J Sun Y Zhou Z . Can systemic inflammatory response syndrome score at admission predict clinical outcome in patients with severe burns? Burns (2019) 45:860–8. doi: 10.1016/j.burns.2018.11.011 Farina JA Rosique MJ Rosique RG . Curbing inflammation in burn patients. Int J Inflam (2013) 2013:1–9. doi: 10.1155/2013/715645 Rani M Nicholson SE Zhang Q Schwacha MG . Damage-associated molecular patterns (DAMPs) released after burn are associated with inflammation and monocyte activation. Burns (2017) 43:297–303. doi: 10.1016/j.burns.2016.10.001 Zhu Z Ding J Tredget EE . The molecular basis of hypertrophic scars. Burn Trauma (2016) 4:2. doi: 10.1186/s41038-015-0026-4 Mulder PPG Vlig M Boekema BKHL Stoop MM Pijpe A van Zuijlen PPM . Persistent systemic inflammation in patients with severe burn injury is accompanied by influx of immature neutrophils and shifts in T cell subsets and cytokine profiles. Front Immunol (2021) 11:621222. doi: 10.3389/fimmu.2020.621222 Kotwal GJ Chien S . Macrophage differentiation in normal and accelerated wound healing. In: Macrophages origin, functions and biointervention (Cham: Springer) (2017) 62:353–64. doi: 10.1007/978-3-319-54090-0_14 Velnar T Bailey T Smrkolj V . The wound healing process: An overview of the cellular and molecular mechanisms. J Int Med Res (2009) 37:1528–42. doi: 10.1177/147323000903700531 Short WD Wang X Keswani SG . The role of T lymphocytes in cutaneous scarring. Adv Wound Care (2022) 11:121–31. doi: 10.1089/wound.2021.0059 Norbury W Herndon DN Tanksley J Jeschke MG Finnerty CC . Infection in burns. Surg Infect (Larchmt) (2016) 17:250–5. doi: 10.1089/sur.2013.134 van Zuijlen PPM Korkmaz HI Sheraton VM Haanstra TM Pijpe A de Vries A . The future of burn care from a complexity science perspective. J Burn Care Res (2022) 43(6):1–10. doi: 10.1093/jbcr/irac029 Mulder PPG Koenen HJPM Vlig M Joosten I de Vries RBM Boekema BKHL . Burn-induced local and systemic immune response: Systematic review and meta-analysis of animal studies. J Invest Dermatol (2022) 142(11):3093–109. doi: 10.1016/j.jid.2022.05.004 Zomer HD Trentin AG . Skin wound healing in humans and mice: Challenges in translational research. J Dermatol Sci (2018) 90:3–12. doi: 10.1016/j.jdermsci.2017.12.009 Orgill DP . Excision and skin grafting of thermal burns. N Engl J Med (2009) 360:893–901. doi: 10.1056/NEJMct0804451 Leelatian N Doxie DB Greenplate AR Sinnaeve J Ihrie RA Irish JM . Preparing viable single cells from human tissue and tumors for cytomic analysis current protocols in molecular biology. Curr Protoc Mol Biol (2018) 118:1–31. doi: 10.1002/cpmb.37.Preparing Boekema BKHL Boekestijn B Breederveld RS . Evaluation of saline, RPMI and DMEM/F12 for storage of split-thickness skin grafts. Burns (2015) 41:848–52. doi: 10.1016/j.burns.2014.10.016 He X de Oliveira VL Keijsers R Joosten I Koenen HJ . Lymphocyte isolation from human skin for phenotypic analysis and &lt;em<Ex vivo&lt;/em< cell culture. J Vis Exp (2016) 110:7–13. doi: 10.3791/52564 Kotecha N Krutzik PO Irish JM . Web-based analysis and publication of flow cytometry experiments. Curr Protoc Cytom (2010) 53:10.17.1–10.17.24. doi: 10.1002/0471142956.cy1017s53 Rodriguez-Rosales YA Langereis JD Gorris MAJ van den Reek JMPA Fasse E Netea MG . Immunomodulatory aged neutrophils are augmented in blood and skin of psoriasis patients. J Allergy Clin Immunol (2021) 148:1030–40. doi: 10.1016/j.jaci.2021.02.041 Metzler KD Goosmann C Lubojemska A Zychlinsky A Papayannopoulos V . Myeloperoxidase-containing complex regulates neutrophil elastase release and actin dynamics during NETosis. Cell Rep (2014) 8:883–96. doi: 10.1016/j.celrep.2014.06.044 Marini O Costa S Bevilacqua D Calzetti F Tamassia N Spina C . Mature CD10+ and immature CD10- neutrophils present in G-CSF-treated donors display opposite effects on T cells. Blood (2017) 129:1343–56. doi: 10.1182/blood-2016-04-713206 Bajnok A Ivanova M Rigó J Toldi G . The distribution of activation markers and selectins on peripheral T lymphocytes in preeclampsia. Mediators Inflammation (2017) 2017:1–7. doi: 10.1155/2017/8045161 Amand M Iserentant G Poli A Sleiman M Fievez V Sanchez IP . Human CD56dimCD16dimCells as an individualized natural killer cell subset. Front Immunol (2017) 8:699. doi: 10.3389/fimmu.2017.00699 Kupper TS Fuhlbrigge RC . Immune surveillance in the skin: Mechanisms and clinical consequences. Nat Rev Immunol (2004) 4:211–22. doi: 10.1038/nri1310 Kim A Lang T Xue M Wijewardana A Jackson C Vandervord J . The role of Th-17 cells and γδ T-cells in modulating the systemic inflammatory response to severe burn injury. Int J Mol Sci (2017) 18:758. doi: 10.3390/ijms18040758 Toth B Alexander M Daniel T Chaudry IH Hubbard WJ Schwacha MG . The role of γδ T cells in the regulation of neutrophil-mediated tissue damage after thermal injury. J Leukoc Biol (2004) 76:545–52. doi: 10.1189/jlb.0404219 Rani M Zhang Q Schwacha MG . Gamma delta T cells regulate wound myeloid cell activity after burn. Shock (2014) 42:133–41. doi: 10.1097/SHK.0000000000000176 Cruz MS Diamond A Russell A Jameson JM . Human αβ and γδ T cells in skin immunity and disease. Front Immunol (2018) 9:1304. doi: 10.3389/fimmu.2018.01304 Menter DG Kopetz S Hawk E Sood AK Loree JM Gresele P . Platelet “first responders” in wound response, cancer, and metastasis. Cancer Metastasis Rev (2017) 36:199–213. doi: 10.1007/s10555-017-9682-0 Widgerow AD King K Tocco-Tussardi I Banyard DA Chiang R Awad A . The burn wound exudate - an under-utilized resource. Burns (2015) 41:11–7. doi: 10.1016/j.burns.2014.06.002 Turner MD Nedjai B Hurst T Pennington DJ . Cytokines and chemokines: At the crossroads of cell signalling and inflammatory disease. Biochim Biophys Acta - Mol Cell Res (2014) 1843:2563–82. doi: 10.1016/j.bbamcr.2014.05.014 Dinarello CA Novick D Kim S Kaplanski G . Interleukin-18 and IL-18 binding protein. Front Immunol (2013) 4:289. doi: 10.3389/fimmu.2013.00289 Foessl I Haudum CW Vidakovic I Prassl R Franz J Mautner SI . MiRNAs as regulators of the early local response to burn injuries. Int J Mol Sci (2021) 22:1–18. doi: 10.3390/ijms22179209 Kotzbeck P Hofmann E Nischwitz SP Kamolz L-P . Differentiating local and systemic inflammatory responses to burn injuries. Burns (2019). doi: 10.1016/j.burns.2018.11.006 Laggner M Lingitz MT Copic D Direder M Klas K Bormann D . Severity of thermal burn injury is associated with systemic neutrophil activation. Sci Rep (2022) 12:1–10. doi: 10.1038/s41598-022-05768-w Jabeen S Clough ECS Thomlinson AM Chadwick SL Ferguson MWJ . Shah m. partial thickness wound: Does mechanism of injury influence healing? Burns (2019) 45:531–42. doi: 10.1016/j.burns.2018.08.010 Mortaz E Alipoor SD Adcock IM Mumby S Koenderman L . Update on neutrophil function in severe inflammation. Front Immunol (2018) 9:2171. doi: 10.3389/fimmu.2018.02171 Manley HR Keightley MC Lieschke GJ . The neutrophil nucleus: An important influence on neutrophil migration and function. Front Immunol (2018) 9:2867. doi: 10.3389/fimmu.2018.02867 Leliefeld PHC Wessels CM Leenen LPH Koenderman L Pillay J . The role of neutrophils in immune dysfunction during severe inflammation. Crit Care (2016) 20:1–9. doi: 10.1186/s13054-016-1250-4 O’Hare FM Watson W O’Neill A Grant T Onwuneme C Donoghue V . Neutrophil and monocyte toll-like receptor 4, CD11b and reactive oxygen intermediates, and neuroimaging outcomes in preterm infants. Pediatr Res (2015) 78:82–90. doi: 10.1038/pr.2015.66 Schmidt T Zündorf J Grüger T Brandenburg K Reiners A-L Zinserling J . CD66b overexpression and homotypic aggregation of human peripheral blood neutrophils after activation by a gram-positive stimulus. J Leukoc Biol (2012) 91:791–802. doi: 10.1189/jlb.0911483 Aratani Y . Myeloperoxidase: Its role for host defense, inflammation, and neutrophil function. Arch Biochem Biophys (2018) 640:47–52. doi: 10.1016/j.abb.2018.01.004 Wang J . Neutrophils in tissue injury and repair. Cell Tissue Res (2018) 371:531–9. doi: 10.1007/s00441-017-2785-7 Italiani P Boraschi D . From monocytes to M1/M2 macrophages: Phenotypical vs. functional differentiation. Front Immunol (2014) 5:514. doi: 10.3389/fimmu.2014.00514 Olingy CE San Emeterio CL Ogle ME Krieger JR Bruce AC Pfau DD . Non-classical monocytes are biased progenitors of wound healing macrophages during soft tissue injury. Sci Rep (2017) 7:1–16. doi: 10.1038/s41598-017-00477-1 Krzyszczyk P Schloss R Palmer A Berthiaume F . The role of macrophages in acute and chronic wound healing and interventions to promote pro-wound healing phenotypes. Front Physiol (2018) 9:419. doi: 10.3389/fphys.2018.00419 Chávez-Galán L Olleros ML Vesin D Garcia I . Much more than M1 and M2 macrophages, there are also CD169+ and TCR+ macrophages. Front Immunol (2015) 6:263. doi: 10.3389/fimmu.2015.00263 Jeschke MG van Baar ME Choudhry MA Chung KK Gibran NS Logsetty S . Burn injury. Nat Rev Dis Prim (2020) 6:1–25. doi: 10.1038/s41572-020-0145-5 Dhaiban S Al-Ani M Elemam NM Maghazachi AA . Targeting chemokines and chemokine receptors in multiple sclerosis and experimental autoimmune encephalomyelitis. J Inflammation Res (2020) 13:619–33. doi: 10.2147/JIR.S270872 Peters JH Tjabringa GS Fasse E de Oliveira VL Schalkwijk J Koenen HJPM . Co-Culture of healthy human keratinocytes and T-cells promotes keratinocyte chemokine production and RORγt-positive IL-17 producing T-cell populations. J Dermatol Sci (2013) 69:44–53. doi: 10.1016/j.jdermsci.2012.10.004 Landén NX Li D Ståhle M . Transition from inflammation to proliferation: a critical step during wound healing. Cell Mol Life Sci (2016) 73:3861–85. doi: 10.1007/s00018-016-2268-0 Matar DY Ng B Darwish O Wu M Orgill DP Panayi AC . Skin inflammation with a focus on wound healing. Adv Wound Care (2022) 12:1–61. doi: 10.1089/wound.2021.0126 Jeschke MG Gauglitz GG Kulp GA Finnerty CC Williams FN Kraft R . Long-term persistance of the pathophysiologic response to severe burn injury. PLoS One (2011) 6:e21245. doi: 10.1371/journal.pone.0021245 Browning JA Cindass R . Burn debridement, grafting, and reconstruction (2021). Available at: http://www.ncbi.nlm.nih.gov/pubmed/31869181. Barayan D Abdullahi A Vinaik R Knuth CM Auger C Jeschke MG . Interleukin-6 blockade, a potential adjunct therapy for post-burn hypermetabolism. FASEB J Off Publ Fed Am Soc Exp Biol (2021) 35:e21596. doi: 10.1096/fj.202100388R