Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 02 November 2022
Sec. Inflammation

Integrating bulk and single-cell sequencing reveals the phenotype-associated cell subpopulations in sepsis-induced acute lung injury

Fuquan Wang,Fuquan Wang1,2Ming Chen,Ming Chen1,2Jiamin Ma,Jiamin Ma1,2Chenchen Wang,Chenchen Wang1,2Jingxu Wang,Jingxu Wang1,2Haifa Xia,Haifa Xia1,2Dingyu Zhang,,*Dingyu Zhang1,2,3*Shanglong Yao,*Shanglong Yao1,2*
  • 1Department of Anesthesiology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
  • 2Department of Anesthesiology, Institute of Anesthesia and Critical Care Medicine, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
  • 3Wuhan Jinyintan Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

The dysfunctional immune response and multiple organ injury in sepsis is a recurrent theme impacting prognosis and mortality, while the lung is the first organ invaded by sepsis. To systematically elucidate the transcriptomic changes in the main constituent cells of sepsis-injured lung tissue, we applied single-cell RNA sequencing to the lung tissue samples from septic and control mice and created a comprehensive cellular landscape with 25044 cells, including 11317 immune and 13727 non-immune cells. Sepsis alters the composition of all cellular compartments, particularly neutrophils, monocytes, T cells, endothelial, and fibroblasts populations. Our study firstly provides a single-cell view of cellular changes in septic lung injury. Furthermore, by integrating bulk sequencing data and single-cell data with the Scissors-method, we identified the cell subpopulations that are most associated with septic lung injury phenotype. The phenotypic-related cell subpopulations identified by Scissors-method were consistent with the cell subpopulations with significant composition changes. The function analysis of the differentially expressed genes (DEGs) and the cell-cell interaction analysis further reveal the important role of these phenotype-related subpopulations in septic lung injury. Our research provides a rich resource for understanding cellular changes and provides insights into the contributions of specific cell types to the biological processes that take place during sepsis-induced lung injury.

Introduction

Sepsis is an organ dysfunction characterized by high morbidity, severity, and mortality due to the dysregulated host response to infection (1), which imposes a huge burden on human health and socioeconomics (2). Although considerable research has been devoted to the pathogenesis and treatment of sepsis, the exact pathogenesis of sepsis has not been elucidated, and the clinical outcomes of sepsis patients are still unsatisfactory. The lung is the first and most commonly damaged organ in sepsis (3). Approximately 50% of patients with sepsis develop acute lung injury (ALI) (4). ALI and its more severe form, acute respiratory distress syndrome(ARDS), are the leading cause of death in sepsis patients (5).

Sepsis induced-ALI is characterized by the overactivated inflammatory response, disruption of the alveolar-capillary membrane barrier, pulmonary edema, and pulmonary dysfunction. Numerous immune cells and non-immune cells are involved in the occurrence and progression of septic lung injury. Therefore, although research contributing to the pathological mechanism of septic lung injury has been extensively studied for decades, there is currently a lack of effective enough understanding of septic lung injury due to the complex nature and the involvement of multiple cell types (6).

Sepsis-induced lung injury comprises a complex microenvironment of interacting lung structural constituent cells and immune cells. However, the current cellular-level studies on the pathological mechanism of septic lung injury are relatively independent, most of the studies only focus on one or two of the lung cell components such as epithelium, endothelium, monocytes, and neutrophils (3, 79), and lack a systematic explanation of the overall cellular changes during lung injury. Studies of septic lung injury that focus only on individual cell types are not enough. In addition, the current research on septic lung injury mainly uses peripheral blood or bronchoalveolar lavage fluid, and there is no systematic elucidation of the damage to the main constituent cells of lung tissue. Therefore, there is an urgent need to elucidate the comprehensive landscape and heterogeneity of immune and non-immune cells of the lung tissue during septic lung injury.

Single-cell RNA sequencing (scRNA-seq) is an efficient tool for exploring cell heterogeneity (10). Therefore, by using the well-established mouse model of sepsis-induced lung injury, we performed single-cell sequencing on mouse lung tissue to develop a deeper and more accurate understanding of the changes and interactions of various cell types in the lung tissue microenvironment during septic lung injury (The main flow of the experiment is shown in Figure 1A). The single-cell RNA sequencing is capable of distinguishing cell types, states, and lineages well. However, the current single-cell data cannot directly link cell clusters to specific phenotypes. The Scissors-method developed by Xia Zheng et al. can identify cell subsets related to a given phenotype by integrating phenotype-associated bulk expression data and single-cell data (11). So we also performed bulk RNA sequencing and used the Scissors-method to identify the cell subpopulations most associated with septic lung injury.

FIGURE 1
www.frontiersin.org

Figure 1 Identifying the infiltrated cell types in septic and non-septic lung tissue. (A) The main workflow of the study. (B) The He staining results of the lung tissue (magnification 200×). (C) The UMAP plot of all scRNA-seq data, shows a total of 14 distinct cell types that were identified. (D) UMAP plot of sequenced cells grouped by color by immune and non-immune cells. (E) Dot plot for cell-type-specific signature genes. (F) The number of immune cells and non-immune cells in each group. (G) The proportion of various types of cells in different groups of immune cells. (H) The proportion of various types of cells in different groups of non-immune cells.

Here, we report extensive profiling of the cellular composition in the sepic mouse lung by generating scRNA-seq profiles of 25044 cells. We provide a detailed cell atlas of normal and septic mouse lung tissue. We also identified the cell subpopulations closely associated with the septic lung injury phenotype with Scissors-method and revealed the specific genes, gene function enrichment, and interactions with other cell clusters for these subpopulations. Our results provide insights into sepsis-associated lung injury and provide an essential resource for the future discovery of effective treatment strategies.

Methods

Ethics statement

The experiment was performed in accordance with the Chinese Animal Research Guidelines and was approved by the Laboratory Animal Management Committee of Tongji Medical College of Huazhong University of Science and Technology.

Animals

Six weeks old male C57BL/6 mice (weight, 20–25 g) were housed in the standard laboratory at 22°C with 40–60% and 12 h/12 day/night cycle. The mice had free access to food and water. The mice were allowed to adapt to the environment 1 week before the experiment. The mice were randomly assigned to three groups: sham group, sepsis group (including CLP group, and CLP-48h group).

We used the Cecal ligation and puncture (CLP), the most commonly used method to induce polymicrobial sepsis, to construct the mouse model of septic lung injury. CLP was performed as described previously (12). The mice in the Sham group and CLP group were sacrificed 24 hours after the operation, while the CLP-48h group was sacrificed 48 hours after the operation. The lung tissue samples were collected for subsequent testing and analysis. We performed bulk RNA sequencing on 3 Sham group specimens and 3 CLP group specimens. Three sham group specimens, three CLP specimens, and three CLP-48h specimens were mixed for single-cell RNA sequencing.

