Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 20 September 2022
Sec. Inflammation
This article is part of the Research Topic Understanding How Myeloid Cell Development and Function Meet Tissue Distinct Metabolic Requirements View all 8 articles

IFN-γ transforms the transcriptomic landscape and triggers myeloid cell hyperresponsiveness to cause lethal lung injury

Atul K. Verma&#x;Atul K. Verma1†Michael McKelvey&#x;Michael McKelvey2†Md Bashir UddinMd Bashir Uddin1Sunil PalaniSunil Palani1Meng NiuMeng Niu3Christopher BauerChristopher Bauer4Shengjun ShaoShengjun Shao1Keer Sun,*Keer Sun1,4*
  • 1Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, TX, United States
  • 2Department of Experimental Pathology, University of Texas Medical Branch, Galveston, TX, United States
  • 3Department of Genetics, Cell Biology & Anatomy, University of Nebraska Medical Center, Omaha, NE, United States
  • 4Department of Pathology and Microbiology, University of Nebraska Medical Center, Omaha, NE, United States

Acute Respiratory Distress Syndrome (ARDS) is an inflammatory disease that is associated with high mortality but no specific treatment. Our understanding of initial events that trigger ARDS pathogenesis is limited. We have developed a mouse model of inflammatory lung injury by influenza and methicillin-resistant Staphylococcus aureus (MRSA) coinfection plus daily antibiotic therapy. Using this pneumonic ARDS model, here we show that IFN-γ receptor signaling drives inflammatory cytokine storm and lung tissue damage. By single-cell RNA sequencing (scRNA-seq) analysis, we demonstrate that IFN-γ signaling induces a transcriptional shift in airway immune cells, particularly by upregulating macrophage and monocyte expression of genes associated with inflammatory diseases. Further evidence from conditional knockout mouse models reveals that IFN-γ receptor signaling in myeloid cells, particularly CD11c+ mononuclear phagocytes, directly promotes TNF-α hyperproduction and inflammatory lung damage. Collectively, the findings from this study, ranging from cell-intrinsic gene expression to overall disease outcome, demonstrate that influenza-induced IFN-γ triggers myeloid cell hyperresponsiveness to MRSA, thereby leading to excessive inflammatory response and lethal lung damage during coinfection.

Introduction

Acute respiratory distress syndrome (ARDS) is characterized by an overwhelming inflammatory response. The limited treatment option, with no effective therapy to prevent lung injury, has led to a mortality rate of approximately 40% in afflicted (1, 2). Hypercytokinemia, also known as “cytokine storm”, plays a central role in ARDS pathogenesis. As such, inflammatory cytokines TNF-α and IL-6 are significant predictors of severity and death of diseases including COVID-19 (3, 4). However, our understanding of initial events that trigger hypercytokinemia and ARDS is very limited.

Viral and bacterial pneumonia that activates excessive and dysregulated lung inflammation is the most frequent cause of acute lung injury (ALI) and ARDS (5, 6). Respiratory illness caused by influenza A virus (IAV) alone is generally self-limited. However, IAV infection predisposes hosts to secondary bacterial pneumonia, which is a major cause of hospitalization and death in influenza pandemics and epidemics (712). Staphylococcus aureus is one such opportunistic pathogen often implicated in post-influenza bacterial pneumonia. The increasing prevalence of methicillin-resistant S. aureus (MRSA) in communities and healthcare settings further complicates the problem. Particularly, IAV/MRSA coinfection can quickly develop into severe, necrotizing pneumonia with over 50% mortality, despite available flu vaccines, anti-viral drugs, and antibiotics (9, 12, 13).

Through post-influenza MRSA infection plus daily antibiotic treatment, we have developed a mouse model of ARDS that mimics the poor outcome of patients with influenza-complicated bacterial pneumonia. We show in this model that antibiotic therapy restores bacterial control but is not sufficient to rescue animals from inflammatory lung damage (14). Using this coinfection-induced ARDS model, we have shown that phagocyte oxidative burst promotes necrotic death of inflammatory cells, thereby resulting in lethal damage to surrounding tissue (15). Furthermore, we have demonstrated that the overproduction of proinflammatory cytokine TNF-α, in response to influenza-induced IFN-γ, drives lethal lung inflammation (16).

In the present study, we have performed single-cell RNA sequencing (scRNA-seq) to profile how IFN-γ signaling instructs lung inflammation. We demonstrate that IFN-γ induces a transcriptional shift in airway immune cells during IAV/MRSA coinfection. While IAV-induced Ifng expression is limited to the T cell cluster, myeloid cells are the main producers of Tnfa and Il1b in response to MRSA super-infection. These findings suggest an IFN-γ-mediated interplay between T cells and inflammatory cells in perpetuating immune dysfunction and inflammatory damage. Using a series of IFN-γ receptor 1 (Ifngr1) conditional knockout mice, we have further investigated the cellular pathways underlying IFN-γ-driven hypercytokinemia and their contributions to ARDS pathogenesis. We found that the proinflammatory cytokine response in the lung is primarily driven by IFN-γ signaling in myeloid cells, with the greatest contribution from mononuclear phagocytes. As such, mice with Ifngr1 deletion in CD11c+ myeloid cells exhibited significantly attenuated lung inflammation and mortality after IAV/MRSA coinfection and gentamicin treatment.

Materials and methods

Murine model of viral and bacterial infection

Specific pathogen-free, C57BL/6 WT, Ifng-/-, Ifngr1-/-, Ifngr1fl/fl, Lyz2Cre, Cd11cCre, Cx3cr1Cre, Mrp8Cre and mTmG (17) mice were purchased from the Jackson Laboratory (Bar Harbor, ME) and bred at the University of Nebraska Medical Center (UNMC) and University of Texas Medical Branch (UTMB) following Animal Care and Use Committee (IACUC) guidelines. All animal experiments were approved by UNMC and UTMB, and all experiments were carried out in accordance with UNMC and UTMB Assurance of Compliance with PHS Policy on Humane Care and Use of Laboratory Animals, which is on file with the Office of Protection from Research Risks, NIH.

