- Open Access
Biomarkers of cell damage, neutrophil and macrophage activation associated with in-hospital mortality in geriatric COVID-19 patients
Immunity & Ageing volume 19, Article number: 65 (2022)
The risk for symptomatic COVID-19 requiring hospitalization is higher in the older population. The course of the disease in hospitalised older patients may show significant variation, from mild to severe illness, ultimately leading to death in the most critical cases. The analysis of circulating biomolecules involved in mechanisms of inflammation, cell damage and innate immunity could lead to identify new biomarkers of COVID-19 severity, aimed to improve the clinical management of subjects at higher risk of severe outcomes.
In a cohort of COVID-19 geriatric patients (n= 156) who required hospitalization we analysed, on-admission, a series of circulating biomarkers related to neutrophil activation (neutrophil elastase, LL-37), macrophage activation (sCD163) and cell damage (nuclear cfDNA, mithocondrial cfDNA and nuclear cfDNA integrity). The above reported biomarkers were tested for their association with in-hospital mortality and with clinical, inflammatory and routine hematological parameters. Aim of the study was to unravel prognostic parameters for risk stratification of COVID-19 patients.
Lower n-cfDNA integrity, higher neutrophil elastase and higher sCD163 levels were significantly associated with an increased risk of in-hospital decease. Median (IQR) values observed in discharged vs. deceased patients were: 0.50 (0.30-0.72) vs. 0.33 (0.22-0.62) for n-cfDNA integrity; 94.0 (47.7-154.0) ng/ml vs. 115.7 (84.2-212.7) ng/ml for neutrophil elastase; 614.0 (370.0-821.0) ng/ml vs. 787.0 (560.0-1304.0) ng/ml for sCD163. The analysis of survival curves in patients stratified for tertiles of each biomarker showed that patients with n-cfDNA integrity < 0.32 or sCD163 in the range 492-811 ng/ml had higher risk of in-hospital decease than, respectively, patients with higher n-cfDNA integrity or lower sCD163. These associations were further confirmed in multivariate models adjusted for age, sex and outcome-related clinical variables. In these models also high levels of neutrophil elastase (>150 ng/ml) appeared to be independent predictor of in-hospital death. An additional analysis of neutrophil elastase in patients stratified for n-cfDNA integrity levels was conducted to better describe the association of the studied parameters with the outcome.
On the whole, biomarkers of cell-free DNA integrity, neutrophil and macrophage activation might provide a valuable contribution to identify geriatric patients with high risk of COVID-19 in-hospital mortality.
The coronavirus disease 2019 (COVID-19) global pandemic caused by the SARS coronavirus 2 (SARS-CoV-2) continues to pose significant challenges to global safety in public health. The pandemic, declared by the World Health Organization (WHO) on March 11, 2020, caused 375.465.073 total cases and 6.110.611 fatalities worldwide, as of 25/03/2022 . Older patients with COVID-19, in particular, present higher risk of death compared to younger age groups. However, the severity of the disease and the clinical trajectories of the older infected cases show significant variation, from mild illness to severe disease requiring hospitalization and, in the most critic cases, the admission to intensive care unit . This variability suggests the importance to find biomarkers of COVID-19 severity that can be used as personalised predictors, especially in the setting of older patients.
An exacerbate inflammatory response to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) with the overproduction of many inflammatory cytokines, known as “cytokine storm” and often displaying the features of macrophage activation syndrome (MAS), is thought to be a major cause of disease severity and death in patients with COVID-19 pneumonia [3, 4]. Among the circulating biomolecules useful to unravel prognostic parameters for risk stratification of COVID-19 patients with the most severe disease, a number of biomarkers of neutrophil and macrophage activation were proposed . Soluble (s)CD163 is the shed form of the corresponding membrane marker of monocyte/macrophage. As a result of the shedding, during inflammation and macrophage activation, the extracellular portion of CD163 circulates in the blood as sCD163 and its levels significantly increase during inflammatory responses, including SARS-CoV-2 infection . In recent years, the CD163 receptor has also been reported to directly bind some pathogenic bacteria and also viruses .
Neutrophil elastase (NE), a serine proteases, and extracellular DNA are biomarkers of neutrophil extracellular traps (NET) formation . The release of NETs, i.e. complexes created by cationic effectors, histones, and decondensed nuclear DNA (e.g., after histone citrullination), has been identified as one of the effector mechanisms utilized by neutrophils to kill microbes [9, 10]. Increasing evidence suggested that NETs can contribute to the inflammatory storm that leads to respiratory failure in many patients with COVID-19 [11,12,13].
In NETs, chromatin fibers are associated to chromatin proteins and other proteins released by neutrophils, including, among others, NE and cathelicidin (LL-37) . The release of NETs, that is, NETosis, plays a major role in entrapping and eliminating pathogens, thus preventing their dissemination, including viruses, as confirmed for influenza [15, 16]. However, the beneficial role of NETs in the innate immune response coexists with its detrimental effects, such as the promotion of tissue damage, thrombosis, and autoimmunity [17, 18]. The inhibition of neutrophils and NETs was proven protective in various models of influenza-associated Acute Respiratory Distress Syndrome (ARDS) [19, 20]. Cell-free DNA (cfDNA) is present at low concentration in plasma in healthy conditions [21, 22]. Circulating cell-free DNA (cfDNA) mostly derives from cells that underwent apoptosis or necrosis and from the formation of NETs [23,24,25]. The relative contribution of the different processes by which cells release their nuclear DNA into body fluids can vary among different individuals and disease conditions, and this variability affects the value of cfDNA integrity, i.e. the ratio of longer to shorter DNA fragments in cfDNA . Recent findings showed that in COVID-19, plasma cfDNA amount correlates with the WHO ordinal scale for disease progression . Finally, cfDNA can be considered a powerful modulator of systemic inflammation, especially in older patients affected by COVID-19 .
In this work we aimed to evaluate the clinical significance and interlink between parameters of cell damages (n-cfDNA, mt-cfDNA, n-cfDNA integrity), neutrophil (NE and LL-37) and macrophage (sCD163) activation, in order to assess their prognostic value for risk stratification of COVID-19 patients with the most severe disease. Considering a cohort of geriatric COVID-19 patients admitted at the INRCA Hospital (Ancona), we analysed plasma/serum samples collected on admission, and verified the association of the selected biomarkers with survival and with clinical parameters.
Clinical and laboratory characteristics of the patients in relation with in-hospital mortality
The clinical characteristics of the enrolled COVID-19 patients, grouped according to the outcome of hospital stay, are reported in Table 1. Patients deceased during hospital stay had higher median age and increased prevalence of stroke, chronic obstructive pulmonary disease (COPD) and chronic kidney disease (CKD). In addition, a significantly higher proportion of CFS (Clinical Frailty Scale) scores ≥ 8 was observed among deceased patients (Table 1). As expected, COVID-19 patients who died in hospital also had, at the admission, increased number of white blood cells (WBC), neutrophil % and neutrophil to lymphocyte ratio (NLR) and decreased counts of the other white blood cells. Among inflammatory factors and cytokines, CRP, IL-6, and IL-10 were significantly higher in patients who deceased during hospital stay. The proportion of plasma samples positive for SARS-CoV-2 RNA was about five-fold higher in patients who deceased compared to patients who survived (21/49 vs. 9/107, p<0.001).
Association of markers of cell damage, neutrophil, and macrophage activation with in-hospital mortality, and comparison of n-cfDNA parameters between infected and uninfected patients
As shown in Table 2, patients who deceased in-hospital had lower n-cfDNA integrity in plasma and higher levels of plasma NE and serum sCD163 with respect to discharged patients. Other tested biomarkers, including plasma concentration of n-cfDNA (Alu 247, Alu 115), mt-cfDNA and LL-37 were not significantly different between deceased and discharged patients. A ROC curve analysis found the following AUC (area under curve) (±std.err.) and p-values for the three biomarkers significantly associated with in-hospital death: for Alu 247/115, (n-cfDNA integrity) 0.6082 (±0.05), p= 0.033; for NE, 0.6409 (±0.046), p = 0.002; for sCD163, 0.6457 (±0.0491), p=0.003.
A comparison of n-cfDNA parameters in non-COVID-19 patients (Supplementary Tables S1 and S2) with those observed in COVID-19 patients found significantly higher values for all the three parameters (Alu 247, Alu 115, Alu 247/115) in COVID-19 patients (considered as a whole, or as subgroups of “discharged” or “deceased” patients) with respect to non-COVID-19 patients, with the partial exception of the Alu 247/115 parameter (not significantly different between non-COVID-19 and deceased COVID-19 patients).
Distribution of plasma concentration of NE in patients with different outcomes and different levels of cfDNA integrity
In order to further evaluate the relationship between the studied biomarkers and the outcome, the association of the NE level with the outcome was tested after stratifying for different values of cfDNA integrity (Fig. 1). Notably, in patients with high cfDNA integrity (Alu247/115 > 0.625, tertile 3), median concentration of NE was nearly three times higher in deceased compared to survived patients [median NE (IQR): 347.9 (203.1-463.9) ng/mL vs. 123.1 (59.6-190.4) ng/mL, p=0.001]. A ROC curve analysis conducted using the same kind of patients’ stratification (based on tertiles of cfDNA integrity) showed that the predictive value of NE for in-hospital death was higher for patients with highest n-cfDNA integrity (AUC ± std.err. for NE, in patients with n-cfDNA integrity in tertile 3, was 0.8354 ± 0.0623, p<0.001) than for patients with lower cf-DNA integrity (AUC ± std.err. 0.6627 ± 0.0758, p=0.032 or 0.5838 ± 0.0859, p=0.329, for NE in patients with n-cfDNA integrity in tertile 1 or in tertile 2 respectively).
Association of markers of cell damage, neutrophil and macrophage activation with SARS-CoV-2 RNAemia
Table 3 shows that patients with detectable presence of SARS-COV-2 RNA in plasma had an increased abundance of n-cfDNA (Alu115), mt-cfDNA and NE. Other parameters, including integrity index of n-cfDNA (Alu247/115), sCD163 and LL-37 were not associated with SARS-CoV-2 RNAemia.
Evaluation of markers of cell damage, neutrophil, and macrophage activation as predictors of in-hospital survival
The three biomarkers showing differences between survived and deceased patients (Alu247/115, NE, sCD163) were then tested for their association with survival. To this purpose, each variable was categorized into tertiles. The Kaplan–Meier estimator was used as univariate analysis to estimate the overall survival associated with different levels of the three biomarkers. The results (Fig. 2A) show that patients with the lowest cfDNA integrity (Alu247/115 ≤ 0.32) have the lowest mean survival time compared to patients with average or higher n-cfDNA integrity: survival time (mean±sd) was 22.1 ± 1.9 days for tertile 1 (Alu247/115 ≤ 0.321); 32.2 ± 3.1 days for tertile 2 (0.321<Alu247/115≤0.625); 40.8 ± 4.3 days for tertile 3 (Alu247/115 > 0.625). On the other hand, patients with low levels of sCD163 (sCD163 < 491) show a longer survival with respect to the two subgroups with intermediate or high levels: 37.6 ± 3.0 days for tertile 1 (sCD163 ≤ 491 ng/ml); 27.3 ± 3.4 days for tertile 2 (491 ng/ml < sCD163 ≤ 811 ng/ml); 29.2 ± 3.8 days for tertile 3 (sCD163 > 811 ng/ml). When classifying patients according to NE, no significant differences were observed: survival time was 34.9±3.4 days for tertile 1 (NE ≤ 74.0 ng/ml); 31.3 ± 3.1 days for tertile 2 (74 ng/ml < NE ≤ 146.7 ng/ml); 31.6 ± 3.8 days for tertile 3 (NE > 146.7 ng/ml).
The three parameters (Alu 247/115, NE and sCD163) were also simultaneously tested for their association with survival in multivariate Cox regression models (Table 4) including sex and age (model 1), or sex, age, comorbidity and frailty (model 2), or all the previous variables and SARS-CoV-2 RNAemia (model 3). cfDNA integrity (Alu 247/115) was associated with the outcome in all the tested multivariate models. In particular, a lower cfDNA integrity index was significantly associated to an increased risk of in-hospital death in the three multivariate models. Higher levels of NE were associated with a higher risk of in-hospital death in multivariate models 2 and 3 (Table 4). Concerning sCD163, plasma concentration in the second tertile was associated with a higher risk of death in all the tested multivariate models (Table 4). The three tested biomarkers maintained their significant association with the outcome if pharmacological treatments were added to the multivariate Cox regression model, while two of them (Alu 247/115 and NE) maintained their significance even when oxygen therapy/ventilation was included as an additional confounding variable (Supplementary Table S3). A graphical representation of the adjusted survival curves for the three different biomarkers, based on the Cox multivariate model 3, is shown in Fig. 2 B.
Correlation of markers of cell damage, neutrophil and macrophage activation with other laboratory parameters
Spearman rank correlation was tested to explore reciprocal associations among different plasma biomarkers of cell damage, neutrophil and macrophage activation (supplementary Table S4) and between these markers and other laboratory parameters (supplementary Table S5). A positive correlation was observed between the plasma concentration of n-cfDNA (Alu 115, Alu 247) and mt-cfDNA (MT-CO3), while the n-cfDNA integrity (Alu 247/115) in plasma was positively correlated with both parameters of n-cfDNA concentration but not with mt-cfDNA plasma concentration (supplementary Table S4). NE was positively correlated with all the cfDNA parameters (Alu 115, Alu 247, mt-cfDNA and Alu 247/115), sCD163 was positively correlated with one of the parameters of n-cfDNA concentration (Alu 115) and with NE, while LL-37 was not correlated with any of the above cited parameters (supplementary Table S4). Scatterplots with the individual data points of the correlation between NE levels and n-cfDNA parameters are shown in Supplementary Figures.
When analysing the relationship between cytokines, or other inflammatory markers, and the studied biomarkers, it resulted that both parameters of n-cfDNA abundance were negatively correlated with serum IFN-α, and one of them (Alu 247) was also positively correlated with serum IL-6 (supplementary Table S5). mt-cfDNA in plasma was positively correlated with serum IL-10. NE was positively correlated with the tested inflammatory parameters with exception of serum IFN-α (negatively correlated), fibrinogen and TNF-α (not correlated). Finally, LL-37 and sCD163 were not significantly correlated to the serum levels of any of the considered inflammatory parameters and cytokines (supplementary Table S5).
We assessed markers of cell damage, neutrophil and macrophage activation in hospitalized COVID-19 geriatric patients, identifying three plasma/serum parameters, i.e., integrity of n-cfDNA, NE and sCD163 levels, as significantly associated with in-hospital death .
In particular, concentrations of NE and sCD163 on-hospital admission were higher in patients who deceased in-hospital, in accordance with previous observations that exuberant neutrophil and macrophage activation is often associated with and could contribute to COVID-19 severity [4, 12, 13]. In addition, a third parameter, i.e. a low n-cfDNA integrity, was found to be associated with increased risk of in-hospital mortality. Different values of n-cfDNA integrity are considered to be representative of different origins of cfDNA: a low integrity should indicate a prevalent apoptotic origin, whereas higher integrity can be interpreted as a marker of necrosis [26, 29] or, as we could suggest, of abundant ongoing NETosis. Indeed, in our cohort of geriatric COVID-19 patients, n-cfDNA integrity was positively correlated with plasma NE, a well-recognized biomarker of NETosis. To explain this correlation, we could consider that NETosis is characterized by the release of webs of nuclear (and mithocondrial) DNA  in absence of DNA fragmentation , hence it is conceivable that massive NETosis involves an abundant release of relatively intact cfDNA in plasma. However, with an apparent paradox, n-cfDNA integrity and NE showed reciprocally direct correlation but opposite association with the outcome in the studied cohort (in-hospital death being associated with increased NE and with decreased n-cfDNA integrity). An analysis of NE in patients with different outcome, conducted after stratifying for different levels of cfDNA integrity, showed that NE was associated with a negative outcome only (or more significantly) in patients with high cfDNA integrity, i.e. in the subgroup of patients at lower risk of death. Conversely, most of the patients who deceased in hospital had low cfDNA integrity, but NE levels not significantly different (or only marginally different) from patients who survived. A possible interpretation for such data is that a combined analysis of n-cfDNA integrity and NE could help in identifying two groups of patients both with relatively high risk of mortality, but with different burden of NETosis: 1) patients with massive NETosis, showing high levels of elastase and high cfDNA integrity due to increased release of non-fragmented DNA from neutrophil; 2) patients without massive NETosis, showing levels of NE around or below the median value, and with low cfDNA integrity, possibly indicating increased apoptosis of immune cells. Considering that literature data show massive NETosis and T lymphocyte apoptosis as frequent characteristics of COVID-19 most severe outcomes,[23, 32] the prevalence of one or the other pathogenetic mechanism in different geriatric patients could help explaining the results obtained in this study. However, some considerations suggest caution when interpreting the present results: 1) as a circulating biomarker of NETosis, NE is not completely specific, since it can be released in plasma by neutrophils also by degranulation ; 2) the method of analysis of n-cfDNA employed in this study is not aimed to ascertain if the cfDNA derives from neutrophils or from other cell types. Hence, our hypothesis of different pathogenic mechanisms contributing to circulating n-cfDNA in different COVID-19 patients should be checked in future studies, possibly including the analysis of very specific biomarkers of NETosis such as the plasma levels of citrullinated histone H3  or circulating myeloperoxidase (MPO)-DNA complexes .
Univariate survival analysis confirmed the association of n-cfDNA integrity and sCD163 with survival while, of note, multivariate analysis showed that n-cfDNA integrity, NE and sCD163 independently predict the outcome even if major prognostic clinical variables (SARS-CoV-2 RNAemia, comorbidities, frailty score, age and sex) are taken in consideration as covariates. Significantly, at least two of the studied biomarkers (n-cfDNA integrity, NE) appear to maintain their association with survival also when different interventions and medications are included in the multivariate model as additional confounding variables. Overall, the tested survival models confirmed that a lower risk of death was associated with higher n-cfDNA integrity, lower levels of NE and lower levels of sCD163.
Furthermore, absolute values of both n-cfDNA (Alu 115) and mt-cfDNA were associated with SARS-CoV-2 RNAemia, and n-cfDNA was also negatively correlated with serum INF-alpha. Considering that on-admission SARS-CoV-2 RNAemia and low serum IFN-alpha are both established markers of COVID-19 severe cases, [36, 37] the above reported observations are in agreement with the association between cfDNA and COVID-19 severity shown in other studies [34, 38].
The significantly increased levels of the most relevant cytokines modulating the inflammatory processes, such as IL-6 and IL-10, further support the notion that increased systemic inflammatory conditions are associated with poorer COVID-19 outcomes, especially in the setting of geriatric patients . The observed positive correlation between IL-6, n-cfDNA abundance (ALU 247) and NE in plasma is in agreement with the results of a recent study showing that circulating markers of NET formation are associated with biomarkers of inflammation . Of, note, the same authors also reported a prognostic value of circulating markers of NET formation including plasma level of NE, cfDNA abundance and citrullinated histone H3,  which is confirmed, in our study, for NE. The different results obtained for the parameters of cfDNA abundance (represented in our study by Alu 115 and Alu 247 values, not significantly associated with prognosis) could be attributed to the use of different methods for cfDNA analysis, or to the different (age) selection of COVID-19 cases. In addition, the elevated plasma levels of NET markers found in COVID-19 patients when compared to uninfected controls in the above cited study , could be considered in overall agreement with our results of an increased abundance of n-cfDNA (Alu 115 and Alu 247 DNA) in infected versus uninfected geriatric patients. Concerning the other main biomarkers analysed in this study, the median plasma concentration of NE and sCD163 observed in our cohort of COVID-19 patients was respectively at least five times higher, and one-point-six times higher, than the values (16-21 ng/ml and 400 ng/ml for NE and sCD163, respectively) reported in recent publications for control groups of 60-70 years old subjects [34, 40, 41]. Our study has some limitations that need to be addressed, most notably its retrospective nature and single-center design. In addition, the method used to extract nucleic acids from plasma samples was not specifically designed to collect highly fragmented DNA, and this could have potentially lead to an overall overestimation of the cfDNA integrity index, albeit without bias between different samples.
On the whole, our data show that a combined analysis of n-cfDNA integrity, NE and sCD163 on hospital admission appear to be capable to discriminate geriatric COVID-19 patients at high risk of in-hospital death. If future studies will confirm the present results, the development of standardized clinical laboratory tests for cfDNA integrity, NE and sCD163 could be advisable, in order to facilitate the possible implementation of these biomarkers in the clinical practice.
Materials and methods
The present study utilizes data and biological samples from the Report-Age COVID project, an observational study conducted at the Italian National Center on Aging (IRCCS INRCA), Italy. The aim of this study is to provide a deeper understanding of COVID-19 disease in older hospitalized patients (age>65 years). All the selected subjects were cases of COVID-19 as confirmed by the positive detection of SARS-CoV-2 RNA in nasal/oro-pharyngeal swabs using real-time reverse transcriptase-polymerase chain reaction assay. The study protocol has been approved by the Ethics Committee of the IRCCS INRCA hospital, Ancona, Italy (reference number CE-INRCA-20008) and registered under the ClinicalTrials.gov database (reference number NCT04348396). Clinical and epidemiological data of hospitalization were gathered in a retrospective manner and anonymized prior to release. All the patients enrolled in the Report-Age COVID received treatment from INRCA hospital from March 1st 2020 to date. Among these patients, we selected 156 subjects from the database based on the availability of biological samples (serum and plasma at hospital admission). These selected patients were admitted to INRCA hospital between October 11th and December 31st 2020. None of them received the anti-COVID-19 vaccine. Most of the patients received corticosteroids during their hospital stays, without significant differences with respect to the outcome. Standard oxygen or nonivasive ventilation (NIV) were used for most of the patients, while few patients were transferred to an intensive care unit (ICU). The initial decision to admit or not a COVID 19 older patient to the ICU was taken after a comprehensive evaluation by a multidisciplinary team, composed by the anesthesiologist, the geriatrician, the cardiologist and the palliativist, as appropriate, which also involved the patient and the family, whenever possible.
In addition to the cohort of COVID-19 patients, and in order to compare values for n-cfDNA parameters in hospitalized older patients with and without SARS-CoV-2 infection, a group of 36 older adults was selected from the Report-Age project, a large-scale ongoing observational study on the health conditions of hospitalized older adults at INRCA Hospital (Trial Registration no. NCT01397682) . The study protocol of the Report-Age project has been approved by the Ethics Committee of the IRCCS INRCA hospital, Ancona, Italy . The control group was composed of older adults who accessed the INRCA hospital for common geriatric disorders in the period 2013-2017 (before the COVID-19 pandemic outbreak), with available biological (plasma) samples and clinical and follow-up data. Patients with evidence of infectious or acute respiratory diseases or with intra-hospital mortality, or with less than two-years of follow-up survival after hospital discharge, were excluded. The patients of this control non-COVID-19 group had comparable sex ratio, median age and prevalence of most of the common age-related diseases as the analysed COVID-19 cohort.
Plasma/serum collection and Nucleic acid extraction
EDTA plasma tubes were gently inverted 8 times and centrifuged. Tubes were centrifuged at 2500 x g at 4°C for 15 minutes. After the centrifugation of blood EDTA tube, three layers were obtained (plasma, buffy coat and erythrocytes). About 2.5 ml of plasma were collected from the upper part of the plasma layer, and placed in a new vial (avoiding collecting plasma from the lower part of the plasma layer near the buffy coat). The collected plasma EDTA volume was then subjected to a second centrifugation, conducted for 8 minutes at 10.000 RCF. After the centrifugation, the upper 80% of the volume was carefully collected (avoiding collecting of the lower 20% volume, enriched in cell debris and platelets) and stored in aliquots of 400-500 microliters, to be immediately frozen at -80°C. Immediately after collection, serum tubes were gently inverted 8 times and left at room temperature for 60 minutes to allow clotting, and then centrifuged at 2500 x g at 4°C for 15 minutes. After centrifugation, the top of the supernatant was carefully aspirated and stored in aliquots of 500-1000 microliters, to be immediately frozen at -80°C. Nucleic acids were extracted from 140 μl of plasma using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany). The kit does not discriminate between RNA and DNA, and it was used to extract plasma cell free DNA (of cellular origin) and viral RNA (for detection of SARS-CoV-2 RNAeamia) from the same plasma samples. The purified nucleic acid sample was conserved at -80°C before the analysis.
PCR amplification of mitochondrial cfDNA (mt-cfDNA) was measured using a real-time quantitative assay for Human cytochrome C oxidase subunit III (MT-CO3) gene. All assays were performed on a Rotor-Gene Q detection system (Qiagen) using a 72-well carousel.
The reaction mixture consisted of 2 μl of nucleic acids from plasma samples and 18 μl master mix, which was composed of 7 μl H2O, 10 μl SYBR qPCR Master Mix (Vazyme, Nanjing, China) and 1 μl of 5 μM forward and reverse primers, respectively. The primers used were: forward 5’-ATGACCCACCAATCACATGC-3’, reverse 5’- ATCACATGGCTAGGCCGGAG-3’ (IDT, Coralville, IA).
PCR conditions were set to: 95 °C for 3 min, followed by 40 cycles of 10 s denaturation at 95 °C, 30 s annealing at 55 °C. Absolute quantification of the target sequence in each biological sample was estimated by comparison to a RT-qPCR standard amplification curve. To generate the standard MT-DNA, a selected region of human purified genomic DNA containing the target sequence for MT-CO3 (forward 5’-ATGACCCACCAATCACATGC-3’, reverse 5’-ATCAATAGATGGAGACATAC-3’) was amplified by PCR under the following conditions: initial denaturation at 95° C for 1 min, then 35 cycles of denaturation at 95° C for 30 s, annealing at 50° C for 30 s and extension at 72°C for 1 min. The generated amplicon had a length of 775 nucleotides and a molecular weight of 478879.57 Da. After amplification, the MT-DNA standard was purified by using MinElute columns (Qiagen) and quantified using a NanoDrop spectrophotometer (ThermoScientific). Serial dilutions were then used to calibrate the RT-qPCR standard curves. Each sample was quantified in duplicate, and triplicates of the standard curve were included in each run.
Determination of abundance and integrity of n-cfDNA
The polymerase chain reaction protocol used for the quantification of nuclear cfDNA (n-cfDNA) and for the evaluation of n-cfDNA integrity was derived from a previous study . Two sets of primers complementary to the consensus sequence of human Alu repeats were used: ALU115 primer forward, 5’-CCTGAGGTCAGGAGTTCGAG-3’; ALU115 primer reverse, 5’-CCCGAGTAGCTGGGATTACA-3’; ALU247 primer forward, 5’-GTGGCTCACGCCTGTAATC-3’; ALU247 primer reverse, 5’-CAGGCTGGAGTGCAGTGG-3’. The first set was used to amplify the 115-bp amplicon (ALU 115), the second for the 247-bp amplicon (ALU 247). The reaction was conducted using 0.2 μM each of forward primer and reverse primer (ALU 115 or ALU 247), 1x Ssofast mix containing enzyme and Syber Green dye (Bio-Rad Laboratories), in a reaction volume of 15 μl. Thermal cycles were conducted on a Qiagen Rotor-Gene instrument, as follows: pre-denaturation 2.5 min at 98°C, followed by 35 cycles of denaturation 15 sec at 95°, annealing/extension 60 sec at 64°C. A melting curve was conducted at the end of each run, to check for the presence of unspecific amplification. A calibration curve constructed by amplifying serial dilutions of a genomic DNA standard sample (0.02 to 200 pg/μl) was present in each assay run, and was used to assess the absolute equivalent concentration of genomic DNA in each plasma sample. Each sample was run in duplicate.
All assays were conducted in blind without knowing the sample identity. The ratio of concentrations calculated for each sample with the two set of primers, i.e. (concentration of Alu 247) / (concentration of Alu 115), hereafter indicated as Alu247/115, was used as index of n-cfDNA integrity.
Measurement of soluble biomarkers
Serum concentration of IL-6, IFN-α, TNF-α and IL-10 were assessed by using four Pro Quantum Immunoassays (Thermo Fisher Scientific) specific for each of the four biomolecules. Each reaction was conducted in duplicate, using 2 μl of serum, and run on an Agilent Aria Mx Real Time PCR instrument at the conditions indicated by the producer.
The plasma levels of NE and LL37 were measured using commercial ELISA kits (ABCAM, AB119553 PMN Elastase Human ELISA Kit and Hycult Biotech, HK321 Human LL-37). The serum levels of CD163 were measured using a commercial ELISA kit (R&DSystems, Human CD163 Immunoassay). All the serum/plasma samples were diluted 1:20 (for sCD163 and LL37 assayes) or 1:50 (for neutrophile elastase assay). All samples were ran in duplicate and samples with an intra-assay coefficient of variation below 10.0% were included in this study.
Assessment of RNAemia
The presence of SARS-CoV-2 RNA in plasma was assessed by analyzing 10 μl of nucleic acids extracted from plasma in a Real Quality RQ-SARS-CoV-2 real time PCR assay (AB Analitica). The assay was specific for the genes RdRp (encoding the RNA-dependent RNA polymerase) gene and S gene (encoding the spike protein). The presence of an endogenous control (targeting the human RNAse P gene) in the same assay was used as internal control to ensure proper sampling and nucleic acids extraction. The assay was used following the conditions indicated by the producer, on an Agilent Aria Mx Real Time PCR instrument. The plasma sample was considered positive for the SARS-CoV-2 RNAemia if both viral genes were detected with a cycle threshold (CT) < 40.
Continuous variables were reported as median and interquartile range after their non-normality have been assessed using Shapiro-Wilk test and comparison of variables between groups was performed by unpaired Mann-Whitney U test. Tertiles of NE, sCD163 and n-cfDNA integrity (Alu247/115) were calculated and, as other categorical variables, were expressed as absolute number.
Spearman rank correlation was conducted to check associations among different biomarkers and between these biomarkers and other parameters.
The Kaplan–Meier survival curves were used to estimate the association of in-hospital mortality risk with different levels of the analysed biomarkers. Cox proportional hazards analysis were used to derive age- and gender-adjusted (Model 1); age- gender- comorbidities- and Clinical Frailty Scale (CFS)-adjusted (Model 2); age- gender- comorbidities- Clinical Frailty Scale (CFS)- and SARS-CoV-2 RNAemia-adjusted (Model 3) and age- gender- comorbidities- Clinical Frailty Scale (CFS)- SARS-CoV-2 RNAemia- and treatments-adjusted (Model 4 and Model 5) hazard ratios (HR) and 95% confidence intervals (95% CI) of the association between all independent variables and study outcome. The choice to consider different Cox proportional hazards models, each containing the same tested variables (biomarkers), but adjusted for an increasing number of confounders, was aimed to verify if a statistically significant relationship between a biomarker and survival maintained its significance after the addition of confounders, and which (if any) confounders could determine a loss of significance in the tested variables. The length of hospital stay was used as the time to failure variable for the model. A two-tailed P value < .05 was considered significant.
The statistical analyses were performed using the IBM SPSS Statistics program (version 27) and the STATA version15.1 Statistical Software Package for Windows (Stata Corp, College Station, TX).
Availability of data and materials
The dataset obtained during the current study are part of a more extensive unpublished work, and thus are not publicly available. However, manuscript data can be made available from the corresponding author upon reasonable request.
COVID-19 Dashboard by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU). https://coronavirus.jhu.edu/map.html. Accessed 25/03/2022.
Prendki V, Tiseo G, Falcone M. ESCMID Study Group for Infections in the Elderly (ESGIE). Caring for older adults during the COVID-19 pandemic. Clin Microbiol Infect. 2022;28(6):785–91. https://doi.org/10.1016/j.cmi.2022.02.040.
Giamarellos-Bourboulis EJ, Netea MG, Rovina N, Akinosoglou K, Antoniadou A, Antonakos N, et al. Complex Immune Dysregulation in COVID-19 Patients with Severe Respiratory Failure. Cell Host Microbe. 2020;27:992–1000.e3. https://doi.org/10.1016/j.chom.2020.04.009.
Ombrello MJ, Schulert GS. COVID-19 and cytokine storm syndrome: are there lessons from macrophage activation syndrome? Transl Res. 2021;232:1–12. https://doi.org/10.1016/j.trsl.2021.03.002.
Sabbatinelli J, Matacchione G, Giuliani A, Ramini D, Rippo MR, Procopio AD, et al. Circulating biomarkers of inflammaging as potential predictors of COVID-19 severe outcomes. Mech Ageing Dev. 2022;204:111667. https://doi.org/10.1016/j.mad.2022.111667.
Zingaropoli MA, Nijhawan P, Carraro A, Pasculli P, Zuccalà P, Perri V, et al. Increased sCD163 and sCD14 Plasmatic Levels and Depletion of Peripheral Blood Pro-Inflammatory Monocytes, Myeloid and Plasmacytoid Dendritic Cells in Patients With Severe COVID-19 Pneumonia. Front Immunol. 2021;12:627548. https://doi.org/10.3389/fimmu.2021.627548.
Sánchez-Torres C, Gómez-Puertas P, Gómez-del-Moral M, Alonso F, Escribano JM, Ezquerra A, et al. Expression of porcine CD163 on monocytes/macrophages correlates with permissiveness to African swine fever infection. Arch Virol. 2003;148:2307–23.
Segal AW. How neutrophils kill microbes. Annu Rev Immunol. 2005;23:197–223.
Brinkmann V, Zychlinsky A. Beneficial suicide: why neutrophils die to make NETs. Nat Rev Microbiol. 2007;5(8):577–82.
Kruger P, Saffarzadeh M, Weber AN, Rieber N, Radsak M, von Bernuth H, et al. Neutrophils: Between Host Defence, Immune Modulation, and Tissue Injury. PLoS Pathog. 2015;11(3):e1004651. https://doi.org/10.1371/journal.ppat.1004651.
Mehta P, McAuley DF, Brown M, Sanchez E, Tattersall RS, Manson JJ, et al. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet. 2020;395(10229):1033–4.
Szturmowicz M, Demkow U. Neutrophil Extracellular Traps (NETs) in Severe SARS-CoV-2 Lung Disease. Int J Mol Sci. 2021;22(16):8854. https://doi.org/10.3390/ijms22168854.
Iliadi V, Konstantinidou I, Aftzoglou K, Iliadis S, Konstantinidis TG, Tsigalou C. The Emerging Role of Neutrophils in the Pathogenesis of Thrombosis in COVID-19. Int J Mol Sci. 2021;22(10):5368. https://doi.org/10.3390/ijms22105368.
Block H, Zarbock A. A Fragile Balance: Does Neutrophil Extracellular Trap Formation Drive Pulmonary Disease Progression? Cells. 2021;10(8):1932. https://doi.org/10.3390/cells10081932.
Narasaraju T, Yang E, Samy RP, Ng HH, Poh WP, Liew AA, et al. Excessive neutrophils and neutrophil extracellular traps contribute to acute lung injury of influenza pneumonitis. Am J Pathol. 2011;179:199–210.
Ashar HK, Mueller NC, Rudd JM, Snider TA, Achanta M, Prasanthi M, et al. The Role of Extracellular Histones in Influenza Virus Pathogenesis. Am J Pathol. 2018;188(1):135–48. https://doi.org/10.1016/j.ajpath.2017.09.014.
Cahilog Z, Zhao H, Wu L, Alam A, Eguchi S, Weng H, et al. The Role of Neutrophil NETosis in Organ Injury: Novel Inflammatory Cell Death Mechanisms. Inflammation. 2020;43(6):2021–32. https://doi.org/10.1007/s10753-020-01294-x.
Grayson PC, Kaplan MJ. At the Bench: Neutrophil extracellular traps (NETs) highlight novel aspects of innate immune system involvement in autoimmune diseases. J Leukoc Biol. 2016;99(2):253–64. https://doi.org/10.1189/jlb.5BT0615-247R.
Moorthy AN, Tan KB, Wang S, Narasaraju T, Chow VT. Effect of high-fat diet on the formation of pulmonary neutrophil extracellular traps during influenza pneumonia in BALB/c mice. Front Immunol. 2016;7:289.
Rudd JM, Pulavendran S, Ashar HK, Ritchey JW, Snider TA, Malayer JR, et al. Neutrophils Induce a Novel Chemokine Receptors Repertoire During Influenza Pneumonia. Front Cell Infect Microbiol. 2019;9:108. https://doi.org/10.3389/fcimb.2019.00108.
Kustanovich A, Schwartz R, Peretz T, Grinshpun A. Life and death of circulating cell-free DNA. Cancer Biol Ther. 2019;20(8):1057–67. https://doi.org/10.1080/15384047.2019.1598759.
Khier S, Lohan L. Kinetics of circulating cell-free DNA for biomedical applications: critical appraisal of the literature. Future Sci OA. 2018;4(4):FSO295. https://doi.org/10.4155/fsoa-2017-0140.
Barnes BJ, Adrover JM, Baxter-Stoltzfus A, Borczuk A, Cools-Lartigue J, Crawford JM, et al. Targeting potential drivers of COVID-19: Neutrophil extracellular traps. J Exp Med. 2020;217(6):e20200652. https://doi.org/10.1084/jem.20200652.
Thierry AR, Roch B. Neutrophil Extracellular Traps and By-Products Play a Key Role in COVID-19: Pathogenesis, Risk Factors, and Therapy. J Clin Med. 2020;9(9):2942. https://doi.org/10.3390/jcm9092942.
Zuo Y, Yalavarthi S, Shi H, Gockman K, Zuo M, Madison JA, et al. Neutrophil extracellular traps in COVID-19. JCI Insight. 2020;5(11):e138999. https://doi.org/10.1172/jci.insight.138999.
Umetani N, Kim J, Hiramatsu S, Reber HA, Hines OJ, Bilchik AJ, et al. Increased integrity of free circulating DNA in sera of patients with colorectal or periampullary cancer: direct quantitative PCR for ALU repeats. Clin Chem. 2006;52(6):1062–9. https://doi.org/10.1373/clinchem.2006.068577.
Cheng AP, Cheng MP, Gu W, Sesing Lenz J, Hsu E, Schurr E, et al. Cell-free DNA tissues of origin by methylation profiling reveals significant cell, tissue, and organ-specific injury related to COVID-19 severity. Med (N Y). 2021;2(4):411–422.e5. https://doi.org/10.1016/j.medj.2021.01.001.
Storci G, Bonifazi F, Garagnani P, Olivieri F, Bonafè M. The role of extracellular DNA in COVID-19: clues from inflamm-aging. Ageing Res Rev. 2021;66:101234. https://doi.org/10.1016/j.arr.2020.101234.
Jiang P, Chan CW, Chan KC, Cheng SH, Wong J, Wong VW, et al. Lengthening and shortening of plasma DNA in hepatocellular carcinoma patients. Proc Natl Acad Sci U S A. 2015;112(11):E1317–25. https://doi.org/10.1073/pnas.1500076112.
Han DSC, Lo YMD. The Nexus of cfDNA and Nuclease Biology. Trends Genet. 2021;37(8):758–70. https://doi.org/10.1016/j.tig.2021.04.005.
Fuchs TA, Abed U, Goosmann C, Hurwitz R, Schulze I, Wahn V, et al. Novel cell death program leads to neutrophil extracellular traps. J Cell Biol. 2007;176(2):231–41. https://doi.org/10.1083/jcb.200606027.
André S, Picard M, Cezar R, Roux-Dalvai F, Alleaume-Butaux A, Soundaramourty C, et al. T cell apoptosis characterizes severe Covid-19 disease. Cell Death Differ. 2022;29(8):1486–99. https://doi.org/10.1038/s41418-022-00936-x.
Parzy E, Bouchaud V, Massot P, Voisin P, Koonjoo N, Moncelet D, et al. Overhauser-enhanced MRI of elastase activity from in vitro human neutrophil degranulation. PLoS One. 2013;8(2):e57946. https://doi.org/10.1371/journal.pone.0057946.
Ng H, Havervall S, Rosell A, Aguilera K, Parv K, von Meijenfeldt FA, et al. Circulating Markers of Neutrophil Extracellular Traps Are of Prognostic Value in Patients With COVID-19. Arterioscler Thromb Vasc Biol. 2021;41(2):988–94. https://doi.org/10.1161/ATVBAHA.120.315267.
Donkel SJ, Wolters FJ, Ikram MA, de Maat MPM. Circulating Myeloperoxidase (MPO)-DNA complexes as marker for Neutrophil Extracellular Traps (NETs) levels and the association with cardiovascular risk factors in the general population. PLoS One. 2021;16(8):e0253698. https://doi.org/10.1371/journal.pone.0253698.
Miki S, Sasaki H, Horiuchi H, Miyata N, Yoshimura Y, Miyazaki K, et al. On-admission SARS-CoV-2 RNAemia as a single potent predictive marker of critical condition development and mortality in COVID-19. PLoS One. 2021;16(7):e0254640. https://doi.org/10.1371/journal.pone.0254640.
Contoli M, Papi A, Tomassetti L, Rizzo P, Vieceli Dalla Sega F, Fortini F, et al. Blood Interferon-α Levels and Severity, Outcomes, and Inflammatory Profiles in Hospitalized COVID-19 Patients. Front Immunol. 2021;12:648004. https://doi.org/10.3389/fimmu.2021.648004.
Hammad R, Eldosoky MAER, Fouad SH, Elgendy A, Tawfeik AM, Alboraie M, et al. Circulating cell-free DNA, peripheral lymphocyte subsets alterations and neutrophil lymphocyte ratio in assessment of COVID-19 severity. Innate Immun. 2021;27(3):240–50. https://doi.org/10.1177/1753425921995577.
Bonafè M, Prattichizzo F, Giuliani A, Storci G, Sabbatinelli J, Olivieri F. Inflamm-aging: Why older men are the most susceptible to SARS-CoV-2 complicated outcomes. Cytokine Growth Factor Rev. 2020;53:33–7. https://doi.org/10.1016/j.cytogfr.2020.04.005.
Smith P, Rosell A, Farm M, Bruzelius M, Aguilera Gatica K, Mackman N, et al. Markers of neutrophil activation and neutrophil extracellular traps in diagnosing patients with acute venous thromboembolism: A feasibility study based on two VTE cohorts. PLoS One. 2022;17(7):e0270865. https://doi.org/10.1371/journal.pone.0270865.
Daftarian N, Zandi S, Piryaie G, Nikougoftar Zarif M, Ranaei Pirmardan E, Yamaguchi M, et al. Peripheral blood CD163(+) monocytes and soluble CD163 in dry and neovascular age-related macular degeneration. FASEB J. 2020;34(6):8001–11. https://doi.org/10.1096/fj.201901902RR.
Bustacchini S, Abbatecola AM, Bonfigli AR, Chiatti C, Corsonello A, Di Stefano G, et al. The Report-AGE project: a permanent epidemiological observatory to identify clinical and biological markers of health outcomes in elderly hospitalized patients in Italy. Aging Clin Exp Res. 2015;27:893–901. https://doi.org/10.1007/s40520-015-0350-3.
The authors wish to thank Belinda Giorgetti, Beatrice Bartozzi, Denise Quinzi, and Elisabetta Farnocchia for their invaluable work in the detection of SARS-CoV-2 infection by performing real-time reverse transcriptase-polymerase chain reaction assays in nasal/oro-pharyngeal swab samples and in plasma samples. The authors wish also to thank Dr. Roberta Galeazzi and the other biologists and laboratory technicians of the Clinical Laboratory of the INRCA hospital for the collection of serum and plasma samples and for the routine clinical laboratory data. We also thank BioGer, the IRCCS INRCA institutional biobank, for providing the biological samples.
This work was supported by Ricerca Corrente to IRCCS INRCA (Italian Ministry of Health).
Ethics approval and consent to participate
The study was conducted according to the guidelines of the Declaration of Helsinki. Patients selected for the study were enrolled according to The Report-Age COVID and Report-Age protocol studies that have been approved by the Ethics Committee of the IRCCS INRCA, Ancona, Italy (reference numbers CE-INRCA-20008 and SC/10/272) and registered under the ClinicalTrials.gov database (reference numbers NCT04348396 and NCT01397682). All the statistical analysis was done on anonymized data. All research was conducted in accordance with relevant guidelines and regulations.
Consent for publication
No details/ images/ videos of an individual person are included in this manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Supplementary Table S1. Characteristics and comorbidities of non-COVID-19 patients.
Additional file 2:
Supplementary Table S2. Parameters of n-cfDNA in non-COVID-19 geriatric patients and comparison with COVID-19 patients.
Additional file 3:
Supplementary Table S3. Cox proportional hazard ratios (HR) for in-hospital survival of COVID-19 patients based on different multivariate models (with pharmacological treatments/oxygen therapy as additional confounding variables).
Additional file 4:
Supplementary Table S4. Correlation among markers of cell damage, neutrophil and macrophage activation.
Additional file 5:
Supplementary Table S5. Correlation of markers of cell damage, neutrophil and macrophage activation with other laboratory parameters.
Additional file 6:
Supplementary Figures. Scatterplots with the individual data points of the correlation between NE levels and n-cfDNA parameters.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Cardelli, M., Pierpaoli, E., Marchegiani, F. et al. Biomarkers of cell damage, neutrophil and macrophage activation associated with in-hospital mortality in geriatric COVID-19 patients. Immun Ageing 19, 65 (2022). https://doi.org/10.1186/s12979-022-00315-7
- Neutrophil elastase
- Prognostic biomarker
- Geriatric patients
- in-hospital mortality