Hematoxylin and eosin staining

Lung tissue was fixed with 4% paraformaldehyde and embedded in paraffin. The paraffin was then cut into 5 µm sections, deparaffinized, rehydrated, and stained with H&E. The pathological changes in lung tissue were observed under a light microscope.

RNA sequencing

The total lung RNAs were used for RNA-seq. The RNA sequencing libraries were constructed and sequenced on the BGISEQ-500 platforms.

Single-cell suspension preparation

After being isolated from the mice, the fresh lung tissue was rinsed twice with the pre-cooled RPMI 1640 + 0.04% BSA medium under sterile conditions. Then the samples were sufficiently cut into small pieces of about 0.5 mm3. Then the samples were incubated with a mixture of 2 ml of trypsin (Gibco, Cat: R001100), 1 ml of collagenase II (Biofrox, Cat:2275MG100), and 100 µl of DNase (Servicebio, Cat: 1121MG010) for 1 h on a 37°C shaker. The digested cell suspension was filtered 1-2 times with a BD 40 μm cell sieve and centrifuged at 300 g for 5 min at 4°C. After the pellet was resuspended with an appropriate amount of medium, an equal volume of red blood cell lysate (MACS, Cat. No. 130-094-183) was added. After mixing, let stand at 4°C for 10 min, the cell suspension was centrifuged at 300g for 5 min, and the supernatant was discarded. The pellet was washed once with medium and centrifuged at 300 g for 5 min, the supernatant was discarded, the cell pellet was resuspended with 100 μl of the medium, and the cell concentration and viability were calculated by Luna cell counter.

For CD45 sorting, the prepared single-cell suspensions were subjected to Magnetic-activated cell sorting (MACS) using Miltenyi mouse CD45 magnetic beads according to the manufacturer’s recommendations. The CD45- and CD45+ cell suspensions were mixed 1:1 for single-cell sequencing.

After adjusting the cell concentration of the single-cell suspension to 700-1200 cells/μl, the single-cell library was prepared following the protocol of the 10×Genomics Chromium Next GEM Single Cell 3′ Reagent Kits v3.1. High-throughput sequencing was performed on the constructed library using the Illumina Nova 6000 PE150 platform.

Single-cell data preprocessing

Cell Ranger (version 5.0.0) was used to process the raw data and generate the raw unique molecular identifier(UMI) count matrix. The merged matrix was transferred into the R statistical environment for further analysis using the Seurat package (v. 4.1.2).

The cells with UMI counts < 1000, mitochondrial UMI counts greater than 10%, or below 500 genes were identified as low-quality cells and removed (Supplementary Figure 1A, All Supplementary Figures are in the “Supplementary Figures” file). The gene expression matrix was normalized using the “NormalizeData” function and 2000 highly variable genes (HVGs) were identified by the “FindVariableFeatures” function. Batch effects were corrected by the “Harmony” function. Principal component analysis (PCA) was performed to reduce the dimension of the scRNA-Seq dataset. By the “ElbowPlot” function in Seurat, we used TOP 50 PCs for downstream analysis (Supplementary Figure 1B).

The cells were clustered and identified by the “FindNeighbors” and “FindClusters” functions with a resolution of 1. The clusters were projected into a two-dimensional plot for visualization using the “RunTSNE” and “RunUMAP” functions.

Cluster marker identification and cell-type annotation

The differentially expressed genes(DEGs) of each cluster were identified by the “FindAllMarkers” function in Seurat, the clusters were annotated according to the expression of classic markers. When identifying the main cell types, the cells that expressed more than one typical cell type marker were considered double cells and excluded from further analysis. Firstly, all cells were divided into 14 major cell types, and then the major cell types were further clustered and annotated separately to detect heterogeneity within each cell type.

“Scissor” method identifies the cell subpopulations associated with sepsis-induced lung injury

“Scissor” integrates the single-cell data and phenotype-relevant bulk sequencing data by first quantifying the similarity between each bulk and each cell sample. The scissor method can then optimize the regression model on the matrix associated with the sample phenotype to identify relevant subpopulations. A detailed description can be found in the literature (11).

Pseudotime analysis

The “Monocle 3” package (13) was used to infer the potential lineage differentiation trajectory.

Functional enrichment analysis

GO enrichment of DEGs of cell subgroups was performed using the “compareCluster” function in the “clusterProfiler” packages (14).

Cell-cell communication analysis

CellPhoneDB2, a Python-based computational analysis tool developed by Roser Vento-Tormo et al. (15), was used to analyze cell-to-cell communication at the molecular level.

Results

Detailed map of cellular composition in the sepsis-induced lung injury

The condition of the lung tissue of the mice was evaluated by HE staining and the results showed that the mice in the CLP group and the CLP-48h group had obvious lung injury (Figure 1B). To create a comprehensive cellular map of the normal and septic lung, we first isolated and sequenced a total of 31532 cells from lung tissue suspensions derived from three groups, including the Sham group (8742 cells), CLP group (12159 cells), CLP-48h group (10631 cells).

After standard data processing and the low-quality cells were filtered out, we performed downstream analysis on a total of 25044 lung cells, including 6825 cells from the Sham group, 9914 cells from the CLP group, and 8305 cells from the CLP-8h group. Based on the expression of the classic markers genes (Figure 1E; Supplementary Figure 1C), we clustered and annotated these cells into 14 major cell types (Figure 1C), including Cycle cell, Alveolar Macrophages, Mono cell, Macrophages, DC cell, neutrophil cell, B cell, NK cell, T cell, endothelial cell, epithelial cell, Fibroblast cell, smooth muscle cell, pericyte cell.

To reveal the cellular composition of normal lungs and septic lungs, we first divided the cells into immune cells and non-immune cells according to the expression of CD45 (Figures 1D, F; Supplementary Figure 1C), and the results showed that the numbers of immune cells and non-immune cells were in the same order of magnitude, indicating that magnetic bead sorting was effective.

To reveal the cellular composition of immune and non-immune cells, we compared the relative percentages of each cell type. Among immune cells (Figure 1G), neutrophil, mono, and T cells were the most three prevalent subsets, while in the non-immune cell fraction (Figure 1H), endothelial cells and fibroblast cells are relatively dominant. Next, we mainly characterized the transcriptomic features of the main cell clusters.

The heterogeneity analysis of monocytes