Lyz2Cre, Cd11cCre, Cx3cr1Cre, and Mrp8Cre reporter mice were generated by crossing mTmG mice with corresponding Cre+ mice. Lyz2Cre_Ifngr1fl/fl, Cd11cCre_Ifngr1fl/fl, Cx3cr1Cre_Ifngr1fl/fl, and Mrp8Cre_Ifngr1fl/fl mice were generated by crossing Ifngr1fl/fl mice with corresponding Cre+ mice.

Viral challenge was performed with a sublethal (0.2LD50) dose of PR8 administered i.n. to anesthetized, sex and age-matched adult mice in 50 µl of sterile PBS. Viral titers were determined by plaque assays on Madin-Darby canine kidney cell monolayers. Bacterial super-challenge was performed seven days later. To induce bacterial pneumonia, anesthetized mice were inoculated i.n. 4~6x108 CFU/mouse of ATCC MRSA strain BAA-1695 in 50 μl of PBS. Bacterial burdens were measured by plating serial 10-fold dilutions of each sample onto blood agar plates.

Bronchoalveolar lavage cell analysis

BAL fluid (BALF) samples were collected by making an incision in the trachea and lavaging the lung twice with 0.8 ml PBS, pH 7.4. Total leukocyte counts were determined using a hemocytometer.

For flow cytometric analysis, BALF cells were incubated with 2.4G2 mAb and stained with allophycocyanin-conjugated anti-CD11c (Biolegend), BUV395-conjugated or BV510-conjugated anti-CD11b (BD Biosciences), FITC-conjugated or PE-Cy7-conjugated anti-Ly6G (clone 1A8, Biolegend), PerCP-Cy5.5-conjugated (eBiosciences) or PE-conjugated anti-Ly6C (Biolegend), and BV421-conjugated or PE-conjugated anti-Siglec-F (Biolegend) mAbs. The stained cells were analyzed on a BD LSRII-green using BD FACSDiva and FlowJo software analysis.

Determination of cytokine production by ELISA

BALF were harvested and assayed for IFN-γ, TNF-α, IL-1β, and IL-6 by ELISA using commercially available kits from BD Biosciences and R&D Systems (Minneapolis, MN).

Single cell RNA sequencing

scRNA-seq was performed with the assistance of the UNMC Genomics Core. The core operates a 10x Genomics Chromium instrument for isolation of single cells and generation of Illumina compatible sequencing libraries (10x Genomics, Pleasanton, CA). The instrument captures individual cells using a gel bead in emulsion (GEM) technology. The droplets contain the reagents for reverse transcription and all transcripts within each cell are uniquely barcoded (all transcripts within a given cell contain the same barcode for cell identification) and Illumina compatible cDNA libraries created are for each cell within the given sample. Prior to beginning the capture, BALF cell suspensions were evaluated by light microscopy with respect to debris and viability. The single‐cell suspensions were diluted to the working concentration required to 8000 target cells captured, with Ifngr1-/- or Ifng-/- cells and corresponding WT controls each divided in half for library preparation. The pooled cells were captured, lysed, RNA reverse transcribed and barcoded using Chromium Single Cell 3’ Reagent Kits v3 reagents. Libraries are quantified by qPCR using the KAPA Library Quant Kit (Illumina) from KAPA Biosystems (Roche, Pleasonton, CA) and were loaded at a concentration of 1.3 pM on Illumina sequencing flowcells. The samples were sequenced on an Illumina NextSeq500 instrument (as recommended for 3’ end sequencing of transcripts from single cells).

ScRNA-seq analyses

The 10x Genomics Cell ranger analysis pipeline (Version 3.0.0) was used for demultiplexing and generating the Feature-Barcode Matrices (18). Briefly, cellranger mkfastq was used to convert the raw base call (BCL) files of scRNA-seq into FASTQ files. Then cellranger count was used for alignment with reference genome mm10, filtering, and processing barcode data to create the feature-barcode matrices for each sample; cellranger aggr was used to normalize the samples to the same sequencing depth and then recomputing the feature-barcode matrices and perform analysis on the combined data including clustering and gene express analysis. We used 10x Genomics Loupe cell browser for differential expression analysis and creating t-SNE plots on which we mapped the marker gene expressions and the clustering results.

Additional bioinformatic analysis was performed using IPA (QIAGEN Inc.). We provided the DEG lists (fold change >4) to IPA for its core analysis and then overlaid with the global molecular network in the Ingenuity pathway knowledge base. As a result, the biological functions, disease annotations, and canonical pathways that were enriched in the DEG lists were identified and analyzed.

Lung histology analysis

Mice were euthanized seven days after MRSA super-infection and the lungs were removed for histological analyses. Paraffin-embedded tissues were sectioned to a thickness of 5 µm and stained with hematoxylin and eosin (H&E) using standard methods. Whole mount H&E-stained lung tissues were scanned using Leica Aperio LV1 scanner and software, and semi quantitatively assessed at low-power (40×) for the proportion of parenchyma with alveoli containing intraluminal material (proteinaceous exudate and/or fibrin) in the background of interstitial expansion and inflammation. Each lung was scored by the relative amount of abnormal tissue as follows: normal-0, 1-25%-1, 26-50%-2, 51-75%-3, >76%-4 (15). Digital images were generated using Leica Biosystems Aperio ImageScope 12.

Evaluation of airway damage

BALF were harvested and assayed for albumin by ELISA using a commercially available kit from Bethyl Laboratories (Montgomery, TX). Total protein levels and lactic acid dehydrogenase (LDH) activities in BALF were analyzed by a Micro bicinchoninic acid assay (BCA) protein assay kit (Thermo Scientific) and a LDH cytotoxicity assay kit (Thermo Scientific), respectively.

Treatment with antibiotics

