Research Paper Volume 15, Issue 7 pp 2395—2417
p21 facilitates chronic lung inflammation via epithelial and endothelial cells
- 1 Department of Molecular Cell Biology, The Weizmann Institute of Science, Rehovot 7610001, Israel
- 2 Department of Life Sciences Core Facilities, The Weizmann Institute of Science, Rehovot 7610001, Israel
- 3 Department of Veterinary Resources, The Weizmann Institute of Science, Rehovot 7610001, Israel
- 4 Department of Immunology, The Weizmann Institute of Science, Rehovot 7610001, Israel
Received: November 8, 2022 Accepted: March 17, 2023 Published: March 30, 2023
https://doi.org/10.18632/aging.204622How to Cite
Copyright: © 2023 Levi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Abstract
Cellular senescence is a stable state of cell cycle arrest that regulates tissue integrity and protects the organism from tumorigenesis. However, the accumulation of senescent cells during aging contributes to age-related pathologies. One such pathology is chronic lung inflammation. p21 (CDKN1A) regulates cellular senescence via inhibition of cyclin-dependent kinases (CDKs). However, its role in chronic lung inflammation and functional impact on chronic lung disease, where senescent cells accumulate, is less understood. To elucidate the role of p21 in chronic lung inflammation, we subjected p21 knockout (p21-/-) mice to repetitive inhalations of lipopolysaccharide (LPS), an exposure that leads to chronic bronchitis and accumulation of senescent cells. p21 knockout led to a reduced presence of senescent cells, alleviated the pathological manifestations of chronic lung inflammation, and improved the fitness of the mice. The expression profiling of the lung cells revealed that resident epithelial and endothelial cells, but not immune cells, play a significant role in mediating the p21-dependent inflammatory response following chronic LPS exposure. Our results implicate p21 as a critical regulator of chronic bronchitis and a driver of chronic airway inflammation and lung destruction.
Introduction
Chronic lung inflammation is the major component of chronic obstructive pulmonary disease (COPD), the third most common cause of death in developed countries, that has become a global epidemic and is increasing worldwide [1, 2]. This disease is associated with chronic inflammation of the peripheral airways and lung parenchyma (chronic bronchitis), the subsequent destruction of the lung alveoli (emphysema), and the accumulation of senescent cells in different lung compartments [1, 3–6]. The presence of senescent cells is associated with chronic inflammation, which is characterized by the recruitment of innate and adaptive immune cells, [7–10]. The peripheral immune cells are attracted to the lung tissue by secretion of pro-inflammatory mediators, including cytokines, chemokines and growth factors, by epithelial [11], endothelial [12] and fibroblast [13] cells in the lungs. Epithelial cells are activated in response to chronic airway irritation, caused by smoke and other pollutants, to produce inflammatory mediators [11, 14–16]. Endothelial cells are also activated following chronic damage to the lung [12, 17] and express cytokines and chemokines, via NF-κB transcriptional activation [18]. Activated endothelial cells also express adhesion molecules which increase vascular barrier permeability, such as intercellular adhesion molecule-1 (ICAM-1), and vascular cell adhesion molecule-1 (VCAM-1), which promotes adherence of immune cells [19–21]. However, the contribution of specific molecular and cellular mechanisms in epithelial and endothelial cells to the chronic inflammation requires further investigation.
Aging is associated with an increased risk of chronic inflammatory lung diseases [22–24]. One of the hallmarks of aging is cellular senescence, and indeed, senescent cells of epithelial, endothelial, fibroblast and immune origin accumulate in lungs with chronic inflammation [25, 26]. Cellular senescence refers to a stable state of cell cycle arrest, accompanied by a profound secretory phenotype [27–29]. On the short term, the presence of senescence cells plays a role in tumor suppression, wound healing, as well as in embryonic development [30–36]. Conversely, accumulation of senescent cells in tissues, a process that occurs during aging, stimulates tumorigenesis and the development of age-related diseases by promoting chronic inflammation [37–39]. This process is mediated by expression and secretion of proteins, members of the senescence associated secretory phenotype (SASP), which leads to chronic activation of immune responses in the lungs [40–45]. However, the relative contribution of senescence mechanisms in different cell types to the chronic inflammation is not sufficiently understood.
The senescence growth arrest is established and maintained by the p53-p21 and p16-Rb pathways [27, 46–48]. One of the cell-cycle inhibitors that is often expressed by senescent cells is the cyclin-dependent kinase inhibitor (CDKI) p21 (also termed CDKN1A). p21 can inhibit the cyclin-dependent kinases (CDKs) CDK2 and CDK4, leading to suppression of retinoblastoma protein (Rb) phosphorylation, and to a subsequent downregulation of genes necessary for cell-cycle progression [42, 49–51]. In addition to its role in senescence, p21 is involved in lung pathologies and p21 ablation attenuates the chronic inflammatory responses and emphysema [52–55]. The contribution of p21 to chronic lung inflammation may rely on its ability to regulate multiple phenotypes of senescent cells. Indeed, one of the characteristics of senescent cells that is necessary for their persistence, resistance to cell death [56, 57] is mediated by p21 and BCL-2 family members [58–60].
While p21 is a key regulator of cellular senescence, its role in specific cell populations in the lung during chronic inflammation and functional output of chronic lung disease requires further investigation. Of note, p53, the upstream mediator of p21 expression, facilitates chronic lung inflammation in bronchial epithelial cells by promoting cellular senescence [26]. We, therefore, utilized a LPS inhalation-induced chronic bronchitis procedure [26] to study the effects of repetitive LPS exposure on p21 knockout (p21-/-) mice. In this model, p21 knockout leads to a decrease in the inflammatory response, reduced accumulation of senescent cells, improved fitness, and alleviation of the pathological manifestations of chronic lung inflammation. Furthermore, we aimed to examine the specific contribution of the epithelial, endothelial and immune compartments to chronic bronchitis pathology. We revealed that resident epithelial and endothelial cells in the lungs, but not immune cells, play a significant role in mediating the p21-dependent inflammatory response following chronic LPS exposure. Overall, these findings suggest that p21 is a crucial regulator of chronic airway inflammation and lung destruction.
Results
Chronic LPS inhalation induces p21-dependent cellular senescence in bronchial epithelia
Repetitive lung injury causes p53-dependent accumulation of senescent cells in the lungs, which contributes to chronic lung inflammation and promotes chronic bronchitis [26]. p21, a downstream mediator of p53, regulates cell cycle arrest and the viability of senescent cells [60]. Here we set to understand the effect of p21 during chronic lung inflammation. To this end, we exposed wild type (WT) and p21 knockout (p21-/-) 8-week-old female mice to chronic LPS-inhalation (0.5 mg/ml) regimen, which causes chronic lung inflammation [26], or to PBS-inhalation regimen as a control, 3 times a week for 10 weeks. In order to understand the effect of p21 knockout on accumulation of senescent cells following exposure to repetitive LPS inhalations, we analyzed the mice lungs for cellular senescence markers. We first evaluated the protein expression levels of the senescence markers p15, p16, p21 and p53 in the lungs of WT and p21-/- mice by immunoblot analysis and the mRNA expression levels of p15, p16 and p21. This analysis revealed increased expression of the senescence markers in the lungs of LPS exposed WT mice, compared to PBS exposed mice, while p21 knockout abolished this effect (Figure 1A, 1B and Supplementary Figure 1A, representative immunoblots).
Figure 1. Accumulation of senescent cells is decreased in the lungs of p21-/- mice. WT and p21-/- mice were exposed to either PBS or aerosolized LPS (0.5 mg/ml), 3 times a week for 10 weeks. At 48 hours following the last LPS exposure, the lungs were harvested and frozen. Alternatively, lungs were harvested, fixed, and analyzed for markers of senescence. (A) Representative immunoblots for senescence-associated proteins p15, p16, p21 and p53 in the mice lungs. (B) mRNA expression levels of senescence markers p15, p16 and p21 in the mice lungs. (C, D) Immunohistochemistry (IHC) of lung sections for cleaved caspase 3, (C) and p21, (D). Scale bar represents 200μm. (E) SA-β-Gal staining of lung sections. Scale bar represents 200μm. (F) Quantification of the number of p21 positive cells per bronchial area, of the lung sections presented in (D). (G) Quantification of SA-β-gal (%) per bronchial area, of the lung sections presented in (E). Data information: Data were analyzed using one-way ANOVA, *P<0.05. **<0.005. ***P<0.0005 (B), and by Student’s t-test, *p<0.05, **p<0.01, and ***p<0.005 (F, G). Data represent mean ±SEM (A, n=3; B, n=3-10; C–G, n=3-6 independent repeats).
Senescent cells exacerbate chronic bronchitis pathology and their elimination limits the disease progression [26]. Furthermore, p21 silencing results in elimination of senescent cells via apoptosis [60]. To determine the effect of p21 knockout on accumulation of senescent bronchial epithelial cells we evaluated lung sections stained for the apoptotic marker cleaved caspase-3, a cellular senescence marker p21, and for SA-β-gal activity (Figure 1C–1E). Neither WT nor p21-/- mice presented caspase-3 cleavage following chronic LPS exposure (Figure 1C and Supplementary Figure 1B - positive control), suggesting that the bronchial epithelial cells do not undergo significant apoptosis in our experimental end-point. The immunoblot analysis of the cleavage caspase-3 did not reveal a significant difference in expression of this marker between the genotypes (Supplementary Figure 1C). Alternatively, immunohistochemical analysis of p21 expression revealed a marked elevation in the amount of p21 expressing cells in LPS exposed WT mice, compared to the PBS exposed mice (Figure 1D, 1F). Furthermore, the levels of SA-β-gal activity in the lungs of LPS exposed WT mice were increased relative to lungs of PBS exposed mice (Figure 1E), while p21 knockout abolished this effect (Figure 1E, 1G). These results suggest that following chronic LPS exposure lung cells undergo senescence via a p21-dependent pathway.
p21-/- diminishes the pathological manifestations of chronic bronchitis
In order to understand the effect of p21 knockout on the progression of chronic bronchitis following chronic LPS exposure, we evaluated the pathology parameters of the mice in our experimental groups. Hematoxylin and eosin (H&E) staining revealed that chronic LPS exposure in WT mice resulted in airway inflammation, as reflected by infiltration of immune cells around the bronchi (Figure 2A), and in a significant increase in bronchial wall thickness, compared to PBS treated mice (Figure 2B). Strikingly, p21-/- mice had not shown accumulation of immune cells around the bronchi following chronic LPS exposure (Figure 2A) and their bronchial wall thickness was not significantly different compared to PBS exposed p21-/- mice (Figure 2B). Therefore, p21 knockout alleviates the histological signs of chronic lung inflammation induced by chronic LPS exposure.
Figure 2. p21-/- improves the pathophysiology triggered by chronic LPS exposure. WT and p21-/- mice were exposed to either PBS or aerosolized LPS (0.5 mg/ml), 3 times a week for 10 weeks. At 48 hours following the last LPS exposure the lungs were harvested, fixed, and stained by H&E. (A) Representative images of lung sections featuring bronchi and alveoli of WT and p21-/- mice. Scale bar represents: 200μm. (B) Airway wall thickness of the lung described in (A). (C) mRNA expression levels of chronic bronchitis associated component Mmp12 in the lungs described in (A). (D) At week 9, mice fitness was assessed by a treadmill assay and was presented as maximal actual distance. (E) Delta maximal actual distance was calculated from the baseline point until week 9 of the mice described in (D). Data information: Data were analyzed using one-way ANOVA, *P<0.05. **<0.005. ***P<0.0005 (B, C), and by Mann Whitney u test (D, E), *p<0.05. Data represent mean ±SEM (A, B, n=5; C, n=3-8; D, E, n=4-5 independent repeats).
Matrix Metalloproteases (MMPs) are a family of proteinases that have the ability to degrade different extracellular matrix components [61]. MMP12, a member of the MMPs family, mediates alveolar destruction during emphysema and it promotes the pathogenesis of chronic bronchitis [62]. Therefore, we set to evaluate MMP12 expression following chronic LPS exposure in WT and p21-/- mice. As expected, WT mice that were exposed to LPS presented increased mRNA expression levels of Mmp12 compared to the PBS exposed mice. However, p21 knockout decreased this effect, and p21-/- mice presented significantly lower levels of Mmp12 then observed in the LPS exposed WT mice (Figure 2C). We then evaluated the accumulation of fibrosis in our model. Sirius Red staining, which stains ECM (Supplementary Figure 2A, 2B) and protein expression of COL1 (Supplementary Figure 2C) revealed no significant changes between PBS and LPS treated WT mice. Therefore, while Mmp12 levels are affected by LPS exposure, no significant accumulation of fibrotic tissue was observed following this treatment.
Since chronic bronchitis causes damages to the lung in a manner that affects fitness [63], we set to examine whether chronic LPS exposure has an effect on the mice fitness. We subjected the mice to a treadmill assay, a behavioral test that measures the ability of mice to sustain running. The assay was performed at 5 time points, specifically, 3 days before the first LPS exposure (hereafter, baseline), 3, 5, 7 and 9 weeks in-parallel with LPS exposure. We observed that following 9 weeks of chronic LPS exposure, the maximal actual distance mice run until exhaustion was significantly shorter among LPS exposed WT mice compared to PBS exposed WT mice (Figure 2D). In contrast, the maximal actual distance that p21-/- mice run until exhaustion did not differ between PBS and LPS exposed mice. Moreover, the delta maximal actual distance that was calculated from the baseline point was significantly higher in the PBS exposed WT mice compared to the LPS exposed WT mice (Figure 2E), but did not differ between p21-/- mice exposed to either PBS or LPS. These results indicate that the LPS exposure regimen reduces the fitness of WT, but not p21-/- mice exposed to LPS at the end of the exposure period. Altogether, p21 knockout limits LPS-induced chronic bronchitis as manifested both at tissue pathology and organismal fitness.
p21 promotes chronic inflammatory responses in the lung airways
A central hallmark of chronic lung inflammation is the recruitment and accumulation of immune cells within the lung parenchyma and the bronchoalveolar lavage (BAL) fluid [9, 64]. Therefore, to determine whether p21-/- affects the chronic inflammatory response in the lungs, we examined the cellular immune response in both the lung parenchyma and the BAL fluid of mice subjected to the chronic LPS regimen. The cellular immune response and infiltration of both innate and adaptive immune cells in the mice lung parenchyma following chronic LPS exposure was assessed by flow cytometry. This analysis revealed a significantly higher influx of CD45+ immune cells (Figure 3A) in the lungs of LPS exposed WT mice compared to the PBS exposed WT mice. By analysing the presence of specific components of the immune system, we revealed a significant accumulation of neutrophils, interstitial macrophages, CD3+, CD4+, and CD8+ T cells (Figure 3B–3F) following chronic LPS exposure. Remarkably, p21 knockout resulted in a significantly lower accumulation of all these immune cell subsets (Figure 3A–3F). The inflammatory response in the lung is accompanied by an elevation in inflammatory secreted components, such as cytokines and chemokines [65–67]. To examine whether p21 knockout altered the inflammatory response in LPS-exposed mice, we analyzed the mRNA expression levels of pro-inflammatory SASP components. We observed a significant increase in the expression of Il-1β, Ccl5, Cxcl1, Cxcl2, Cxcl9, Cxcl10, Cxcl11 and Tnf-α in LPS exposed WT mice compared with PBS exposed mice, but not in p21-/- mice (Figure 3G). In both WT and p21-/- mice, however, the LPS-triggered increase in Ifn-γ, Il-6, Kc and Cxcl5 remained the same.
Figure 3. p21-/- decreases chronic inflammatory responses and iBALT formation caused by exposure to aerosolized LPS. WT and p21-/- mice were exposed to either PBS or aerosolized LPS (0.5 mg/ml), 3 times a week for 10 weeks. (A–F) At 48 hours following the last LPS exposure, whole lungs were dissociated into single cell suspensions and analyzed by flow cytometry to determine: (A) numbers of immune cells (CD45+), (B) numbers of neutrophils (CD45+/Ly6G+/CD11b+), (C) numbers of interstitial macrophages (CD45+/CD11c+/SiglecF-/CD11b+/CD24+), (D) numbers of CD3+ T cells (CD45+/CD3+), (E) numbers of CD4+ T cells (CD45+/CD3+/CD4+), and (F) numbers of CD8+ T cells (CD45+/CD3+/CD8+). (G) mRNA expression levels of the indicated SASP factors in the mice lungs. (H) Representative images of lungs stained for CD3+ (red) and B220 (green) depict accumulation of iBALTs in LPS exposed mice. Scale bar represents 200μm. (I, J) Numbers (I) and sizes (J) of iBALTs in the lungs of mice exposed to aerosolized LPS. Data information: Data were analyzed using one-way ANOVA, *P<0.05. **<0.005. ***P<0.0005 (A–G), and by Student’s t-test, *p<0.05, **p<0.01, and ***p<0.005 (I, J). Data represent mean ±SEM (A–F, n=9-12; G, n=3-6; H–J; n=4-6 independent repeats).
To understand the inflammatory response in the lung parenchyma, we explored in-situ the location of the adaptive immune populations of T and B cells, which aggregate in the LPS-exposed mice lungs. These aggregations present histological features of inducible bronchus-associated lymphoid tissue (iBALT). Using immunofluorescence co-staining we detected the presence of iBALTs in the chronic LPS exposed mice lungs (Figure 3H). Notably, in p21-/- LPS exposed mice, the iBALTs number was 4-fold lower than in WT LPS exposed mice (Figure 3I). In addition, the iBALTs size was smaller in the p21-/- mice compared to the WT mice (Figure 3J, p=0.06). Since p21-/- mice presented a low number of iBALTs in comparison to WT mice, fewer iBALTs were measured, thus negatively affecting the ability to identify the difference between the genotypes. However, it is evident that p21-/- leads to a reduction in the number and the size of iBALTs, indicating lower level of immune cell accumulation in these mice.
Then, we evaluated the cellular immune response in the mice BAL fluid by flow cytometry. As in the lung parenchyma, we observed a significantly higher influx of neutrophils, B cells, CD3+, CD4+ and CD8+ T cells (Figure 4A, 4C–4F), accompanied by an elevated influx of alveolar macrophages (Figure 4B) in WT LPS exposed mice compared to the PBS exposed WT mice. In contrast, a 2-fold decrease on average in the influx of these immune cell populations was observed in p21-/- LPS exposed mice (Figure 4A–4F). It, therefore, appears that p21 knockout alleviates the accumulation of immune cells in lungs, including specifically in iBALTs, following chronic LPS exposure.
Figure 4. p21-/- protects lung airways from chronic inflammatory responses triggered by chronic LPS exposure. WT and p21-/- mice were exposed to either PBS or aerosolized LPS (0.5 mg/ml), 3 times a week for 10 weeks. (A–F) At 48 hours following the last LPS exposure BAL fluid was collected and flow cytometry was used to determine: (A) numbers of neutrophils (CD45+/Ly6G+/CD11b+), (B) numbers of alveolar macrophages (CD45+/CD11c+/SiglecF+), (C) numbers of B cells (CD45+/B220+), (D) numbers of CD3+ T cells (CD45+/CD3+), (E) numbers of CD4+ T cells (CD45+/CD3+/CD4+), and (F) numbers of CD8+ T cells (CD45+/CD3+/CD8+). Data information: Data were analyzed using one-way ANOVA. ***P<0.0005. Data represent mean ±SEM (A–F, n=8-11 independent repeats).
Epithelial and endothelial cell populations mediate p21-dependent inflammatory responses following chronic LPS inhalation
Tissue physiology and pathophysiology are dictated by the presence and interaction of different cell types. Of note, LPS treatment results in profound changes in inflammation and immune cells recruitment to the lungs. To better understand the role of p21 in different cell types in the inflammatory response during chronic lung inflammation, we performed a transcriptome analysis for the lung’s three most abundant cell populations, namely, epithelial, endothelial and immune cells. Due to the absence of fibrosis in this experimental system (Supplementary Figure 2A–2C), the amount of fibroblasts in the lung is relatively small. Therefore, this cell population was not evaluated. The cells were FACS-sorted from WT and p21-/- mice lungs and subjected to mRNA sequencing. We analyzed the data for each cell type separately and identified the differentially-expressed genes (DEGs) affected by LPS treatment and/or p21 knockout. We performed three distinct comparisons: LPS WT vs. PBS WT, LPS p21-/- vs. PBS p21-/- and LPS p21-/- vs. LPS WT. Out of these comparisons, we selected the DEGs that passed the threshold of |log2FoldChange|≥1, average mean normalized reads |≥5 and adjusted p value≤0.05 in at least one of the comparisons. 13754, 12947 and 12806 genes were identified in epithelial, endothelial and immune cells, respectively. Out of the identified genes, 797, 228 and 238 DEGs were identified in epithelial, endothelial and immune cells, respectively. We performed partitioning clustering for each cell type and presented the data as a heatmap (Figure 5A). For each cell population we examined the clusters and observed that following chronic LPS exposure WT mice presented clusters of DEGs that were upregulated in epithelial (5 clusters), endothelial (3 clusters) and immune cell (2 clusters) populations, compared to PBS exposed mice. However, p21 knockout resulted in a general downregulation in the expression of DEGs in the p21-affected clusters. Then, we performed a functional analysis of the genes of the p21-affected clusters to uncover the enriched Gene Ontology (GO) terms of each specific cell population (Figure 5B). The “inflammatory response” term was enriched in all cell populations, reflecting the severe response to chronic LPS exposure. Surprisingly, however, the epithelial and endothelial cell populations showed enrichment in 8 and 11 GO terms, respectively, related to inflammatory pathways, while the immune cells were enriched only in 3 such GO terms (Figure 5B).
Figure 5. Epithelial and endothelial cell populations mediate p21-dependent inflammatory responses following chronic LPS inhalation. WT and p21-/- mice were exposed to either PBS or aerosolized LPS (0.5 mg/ml), 3 times a week for 10 weeks. At 48 hours following the last LPS exposure, whole lungs were dissected and dissociated into single cell suspensions and analyzed by flow cytometry and sorting for subsequent RNA sequencing. Data was analyzed as follows: (A) Hierarchical clustering heatmap of differentially expressed genes (DEGs) in epithelial, endothelial and immune cells populations in the mice lungs. One sample of the immune cells population, in the PBS p21-/- mice group, was excluded from the analysis. (B) Functional analysis of the DEGs in epithelial, endothelial and immune cells populations in the mice lungs. Abbreviations: antigen processing and presentation - antigen presentation, cytokine-mediated signaling pathway - cytokine-mediated signaling, external side of plasma membrane - external plasma membrane, glomerulus vasculature morphogenesis - glomerulus vasculature morphogenesis, positive regulation of immune response - immune response regulation, positive regulation of MAPK cascade - regulation of MAPK cascade, regulation of cytokine production - cytokine production, regulation of tumor necrosis factor superfamily cytokine production - TNF cytokine production. (C) Heatmap based on functional analysis of the upregulated DEGs in epithelial, endothelial and immune cells populations in the mice lungs. (D) Number of upregulated LPS responsive DEGs in epithelial, endothelial and immune cells populations in the mice lungs. (E) Number of upregulated LPS responsive DEGs in p21-/- affected clusters in epithelial, endothelial and immune cells populations in the mice lungs. Data information: Data were analyzed using Chi-squared test. ***P<0.0005. (A–E), n=4-5 independent repeats.
To reveal the possible functional impact of the upregulated DEGs, we performed a functional analysis to compare the enrichment of the most significant GO terms between epithelial, endothelial and immune cells populations. We observed that both epithelial and endothelial cell populations presented enrichment of inflammatory response related terms such as “superoxide anion generation”, “antigen processing and presentation of exogenous peptide antigen”, “positive regulation of immune response” and “leukocyte proliferation”, which, surprisingly, were not enriched in the immune cells population (Figure 5C). Interestingly, we observed that the terms “axoneme assembly” and “epithelial cilium movement”, which are functionally related to epithelia, were enriched only in the epithelial cell population.
To further examine which genes are affected by p21 following LPS treatment in each cell population, we compared the upregulated DEGs following LPS treatment in both WT (LPS WT vs. PBS WT) and p21-/- (LPS p21-/- vs. PBS p21-/-) mice. First, a KEGG pathway analysis in epithelial, endothelial and immune cells revealed no differences in gene expression when comparing between PBS WT and PBS p21-/- mice. In contrast, we observed that chronic LPS exposure resulted in a significant upregulation of DEGs in WT mice compared with p21-/- mice, in the epithelial (p<0.001) and endothelial (p<0.001) cell populations, but not in the immune cells population (p=0.543) (Figure 5D). Furthermore, a KEGG pathway analysis revealed that chronic LPS exposure resulted in changes in the inflammatory response, the p53 signaling pathway and cellular senescence pathways in WT mice compared with p21-/- mice in endothelial cells (Supplementary Figure 3A).
Then, in order to determine whether the differences observed in the epithelial and endothelial cell populations following chronic LPS exposure are p21-dependent, we compared the amount of upregulated DEGs following LPS treatment specifically in the p21 affected clusters in both WT (LPS WT vs. PBS WT) and p21-/- (LPS p21-/- vs. PBS p21-/-) mice. We observed that chronic LPS exposure resulted in a significant upregulation of DEGs in the p21 affected clusters in WT mice compared with p21-/- mice, in the epithelial (p<0.001) and endothelial (p<0.001) cell populations, but not in the immune cells population (p=0.644) (Figure 5E). Therefore, the effect of p21 knockout is prominent in epithelial and endothelial but not in immune cells. Overall, this analysis suggests that resident epithelial and endothelial cells in the lungs, but not immune cells, play a significant role in mediating the p21-dependent regulation of inflammatory response following chronic LPS exposure.
Discussion
The long-term persistence of senescent cells is deleterious and can promote the development of age-related pathologies [28, 29, 68]. Here, we show that senescent cells accumulate during chronic bronchitis in mice, and p21 knockout decreases their presence. We observed a significant decrease in the main characteristics of chronic bronchitis in p21 knockout, compared to WT mice. The main characteristics affected were the influx of immune cells into the BAL fluid and lung parenchyma, inflammatory cytokines in the broncho-alveolar spaces, and the pathological manifestations of chronic bronchitis, including iBALTs. These results show that p21 knockout reduces the long-term immunological response to LPS. However, the signs of inflammation, including CD3+ and CD4+ T cells in the BAL fluid and the neutrophils and interstitial macrophages in the lung parenchyma, were still present in LPS treated p21 knockout mice. Furthermore, our study illustrates that resident epithelial and endothelial cells in the lungs, but not immune cells, play a significant role in mediating the p21-dependent inflammatory response during chronic bronchitis. These findings suggest that p21 is a crucial regulator of chronic airway inflammation and lung destruction and a significant driver of chronic bronchitis. These effects are mediated through epithelial and endothelial cells.
The increase in p21 levels maintains the viability of senescent cells following the stress they might undergo in pathological conditions [60]. Senescent cells expressing p21 contribute to the progression of lung pathologies, such as lung fibrosis and chronic lung inflammation [54, 68–72]. Moreover, p21 deficiency reduces senescent cells accumulation and protects mice from liver fibrosis [60], renal fibrosis [73], osteopenia [74] and chronic pancreatitis [75]. Emerging evidence suggests that various types of cells undergo senescence and upregulate p21 in lungs affected by chronic bronchitis, including in alveolar cells, bronchial epithelial cells, fibroblasts, smooth muscle cells, endothelial cells, and leukocytes [72, 76–78]. These cells can contribute to the inflammatory response following chronic LPS exposure through secretion of SASP in an autocrine and paracrine manner. However, the specific contribution of each cell type population to the p21-dependent inflammatory response has yet been examined.
Senescent cells of different cell types might affect the progression of chronic bronchitis differently. Senescence of both adaptive and innate immune cells (immunosenescence) may result in impaired immune responses and aggravate lung inflammation [79–81]. However, here we showed that the p21-dependent inflammatory response following chronic LPS exposure is not mediated by the immune cells, but rather by resident epithelial and endothelial cells in the lungs. Indeed, epithelial and endothelial cells accumulate in the lungs during chronic bronchitis [76, 82]. These cells may contribute differently to the progression of chronic bronchitis pathology. Accumulation of p21-positive senescent airway epithelial cells during chronic bronchitis results in impaired lung epithelial regeneration, in addition to accelerated p38 MAPK-dependent airway inflammation [83]. Deleting p21 in epithelial cells also attenuates DNA damage, inflammatory responses, and oxidative stress induced by chronic bronchitis [54, 84]. p21 expression in senescent endothelial cells contributes to the inflammatory response by secretion of SASP [77, 85], while promoting immune cells adhesion by expression of ICAM-1 [77]. Overall, p21 in epithelial and endothelial cells potentiates the pro-inflammatory response in chronic bronchitis following repeated injury by various mechanisms. Ablation of p21 disrupts the accumulation of senescent cells and alleviates chronic lung inflammation.
Elimination of senescent cells may allow for better recovery of tissues from damage. Indeed, senescent cells elimination via transgenic techniques improves aging-associated phenotypes [31, 38, 86]. Interestingly, p21 knockout supports tissue regeneration and limits aging phenotypes in several mouse models. For example, p21 prolongs the lifespan of telomerase deficient mice by rescuing the proliferation of intestinal progenitor cells and improving the repopulation capacity of hematopoietic stem cells [87]. Additionally, p21 deficiency supports tissue regeneration, including regeneration of articular cartilage [88, 89], bone regeneration following injury [90] and skin wound healing [91]. Lastly, p21 attenuates the phenotype of aged BubR1 hypomorphic progeroid mice [92]. Therefore, we suggest that p21-dependent elimination of senescent cells may limit the damage induced by the pro-inflammatory presence of senescent cells, but also promote tissue regeneration. Therefore, inhibition of p21 represents a promising strategy for limiting age-related inflammatory disorders in general and obstructive lung diseases in particular.
Materials and Methods
Contact for reagent and resource sharing
Further information and requests for reagents may be directed to, and will be fulfilled by the corresponding author, Dr. Valery Krizhanovsky (
Methods and protocols
Mice
Female C57BL/6J mice were purchased from Harlan Laboratories. p21-/- mice [B6.129S6(Cg)-Cdkn1a
Chronic lung inflammation induction
Chronic LPS exposure regime was performed to age-matched C57BL/6J and p21−/− mice that were exposed to PBS or PBS containing Escherichia coli LPS (L2630, Sigma-Aldrich) aerosol, as previously described by us [26]. Following 48 hours from the last exposure all mice were harvested. BAL fluid and right lung lobes were subjected to fluorescence-activated cell sorting (FACS) analysis. Left lung lobes were taken for cell sorting, RNA and protein extraction, fixation using 4% paraformaldehyde (PFA) for immunohistology or embedded in OCT solution for cryosectioning and SA-β-gal stains.
Single-cell lung homogenate preparation, FACS analysis and sorting
The lungs of euthanized mice were removed, washed in RPMI medium (11875-093, Thermo Fisher Scientific), minced, and incubated at 37C for 45min in RPMI medium containing 1mg/ml collagenase type 4 (LS004189, Worthington), and 0.02mg/ml DNaseI (04536282001, Roche). Lung cell suspensions were pushed through a 100-μm cell strainer and spun, and red blood cells were lysed with red blood cell lysis buffer (R7757, Sigma-Aldrich). Cells from whole lungs were collected, as well as BAL fluid, washed twice with FACS buffer, and immunolabeled with antibodies against B220 (PerCP/Cy5.5, #103236), CD4 (APC, #100412), CD8a (FITC, #100706), CD11b (PerCP/Cy5.5, #101228), CD11c (APC, #117310), CD24 (PE-CY7, #101822), CD45 (PB, #103126 or FITC, #103108), CD103 (PE, #121406), Ly6G (APC, #127614), MHC-II (APC-CY7, #107627), NKp46 (PE, #137604), TCR-b (PE, #109208) (all from BioLegend) and Siglec-F (PB, #562681, BD Biosciences). The cells were run in a LSR II Flow Cytometer (BD Biosciences) and analyzed using the FlowJo v10 software (BD Biosciences). The gating strategy for the immune subsets was performed as follows: neutrophils (CD45+/Ly6G+/CD11b+), NK cells (CD45+/NKp46+), CD3 (CD45+/TCRb+), CD4 (CD45+/TCRb+/CD4+), CD8 (CD45+/TCRb+/CD8+), B cells (CD45+/B220+), alveolar macrophages (CD45+/CD11c+/SiglecF+) and interstitial macrophages (CD45+/MHCII+, CD11c+, Siglec-F-, CD11b+, CD24-).
For sorting, the lungs of euthanized mice were perfused with cold PBS and then removed. The lungs were then washed in DMEM/F12 media (11330-032, Invitrogen), minced, and incubated at 37C for 20min in DMEM/F12 medium containing 3U/ml elastase (LS002279, Worthington), and 0.33U/ml DNaseI (04536282001, Roche). Lung cell suspensions were then washed with DMEM/F12 supplemented with 100units/ml of penicillin, 100mg/ml of streptomycin (03-031-1B, Biological Industries) and 10% fetal bovine serum (FBS) (Thermo Fisher Scientific), pushed through a 100-um cell strainer and spun. Red blood cells were lysed with ACK lysis buffer (A1049201, Gibco). The cells were then centrifuged, re-suspended in FACS buffer and immunolabeled with antibodies against CD31 (BV605, #102427), CD45 (PE, #103106), EpCam (Alexa Fluor 488, #118210) (all from BioLegend) and TER-119 (eFluor450, #48-5921-82, eBioscience). Cells were stained with SytoxBlue (S34857, Thermo Fisher Scientific) and were sorted using FACS Aria Fusion (BD Biosciences) directly into RLT buffer (74104, QIAGEN) for qPCR or Lysis/Binding buffer (61012, Thermo Fisher Scientific) for subsequent RNA sequencing. The gating strategy for the immune subsets was performed as follows: epithelial cells (CD45-/CD31-/EpCam+), endothelial cells (CD45-/CD31+/EpCam-) and immune cells (CD45+/CD31-/EpCam-).
Quantitative RT-PCR
Total RNA from mice lungs was extracted using NucleoSpin RNA Mini kit (740955.50, Macherey-Nagel). cDNA was produced using random hexamers (N8080127, Thermo-Fisher). The cDNA samples were amplified using Platinum SYBR Green qPCR SuperMix (11744-500, Life Technologies) in a StepOnePlus Real-Time PCR System (Applied Biosystems). Relative expression was normalized using the expression levels of GAPDH. Primer sequences can be found in Supplementary Table 1.
Immunoblotting
Cell lysates (15–30 mg of protein) were resolved by 12.5% SDS–PAGE and transferred onto ImmobilonP membranes (IPVH00010, Millipore). After blocking of the membranes with 5% bovine serum albumin (BSA) in TBST (Tris-buffered saline with 0.01% Tween-20) for 1 hour, they were probed with antibodies against p15 (ab53035), p16 (ab108349) (both from Abcam), cleaved-caspase 3 (#9661, Cell Signaling Technology), GAPDH (MAB374, Millipore), collagen-1 (600-401-103-0.5, Rockland Immunochemicals), p21 (sc-6246, Santa Cruz) and p53 (mix of DO-1 and PAb1801, kindly provided by M. Oren, Weizmann Institute of Science). The blots were developed using either SuperSignal West Pico PLUS chemiluminescent substrate (#34579) or SuperSignal West Femto maximum sensitivity substrate (#34095) (both from Thermo Fisher Scientific). The blots were analyzed using the Image Lab software (Bio-Rad Laboratories).
Histological analysis
Immunofluorescence (IF) was performed on 4-μm paraffin sections and immunohistochemistry (IHC) was performed on 2-μm paraffin sections, according to standard procedures. Primary antibodies recognizing CD3 (1:50, MCA500G, Bio-Rad), CD45R (#14-0452-82, eBioscience, 1:50), cleaved caspase 3 (1:50, #9661, Cell Signaling Technology) and p21 (1:50, 556431, BD Biosciences) were applied overnight at 4C. Cleaved caspase 3 and p21 sections were developed using DAB (SK-4100, Vector Laboratories) followed by hematoxylin counterstaining. CD3 and CD45R sections were incubated with Cy3 anti-rabbit antibody (111-165-144, Jackson ImmunoResearch) for 90min in a humidity chamber. Sections were counterstained by DAPI (D9542, Sigma-Aldrich) and then mounted using Fluormount (0100-01, Southern Biotech). Cleaved caspase 3 and p21 stained sections were examined and photographed with a bright field microscope (Olympus IX81) and images were visualized using the CellP software (Diagnostic Instruments). CD3 and CD45R stained sections were visualized and photographed with a fluorescence microscope (Eclipse Ni-U; Nikon) equipped with Plan Fluor objectives (20x; 40x; 60x) connected to a color camera (DS-Ri1, Nikon) microscope. Images were analyzed using the Image Pro+ software (Nikon).
Paraffin-embedded tissue sections (4-μm) were stained with hematoxylin–eosin (H&E) for routine examination. From these stained slides bronchial epithelial cell wall width were calculated. 10 random points were chosen and the average of the cell wall width was measured. The sections were also stained for with Sirius Red for visualization of fibrotic deposition. The sections were examined and photographed with a bright-field microscope (Olympus). These images were then quantified using the NIH ImageJ software (http://rsb.info.nih.gov/ij/).
For SA-β-gal staining, 14-μm cryosections of OCT embedded mouse lungs were fixed in 0.5% glutaraldehyde (50-261-96, Electron Microscopy Sciences) for 15min, washed with PBS supplemented with 1mM MgCl2 in PBS at pH 5.5, and incubated for 6 hours in X-Gal staining solution (1mg/mL X-Gal (1758-0300, Inalco Pharmaceuticals), 5mM potassium ferrocyanide (104973, Merck Millipore), 5mM potassium ferricyanide (P3289, Sigma-Aldrich) and 1mM MgCl2 in PBS at pH 5.5). Sections were counterstained with Nuclear Fast Red (N3020, Sigma-Aldrich), dehydrated and mounted with Eukitt mounting solution (03989, Sigma-Aldrich). Sections were visualized using an Olympus microscope, and images were analyzed using the CellP software (Diagnostic Instruments).
Treadmill assay
Mice were tested before and during chronic LPS exposure. The assay was performed on a treadmill (Panlab Mouse 5-Lane Treadmill; model#: 760309; HARVARD APPARATUS, USA) at 5 time points: baseline point (3 days before the first LPS exposure) and every other week from week 3 until week 9. Training protocol was described [93] and modified as follows: mice run on the treadmill with no incline, starting at a speed of 15 meters/min for 10min. Afterwards, the speed was increased gradually, over 10min, to a final speed of 35 meters/min. Finally, mice run at a speed of 35 meters/min for another 10min. The maximal actual distance was determined from the starting point until exhaustion (defined as 2 consecutive 30sec intervals in which the mouse steps off the lane at least 5 times).
MARS-Seq sequencing and data analysis
Libraries were prepared from pooled samples of the same cell type (10000 cells per sample) according to a bulk variation of MARS-Seq [94] and were sequenced on Illumina NextSeq500 (Illumina). Data was analyzed using a User-friendly Transcriptome Analysis Pipeline (UTAP) [95]. Only reads with unique mapping to the 3’ of RefSeq annotated genes (mm10, NCBI Mus musculus Annotation Release 109) were analyzed. Gene expression was calculated and normalized using DESeq2 with batch correction. A gene was defined differentially-expressed (DE) according to the following thresholds of adjusted p value ≤0.05, average normalized reads one of all samples ≥5, and absolute value of fold-change ≥2. We selected differentially-expressed genes (DEGs) from the following comparisons: LPS WT- PBS WT, LPS p21−/−- PBS p21−/− and LPS p21−/−-LPS WT. We performed partitioning clustering using the k-means algorithm (Euclidian distance) using Partek Genomics Suite® software (version 7.0). Gene ontology term enrichment analysis was performed using the Metascape software [96].
Cytoscape pathways analysis
Pathway analysis was performed using Cytoscape [97] and STRING application [98] for exploring protein-protein interactions and Enrichment Map application to postulate pathways. Differentially expressed genes having adjusted p-value≤0.05 and Log2 fold change≥0.5 were used as input to project differentially expressed pathways in each comparison. Pathways having FDR value smaller than 0.05 are displayed.
Statistical analysis
Data are presented as mean ±SEM. Statistical significance was performed by the Prism software (GraphPad) and determined using Student’s t-test for pairwise comparisons or one-way ANOVA for comparing multiple groups, as well as by the R software (version 4.0.3) to determine chi-squared test for comparing the number of DEGs. Mann Whitney u test was performed in order to test the statistical significance of the physical activity and fitness of the mice. A P-value of < 0.05 was considered significant.
Data availability statement
The data that support the findings of this study is available from the corresponding author, V.K., on reasonable request.
Abbreviations
BAL: Bronchoalveolar lavage; CDK: Cyclin-dependent kinases; CDKI: Cyclin-dependent kinase inhibitor; CDKN1A: Cyclin-dependent kinase inhibitor 1A (p21); COPD: Chronic obstructive pulmonary disease; DEG: Differentially-expressed genes; GO: Gene Ontology; H&E: Hematoxylin and eosin; iBALT: Inducible bronchus-associated lymphoid tissue; LPS: Lipopolysaccharide; MMPs: Matrix Metalloproteases; NF-κB: Nuclear factor kappa-light-chain-enhancer of activated B cells; p21-/-: p21 knockout; Rb: retinoblastoma protein; SASP: Senescence-associated secretory phenotype; SA-β-gal: Senescence-associated beta-galactosidase; WT: Wild type.
Author Contributions
N.L. and N.P. designed and performed experiments and analysed the data. J.M, L.R., Y.O. and O.R. performed experiments. N.W. and H.A. performed bioinformatics analysis. R.E. conducted and analyzed the immunohistochemical experiments. M.T. participated in the design and analysis of the mouse LPS-inhalation experimental settings and treadmill assay. R.R. performed statistical analysis. N.L., N.P. and V.K. analysed the experiments and wrote the manuscript. V.K. supervised the project.
Acknowledgments
We thank all the members of the Krizhanovsky laboratory for stimulating discussions.
Conflicts of Interest
V.K. is an author of patents on senolytics and senolytic approaches and consultant for Sentaur Bio. None of these influenced the work conducted or analysis of data presented in this manuscript.
Ethical Statement
All mice were housed and maintained in the Animal Facility of the Weizmann Institute of Science during the experiments. All mice experiments were approved by the Weizmann Institute Animal Care and Use Committee (IACUC).
Funding
This work was supported by grants to V.K. from the European Research Council under the European Union’s H2020 (856487), from the Israel Science Foundation (1626/20), DFG (CRC 1506), Israel Ministry of Health, Minerva Center “Aging, from Physical Materials to Human Tissues”, Quinquin Foundation, Belle S. and Irving E. Meller Center for the Biology of Aging and Sagol Institute for Longevity Research. M.T. is the incumbent of the Carolito Stiftung Research Fellow Chair in Neurodegenerative Diseases. V.K. is an incumbent of The Georg F. Duckwitz Professorial Chair and Shimon and Golde Picker – Weizmann Award.
References
- 1. Barnes PJ, Burney PG, Silverman EK, Celli BR, Vestbo J, Wedzicha JA, Wouters EF. Chronic obstructive pulmonary disease. Nat Rev Dis Primers. 2015; 1:15076. https://doi.org/10.1038/nrdp.2015.76 [PubMed]
- 2. WHO. Global Health Estimates: Life expectancy and leading causes of death and disability. (https://www.who.int/data/gho/data/themes/mortality-and-global-health-estimates: World Health Organization). 2019.
- 3. Vogelmeier CF, Criner GJ, Martínez FJ, Anzueto A, Barnes PJ, Bourbeau J, Celli BR, Chen R, Decramer M, Fabbri LM, Frith P, Halpin DM, López Varela MV, et al. Global Strategy for the Diagnosis, Management, and Prevention of Chronic Obstructive Lung Disease 2017 Report: GOLD Executive Summary. Arch Bronconeumol. 2017; 53:128–49. https://doi.org/10.1016/j.arbres.2017.02.001 [PubMed]
- 4. Higham A, Quinn AM, Cançado JED, Singh D. The pathology of small airways disease in COPD: historical aspects and future directions. Respir Res. 2019; 20:49. https://doi.org/10.1186/s12931-019-1017-y [PubMed]
- 5. Barnes PJ. Targeting cellular senescence as a new approach to chronic obstructive pulmonary disease therapy. Curr Opin Pharmacol. 2021; 56:68–73. https://doi.org/10.1016/j.coph.2020.11.004 [PubMed]
- 6. Meiners S, Eickelberg O, Königshoff M. Hallmarks of the ageing lung. Eur Respir J. 2015; 45:807–27. https://doi.org/10.1183/09031936.00186914 [PubMed]
- 7. Freeman CM, Curtis JL. Lung Dendritic Cells: Shaping Immune Responses throughout Chronic Obstructive Pulmonary Disease Progression. Am J Respir Cell Mol Biol. 2017; 56:152–9. https://doi.org/10.1165/rcmb.2016-0272TR [PubMed]
- 8. Givi ME, Folkerts G, Wagenaar GT, Redegeld FA, Mortaz E. Cigarette smoke differentially modulates dendritic cell maturation and function in time. Respir Res. 2015; 16:131. https://doi.org/10.1186/s12931-015-0291-6 [PubMed]
- 9. Wang Y, Xu J, Meng Y, Adcock IM, Yao X. Role of inflammatory cells in airway remodeling in COPD. Int J Chron Obstruct Pulmon Dis. 2018; 13:3341–8. https://doi.org/10.2147/COPD.S176122 [PubMed]
- 10. Bu T, Wang LF, Yin YQ. How Do Innate Immune Cells Contribute to Airway Remodeling in COPD Progression? Int J Chron Obstruct Pulmon Dis. 2020; 15:107–16. https://doi.org/10.2147/COPD.S235054 [PubMed]
- 11. Gao W, Li L, Wang Y, Zhang S, Adcock IM, Barnes PJ, Huang M, Yao X. Bronchial epithelial cells: The key effector cells in the pathogenesis of chronic obstructive pulmonary disease? Respirology. 2015; 20:722–9. https://doi.org/10.1111/resp.12542 [PubMed]
- 12. Szucs B, Szucs C, Petrekanits M, Varga JT. Molecular Characteristics and Treatment of Endothelial Dysfunction in Patients with COPD: A Review Article. Int J Mol Sci. 2019; 20:4329. https://doi.org/10.3390/ijms20184329 [PubMed]
- 13. Barnes PJ. Small airway fibrosis in COPD. Int J Biochem Cell Biol. 2019; 116:105598. https://doi.org/10.1016/j.biocel.2019.105598 [PubMed]
-
14.
Hiemstra PS, McCray PB
Jr , Bals R. The innate immune function of airway epithelial cells in inflammatory lung disease. Eur Respir J. 2015; 45:1150–62. https://doi.org/10.1183/09031936.00141514 [PubMed] - 15. Inui T, Watanabe M, Nakamoto K, Sada M, Hirata A, Nakamura M, Honda K, Ogawa Y, Takata S, Yokoyama T, Saraya T, Kurai D, Wada H, et al. Bronchial epithelial cells produce CXCL1 in response to LPS and TNFα: A potential role in the pathogenesis of COPD. Exp Lung Res. 2018; 44:323–31. https://doi.org/10.1080/01902148.2018.1520936 [PubMed]
- 16. Aghapour M, Raee P, Moghaddam SJ, Hiemstra PS, Heijink IH. Airway Epithelial Barrier Dysfunction in Chronic Obstructive Pulmonary Disease: Role of Cigarette Smoke Exposure. Am J Respir Cell Mol Biol. 2018; 58:157–69. https://doi.org/10.1165/rcmb.2017-0200TR [PubMed]
- 17. Lu Q, Gottlieb E, Rounds S. Effects of cigarette smoke on pulmonary endothelial cells. Am J Physiol Lung Cell Mol Physiol. 2018; 314:L743–56. https://doi.org/10.1152/ajplung.00373.2017 [PubMed]
- 18. Orosz Z, Csiszar A, Labinskyy N, Smith K, Kaminski PM, Ferdinandy P, Wolin MS, Rivera A, Ungvari Z. Cigarette smoke-induced proinflammatory alterations in the endothelial phenotype: role of NAD(P)H oxidase activation. Am J Physiol Heart Circ Physiol. 2007; 292:H130–9. https://doi.org/10.1152/ajpheart.00599.2006 [PubMed]
- 19. Shen Y, Rattan V, Sultana C, Kalra VK. Cigarette smoke condensate-induced adhesion molecule expression and transendothelial migration of monocytes. Am J Physiol. 1996; 270:H1624–33. https://doi.org/10.1152/ajpheart.1996.270.5.H1624 [PubMed]
- 20. Overbeek SA, Braber S, Henricks PA, Kleinjan M, Kamp VM, Georgiou NA, Garssen J, Kraneveld AD, Folkerts G. Cigarette smoke induces β2-integrin-dependent neutrophil migration across human endothelium. Respir Res. 2011; 12:75. https://doi.org/10.1186/1465-9921-12-75 [PubMed]
- 21. Sharma J, Young DM, Marentette JO, Rastogi P, Turk J, McHowat J. Lung endothelial cell platelet-activating factor production and inflammatory cell adherence are increased in response to cigarette smoke component exposure. Am J Physiol Lung Cell Mol Physiol. 2012; 302:L47–55. https://doi.org/10.1152/ajplung.00179.2011 [PubMed]
- 22. Yoon YS, Jin M, Sin DD. Accelerated lung aging and chronic obstructive pulmonary disease. Expert Rev Respir Med. 2019; 13:369–80. https://doi.org/10.1080/17476348.2019.1580576 [PubMed]
- 23. Easter M, Bollenbecker S, Barnes JW, Krick S. Targeting Aging Pathways in Chronic Obstructive Pulmonary Disease. Int J Mol Sci. 2020; 21:6924. https://doi.org/10.3390/ijms21186924 [PubMed]
- 24. Kukrety SP, Parekh JD, Bailey KL. Chronic obstructive pulmonary disease and the hallmarks of aging. Lung India. 2018; 35:321–7. https://doi.org/10.4103/lungindia.lungindia_266_17 [PubMed]
- 25. Baker JR, Donnelly LE, Barnes PJ. Senotherapy: A New Horizon for COPD Therapy. Chest. 2020; 158:562–70. https://doi.org/10.1016/j.chest.2020.01.027 [PubMed]
- 26. Sagiv A, Bar-Shai A, Levi N, Hatzav M, Zada L, Ovadya Y, Roitman L, Manella G, Regev O, Majewska J, Vadai E, Eilam R, Feigelson SW, et al. p53 in Bronchial Club Cells Facilitates Chronic Lung Inflammation by Promoting Senescence. Cell Rep. 2018; 22:3468–79. https://doi.org/10.1016/j.celrep.2018.03.009 [PubMed]
- 27. Childs BG, Gluscevic M, Baker DJ, Laberge RM, Marquess D, Dananberg J, van Deursen JM. Senescent cells: an emerging target for diseases of ageing. Nat Rev Drug Discov. 2017; 16:718–35. https://doi.org/10.1038/nrd.2017.116 [PubMed]
- 28. 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:75–95. https://doi.org/10.1038/s41580-020-00314-w [PubMed]
- 29. Gorgoulis V, Adams PD, Alimonti A, Bennett DC, Bischof O, Bishop C, Campisi J, Collado M, Evangelou K, Ferbeyre G, Gil J, Hara E, Krizhanovsky V, et al. Cellular Senescence: Defining a Path Forward. Cell. 2019; 179:813–27. https://doi.org/10.1016/j.cell.2019.10.005 [PubMed]
- 30. Krizhanovsky V, Yon M, Dickins RA, Hearn S, Simon J, Miething C, Yee H, Zender L, Lowe SW. Senescence of activated stellate cells limits liver fibrosis. Cell. 2008; 134:657–67. https://doi.org/10.1016/j.cell.2008.06.049 [PubMed]
- 31. Demaria M, Ohtani N, Youssef SA, Rodier F, Toussaint W, Mitchell JR, Laberge RM, Vijg J, Van Steeg H, Dollé ME, Hoeijmakers JH, de Bruin A, Hara E, Campisi J. An essential role for senescent cells in optimal wound healing through secretion of PDGF-AA. Dev Cell. 2014; 31:722–33. https://doi.org/10.1016/j.devcel.2014.11.012 [PubMed]
- 32. Serrano M, Lin AW, McCurrach ME, Beach D, Lowe SW. Oncogenic ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a. Cell. 1997; 88:593–602. https://doi.org/10.1016/s0092-8674(00)81902-9 [PubMed]
- 33. Gal H, Lysenko M, Stroganov S, Vadai E, Youssef SA, Tzadikevitch-Geffen K, Rotkopf R, Biron-Shental T, de Bruin A, Neeman M, Krizhanovsky V. Molecular pathways of senescence regulate placental structure and function. EMBO J. 2019; 38:e100849. https://doi.org/10.15252/embj.2018100849 [PubMed]
- 34. Chuprin A, Gal H, Biron-Shental T, Biran A, Amiel A, Rozenblatt S, Krizhanovsky V. Cell fusion induced by ERVWE1 or measles virus causes cellular senescence. Genes Dev. 2013; 27:2356–66. https://doi.org/10.1101/gad.227512.113 [PubMed]
- 35. Muñoz-Espín D, Serrano M. Cellular senescence: from physiology to pathology. Nat Rev Mol Cell Biol. 2014; 15:482–96. https://doi.org/10.1038/nrm3823 [PubMed]
- 36. Storer M, Mas A, Robert-Moreno A, Pecoraro M, Ortells MC, Di Giacomo V, Yosef R, Pilpel N, Krizhanovsky V, Sharpe J, Keyes WM. Senescence is a developmental mechanism that contributes to embryonic growth and patterning. Cell. 2013; 155:1119–30. https://doi.org/10.1016/j.cell.2013.10.041 [PubMed]
- 37. Baker DJ, Wijshake T, Tchkonia T, LeBrasseur NK, Childs BG, van de Sluis B, Kirkland JL, van Deursen JM. Clearance of p16Ink4a-positive senescent cells delays ageing-associated disorders. Nature. 2011; 479:232–6. https://doi.org/10.1038/nature10600 [PubMed]
- 38. Baker DJ, Childs BG, Durik M, Wijers ME, Sieben CJ, Zhong J, Saltness RA, Jeganathan KB, Verzosa GC, Pezeshki A, Khazaie K, Miller JD, van Deursen JM. Naturally occurring p16(Ink4a)-positive cells shorten healthy lifespan. Nature. 2016; 530:184–9. https://doi.org/10.1038/nature16932 [PubMed]
- 39. Ovadya Y, Landsberger T, Leins H, Vadai E, Gal H, Biran A, Yosef R, Sagiv A, Agrawal A, Shapira A, Windheim J, Tsoory M, Schirmbeck R, et al. Impaired immune surveillance accelerates accumulation of senescent cells and aging. Nat Commun. 2018; 9:5435. https://doi.org/10.1038/s41467-018-07825-3 [PubMed]
- 40. Acosta JC, Banito A, Wuestefeld T, Georgilis A, Janich P, Morton JP, Athineos D, Kang TW, Lasitschka F, Andrulis M, Pascual G, Morris KJ, Khan S, et al. A complex secretory program orchestrated by the inflammasome controls paracrine senescence. Nat Cell Biol. 2013; 15:978–90. https://doi.org/10.1038/ncb2784 [PubMed]
- 41. Ovadya Y, Krizhanovsky V. Senescent cells: SASPected drivers of age-related pathologies. Biogerontology. 2014; 15:627–42. https://doi.org/10.1007/s10522-014-9529-9 [PubMed]
- 42. Kumari R, Jat P. Mechanisms of Cellular Senescence: Cell Cycle Arrest and Senescence Associated Secretory Phenotype. Front Cell Dev Biol. 2021; 9:645593. https://doi.org/10.3389/fcell.2021.645593 [PubMed]
- 43. Lopes-Paciencia S, Saint-Germain E, Rowell MC, Ruiz AF, Kalegari P, Ferbeyre G. The senescence-associated secretory phenotype and its regulation. Cytokine. 2019; 117:15–22. https://doi.org/10.1016/j.cyto.2019.01.013 [PubMed]
- 44. Coppé JP, Desprez PY, Krtolica A, Campisi J. The senescence-associated secretory phenotype: the dark side of tumor suppression. Annu Rev Pathol. 2010; 5:99–118. https://doi.org/10.1146/annurev-pathol-121808-102144 [PubMed]
- 45. Kale A, Sharma A, Stolzing A, Desprez PY, Campisi J. Role of immune cells in the removal of deleterious senescent cells. Immun Ageing. 2020; 17:16. https://doi.org/10.1186/s12979-020-00187-9 [PubMed]
- 46. Rodier F, Campisi J. Four faces of cellular senescence. J Cell Biol. 2011; 192:547–56. https://doi.org/10.1083/jcb.201009094 [PubMed]
- 47. Sharpless NE, Sherr CJ. Forging a signature of in vivo senescence. Nat Rev Cancer. 2015; 15:397–408. https://doi.org/10.1038/nrc3960 [PubMed]
- 48. Liu X, Wan M. A tale of the good and bad: Cell senescence in bone homeostasis and disease. Int Rev Cell Mol Biol. 2019; 346:97–128. https://doi.org/10.1016/bs.ircmb.2019.03.005 [PubMed]
- 49. Salama R, Sadaie M, Hoare M, Narita M. Cellular senescence and its effector programs. Genes Dev. 2014; 28:99–114. https://doi.org/10.1101/gad.235184.113 [PubMed]
- 50. Fischer M, Müller GA. Cell cycle transcription control: DREAM/MuvB and RB-E2F complexes. Crit Rev Biochem Mol Biol. 2017; 52:638–62. https://doi.org/10.1080/10409238.2017.1360836 [PubMed]
- 51. Otto T, Sicinski P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer. 2017; 17:93–115. https://doi.org/10.1038/nrc.2016.138 [PubMed]
- 52. Chiappara G, Gjomarkaj M, Virzì A, Sciarrino S, Ferraro M, Bruno A, Montalbano AM, Vitulo P, Minervini MI, Pipitone L, Pace E. The role of p21 Waf1/Cip1 in large airway epithelium in smokers with and without COPD. Biochim Biophys Acta. 2013; 1832:1473–81. https://doi.org/10.1016/j.bbadis.2013.04.022 [PubMed]
- 53. Lv X, Liu C, Liu S, Li Y, Wang W, Li K, Hua F, Cui B, Zhang X, Yu J, Yu J, Hu Z. The cell cycle inhibitor P21 promotes the development of pulmonary fibrosis by suppressing lung alveolar regeneration. Acta Pharm Sin B. 2022; 12:735–46. https://doi.org/10.1016/j.apsb.2021.07.015 [PubMed]
- 54. Yao H, Yang SR, Edirisinghe I, Rajendrasozhan S, Caito S, Adenuga D, O’Reilly MA, Rahman I. Disruption of p21 attenuates lung inflammation induced by cigarette smoke, LPS, and fMLP in mice. Am J Respir Cell Mol Biol. 2008; 39:7–18. https://doi.org/10.1165/rcmb.2007-0342OC [PubMed]
- 55. Yao H, Chung S, Hwang JW, Rajendrasozhan S, Sundar IK, Dean DA, McBurney MW, Guarente L, Gu W, Rönty M, Kinnula VL, Rahman I. SIRT1 protects against emphysema via FOXO3-mediated reduction of premature senescence in mice. J Clin Invest. 2012; 122:2032–45. https://doi.org/10.1172/JCI60132 [PubMed]
- 56. Wang E. Senescent human fibroblasts resist programmed cell death, and failure to suppress bcl2 is involved. Cancer Res. 1995; 55:2284–92. [PubMed]
- 57. Sanders YY, Liu H, Zhang X, Hecker L, Bernard K, Desai L, Liu G, Thannickal VJ. Histone modifications in senescence-associated resistance to apoptosis by oxidative stress. Redox Biol. 2013; 1:8–16. https://doi.org/10.1016/j.redox.2012.11.004 [PubMed]
- 58. Chang J, Wang Y, Shao L, Laberge RM, Demaria M, Campisi J, Janakiraman K, Sharpless NE, Ding S, Feng W, Luo Y, Wang X, Aykin-Burns N, et al. Clearance of senescent cells by ABT263 rejuvenates aged hematopoietic stem cells in mice. Nat Med. 2016; 22:78–83. https://doi.org/10.1038/nm.4010 [PubMed]
- 59. Yosef R, Pilpel N, Tokarsky-Amiel R, Biran A, Ovadya Y, Cohen S, Vadai E, Dassa L, Shahar E, Condiotti R, Ben-Porath I, Krizhanovsky V. Directed elimination of senescent cells by inhibition of BCL-W and BCL-XL. Nat Commun. 2016; 7:11190. https://doi.org/10.1038/ncomms11190 [PubMed]
- 60. Yosef R, Pilpel N, Papismadov N, Gal H, Ovadya Y, Vadai E, Miller S, Porat Z, Ben-Dor S, Krizhanovsky V. p21 maintains senescent cell viability under persistent DNA damage response by restraining JNK and caspase signaling. EMBO J. 2017; 36:2280–95. https://doi.org/10.15252/embj.201695553 [PubMed]
- 61. Levi N, Papismadov N, Solomonov I, Sagi I, Krizhanovsky V. The ECM path of senescence in aging: components and modifiers. FEBS J. 2020; 287:2636–46. https://doi.org/10.1111/febs.15282 [PubMed]
- 62. Baggio C, Velazquez JV, Fragai M, Nordgren TM, Pellecchia M. Therapeutic Targeting of MMP-12 for the Treatment of Chronic Obstructive Pulmonary Disease. J Med Chem. 2020; 63:12911–20. https://doi.org/10.1021/acs.jmedchem.0c01285 [PubMed]
- 63. Gelinas J, Harper M, Sasso J, Wright S, Melzer B, Agar G, Guenette J, duManoir G, Roman M, Rolf JD, Eves N. Phenotyping Cardiopulmonary Exercise Limitations in Chronic Obstructive Pulmonary Disease. Front Physiol. 2022; 13:816586. https://doi.org/10.3389/fphys.2022.816586 [PubMed]
- 64. King PT. Inflammation in chronic obstructive pulmonary disease and its role in cardiovascular disease and lung cancer. Clin Transl Med. 2015; 4:68. https://doi.org/10.1186/s40169-015-0068-z [PubMed]
- 65. Barnes PJ, Chowdhury B, Kharitonov SA, Magnussen H, Page CP, Postma D, Saetta M. Pulmonary biomarkers in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2006; 174:6–14. https://doi.org/10.1164/rccm.200510-1659PP [PubMed]
- 66. Silver JS, Kearley J, Copenhaver AM, Sanden C, Mori M, Yu L, Pritchard GH, Berlin AA, Hunter CA, Bowler R, Erjefalt JS, Kolbeck R, Humbles AA. Inflammatory triggers associated with exacerbations of COPD orchestrate plasticity of group 2 innate lymphoid cells in the lungs. Nat Immunol. 2016; 17:626–35. https://doi.org/10.1038/ni.3443 [PubMed]
- 67. Ladjemi MZ, Martin C, Lecocq M, Detry B, Nana FA, Moulin C, Weynand B, Fregimilicka C, Bouzin C, Thurion P, Carlier F, Serré J, Gayan-Ramirez G, et al. Increased IgA Expression in Lung Lymphoid Follicles in Severe Chronic Obstructive Pulmonary Disease. Am J Respir Crit Care Med. 2019; 199:592–602. https://doi.org/10.1164/rccm.201802-0352OC [PubMed]
- 68. Schafer MJ, White TA, Iijima K, Haak AJ, Ligresti G, Atkinson EJ, Oberg AL, Birch J, Salmonowicz H, Zhu Y, Mazula DL, Brooks RW, Fuhrmann-Stroissnigg H, et al. Cellular senescence mediates fibrotic pulmonary disease. Nat Commun. 2017; 8:14532. https://doi.org/10.1038/ncomms14532 [PubMed]
- 69. Liu RM, Liu G. Cell senescence and fibrotic lung diseases. Exp Gerontol. 2020; 132:110836. https://doi.org/10.1016/j.exger.2020.110836 [PubMed]
- 70. Kuwano K, Araya J, Hara H, Minagawa S, Takasaka N, Ito S, Kobayashi K, Nakayama K. Cellular senescence and autophagy in the pathogenesis of chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF). Respir Investig. 2016; 54:397–406. https://doi.org/10.1016/j.resinv.2016.03.010 [PubMed]
- 71. 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:556–64. https://doi.org/10.1164/rccm.201810-1975TR [PubMed]
- 72. Takasaka N, Araya J, Hara H, Ito S, Kobayashi K, Kurita Y, Wakui H, Yoshii Y, Yumino Y, Fujii S, Minagawa S, Tsurushige C, Kojima J, et al. Autophagy induction by SIRT6 through attenuation of insulin-like growth factor signaling is involved in the regulation of human bronchial epithelial cell senescence. J Immunol. 2014; 192:958–68. https://doi.org/10.4049/jimmunol.1302341 [PubMed]
- 73. Megyesi J, Tarcsafalvi A, Li S, Hodeify R, Seng NS, Portilla D, Price PM. Increased expression of p21WAF1/CIP1 in kidney proximal tubules mediates fibrosis. Am J Physiol Renal Physiol. 2015; 308:F122–30. https://doi.org/10.1152/ajprenal.00489.2014 [PubMed]
- 74. Premnath P, Ferrie L, Louie D, Boyd S, Krawetz R. Absence of p21(WAF1/CIP1/SDI1) protects against osteopenia and minimizes bone loss after ovariectomy in a mouse model. PLoS One. 2019; 14:e0215018. https://doi.org/10.1371/journal.pone.0215018 [PubMed]
- 75. Seleznik GM, Reding T, Peter L, Gupta A, Steiner SG, Sonda S, Verbeke CS, Dejardin E, Khatkov I, Segerer S, Heikenwalder M, Graf R. Development of autoimmune pancreatitis is independent of CDKN1A/p21-mediated pancreatic inflammation. Gut. 2018; 67:1663–73. https://doi.org/10.1136/gutjnl-2016-313458 [PubMed]
- 76. Rutten EP, Gopal P, Wouters EF, Franssen FM, Hageman GJ, Vanfleteren LE, Spruit MA, Reynaert NL. Various Mechanistic Pathways Representing the Aging Process Are Altered in COPD. Chest. 2016; 149:53–61. https://doi.org/10.1378/chest.15-0645 [PubMed]
- 77. Amsellem V, Gary-Bobo G, Marcos E, Maitre B, Chaar V, Validire P, Stern JB, Noureddine H, Sapin E, Rideau D, Hue S, Le Corvoisier P, Le Gouvello S, et al. Telomere dysfunction causes sustained inflammation in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2011; 184:1358–66. https://doi.org/10.1164/rccm.201105-0802OC [PubMed]
- 78. Ahmad T, Sundar IK, Tormos AM, Lerner CA, Gerloff J, Yao H, Rahman I. Shelterin Telomere Protection Protein 1 Reduction Causes Telomere Attrition and Cellular Senescence via Sirtuin 1 Deacetylase in Chronic Obstructive Pulmonary Disease. Am J Respir Cell Mol Biol. 2017; 56:38–49. https://doi.org/10.1165/rcmb.2016-0198OC [PubMed]
- 79. Sadighi Akha AA. Aging and the immune system: An overview. J Immunol Methods. 2018; 463:21–6. https://doi.org/10.1016/j.jim.2018.08.005 [PubMed]
- 80. Maté I, Martínez de Toda I, Arranz L, Álvarez-Sala JL, De la Fuente M. Accelerated immunosenescence, oxidation and inflammation lead to a higher biological age in COPD patients. Exp Gerontol. 2021; 154:111551. https://doi.org/10.1016/j.exger.2021.111551 [PubMed]
- 81. Cho WK, Lee CG, Kim LK. COPD as a Disease of Immunosenescence. Yonsei Med J. 2019; 60:407–13. https://doi.org/10.3349/ymj.2019.60.5.407 [PubMed]
- 82. Tsuji T, Aoshiba K, Nagai A. Alveolar cell senescence in patients with pulmonary emphysema. Am J Respir Crit Care Med. 2006; 174:886–93. https://doi.org/10.1164/rccm.200509-1374OC [PubMed]
- 83. Zhou F, Onizawa S, Nagai A, Aoshiba K. Epithelial cell senescence impairs repair process and exacerbates inflammation after airway injury. Respir Res. 2011; 12:78. https://doi.org/10.1186/1465-9921-12-78 [PubMed]
- 84. Yao H, Sundar IK, Gorbunova V, Rahman I. P21-PARP-1 pathway is involved in cigarette smoke-induced lung DNA damage and cellular senescence. PLoS One. 2013; 8:e80007. https://doi.org/10.1371/journal.pone.0080007 [PubMed]
- 85. Paschalaki K, Rossios C, Pericleous C, MacLeod M, Rothery S, Donaldson GC, Wedzicha JA, Gorgoulis V, Randi AM, Barnes PJ. Inhaled corticosteroids reduce senescence in endothelial progenitor cells from patients with COPD. Thorax. 2022; 77:616–20. https://doi.org/10.1136/thoraxjnl-2020-216807 [PubMed]
- 86. Folgueras AR, Freitas-Rodríguez S, Velasco G, López-Otín C. Mouse Models to Disentangle the Hallmarks of Human Aging. Circ Res. 2018; 123:905–24. https://doi.org/10.1161/CIRCRESAHA.118.312204 [PubMed]
- 87. Choudhury AR, Ju Z, Djojosubroto MW, Schienke A, Lechel A, Schaetzlein S, Jiang H, Stepczynska A, Wang C, Buer J, Lee HW, von Zglinicki T, Ganser A, et al. Cdkn1a deletion improves stem cell function and lifespan of mice with dysfunctional telomeres without accelerating cancer formation. Nat Genet. 2007; 39:99–105. https://doi.org/10.1038/ng1937 [PubMed] Erratum in: Nat Genet. 2020; 52:548. https://doi.org/10.1038/s41588-020-0593-6 [PubMed]
- 88. Ibaraki K, Hayashi S, Kanzaki N, Hashimoto S, Kihara S, Haneda M, Takeuchi K, Niikura T, Kuroda R. Deletion of p21 expression accelerates cartilage tissue repair via chondrocyte proliferation. Mol Med Rep. 2020; 21:2236–42. https://doi.org/10.3892/mmr.2020.11028 [PubMed]
- 89. Jablonski CL, Besler BA, Ali J, Krawetz RJ. p21-/- Mice Exhibit Spontaneous Articular Cartilage Regeneration Post-Injury. Cartilage. 2021; 13:1608S–17S. https://doi.org/10.1177/1947603519876348 [PubMed]
- 90. Premnath P, Jorgenson B, Hess R, Tailor P, Louie D, Taiani J, Boyd S, Krawetz R. p21-/- mice exhibit enhanced bone regeneration after injury. BMC Musculoskelet Disord. 2017; 18:435. https://doi.org/10.1186/s12891-017-1790-z [PubMed]
- 91. Jiang D, de Vries JC, Muschhammer J, Schatz S, Ye H, Hein T, Fidan M, Romanov VS, Rinkevich Y, Scharffetter-Kochanek K. Local and transient inhibition of p21 expression ameliorates age-related delayed wound healing. Wound Repair Regen. 2020; 28:49–60. https://doi.org/10.1111/wrr.12763 [PubMed]
- 92. Baker DJ, Weaver RL, van Deursen JM. p21 both attenuates and drives senescence and aging in BubR1 progeroid mice. Cell Rep. 2013; 3:1164–74. https://doi.org/10.1016/j.celrep.2013.03.028 [PubMed]
- 93. Umansky KB, Gruenbaum-Cohen Y, Tsoory M, Feldmesser E, Goldenberg D, Brenner O, Groner Y. Runx1 Transcription Factor Is Required for Myoblasts Proliferation during Muscle Regeneration. PLoS Genet. 2015; 11:e1005457. https://doi.org/10.1371/journal.pgen.1005457 [PubMed]
- 94. Jaitin DA, Kenigsberg E, Keren-Shaul H, Elefant N, Paul F, Zaretsky I, Mildner A, Cohen N, Jung S, Tanay A, Amit I. Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science. 2014; 343:776–9. https://doi.org/10.1126/science.1247651 [PubMed]
- 95. Kohen R, Barlev J, Hornung G, Stelzer G, Feldmesser E, Kogan K, Safran M, Leshkowitz D. UTAP: User-friendly Transcriptome Analysis Pipeline. BMC Bioinformatics. 2019; 20:154. https://doi.org/10.1186/s12859-019-2728-2 [PubMed]
- 96. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019; 10:1523. https://doi.org/10.1038/s41467-019-09234-6 [PubMed]
- 97. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13:2498–504. https://doi.org/10.1101/gr.1239303 [PubMed]
- 98. Szklarczyk D, Kirsch R, Koutrouli M, Nastou K, Mehryary F, Hachilif R, Gable AL, Fang T, Doncheva NT, Pyysalo S, Bork P, Jensen LJ, von Mering C. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023; 51:D638–46. https://doi.org/10.1093/nar/gkac1000 [PubMed]