Monocytes represent the first line of defense against infection. In the early stages of infection and sepsis, monocytes undergo activation, and functional and morphological changes following infectious stimulation, showing a high degree of heterogeneity. Based on the differences in marker genes (Figure 2B), we clustered the monocytes into six subgroups (Figure 2A), including Alveolar Macrophages, Mono0, MacroCxcl3, MacroMrc1, MacroGrk3, and Mono1.

FIGURE 2
www.frontiersin.org

Figure 2 Cellular composition of monocytes in normal and septic lung tissue. (A) a total of five clusters of monocytes were identified in lung tissue. Cell subpopulations are colored as shown in the legend. (B) Dot plot for principal identifiers of different subclusters. (C) The proportions of different monocyte subclusters. (D) The UMAP visualization of the Scissor-medthod selected cells. The red dots are cells associated with the sepsis phenotypes. (E) The prediction of monocyte differentiation trajectories with Monocle. (F) Heatmap of DEGs of the Mono1 subcluster. (G) Enriched GO functions of upregulated genes in Mono1.

We observed that Mono1 was specifically identified in septic lung injury, the proportion of which in the CLP group was higher than that in the CLP-48h group (Figure 2C). The scissors method suggested that Mono1 cluster were the subclusters that were most associated with the sepsis-induced lung injury phenotype (Figure 2D). Of note, the percentage of Mono0 was decreased in the CLP group compared with the Sham group, which was even lower than the CLP-48h group, while the proportion of Macro-Cxcl3 showed an obvious opposite trend. To determine if there is a link between Mono0 and Mono1, we performed the pseudotime analysis. The Monocle3 pseudotime analysis suggested (Figure 2E) that Mono0 was differentiated from Mono1 clusters. The mono1 population may be a specific cell cluster migrating from peripheral blood under the condition of septic lung injury. Interestingly, the proportion of Mono1 in the CLP-48h group is lower than that in the CLP group, indicating that Mono1 may be closely related to the acute phase of sepsis-induced lung injury. We further identified the DEGs of Mono1, as expected, the DEGs of Mono1 were shown to possess the characteristics of promoting inflammatory response, including S100A8/A9, Hmgb1, Mmp8, etc. (Figure 2F). The Go enrichment analyses of the Mono1 DEGs were predominantly involved in the inflammatory signals, including chemotaxis of neutrophils and lymphocytes, regulation of cytokine production, etc. (Figure 2G). Therefore, Mono1 is more likely to be the most important subgroup that aims to promote the inflammatory response in sepsis acute lung injury.

Among the DEGs, the most noteworthy gene is Hmgb1 and S100A8/A9. A considerable number of studies have been published to confirm the important roles of S100A8/A9 in inflammatory-related diseases (16, 17). Previous studies have indicated that the upregulation of S100A8/A9 monocytes may contribute to cytokine storms in severe patients (18). Our study shows that S100A8/A9+ monocytes are also closely associated with septic lung injury. It is worth noting that the expression level of S100A8/A9 in the CLP-48h group was higher than that in the CLP group, and began to be significantly expressed in the Mono0 subgroup (Supplementary Figure 2A).

As a ubiquitous nuclear protein that can be secreted by activated monocytes/macrophages, the role of HMGB1 in orchestrating inflammatory responses has been demonstrated in previous studies (19, 20). Furthermore, clinical evidence suggests that the levels of circulating Hmgb1 are positively associated with severity and mortality in septic patients (21). Compared with S100a8/a9, the expression of Hmgb1 showed more obvious subgroup specificity, our study showed that in septic lung injury, Hmgb1 is not indiscriminately elevated in all monocytes/macrophage populations, but is concentrated in the Mono1 subclusters, especially the Mono1 subpopulation in the CLP group.

We next investigated the network of cell-cell interactions among Mono1 and other immune cells identified in our current work by using CellphoneDB2.0 (Supplementary Figure 2B). Notably, CellphoneDB is only applicable to human genes, and it is necessary to convert mouse genes into human genes for relevant analysis. Mono1 expresses higher levels of CCR1, the receptors of which includes CCC15, CCL18, and CCL23 (The corresponding relationships in cellphone are: CCL15-Ccl6, CCL18-Ccl3, CCL23-Ccl6). Interestingly, CCC15, CCL18, and CCL23 of Mono1 cells were also highly expressed. In addition, the above-mentioned Ccr1 receptor-ligand pair is mainly expressed in some subsets of monocytes/macrophages and neutrophils. This suggests that the chemotaxis of Mono1 is regulated by other immune cells and that Mono1 is also responsible for the infiltration of certain other immune cells.

The heterogeneity analysis of neutrophils

Excessive accumulation and activation of neutrophils in lung tissue is the main pathological feature of septic lung injury (22). Previous studies have suggested that neutrophil populations are not homogenous (23). We identified four distinct neutrophil clusters (Figure 3A) based on their expression profiles (Figure 3B), including Neutrophil0-Ccl6, Neutrophil1-Krt83, Neutrophil2-Csta2, Neutrophil3-Cxcl3. The analysis of the composition of neutrophils in each group showed that the neutrophils in sepsis-related lung tissue were significantly different from the normal lung tissue (Figures 3C, D). Among them, the proportion of Neutrophil3-Cxcl3 subclusters in the CLP group was significantly increased.

FIGURE 3
www.frontiersin.org

Figure 3 Cellular composition of neutrophils in normal and septic lung tissue. (A) a total of four subclusters of neutrophils were identified in lung tissue. Cell subpopulations are colored as shown in the legend. (B) Dot plot for principal identifiers of different subclusters. (C) The proportions of different monocyte subclusters. (D) The prediction of neutrophils differentiation trajectories with Monocle. (E) The UMAP visualization of the Scissor-method selected cells. The red dots are cells associated with the sepsis phenotype. (F) Violin plots of DEGs for Neutrophil3-Cxcl3. (G) Enriched GO functions of upregulated genes in Neutrophil3-Cxcl3. (H) Dot plot show ligand-receptor pairs of cytokines between Neutrophil3-Cxcl3 and other immune cell groups.

To demonstrate the developmental progression among different neutrophil clusters described above, we ordered these clusters by performing the trajectory analysis using Monocle 3 (Figure 3E). Mme is an important marker of neutrophil maturation. So we also analyzed the expression of Mme in different neutrophil subpopulations. The results showed that although the overall expression level was very low, the expression level of Mme in the Neutrophil3-Cxcl3 subpopulation was still significantly lower than that of other subgroups(As shown in Supplementary Figure 3A). Our results suggested that the Neutrophil3-Cxcl3 subpopulation is more likely to be derived from peripheral blood, then differentiate into Neutrophil0-Ccl6, and then differentiate into Neutrophil1-Krt83 and Neutrophil2-Csta2. The scissors method suggested that almost all neutrophils in the CLP and CLP-48h groups were associated with the phenotype of septic lung injury. Combining the above findings, we consider the Neutrophil3-Cxcl3 subclusters to be the main driver of lung injury in the CLP group.