Mice were intraperitoneally (i.p.) injected with a therapeutic dose of vancomycin (300 mg/kg) or gentamicin (100 mg/kg) beginning 2 or 4 h after MRSA infection and then followed by 150 mg/kg/day for vancomycin or 50 mg/kg/day for gentamicin, respectively. All antibiotic treatment injections continued through day seven after MRSA infection. Minimum inhibitory concentration (MIC) values of antibiotics for MRSA BAA-1695 were ≤ 12.5 µg/ml for gentamicin and ≤ 1.17 µg/ml for vancomycin in vitro (14).

Statistics

Significant differences between experimental groups were determined using a two-tailed Student t-test (to compare two samples), or an ANOVA analysis followed by Tukey’s multiple comparisons test (to compare multiple samples) in GraphPad Prism 9 (La Jolla, CA). Survival analyses were performed using the log-rank test. For all analyses, a P value <0.05 was considered to be significant.

Results

IFN-γ receptor signaling promotes hypercytokinemia during post-influenza MRSA pneumonia

We have developed a mouse model of ALI induced by post-influenza MRSA infection (15). This secondary bacterial pneumonia model replicates the poor outcome of patients with pneumonic ARDS, despite appropriate antibiotic treatment (8, 13, 1921). Specifically, C57BL/6 mice were challenged with IAV strain A/PR8/1934 followed by intranasal (i.n.) super-infection of MRSA seven days later (15). IAV/MRSA coinfected mice were treated daily with antibiotic vancomycin or gentamicin for effective bacterial control. However, antibiotic therapy alone is not sufficient to improve animal survival (14, 16).

We have recently shown in this model that combining antibiotic treatment with IFN-γ neutralizing antibodies prevents hypercytokinemia and thereby improves protection against ALI and animal mortality (16). This IFN-γ-driven excessive lung inflammation was further verified in the global Ifngr1 gene-deficient (Ifngr1−/−) mouse model (Figure 1A). Vancomycin-treated WT and Ifngr1-/- mice exhibited similar bacterial burdens after MRSA infection alone and IAV/MRSA coinfection (Figure 1B). Interestingly, compared with respective WT controls, Ifngr1-/- mice exhibited increased airway neutrophil accumulation after coinfection (Figure 1C). In contrast, acute proinflammatory cytokines TNF-α, IL-1β and IL-6 were diminished in coinfected Ifngr1-/- mice (Figure 1D), while their IFN-γ production was not significantly different from coinfected WT controls (Figure 1E). Importantly, acute lung damage was attenuated in Ifngr1-/- mice, as indicated by significantly reduced BALF protein and LDH levels 7 days post-MRSA super-infection (dps) (Figures 1F-G). Together, these results establish that IFN-γ receptor signaling promotes hypercytokinemia and lung damage during secondary MRSA pneumonia.

FIGURE 1
www.frontiersin.org

Figure 1 IFN-γ signaling induces hypercytokinemia and acute lung damage. (A) Experimental scheme. Mice were infected with PR8 virus and then challenged with MRSA on day seven post-influenza infection. All PR8-infected (Flu), MRSA-infected (MRSA), or coinfected mice (Flu+MRSA) were treated with vancomycin 2 h after MRSA or PBS inoculation. (B) Lung bacterial burdens, (C) BALF levels of (mean±SD, n≥4) AM, monocytes (Mo) and neutrophils (PMN), (D) TNF-α, IL-1β and IL-6, and (E) IFN-γ at 24 h after PR8 and/or MRSA infection. (F) BALF levels of protein and (G) lactate dehydrogenase (LDH) at 7 dps in WT and Ifngr1-/- mice (mean±SEM, n≥5) after PR8/MRSA coinfection. Data in (F) represented protein levels relative to naïve mice. AU, arbitrary unit. *P<0.05, **P<0.01, ***P<0.001, t-test. Data shown are representative of at least two independent experiments.

IFN-γ signaling transforms the transcriptomic landscape of airway immune cells

To determine the impact of IFN-γ signaling on the lung inflammatory profile, we performed scRNA-seq of bronchoalveolar lavage fluid (BALF) cells 24 h after IAV/MRSA coinfection of Ifngr1-/- and WT mice side-by-side (Figure 2A). The clustering analysis revealed a clear segregation of Ifngr1-/- cell populations from WT controls (Figure 2B). We next compared gene expressions between WT and Ifngr1-/- cells. A total of 4676 differentially expressed genes (DEGs, P<0.05), including 1954 upregulated and 2722 downregulated, were identified in WT cells (Figure 2C).

FIGURE 2
www.frontiersin.org

Figure 2 IFN-γ induces fundamental changes in the airway transcriptomic landscape. (A) Experimental scheme for scRNA-seq analysis. (B) t-SNE projections distinguishing Ifngr1-/- cells (red) from WT controls (brown) in the same experiment. Each dot represents a cell. (C) Volcano plot for DEGs in WT cells as compared to Ifngr1-/- datasets. (D) t-SNE projections and (E) volcano plot for DEGs in WT cells (blue) as compared to Ifng-/- datasets (green). (F) t-SNE projections showing overlap among Ifngr1-/-, Ifng-/- and WT cells. Four datasets were combined from two independent experiments. (G) Diagram showing top DEG (fold change >4) associated with inflammatory disease conditions, overlapping with genes related to 1respiratory disease, 2cell death and survival, and 3free radical scavenging in WT mice as compared to corresponding Ifngr1-/- animals. Red: upregulated; blue: downregulated in WT cells. (H) Representative classification of diseases and functions mediated by IFN-γ receptor signaling. Categories are shown in terms of the –log10 (P-value), as represented by the left y-axis, and the number of DEGs counted, represented by the right y-axis.

To determine whether IFN-γ signaling is solely responsible for the fundamental changes in the transcriptomic landscape, we further performed scRNA-seq of BALF cells from Ifng-/- mice as compared with corresponding WT controls (Figures 2D-E). Like Ifngr1-/- mice, there was very limited overlap between Ifng-/- and corresponding WT landscapes on the t-SNE plot. Importantly, the combined analysis of four datasets revealed that Ifngr1-/- and Ifng-/- cell clusters from two independent experiments were more intermixed than with their respective WT controls (Figure 2F). These results establish that IFN-γ signaling induces fundamental changes in the airway transcriptome after IAV/MRSA coinfection.

