Abstract
Background SARS-CoV-2 infection of the respiratory system can progress to a multi-systemic disease with aberrant inflammatory response. Cellular senescence promotes chronic inflammation, named as senescence-associated secretory phenotype (SASP). We investigated whether COVID-19 disease is associated with cellular senescence and SASP.
Methods Autopsy lung tissue samples from 11 COVID-19 patients and 43 age-matched non-COVID controls with similar comorbidities were analysed by immunohistochemistry for SARS-CoV-2, markers of senescence and key SASP cytokines. Virally-induced senescence was functionally recapitulated in vitro, by infecting epithelial Vero-E6 cells and a three-dimensional alveosphere system of alveolar type 2 (AT2) cells with SARS-CoV-2 strains isolated from COVID-19 patients.
Results SARS-CoV-2 was detected by immunocytochemistry and electron microscopy predominantly in AT2 cells. Infected AT2 cells expressed the angiotensin-converting-enzyme 2 (ACE2) and exhibited increased senescence (p16INK4A and SenTraGorTM positivity) and IL-1β and IL-6 expression. In vitro, infection of Vero-E6 cells with SARS-CoV-2 induced senescence (SenTraGorTM), DNA damage (γ-H2AX) and increased cytokine (IL-1β, IL-6, CXCL8) and Apolipoprotein B mRNA-editing (APOBEC) enzyme expression. Next-generation-sequencing analysis of progenies obtained from infected/senescent Vero-E6 cells demonstrated APOBEC-mediated SARS-CoV-2 mutations. Dissemination of the SARS-CoV-2-infection and senescence was confirmed in extra-pulmonary sites (kidney and liver) of a COVID-19 patient.
Conclusions We demonstrate that in severe COVID-19, AT2 cells infected by SARS-CoV-2 exhibit senescence and a proinflammatory phenotype. In vitro, SARS-CoV-2 infection induces senescence and inflammation. Importantly, infected senescent cells may act as a source of SARS-CoV-2 mutagenesis mediated by APOBEC enzymes. Therefore, SARS-CoV-2-induced senescence may be an important molecular mechanism of severe COVID-19, disease persistence and mutagenesis.
Introduction
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19) that primarily affects the respiratory system. The clinical course of the patients ranges from asymptomatic to a life-threatening respiratory failure accompanied by a multi-systemic inflammatory disease [1, 2]. Systemic disease may occur through a virally-mediated inflammatory response that consists of a variety of proinflammatory cytokines and chemokines, including interleukin (IL)-1β, IL-6, IL-12, IFN-γ, TNF-α, CXCL-8, CXCL-10 and CCL-2 [3, 4]. The link between viral infection of cells and development of severe lung disease and systemic manifestations is still poorly understood. Viral infection results in the activation of complex innate and adaptive immune responses that are orchestrated sequentially, involving several cell types and inflammatory mediators [5, 6]. At the cellular level, intrinsic defence mechanisms are activated and outcomes range from complete recovery to cell death [7–11]. An “intermediate” and essential cellular state that is overlooked, due to lack of efficient methodological tools, is cellular senescence [12, 13].
Cellular senescence is a stress response mechanism that preserves homeostasis. Senescent cells are characterized by prolonged and generally irreversible cell-cycle arrest and resistance to apoptosis [12, 14]. Additionally, they also exhibit secretory features, collectively described as the senescence-associated secretory phenotype (SASP) [12]. SASP includes a variety of cytokines, chemokines, growth factors, proteases and other molecules, depending on the type of senescence and the cells involved [12, 15]. These inflammatory proteins are released in the extracellular space as soluble factors, transmembrane proteins following ectodomain shedding, or as molecules engulfed within extracellular vesicles [16–18]. Cellular senescence plays a key role in several lung diseases, involving the senescence of several cell types in the lung [19, 20]. Under physiological conditions, senescence is transiently activated and SASP mediates the recruitment of immune cells for senescent cell clearance. In addition, other SASP factors promote tissue regeneration and repair, overall ensuring cellular/tissue homeostasis. On the contrary, persistence of senescent cells exerts harmful properties promoting tissue dysfunction and the maintenance of chronic inflammation, via paracrine and systemic SASP [12, 15].
There is little published evidence linking viral infection to cellular senescence [21–26]. Given the implication of an inflammatory response in the progression of COVID-19 and the SASP secretion by senescent cells, we investigated whether cellular senescence occurs in COVID-19. Since individual markers are not adequate to unequivocally detect senescence, especially in vivo, as they may also be present in non senescent conditions, we followed in a clinical setting a detailed multi-marker algorithmic approach that we and others recently published and was approved by the senescence community [12, 27]. We were the first to provide evidence supporting the proof for senescence in COVID-19 infected lung tissue by applying this algorithm in a previous, preprint (bioRvix), version of the current manuscript [28], which subsequently has been confirmed by others [29]. This is the largest clinical study verifying viral-induced senescence by SARS-CoV-2 infection linked to a proinflammatory phenotype, which may contribute to acute and chronic clinical manifestations in severe COVID-19. Importantly, we also provide novel evidence for the generation of viral mutations promoted by SARS-CoV-2 persistence in senescent cells.
Materials and methods
Tissues
Formalin-fixed and paraffin-embedded autopsy lung tissue samples were obtained from 11 patients that died from severe COVID-19 (confirmed by RT-qPCR) (table 1 and Table S1A). Lung parenchyma displaying analogous lesions (atelectasia, fibrosis and infiltration by immune cells) from 43 age-matched patients with similar comorbidities from previously published cohort [30] and new cases, resected prior to the COVID-19 outbreak, were used for comparison with COVID-19 samples for all experiments detailed below (negative controls) (table 1 and Table S1B). A non-COVID cohort of 60 age-matched lung samples was used as biological negative controls for validation of the anti-SARS-CoV-2 (G2) monoclonal antibody (online supplement). Moreover, renal and liver tissue from one of the eleven COVID-19 cases and five non-COVID-19 patients were also available and analysed for the presence of the virus and cellular senescence. Clinical sample collection and their experimental use were approved by the Commission Cantonale D'éthique de la Recherche, University of Lausanne, Switzerland (2020–01257) and the Bio-Ethics Committee of the University of Athens Medical School, Greece.
Summary of clinical characteristics of patients
Cells and SARS-CoV-2 culture
SARS-CoV-2 (isolate 30–287, B.1.222 and B.1 strain (accession No MT459880.1) was obtained through culture in Vero E6 cells (ATCC® CRL-1586), from two different infected (COVID-19) patients. The virus was recovered from a nasopharyngeal swab, rinsed in 1 mLsaline and filtered twice through a 0.22 nm filter. Virus stock was prepared by infecting fully confluent Vero-E6 cells in DMEM, 10% fetal bovine serum (FBS) with antibiotics, at 37°C and 5% CO2. Virus stock was collected four days after inoculation, sequenced by NGS (online supplement) and the supernatant was frozen (−80°C) until use. Infections were carried out in 24-well plates, using SARS-CoV-2 at a 0.01 MOI. Cells were either fixed with 4% paraformaldehyde or lysed with NucleoZOL (Macherey-Nagel) 17 days post infection. Administration of the specific ataxia-telangiectasia-mutated (ATM) protein kinase inhibitor KU-55933 was carried out as described elsewhere [31]. Three independent experiments were performed. All manipulations were carried out in a Biosafety level 3 facility.
SARS-CoV-2 primary human lung alveosphere system
Formalin fixed and paraffin embedded sections from a three-dimensional alveosphere system infected with SARS-CoV-2 and corresponding controls (non-infected) were analysed. The alveospheres consist of self-renewing AT2 cells that express the well-established AT2 cell markers, HTII-280 and surfactant protein C as well as of some cells expressing the AT1 cell marker HTI 56 [32]. Alveospheres were gently opened using pipetting to allow infection and to avoid cellular stress [32].
Anti-SARS-COV-2 (G2) antibody generation
Mouse immunization and antibody collection, selection and specificity determination are described in the online supplement. Transcriptome analysis of hybridomas and amino acid determination of selected clones are also provided. Four clones, (479-S1, 480-S2, 481-S3 and 482-S4) are under patent application (Gorgoulis VG, Vassilakos D, Kastrinakis N. GR-patent application no: 22–0003846810).
RNA extraction and Reverse-Transcription real-time PCR (RT-qPCR) detection were performed as previously described (online supplement) [33].
Next generation sequencing (NGS)
NGS was performed as previously described (online supplement) [34].
Protein extraction and immunoblot analysis
Protein extraction and immunoblot analysis was performed as described elsewhere [31]. For details see online supplement. Horse Radish Peroxidase conjugated anti-goat, anti-mouse and anti-rabbit secondary antibodies (1:1000 dilution) (Cell Signaling) were used. Primary antibodies were: anti-APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) 3G/A3G (Abcam ab109727 and ab172694), anti–APOBEC3H (LSBio LS–C151868) and anti-GAPDH (Cell Signaling).
Immunocytochemistry and immunohistochemistry
Immunocytochemistry and immunohistochemistry were performed according to published protocols [34]. The following primary antibodies were applied overnight at 4°C: anti-SARS-CoV-2 (G2) monoclonal antibody (at a dilution 1:300), anti-SARS-CoV-2 monoclonal antibody (1A9 clone, Genetek), anti-ACE-2 (Abcam), anti-thyroid transcription factor (TTF)-1 (Dako), anti-Surfactant Protein B (SP-B) (1B9 clone, ZETA Corporation), anti-p16INK4A (Santa Cruz), anti-IL-1β (Abcam), anti-IL-6 (R&D systems), anti-phospho-histone (Ser139) H2AX (γΗ2ΑΧ) (Cell Signaling), anti-Ki67 (Abcam) and anti-p21Waf1/Cip1(Cell Signaling).
Senescence-SenTraGorTM (GL13) methodology
SenTraGorTM (trademark of GL13 compound) staining and double staining experiments were performed and evaluated as previously described [27, 35]. The mean percentage of SenTraGor positive alveolar cells in at least 10 high power fields (×400) per patient was quantified.
Electron microscopy
Representative area from hematoxylin and eosin-stained paraffin sections of the lung autopsy of COVID-19 patients and non COVID-19 controls were chosen under the light microscope and marked. Paraffin-embedded tissue was deparaffinized, rehydrated and fixed in 2.5% glutaraldehyde in PBS for 24 h and post-fixed in 1% aqueous osmium.
Bioinformatic analysis for identification of mutational signatures in the SARS-CoV-2 genome
As described in detail in online supplement, candidacy for APOBEC deamination of C→U mutational substitutions in the SARS-CoV-2 genome of strains available in GISAID database [36, 37] and in those obtained after cell culture infection, was examined and verified against experimentally validated APOBEC motifs (Figure S7).
Statistical analysis
Data are expressed as mean±sd. Comparisons were performed with unpaired non-parametric Mann-Whitney U test (comparison between two groups) or Kruskal Wallis test followed by Dunn's post-hoc analysis (comparison between three groups). The Wilcoxon paired non-parametric test was used to compare values of infected cells in vitro.
Results
Detection of SARS-CoV-2 in lung cells
In order to detect SARS-CoV-2 in lung tissue we developed monoclonal antibodies which react against receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 and identified a high affinity antibody (G2), the validity of which was recently verified [38] (Figure S1, Table S1). SARS-CoV-2 was detected (using both our “in house” G2 clone and a commercially available antibody from GeneTek) predominantly in alveolar type 2 (AT2) cells, which were identified by TTF-1 and Surfactant Protein B (SP-B) positivity, and in sparse inflammatory cells (alveolar and tissue macrophages) in all COVID-19 patients, ranging from <5 cells/4mm2 tissue to >50 cells/4mm2 tissue (figure 1a, Table S1A). Surface epithelial cells in small peripheral airways stained also positive in certain cases (Figure S2). SARS-CoV-2 infected AT2 cells were occasionally large and appeared isolated (denuded or syncytial) or clustered (hyperplasia), exhibiting a variety of topological distribution (figure 1a). These cells co-expressed the angiotensin-converting enzyme 2 (ACE2) receptor (figure 1b), supporting SARS-CoV-2 infection being mediated by the ACE2 receptor [39]. In addition, electron microscopy analysis in representative COVID-19 cases confirmed the presence of virus within AT2 cells (figure 1c) and high magnification revealed virions in the proximity of the endoplasmic reticulum, indicating their likely assembly and budding, as well as virions residing in cytoplasmic vesicles, implying their transfer and release into the extracellular space.
Detection of SARS-CoV-2 in lung cells from patients who died of COVID-19. a. Representative images of SARS-CoV-2 (G2 mab) immunohistochemistry (IHC) staining in COVID-19 autopsy lung tissue. Competition with anti-peptide (S protein) showing specificity of the G2 immunostaining. Verification of G2 immunopositivity with a commercial anti-S antibody (clone 1A9 clone, Genetek). Representative negative control immunostaining in non-COVID-19 lung tissues. Graph shows quantification of SARS-CoV-2 staining in the lung samples. b. Detection of SARS-CoV-2 in AT2 cells (confirmed by positive TTF-1 and Surfactant Protein B (SP-B) staining) and in ACE2-expressing cells. Double immunostaining is shown for SARS-CoV-2 and TTF-1. C. Detection of SARS-CoV-2 by transmission electron microscopy (EM) in a representative COVID-19 patient. Presence of SARS-CoV-2 within AT2 cells (i, ii) and of virions in the proximity of the endoplasmic reticulum (iii, iv) as well as in cytoplasmic vesicles (iii, v-vi). Corresponding scale bars are depicted. AT2: alveolar type 2 cells; ACE2: Angiotensin-Converting Enzyme 2; ****Statistical significance: p<0.0001.
Senescence in SARS-CoV-2 infected cells
A proportion of SARS-CoV-2 infected AT2 cells (range 8–21%) displayed a senescent phenotype, with positive staining for SenTraGor™ and p16INK4A (figure 2a-c, Figure S3) [12, 29, 30]. SenTraGor™ is an established marker of senescence and has the unique property of detecting senescent cells in any setting, including archival material [formalin-fixed paraffin embedded (FFPE) tissue)] [12, 27, 35]. By contrast, lung tissues exhibiting fibrosis, atelectasia and inflammatory infiltrates from age-matched non-COVID-19 cases with similar co-morbidities (table 1 and Table S1), including a separate control cohort of patients with acute pneumonia (aspiration pneumonia), showed significantly lower senescence (range 1–2%, p<0.0001) (figure 2a-c, Figure S4), suggesting that SARS-CoV-2 infection may induce senescence.
Senescence in SARS-CoV-2 infected cells. a. Representative images of SARS-CoV-2, SenTraGor™ (senescence) and ACE-2 staining in serial sections of COVID-19 lung tissue (1). Double-immunostaining analysis for SARS-CoV-2, SenTraGor™ (senescence), p16INK4A, TTF-1 and Surfactant Protein B (SP-B) in COVID-19 lung tissue (2). b. Representative results from serial staining for SARS-CoV-2, SenTraGor™ (senescence) and ACE-2, and double-staining for SARS-CoV-2 and p16INK4A in non-COVID-19 lung tissue. C. Graphs showing increased expression of SenTraGor™ and p16INK4A in COVID-19 lung tissue. Corresponding scale bars are depicted. ****Statistical significance: p<0.0001.
In order to demonstrate that SARS-CoV-2 can induce cellular senescence, we infected Vero cells with a viral strain isolated from a COVID-19 patient. Vero-E6 cells is an established cellular system for viral propagation and studies, as apart from their high infectivity to SARS-CoV-2 they are among the few cell lines demonstrating SARS-CoV-2-mediated cytopathic effects, an essential aspect in diagnostics [40, 41]. Infection was carried out at a low multiplicity of infection (MOI) to mimic natural coronavirus infection [42]. In line with our hypothesis, the infected cells following an initial surge of cell death reached an equilibrium demonstrating clear evidence of senescence (increased SenTraGor positivity, increased p21WAF1/CIP1 and reduced Ki-67), as compared to the non-infected control cells at 17 days post infection (figure 3 and Fig. S5A). As Vero cells lack p16INK4A[43], the most likely trigger of senescence is DNA damage, as previously reported [12, 13]. DNA damage measured by γ-H2AX immunostaining was evident in SARS-CoV-2 infected cells (figure 3A4). Treatment of infected cells with a selective ATM inhibitor (KU-55933) dramatically decreased senescence assessed by SenTraGor staining (Figure S6). ATM has been previously reported as a key driver of NF-κB-dependent DNA-damage-induced senescence [44]. It appears that genotoxic stress results from a vicious cycle imposed by the virus in host cells as it hijacks most intracellular protein machineries [11, 45]. Likewise, cellular senescence was identified in infected alveolar cells/alveospheres of an established primary lung alveolar three-dimensional model, confirming the in vivo observations and the findings in Vero cells (figure 3b, Figure S5B).
SARS-CoV-2 induced senescence in Vero cells and alveospheres. a. SARS-CoV-2 presence (1), senescence induction (2), cellular proliferation (3) and DNA damage activation (4), with corresponding quantitative histograms, in Vero cells with and without SARS-CoV-2 infection (17 days post infection). b. Senescence induction (assessed by SenTraGor and p21WAF1/CIP1 staining) in alveolar cells (alveospheres) of a primary lung alveolar three-dimensional model following SARS-CoV-2 infection (black arrows) in comparison to non senescent ones (yellow arrows). Absence of senescence is clearly evident in non infected alveolar cells. C. Graph showing induction of SASP related cytokines following SARS-CoV-2 infection of Vero cells (mRNA expression–normalized against β2-microglobulin mRNA expression) **Statistical significance: p<0.01.
Senescence-associated secretory phenotype
We found very high expression of both IL-1β and IL-6 by senescent AT2 cells in the lungs of COVID-19 patients, while in the non COVID-19 control cases, including samples from patients with acute pneumonia, expression was very low in the few senescent AT2 cells detected in lung tissue with fibrosis, atelectasia and inflammatory infiltrates (p<0.0001) (figure 4a-c, Table S1). As both cytokines are key components of the inflammatory response to SARS-CoV-2, our findings may imply this pattern of inflammation may be due to the SASP as a result of cellular senescence in COVID-19 patients. Likewise, SARS-CoV-2 senescent Vero cells displayed expression of SASP-related cytokines, as assessed by our algorithmic assessment of senescence, supporting our in vivo findings (figure 3c).
Senescence associated secretory phenotype (SASP) in COVID-19 lung tissues. a. Representative staining results (at low and high magnification) of SenTraGorTM, IL-6, IL-1β and TTF-1 and/or Surfactant Protein B (SP-B) in corresponding serial sections (1) and as double immunostaining analysis (2) of COVID-19 lung tissue. Original magnification: 400×. b. Representative staining results showing absence or minimal levels of SenTraGor™, IL-6 and IL-1β in age-matched non-COVID-19 control samples. Corresponding scale bars are depicted. IL-1 β: Interleukin 1β; IL-6: Interleukin-6; **** Statistical significance: p<0.001.
SARS CoV-2 extra pulmonary dissemination
SARS-CoV-2 immunoreactivity was additionally detected by applying both our “in house” G2 clone and a commercial one from GeneTek in serial sections of available kidney and liver tissue from one COVID-19 patient (figure 5a). A clear cytoplasmic signal was evident in a number of renal tubules as well as in areas of the liver parenchyma. In contrast, in tissues from non-COVID-19 cases the signal was absent, as expected (figure 5a). Interestingly, by serial section analysis, areas harboring the virus also exhibited senescence (SenTraGor positivity), implying occurrence of SARS-CoV-2 induced senescence in extra pulmonary sites (figure 5b).
SARS-CoV-2 dissemination in extra pulmonary sites. a. Representative staining results of SARS-CoV-2 in the kidney and the liver of a COVID-19 patient. Absence of the immunohistochemical signal is evident in the corresponding tissues of a representative non COVID-19 case. b. Serial section analysis the representative COVID-19 case (see panel A) depicting concurrent positivity for SARS-CoV-2 and the senescence marker (SenTraGor™) in cells of the liver and kidney.
Infected senescent cells and SARS-CoV-2 mutagenesis
We recently hypothesized that senescent cells could represent a source of SARS-CoV-2 quasispecies generation [46]. This assumption was based on the fact that the apoptotic tolerant nature of senescent cells allows the virus to be hosted for longer periods compared to other cells with higher cell turnover [47, 48], thus rendering the SARS-CoV-2 virome more susceptible to host mediated editing. APOBEC enzymes are well known to participate in viral genome editing [12, 33, 49]. Notably, we recently showed these enzymes to be highly expressed in senescent cells [49].
Updating our initial bioinformatic investigation [46], we screened ∼ 4,500,000 SARS-CoV-2 strains from the GISAID database [36, 37]. This analysis mostly indicated that the predominant signature of SARS-CoV-2 variants is APOBEC-mediated (figure 6), in line with recent literature [50]. To validate our hypothesis, we proceeded to cell culture infection with two strains isolated from patients and obtained progenies from senescent cells (figure 6b). Next generation sequencing (NGS) analysis demonstrated that the collected viral genomes acquired mutations. Importantly, the predominant pattern was APOBEC-driven, in agreement with the bioinformatic analysis of the GISAID database. Moreover, and in line with our presumption and our previous study [49], the SARS-CoV-2 induced senescent cells demonstrated significantly higher mRNA and protein levels of the APOBEC enzymes, particularly G and H (RNA editing cytoplasmic variants), which are reported to play a pivotal role in viral RNA editing [12, 33, 49]. Other mutation signatures were also found both in the cell culture isolated virome genomes and in the strains from GISAID, and maybe the outcome of the oxidative stress present in the senescent cells [51].
Infected senescent cells as a putative source for SARS-CoV-2 quasi-species variant generation. a. Schematic layout presenting the identification of APOBEC-mediated mutations in the genome of 4,672,296 SARS-CoV-2 strains available in GISAID database. For algorithmic assessment see Figure S7 and corresponding section in Materials & Methods. b. Identification of SARS-CoV-2 mutations acquired in senescent cells in two strains (i: B.1.222 and ii: B.1) isolated from two different patients, following cell culture infection for the depicted time periods. Yellow bars show the frequency of C→U substitutions observed in the GISAID database in the NGS read counts (green pileups) and red bars depict the C→U frequency of mutation observed, post infection of cells in the current study, in the genome of the two employed SARS-CoV-2 strains. These sites are ranked as highest relatively to the C→U counts as observed from the GISAID database. The genomic sequences flanking each C→U are also presented, while their co-ordinates are shown relative to the superimposed SARS-CoV-2 genome. NGS (Next Generation Sequencing) reads were confirmed in triplicate reads.
Discussion
We have demonstrated the presence of SARS-CoV-2 in AT2 cells of patients who died from COVID-19 using our novel and a commercial anti-SARS-CoV-2 antibody against its spike protein and by electron microscopy. We have shown for the first time following the senescence detecting algorithm in vivo that a proportion of SARS-CoV-2-infected AT2 cells exhibit features of cellular senescence (as demonstrated by significantly increased staining with the novel senescence marker SenTraGor™ and also of p16INK4A) [12, 27, 35]. In this respect, our investigation encompassing serial section analysis and co-staining is the first to demonstrate in vivo which cells are truly senescent and infected. The finding that in age-matched non-COVID-19 controls, including lung tissue samples from a cohort of patients with acute pneumonia, the percentage of senescent cells was much lower (1–2%) than that of the COVID-19 patients (8–21%), strongly indicates that SARS-CoV-2 triggers senescence. To confirm that SARS-CoV-2 induces cellular senescence per se, we demonstrated that infection of epithelial cells with SARS-CoV-2 virus (B.1.222 strain) in vitro increased SenTraGor™ staining and induced DNA damage, measured by increased γ-H2AX expression. Importantly, we found that inhibition of ATM (ataxia-telangiectasia-mutated protein kinase), an apical orchestrator of the DNA damage response pathway, dramatically reduced senescence in infected samples, suggesting SARS-CoV-2 induced senescence being mediated by DNA damage and activation of the DNA damage response pathway. In line with our in vivo and in vitro findings, SARS-CoV-2 induced senescence was identified in a previously established primary lung alveolar three-dimensional model [32]. The fact that alveospheres underwent gentle mechanical manipulations is unlikely to have contributed to the occurrence of senescence, given that such interventions are frequently applied on 3D organoid systems and no senescence phenotypes have been reported nor are evident from the post-infection proliferation rates and culture [52].
We also demonstrated that cells infected with SARS-CoV-2 exhibit high expression of IL-1β and IL-6, both components of the SASP and implicated in systemic features of severe COVID-19 [3, 4]. Therefore, our in vitro and in vivo findings suggest that SARS-CoV-2 attaches to AT2 cells via ACE2 to infect these cells and through activation of DNA damage response signalling roots, induces cellular senescence and associated proinflammatory phenotype (SASP). In line with our observations, Lee et al. recently demonstrated that SARS-CoV-2-induced senescence exhibited enhanced γ-H2AX DNA damage foci that were abrogated by the reactive oxygen species (ROS) scavenger N-acetylcysteine (NAC) [29], suggesting implication of ROS and oxidative stress (Figure S8). Moreover, infected cells displayed SASP-positive senescence mediated by p53 and activation of cGAS/STING signaling pathway [29]. They also demonstrated that senolytics reduced COVID-19 lung disease and inflammation in infected hamsters and mice [29], further suggesting SASP as an outcome of viral induced senescence. Camell et al. also reported that treatment with senolytics of old mice infected with a SARS-CoV-2–related mouse β-coronavirus reduced inflammation [53]. Based on our findings and recent evidence a potential mechanistic scenario could be the following (Figure S8): upon entrance, the virus highjacks certain energy consuming functions related to RNA processing, translation and the ER [11]. Thus, an energy shift occurs resulting in increased ROS production, oxidative stress and DNA damage/DNA Damage Response (DDR) pathway activation. In addition, DDR and the subsequent cell cycle arrest may be driven by an interaction of Coronavirus nsp13 protein and DNA polymerase δ [54] and SASP via the cGAS/STING and other DNA damage dependent pathways (Figure S8). Therefore, our previous work showing that DNA damage/DDR is a potent inducer of senescence, signifies its importance in the establishment of the viral induced senescence phenotype [31].
Senescent cells are in a state of cell cycle arrest but remain metabolically active and secrete a typical profile of inflammatory proteins known as the SASP. SASP components include the proinflammatory cytokines IL-1β and IL-6, which are elevated in plasma of COVID-19 patients that have acute respiratory distress syndrome (ARDS) or systemic inflammatory features. We demonstrate that in COVID-19, senescence alters the properties and function of AT2 cells, producing cytokines that could be released into the systemic circulation, amplifying and perpetuating chronic inflammation. Of course, the availability of tissue samples and the absence of more refined COVID-19 in vitro and in vivo models make it difficult to draw robust conclusions. It is likely that SARS-CoV-2 spreads from epithelial cells in the lower airways to infect AT2 cells in the alveolar walls, which express ACE2, and causes local senescence and inflammation via SASP in the lung. In line with this notion, it has been previously reported in vitro that cellular senescence has also a role as a natural antiviral defence mechanism and that SASP acts as the major contributor of this response, activating and recruiting the immune system to clear out the infection [21, 22]. Subsequently, other cellular compartments (stem/progenitor, endothelial and inflammatory cells) may be affected via different cell and non-cell autonomous mechanisms, leading to aberrant immune responses, chronic inflammation, tissue dysfunction and/or fibrosis and eventually to lung damage and failure [26, 55]. In particular, senescence in stem/progenitor cells impairs lung regenerative capacity. Immunosenescence results in elevated neutrophils-to-lymphocytes ratio (NLR) and high IL-6 production. In addition, IL-1, IL-6, CXCL8 and TNF-α might be crucial for maintaining chronic inflammation and TGF-β, PAI-1, and MMPs might favour fibrosis via the expression of fibrotic genes (ACTA2, COL1A1, COL1A2, and FN1) in the surrounding microenvironment. The virus may then enter the circulation and senescence may subsequently spread systemically to affect other organs (figure 5), leading to multi-organ failure/multi-morbidity and death in the acute phase [1, 25], or leading to post-acute sequelae of COVID-19 (PASC or long-COVID), an evolving syndrome with long-term complications [26, 56]. In this manner, senescence could also predispose to disease severity upon re-infection by SARS-CoV-2 or infection by other viruses and even negatively impact vaccine efficacy, similarly to what has been proposed for aged individuals [26].
An additional implication relates to the prolonged survival of senescent cells that are infected with the virus, as senescent cells are resistant to apoptosis and clearance by efferocytosis [12, 14]. Such a context can provide an extended time “window” for virus replication, therefore exposing its genome to host-mediated editing. Along this vein we showed that the senescent cell compartment acts as a “fertile” environment for mutational evolution of the virus, as it is more susceptible to APOBEC-3G and 3H RNA editing (figure 6) [33, 49]. Of course, other mechanisms extending survival and replication of SARS-CoV-2 may take place. Coronaviruses code for an important multifunctional enzyme termed papain-like protease (PLP) that exerts intrinsic deubiquitinating and deISGylating activities. The latter is related to interferon-stimulated genes (ISG) up-regulation and can serve to antagonize the host's immune response that would otherwise hinder infection [57]. In this manner increased and prolonged viral replication and exposure of its genome to host mediated editing and quasispecies generation could also occur. Taking into account that interferons are known senescence inducers, the picture gets even more complicated, rendering difficult to discriminate IGS-mediated mechanisms from that involving virally-induced senescence [58]. Of course, additional studies are needed and the guideline algorithmic approach for in vivo senescence assessment, followed in the current work, is anticipated to clarify in the future a lot of unanswered questions in clinical settings [12, 27]. Notably, in the recently identified Omicron (B.1.1.529) SARS-CoV-2 variant/strain the APOBEC signature was identified as the predominant mutational profile (Figure S9). To better understand the role of APOBEC in the generation of viral mutations, inhibition of APOBEC proteins by RNA interference would be valuable. However, this is very challenging, since APOBEC suppression induces DNA damage and sensitises cells to stress induced death [59–61]. Interestingly and in line with our concept linking cellular senescence with viral mutagenesis, by conducting a bioinformatic analysis, we found that viruses bypassing cellular senescence (oncogenic) exhibit a significant lower mutation burden compared to viruses inducing senescence (Figure S10). Given that mutations can affect the effectiveness of vaccines, our findings could imply a potential implication of viral induced senescence in vaccination strategies. Further studies are required towards this root.
Despite the fact that up to the present this is the largest clinical study demonstrating virally-induced senescence by SARS-CoV-2, a limitation of our study is the relatively small sample size of the COVID-19 lung autopsies, due to difficulty of accessing this rare material, given that the autopsies are limited and were performed only in the initial phase of the COVID-19 outbreak in order to identify the pathological basis of this new entity. Another limitation due to the nature of the disease is that pathological features, such as cellular senescence in lungs, can only be investigated in cadaverous material, which represents the most severe outcome of the spectrum of COVID-19 clinical manifestations. Therefore, evaluation of senescence in less severe conditions is not currently feasible. Moreover, in our study we examined 18 cases of aspiration pneumonia, which displayed similar degrees of senescence to other non-infected control cases. It would be interesting to extend these studies in other types of pneumonia. However, material from patients with an acute but not COVID-19 related inflammatory response is rare, given that these patients commonly recover following treatment and in the infrequent cases of fatal outcome autopsy is rarely performed.
Overall, SARS-CoV-2 induced senescence justifies the therapeutic application of senotherapeutics for the treatment of COVID-19 patients and possibly for long-COVID syndromes [56]. In line, quercetin was recently tested against standard care in two randomized clinical trials (ClinicalTrials.gov identifiers NCT04578158 and NCT04861298) in patients with confirmed SARS-CoV-2 infection and mild COVID-19-associated symptoms. In both trials, senolytic intervention led to clinical improvement [62, 63]. Overall, these data imply the beneficial role of senolytics in significantly improving the outcome during SARS-CoV-2 infection.
Acknowledgements
We would like to thank Konstantinos Ntostoglou for his valuable help in preparing the material. We acknowledge support in RNA sequencing by the Greek Research Infrastructure for Personalised Medicine (pMED-GR) (MIS 5002802), which is implemented under the Action “Reinforcement of the Research and Innovation Infrastructure”, funded by the Operational Programme “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014–2020) and co-financed by Greece and the European Union (European Regional Development Fund). This work was supported by the: National Public Investment Program of the Ministry of Development and Investment / General Secretariat for Research and Technology, in the framework of the Flagship Initiative to address SARS-CoV-2 (2020ΣΕ01300001); Horizon 2020 Marie Sklodowska-Curie training program no. 722729 (SYNTRAIN); Welfare Foundation for Social & Cultural Sciences, Athens, Greece (KIKPE); H. Pappas donation; Hellenic Foundation for Research and Innovation (HFRI) grants no. 775 and 3782 and NKUA-SARG grant 70/3/8916. KP is financially supported by National Heart & Lung Institute – Imperial College London and a British Heart Foundation Project Grant (PG/19/75/34686).
Footnotes
Author contributions: K.E., D.V.: immunocytochemistry, SenTraGorTM staining; K.E., P.G.F., C.K., L.DL, V.G.G.: histopathological examination and evaluation of the stainings; M.D., I.K.: cell culture and viral manipulations; B.S.: 3D alveosphere system and viral manipulations; A.K., An.P.: RNA and protein extraction, RT-qPCR, NGS, western blot; N.L., S.H.: electron microscopy; O.H., Ai.P.: Statistical analysis, bioinformatic analysis and graphs generation; D.Va.: anti-SARS-CoV-2 antibody generation; K.E., A.G.T, K.P., S.T., V.G.G.: data analysis and interpretation, manuscript preparation; S.T., Ar. P., G.B., P.J.B., V.G.G.: experimental design, guidance, manuscript writing with input for all co-authors.
Funding: This work was supported by the: National Public Investment Program of the Ministry of Development and Investment / General Secretariat for Research and Technology, in the framework of the Flagship Initiative to address SARS-CoV-2 (2020ΣΕ01300001); Horizon 2020 Marie Sklodowska-Curie training program no. 722729 (SYNTRAIN); Welfare Foundation for Social & Cultural Sciences, Athens, Greece (KIKPE); H. Pappas donation; Hellenic Foundation for Research and Innovation (HFRI) grants no. 775 and 3782 and NKUA-SARG grant 70/3/8916. National Public Investment Program of the Ministry of Development and Investment / General Secretariat for Research and Technology, in the framework of the Flagship Initiative to address SARS-CoV-2 ; Grant: 2020ΣΕ01300001; H2020 Marie Skłodowska-Curie Actions; DOI: http://dx.doi.org/10.13039/100010665; Grant: 722729 (SYNTRAIN); Welfare Foundation for Social & Cultural Sciences, Athens, Greece (KIKPE); Hellenic Foundation for Research and Innovation; DOI: http://dx.doi.org/10.13039/501100013209; Grant: 3782, 70/3/8916, 775, NKUA-SARG grant; H. Pappas.
Conflict of Interest: PB receives research funding from AstraZeneca and Boehringer-Ingelheim and is a scientific advisor to AstraZeneca, Boehringer-Ingelheim, Epi-Endo, Novartis, Pieris and Teva. The other authors wish to declare no conflict of interest.
- Received August 27, 2021.
- Accepted December 24, 2021.
- Copyright ©The authors 2022.
This version is distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0. For commercial reproduction rights and permissions contact permissions{at}ersnet.org