The highly expressed genes in Neutrophil3-Cxcl3 (Figure 3F), such as Cxcl2, Cxcl3, and Ccl3, have been reported to involve in neutrophil and other immune cell migration (2426). The GO analyses (Figure 3G) also suggested the DEGs in Neutrophil3-Cxcl3 are mainly enriched in the chemotaxis and migration of neutrophils, lymphocytes, monocytes, etc. So we next explored the cell-cell interaction network between Neutrophil3-Cxcl3 and other immune cell types(Figure). CellphoneDB results (Figure 3H) showed that Neutrophil3-Cxcl3 mainly expressed CCL18, the receptor of which was mainly expressed on certain neutrophils, monocytes/macrophages, and DC cell subsets. Neutrophil3-Cxcl3 also highly expresses the chemokine receptor Ccr1, the ligand of which is expressed in a variety of immune cells, indicating that the migration and infiltration of C during septic lung injury is regulated by a variety of immune cells. Our studies suggested that the chemotaxis of Neutrophil3-Cxcl3 may play a positive feedback effect on inflammatory cell recruitment. Furthermore, the major chemokines and chemokine receptors of Neutrophil3-Cxcl3, such as Ccl3, and Cxcr1, may be potential main intervention targets for inhibiting neutrophil-sustained inflammatory conditions in septic lung injury.

The heterogeneity analysis of lymphocytes

We identified 5 subgroups of T cells based on the expression of classic T cell markers (Figure 4B), including “CD8 Tcell”, “Th17 cell”, “CD8 TCM”, “CD4 Tcell”, and “ILC2” (Figure 4A). We originally selected T lymphocytes for further analysis, but during the subpopulation annotation process, we found a subcluster that negatively expressed Cd3d, Cd4, and CD8a, but was positive for Gata3(Supplementary Figure 4A), which was identified as ILC2.

FIGURE 4
www.frontiersin.org

Figure 4 Cellular composition of lymphocytes in normal and septic lung tissue. (A) a total of five subclusters of lymphocytes were identified in lung tissue. Cell subpopulations are colored as shown in the legend. (B) Dot plot for principal identifiers of different subclusters. (C) The proportions of different lymphocyte subclusters. (D) The UMAP visualization of the Scissor-method selected cells. The red and blue dots are cells associated with sepsis and normal phenotypes (E) Heatmap of DEGs of the ILC2 subcluster. (F) Enriched GO functions of upregulated genes in ILC2. (G) Violin plot of the expression levels of Areg in different lymphocyte subclusters.

We observed that the proportion of CD8 Tcell was significantly reduced and the ILC2 cell was significantly increased in the sepsis group compared with the Sham group (Figure 4C). The scissors method suggested that ILC2 was the subgroup most associated with the sepsis phenotype (Figure 4D), so we focus on the ILC2 subclusters for further analysis. Considering that ILC2 is mainly expressed in the CLP group, we directly compared the DEGs between ILC2 and other subclusters. The top10 differential genes are shown in Figure 4E.

Previous studies have shown that ILC2 plays an important role in regulating adaptive immunity, macrophage polarization, etc (27, 28). The results of the Go enrichment analysis (Figure 4F) showed that DEGs of ILC2 were mainly focused on regulating innate immune cells, including regulating the myeloid cell differentiation, and positive regulation of cytokine (18). As a member of the epidermal growth factor (EGF)family, amphiregulin (Areg) has been implicated in exerting an important role in regulating tissue repair in acute epithelial injury and asthma (29, 30). ILC2 is an important source of Areg in the lung (31), however, the expression of Areg in ILC2 is significantly inhibited in septic lung injury (Figure 4G; Supplementary Figure 4B), which may be one of the important mechanisms of septic lung injury.

The heterogeneity analysis of endothelial cells

Inflammation and dysfunction of endothelial function play a central role in sepsis-induced lung injury. We identified endothelial cells (ECs) into six clusters of distinct expression profiles (Figure 5B), including Capillary ECs-Sema3c, Capillary ECs-Ednrb, lymphatic ECs, Arterial ECs, Vein ECs, Capillary ECs (Figure 5A).

FIGURE 5
www.frontiersin.org

Figure 5 Cellular composition of endothelial cells in normal and septic lung tissue (A) a total of six subclusters of endothelial cells were identified in lung tissue. Cell subpopulations are colored as shown in the legend. (B) Dot plot for principal identifiers of different subclusters. (C) The proportions of different lymphocyte subclusters. (D) The UMAP visualization of the Scissor-method selected cells. The red and blue dots are cells associated with sepsis and normal phenotypes. (E) Enriched GO functions upregulated genes of Capillary ECs-Sema3c and Capillary Ecs in the Sham group. (F) Enriched GO functions upregulated genes of lymphatic ECs in the sepsis group. (G) Heatmap of DEGs of Capillary ECs-Sema3c and Capillary ECs subpopulations in the Sham group. (H) Heatmap of DEGs of lymphatic ECs in the sepsis group.

The percentage of endothelial cell composition in each group showed the cell subsets with the most obvious percentage changes were Capillary ECs-Ednrb and Capillary ECs cells (Figure 5C). The scissors method showed that populations of Capillary ECs-Sema3c and Capillary Ecs are the workhorses of endothelial cells against septic lung injury, while the lymphatic ECs are closely related to the progression of sepsis(Figure 5D). So we selected the three subclusters mentioned above for further analysis. The top10 DEGs of the Capillary ECs-Sema3c and Capillary ECs of the Sham group are shown in Figure 5G. We found that the DEGs of Capillary ECs-Sema3c and Capillary ECs in the Sham group were mainly enriched in the function of suppressing inflammatory responses, including in the response to interferon γ, promoting the activation of T cells and inhibiting other immune cells (Figure 5E). Immune interferon-gamma (IFN-γ) plays an important role in fighting infection and exerts its effects through an extensive transcriptional program involving approximately 2000 genes (32, 33), of which the Guanylate-binding proteins (GBP) family is an important member. GBPs are a family of dynamin-related large GTPases which are expressed in response to interferons. The antibacterial function of GBPs has been confirmed by multiple studies (34, 35), consistent with the GO enrichment showing that the Sham group had an increased response to interferon-γ. The gene expression of GBP2, GBP6, and GBP7 in the Sham group was significantly higher than that in the sepsis group.