Using the ingenuity pathways analysis (IPA) system with a threshold of −log (P−value) > 4 and DEG counts ≥ 10, we found that IFN-γ-activated molecular and cellular functions were associated with cell death and survival, and inflammatory diseases (Figure 2G). Further disease and function analysis indicates that apart from inflammatory response, IFN-γ receptor signaling is also involved in multiple organ/tissue disorders, organismal injury and abnormalities, as well as organismal survival (Figure 2H). These fundamental changes in the transcriptomic profile suggest that IFN-γ plays a central role in causing immune cell dysfunction and inflammatory tissue damage during post-influenza MRSA pneumonia.

IFN-γ upregulates monocyte and macrophage expression of genes associated with inflammatory diseases

To further differentiate the effect of IFN-γ signaling on airway immune cell subsets, we performed pairwise comparison of Ifngr1-/- and Ifng-/- transcription profiles with corresponding WT controls within three major cell clusters, representing neutrophils (Lyz2+Cxcr2+), T cells (Cd3+), and monocytes/macrophages (Mo/Mψ, Lyz2+Ccr2+) (Figures 3A-B). Of note, Cd11c+Marco+ alveolar macrophages (AM) contributed to the Mo/Mψ cluster but only a minor proportion (Figure 3A).

FIGURE 3
www.frontiersin.org

Figure 3 IFN-γ induces robust upregulation of inflammatory disease genes in macrophages and monocytes. (A) Visualization of expression of cell type signature genes and (B) t-SNE plot and the percentages of three major clusters in scRNA-seq datasets combined of WT, Ifngr1-/- and Ifng-/- BALF cells. (C) Single-cell expression of six selected inflammatory genes. (D) Diagram showing top DEGs (fold change > 4) associated with inflammatory disease conditions in WT cell clusters as compared with corresponding Ifngr1-/- cells, overlapping with genes related to 1respiratory disease, 2cell death and survival, and 3free radical scavenging. Red: upregulated; blue: downregulated in WT cells.

In line with IFN-γ-induced fundamental changes in the airway transcriptome, we detected broad Ifngr1 expression in these immune cells (Figure 3C). As expected, Ifng expression was limited to the T cell cluster. On the other hand, inflammatory cells, as indicated by their Lyz2 and Nox2 expression, appeared to be the predominant producers of Tnfa and Il1b. Surprisingly, Il6 mRNA was detected only in a small number of cells in the Mo/Mψ cluster. Thus, it is possible that IL-6 is primarily produced by lung stromal cells rather than immune cells. Together, these findings imply an IFN-γ-initiated interplay of T cells, inflammatory cells, and presumably, lung tissue cells in driving hypercytokinemia after IAV/MRSA coinfection.

A further IPA analysis of these IFN-γ-induced DEGs revealed an enrichment of top genes, i.e., Log (P−value) > 4 and unique molecular identifier (UMI) > 0.1, associated with inflammatory response and organismal injury, overlapping with genes related to respiratory disease, cell death and survival, and free radical scavenging (Figure 3D). It is noteworthy that these top DEGs associated with respiratory inflammation, e.g. > 50-fold increased Cxcl9 and Il12b and >15-fold decreased Arg1 in WT cells, appeared to be most active within myeloid cells, particularly the Mo/Mψ cluster (Figure 4). This suggests that in response to IAV/MRSA coinfection, T cell-derived IFN-γ may directly activate myeloid cells to enhance their proinflammatory properties.

FIGURE 4
www.frontiersin.org

Figure 4 IFN-γ primarily induces transcriptomic changes in myeloid cells. ScRNA-seq analysis of DEGs (P<0.05) in WT (A) neutrophil, (B) Mo/Mψ, and (C) T cell clusters as compared with corresponding Ifngr1-/- (left) and Ifng-/- (right) subset in two independent experiments. Volcano plots show the logarithm of fold change (log2) on x-axis and UMI of each transcript on y-axis.

IFN-γ signaling in myeloid cells drives hypercytokinemia and lethal IAV/MRSA coinfection

To determine the contribution of different myeloid cell subsets to coinfection pathogenesis, we developed cre-lox mouse models with Ifngr1 deletion (22) in myeloid cells (Lyz2Cre) (23), Mo/Mψ (Cx3cr1Cre) (24), and neutrophils (Mrp8Cre) (25) (Figure 5A). These conditional Ifngr1 knockout mice displayed similar lung bacterial burdens as WT controls after IAV/MRSA coinfection and vancomycin treatment (Figure 5B). Nonetheless, deletion of Ifngr1 in myeloid cells (Lyz2Cre_Ifngr1fl/fl) led to nearly 6- and 10-fold decreased TNF-α and IL-6 production, respectively (Figures 5C-D). These results indicate that IFN-γ-activated myeloid cells are responsible for driving hypercytokinemia after coinfection. Mice carrying Ifngr1 deficiency in Mo/Mψ (Cx3cr1Cre_Ifngr1fl/fl) also exhibited reduced TNF-α production; their IL-6 production, however, was comparable to Ifngr1fl/fl controls. In contrast, deletion of Ifngr1 primarily in neutrophils (Mrp8Cre_Ifngr1fl/fl) had no significant effect on TNF-α or IL-6 response. These results suggest that IFN-γ-activated mononuclear phagocytes play a more key role in TNF-α hyperproduction.

FIGURE 5
www.frontiersin.org

Figure 5 IFN-γ signaling in myeloid cells drives TNF-α and IL-6 hyperproduction and animal mortality. (A) Flow cytometry analysis of percentages (mean of ‗3 mice/group) of EGFP-expressing AM, monocytes (Mo), neutrophils (PMN) in Lyz2Cre_mTmG, Cx3cr1Cre_mTmG, and Mrp8Cre_mTmG mice (black line), and mTmG controls (solid grey) 24 h after PR8/MRSA coinfection. (B) Lung bacterial burdens, (C) BALF TNF-α and (D) IL-6 levels at 24 h, and (E) animal survival of WT, Lyz2Cre, Ifngr1fl/fl, Lyz2Cre_Ifngr1fl/fl, Cx3cr1Cre_Ifngr1fl/fl, and Mrp8Cre_Ifngr1fl/fl mice after PR8/MRSA coinfection. All mice were treated with vancomycin starting 2 h after MRSA infection. *P<0.05, **P<0.01, ***P<0.001, Tukey’s multiple comparisons test (B-D) or log-rank test (E). Data shown were combined from more than two independent experiments.

