PPARγ Attenuates Cellular Senescence of Alveolar Macrophages in Asthma- COPD Overlap

PPARγ Attenuates Cellular Senescence of Alveolar Macrophages in Asthma- COPD Overlap | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article PPARγ Attenuates Cellular Senescence of Alveolar Macrophages in Asthma- COPD Overlap Rongjun Wan, Prakhyath Srikaram, Shaobing Xie, Qiong Chen, Chengping Hu, and 3 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4009724/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 20 Apr, 2024 Read the published version in Respiratory Research → Version 1 posted 7 You are reading this latest preprint version Abstract Asthma-chronic obstructive pulmonary disease (COPD) overlap (ACO) represents a complex condition characterized by shared clinical and pathophysiological features of asthma and COPD in older individuals. However, the pathophysiology of ACO remains unexplored. We aimed to identify the major inflammatory cells in ACO, examine senescence within these cells, and elucidate the genes responsible for regulating senescence. Bioinformatic analyses were performed to investigate major cell types and cellular senescence signatures in a public single-cell RNA sequencing (scRNA-Seq) dataset derived from the lung tissues of patients with ACO. Similar analyses were carried out in an independent cohort study Immune Mechanisms Severe Asthma (IMSA), which included bulk RNA-Seq and CyTOF data from bronchoalveolar lavage fluid (BALF) samples. The analysis of the scRNA-Seq data revealed that monocytes/ macrophages were the predominant cell type in the lung tissues of ACO patients, constituting more than 50% of the cells analyzed. Lung monocytes/macrophages from patients with ACO exhibited a lower prevalence of senescence as defined by lower enrichment scores of SenMayo and expression levels of cellular senescence markers. Intriguingly, analysis of the IMSA dataset showed similar results in patients with severe asthma. They also exhibited a lower prevalence of senescence, particularly in airway CD206 + macrophages, along with increased cytokine expression (e.g., IL-4, IL-13 , and IL-22 ). Further exploration identified alveolar macrophages as a major subtype of monocytes/macrophages driving cellular senescence in ACO. Differentially expressed genes related to oxidation-reduction, cytokines, and growth factors were implicated in regulating senescence in alveolar macrophages. PPARγ (Peroxisome Proliferator-Activated Receptor Gamma) emerged as one of the predominant regulators modulating the senescent signature of alveolar macrophages in ACO. Collectively, the findings suggest that senescence in macrophages, particularly alveolar macrophages, plays a crucial role in the pathophysiology of ACO. Furthermore, PPARγ may represent a potential therapeutic target for interventions aimed at modulating senescence-associated processes in ACO. ACO Asthma COPD Macrophages Senescence PPARγ Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 INTRODUCTON Asthma and chronic obstructive pulmonary disease (COPD) are two of the most prevalent obstructive lung diseases. Asthma-COPD overlap (ACO) is a term used to describe a complex condition characterized by shared clinical and pathophysiological features of both diseases, typically observed in older individuals with a long history of disease and exposure to environmental factors 1 – 3 . ACO is estimated to affect approximately 10–30% of individuals with asthma and around 25% of those with COPD 4 . Cumulative environmental exposures, such as exposure to allergens and toxic particles or gases (e.g., smoking and indoor biomass), are believed to play a crucial role in the development and progression of ACO 1 , 5 . Additionally, asthma and atopy have been identified as potential risk factors for the development of COPD and, consequently, ACO 6 . More importantly, patients with ACO often experience increased rates of exacerbations and clinical symptoms, posing a significant challenge for clinical treatment and management 7 , 8 . However, managing ACO patients can be intricate due to the overlapping nature of these diseases. The inflammatory mechanisms in asthma and COPD are distinct. In asthma, inflammation is primarily driven by Th2 cell-mediated responses, leading to eosinophilic inflammation and an increase in cytokines like IL-4, IL-5, and IL-13 9 . On the other hand, COPD is characterized by inflammation with a higher presence of neutrophils and a dominance of Th1 and Th17 cell-mediated responses 10 . In both conditions, macrophages, the primary immune cells in the lungs, play crucial roles in the immune response, defense against infections, tissue homeostasis, and inflammation resolution 11 . Macrophages in asthma and COPD can contribute to airway inflammation and remodeling by releasing pro-inflammatory cytokines, proteases, and reactive oxygen species 12 , 13 . Therefore, it is reasonable to believe that macrophages also play a vital role in the inflammatory processes and pathophysiology of ACO. ACO is more commonly diagnosed in older individuals, and its prevalence tends to increase with age, as supported by various studies 14 , 15 . Age has been identified as a significant risk factor for ACO in asthma patients 16 . Consequently, the aging process and senescence appear to play substantial roles in the development of lung inflammation in individuals with ACO. Senescence represents a complex cellular state characterized by cellular stress, DNA damage, cell cycle arrest, and the release of senescence-associated secretory phenotype (SASP) factors 17 – 19 . These SASP factors encompass chemokines, cytokines, growth factors, adhesion molecules, and lipid components that can contribute to multiple age-related disorders with both local and systemic consequences 20 – 27 . It's worth noting that cellular senescence has been associated with both asthma 28 , 29 and COPD 30 – 32 . Consequently, senescence may exert a significant influence on the development and management of ACO. However, it is essential to recognize that senescence is a complex process influenced by factors such as cell type, age, and specific diseases 29 , 33 . Therefore, it is crucial to investigate whether senescence in lung macrophages contributes to the inflammatory processes and pathophysiology unique to ACO. In our current research, we conducted an analysis using publicly available single-cell RNA sequencing (scRNA-Seq) data obtained from human lung tissues, comparing individuals with and without ACO. Our primary focus was on identifying the major cell types and identified monocytes/macrophages as the predominant cell types in patients with ACO. We delved into gene signatures associated with senescence, specifically within monocytes/macrophages. The findings were further validated through another independent cohort (IMSA) with a specific emphasis on investigating the relationship between cellular senescence features and the severity of asthma. Additionally, we explored differentially expressed genes and pathways within alveolar macrophages of individuals with ACO and identified PPARγ as a key regulatory factor in driving cellular senescence in alveolar macrophages. METHODS Data Source The scRNA-Seq data was publicly available and was generated from human lung tissues of a patient with ACO who died with exacerbation and two transplant donors 34 . The Bulk RNA-Seq and mass cytometry (cytometry by time of flight, CyTOF) data were derived from bronchoalveolar lavage fluid (BALF) samples and downloaded from the Gene Expression Omnibus (GEO) database with the accession number GSE136587, and the Flow Repository (FR) database with the identifier FR-FCM-Z395 35 . The data were collected as a part of the Immune Mechanisms Severe Asthma (IMSA) study and involved the analysis of bronchoscopically obtained distal airway and alveolar cells. The study included a total of 39 subjects consisting of 6 healthy controls, 17 mild/moderate asthma patients, and 16 severe asthma patients. Single Cell RNA-Seq Processing and Analyzing Quality Control The raw counts matrixes downloaded were loaded to R workspace to initiate a Seurat workflow. A series of critical preprocessing steps were performed to ensure data quality before the analysis. Specifically, cells with a count of unique molecular identifiers (UMIs) below 500, cells expressing fewer than 200 genes with a log10GenesPerUMI value below 0.75, and cells with mitochondrial gene expression contributing to over 20% of the total gene expression were removed. Dimension Reduction and Cell Annotation Following the quality control procedures, the processed data underwent integration using distinct strategies: the complete dataset was integrated utilizing Seurat 36 , while the Monocytes/Macrophages subset was integrated using the Harmony approach 37 . The dimensionality reduction and visualization were achieved using the Uniform Manifold Approximation and Projection (UMAP) algorithm. UMAP not only captures inter-cellular distances but also provides a holistic view of the data's global structure 38 . Cell clusters were identified by Leiden algorithm, an automated algorithm tailored for effectively clustering cells in scRNA-seq data 39 . Cell markers were these with differentially expressed genes (DEGs) identified by the FindAllMarkers function from the Seurat package. Cellular annotations were identified by SingleR or manually based on specific marker expression, which allows us to get the composition and dynamics of cells within the dataset 40 . Differential Expressing Analyses Differential expressing analyses were conducted using the Libra R package (Version 1.0.0). Within this framework, the "run_de" function was employed, applying Wilcoxon rank sum tests to assess statistical differences between groups. Genes with an adjusted P value 0.5 were termed as statistical significance. Moreover, the Seurat package's "FindMarkers" and "FindAllMarkers" functions were also utilized to identify markers that distinguish cell clusters with Wilcoxon rank sum tests. Enrichment Analyses The enrichment analyses included geneset scoring through the AUCell algorithm, and over-representation analysis (ORA). The AUCell algorithm was employed to gauge cellular senescence enrichment. This algorithm calculated an enrichment score by juxtaposing the input data against the SenMayo senescence geneset, an optimal choice for senescence screening since it consistently aligns with results from biological validation experiments 41 . The ORA analysis function RunEnricher was embedded within the High-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA) R package. It was employed to investigate the biological characteristics of selected co-expression modules. Trajectory Analyses A partition-based graph abstraction (PAGA) based trajectory analysis was used to investigate the trajectory of monocytes/macrophages within single-cell RNA sequencing (scRNA-seq) data. This procedure was performed using Scanpy version 1.9.1 in Python 3.10.9., which revealed the intricate relationships and transitions among different cell states. This workflow simplified the dataset's graph representation by clustering similar cells, enabling the identification of connected and disconnected regions. By combining well-defined paths rooted in high-confidence connections and leveraging a random-walk-based distance metric, cells are systematically ordered within each partition according to their proximity to a designated reference cell 42 . High-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA) To explore the hub genes and the underlying expression network associated with cellular senescence and related biological processes, we harnessed the power of hdWGCNA (Version 0.2.18). This algorithm is a refinement of the traditional WGCNA designed specifically for scRNA-Seq data. It operates by computing a sparsity-reduced expression matrix termed metacell, akin to the concept of pseudobulk to mitigate the influence of undetected genes and facilitate the identification of gene co-expression modules with clinical or biological information 43 , 44 Transcription Factor Analyses Transcription factors are crucial for the initial stage of decoding DNA sequencing and regulation of gene expression 45 . To initiate this process, we employed PySCENIC (Version 0.12.1) 46 to identify cell cluster-specific transcription factors. Moreover, we extended our investigation by extracting a gene list from the cellular senescence-related gene module identified using hdWGCNA. Subsequently, this gene list was subjected to transcription factor enrichment analysis using the Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining (TRRUST) database on the metascape website (metascape.org) 47 , 48 . Bulk RNA-Seq Processing and Analyzing Data Preprocessing The expression matrix encompassing raw counts and pertinent phenotype information was procured from the GEO website, facilitated by the GEOQuery R package (Version 2.66.0) 49 . The raw counts with batch details were downloaded for differential expressing analysis. Furthermore, the raw counts were normalized and then transformed into log2-transformed transcripts per million (log2TPM) values. To mitigate the batch-induced discrepancies, the Combat algorithm embedded within the sva R package (Version 3.44.0) 50 . Unsupervised Clustering and Differential Expressing Analyses Due to heterogeneity of asthma, we did not make direct comparisons for senescence signature between different asthma groups. Instead, we employed an unsupervised clustering approach that focused on the expression profiles of genes associated with cellular senescence 51 . For the clustering process, we utilized the k-means algorithm, which groups data points based on their similarity as measured by the Euclidean distance with the ConsensusClusterPlus R package (Version 1.60.0) 52 , 53 . The differential expressing analysis was performed to explore the variations in gene expression across these clusters by following the limma-voom approach with batch information adjusted in the model 54 . CyTOF Processing and Analyzing Data Preprocessing An outlier sample displaying an unusually high intensity value of 3 in the live_dead channels was identified and subsequently excluded from the CyTOF analysis. Cells exhibiting extremely low CD45 expression were also filtered out using Flowjo. The resultant fcs data was then transferred to an R workspace for subsequent analysis. To ensure accurate analysis, the data underwent a bead-based normalization process utilizing the "normCytof" function within the CATALYST R package (Version 1.20.1) 55 . Dimension Reduction and Cell Annotation Similar to scRNA-Seq, dimension reduction for visualization and cell annotation were required for CyTOF core analyses. Especially, the t-Distributed Stochastic Neighbor Embedding (t-SNE) algorithm was employed 56 . For automated cluster identification, the FlowSOM 57 was used. Cell type annotations were based on expression profiles of major cell type markers similar to that described in relevant publications 35 . Differential Abundance (DA) and Differential States (DS) Analyses Based on CyTOF Data Differential cell proportion (abundance) and cytokine expression (states) between cellular senescence clusters identified by unsupervised clustering were calculated with diffcyt R package (1.16.0). Both limma based DA (diffcyt-DS-voom), and DS (diffcyt-DA-trend) parameters were chosen for differential expression analyses 54 , 58 . Statistical Analysis and Data Visualization All analysis and visualizations were processed with R (Version 4.1.2 for scRNA-Seq and Version 4.2.1 for CyTOF and Bulk RNA-Seq), Python (Version 3.10.9) and GraphPad (Version 9.3). Wilcoxon rank sum tests were employed for continuous variables between two groups; Chi-Square tests were performed for the comparisons of categorical variables. Correlation analyses (Pearson) were based on Single-Cell Variational Inference (scVI-tools, Version 0.20.2) imputed expression data to avoid sparsity influence 59 . Data visualizations were performed with ggplot2 (Version 3.4.1), genekitr (Version 1.0.5), ggVennDiagram (Version 1.2.2) and ComplexHeatmap (Version 2.10.0) R packages. RESULTS Alveolar Macrophages were the Major Cell types in Patients with ACO To better understand the pathogenesis of ACO, we analyzed the existing scRNA-Seq dataset generated from human lung tissues of a patient with ACO and two transplant donors without ACO 34 . By using the Leiden algorithm, an automated algorithm designed to effectively cluster cells in scRNA-seq data 39 , we identified a total of 12 cell clusters within dataset (Fig. 1 , A). These cell clusters were further annotated as different cell types by using the SingleR algorithm with manual adjustments. Cell markers were these differentially expressed across different cell types identified by the FindAllMarkers function from the Seurat package 36 . A total of nine different cell types were identified within the dataset, which included monocytes/macrophages, T cells, NK cells, AT2 alveolar type II cells, endothelial cells, airway epithelial cells, B cells, fibroblasts, and mast cells (Fig. 1 , B). The top-ranking differentially expressed genes for each cell type were presented in Fig. 1 , C and Figure E1). For example, complement-related genes ( C1QA, C1QB ), CD68 , and APOC1 were highly expressed in monocytes/macrophages (Figure E1, A), while genes like SFTPB, SFTPA1 , and SFTPA2 were uniquely and significantly expressed in AT2 airway epithelial cells ( Figure E1, B). Tryptase genes ( TPSAB1, TPSB2 ) were highly expressed in mast cells ( Figure E1, C). To identify the distribution of different cell types in the samples from ACO and control group, we performed cell proportion analysis as defined by relative percentages (Fig. 1 , D). Of these, monocytes/macrophages were a predominant cell type, constituting more than 50% of the cells in the different groups being analyzed. Furthermore, the proportion of monocytes/macrophages was significantly higher in ACO patient compared to control group ( p < 0.001). Next, we focused on those monocytes/macrophages specifically and re-clustered them by using the Leiden algorithm. A total of seven clusters were identified (Fig. 1 , E). These cell clusters were further annotated as different sub-types of monocytes/macrophages (Fig. 1 , F), including alveolar macrophages (AM), cycling cells, interstitial macrophages (IM), and monocytes. The top-ranking differentially expressed genes for each cell type were presented in Fig. 1 , G and Figure E2). Similar to monocytes/macrophages, the subtype AMs express the complement-related genes C1QA, C1QB ( Figure E2, A). Additionally, several other genes RBP4, CD9, SERPING1 , and CES1 are also significantly expressed in AMs. The proliferation markers PCLAF, TOP2A and MKI67 were highly expressed in cycling cells ( Figure E2, B). IMs were in an intermediate state in the UMAP plot and expressed high levels of molecules such as LGMN , RNASE1 , and CCL2 ( Figure E2, C). While CFP, FCN1 , and S100A8 ,were highly expressed in monocytes ( Figure E2, D). Taken together, the results suggest that monocytes/macrophages in ACO are major cells that may drive airway inflammation. Decreased cellular Senescence in Monocytes/Macrophages of Patients with ACO Cellular senescence has gained considerable attention across various diseases, including respiratory conditions such as asthma 29 . To determine whether senescence contributes to the development of ACO, we specifically assessed the enrichment of genes associated with senescence within lung monocytes/macrophages. Employing the AUCell algorithm, we computed enrichment scores against a senescence-related gene set known as SenMayo 41 . Among all the subtypes, AMs/monocytes showed the higher enrichment sore as assessed by the density of AUCell enriched SenMayo senescence scores (Fig. 2 , A). Next, we investigated the distribution of ACO and control group among monocytes/macrophages (Fig. 2 , B) and enrichment scores of SenMayo senescence in ACO and control group (Fig. 2 , C). As illustrated in Fig. 2 , C, the senescent signatures are enriched in the control group as compared to ACO group. Among them, AMs were a predominant cell type, constituting more than 50% of the cells among all these subtypes (Fig. 2 , D). The proportion of AMs was significantly higher in ACO patients as compared to control group ( p < 0.001). Additionally, cell cycle arrested in the G1 phase is one of the important features of cellular senescence 60 . We found that monocytes/macrophages that were in the G1 phase of the cell cycle were more abundant in the control group compared with ACO group (Fig. 2 , E). These findings were further supported by the expression of cellular senescence markers, CDKN1A (p21) and CDKN2A (p16) 61 . The expression of CDKN1A was much lower in the ACO group as compared to the control group (Fig. 2 , F). In contrast, no significant change was observed for CDKN2A (Fig. 2 , G). These findings suggest a lower prevalence of cellular senescence within lung monocytes/macrophages among patients with ACO. A lower prevalence of senescence observed in patients with severe asthma We further validated the results through another independent cohort (IMSA) 35 and analyzed the relationship between cellular senescence features and the severity of asthma. After stringent quality control, a total of 39 subjects were finally included for analysis: 6 healthy controls (HC), 17 mild/moderate asthma patients (MMA), and 16 severe asthma patients (SA). Definition of these subjects has been previously reported 9 . Consensus clustering analysis with the Consensus Cumulative Distribution Function (CDF) identified cellular senescence clusters for immune cells from IMSA. Consensus clustering analysis involves repeated clustering of the data with different values of k (from 2 to 5) within a dataset to identify stable and meaningful clusters (Fig. 3 , A). Based on the consensus CDF (Fig. 3 , B), it is evident that cellular senescence clusters 2 provides the most consistent and stable clustering solution for the IMSA dataset. According to the ssGSEA SenMayo Senescence score, the cluster 2 was divided into senescence clustering low and high group ( P = 0.002, Fig. 3 , C). Senescence clustering low group showed lower expression of senescence markers CDKN1A as compared to senescence high group ( P = 0.011, Fig. 3 , D). No difference was noted for CDKN2A between senescence clustering low and high group ( P = 0.29, Fig. 3 , E). In contrast, senescence low group showed higher expression of cell proliferation marker MKI67 ( P = 0.0015, Fig. 3 , F). Intriguingly, among the analyzed subjects consisting of HC, MMA, and SA participants, the large proportion of severe asthmatic patients was observed in senescence low group as compared to senescence high group ( P < 0.05, Fig. 3 , G). These findings indicate that cellular senescence may be negatively associated with the severity of asthma. Airway CD206 + macrophages show a lower prevalence of senescence and increased cytokines To further explore the relationship between cellular senescence and airway inflammation in asthma, we focused on the CyTOF dataset that targets lineage markers for the adaptive and innate immune systems 35 . Unsupervised cell clustering was performed in BAL fluids from HC, MMA, and SA participants using the FlowSOM algorithm 62 , 63 . Cell types in those clusters were determined by the surface marker staining intensities across t-SNE spaces. A total of 7 different cell clusters were annotated using a combination of surface marker genes, including B lymphocytes (CD19 + ), CD206 − macrophages (CD11C + CD206 − ), CD206 + macrophages (CD11C + CD206 − ), CD4 T lymphocytes (CD3 + CD4 + ), CD8 T lymphocytes (CD3 + CD8 + ), γδ T lymphocytes (TCRγδ + ), and NK cells (CD56 + ) (Fig. 4 , A, B, and Figure E3). To identify the link between these airway immune cells and cellular senescence, the percentages of those identified cell populations were compared in senescence high and low groups. Among all these populations, a significant difference was found for CD206 + macrophages, which showed a higher proportion in senescence low in relative to senescence high group (Fig. 4 , C). Interestingly, the senescence low group showed elevated levels of cytokines IL-5, IL-13, IL-10, and IL-22 in BAL fluids of the IMSA cohort (Fig. 4 , D). Particularly, CD206 + macrophages showed increased expression of IL-4, IL-13, and IL-22 in senescence low group as compared those in senescence high group (Fig. 4 , E). Taken together, these data suggest that CD206 + macrophages are major airway immune cells associated with cellular senescence and cellular senescence is negatively correlated with cytokines in airway CD206 + Macrophages. Differentially expressed genes and pathways in Alveolar macrophages of ACO Given the significance of CD206 + macrophages in cellular senescence, we re-visited the ACO dataset and determined whether the specific subtype of monocytes/macrophages drive the difference in cellular senescence between ACO and control group. To explore this, we used trajectory analysis to infer the relationships and distances between different subtypes of monocytes/macrophages. Analyses were started with cycling cells, we identified a specific cluster that has the highest distance (pseudo-time) and is in the upper right corner of a two-dimensional scatter plot (Fig. 5 , A). This cluster was mostly composed of cells from Leiden cluster 1 (Fig. 5 , B) and was predominantly found in the control group (Fig. 5 , C), suggesting that cells in Leiden cluster 1 may contribute to the difference in biological processes between ACO and control group. To further characterize the Leiden cluster 1, we employed the FindAllMarkers to investigate the highly expressed marker genes for cluster 1 against other cell clusters. A total of 47 genes were identified specifically for the Leiden cluster 1, such as SERPING1, CES1, HLA-DQA1, INHBA, FABP4, LGALS3BP, AKR1C3, CITED2, TERM1, FN1, EVL, ALDH1A1 , and PDLIM1 (Fig. 5 , D and see Table E2). Additionally, we checked the lineage markers that are commonly used for monocytes/macrophages (Fig. 5 , E). we found that PPARγ, FCGR1A, FCGR3A, MSR1 (CD205), and MRC1 (CD206) were enriched in Leiden cluster 1. Furthermore, we used AUCell scoring to analyze the biological characteristics of Leiden cluster 1 in comparison to other clusters. Especially, the KEGG enrichment analysis (Fig. 5 , F and see Table E3) indicated that Leiden cluster 1 had significant enrichments in several pathways, such as antigen processing and presentation, Peroxisome Proliferator-Activated Receptor (PPAR) signaling pathway. Next, we investigated the differences in cellular senescence between ACO and control groups specially in Leiden cluster 1. Consistent with our previous analyses among all the monocytes/macrophages, the SenMayo senescence score was remarkably lower in ACO groups compared with control group (Fig. 5 , G). This was further supported by the expression of cellular senescence marker CDKN1A (Fig. 5 , H), which also showed reduced expression in the ACO group. No statistical difference was observed for CDKN2A (Fig. 5 , I). Additionally, to further display the characteristics of cellular senescence, we identified differentially expressed genes related to oxidation-reduction, cytokines, and growth factors between ACO and control groups in Leiden cluster 1. Of those, TGFB1 was increased but several chemokines such as CXCL2, CXCL3, CXCL8, and CCL18 were decreased in ACO group compared with control group (Fig. 5 , J and see Table E4). These results suggest correlations between cellular senescence and cytokine expression in alveolar macrophages of Patients with ACO. Furthermore, strong correlations were identified for CDKN1A expression levels and those identified genes in Leiden cluster 1 and traditional lineage markers for monocytes/macrophages (Fig. 5 , K and see Table E5 in this article’s Online Repository at www.jaci.org ). For example, CDNK1A expression was correlated with genes with either positive regulation (e.g., CD14, CD68, CD86, CD163, FCGR3A, MSR1, SERPING1, HLA-DQA1, INHBA, LGALS3BP, AKR1C, TREM1, FN1 , and ALDH1A1 ) or negative regulation (e.g., MARCO, FCGR1A, ITGB2 , and MRC1 ). Collectively, the results indicate that these identified genes and biological pathways may be involved in regulating cellular senescence in alveolar macrophages of patients with ACO. PPARγ, a Key Regulating Factor of Cellular Senescence in Alveolar Macrophages To delve deeper into the underlying regulatory mechanisms of biological changes, we first utilized the hdWGCNA analysis focused on Leiden cluster 1 to explore the expression network associated with senescence signatures. The algorithm discerned six gene modules: green, blue, red, brown, cyan, and tan (Fig. 6 , A). Within these, CDKN1A (p21) was pinpointed in the brown module, and the module feature gene score was positively correlated with the difference in cellular senescence between ACO and control groups (R = 0.669, p < 0.0001, Fig. 6 , B). Prominently, the top 30 hub genes with high correlation in this module encompassed those related to growth arrest and proliferation (e.g., GADD45B, PPP1R15A , and DUSP2 ), TP53-mediated cell senescence-associated heat shock proteins ( DNAJA1, DNAJB1 , and HSPA5 ), DNA damage and repair (e.g., MYL12A), and energy, protein, and lipid metabolism (e.g., GLUL, PLIN2, RPS4Y1, PDIA3 , and SOD2 ). Furthermore, genes associated with SASP and its upstream regulators (e.g., CCL20, CXCL3, CXCL5, CXCL8, JUN, NFKB1A , and NFKBIZ ) and those tied to antigen presentation, antibody response, and macrophage activation (e.g., C83, CD68, FCGR3A , and B2M ) were also included in these hub genes (Fig. 6 , C). These genes or biological processes are directly or indirectly involved in senescence processes of alveolar macrophages. Subsequent enrichment analyses of entire module genes showed that the brown module was mainly associated with macrophage markers and cellular senescence (Fig. 6 , D), cytokine signaling in immune system (Fig. 6 , E), and antigen processing and presentation and PPAR signaling (Fig. 6 , F). These biological characteristics are closely related and mutually regulated. To pinpoint the central regulator of senescence in alveolar macrophages, we leveraged SCENIC to discern the most specific transcription factors for each Leiden cluster. Notably, for Leiden cluster 1, PPARγ emerged as the predominant regulator (Fig. 6 , G). This was further supported by the transcription factor enrichment analysis with the brown module against the TRRUST database. PPARγ was identified to be one of the top 10 pivotal transcription factors modulating gene expression within this module (Fig. 6 , H). Intriguingly, each of these transcription factors had the potential to regulate CDKN1A expression. Taken together, these findings indicate that the brown module had the core impact on cellular senescence, and PPARγ is one of the predominant regulators modulating the senescent signature of alveolar macrophages. DISCUSSION Both asthma and COPD can co-exist in older individuals with clinically overlapping phenotypes 1 – 3 . ACO has recently been a focus of interest because patients with ACO often experience poor health-related quality of life, increased rates of exacerbations, and severe clinical symptoms 7 , 8 . While the exact pathophysiology of ACO is not fully understood, it is thought to result from various features of both asthma and COPD 64 , 65 . Both asthma and COPD involve chronic airway inflammation, but they have different underlying mechanisms. In asthma, there is typically eosinophilic inflammation driven by a Th2 immune response 65 , while COPD is characterized by neutrophilic inflammation associated with Th1 and Th17-driven inflammation 10 , 66 , 67 . In both conditions, macrophages play a central role in the recruitment and activation of neutrophils and eosinophils 11 – 13 , 68 . Indeed, a very recent cross-sectional exploratory study was performed to investigate cellular changes (e.g., macrophages, neutrophils, eosinophils, mast cells, CD8 + , and CD4 + T lymphocytes) in the airway wall of ACO compared with asthma, COPD current smokers, and ex-smokers, normal lung function smokers, and non-smoker controls. They found that the ACO airway tissue inflammatory cellular profile differed from the contributing diseases of asthma and COPD with a predominance of macrophages 65 . The results were supported by our analyses on an existing scRNA-Seq dataset generated from human lung tissues of patients with ACO. Monocytes/macrophages were a predominant cell type among all cell type to be analyzed (monocytes/macrophages, T cells, NK cells, AT2 alveolar type II cells, endothelial cells, airway epithelial cells, B cells, fibroblasts, and mast cells), constituting more than 50% of the total cells, and the proportion of monocytes/macrophages was significantly higher in patients with ACO in relative to control group. Importantly, our further study identified several sub-types of monocytes/macrophages, including alveolar macrophages, cycling cells, interstitial macrophages, and monocytes. Of these, alveolar macrophage was the most predominant cell type with the increased expression of complement-related genes C1QA, C1QB, and several other genes RBP4, CD9, SERPING1, and CES1. Thus, our findings provide further evidence that monocytes/macrophages, particularly alveolar macrophages, are major cells that may contribute to airway inflammation in patients with ACO. Next, we explored how monocytes/macrophages drive airway inflammation with the emphasis on cellular senescence that has been implicated in the pathophysiology of various diseases, including asthma 29 . Cellular senescence is an extremely complex and dynamic biological process induced by several factors, including aging, oxidative stress, DNA damage, mitochondrial dysfunction, epigenetic modifications, and telomere shortening 29 , 69 , 70 . Cellular senescence has been associated with both asthma 29 and COPD 30 – 32 . Several stimuli causing cellular senescence have been associated with asthma, such as telomere shortening, oxidative stress, inflammation, and autophagy/mitophagy. Therefore, senescence may exert a significant influence on the function and activation of ACO-associated target cells and subsequently development and management of ACO. Here we investigated the relationship between senescence and ACO by analyzing gene signatures associated with senescence in all the subtypes of monocytes/macrophages. We demonstrated for the first time a lower prevalence of cellular senescence in alveolar macrophages of patients with ACO. Furthermore, CDKN1A (p21) but not CDKN2A (p16) showed lower expression in alveolar macrophages of patients with ACO, suggesting that CDKN1A (p21) pathway may play a major role in mediating senescence of alveolar macrophages in ACO. To strengthen our findings on the role of cellular senescence in ACO, we analyzed the relationship of cellular senescence features with the severity of asthma in a total of 39 individuals from IMSA 35 . Specifically, senescence in those individuals were clustered into senescence low and high group and then the distribution of healthy controls, mild/moderate asthma, and severe asthma was analyzed in these two groups. Intriguingly, the large proportion of severe asthmatic patients was observed in senescence clustering low group compared with senescence clustering high group. Furthermore, we focused on the CyTOF dataset from BAL fluids that targets lineage markers and identified a total of 7 different cell clusters, including B lymphocytes, CD206 − macrophages, CD206 + macrophages, CD4 + T, CD8 + T, γδ T, and NK cells. Of these, CD206 + macrophages were distributed in senescence clustering low group and showed increased expression of Th2 cytokines IL-4, IL-13 and Th17 cytokine IL-22, suggesting that, CD206 + macrophages may drive both Th2 and Th17 populations commonly seen in patients with ACO 64 , 65 . CD206 expression is a characteristic feature of some tissue-resident macrophages, including alveolar macrophages 11 , 71 , and CD206 + macrophages represent a large proportion of alveolar macrophages. Collectively, these independent data provide evidence that patients with severe asthma have a lower prevalence of cellular senescence, and CD206 + macrophages may contribute to the severity of asthma, possibly ACO, by releasing inflammatory cytokines. We also explored the underly mechanisms as for how cellular senescence is regulated in macrophages. We investigated the differentially expressed genes and pathways in alveolar macrophages in Leiden cluster 1, a population mostly likely contributing to the difference in cellular senescence between ACO and control group. Several pathways were also identified, such as antigen processing and presentation, PPARγ signaling pathway. Additionally, several markers commonly used for monocytes/macrophages also showed different expression, including PPARγ, FCGR1A, FCGR3A, CD205, and CD206. As expected, most of these differentially expressed genes were either positively or negatively associated with senescence as defined by the expression of CDNK1A and enrichment scores of SenMayo senescence (data not shown). To further delve deeper into the underlying regulatory mechanisms of biological changes, we utilized the hdWGCNA analysis focused on alveolar macrophages and identified several expression networks associated with senescence signatures, including growth arrest and proliferation, TP53-mediated cell senescence-associated heat shock proteins, DNA damage and repair, and energy, protein, and lipid metabolism. These genes or biological processes are directly or indirectly involved in senescence processes of alveolar macrophages. Importantly, we identified PPARγ as one of the top 10 pivotal transcription factors modulating senescent signatures of alveolar macrophages. Interestingly, PPARγ as a transcription factor has been shown to play a crucial role in various physiological processes, including adipogenesis 72 , 73 , glucose homeostasis 74 , inflammation 75 , 76 , and metabolism and function of macrophages 77 . Studies have implicated that PPARγ regulates senescence in various cell types, including fibroblasts and endothelial cells 78 , 79 , and macrophages 77 . PPARγ activation can promote senescence by upregulating the expression of specific senescence-associated genes and inducing cell cycle arrest 79 . PPARγ activation can also have anti-senescent effects by reducing senescence and improve metabolic function 79 – 82 Taken together, we present groundbreaking findings, demonstrating for the first time that monocytes/macrophages, particularly alveolar macrophages, are the predominant cell types in patients with ACO. Our research reveals a lower prevalence of senescence within alveolar macrophages in both ACO patients and those with severe asthma. Mechanistically, our in-depth exploration of differentially expressed genes within alveolar macrophages has identified PPARγ as a crucial regulatory factor driving senescence in ACO. Notably, there are several limitations, including small sample size, lack of mouse model to mimic patients with ACO, and complicated and context-dependent role of PPARγ in cellular senescence. However, our observations in this study provide a solid foundation for future investigations aimed at a comprehensive understanding of the mechanisms through which PPARγ regulates senescence. Furthermore, it opens doors for exploring PPARγ as a potential therapeutic target for interventions aimed at modulating senescence-associated processes in ACO -a condition marked by mixed features of both asthma and COPD, known for its increased severity. Abbreviations IL Interleukins AM Alveolar Macrophages IM Interstitial Macrophages SASP Senescence Associated Secretory Phenotype COPD Chronic Obstructive Pulmonary Diseases ACO Asthma-COPD Overlap CyTOF Cytometry by Time of Flight BALF Bronchoalveolar Lavage Fluid IMSA Immune Mechanisms Severe Asthma UMAP Uniform Manifold Approximation and Projection DEGs Differential Expressed Genes H&E Hematoxylin and Eosin SCENIC Single-Cell rEgulatory Network Inference and Clustering hdWGCNA High-Dimensional Weighted Gene Co-epression Network Analysis PAGA Partition-Based Graph Abstraction GEO Gene Expression Omnibus GSEA Geneset Enrichment Analysis TRRUST: Transcriptional Regulatory Relationships Unraveled by Sentence-based Text ` mining ssGSEA Single Sample GSEA CRE Cockroach Extract KEGG Kyoto Encyclopedia of Genes and Genomes GO Gene Ontology PPARγ Peroxisome Proliferator Activated Receptor Gamma Declarations Disclosure of potential conflict of interest: The authors have declared that no conflict of interest exists Author Contribution analyzed data, performed experiments, and review the manuscript. R.W., M. W., P. G. designed and wrote the manuscript. Q. C., C. H., M. W., P. G., supervised the study. R.W., P. G. designed and wrote the manuscript. All authors read and approved the final version of the manuscript. ACKNOWLEDGEMENTS This work was supported by grants from the US National Institutes of Health (NIH) (1R01AI153331 and R01AI141642 to Gao). Availability of Data All bioinformatical data can be downloaded from figshare database (scRNA-Seq from fatal ACO and two health control: https://doi.org/10.6084/m9.figshare . 14782377.v1.), GEO (Bulk RNA-Seq data from IMSA BALF: GSE136587) database and flow repository database (CyTOF data from IMSA BALF: FR-FCM-Z395). References Leung JM, Sin DD. Asthma-COPD overlap syndrome: pathogenesis, clinical features, and therapeutic targets. BMJ. 2017;358:j3772. Mekov E, Nuñez A, Sin DD, et al. Update on Asthma-COPD Overlap (ACO): A Narrative Review. Int J Chron Obstruct Pulmon Dis. 2021;16:1783–99. Barczyk A, Maskey-Warzęchowska M, Górska K, et al. Asthma-COPD Overlap-A Discordance Between Patient Populations Defined by Different Diagnostic Criteria. J Allergy Clin Immunol Pract. 2019;7(7):2326–36e5. Alshabanat A, Zafari Z, Albanyan O, Dairi M, FitzGerald JM. Asthma and COPD Overlap Syndrome (ACOS): A Systematic Review and Meta Analysis. PLoS ONE. 2015;10(9):e0136065. Hikichi M, Hashimoto S, Gon Y. Asthma and COPD overlap pathophysiology of ACO. Allergol Int. 2018;67(2):179–86. Gottlieb DJ, Sparrow D, O'Connor GT, Weiss ST. Skin test reactivity to common aeroallergens and decline of lung function. The Normative Aging Study. Am J Respir Crit Care Med. 1996;153(2):561–6. Caillaud D, Chanez P, Escamilla R et al. Asthma-COPD overlap syndrome (ACOS) vs 'pure' COPD: a distinct phenotype? Allergy 2017; 72(1): 137 – 45. Sánchez Castillo S, Smith L, Díaz Suárez A, López Sánchez GF. Limitations in Activities of Daily Living among Older Adults with COPD, Asthma, or Asthma-COPD Overlap Residing in Spain. Int J Environ Res Public Health 2023; 20(4). Woodruff PG, Modrek B, Choy DF, et al. T-helper type 2-driven inflammation defines major subphenotypes of asthma. Am J Respir Crit Care Med. 2009;180(5):388–95. Qin K, Xu B, Pang M, Wang H, Yu B. The functions of CD4 T-helper lymphocytes in chronic obstructive pulmonary disease. Acta Biochim Biophys Sin (Shanghai). 2022;54(2):173–8. Saradna A, Do DC, Kumar S, Fu QL, Gao P. Macrophage polarization and allergic asthma. Transl Res. 2018;191:1–14. Cai Y, Sugimoto C, Arainga M, Alvarez X, Didier ES, Kuroda MJ. In vivo characterization of alveolar and interstitial lung macrophages in rhesus macaques: implications for understanding lung disease in humans. J Immunol. 2014;192(6):2821–9. Draijer C, Peters-Golden M. Alveolar Macrophages in Allergic Asthma: the Forgotten Cell Awakes. Curr Allergy Asthma Rep. 2017;17(2):12. Baptist AP, Busse PJ. Asthma Over the Age of 65: All's Well That Ends Well. J Allergy Clin Immunol Pract. 2018;6(3):764–73. de Marco R, Pesce G, Marcon A, et al. The coexistence of asthma and chronic obstructive pulmonary disease (COPD): prevalence and risk factors in young, middle-aged and elderly people from the general population. PLoS ONE. 2013;8(5):e62985. Mendy A, Forno E, Niyonsenga T, Carnahan R, Gasana J. Prevalence and features of asthma-COPD overlap in the United States 2007–2012. Clin Respir J. 2018;12(8):2369–77. Di Micco R, Krizhanovsky V, Baker D. d'Adda di Fagagna F. Cellular senescence in ageing: from mechanisms to therapeutic opportunities. Nat Rev Mol Cell Biol. 2021;22(2):75–95. Santoro A, Bientinesi E, Monti D. Immunosenescence and inflammaging in the aging process: age-related diseases or longevity? Ageing Res Rev. 2021;71:101422. Shvedova M, Samdavid Thanapaul RJR, Thompson EL, Niedernhofer LJ, Roh DS. Cellular Senescence in Aging, Tissue Repair, and Regeneration. Plast Reconstr Surg. 2022;150:4s–11s. Muñoz-Espín D, Cañamero M, Maraver A, et al. Programmed cell senescence during mammalian embryonic development. Cell. 2013;155(5):1104–18. Wang X, Qu M, Li J, Danielson P, Yang L, Zhou Q. Induction of Fibroblast Senescence During Mouse Corneal Wound Healing. Investig Ophthalmol Vis Sci. 2019;60(10):3669–79. Calcinotto A, Kohli J, Zagato E, Pellegrini L, Demaria M, Alimonti A. Cellular Senescence: Aging, Cancer, and Injury. Physiol Rev. 2019;99(2):1047–78. Vicente R, Mausset-Bonnefont AL, Jorgensen C, Louis-Plence P, Brondello JM. Cellular senescence impact on immune cell fate and function. Aging Cell. 2016;15(3):400–6. Kang TW, Yevsa T, Woller N, et al. Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. Nature. 2011;479(7374):547–51. Lecot P, Alimirah F, Desprez PY, Campisi J, Wiley C. Context-dependent effects of cellular senescence in cancer development. Br J Cancer. 2016;114(11):1180–4. Yamamoto-Imoto H, Minami S, Shioda T, et al. Age-associated decline of MondoA drives cellular senescence through impaired autophagy and mitochondrial homeostasis. Cell Rep. 2022;38(9):110444. Moiseeva V, Cisneros A, Sica V, et al. Senescence atlas reveals an aged-like inflamed niche that blunts muscle regeneration. Nature. 2023;613(7942):169–78. Soma T, Nagata M. Immunosenescence, Inflammaging, and Lung Senescence in Asthma in the Elderly. Biomolecules 2022; 12(10). Wan R, Srikaram P, Guntupalli V, Hu C, Chen Q, Gao P. Cellular senescence in asthma: from pathogenesis to therapeutic challenges. EBioMedicine. 2023;94:104717. Woldhuis RR, Heijink IH, van den Berge M, et al. COPD-derived fibroblasts secrete higher levels of senescence-associated secretory phenotype proteins. Thorax. 2021;76(5):508–11. Woldhuis RR, de Vries M, Timens W, et al. Link between increased cellular senescence and extracellular matrix changes in COPD. Am J Physiol Lung Cell Mol Physiol. 2020;319(1):L48–L60. Adams TS, Schupp JC, Poli S et al. Single-cell RNA-seq reveals ectopic and aberrant lung-resident cell populations in idiopathic pulmonary fibrosis. Sci Adv 2020; 6(28): eaba1983. Schuliga M, Read J, Knight DA. Ageing mechanisms that contribute to tissue remodeling in lung disease. Ageing Res Rev. 2021;70:101405. Moghbeli K, Valenzi E, Naramore R, et al. β-Agonist exposure preferentially impacts lung macrophage cyclic AMP-related gene expression in asthma and asthma COPD overlap syndrome. Am J Physiol Lung Cell Mol Physiol. 2021;321(5):L837–l43. Camiolo MJ, Zhou X, Oriss TB, et al. High-dimensional profiling clusters asthma severity by lymphoid and non-lymphoid status. Cell Rep. 2021;35(2):108974. Hao Y, Hao S, Andersen-Nissen E, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184(13):3573–87e29. Korsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods. 2019;16(12):1289–96. McInnes L, Healy J, Melville J, Umap. Uniform manifold approximation and projection for dimension reduction. arXiv preprint arXiv:180203426 2018. Traag VA, Waltman L, van Eck NJ. From Louvain to Leiden: guaranteeing well-connected communities. Sci Rep. 2019;9(1):5233. Aran D, Looney AP, Liu L, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163–72. Saul D, Kosinsky RL, Atkinson EJ, et al. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat Commun. 2022;13(1):4827. Wolf FA, Hamey FK, Plass M, et al. PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells. Genome Biol. 2019;20(1):59. Morabito S, Reese F, Rahimzadeh N, Miyoshi E, Swarup V. High dimensional co-expression networks enable discovery of transcriptomic drivers in complex biological systems. bioRxiv 2022: 2022.09.22.509094. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559. Lambert SA, Jolma A, Campitelli LF, et al. Hum Transcription Factors Cell. 2018;172(4):650–65. Aibar S, Gonzalez-Blas CB, Moerman T, et al. SCENIC: single-cell regulatory network inference and clustering. Nat Methods. 2017;14(11):1083–6. Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523. Han H, Cho JW, Lee S, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380–d6. Davis S, Meltzer PS. GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics. 2007;23(14):1846–7. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3. Jassal B, Matthews L, Viteri G, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48(D1):D498–d503. Hartigan JA, Wong MA, Algorithm AS. A k-means clustering algorithm. J royal Stat Soc Ser c (applied statistics). 1979;136(1):100–8. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572–3. Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. Nowicka M, Krieg C, Crowell HL, et al. CyTOF workflow: differential discovery in high-throughput high-dimensional cytometry datasets. F1000Research. 2019;6:748. Van der Maaten L, Hinton G. Visualizing data using t-SNE. J Mach Learn Res 2008; 9(11). Quintelier K, Couckuyt A, Emmaneel A, Aerts J, Saeys Y, Van Gassen S. Analyzing high-dimensional cytometry data using FlowSOM. Nat Protoc. 2021;16(8):3775–801. Weber LM, Nowicka M, Soneson C, Robinson MD. diffcyt: Differential discovery in high-dimensional cytometry via high-resolution clustering. Commun Biol. 2019;2:183. Lopez R, Regier J, Cole MB, Jordan MI, Yosef N. Deep generative modeling for single-cell transcriptomics. Nat Methods. 2018;15(12):1053–8. Martinez-Zamudio RI, Robinson L, Roux PF, Bischof O, SnapShot. Cell Senescence Pathways Cell. 2017;170(4):816–e1. Chandra A, Lagnado AB, Farr JN, et al. Targeted clearance of p21- but not p16-positive senescent cells prevents radiation-induced osteoporosis and increased marrow adiposity. Aging Cell. 2022;21(5):e13602. Kohonen T. Essentials of the self-organizing map. Neural Netw. 2013;37:52–65. Van Gassen S, Callebaut B, Van Helden MJ, et al. FlowSOM: Using self-organizing maps for visualization and interpretation of cytometry data. Cytometry A. 2015;87(7):636–45. Dey S, Eapen MS, Chia C, Gaikwad AV, Wark PAB, Sohal SS. Pathogenesis, clinical features of asthma COPD overlap, and therapeutic modalities. Am J Physiol Lung Cell Mol Physiol. 2022;322(1):L64–L83. Dey S, Lu W, Haug G, et al. Airway inflammatory changes in the lungs of patients with asthma-COPD overlap (ACO): a bronchoscopy endobronchial biopsy study. Respir Res. 2023;24(1):221. Jeffery PK, Laitinen A, Venge P. Biopsy markers of airway inflammation and remodelling. Respiratory medicine 2000; 94 Suppl F: S9-15. Eapen MS, Myers S, Walters EH, Sohal SS. Airway inflammation in chronic obstructive pulmonary disease (COPD): a true paradox. Expert Rev Respir Med. 2017;11(10):827–39. Eapen MS, Hansbro PM, McAlinden K, et al. Abnormal M1/M2 macrophage phenotype profiles in the small airway wall and lumen in smokers and chronic obstructive pulmonary disease (COPD). Sci Rep. 2017;7(1):13392. Gorgoulis V, Adams PD, Alimonti A, et al. Cellular Senescence: Defining a Path Forward. Cell. 2019;179(4):813–27. Barnes PJ, Baker J, Donnelly LE. Cellular Senescence as a Mechanism and Target in Chronic Lung Diseases. Am J Respir Crit Care Med. 2019;200(5):556–64. Mitsi E, Kamng'ona R, Rylance J, et al. Human alveolar macrophages predominately express combined classical M1 and M2 surface markers in steady state. Respir Res. 2018;19(1):66. Park HS, Ju UI, Park JW, et al. PPARgamma neddylation essential for adipogenesis is a potential target for treating obesity. Cell Death Differ. 2016;23(8):1296–311. Chen Z, Ghavimi SAA, Wu M et al. PPARgamma agonist treatment reduces fibroadipose tissue in secondary lymphedema by exhausting fibroadipogenic PDGFRalpha + mesenchymal cells. JCI Insight 2023; 8(24). Mosana M, Ayeleso A, Nyakudya T, Erlwanger K, Mukwevho E. Potential Protective Effects of Neonatal Supplementation with Oleanolic Acid on Peroxisome Proliferator-Activated Receptor Gamma (PPARγ)-Ligand Dependent Regulation of Glucose Homeostasis in High Fructose-Fed Rats. Nat Prod Commun 2020; 15(3). Stark JM, Coquet JM, Tibbitt CA. The Role of PPAR-gamma in Allergic Disease. Curr Allergy Asthma Rep. 2021;21(11):45. Wei R, Han X, Li M, et al. The nuclear cytokine IL-37a controls lethal cytokine storms primarily via IL-1R8-independent transcriptional upregulation of PPARgamma. Cell Mol Immunol. 2023;20(12):1428–44. Yu L, Gao Y, Aaron N, Qiang L. A glimpse of the connection between PPARgamma and macrophage. Front Pharmacol. 2023;14:1254317. Liu X, Gu Y, Kumar S, et al. Oxylipin-PPARgamma-initiated adipocyte senescence propagates secondary senescence in the bone marrow. Cell Metab. 2023;35(4):667–84. e6. Toobian D, Ghosh P, Katkar GD. Parsing the Role of PPARs in Macrophage Processes. Front Immunol. 2021;12:783780. Yang HC, Deleuze S, Zuo Y, Potthoff SA, Ma LJ, Fogo AB. The PPARgamma agonist pioglitazone ameliorates aging-related progressive renal injury. J Am Soc Nephrol. 2009;20(11):2380–8. Villapol S. Roles of Peroxisome Proliferator-Activated Receptor Gamma on Brain and Peripheral Inflammation. Cell Mol Neurobiol. 2018;38(1):121–32. Daniel B, Nagy G, Czimmerer Z et al. The Nuclear Receptor PPARgamma Controls Progressive Macrophage Polarization as a Ligand-Insensitive Epigenomic Ratchet of Transcriptional Memory. Immunity 2018; 49(4): 615 – 26 e6. Additional Declarations No competing interests reported. Supplementary Files FigureS1.tif FigureS2.tif FigureS3.tif WanetalSupplementalTables12192023.xlsx Cite Share Download PDF Status: Published Journal Publication published 20 Apr, 2024 Read the published version in Respiratory Research → Version 1 posted Editorial decision: Revision requested 19 Mar, 2024 Reviews received at journal 17 Mar, 2024 Reviewers agreed at journal 04 Mar, 2024 Reviewers invited by journal 04 Mar, 2024 Editor assigned by journal 04 Mar, 2024 Submission checks completed at journal 04 Mar, 2024 First submitted to journal 03 Mar, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4009724","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":276338234,"identity":"ed7b1d97-4897-44c6-9472-84992f92f497","order_by":0,"name":"Rongjun Wan","email":"","orcid":"","institution":"Johns Hopkins University School of Medicine","correspondingAuthor":false,"prefix":"","firstName":"Rongjun","middleName":"","lastName":"Wan","suffix":""},{"id":276338235,"identity":"979a9986-6795-4745-85be-81ff3ef77148","order_by":1,"name":"Prakhyath Srikaram","email":"","orcid":"","institution":"Johns Hopkins University School of Medicine","correspondingAuthor":false,"prefix":"","firstName":"Prakhyath","middleName":"","lastName":"Srikaram","suffix":""},{"id":276338236,"identity":"15c32370-2865-49c3-8e12-aa3151608681","order_by":2,"name":"Shaobing Xie","email":"","orcid":"","institution":"Johns Hopkins University School of Medicine","correspondingAuthor":false,"prefix":"","firstName":"Shaobing","middleName":"","lastName":"Xie","suffix":""},{"id":276338237,"identity":"f1f6eb29-b205-4c66-a6b1-bd4a10073c86","order_by":3,"name":"Qiong Chen","email":"","orcid":"","institution":"Central South University","correspondingAuthor":false,"prefix":"","firstName":"Qiong","middleName":"","lastName":"Chen","suffix":""},{"id":276338238,"identity":"3696eba0-e78d-42f6-ae0d-2e1800339073","order_by":4,"name":"Chengping Hu","email":"","orcid":"","institution":"Central South University","correspondingAuthor":false,"prefix":"","firstName":"Chengping","middleName":"","lastName":"Hu","suffix":""},{"id":276338239,"identity":"fb56dd89-83d6-4f63-957c-8241860912c9","order_by":5,"name":"Mei Wan","email":"","orcid":"","institution":"Johns Hopkins University School of Medicine","correspondingAuthor":false,"prefix":"","firstName":"Mei","middleName":"","lastName":"Wan","suffix":""},{"id":276338240,"identity":"fb2531c1-76e0-4893-a6a6-4be73bda00c9","order_by":6,"name":"Yuanyaun Li","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAyUlEQVRIiWNgGAWjYBACPijN2A+hmQlrYYNpmdlAspYNB4jWwt587MPHtjuym48ffybBUGGd2MB+9gB+LTzHkmfObHtmvO1MjpkEw5n0xAaevAT8WiRyjJl52w4nbrvBwybBCGQ0SPAY4Nci//4zWMvmGezPJBj/EaNFgocZrGWDBIOZBGMDMVp40owZZ5w7bDzjTI6xRcKxdOM2nhz8WvjZDz9m+FB2WLa//fjDGx9qrGX72c/g14IKEhgQMTUKRsEoGAWjgAIAAFZXQENOz5wIAAAAAElFTkSuQmCC","orcid":"","institution":"Central South University","correspondingAuthor":true,"prefix":"","firstName":"Yuanyaun","middleName":"","lastName":"Li","suffix":""},{"id":276338241,"identity":"c040970f-5d7e-4ead-ab54-5873bcdc59d0","order_by":7,"name":"Peisong Gao","email":"","orcid":"","institution":"Johns Hopkins University School of Medicine","correspondingAuthor":false,"prefix":"","firstName":"Peisong","middleName":"","lastName":"Gao","suffix":""}],"badges":[],"createdAt":"2024-03-03 21:59:58","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4009724/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4009724/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12931-024-02790-6","type":"published","date":"2024-04-20T23:00:23+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":52036654,"identity":"bc0b3702-08e4-446a-92e3-7920b9473abd","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":128716,"visible":true,"origin":"","legend":"\u003cp\u003eAlveolar macrophages were the major cell types in patients with ACO\u003cem\u003e. \u003c/em\u003e\u003cstrong\u003eA\u003c/strong\u003e, A total of 12 cell clusters were identified within the scRNA-Seq dataset generated from human lung tissues of patients with ACO (n=1) and without ACO (n=2). \u003cstrong\u003eB\u003c/strong\u003e, Cell clusters in (\u003cstrong\u003eA\u003c/strong\u003e) were further annotated as different cell types by using the SingleR algorithm. \u003cstrong\u003eC\u003c/strong\u003e, Heatmap represents the top-ranking differentially expressed genes for each cell type in (\u003cstrong\u003eB\u003c/strong\u003e). \u003cstrong\u003eD\u003c/strong\u003e, Relative percentage of each cell type in ACO and controls. \u003cstrong\u003eE\u003c/strong\u003e, A total of seven clusters were identified within monocytes/macrophages by using the Leiden algorithm. \u003cstrong\u003eF\u003c/strong\u003e, Cell clusters in (\u003cstrong\u003eE\u003c/strong\u003e) were further annotated as different sub-types of monocytes/macrophages. \u003cstrong\u003eG\u003c/strong\u003e, Heatmap represents the top-ranking differentially expressed genes for each cell type in (\u003cstrong\u003eF\u003c/strong\u003e). ACO: asthma-COPD overlap. Chi-Square tests were performed for the comparisons of category variables.\u003c/p\u003e","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/5bbcf2faba23fa92fcb757fd.png"},{"id":52036659,"identity":"4c208ccd-3eaa-4800-b6a2-913d108ff3bb","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":303509,"visible":true,"origin":"","legend":"\u003cp\u003eDeceased cellular senescence in monocytes/macrophages of patients with ACO. \u003cstrong\u003eA\u003c/strong\u003e, Density of AUCell enriched SenMayo senescence score within lung monocytes/macrophages. \u003cstrong\u003eB\u003c/strong\u003e, Distribution of ACO and control group among monocytes/macrophages. \u003cstrong\u003eC\u003c/strong\u003e, Enrichment scores of SenMayo senescence in ACO and control group. \u003cstrong\u003eD\u003c/strong\u003e, Relative percentage of each monocyte/macrophage subtype in ACO and controls. \u003cstrong\u003eE\u003c/strong\u003e, Relative percentage of each cell cycle phases in ACO and controls. \u003cstrong\u003eF-G\u003c/strong\u003e, Expression of cellular senescence markers CDKN1A (p21, \u003cstrong\u003eF\u003c/strong\u003e) and CDKN2A (p16, \u003cstrong\u003eG\u003c/strong\u003e) in ACO and controls. Chi-Square tests were performed for the comparisons of category variables. Wilcox rank sum tests were used for the comparison of numerical variables. A p \u0026lt; 0.05 was considered as statistical significance.\u003c/p\u003e","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/210a84ec31ee3b98b897c5d4.png"},{"id":52036655,"identity":"2964c45a-1641-4ade-bcbd-5f9516720d97","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":45865,"visible":true,"origin":"","legend":"\u003cp\u003eA lower prevalence of senescence in patients with severe asthma in a IMSA dataset. \u003cstrong\u003eA\u003c/strong\u003e, Identification of stable and meaningful clusters by repeated consensus clustering analysis with different values of k (from 2 to 5) within a IMSA dataset. \u003cstrong\u003eB\u003c/strong\u003e, Consistent and stable cellular senescence clusters identified by Cumulative Distribution Function (CDF). \u003cstrong\u003eC\u003c/strong\u003e, Two clusters were assigned senescence low and high group according to the ssGSEA SenMayo Senescence score. \u003cstrong\u003eD-F\u003c/strong\u003e, Relative expression of senescence markers CDKN1A (p21, \u003cstrong\u003eD\u003c/strong\u003e), CDKN2A (p16, \u003cstrong\u003eE\u003c/strong\u003e), and proliferation marker MKI67\u003cem\u003e \u003c/em\u003e(\u003cstrong\u003eF\u003c/strong\u003e) in senescence clustering low and high group. \u003cstrong\u003eG\u003c/strong\u003e, Relative percentage of health control (HC, n=5), mild/moderate asthma patients (MMA, n=17), and severe asthma (SA, n=16) participants in senescence clustering low and high group. Fisher’s Exact test was performed for the comparisons of category variables. Wilcox rank sum tests were used for the comparison of numerical variables between 2 groups. A p \u0026lt; 0.05 was considered as statistical significance.\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/4612947a97309c49d39058db.png"},{"id":52036658,"identity":"fd9fa1ba-a498-4df8-a13e-c5a027abdaa2","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":196190,"visible":true,"origin":"","legend":"\u003cp\u003eAirway CD206\u003csup\u003e+\u003c/sup\u003emacrophages show a lower prevalence of senescence but increased cytokines. \u003cstrong\u003eA\u003c/strong\u003e, Density scatter plots showed the distribution of relative expression of markers for cell type annotation. \u003cstrong\u003eB\u003c/strong\u003e, A total of 7 different cell clusters annotated using a combination of surface marker genes. \u003cstrong\u003eC\u003c/strong\u003e, Relative percentage of each cell type in senescence clustering low and high group. \u003cstrong\u003eD\u003c/strong\u003e, Levels of cytokines in BAL fluids in senescence clustering low and high group of the IMSA cohort. \u003cstrong\u003eE\u003c/strong\u003e, Heatmap represents cytokine expression in CD206\u003csup\u003e+\u003c/sup\u003e macrophages of senescence clustering low and high group. Variables. Limma based multiple linear regressions were performed for abundance and expression comparisons between senescence high and low groups. A p \u0026lt; 0.05 was considered as statistical significance.\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/829ad394e790a3216f49d4ec.png"},{"id":52036968,"identity":"0ba5bb77-dfd0-46dc-a6ef-1d790e85de20","added_by":"auto","created_at":"2024-03-05 17:16:43","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":960582,"visible":true,"origin":"","legend":"\u003cp\u003eDifferentially expressed genes and pathways in alveolar macrophages of ACO. \u003cstrong\u003eA\u003c/strong\u003e, A two-dimensional scatter plot of trajectory estimation of all cells from monocytes/macrophages. \u003cstrong\u003eB\u003c/strong\u003e, Distribution of cell clusters identified by Leiden matching trajectory analysis. \u003cstrong\u003eC\u003c/strong\u003e, Distribution of cells strata by groups matching trajectory analysis. \u003cstrong\u003eD\u003c/strong\u003e, Relative expression of lineage markers commonly used for monocytes/macrophages\u003cstrong\u003e. E\u003c/strong\u003e, Top-ranking high expressed marker genes across 7 cell clusters identified by FindMarkers. \u003cstrong\u003eF\u003c/strong\u003e, Difference of AUCell KEGG enrichment scores between Leiden cluster 1 and other clusters. \u003cstrong\u003eG\u003c/strong\u003e, Difference in SenMayo senescence score between ACO and control group. \u003cstrong\u003eH-I\u003c/strong\u003e, Relative expression of senescence markers CDKN1A (p21, \u003cstrong\u003eH\u003c/strong\u003e) and CDKN2A (p16, \u003cstrong\u003eI\u003c/strong\u003e) in ACO and control group. \u003cstrong\u003eJ\u003c/strong\u003e, Heatmap representing the differential expression of genes linked to oxidoreductases, cytokines, and growth factors between ACO and control group. \u003cstrong\u003eK\u003c/strong\u003e, Correlations between CDKN1A expression and differentially expressed genes in Leiden cluster 1 and lineage markers for monocytes/macrophages. Wilcox rank sum tests were used for the comparison of numerical variables between 2 groups; Pearson’s correlations were used for computing coefficients between imputed expression values. A p \u0026lt; 0.05 was considered as statistical significance.\u003c/p\u003e","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/248bcc998a7b98cbf22ebdfc.png"},{"id":52036656,"identity":"96d21f66-6b21-4ab8-9fb8-3975e3a571be","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":157468,"visible":true,"origin":"","legend":"\u003cp\u003ePPARγ as a key regulating factor of cellular senescence in alveolar macrophages. \u003cstrong\u003eA\u003c/strong\u003e, Six gene modules identified by the hdWGCNA analysis focused on Leiden cluster 1. \u003cstrong\u003eB\u003c/strong\u003e, Brown module was positively correlated with the difference in cellular senescence between ACO and control groups. C, Top 30 hub genes withhigh correlations in the brown module. \u003cstrong\u003eD-F\u003c/strong\u003e, Different pathways in the brown module were identified by WikiPathways (\u003cstrong\u003eD\u003c/strong\u003e), Reactome (\u003cstrong\u003eE)\u003c/strong\u003eand KCGG (\u003cstrong\u003eF\u003c/strong\u003e). \u003cstrong\u003eG\u003c/strong\u003e, Top 10 specific regulons for Leiden cluster 1 ranked by score value generated from PySCENIC. \u003cstrong\u003eH\u003c/strong\u003e, Transcription factor enrichment analysis within the brown module against the TRRUST database.\u003c/p\u003e","description":"","filename":"Figure6.png","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/4168faef00eec2226e9c3254.png"},{"id":55691196,"identity":"2e355a56-4246-41a9-95e4-e4c89a5d8d62","added_by":"auto","created_at":"2024-05-01 23:08:50","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":3100494,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/3622f6ea-139c-4141-95fa-8015ce73d87a.pdf"},{"id":52036662,"identity":"bc15c45c-a902-451c-b0db-535965493346","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"tif","order_by":8,"title":"","display":"","copyAsset":false,"role":"supplement","size":331158,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS1.tif","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/f73537766fb5451ecb554c6f.tif"},{"id":52036664,"identity":"159bf2b3-1b61-414d-b88a-32a9e9bebfc5","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"tif","order_by":9,"title":"","display":"","copyAsset":false,"role":"supplement","size":582260,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS2.tif","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/87cdb2a017801b32f36faf76.tif"},{"id":52036969,"identity":"4cf16421-2b9a-42a1-9fc4-3a6563456591","added_by":"auto","created_at":"2024-03-05 17:16:43","extension":"tif","order_by":10,"title":"","display":"","copyAsset":false,"role":"supplement","size":702242,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS3.tif","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/859fecc355fcf6a3eb4752be.tif"},{"id":52036660,"identity":"e01ba548-3977-45f9-b5f9-a7b66a9fe38a","added_by":"auto","created_at":"2024-03-05 17:08:43","extension":"xlsx","order_by":11,"title":"","display":"","copyAsset":false,"role":"supplement","size":71503,"visible":true,"origin":"","legend":"","description":"","filename":"WanetalSupplementalTables12192023.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-4009724/v1/dff6f68f7919afb96c2c1d8c.xlsx"}],"financialInterests":"No competing interests reported.","formattedTitle":"PPARγ Attenuates Cellular Senescence of Alveolar Macrophages in Asthma- COPD Overlap","fulltext":[{"header":"INTRODUCTON","content":"\u003cp\u003eAsthma and chronic obstructive pulmonary disease (COPD) are two of the most prevalent obstructive lung diseases. Asthma-COPD overlap (ACO) is a term used to describe a complex condition characterized by shared clinical and pathophysiological features of both diseases, typically observed in older individuals with a long history of disease and exposure to environmental factors\u003csup\u003e\u003cspan additionalcitationids=\"CR2\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e. ACO is estimated to affect approximately 10\u0026ndash;30% of individuals with asthma and around 25% of those with COPD\u003csup\u003e\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u003c/sup\u003e. Cumulative environmental exposures, such as exposure to allergens and toxic particles or gases (e.g., smoking and indoor biomass), are believed to play a crucial role in the development and progression of ACO\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e. Additionally, asthma and atopy have been identified as potential risk factors for the development of COPD and, consequently, ACO\u003csup\u003e\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e\u003c/sup\u003e. More importantly, patients with ACO often experience increased rates of exacerbations and clinical symptoms, posing a significant challenge for clinical treatment and management\u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e,\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. However, managing ACO patients can be intricate due to the overlapping nature of these diseases.\u003c/p\u003e \u003cp\u003eThe inflammatory mechanisms in asthma and COPD are distinct. In asthma, inflammation is primarily driven by Th2 cell-mediated responses, leading to eosinophilic inflammation and an increase in cytokines like IL-4, IL-5, and IL-13\u003csup\u003e9\u003c/sup\u003e. On the other hand, COPD is characterized by inflammation with a higher presence of neutrophils and a dominance of Th1 and Th17 cell-mediated responses \u003csup\u003e\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e. In both conditions, macrophages, the primary immune cells in the lungs, play crucial roles in the immune response, defense against infections, tissue homeostasis, and inflammation resolution\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. Macrophages in asthma and COPD can contribute to airway inflammation and remodeling by releasing pro-inflammatory cytokines, proteases, and reactive oxygen species\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. Therefore, it is reasonable to believe that macrophages also play a vital role in the inflammatory processes and pathophysiology of ACO.\u003c/p\u003e \u003cp\u003eACO is more commonly diagnosed in older individuals, and its prevalence tends to increase with age, as supported by various studies\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e,\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e. Age has been identified as a significant risk factor for ACO in asthma patients\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e. Consequently, the aging process and senescence appear to play substantial roles in the development of lung inflammation in individuals with ACO. Senescence represents a complex cellular state characterized by cellular stress, DNA damage, cell cycle arrest, and the release of senescence-associated secretory phenotype (SASP) factors\u003csup\u003e\u003cspan additionalcitationids=\"CR18\" citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e. These SASP factors encompass chemokines, cytokines, growth factors, adhesion molecules, and lipid components that can contribute to multiple age-related disorders with both local and systemic consequences\u003csup\u003e\u003cspan additionalcitationids=\"CR21 CR22 CR23 CR24 CR25 CR26\" citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. It's worth noting that cellular senescence has been associated with both asthma\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e,\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e and COPD\u003csup\u003e\u003cspan additionalcitationids=\"CR31\" citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Consequently, senescence may exert a significant influence on the development and management of ACO. However, it is essential to recognize that senescence is a complex process influenced by factors such as cell type, age, and specific diseases\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e,\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. Therefore, it is crucial to investigate whether senescence in lung macrophages contributes to the inflammatory processes and pathophysiology unique to ACO.\u003c/p\u003e \u003cp\u003eIn our current research, we conducted an analysis using publicly available single-cell RNA sequencing (scRNA-Seq) data obtained from human lung tissues, comparing individuals with and without ACO. Our primary focus was on identifying the major cell types and identified monocytes/macrophages as the predominant cell types in patients with ACO. We delved into gene signatures associated with senescence, specifically within monocytes/macrophages. The findings were further validated through another independent cohort (IMSA) with a specific emphasis on investigating the relationship between cellular senescence features and the severity of asthma. Additionally, we explored differentially expressed genes and pathways within alveolar macrophages of individuals with ACO and identified PPARγ as a key regulatory factor in driving cellular senescence in alveolar macrophages.\u003c/p\u003e"},{"header":"METHODS","content":"\u003cp\u003e \u003cb\u003eData Source\u003c/b\u003e \u003c/p\u003e \u003cp\u003eThe scRNA-Seq data was publicly available and was generated from human lung tissues of a patient with ACO who died with exacerbation and two transplant donors\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. The Bulk RNA-Seq and mass cytometry (cytometry by time of flight, CyTOF) data were derived from bronchoalveolar lavage fluid (BALF) samples and downloaded from the Gene Expression Omnibus (GEO) database with the accession number GSE136587, and the Flow Repository (FR) database with the identifier FR-FCM-Z395\u003csup\u003e35\u003c/sup\u003e. The data were collected as a part of the Immune Mechanisms Severe Asthma (IMSA) study and involved the analysis of bronchoscopically obtained distal airway and alveolar cells. The study included a total of 39 subjects consisting of 6 healthy controls, 17 mild/moderate asthma patients, and 16 severe asthma patients.\u003c/p\u003e \u003cp\u003e \u003cb\u003eSingle Cell RNA-Seq Processing and Analyzing\u003c/b\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eQuality Control\u003c/b\u003e \u003c/p\u003e \u003cp\u003eThe raw counts matrixes downloaded were loaded to R workspace to initiate a Seurat workflow. A series of critical preprocessing steps were performed to ensure data quality before the analysis. Specifically, cells with a count of unique molecular identifiers (UMIs) below 500, cells expressing fewer than 200 genes with a log10GenesPerUMI value below 0.75, and cells with mitochondrial gene expression contributing to over 20% of the total gene expression were removed.\u003c/p\u003e \u003cp\u003e \u003cb\u003eDimension Reduction and Cell Annotation\u003c/b\u003e \u003c/p\u003e \u003cp\u003eFollowing the quality control procedures, the processed data underwent integration using distinct strategies: the complete dataset was integrated utilizing Seurat\u003csup\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e, while the Monocytes/Macrophages subset was integrated using the Harmony approach\u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e\u003c/sup\u003e. The dimensionality reduction and visualization were achieved using the Uniform Manifold Approximation and Projection (UMAP) algorithm. UMAP not only captures inter-cellular distances but also provides a holistic view of the data's global structure\u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e. Cell clusters were identified by Leiden algorithm, an automated algorithm tailored for effectively clustering cells in scRNA-seq data\u003csup\u003e\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e. Cell markers were these with differentially expressed genes (DEGs) identified by the FindAllMarkers function from the Seurat package. Cellular annotations were identified by SingleR or manually based on specific marker expression, which allows us to get the composition and dynamics of cells within the dataset\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eDifferential Expressing Analyses\u003c/b\u003e \u003c/p\u003e \u003cp\u003eDifferential expressing analyses were conducted using the Libra R package (Version 1.0.0). Within this framework, the \"run_de\" function was employed, applying Wilcoxon rank sum tests to assess statistical differences between groups. Genes with an adjusted P value\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and absolute log2-fold change (|log2FC|)\u0026thinsp;\u0026gt;\u0026thinsp;0.5 were termed as statistical significance. Moreover, the Seurat package's \"FindMarkers\" and \"FindAllMarkers\" functions were also utilized to identify markers that distinguish cell clusters with Wilcoxon rank sum tests.\u003c/p\u003e \u003cp\u003e \u003cb\u003eEnrichment Analyses\u003c/b\u003e \u003c/p\u003e \u003cp\u003eThe enrichment analyses included geneset scoring through the AUCell algorithm, and over-representation analysis (ORA). The AUCell algorithm was employed to gauge cellular senescence enrichment. This algorithm calculated an enrichment score by juxtaposing the input data against the SenMayo senescence geneset, an optimal choice for senescence screening since it consistently aligns with results from biological validation experiments\u003csup\u003e\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. The ORA analysis function RunEnricher was embedded within the High-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA) R package. It was employed to investigate the biological characteristics of selected co-expression modules.\u003c/p\u003e \u003cp\u003e \u003cb\u003eTrajectory Analyses\u003c/b\u003e \u003c/p\u003e \u003cp\u003eA partition-based graph abstraction (PAGA) based trajectory analysis was used to investigate the trajectory of monocytes/macrophages within single-cell RNA sequencing (scRNA-seq) data. This procedure was performed using Scanpy version 1.9.1 in Python 3.10.9., which revealed the intricate relationships and transitions among different cell states. This workflow simplified the dataset's graph representation by clustering similar cells, enabling the identification of connected and disconnected regions. By combining well-defined paths rooted in high-confidence connections and leveraging a random-walk-based distance metric, cells are systematically ordered within each partition according to their proximity to a designated reference cell\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eHigh-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA)\u003c/b\u003e \u003c/p\u003e \u003cp\u003eTo explore the hub genes and the underlying expression network associated with cellular senescence and related biological processes, we harnessed the power of hdWGCNA (Version 0.2.18). This algorithm is a refinement of the traditional WGCNA designed specifically for scRNA-Seq data. It operates by computing a sparsity-reduced expression matrix termed metacell, akin to the concept of pseudobulk to mitigate the influence of undetected genes and facilitate the identification of gene co-expression modules with clinical or biological information\u003csup\u003e\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e,\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003cp\u003e \u003cb\u003eTranscription Factor Analyses\u003c/b\u003e \u003c/p\u003e \u003cp\u003eTranscription factors are crucial for the initial stage of decoding DNA sequencing and regulation of gene expression\u003csup\u003e\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e. To initiate this process, we employed PySCENIC (Version 0.12.1)\u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e to identify cell cluster-specific transcription factors. Moreover, we extended our investigation by extracting a gene list from the cellular senescence-related gene module identified using hdWGCNA. Subsequently, this gene list was subjected to transcription factor enrichment analysis using the Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining (TRRUST) database on the metascape website (metascape.org)\u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e,\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eBulk RNA-Seq Processing and Analyzing\u003c/b\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eData Preprocessing\u003c/b\u003e \u003c/p\u003e \u003cp\u003eThe expression matrix encompassing raw counts and pertinent phenotype information was procured from the GEO website, facilitated by the GEOQuery R package (Version 2.66.0)\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. The raw counts with batch details were downloaded for differential expressing analysis. Furthermore, the raw counts were normalized and then transformed into log2-transformed transcripts per million (log2TPM) values. To mitigate the batch-induced discrepancies, the Combat algorithm embedded within the sva R package (Version 3.44.0)\u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eUnsupervised Clustering and Differential Expressing Analyses\u003c/b\u003e \u003c/p\u003e \u003cp\u003eDue to heterogeneity of asthma, we did not make direct comparisons for senescence signature between different asthma groups. Instead, we employed an unsupervised clustering approach that focused on the expression profiles of genes associated with cellular senescence\u003csup\u003e\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e\u003c/sup\u003e. For the clustering process, we utilized the k-means algorithm, which groups data points based on their similarity as measured by the Euclidean distance with the ConsensusClusterPlus R package (Version 1.60.0)\u003csup\u003e\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e,\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e\u003c/sup\u003e. The differential expressing analysis was performed to explore the variations in gene expression across these clusters by following the limma-voom approach with batch information adjusted in the model\u003csup\u003e\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eCyTOF Processing and Analyzing\u003c/b\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eData Preprocessing\u003c/b\u003e \u003c/p\u003e \u003cp\u003eAn outlier sample displaying an unusually high intensity value of 3 in the live_dead channels was identified and subsequently excluded from the CyTOF analysis. Cells exhibiting extremely low CD45 expression were also filtered out using Flowjo. The resultant fcs data was then transferred to an R workspace for subsequent analysis. To ensure accurate analysis, the data underwent a bead-based normalization process utilizing the \"normCytof\" function within the CATALYST R package (Version 1.20.1)\u003csup\u003e\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eDimension Reduction and Cell Annotation\u003c/b\u003e \u003c/p\u003e \u003cp\u003eSimilar to scRNA-Seq, dimension reduction for visualization and cell annotation were required for CyTOF core analyses. Especially, the t-Distributed Stochastic Neighbor Embedding (t-SNE) algorithm was employed\u003csup\u003e\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e\u003c/sup\u003e. For automated cluster identification, the FlowSOM\u003csup\u003e\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e\u003c/sup\u003e was used. Cell type annotations were based on expression profiles of major cell type markers similar to that described in relevant publications\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eDifferential Abundance (DA) and Differential States (DS) Analyses Based on CyTOF Data\u003c/b\u003e \u003c/p\u003e \u003cp\u003eDifferential cell proportion (abundance) and cytokine expression (states) between cellular senescence clusters identified by unsupervised clustering were calculated with diffcyt R package (1.16.0). Both limma based DA (diffcyt-DS-voom), and DS (diffcyt-DA-trend) parameters were chosen for differential expression analyses\u003csup\u003e\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e,\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eStatistical Analysis and Data Visualization\u003c/b\u003e \u003c/p\u003e \u003cp\u003eAll analysis and visualizations were processed with R (Version 4.1.2 for scRNA-Seq and Version 4.2.1 for CyTOF and Bulk RNA-Seq), Python (Version 3.10.9) and GraphPad (Version 9.3). Wilcoxon rank sum tests were employed for continuous variables between two groups; Chi-Square tests were performed for the comparisons of categorical variables. Correlation analyses (Pearson) were based on Single-Cell Variational Inference (scVI-tools, Version 0.20.2) imputed expression data to avoid sparsity influence\u003csup\u003e\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e\u003c/sup\u003e. Data visualizations were performed with ggplot2 (Version 3.4.1), genekitr (Version 1.0.5), ggVennDiagram (Version 1.2.2) and ComplexHeatmap (Version 2.10.0) R packages.\u003c/p\u003e"},{"header":"RESULTS","content":"\u003cp\u003e \u003cem\u003eAlveolar Macrophages were the Major Cell types in Patients with ACO\u003c/em\u003e \u003c/p\u003e \u003cp\u003eTo better understand the pathogenesis of ACO, we analyzed the existing scRNA-Seq dataset generated from human lung tissues of a patient with ACO and two transplant donors without ACO\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. By using the Leiden algorithm, an automated algorithm designed to effectively cluster cells in scRNA-seq data\u003csup\u003e\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e, we identified a total of 12 cell clusters within dataset (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, A). These cell clusters were further annotated as different cell types by using the SingleR algorithm with manual adjustments. Cell markers were these differentially expressed across different cell types identified by the FindAllMarkers function from the Seurat package\u003csup\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e. A total of nine different cell types were identified within the dataset, which included monocytes/macrophages, T cells, NK cells, AT2 alveolar type II cells, endothelial cells, airway epithelial cells, B cells, fibroblasts, and mast cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, B). The top-ranking differentially expressed genes for each cell type were presented in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, C and Figure E1). For example, complement-related genes (\u003cem\u003eC1QA, C1QB\u003c/em\u003e), \u003cem\u003eCD68\u003c/em\u003e, and \u003cem\u003eAPOC1\u003c/em\u003e were highly expressed in monocytes/macrophages (Figure E1, A), while genes like \u003cem\u003eSFTPB, SFTPA1\u003c/em\u003e, and \u003cem\u003eSFTPA2\u003c/em\u003e were uniquely and significantly expressed in AT2 airway epithelial cells ( Figure E1, B). Tryptase genes (\u003cem\u003eTPSAB1, TPSB2\u003c/em\u003e) were highly expressed in mast cells ( Figure E1, C). To identify the distribution of different cell types in the samples from ACO and control group, we performed cell proportion analysis as defined by relative percentages (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, D). Of these, monocytes/macrophages were a predominant cell type, constituting more than 50% of the cells in the different groups being analyzed. Furthermore, the proportion of monocytes/macrophages was significantly higher in ACO patient compared to control group (\u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001). Next, we focused on those monocytes/macrophages specifically and re-clustered them by using the Leiden algorithm. A total of seven clusters were identified (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, E). These cell clusters were further annotated as different sub-types of monocytes/macrophages (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, F), including alveolar macrophages (AM), cycling cells, interstitial macrophages (IM), and monocytes. The top-ranking differentially expressed genes for each cell type were presented in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, G and Figure E2). Similar to monocytes/macrophages, the subtype AMs express the complement-related genes \u003cem\u003eC1QA, C1QB\u003c/em\u003e ( Figure E2, A). Additionally, several other genes \u003cem\u003eRBP4, CD9, SERPING1\u003c/em\u003e, and \u003cem\u003eCES1\u003c/em\u003e are also significantly expressed in AMs. The proliferation markers \u003cem\u003ePCLAF, TOP2A\u003c/em\u003e and \u003cem\u003eMKI67\u003c/em\u003e were highly expressed in cycling cells ( Figure E2, B). IMs were in an intermediate state in the UMAP plot and expressed high levels of molecules such as \u003cem\u003eLGMN\u003c/em\u003e, \u003cem\u003eRNASE1\u003c/em\u003e, and \u003cem\u003eCCL2\u003c/em\u003e ( Figure E2, C). While \u003cem\u003eCFP, FCN1\u003c/em\u003e, and \u003cem\u003eS100A8\u003c/em\u003e,were highly expressed in monocytes ( Figure E2, D). Taken together, the results suggest that monocytes/macrophages in ACO are major cells that may drive airway inflammation.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eDecreased cellular Senescence in Monocytes/Macrophages of Patients with ACO\u003c/em\u003e \u003c/p\u003e \u003cp\u003eCellular senescence has gained considerable attention across various diseases, including respiratory conditions such as asthma\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. To determine whether senescence contributes to the development of ACO, we specifically assessed the enrichment of genes associated with senescence within lung monocytes/macrophages. Employing the AUCell algorithm, we computed enrichment scores against a senescence-related gene set known as SenMayo\u003csup\u003e\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. Among all the subtypes, AMs/monocytes showed the higher enrichment sore as assessed by the density of AUCell enriched SenMayo senescence scores (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, A). Next, we investigated the distribution of ACO and control group among monocytes/macrophages (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, B) and enrichment scores of SenMayo senescence in ACO and control group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, C). As illustrated in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, C, the senescent signatures are enriched in the control group as compared to ACO group. Among them, AMs were a predominant cell type, constituting more than 50% of the cells among all these subtypes (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, D). The proportion of AMs was significantly higher in ACO patients as compared to control group (\u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001). Additionally, cell cycle arrested in the G1 phase is one of the important features of cellular senescence\u003csup\u003e\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e\u003c/sup\u003e. We found that monocytes/macrophages that were in the G1 phase of the cell cycle were more abundant in the control group compared with ACO group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, E). These findings were further supported by the expression of cellular senescence markers, CDKN1A (p21) and CDKN2A (p16)\u003csup\u003e\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e\u003c/sup\u003e. The expression of CDKN1A was much lower in the ACO group as compared to the control group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, F). In contrast, no significant change was observed for CDKN2A (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, G). These findings suggest a lower prevalence of cellular senescence within lung monocytes/macrophages among patients with ACO.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eA lower prevalence of senescence observed in patients with severe asthma\u003c/em\u003e \u003c/p\u003e \u003cp\u003eWe further validated the results through another independent cohort (IMSA)\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e and analyzed the relationship between cellular senescence features and the severity of asthma. After stringent quality control, a total of 39 subjects were finally included for analysis: 6 healthy controls (HC), 17 mild/moderate asthma patients (MMA), and 16 severe asthma patients (SA). Definition of these subjects has been previously reported\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u003c/sup\u003e. Consensus clustering analysis with the Consensus Cumulative Distribution Function (CDF) identified cellular senescence clusters for immune cells from IMSA. Consensus clustering analysis involves repeated clustering of the data with different values of k (from 2 to 5) within a dataset to identify stable and meaningful clusters (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, A). Based on the consensus CDF (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, B), it is evident that cellular senescence clusters 2 provides the most consistent and stable clustering solution for the IMSA dataset. According to the ssGSEA SenMayo Senescence score, the cluster 2 was divided into senescence clustering low and high group (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.002, Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, C). Senescence clustering low group showed lower expression of senescence markers CDKN1A as compared to senescence high group (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.011, Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, D). No difference was noted for CDKN2A between senescence clustering low and high group (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.29, Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, E). In contrast, senescence low group showed higher expression of cell proliferation marker MKI67 (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.0015, Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, F). Intriguingly, among the analyzed subjects consisting of HC, MMA, and SA participants, the large proportion of severe asthmatic patients was observed in senescence low group as compared to senescence high group (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05, Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, G). These findings indicate that cellular senescence may be negatively associated with the severity of asthma.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eAirway CD206\u003c/em\u003e \u003csup\u003e \u003cem\u003e+\u003c/em\u003e \u003c/sup\u003e \u003cem\u003emacrophages show a lower prevalence of senescence and increased cytokines\u003c/em\u003e \u003c/p\u003e \u003cp\u003eTo further explore the relationship between cellular senescence and airway inflammation in asthma, we focused on the CyTOF dataset that targets lineage markers for the adaptive and innate immune systems\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e. Unsupervised cell clustering was performed in BAL fluids from HC, MMA, and SA participants using the FlowSOM algorithm\u003csup\u003e\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e,\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e\u003c/sup\u003e. Cell types in those clusters were determined by the surface marker staining intensities across t-SNE spaces. A total of 7 different cell clusters were annotated using a combination of surface marker genes, including B lymphocytes (CD19\u003csup\u003e+\u003c/sup\u003e), CD206\u003csup\u003e\u0026minus;\u003c/sup\u003e macrophages (CD11C\u003csup\u003e+\u003c/sup\u003eCD206\u003csup\u003e\u0026minus;\u003c/sup\u003e), CD206\u003csup\u003e+\u003c/sup\u003e macrophages (CD11C\u003csup\u003e+\u003c/sup\u003eCD206\u003csup\u003e\u0026minus;\u003c/sup\u003e), CD4 T lymphocytes (CD3\u003csup\u003e+\u003c/sup\u003eCD4\u003csup\u003e+\u003c/sup\u003e), CD8 T lymphocytes (CD3\u003csup\u003e+\u003c/sup\u003eCD8\u003csup\u003e+\u003c/sup\u003e), γδ T lymphocytes (TCRγδ\u003csup\u003e+\u003c/sup\u003e), and NK cells (CD56\u003csup\u003e+\u003c/sup\u003e) (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e, A, B, and Figure E3). To identify the link between these airway immune cells and cellular senescence, the percentages of those identified cell populations were compared in senescence high and low groups. Among all these populations, a significant difference was found for CD206\u003csup\u003e+\u003c/sup\u003e macrophages, which showed a higher proportion in senescence low in relative to senescence high group (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e, C). Interestingly, the senescence low group showed elevated levels of cytokines IL-5, IL-13, IL-10, and IL-22 in BAL fluids of the IMSA cohort (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e, D). Particularly, CD206\u003csup\u003e+\u003c/sup\u003e macrophages showed increased expression of IL-4, IL-13, and IL-22 in senescence low group as compared those in senescence high group (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e, E). Taken together, these data suggest that CD206\u003csup\u003e+\u003c/sup\u003emacrophages are major airway immune cells associated with cellular senescence and cellular senescence is negatively correlated with cytokines in airway CD206\u003csup\u003e+\u003c/sup\u003eMacrophages.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eDifferentially expressed genes and pathways in Alveolar macrophages of\u003c/em\u003e ACO\u003c/p\u003e \u003cp\u003eGiven the significance of CD206\u003csup\u003e+\u003c/sup\u003emacrophages in cellular senescence, we re-visited the ACO dataset and determined whether the specific subtype of monocytes/macrophages drive the difference in cellular senescence between ACO and control group. To explore this, we used trajectory analysis to infer the relationships and distances between different subtypes of monocytes/macrophages. Analyses were started with cycling cells, we identified a specific cluster that has the highest distance (pseudo-time) and is in the upper right corner of a two-dimensional scatter plot (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, A). This cluster was mostly composed of cells from Leiden cluster 1 (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, B) and was predominantly found in the control group (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, C), suggesting that cells in Leiden cluster 1 may contribute to the difference in biological processes between ACO and control group. To further characterize the Leiden cluster 1, we employed the FindAllMarkers to investigate the highly expressed marker genes for cluster 1 against other cell clusters. A total of 47 genes were identified specifically for the Leiden cluster 1, such as \u003cem\u003eSERPING1, CES1, HLA-DQA1, INHBA, FABP4, LGALS3BP, AKR1C3, CITED2, TERM1, FN1, EVL, ALDH1A1\u003c/em\u003e, and \u003cem\u003ePDLIM1\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, D and see Table E2). Additionally, we checked the lineage markers that are commonly used for monocytes/macrophages (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, E). we found that PPARγ, FCGR1A, FCGR3A, MSR1 (CD205), and MRC1 (CD206) were enriched in Leiden cluster 1. Furthermore, we used AUCell scoring to analyze the biological characteristics of Leiden cluster 1 in comparison to other clusters. Especially, the KEGG enrichment analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, F and see Table E3) indicated that Leiden cluster 1 had significant enrichments in several pathways, such as antigen processing and presentation, Peroxisome Proliferator-Activated Receptor (PPAR) signaling pathway. Next, we investigated the differences in cellular senescence between ACO and control groups specially in Leiden cluster 1. Consistent with our previous analyses among all the monocytes/macrophages, the SenMayo senescence score was remarkably lower in ACO groups compared with control group (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, G). This was further supported by the expression of cellular senescence marker CDKN1A (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, H), which also showed reduced expression in the ACO group. No statistical difference was observed for CDKN2A (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, I). Additionally, to further display the characteristics of cellular senescence, we identified differentially expressed genes related to oxidation-reduction, cytokines, and growth factors between ACO and control groups in Leiden cluster 1. Of those, TGFB1 was increased but several chemokines such as CXCL2, CXCL3, CXCL8, and CCL18 were decreased in ACO group compared with control group (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, J and see Table E4). These results suggest correlations between cellular senescence and cytokine expression in alveolar macrophages of Patients with ACO. Furthermore, strong correlations were identified for CDKN1A expression levels and those identified genes in Leiden cluster 1 and traditional lineage markers for monocytes/macrophages (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, K and see Table E5 in this article\u0026rsquo;s Online Repository at \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e\u003ca href=\"http://www.jaci.org\" target=\"_blank\"\u003ewww.jaci.org\u003c/a\u003e\u003c/span\u003e\u003cspan address=\"http://www.jaci.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). For example, CDNK1A expression was correlated with genes with either positive regulation (e.g., \u003cem\u003eCD14, CD68, CD86, CD163, FCGR3A, MSR1, SERPING1, HLA-DQA1, INHBA, LGALS3BP, AKR1C, TREM1, FN1\u003c/em\u003e, and \u003cem\u003eALDH1A1\u003c/em\u003e) or negative regulation (e.g., \u003cem\u003eMARCO, FCGR1A, ITGB2\u003c/em\u003e, and \u003cem\u003eMRC1\u003c/em\u003e). Collectively, the results indicate that these identified genes and biological pathways may be involved in regulating cellular senescence in alveolar macrophages of patients with ACO.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003ePPARγ, a Key Regulating Factor of Cellular Senescence in Alveolar Macrophages\u003c/em\u003e \u003c/p\u003e \u003cp\u003eTo delve deeper into the underlying regulatory mechanisms of biological changes, we first utilized the hdWGCNA analysis focused on Leiden cluster 1 to explore the expression network associated with senescence signatures. The algorithm discerned six gene modules: green, blue, red, brown, cyan, and tan (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, A). Within these, CDKN1A (p21) was pinpointed in the brown module, and the module feature gene score was positively correlated with the difference in cellular senescence between ACO and control groups (R\u0026thinsp;=\u0026thinsp;0.669, p\u0026thinsp;\u0026lt;\u0026thinsp;0.0001, Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, B). Prominently, the top 30 hub genes with high correlation in this module encompassed those related to growth arrest and proliferation (e.g., \u003cem\u003eGADD45B, PPP1R15A\u003c/em\u003e, and \u003cem\u003eDUSP2\u003c/em\u003e), TP53-mediated cell senescence-associated heat shock proteins (\u003cem\u003eDNAJA1, DNAJB1\u003c/em\u003e, and \u003cem\u003eHSPA5\u003c/em\u003e), DNA damage and repair (e.g., MYL12A), and energy, protein, and lipid metabolism (e.g., \u003cem\u003eGLUL, PLIN2, RPS4Y1, PDIA3\u003c/em\u003e, and \u003cem\u003eSOD2\u003c/em\u003e). Furthermore, genes associated with SASP and its upstream regulators (e.g., \u003cem\u003eCCL20, CXCL3, CXCL5, CXCL8, JUN, NFKB1A\u003c/em\u003e, and \u003cem\u003eNFKBIZ\u003c/em\u003e) and those tied to antigen presentation, antibody response, and macrophage activation (e.g., \u003cem\u003eC83, CD68, FCGR3A\u003c/em\u003e, and \u003cem\u003eB2M\u003c/em\u003e) were also included in these hub genes (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, C). These genes or biological processes are directly or indirectly involved in senescence processes of alveolar macrophages. Subsequent enrichment analyses of entire module genes showed that the brown module was mainly associated with macrophage markers and cellular senescence (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, D), cytokine signaling in immune system (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, E), and antigen processing and presentation and PPAR signaling (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, F). These biological characteristics are closely related and mutually regulated. To pinpoint the central regulator of senescence in alveolar macrophages, we leveraged SCENIC to discern the most specific transcription factors for each Leiden cluster. Notably, for Leiden cluster 1, PPARγ emerged as the predominant regulator (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, G). This was further supported by the transcription factor enrichment analysis with the brown module against the TRRUST database. PPARγ was identified to be one of the top 10 pivotal transcription factors modulating gene expression within this module (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, H). Intriguingly, each of these transcription factors had the potential to regulate CDKN1A expression. Taken together, these findings indicate that the brown module had the core impact on cellular senescence, and PPARγ is one of the predominant regulators modulating the senescent signature of alveolar macrophages.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e"},{"header":"DISCUSSION","content":"\u003cp\u003eBoth asthma and COPD can co-exist in older individuals with clinically overlapping phenotypes\u003csup\u003e\u003cspan additionalcitationids=\"CR2\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e. ACO has recently been a focus of interest because patients with ACO often experience poor health-related quality of life, increased rates of exacerbations, and severe clinical symptoms\u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e,\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. While the exact pathophysiology of ACO is not fully understood, it is thought to result from various features of both asthma and COPD \u003csup\u003e\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e,\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u003c/sup\u003e. Both asthma and COPD involve chronic airway inflammation, but they have different underlying mechanisms. In asthma, there is typically eosinophilic inflammation driven by a Th2 immune response\u003csup\u003e\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u003c/sup\u003e, while COPD is characterized by neutrophilic inflammation associated with Th1 and Th17-driven inflammation\u003csup\u003e\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e,\u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e66\u003c/span\u003e,\u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e67\u003c/span\u003e\u003c/sup\u003e. In both conditions, macrophages play a central role in the recruitment and activation of neutrophils and eosinophils\u003csup\u003e\u003cspan additionalcitationids=\"CR12\" citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e68\u003c/span\u003e\u003c/sup\u003e. Indeed, a very recent cross-sectional exploratory study was performed to investigate cellular changes (e.g., macrophages, neutrophils, eosinophils, mast cells, CD8\u003csup\u003e+\u003c/sup\u003e, and CD4\u003csup\u003e+\u003c/sup\u003eT lymphocytes) in the airway wall of ACO compared with asthma, COPD current smokers, and ex-smokers, normal lung function smokers, and non-smoker controls. They found that the ACO airway tissue inflammatory cellular profile differed from the contributing diseases of asthma and COPD with a predominance of macrophages\u003csup\u003e\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u003c/sup\u003e. The results were supported by our analyses on an existing scRNA-Seq dataset generated from human lung tissues of patients with ACO. Monocytes/macrophages were a predominant cell type among all cell type to be analyzed (monocytes/macrophages, T cells, NK cells, AT2 alveolar type II cells, endothelial cells, airway epithelial cells, B cells, fibroblasts, and mast cells), constituting more than 50% of the total cells, and the proportion of monocytes/macrophages was significantly higher in patients with ACO in relative to control group. Importantly, our further study identified several sub-types of monocytes/macrophages, including alveolar macrophages, cycling cells, interstitial macrophages, and monocytes. Of these, alveolar macrophage was the most predominant cell type with the increased expression of complement-related genes C1QA, C1QB, and several other genes RBP4, CD9, SERPING1, and CES1. Thus, our findings provide further evidence that monocytes/macrophages, particularly alveolar macrophages, are major cells that may contribute to airway inflammation in patients with ACO.\u003c/p\u003e \u003cp\u003eNext, we explored how monocytes/macrophages drive airway inflammation with the emphasis on cellular senescence that has been implicated in the pathophysiology of various diseases, including asthma\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. Cellular senescence is an extremely complex and dynamic biological process induced by several factors, including aging, oxidative stress, DNA damage, mitochondrial dysfunction, epigenetic modifications, and telomere shortening\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e,\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e,\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e\u003c/sup\u003e. Cellular senescence has been associated with both asthma\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e and COPD\u003csup\u003e\u003cspan additionalcitationids=\"CR31\" citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Several stimuli causing cellular senescence have been associated with asthma, such as telomere shortening, oxidative stress, inflammation, and autophagy/mitophagy. Therefore, senescence may exert a significant influence on the function and activation of ACO-associated target cells and subsequently development and management of ACO. Here we investigated the relationship between senescence and ACO by analyzing gene signatures associated with senescence in all the subtypes of monocytes/macrophages. We demonstrated for the first time a lower prevalence of cellular senescence in alveolar macrophages of patients with ACO. Furthermore, CDKN1A (p21) but not CDKN2A (p16) showed lower expression in alveolar macrophages of patients with ACO, suggesting that CDKN1A (p21) pathway may play a major role in mediating senescence of alveolar macrophages in ACO.\u003c/p\u003e \u003cp\u003eTo strengthen our findings on the role of cellular senescence in ACO, we analyzed the relationship of cellular senescence features with the severity of asthma in a total of 39 individuals from IMSA\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e. Specifically, senescence in those individuals were clustered into senescence low and high group and then the distribution of healthy controls, mild/moderate asthma, and severe asthma was analyzed in these two groups. Intriguingly, the large proportion of severe asthmatic patients was observed in senescence clustering low group compared with senescence clustering high group. Furthermore, we focused on the CyTOF dataset from BAL fluids that targets lineage markers and identified a total of 7 different cell clusters, including B lymphocytes, CD206\u003csup\u003e\u0026minus;\u003c/sup\u003e macrophages, CD206\u003csup\u003e+\u003c/sup\u003e macrophages, CD4\u003csup\u003e+\u003c/sup\u003e T, CD8\u003csup\u003e+\u003c/sup\u003e T, γδ T, and NK cells. Of these, CD206\u003csup\u003e+\u003c/sup\u003e macrophages were distributed in senescence clustering low group and showed increased expression of Th2 cytokines IL-4, IL-13 and Th17 cytokine IL-22, suggesting that, CD206\u003csup\u003e+\u003c/sup\u003e macrophages may drive both Th2 and Th17 populations commonly seen in patients with ACO\u003csup\u003e\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e,\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u003c/sup\u003e. CD206 expression is a characteristic feature of some tissue-resident macrophages, including alveolar macrophages\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e,\u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e\u003c/sup\u003e, and CD206\u003csup\u003e+\u003c/sup\u003e macrophages represent a large proportion of alveolar macrophages. Collectively, these independent data provide evidence that patients with severe asthma have a lower prevalence of cellular senescence, and CD206\u003csup\u003e+\u003c/sup\u003emacrophages may contribute to the severity of asthma, possibly ACO, by releasing inflammatory cytokines.\u003c/p\u003e \u003cp\u003eWe also explored the underly mechanisms as for how cellular senescence is regulated in macrophages. We investigated the differentially expressed genes and pathways in alveolar macrophages in Leiden cluster 1, a population mostly likely contributing to the difference in cellular senescence between ACO and control group. Several pathways were also identified, such as antigen processing and presentation, PPARγ signaling pathway. Additionally, several markers commonly used for monocytes/macrophages also showed different expression, including PPARγ, FCGR1A, FCGR3A, CD205, and CD206. As expected, most of these differentially expressed genes were either positively or negatively associated with senescence as defined by the expression of CDNK1A and enrichment scores of SenMayo senescence (data not shown). To further delve deeper into the underlying regulatory mechanisms of biological changes, we utilized the hdWGCNA analysis focused on alveolar macrophages and identified several expression networks associated with senescence signatures, including growth arrest and proliferation, TP53-mediated cell senescence-associated heat shock proteins, DNA damage and repair, and energy, protein, and lipid metabolism. These genes or biological processes are directly or indirectly involved in senescence processes of alveolar macrophages.\u003c/p\u003e \u003cp\u003eImportantly, we identified PPARγ as one of the top 10 pivotal transcription factors modulating senescent signatures of alveolar macrophages. Interestingly, PPARγ as a transcription factor has been shown to play a crucial role in various physiological processes, including adipogenesis\u003csup\u003e\u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e72\u003c/span\u003e,\u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e73\u003c/span\u003e\u003c/sup\u003e, glucose homeostasis\u003csup\u003e\u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e74\u003c/span\u003e\u003c/sup\u003e, inflammation\u003csup\u003e\u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e,\u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e76\u003c/span\u003e\u003c/sup\u003e, and metabolism and function of macrophages\u003csup\u003e\u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e77\u003c/span\u003e\u003c/sup\u003e. Studies have implicated that PPARγ regulates senescence in various cell types, including fibroblasts and endothelial cells\u003csup\u003e\u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e78\u003c/span\u003e,\u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e79\u003c/span\u003e\u003c/sup\u003e, and macrophages\u003csup\u003e\u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e77\u003c/span\u003e\u003c/sup\u003e. PPARγ activation can promote senescence by upregulating the expression of specific senescence-associated genes and inducing cell cycle arrest\u003csup\u003e\u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e79\u003c/span\u003e\u003c/sup\u003e. PPARγ activation can also have anti-senescent effects by reducing senescence and improve metabolic function\u003csup\u003e\u003cspan additionalcitationids=\"CR80 CR81\" citationid=\"CR79\" class=\"CitationRef\"\u003e79\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e82\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eTaken together, we present groundbreaking findings, demonstrating for the first time that monocytes/macrophages, particularly alveolar macrophages, are the predominant cell types in patients with ACO. Our research reveals a lower prevalence of senescence within alveolar macrophages in both ACO patients and those with severe asthma. Mechanistically, our in-depth exploration of differentially expressed genes within alveolar macrophages has identified PPARγ as a crucial regulatory factor driving senescence in ACO. Notably, there are several limitations, including small sample size, lack of mouse model to mimic patients with ACO, and complicated and context-dependent role of PPARγ in cellular senescence. However, our observations in this study provide a solid foundation for future investigations aimed at a comprehensive understanding of the mechanisms through which PPARγ regulates senescence. Furthermore, it opens doors for exploring PPARγ as a potential therapeutic target for interventions aimed at modulating senescence-associated processes in ACO -a condition marked by mixed features of both asthma and COPD, known for its increased severity.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eIL\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Interleukins\u003c/p\u003e\n\u003cp\u003eAM\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Alveolar Macrophages\u003c/p\u003e\n\u003cp\u003eIM\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;Interstitial Macrophages\u003c/p\u003e\n\u003cp\u003eSASP\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;Senescence Associated Secretory Phenotype\u003c/p\u003e\n\u003cp\u003eCOPD\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Chronic Obstructive Pulmonary Diseases\u003c/p\u003e\n\u003cp\u003eACO\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Asthma-COPD Overlap\u003c/p\u003e\n\u003cp\u003eCyTOF\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Cytometry by Time of Flight\u003c/p\u003e\n\u003cp\u003eBALF\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Bronchoalveolar Lavage Fluid\u003c/p\u003e\n\u003cp\u003eIMSA\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Immune Mechanisms Severe Asthma\u003c/p\u003e\n\u003cp\u003eUMAP\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;Uniform Manifold Approximation and Projection\u003c/p\u003e\n\u003cp\u003eDEGs\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Differential Expressed Genes\u003c/p\u003e\n\u003cp\u003eH\u0026amp;E\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Hematoxylin and Eosin\u003c/p\u003e\n\u003cp\u003eSCENIC\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Single-Cell rEgulatory Network Inference and Clustering\u003c/p\u003e\n\u003cp\u003ehdWGCNA\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;High-Dimensional Weighted Gene Co-epression Network Analysis\u003c/p\u003e\n\u003cp\u003ePAGA\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Partition-Based Graph Abstraction\u003c/p\u003e\n\u003cp\u003eGEO\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Gene Expression Omnibus\u003c/p\u003e\n\u003cp\u003eGSEA\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Geneset Enrichment Analysis\u003c/p\u003e\n\u003cp\u003eTRRUST: \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;Transcriptional Regulatory Relationships Unraveled by Sentence-based Text `\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;mining\u003c/p\u003e\n\u003cp\u003essGSEA\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;Single Sample GSEA\u003c/p\u003e\n\u003cp\u003eCRE\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Cockroach Extract\u003c/p\u003e\n\u003cp\u003eKEGG\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;Kyoto Encyclopedia of Genes and Genomes\u003c/p\u003e\n\u003cp\u003eGO\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;Gene Ontology\u003c/p\u003e\n\u003cp\u003ePPAR\u0026gamma; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; Peroxisome Proliferator Activated Receptor Gamma\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eDisclosure of potential conflict of interest:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors have declared that no conflict of interest exists\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eanalyzed data, performed experiments, and review the manuscript. R.W., M. W., P. G. designed and wrote the manuscript. Q. C., C. H., M. W., P. G., supervised the study. R.W., P. G. designed and wrote the manuscript. All authors read and approved the final version of the manuscript.\u003c/p\u003e\u003ch2\u003eACKNOWLEDGEMENTS\u003c/h2\u003e \u003cp\u003eThis work was supported by grants from the US National Institutes of Health (NIH) (1R01AI153331 and R01AI141642 to Gao).\u003c/p\u003e\u003ch2\u003eAvailability of Data\u003c/h2\u003e \u003cp\u003eAll bioinformatical data can be downloaded from figshare database (scRNA-Seq from fatal ACO and two health control: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.6084/m9.figshare\u003c/span\u003e\u003cspan address=\"10.6084/m9.figshare\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e. 14782377.v1.), GEO (Bulk RNA-Seq data from IMSA BALF: GSE136587) database and flow repository database (CyTOF data from IMSA BALF: FR-FCM-Z395).\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eLeung JM, Sin DD. Asthma-COPD overlap syndrome: pathogenesis, clinical features, and therapeutic targets. BMJ. 2017;358:j3772.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMekov E, Nu\u0026ntilde;ez A, Sin DD, et al. Update on Asthma-COPD Overlap (ACO): A Narrative Review. Int J Chron Obstruct Pulmon Dis. 2021;16:1783\u0026ndash;99.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBarczyk A, Maskey-Warzęchowska M, G\u0026oacute;rska K, et al. Asthma-COPD Overlap-A Discordance Between Patient Populations Defined by Different Diagnostic Criteria. J Allergy Clin Immunol Pract. 2019;7(7):2326\u0026ndash;36e5.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAlshabanat A, Zafari Z, Albanyan O, Dairi M, FitzGerald JM. Asthma and COPD Overlap Syndrome (ACOS): A Systematic Review and Meta Analysis. PLoS ONE. 2015;10(9):e0136065.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHikichi M, Hashimoto S, Gon Y. Asthma and COPD overlap pathophysiology of ACO. Allergol Int. 2018;67(2):179\u0026ndash;86.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGottlieb DJ, Sparrow D, O'Connor GT, Weiss ST. Skin test reactivity to common aeroallergens and decline of lung function. The Normative Aging Study. Am J Respir Crit Care Med. 1996;153(2):561\u0026ndash;6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCaillaud D, Chanez P, Escamilla R et al. Asthma-COPD overlap syndrome (ACOS) vs 'pure' COPD: a distinct phenotype? \u003cem\u003eAllergy\u003c/em\u003e 2017; 72(1): 137\u0026thinsp;\u0026ndash;\u0026thinsp;45.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eS\u0026aacute;nchez Castillo S, Smith L, D\u0026iacute;az Su\u0026aacute;rez A, L\u0026oacute;pez S\u0026aacute;nchez GF. Limitations in Activities of Daily Living among Older Adults with COPD, Asthma, or Asthma-COPD Overlap Residing in Spain. Int J Environ Res Public Health 2023; 20(4).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWoodruff PG, Modrek B, Choy DF, et al. T-helper type 2-driven inflammation defines major subphenotypes of asthma. Am J Respir Crit Care Med. 2009;180(5):388\u0026ndash;95.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eQin K, Xu B, Pang M, Wang H, Yu B. The functions of CD4 T-helper lymphocytes in chronic obstructive pulmonary disease. Acta Biochim Biophys Sin (Shanghai). 2022;54(2):173\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaradna A, Do DC, Kumar S, Fu QL, Gao P. Macrophage polarization and allergic asthma. Transl Res. 2018;191:1\u0026ndash;14.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCai Y, Sugimoto C, Arainga M, Alvarez X, Didier ES, Kuroda MJ. In vivo characterization of alveolar and interstitial lung macrophages in rhesus macaques: implications for understanding lung disease in humans. J Immunol. 2014;192(6):2821\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDraijer C, Peters-Golden M. Alveolar Macrophages in Allergic Asthma: the Forgotten Cell Awakes. Curr Allergy Asthma Rep. 2017;17(2):12.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBaptist AP, Busse PJ. Asthma Over the Age of 65: All's Well That Ends Well. J Allergy Clin Immunol Pract. 2018;6(3):764\u0026ndash;73.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ede Marco R, Pesce G, Marcon A, et al. The coexistence of asthma and chronic obstructive pulmonary disease (COPD): prevalence and risk factors in young, middle-aged and elderly people from the general population. PLoS ONE. 2013;8(5):e62985.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMendy A, Forno E, Niyonsenga T, Carnahan R, Gasana J. Prevalence and features of asthma-COPD overlap in the United States 2007\u0026ndash;2012. Clin Respir J. 2018;12(8):2369\u0026ndash;77.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDi Micco R, Krizhanovsky V, Baker D. d'Adda di Fagagna F. Cellular senescence in ageing: from mechanisms to therapeutic opportunities. Nat Rev Mol Cell Biol. 2021;22(2):75\u0026ndash;95.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSantoro A, Bientinesi E, Monti D. Immunosenescence and inflammaging in the aging process: age-related diseases or longevity? Ageing Res Rev. 2021;71:101422.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShvedova M, Samdavid Thanapaul RJR, Thompson EL, Niedernhofer LJ, Roh DS. Cellular Senescence in Aging, Tissue Repair, and Regeneration. Plast Reconstr Surg. 2022;150:4s\u0026ndash;11s.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMu\u0026ntilde;oz-Esp\u0026iacute;n D, Ca\u0026ntilde;amero M, Maraver A, et al. Programmed cell senescence during mammalian embryonic development. Cell. 2013;155(5):1104\u0026ndash;18.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang X, Qu M, Li J, Danielson P, Yang L, Zhou Q. Induction of Fibroblast Senescence During Mouse Corneal Wound Healing. Investig Ophthalmol Vis Sci. 2019;60(10):3669\u0026ndash;79.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCalcinotto A, Kohli J, Zagato E, Pellegrini L, Demaria M, Alimonti A. Cellular Senescence: Aging, Cancer, and Injury. Physiol Rev. 2019;99(2):1047\u0026ndash;78.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVicente R, Mausset-Bonnefont AL, Jorgensen C, Louis-Plence P, Brondello JM. Cellular senescence impact on immune cell fate and function. Aging Cell. 2016;15(3):400\u0026ndash;6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKang TW, Yevsa T, Woller N, et al. Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. Nature. 2011;479(7374):547\u0026ndash;51.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLecot P, Alimirah F, Desprez PY, Campisi J, Wiley C. Context-dependent effects of cellular senescence in cancer development. Br J Cancer. 2016;114(11):1180\u0026ndash;4.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYamamoto-Imoto H, Minami S, Shioda T, et al. Age-associated decline of MondoA drives cellular senescence through impaired autophagy and mitochondrial homeostasis. Cell Rep. 2022;38(9):110444.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMoiseeva V, Cisneros A, Sica V, et al. Senescence atlas reveals an aged-like inflamed niche that blunts muscle regeneration. Nature. 2023;613(7942):169\u0026ndash;78.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSoma T, Nagata M. Immunosenescence, Inflammaging, and Lung Senescence in Asthma in the Elderly. Biomolecules 2022; 12(10).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWan R, Srikaram P, Guntupalli V, Hu C, Chen Q, Gao P. Cellular senescence in asthma: from pathogenesis to therapeutic challenges. EBioMedicine. 2023;94:104717.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWoldhuis RR, Heijink IH, van den Berge M, et al. COPD-derived fibroblasts secrete higher levels of senescence-associated secretory phenotype proteins. Thorax. 2021;76(5):508\u0026ndash;11.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWoldhuis RR, de Vries M, Timens W, et al. Link between increased cellular senescence and extracellular matrix changes in COPD. Am J Physiol Lung Cell Mol Physiol. 2020;319(1):L48\u0026ndash;L60.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAdams TS, Schupp JC, Poli S et al. Single-cell RNA-seq reveals ectopic and aberrant lung-resident cell populations in idiopathic pulmonary fibrosis. \u003cem\u003eSci Adv\u003c/em\u003e 2020; 6(28): eaba1983.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchuliga M, Read J, Knight DA. Ageing mechanisms that contribute to tissue remodeling in lung disease. Ageing Res Rev. 2021;70:101405.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMoghbeli K, Valenzi E, Naramore R, et al. β-Agonist exposure preferentially impacts lung macrophage cyclic AMP-related gene expression in asthma and asthma COPD overlap syndrome. Am J Physiol Lung Cell Mol Physiol. 2021;321(5):L837\u0026ndash;l43.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCamiolo MJ, Zhou X, Oriss TB, et al. High-dimensional profiling clusters asthma severity by lymphoid and non-lymphoid status. Cell Rep. 2021;35(2):108974.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHao Y, Hao S, Andersen-Nissen E, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184(13):3573\u0026ndash;87e29.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKorsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods. 2019;16(12):1289\u0026ndash;96.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcInnes L, Healy J, Melville J, Umap. Uniform manifold approximation and projection for dimension reduction. \u003cem\u003earXiv preprint arXiv:180203426\u003c/em\u003e 2018.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTraag VA, Waltman L, van Eck NJ. From Louvain to Leiden: guaranteeing well-connected communities. Sci Rep. 2019;9(1):5233.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAran D, Looney AP, Liu L, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163\u0026ndash;72.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaul D, Kosinsky RL, Atkinson EJ, et al. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat Commun. 2022;13(1):4827.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWolf FA, Hamey FK, Plass M, et al. PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells. Genome Biol. 2019;20(1):59.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMorabito S, Reese F, Rahimzadeh N, Miyoshi E, Swarup V. High dimensional co-expression networks enable discovery of transcriptomic drivers in complex biological systems. bioRxiv 2022: 2022.09.22.509094.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLangfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLambert SA, Jolma A, Campitelli LF, et al. Hum Transcription Factors Cell. 2018;172(4):650\u0026ndash;65.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAibar S, Gonzalez-Blas CB, Moerman T, et al. SCENIC: single-cell regulatory network inference and clustering. Nat Methods. 2017;14(11):1083\u0026ndash;6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHan H, Cho JW, Lee S, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380\u0026ndash;d6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDavis S, Meltzer PS. GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics. 2007;23(14):1846\u0026ndash;7.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLeek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882\u0026ndash;3.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJassal B, Matthews L, Viteri G, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48(D1):D498\u0026ndash;d503.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHartigan JA, Wong MA, Algorithm AS. A k-means clustering algorithm. J royal Stat Soc Ser c (applied statistics). 1979;136(1):100\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572\u0026ndash;3.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRitchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNowicka M, Krieg C, Crowell HL, et al. CyTOF workflow: differential discovery in high-throughput high-dimensional cytometry datasets. F1000Research. 2019;6:748.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVan der Maaten L, Hinton G. Visualizing data using t-SNE. J Mach Learn Res 2008; 9(11).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eQuintelier K, Couckuyt A, Emmaneel A, Aerts J, Saeys Y, Van Gassen S. Analyzing high-dimensional cytometry data using FlowSOM. Nat Protoc. 2021;16(8):3775\u0026ndash;801.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWeber LM, Nowicka M, Soneson C, Robinson MD. diffcyt: Differential discovery in high-dimensional cytometry via high-resolution clustering. Commun Biol. 2019;2:183.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLopez R, Regier J, Cole MB, Jordan MI, Yosef N. Deep generative modeling for single-cell transcriptomics. Nat Methods. 2018;15(12):1053\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMartinez-Zamudio RI, Robinson L, Roux PF, Bischof O, SnapShot. Cell Senescence Pathways Cell. 2017;170(4):816\u0026ndash;e1.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChandra A, Lagnado AB, Farr JN, et al. Targeted clearance of p21- but not p16-positive senescent cells prevents radiation-induced osteoporosis and increased marrow adiposity. Aging Cell. 2022;21(5):e13602.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKohonen T. Essentials of the self-organizing map. Neural Netw. 2013;37:52\u0026ndash;65.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVan Gassen S, Callebaut B, Van Helden MJ, et al. FlowSOM: Using self-organizing maps for visualization and interpretation of cytometry data. Cytometry A. 2015;87(7):636\u0026ndash;45.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDey S, Eapen MS, Chia C, Gaikwad AV, Wark PAB, Sohal SS. Pathogenesis, clinical features of asthma COPD overlap, and therapeutic modalities. Am J Physiol Lung Cell Mol Physiol. 2022;322(1):L64\u0026ndash;L83.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDey S, Lu W, Haug G, et al. Airway inflammatory changes in the lungs of patients with asthma-COPD overlap (ACO): a bronchoscopy endobronchial biopsy study. Respir Res. 2023;24(1):221.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJeffery PK, Laitinen A, Venge P. Biopsy markers of airway inflammation and remodelling. \u003cem\u003eRespiratory medicine\u003c/em\u003e 2000; 94 Suppl F: S9-15.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEapen MS, Myers S, Walters EH, Sohal SS. Airway inflammation in chronic obstructive pulmonary disease (COPD): a true paradox. Expert Rev Respir Med. 2017;11(10):827\u0026ndash;39.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEapen MS, Hansbro PM, McAlinden K, et al. Abnormal M1/M2 macrophage phenotype profiles in the small airway wall and lumen in smokers and chronic obstructive pulmonary disease (COPD). Sci Rep. 2017;7(1):13392.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGorgoulis V, Adams PD, Alimonti A, et al. Cellular Senescence: Defining a Path Forward. Cell. 2019;179(4):813\u0026ndash;27.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBarnes PJ, Baker J, Donnelly LE. Cellular Senescence as a Mechanism and Target in Chronic Lung Diseases. Am J Respir Crit Care Med. 2019;200(5):556\u0026ndash;64.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMitsi E, Kamng'ona R, Rylance J, et al. Human alveolar macrophages predominately express combined classical M1 and M2 surface markers in steady state. Respir Res. 2018;19(1):66.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePark HS, Ju UI, Park JW, et al. PPARgamma neddylation essential for adipogenesis is a potential target for treating obesity. Cell Death Differ. 2016;23(8):1296\u0026ndash;311.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen Z, Ghavimi SAA, Wu M et al. PPARgamma agonist treatment reduces fibroadipose tissue in secondary lymphedema by exhausting fibroadipogenic PDGFRalpha\u0026thinsp;+\u0026thinsp;mesenchymal cells. JCI Insight 2023; 8(24).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMosana M, Ayeleso A, Nyakudya T, Erlwanger K, Mukwevho E. Potential Protective Effects of Neonatal Supplementation with Oleanolic Acid on Peroxisome Proliferator-Activated Receptor Gamma (PPARγ)-Ligand Dependent Regulation of Glucose Homeostasis in High Fructose-Fed Rats. Nat Prod Commun 2020; 15(3).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStark JM, Coquet JM, Tibbitt CA. The Role of PPAR-gamma in Allergic Disease. Curr Allergy Asthma Rep. 2021;21(11):45.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWei R, Han X, Li M, et al. The nuclear cytokine IL-37a controls lethal cytokine storms primarily via IL-1R8-independent transcriptional upregulation of PPARgamma. Cell Mol Immunol. 2023;20(12):1428\u0026ndash;44.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYu L, Gao Y, Aaron N, Qiang L. A glimpse of the connection between PPARgamma and macrophage. Front Pharmacol. 2023;14:1254317.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu X, Gu Y, Kumar S, et al. Oxylipin-PPARgamma-initiated adipocyte senescence propagates secondary senescence in the bone marrow. Cell Metab. 2023;35(4):667\u0026ndash;84. e6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eToobian D, Ghosh P, Katkar GD. Parsing the Role of PPARs in Macrophage Processes. Front Immunol. 2021;12:783780.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang HC, Deleuze S, Zuo Y, Potthoff SA, Ma LJ, Fogo AB. The PPARgamma agonist pioglitazone ameliorates aging-related progressive renal injury. J Am Soc Nephrol. 2009;20(11):2380\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVillapol S. Roles of Peroxisome Proliferator-Activated Receptor Gamma on Brain and Peripheral Inflammation. Cell Mol Neurobiol. 2018;38(1):121\u0026ndash;32.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDaniel B, Nagy G, Czimmerer Z et al. The Nuclear Receptor PPARgamma Controls Progressive Macrophage Polarization as a Ligand-Insensitive Epigenomic Ratchet of Transcriptional Memory. Immunity 2018; 49(4): 615\u0026thinsp;\u0026ndash;\u0026thinsp;26 e6.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"respiratory-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"rere","sideBox":"Learn more about [Respiratory Research](http://respiratory-research.biomedcentral.com/)","snPcode":"12931","submissionUrl":"https://submission.nature.com/new-submission/12931/3","title":"Respiratory Research","twitterHandle":"@RespiratoryBMC","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"ACO, Asthma, COPD, Macrophages, Senescence, PPARγ","lastPublishedDoi":"10.21203/rs.3.rs-4009724/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4009724/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eAsthma-chronic obstructive pulmonary disease (COPD) overlap (ACO) represents a complex condition characterized by shared clinical and pathophysiological features of asthma and COPD in older individuals. However, the pathophysiology of ACO remains unexplored. We aimed to identify the major inflammatory cells in ACO, examine senescence within these cells, and elucidate the genes responsible for regulating senescence. Bioinformatic analyses were performed to investigate major cell types and cellular senescence signatures in a public single-cell RNA sequencing (scRNA-Seq) dataset derived from the lung tissues of patients with ACO. Similar analyses were carried out in an independent cohort study Immune Mechanisms Severe Asthma (IMSA), which included bulk RNA-Seq and CyTOF data from bronchoalveolar lavage fluid (BALF) samples. The analysis of the scRNA-Seq data revealed that monocytes/ macrophages were the predominant cell type in the lung tissues of ACO patients, constituting more than 50% of the cells analyzed. Lung monocytes/macrophages from patients with ACO exhibited a lower prevalence of senescence as defined by lower enrichment scores of SenMayo and expression levels of cellular senescence markers. Intriguingly, analysis of the IMSA dataset showed similar results in patients with severe asthma. They also exhibited a lower prevalence of senescence, particularly in airway CD206\u0026thinsp;+\u0026thinsp;macrophages, along with increased cytokine expression (e.g., \u003cem\u003eIL-4, IL-13\u003c/em\u003e, and \u003cem\u003eIL-22\u003c/em\u003e). Further exploration identified alveolar macrophages as a major subtype of monocytes/macrophages driving cellular senescence in ACO. Differentially expressed genes related to oxidation-reduction, cytokines, and growth factors were implicated in regulating senescence in alveolar macrophages. PPARγ (Peroxisome Proliferator-Activated Receptor Gamma) emerged as one of the predominant regulators modulating the senescent signature of alveolar macrophages in ACO. Collectively, the findings suggest that senescence in macrophages, particularly alveolar macrophages, plays a crucial role in the pathophysiology of ACO. Furthermore, PPARγ may represent a potential therapeutic target for interventions aimed at modulating senescence-associated processes in ACO.\u003c/p\u003e","manuscriptTitle":"PPARγ Attenuates Cellular Senescence of Alveolar Macrophages in Asthma- COPD Overlap","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-03-05 17:08:38","doi":"10.21203/rs.3.rs-4009724/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-03-19T19:13:38+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-03-17T11:39:59+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"872a63e1-9146-492e-a106-998939c00d0f","date":"2024-03-04T17:15:08+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-03-04T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-03-04T11:55:44+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-03-04T07:23:36+00:00","index":"","fulltext":""},{"type":"submitted","content":"Respiratory Research","date":"2024-03-03T21:36:07+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"respiratory-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"rere","sideBox":"Learn more about [Respiratory Research](http://respiratory-research.biomedcentral.com/)","snPcode":"12931","submissionUrl":"https://submission.nature.com/new-submission/12931/3","title":"Respiratory Research","twitterHandle":"@RespiratoryBMC","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"6eb99414-9013-48b6-8219-c551af708ca4","owner":[],"postedDate":"March 5th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-05-01T23:00:23+00:00","versionOfRecord":{"articleIdentity":"rs-4009724","link":"https://doi.org/10.1186/s12931-024-02790-6","journal":{"identity":"respiratory-research","isVorOnly":false,"title":"Respiratory Research"},"publishedOn":"2024-04-20 23:00:23","publishedOnDateReadable":"April 20th, 2024"},"versionCreatedAt":"2024-03-05 17:08:38","video":"","vorDoi":"10.1186/s12931-024-02790-6","vorDoiUrl":"https://doi.org/10.1186/s12931-024-02790-6","workflowStages":[]},"version":"v1","identity":"rs-4009724","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4009724","identity":"rs-4009724","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}