CD74, the major histocompatibility complex II-associated invariant chain, mainly expresses antigen-presenting cells (APC), which is related to the function of antigen presentation. Studies have also found that CD74 molecules are also expressed in endothelial cells (36) which are related to the occurrence and development of corresponding diseases (37). The expression of CD74 was significantly decreased in the sepsis-associated group, which might be associated with the excessive inflammatory response (38). Moreover, CD74 is an indispensable part of the receptor complex of macrophage migration inhibitory factor (MIF), and macrophages can produce IL-8 that chemotactic leukocytes to the site of infection (39). The expression of CD74 in the sepsis-related group was significantly inhibited, which may be one of the reasons for the excessive activation of inflammatory cells.

The scissor method showed that the ECs clusters most associated with the septic lung injury phenotype were lymphatic ECs. GO enrichment analysis of lymphatic ECs showed that the DECs (Figure 5H), were enriched in the regulation of inflammatory response and regulatory peptide synthesis, transport, secretion, and other related pathways (Figure 5F). Previous studies have shown that endothelial cells can be involved in sepsis progression by regulating the expression of related peptides, which is consistent with our findings.

Our study also found that lymphatic ECs DEGs were associated with the migration of granulocytes, so we next explored the cell-cell interactions of lymphatic ECs and other immune cells in the CLP group. CellphoneDB results (Supplementary Figure 5A) showed that the lymphatic ECs of the CLP group mainly expressed Ccl7 (Supplementary Figure 5B), the receptor of which includes Ccr1 and Ccr2, thereby exerting chemotactic effects on certain subsets of neutrophils and monocytes. Notably, there were significant interactions in the CLP group between lymphatic ECs and the subpopulations that are closely associated with the septic lung injury phenotype, such as Mono1, and Neutrophil-Cxcl3.

The heterogeneity analysis of fibroblasts

In the current study, our studies reveal the diversity of fibroblast cells with seven cell clusters (Figure 6A) based on differences in marker genes (Figure 6B). There were significant differences in the fibrocystic constituents of the lung tissue of normal mice and septic mice (Figure 6C).

FIGURE 6
www.frontiersin.org

Figure 6 Cellular composition of fibroblasts in normal and septic lung tissue. (A) a total of seven subclusters of fibroblasts were identified in lung tissue. Cell subpopulations are colored as shown in the legend. (B) Dot plot for principal identifiers of different subclusters. (C) The proportions of different fibroblasts subclusters. (D) The UMAP visualization of the Scissor-method selected cells. The red and blue dots are cells associated with sepsis and normal phenotypes. (E) The prediction of fibroblast differentiation trajectories with Monocle. (F) Heatmap of top10 DEGs of different groups of fibroblasts. (G) Dot plot show ligand-receptor pairs of cytokines between fibroblasts and other immune cell groups in the CLP group.

The Scissors-method showed that most of the fibroblasts in the CLP group were not innocent of septic lung injury, especially the constituent cells of the Fib0-Saa1 and Fib2-Ces1d subpopulations (Figure 6D). In addition, the Fib5-Dmkn cluster was only present in the sepsis-related group. The results of Monocle3 pseudotime analysis (Figure 6E) showed that Fib0-Saa1 and Fib2Ces1d subpopulations are differentiated from Fib1-Pi16.

We directly performed differential gene analysis on the fibroblasts of the three groups, and the top 10 differential genes were shown in Figure 6F. The results of the Go enrichment analysis (Supplementary Figure 6A) showed that the differential genes of fibroblasts in the CLP group were mainly concentrated in the acute inflammatory response, lipoprotein-related pathways, and reactive oxygen pathways. The Serum-Amyloid A(SAA) family protein genes SAA1 and SAA3 were the main players involved in the GO pathway. As one of the typical positive acute phase response proteins, SAA is closely related to infectious diseases. In the process of inflammation, plasma concentrations of SAA (mainly SAA1) increase dramatically by up to 1,000-fold (40). SAA could serve as a major diagnostic marker in the acute phase of inflammation in animals and humans. The chemotactic properties of SAA protein for monocytes and neutrophils have been confirmed (41, 42). We found that the expression of SAA1 and SAA3 was significantly increased in the CLP group, especially in the Fib0-Saa1 and Fib2-Ces1d subpopulations(Supplementary Figure 6B).

Previous studies have shown that damage-responsive fibroblasts can alter the lung microenvironment to promote robust infiltration of immune cells at the expense of lung function (43). So we analyzed the interaction between fibroblasts and immune cells in the CLP group. CellphoneDB results (Figure 6G) showed that the fibroblast of the CLP group mainly expressed Ccl7 (Supplementary Figure 6C). Notably, fibroblasts in the CLP group, similar to lymphatic ECs, both exhibited significant chemotaxis to cell subclusters associated with the septic lung injury phenotype, including Mono1 and Neutrophil-Cxcl3.

Discussion

Despite intense research efforts aimed at exploring the pathogenesis and treatment of sepsis-induced lung injury, there is still a lack of effective means to address this clinical challenge. Surviving severe sepsis-induced lung injury depends on a careful balance between building an immune response sufficient to clear the infection and maintaining lung function in the presence of immune-induced tissue damage. There is an urgent need to clarify the precise pathogenic mechanisms of sepsis-induced ALI to identify new therapeutic targets.

The host inflammatory response dominated by immune cells can effectively limit the spread of pathogens and eventually clear the pathogens, but immunopathology is also a major factor in tissue damage and ARDS (5, 44). Non-immune lung cells, including fibroblasts, endothelial cells, and epithelial cells are also involved in coordinating immune responses and are directly involved in the regulation of lung function (45). Therefore, the coordination of the functions of immune cells and non-immune cells of the system is the key to solving the lung injury caused by sepsis. However, there is currently no systematic study of immune and non-immune cells of the lung cell lineages in the context of sepsis-induced lung injury. Therefore, we depicted the high-resolution visualization of the cellular landscape in the septic lung via scRNA-seq technology, which will facilitate a better understanding of the pathogenesis of sepsis-induced lung injury.

We firstly provide extensive profiling of cellular composition in sepsis-induced lung injury by using multiplexed scRNA-seq to assess the expression profile of 25044 cells. Monocytes, neutrophils, T cells, endothelial cells, and fibroblasts were then selected as major cellular components for further analysis. Our study provides insight into the pathogenesis of sepsis-induced lung injury by characterizing several pathological cell subpopulations with distinct molecular expression profiles.

We used the Scissors-method to identify the cell subpopulations most associated with the septic lung injury phenotype. During the analysis of the major subpopulations, we found that the most phenotypically correlated cell subclusters identified by the scissors method were consistent with the subclusters, the proportions of which varied significantly across the three groups. This also proves to a certain extent that the cell subsets we focus on are closely related to septic lung injury.