In agreement with their diminished proinflammatory cytokine response, coinfected Lyz2Cre_Ifngr1fl/fl mice exhibited significantly improved survival with vancomycin treatment, as compared with mice carrying Ifngr1 deletion limited to either Mo/Mψ (Cx3cr1Cre_Ifngr1fl/fl) or neutrophils (Mrp8Cre_Ifngr1fl/fl) (Figure 5E). These results suggest that IFN-γ-activated myeloid cells are responsible for driving hypercytokinemia and lethal IAV/MRSA coinfection.

IFN-γ directly activates mononuclear phagocytes to promote TNF-α hyperproduction and inflammatory lung damage

Our results indicate that IFN-γ-activated mononuclear phagocytes play a dominant role in driving TNF-α hyperproduction (Figure 5C). Considering that CX3CR1 is also expressed in lymphocytes such as NK and T cells (24), we further developed Cd11cCre_Ifngr1fl/fl mouse model (26) to induce Ifngr1 deletion specifically in CD11c+ mononuclear phagocytes (Figure 6). We have previously shown that compared with bacteriolytic vancomycin, treatments with protein synthesis inhibitor antibiotics linezolid and gentamicin provide better protection against S. aureus α-toxin-induced direct lung injury during IAV/MRSA coinfection (14, 16). Thus, we utilized this more sensitive lung inflammatory damage model, i.e., IAV/MRSA coinfection plus gentamicin treatment, to further differentiate the role of myeloid cell subsets in lung injury (Figure 6A).

FIGURE 6
www.frontiersin.org

Figure 6 IFN-γ directly activates mononuclear phagocytes to promote TNF-α production. (A) Experimental scheme. All mice were treated with gentamicin starting 4 h after MRSA super-infection. (B) Flow cytometry analysis of percentages (mean of ‗3 mice/group) of EGFP-expressing AMs, monocytes (Mo), and neutrophils (PMN) in Cd11cCre_mTmG reporter mice (black line), and mTmG controls (solid grey) 24 h after PR8/MRSA coinfection. (C-D) Lung viral and bacterial burdens and (E) TNF-α and IL-6 levels 24 h after PR8/MRSA coinfection of Ifngr1fl/fl, Lyz2Cre_Ifngr1fl/fl, Cx3cr1Cre_Ifngr1fl/fl, and Cd11cCre_Ifngr1fl/fl mice. *P<0.05, **P<0.01, Tukey’s multiple comparisons test. Data shown are representative of two independent experiments.

Deletion of Ifngr1 in CD11c+ myeloid cells had no significant impact on lung viral and bacterial control (Figures 6C-D). Conversely, BALF level of TNF-α was significantly reduced in Cd11cCre_Ifngr1fl/fl animals, in a pattern similar to Lyz2Cre_Ifngr1fl/fl and Cx3cr1Cre_Ifngr1fl/fl mice (Figure 6E). This finding suggests that IFN-γ signaling in monocytes and macrophages drives TNF-α production. In line with that, Cx3cr1Cre_Ifngr1fl/fl mice exhibited reduced histopathological scores at 7 dps, suggesting that IFN-γ signaling in mononuclear phagocytes are responsible for triggering severe lung damage in Ifngr1fl/fl controls (Figure 7A). Importantly, gentamicin treatment significantly improved the survival of Cd11cCre_Ifngr1fl/fl and Cx3cr1Cre_Ifngr1fl/fl mice after IAV/MRSA coinfection (Figure 7B). These findings indicate that IFN-γ directly activates monocytes and macrophages to elicit TNF-α production and lethal lung inflammation during post-influenza MRSA pneumonia.

FIGURE 7
www.frontiersin.org

Figure 7 IFN-γ directly activates mononuclear phagocytes to promote acute lung damage. (A) Experimental scheme, and lung histopathology (H&E) and histopathologic scores (each symbol represents one mouse) at 7 dps. (B) Survival of Ifngr1fl/fl, Cx3cr1Cre_Ifngr1fl/fl, Cd11cCre_Ifngr1fl/fl, Lyz2Cre_Ifngr1fl/fl and Mrp8Cre_Ifngr1fl/fl mice after PR8/MRSA coinfection. *P<0.05, **P<0.01, t-test (A), or log-rank test (B). Data shown are representative of at least two independent experiments.

Discussion

In the present study, we have investigated the signaling pathway underlying IFN-γ-driven lethal lung inflammation during IAV/MRSA coinfection. IAV infection is known to induce prominent IFN-γ production by T cells. Here we show that IFN-γ signaling induces fundamental changes in the airway transcriptome in response to MRSA infection, characterized by upregulations of genes associated with immune cell dysfunction and inflammatory tissue damage. Using a series of Ifngr1 conditional knockout mouse models, we further demonstrate that IFN-γ signaling in myeloid cells, particularly mononuclear phagocytes, directly activates TNF-α production and lethal lung damage. Collectively, these findings, ranging from cell-intrinsic gene regulation to overall impact on disease outcome, suggest that IFN-γ induces hyperresponsiveness of monocytes and macrophages to trigger ARDS pathogenesis during influenza-complicated MRSA infection.

Secondary bacterial infection following influenza is known to cause severe pneumonia. However, the reported studies generally focus on how IAV infection impairs antibacterial immunity and thereby increases host susceptibility to secondary bacterial infection (2734). In the current model, IAV/MRSA-coinfected mice were treated with adequate antibiotics to mediate bacterial killing. Thus, our model mimics the clinical scenario where critically ill patients develop ARDS and even death, despite appropriate antimicrobial therapies. This unique model allows us to investigate how dysregulated inflammatory response contributes to ARDS pathogenesis, independent of defective bacterial control.