It is generally believed that inflammatory storm is an important pathological mechanism of acute lung injury in sepsis, and innate immune cells play a major role in this process. As the first line of defense against the invasion of foreign microorganisms, neutrophils were once considered to be the same type of terminally differentiated white blood cells, with single functional characteristics, short lifespan, and low transcriptional activity. However, recent studies have shown that neutrophils have phenotypic and functional heterogeneity under different physiological or pathological conditions (46). Neutrophils and monocytes are key components of the innate immune system, the excessive activation of which is known to cause organ dysfunction in sepsis (47). However, our study shows that not all immune cells are murderers. In the case of monocytes, the Scissors-method suggested that only the Mono0 and Mono1 subclusters were strongly associated with the septic lung injury phenotype. While in the neutrophils, nearly all of them showed phenotype-related transcriptional changes during septic lung injury.

Many past studies have attempted to address acute lung injury from sepsis by focusing on attenuating epithelial and endothelial cell damage (22, 48, 49). However, not all non-immune cells, including epithelial and endothelial cells, which are more often considered “victim” roles, are innocent during septic lung injury. Previous studies suggested that the physiological function of endothelial cells is severely disrupted during sepsis, and not only functions as a physical barrier (50, 51) but actively participates in the regulation of physiological and pathophysiological processes (52). Our study shows that in the process of septic lung injury, mainly lymphatic endothelium plays a negative regulatory role in an endothelial cell population, and in fibroblasts, mainly Fib0-Saa1 and Fib2-Ces1d are involved in the lung injury process. These non-immune cell subsets also cannot be ignored in the exploration of septic lung injury.

No single subpopulation can function alone during septic lung injury, so we also performed cellular communication analysis for subpopulations with significant heterogeneity, primarily those associated with the lung injury phenotype. Our results suggested that cell subclusters more associated with the septic lung injury phenotype, such as neutrophil3-Cxcl3 and Mono1 subclusters, are more strongly associated with endothelial cells, fibroblasts, and many other immune cells. These cells not only expressed high levels of chemokine receptors but also highly expressed chemokine ligands, which may be associated with the generation of inflammatory storms.

Our study also revealed many possible therapeutic targets, such as our study showed that ILC2 highly expressed Areg, while the expression of Areg of ILC2 in the CLP group was significantly inhibited, which may be an intervention target for promoting airway epithelial integrity and lung tissue homeostasis. For another example, both lymphocytes and endothelial cells associated with the phenotype of septic lung injury highly express the chemokine Ccl7, which plays a chemotactic role for some phenotype-related immune cells. This suggested that Ccl7 in non-immune cells may also be a clue for the intervention of septic lung injury.

There are several limitations in our study, for example, the total number of cells involved in our study is not very large, resulting in the annotation of some cell subclusters may not be detailed enough. Furthermore, because our study was more focused on the comprehensive landscape of immune and non-immune cells in lung tissue during septic lung injury, mainly at 24 h, deep functional validation of individual subpopulations was lacking. In our opinion, the effective functional verification of the identified subpopulation most related to the sepsis lung injury phenotype is to use gene knockout animals, preferably lung tissue conditional knockout mice. Knock out specific functional genes of relevant subsets respectively to explore whether the symptoms and prognosis of sepsis-induced lung injury can be improved. We will continue to explore relevant work in depth in the future. Nevertheless, our study may provide foundations for further exploration of the mechanism of septic lung injury to a certain extent.

In conclusion, our study is the first to systematically provide a cellular transcriptional landscape of lung tissue in sepsis-induced lung injury. We also identified the subpopulations of cells most associated with the septic lung injury phenotype and revealed the transcriptional signature of these subpopulations. Our study will advance a better understanding of the pathogenic roles of different lung cell clusters during sepsis-induced lung injury.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE207651.

Ethics statement

The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Huazhong University of Science and Technology(IACUC Number:2629).

Author contributions

FW, MC, and JW established the animal models and collected samples. FW, JM, and CW finished the bioinformatics analysis. JW and HX prepared the figures. FW wrote this article. SY. and DZ designed the experiments and revised the article. All authors contributed to the article and approved the submitted version.

Funding

The study was supported by the Major Technological Innovation Special Project of Hubei Province of China (2019ACA167).

Acknowledgments

We would like to thank OE Biotech Co., Ltd (Shanghai, China) for providing single-cell RNA-seq, Dr. Junyi Hu (Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology) for the guidance of bioinformatics analysis.

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.981784/full#supplementary-material

References

1. Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA (2016) 315(8):801–10. doi: 10.1001/jama.2016.0287

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Rudd KE, Johnson SC, Agesa KM, Shackelford KA, Tsoi D, Kievlan DR, et al. Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the global burden of disease study. Lancet (London England) (2020) 395(10219):200–11. doi: 10.1016/S0140-6736(19)32989-7

CrossRef Full Text | Google Scholar

3. Park I, Kim M, Choe K, Song E, Seo H, Hwang Y, et al. Neutrophils disturb pulmonary microcirculation in sepsis-induced acute lung injury. Eur Respir J (2019) 53(3):1800786. doi: 10.1183/13993003.00786-2018

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Sevransky JE, Levy MM, Marini JJ. Mechanical ventilation in sepsis-induced acute lung injury/acute respiratory distress syndrome: an evidence-based review. Crit Care Med (2004) 32(11 Suppl):S548–53. doi: 10.1097/01.CCM.0000145947.19077.25

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

6. De Freitas Caires N, Gaudet A, Portier L, Tsicopoulos A, Mathieu D, Lassalle P. Endocan, sepsis, pneumonia, and acute respiratory distress syndrome. Crit Care (London England) (2018) 22(1):280. doi: 10.1186/s13054-018-2222-7

CrossRef Full Text | Google Scholar

7. He R, Liu B, Xiong R, Geng B, Meng H, Lin W, et al. Itaconate inhibits ferroptosis of macrophage via Nrf2 pathways against sepsis-induced acute lung injury. Cell Death discov (2022) 8(1):43. doi: 10.1038/s41420-021-00807-3

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Gao M, Yu T, Liu D, Shi Y, Yang P, Zhang J, et al. Sepsis plasma-derived exosomal miR-1-3p induces endothelial cell dysfunction by targeting SERP1. Clin Sci (London Engl 1979) (2021) 135(2):347–65. doi: 10.1042/CS20200573

CrossRef Full Text | Google Scholar

9. Xu B, Wang H, Chen Z. Puerarin inhibits ferroptosis and inflammation of lung injury caused by sepsis in LPS induced lung epithelial cells. Front Pediatr (2021) 9:706327. doi: 10.3389/fped.2021.706327

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Stubbington MJT, Rozenblatt-Rosen O, Regev A, Teichmann SA. Single-cell transcriptomics to explore the immune system in health and disease. Science (2017) 358(6359):58–63. doi: 10.1126/science.aan6828

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sun D, Guan X, Moran AE, Wu LY, Qian DZ, Schedin P, et al. Identifying phenotype-associated subpopulations by integrating bulk and single-cell sequencing data. Nat Biotechnol (2022) 40(4):527–38. doi: 10.1038/s41587-021-01091-3

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sun S, Wang J, Wang J, Wang F, Yao S, Xia H. Maresin 1 mitigates sepsis-associated acute kidney injury in mice via inhibition of the NF-κB/STAT3/MAPK pathways. Front Pharmacol (2019) 10:1323. doi: 10.3389/fphar.2019.01323

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Cao J, Spielmann M, Qiu X, Huang X, Ibrahim DM, Hill AJ, et al. The single-cell transcriptional landscape of mammalian organogenesis. Nature (2019) 566(7745):496–502. doi: 10.1038/s41586-019-0969-x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an r package for comparing biological themes among gene clusters. Omics J Integr Biol (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

CrossRef Full Text | Google Scholar

15. Vento-Tormo R, Efremova M, Botting RA, Turco MY, Vento-Tormo M, Meyer KB, et al. Single-cell reconstruction of the early maternal-fetal interface in humans. Nature (2018) 563(7731):347–53. doi: 10.1038/s41586-018-0698-6

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Marta-Enguita J, Navarro-Oviedo M, Rubio-Baines I, Aymerich N, Herrera M, Zandio B, et al. Association of calprotectin with other inflammatory parameters in the prediction of mortality for ischemic stroke. J Neuroinflamm (2021) 18(1):3. doi: 10.1186/s12974-020-02047-1

CrossRef Full Text | Google Scholar

17. de Korte AM, Aquarius R, Vogl T, Roth J, Bartels R, Boogaarts HD, et al. Elevation of inflammatory S100A8/S100A9 complexes in intracranial aneurysms. J neurointerv Surg (2020) 12(11):1117–21. doi: 10.1136/neurintsurg-2019-015753

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ren X, Wen W, Fan X, Hou W, Su B, Cai P, et al. COVID-19 immune features revealed by a large-scale single-cell transcriptome atlas. Cell (2021) 184(7):1895–913.e19. doi: 10.1016/j.cell.2021.01.05310.1016/j.cell.2021.01.053

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yang H, Wang H, Andersson U. Targeting inflammation driven by HMGB1. Front Immunol (2020) 11:484. doi: 10.3389/fimmu.2020.00484

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Deng M, Tang Y, Li W, Wang X, Zhang R, Zhang X, et al. The endotoxin delivery protein HMGB1 mediates caspase-11-Dependent lethality in sepsis. Immunity (2018) 49(4):740–53.e7. doi: 10.1016/j.immuni.2018.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Karlsson S, Pettilä V, Tenhunen J, Laru-Sompa R, Hynninen M, Ruokonen E. HMGB1 as a predictor of organ dysfunction and outcome in patients with severe sepsis. Intensive Care Med (2008) 34(6):1046–53. doi: 10.1007/s00134-008-1032-9

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Matthay MA, Zemans RL. The acute respiratory distress syndrome: pathogenesis and treatment. Annu Rev Pathol (2011) 6:147–63. doi: 10.1146/annurev-pathol-011110-130158

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Xie X, Shi Q, Wu P, Zhang X, Kambara H, Su J, et al. Single-cell transcriptome profiling reveals neutrophil heterogeneity in homeostasis and infection. Nat Immunol (2020) 21(9):1119–33. doi: 10.1038/s41590-020-0736-z

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Breda LCD, Breda CNS, de Almeida JRF, Paulo LNM, Jannuzzi GP, Menezes IG, et al. Fonsecaea pedrosoi conidia and hyphae activate neutrophils distinctly: Requirement of TLR-2 and TLR-4 in neutrophil effector functions. Front Immunol (2020) 11:540064. doi: 10.3389/fimmu.2020.540064

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Muntjewerff EM, Parv K, Mahata SK, van Riessen NK, Phillipson M, Christoffersson G, et al. The anti-inflammatory peptide catestatin blocks chemotaxis. J leukocyte Biol (2021) 112(2):273–8. doi: 10.1101/2020.11.23.393934

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sokulsky LA, Garcia-Netto K, Nguyen TH, Girkin JLN, Collison A, Mattes J, et al. A critical role for the CXCL3/CXCL5/CXCR2 neutrophilic chemotactic axis in the regulation of type 2 responses in a model of rhinoviral-induced asthma exacerbation. J Immunol (Baltimore Md 1950) (2020) 205(9):2468–78. doi: 10.4049/jimmunol.1901350

CrossRef Full Text | Google Scholar

27. Akama Y, Satoh-Takayama N, Kawamoto E, Ito A, Gaowa A, Park EJ, et al. The role of innate lymphoid cells in the regulation of immune homeostasis in sepsis-mediated lung inflammation. Diagnostics (Basel Switzerland) (2020) 10(10):808. doi: 10.3390/diagnostics10100808

CrossRef Full Text | Google Scholar

28. Nascimento DC, Melo PH, Piñeros AR, Ferreira RG, Colón DF, Donate PB, et al. IL-33 contributes to sepsis-induced long-term immunosuppression by expanding the regulatory T cell population. Nat Commun (2017) 8:14919. doi: 10.1038/ncomms14919

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Fukumoto J, Harada C, Kawaguchi T, Suetsugu S, Maeyama T, Inoshima I, et al. Amphiregulin attenuates bleomycin-induced pneumopathy in mice. Am J Physiol Lung Cell Mol Physiol (2010) 298(2):L131–8. doi: 10.1152/ajplung.90576.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Enomoto Y, Orihara K, Takamasu T, Matsuda A, Gon Y, Saito H, et al. Tissue remodeling induced by hypersecreted epidermal growth factor and amphiregulin in the airway after an acute asthma attack. J Allergy Clin Immunol (2009) 124(5):913–20.e1-7. doi: 10.1016/j.jaci.2009.08.044

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Monticelli LA, Sonnenberg GF, Abt MC, Alenghat T, Ziegler CG, Doering TA, et al. Innate lymphoid cells promote lung-tissue homeostasis after infection with influenza virus. Nat Immunol (2011) 12(11):1045–54. doi: 10.1038/ni.2131

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhang SY, Boisson-Dupuis S, Chapgier A, Yang K, Bustamante J, Puel A, et al. Inborn errors of interferon (IFN)-mediated immunity in humans: insights into the respective roles of IFN-alpha/beta, IFN-gamma, and IFN-lambda in host defense. Immunol Rev (2008) 226:29–40. doi: 10.1111/j.1600-065X.2008.00698.x

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Martens S, Howard J. The interferon-inducible GTPases. Annu Rev Cell Dev Biol (2006) 22:559–89. doi: 10.1146/annurev.cellbio.22.010305.104619

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Praefcke GJK. Regulation of innate immune functions by guanylate-binding proteins. Int J Med Microbiol IJMM (2018) 308(1):237–45. doi: 10.1016/j.ijmm.2017.10.013

CrossRef Full Text | Google Scholar

35. Kim BH, Shenoy AR, Kumar P, Das R, Tiwari S, MacMicking JD. A family of IFN-γ-inducible 65-kD GTPases protects against bacterial infection. Sci (New York NY) (2011) 332(6030):717–21. doi: 10.1126/science.1201711

CrossRef Full Text | Google Scholar

36. Rodor J, Chen SH, Scanlon JP, Monteiro JP, Caudrillier A, Sweta S, et al. Single-cell RNA-seq profiling of mouse endothelial cells in response to pulmonary arterial hypertension. Cardiovasc Res (2021) 118(11):2519–34. doi: 10.1093/cvr/cvab296

CrossRef Full Text | Google Scholar

37. Stumptner-Cuvelette P, Benaroch P. Multiple roles of the invariant chain in MHC class II function. Biochim Biophys Acta (2002) 1542(1-3):1–13. doi: 10.1016/S0167-4889(01)00166-5

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Baños-Hernández CJ, Bucala R, Hernández-Bello J, Navarro-Zarza JE, Villanueva-Pérez MA, Godínez-Rubí M, et al. Expression of macrophage migration inhibitory factor and its receptor CD74 in systemic sclerosis. Central-European J Immunol (2021) 46(3):375–83. doi: 10.5114/ceji.2021.109756

CrossRef Full Text | Google Scholar

39. Shi X, Leng L, Wang T, Wang W, Du X, Li J, et al. CD44 is the signaling component of the macrophage migration inhibitory factor-CD74 receptor complex. Immunity (2006) 25(4):595–606. doi: 10.1016/j.immuni.2006.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Obici L, Merlini G. AA Amyloidosis: basic knowledge, unmet needs and future treatments. Swiss Med weekly (2012) 142:w13580. doi: 10.4414/smw.2012.13580

CrossRef Full Text | Google Scholar

41. Baranova IN, Vishnyakova TG, Bocharov AV, Kurlander R, Chen Z, Kimelman ML, et al. Serum amyloid a binding to CLA-1 (CD36 and LIMPII analogous-1) mediates serum amyloid a protein-induced activation of ERK1/2 and p38 mitogen-activated protein kinases. J Biol Chem (2005) 280(9):8031–40. doi: 10.1074/jbc.M405009200

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Herold K, Moser B, Chen Y, Zeng S, Yan SF, Ramasamy R, et al. Receptor for advanced glycation end products (RAGE) in a dash to the rescue: inflammatory signals gone awry in the primal response to stress. J leukocyte Biol (2007) 82(2):204–12. doi: 10.1189/jlb.1206751

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Boyd DF, Allen EK, Randolph AG, Guo XJ, Weng Y, Sanders CJ, et al. Exuberant fibroblast activity compromises lung function via ADAMTS4. Nature (2020) 587(7834):466–71. doi: 10.1038/s41586-020-2877-5

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Soares MP, Teixeira L, Moita LF. Disease tolerance and immunity in host protection against infection. Nat Rev Immunol (2017) 17(2):83–96. doi: 10.1038/nri.2016.136

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Iwasaki A, Pillai PS. Innate immunity to influenza virus infection. Nat Rev Immunol (2014) 14(5):315–28. doi: 10.1038/nri3665

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Beyrau M, Bodkin JV, Nourshargh S. Neutrophil heterogeneity in health and disease: a revitalized avenue in inflammation and immunity. Open Biol (2012) 2(11):120134. doi: 10.1098/rsob.120134

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Cohen J. The immunopathogenesis of sepsis. Nature (2002) 420(6917):885–91. doi: 10.1038/nature01326

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Cheng KT, Xiong S, Ye Z, Hong Z, Di A, Tsang KM, et al. Caspase-11-mediated endothelial pyroptosis underlies endotoxemia-induced lung injury. J Clin Invest (2017) 127(11):4124–35. doi: 10.1172/JCI94495

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Li X, Jamal M, Guo P, Jin Z, Zheng F, Song X, et al. Irisin alleviates pulmonary epithelial barrier dysfunction in sepsis-induced acute lung injury via activation of AMPK/SIRT1 pathways. Biomedicine pharmacother = Biomed pharmacotherapie (2019) 118:109363. doi: 10.1016/j.biopha.2019.109363

CrossRef Full Text | Google Scholar

50. Lee WL, Slutsky AS. Sepsis and endothelial permeability. New Engl J Med (2010) 363(7):689–91. doi: 10.1056/NEJMcibr1007320

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Opal SM, van der Poll T. Endothelial barrier dysfunction in septic shock. J Internal Med (2015) 277(3):277–93. doi: 10.1111/joim.12331

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Ince C, Mayeux PR, Nguyen T, Gomez H, Kellum JA, Ospina-Tascón GA, et al. The endothelium in sepsis. Shock (Augusta Ga) (2016) 45(3):259–70. doi: 10.1097/SHK.0000000000000473

CrossRef Full Text | Google Scholar

Keywords: sepsis, lung injury, single-cell sequencing, Scissors-method, cellular landscape

Citation: Wang F, Chen M, Ma J, Wang C, Wang J, Xia H, Zhang D and Yao S (2022) Integrating bulk and single-cell sequencing reveals the phenotype-associated cell subpopulations in sepsis-induced acute lung injury. Front. Immunol. 13:981784. doi: 10.3389/fimmu.2022.981784

Received: 30 June 2022; Accepted: 21 October 2022;
Published: 02 November 2022.

Edited by:

Yongqing Li, School of Medicine, University of Michigan, United States

Reviewed by:

Lu Liu, Affiliated Hospital of Jiangsu University, China
Xianglin Meng, First Affiliated Hospital of Harbin Medical University, China

Copyright © 2022 Wang, Chen, Ma, Wang, Wang, Xia, Zhang and Yao. 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: Dingyu Zhang, 1813886398@qq.com; Shanglong Yao, shanglong_yao@sina.com

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.