Using Ifngr1 conditional knockout mouse models, we show that IFN-γ signaling in myeloid cells is sufficient to trigger hypercytokinemia and ARDS pathogenesis. In agreement, our scRNA-seq data demonstrate that IFN-γ-induced gene upregulation in WT myeloid cells presents an inflammatory signature, e.g., Cxcl9 in monocytes and Cxcl10 in neutrophils. Of particular interest, CXCL10-CXCR3 signaling has been shown to promote oxidative burst in neutrophils, thereby leading to fulminant lung inflammation (3537). In contrast, genes upregulated in Ifngr1-/- and Ifng-/- cells, such as Arg1 and Cd163 genes, are associated with immunoregulatory function and tissue repair (38).

The hallmark of ARDS is inflammatory cell accumulation in the lung. Influenza triggers intensive recruitment of inflammatory monocytes (39), while MRSA infection induces heightened neutrophil infiltration. Interestingly, IAV/MRSA coinfection-induced neutrophilic infiltration increases in the global Ifngr1 knockout mice (Figure 1B). Nonetheless, deletion of Ifngr1 only in myeloid cells did not appear to have a significant impact on the airway neutrophil accumulation (data not shown). These observations suggest that rather than promoting inflammatory infiltration, dysregulation of inflammatory cell function is the key mechanism in IFN-γ-driven alveolar inflammation and animal mortality.

Mononuclear phagocytes, particularly macrophages and monocytes, are known for their functional plasticity in response to cytokine milieu (40). In agreement, here we show that mononuclear phagocytes play a dominant role in driving TNF-α production and ARDS pathogenesis. Nonetheless, future studies are necessary to establish whether IFN-γ signaling promotes TNF-α production directly by this inflammatory cell subset or indirectly by epithelial cells. In contrast, deletion of Ifngr1 only in neutrophils (Mrp8Cre_Ifngr1fl/fl) has no significant effect on inflammatory cytokine response and animal mortality (Figure 6). These results suggest that the regulatory effect of IFN-γ on neutrophils is largely indirect, again highlighting the significance of IFN-γ-activated mononuclear phagocyte subset in ARDS pathogenesis.

Compared to the magnitude of inflammatory cell influx, necrotic cell death is more important in causing acute lung injury. Our previous study has demonstrated that oxidative stress derived from phagocyte NADPH activity promotes necrotic cell death (15). Inflammatory cytokines are responsible for activation of inflammatory cell recruitment and oxidative burst. Interestingly, MRSA infection alone induces intensive production of proinflammatory cytokines, independent of IFN-γ signaling (Figure 1D-E). In this case, however, antibiotic therapy is effective in improving animal survival (14). It has been reported that IFN-γ, in synergy with TNF-α, triggers inflammatory cell death and directly contributes to lethal tissue damage in SARS-CoV-2 infection and cytokine storm syndromes. Furthermore, blocking both IFN-γ and TNF-α with specific antibodies protected against inflammatory cell death and animal mortality from cytokine shock, sepsis, and hemophagocytic lymphohistiocytosis, in addition to SARS-CoV-2 infection (41). Thus, in addition to stimulating TNF-α production, IFN-γ is likely involved in other immunopathogenic aspects of IAV/MRSA coinfection, such as promoting inflammatory cell death.

In summary, our findings from the current study, ranging from cell-specific gene regulation to overall disease outcome, indicate that IFN-γ directly activates mononuclear phagocytes to trigger TNF-α hyperproduction and ARDS pathogenesis during IAV/MRSA coinfection. Considering the prominent IFN-γ-associated cytokine storm in various infectious and autoinflammatory conditions, strategies to target these early events represent potential therapeutic intervention for attenuating inflammatory lung damage in ARDS-prone diseases.

Data availability statement

The data presented in the study are deposited in the GEO repository, accession number GSE212402.

Ethics statement

The animal study was reviewed and approved by University of Nebraska Medical Center (UNMC) and University of Texas Medical Branch (UTMB) Animal Care and Use Committee (IACUC).

Author contributions

Conceived and designed the experiments: KS. Performed the experiments: AV, MM, MU, SP, CB and SS. Analyzed the data: KS and MN. Wrote the paper: KS and MM. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by NIH grant R01 HL118408 to KS. We thank Dr. James Eudy for advice on single-cell sequencing. The UNMC DNA Sequencing Core receives partial support from P20GM103427, P30GM110768, and P30CA036727. We also thank the Bioinformatics and Systems Biology Core at UNMC for providing “NAME” data analysis services, which receives support from Nebraska Research Initiative (NRI) and NIH (2P20GM103427 and 5P30CA036727). The authors also thank UTMB Anatomic Pathology.

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.

References

1. Luyt CE, Bouadma L, Morris AC, Dhanani JA, Kollef M, Lipman J, et al. Pulmonary infections complicating ARDS. Intensive Care Med (2020) 46:2168–83. doi: 10.1007/s00134-020-06292-z

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Menk M, Estenssoro E, Sahetya SK, Neto AS, Sinha P, Slutsky AS, et al. Current and evolving standards of care for patients with ARDS. Intensive Care Med (2020) 46:2157–67. doi: 10.1007/s00134-020-06299-6

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Del Valle DM, Kim-Schulze S, Huang HH, Beckmann ND, Nirenberg S, Wang B, et al. An inflammatory cytokine signature predicts COVID-19 severity and survival. Nat Med (2020) 26:1636–43. doi: 10.1038/s41591-020-1051-9

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Abou-Arab O, Bennis Y, Gauthier P, Boudot C, Bourdenet G, Gubler B, et al. Association between inflammation, angiopoietins, and disease severity in critically ill COVID-19 patients: a prospective study. Br J Anaesth (2021) 126:e127–30. doi: 10.1016/j.bja.2020.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Bauer TT, Ewig S, Rodloff AC, Muller EE. Acute respiratory distress syndrome and pneumonia: a comprehensive review of clinical data. Clin Infect Dis (2006) 43:748–56. doi: 10.1086/506430

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Matthay MA, Zemans RL, Zimmerman GA, Arabi YM, Beitler JR, Mercat A, et al. Acute respiratory distress syndrome. Nat Rev Dis Primers (2019) 5:18. doi: 10.1038/s41572-019-0069-0

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Louria DB, Blumenfeld HL, Ellis JT, Kilbourne ED, Rogers DE. Studies on influenza in the pandemic of 1957-1958. II. pulmonary complications of influenza. J Clin Invest (1959) 38:213–65. doi: 10.1172/JCI103791

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Finelli L, Fiore A, Dhara R, Brammer L, Shay DK, Kamimoto L, et al. Influenza-associated pediatric mortality in the united states: increase of staphylococcus aureus coinfection. Pediatrics (2008) 122:805–11. doi: 10.1542/peds.2008-1336

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Randolph AG, Vaughn F, Sullivan R, Rubinson L, Thompson BT, Yoon G, et al. Critically ill children during the 2009-2010 influenza pandemic in the united states. Pediatrics (2011) 128:e1450–8. doi: 10.1542/peds.2011-0774

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Shah NS, Greenberg JA, McNulty MC, Gregg KS, J.t. Riddell JE, Weber DM, et al. Bacterial and viral co-infections complicating severe influenza: Incidence and impact among 507 U.S. patients, 2013-14. J Clin Virol Off Publ Pan Am Soc Clin Virol (2016) 80:12–9. doi: 10.1016/j.jcv.2016.04.008

CrossRef Full Text | Google Scholar

11. Williams DJ, Hall M, Brogan TV, Farris RW, Myers AL, Newland JG, et al. Influenza coinfection and outcomes in children with complicated pneumonia. Arch Pediatr Adolesc Med (2011) 165:506–12. doi: 10.1001/archpediatrics.2010.295

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Rice TW, Rubinson L, Uyeki TM, Vaughn FL, John BB, Miller RR 3rd, et al. Critical illness from 2009 pandemic influenza a virus and bacterial coinfection in the united states. Crit Care Med (2012) 40:1487–98. doi: 10.1097/CCM.0b013e3182416f23

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Vardakas KZ, Matthaiou DK, Falagas ME. Incidence, characteristics and outcomes of patients with severe community acquired-MRSA pneumonia. Eur Respir J Off J Eur Soc Clin Respir Physiol (2009) 34:1148–58. doi: 10.1183/09031936.00041009

CrossRef Full Text | Google Scholar

14. Verma AK, Bauer C, Yajjala VK, Bansal S, Sun K. Linezolid attenuates lethal lung damage during post-influenza methicillin-resistant staphylococcus aureus pneumonia. Infection Immun (2019) 87(10):e00538–19. doi: 10.1128/IAI.00538-19

CrossRef Full Text | Google Scholar

15. Sun K, Yajjala VK, Bauer C, Talmon GA, Fischer KJ, Kielian T, et al. Nox2-derived oxidative stress results in inefficacy of antibiotics against post-influenza s. aureus pneumonia. J Exp Med (2016) 213:1851–64. doi: 10.1084/jem.20150514

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Verma AK, Bauer C, Palani S, Metzger DW, Sun K. IFN-gamma drives TNF-alpha hyperproduction and lethal lung inflammation during antibiotic treatment of postinfluenza staphylococcus aureus pneumonia. J Immunol (2021) 207:1371–6. doi: 10.4049/jimmunol.2100328

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Verma AK, Bansal S, Bauer C, Muralidharan A, Sun K. Influenza infection induces alveolar macrophage dysfunction and thereby enables noninvasive streptococcus pneumoniae to cause deadly pneumonia. J Immunol (2020) 205:1601–7. doi: 10.4049/jimmunol.2000094

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Niu M, Morsey B, Lamberty BG, Emanuel K, Yu F, Leon-Rivera R, et al. Methamphetamine increases the proportion of SIV-infected Microglia/Macrophages, alters metabolic pathways, and elevates cell death pathways: A single-cell analysis. Viruses 12 (2020) 12(11):1297. doi: 10.3390/v12111297

CrossRef Full Text | Google Scholar

19. McDanel JS, Perencevich EN, Storm J, Diekema DJ, Herwaldt L, Johnson JK, et al. Increased mortality rates associated with staphylococcus aureus and influenza Co-infection, Maryland and Iowa, USA(1). Emerg Infect Dis (2016) 22:1253–6. doi: 10.3201/eid2207.151319

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Metersky ML, Masterton RG, Lode H, File TM Jr., Babinchak T. Epidemiology, microbiology, and treatment considerations for bacterial pneumonia complicating influenza. Int J Infect Dis IJID Off Publ Int Soc Infect Dis (2012) 16:e321–31. doi: 10.1016/j.ijid.2012.01.003

CrossRef Full Text | Google Scholar

21. Centers for Disease, and Prevention. Severe coinfection with seasonal influenza a (H3N2) virus and staphylococcus aureus - Maryland, February-march 2012. MMWR. Morbidity mortality weekly Rep (2012) 61:289–91.

Google Scholar

22. Lee SH, Carrero JA, Uppaluri R, White JM, Archambault JM, Lai KS, et al. Identifying the initiating events of anti-listeria responses using mice with conditional loss of IFN-gamma receptor subunit 1 (IFNGR1). J Immunol (2013) 191:4223–34. doi: 10.4049/jimmunol.1300910

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Clausen BE, Burkhardt C, Reith W, Renkawitz R, Forster I. Conditional gene targeting in macrophages and granulocytes using LysMcre mice. Transgenic Res (1999) 8:265–77. doi: 10.1023/A:1008942828960

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Yona S, Kim KW, Wolf Y, Mildner A, Varol D, Breker M, et al. Fate mapping reveals origins and dynamics of monocytes and tissue macrophages under homeostasis. Immunity (2013) 38:79–91. doi: 10.1016/j.immuni.2012.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Passegue E, Wagner EF, Weissman IL. JunB deficiency leads to a myeloproliferative disorder arising from hematopoietic stem cells. Cell (2004) 119:431–43. doi: 10.1016/j.cell.2004.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Caton ML, Smith-Raska MR, Reizis B. Notch-RBP-J signaling controls the homeostasis of CD8- dendritic cells in the spleen. J Exp Med (2007) 204:1653–64. doi: 10.1084/jem.20062648

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lee MH, Arrecubieta C, Martin FJ, Prince A, Borczuk AC, Lowy FD. A postinfluenza model of staphylococcus aureus pneumonia. J Infect Dis (2010) 201:508–15. doi: 10.1086/650204

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Iverson AR, Boyd KL, McAuley JL, Plano LR, Hart ME, McCullers JA. Influenza virus primes mice for pneumonia from staphylococcus aureus. J Infect Dis (2011) 203:880–8. doi: 10.1093/infdis/jiq113

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Small CL, Shaler CR, McCormick S, Jeyanathan M, Damjanovic D, Brown EG, et al. Influenza infection leads to increased susceptibility to subsequent bacterial superinfection by impairing NK cell responses in the lung. J Immunol (2010) 184:2048–56. doi: 10.4049/jimmunol.0902772

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kudva A, Scheller EV, Robinson KM, Crowe CR, Choi SM, Slight SR, et al. Influenza a inhibits Th17-mediated host defense against bacterial pneumonia in mice. J Immunol (2011) 186:1666–74. doi: 10.4049/jimmunol.1002194

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Robinson KM, Lee B, Scheller EV, Mandalapu S, Enelow RI, Kolls JK, et al. The role of IL-27 in susceptibility to post-influenza staphylococcus aureus pneumonia. Respir Res (2015) 16:10. doi: 10.1186/s12931-015-0168-8

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Robinson KM, McHugh KJ, Mandalapu S, Clay ME, Lee B, Scheller EV, et al. Influenza a virus exacerbates staphylococcus aureus pneumonia in mice by attenuating antimicrobial peptide production. J Infect Dis (2014) 209:865–75. doi: 10.1093/infdis/jit527

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Robinson KM, Choi SM, McHugh KJ, Mandalapu S, Enelow RI, Kolls JK, et al. Influenza a exacerbates staphylococcus aureus pneumonia by attenuating IL-1beta production in mice. J Immunol (2013) 191:5153–9. doi: 10.4049/jimmunol.1301237

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Martinez-Colon GJ, Warheit-Niemi H, Gurczynski SJ, Taylor QM, Wilke CA, Podsiad AB, et al. Influenza-induced immune suppression to methicillin-resistant staphylococcus aureus is mediated by TLR9. PloS Pathog (2019) 15:e1007560. doi: 10.1371/journal.ppat.1007560

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Ichikawa A, Kuba K, Morita M, Chida S, Tezuka H, Hara H, et al. CXCL10-CXCR3 enhances the development of neutrophil-mediated fulminant lung injury of viral and nonviral origin. Am J Respir Crit Care Med (2013) 187:65–77. doi: 10.1164/rccm.201203-0508OC

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Lang S, Li L, Wang X, Sun J, Xue X, Xiao Y, et al. CXCL10/IP-10 neutralization can ameliorate lipopolysaccharide-induced acute respiratory distress syndrome in rats. PloS One (2017) 12:e0169100. doi: 10.1371/journal.pone.0169100

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zeng X, Moore TA, Newstead MW, Deng JC, Kunkel SL, Luster AD, et al. Interferon-inducible protein 10, but not monokine induced by gamma interferon, promotes protective type 1 immunity in murine klebsiella pneumoniae pneumonia. Infection Immun (2005) 73:8226–36. doi: 10.1128/IAI.73.12.8226-8236.2005

CrossRef Full Text | Google Scholar

38. Chen X, Yang B, Tian J, Hong H, Du Y, Li K, et al. Dental follicle stem cells ameliorate lipopolysaccharide-induced inflammation by secreting TGF-beta3 and TSP-1 to elicit macrophage M2 polarization. Cell Physiol Biochem (2018) 51:2290–308. doi: 10.1159/000495873

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Stifter SA, Bhattacharyya N, Pillay R, Florido M, Triccas JA, Britton WJ, et al. Functional interplay between type I and II interferons is essential to limit influenza a virus-induced tissue inflammation. PloS Pathog (2016) 12:e1005378. doi: 10.1371/journal.ppat.1005378

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Schmit T, Guo K, Tripathi JK, Wang Z, McGregor B, Klomp M, et al. Interferon-gamma promotes monocyte-mediated lung injury during influenza infection. Cell Rep (2022) 38:110456. doi: 10.1016/j.celrep.2022.110456

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Karki R, Sharma BR, Tuladhar S, Williams EP, Zalduondo L, Samir P, et al. Synergism of TNF-alpha and IFN-gamma triggers inflammatory cell death, tissue damage, and mortality in SARS-CoV-2 infection and cytokine shock syndromes. Cell (2021) 184:149–168 e17. doi: 10.1016/j.cell.2020.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Acute lung injury, cytokine storm, influenza A virus, methicillin-resistant Staphylococcus aureus, coinfection

Citation: Verma AK, McKelvey M, Uddin MB, Palani S, Niu M, Bauer C, Shao S and Sun K (2022) IFN-γ transforms the transcriptomic landscape and triggers myeloid cell hyperresponsiveness to cause lethal lung injury. Front. Immunol. 13:1011132. doi: 10.3389/fimmu.2022.1011132

Received: 03 August 2022; Accepted: 01 September 2022;
Published: 20 September 2022.

Edited by:

Yi Shi, Fudan University, China

Reviewed by:

Victor C. Huber, University of South Dakota, United States
Adam D. Hoppe, South Dakota State University, United States

Copyright © 2022 Verma, McKelvey, Uddin, Palani, Niu, Bauer, Shao and Sun. 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.

*Correspondence: Keer Sun, a2VzdW5AdXRtYi5lZHU=

These authors have contributed equally to this work and share first authorship

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