1. Introduction
Clinical features of SARS-CoV-2 infection are various, from mild to moderate symptoms, characterized by spontaneous regression, to severe cases, showing conditions leading to ICU admission or death[
1,
2]. Despite virus transmission and clinical presentation are well known, scanty information is available about the pathophysiology of disease worsening.
The severe respiratory failure caused by COVID-19 is primarily manifested as acute respiratory distress syndrome (ARDS). Its related pathological findings were often documented in most COVID-19 autoptic studies[
3]. ARDS nosography has been substantially modified after COVID-19 pandemic, since several studies have shown that COVID-19-related ARDS is atypical, although often falling within the Berlin definition of ARDS[
4,
5].
As a matter of fact, inflammation plays a key role in the disease worsening, if any, of COVID-19 patients. Several studies have focused on the cytokine storm and the lymphocyte fall observed in these patients[
6,
7,
8]. The relationship between the clinical presentation of severe cases of COVID-19 and the related immuno-inflammatory factors have not so far well understood[
9].
The early identification of cases with greater clinical severity is mandatory to ensure the best possible care. The detection of laboratory bio-markers or indices representing the unbalanced innate and adaptive immune responses allows a rapid characterization of patients in a hospital-setting. Venous and arterial samplings, together with a chest computer tomography (CT) could quickly provide enough information to recognize hyper-inflammatory state and respiratory distress, to assess clinical severity and predict outcome.
Several biochemical parameters were under investigation for an early assessment of disease severity[
10,
11]. Among them, C-reactive protein (CRP) and Neutrophil-to-lymphocyte ratio (NLR) are rapid and widely available inflammatory indices, recently under debate as predictors of bad outcome in a large variety of pathological conditions characterized by exaggerated systemic inflammation[
12,
13]. In our previous study we found that both NLR and CRP are reliable prognostic predictors in COVID-19 patients. [
14].
Neutrophilia, lymphopenia and high CRP levels are closely linked with the pathophysiology of COVID-19. Cytokine storm determines a massive recruitment of circulating cells resulting in a continuous activation of neutrophils, incremental consumption and negative counter-regulation of lymphocytes and a strengthened inflammatory response, with an over-production of CRP, that is a protein of acute phase.
In the present study, we aimed to evaluate the time-course of immuno-inflammatory pathways involvement in the development of respiratory failure. We especially evaluated the dynamic changes of neutrophils, lymphocytes and CRP during hospitalization. Finally, we also assessed the link between CRP, neutrophils and lymphocytes with P/F ratio, to better understand the individual role of mediators leading to respiratory failure.
2. Materials and Methods
In this retrospective, multicentre, observational study, a total of 963 patients were consecutively candidates among patients admitted to the Emergency Unit of San Marco Hospital and the COVID-19 Internal Medicine Unit at Cannizzaro Hospital, in the city of Catania, Italy, between October 2020 and September 2022. Within this cohort, 764 patients were selected, based on predefined inclusion criteria (
Figure 1):
- 1)
a reliable diagnosis of SARS-CoV2 infection obtained by RT-PCR molecular swab testing;
- 2)
no history of pharmacological treatments responsible for alterations in the leukocyte count and/or CRP upon admission;
- 3)
no current or past history of conditions responsible for alterations in the leukocyte count and/or CRP;
- 4)
availability of at least three blood tests and blood gas analyses during hospitalization, and a hospitalization period not less than 48 hours.
Routine biomarkers were measured using standard techniques using auto-analyzers (Beckman Coulter DxH 800; Danaher Corporation, Miami, FL, USA, Beckman DxC 700 AU; Danaher Corporation, Miami, FL, USA).
Total and differential leukocyte count and CRP levels were measured three times: at baseline, on the median day of hospitalization, and at discharge (defined as death, transfer to the Intensive Care Unit, ICU, or routine discharge).
The severity of respiratory failure was identified using P/F ratio obtained by blood gas test reports, using PaO2/FiO2 ratio, namely the ratio between. Partial arterial pressure of oxygen/Fraction of inspired oxygen.
The neutrophil-to-lymphocyte ratio (NLR) was calculated as follows: NLR = number of neutrophils/number of lymphocytes.
Relevant information was obtained by reviewing the medical records of the participants, including their demographic and anamnestic data, clinical and laboratory characteristics, treatment regimen, as well as outcome.
Appropriate measures were taken to display information before performing statistical analysis.
The present study exclusively focused on death following SARS-CoV2 infection, so excluding those caused by any other factors.
2.1. Statistical Analysis
The statistical analysis was performed using the IBM-SPSS and R statistical software packages.
Categorical variables were described in terms of absolute frequency and percentage prevalence, while continuous variables were further divided into two groups, after performing the Kolmogorov-Smirnov test to evaluate their distribution. Continuous variables were expressed as mean ± standard deviation in presence of normal distribution and median and interquartile range in presence of non- Gaussian distribution.
The sample was initially divided into three groups, according to outcome. Differences between outcome were evaluated using the chi-square test with Fisher’s correction for categorical variables, one-way ANOVA for normally distributed continuous variables, and the Kruskal-Wallis test for non-normally distributed continuous variables.
Non-parametric tests, such as Spearman’s correlation coefficient and Kendall’s Tau, were used to assess the correlations between the variables of interest, considering their distribution.
Where appropriate, we used the z-scores for variables characterized by extreme kurtosis and asymmetry.
Firstly, to evaluate the relationship between NLR and CRP as continuous variables with Pa02/FiO2 ratio, univariate and multivariate logistic regression models were used, and odds ratios (ORs), both corrected and uncorrected, as well as 95% confidence intervals (CIs), were calculated.
All these relationships were assessed in the whole sample and also according to single outcome.
Multiple logistic and Cox regression models were also built to verify the associations between biomarkers and outcome.
Repeated measures ANOVA was performed, whenever necessary, to test the differences over time in the mean levels of variables.
Finally, to better understand the existing relationships between inflammatory biomarkers and P/F ratio, we performed a mediation analysis, to verify if these relationships resulted mediated by an external variable. A moderation analysis was also done to estimate if the effect of such variable could influence the strength and direction of these relationships, after verifying that the required assumptions were met.
3. Results
Figure 1 depicts the flow chart of patient recruitment, as well as different outcomes.
Median age of the whole sample was 74 years and 54% were men. Median period of hospitalization was 9 days.
Compared to survivors, patients deceased or admitted to ICU showed higher median levels of NLR, CRP, leukocyte and neutrophil count, lower values of lymphocyte count and P/F ratio. Regarding age, significant differences were only observed between deceased and survivors.
As regards comorbidities, no statistically differences were observed between the various outcomes. All demographic and clinical characteristics of the patients are shown in
Table 1.
First, we analysed the temporal trend shown by NLR and CRP during hospitalization, and assessed the differences by repeated measure ANOVA, in the whole sample and for each subgroup (
Figure 2A,B).
While whole sample, survivors and ICU-admitted patients showed a flat pattern of NLR over time, deceased patients, by contrast, showed an increase of NLR from baseline over time. At variance, CRP tended to fall from baseline to median day of hospitalization in all four subgroups, but steeply increased at the end of hospitalization only in ICU-admitted patients.
As shown in
Figure 2, admission to ICU occurred within an interval corresponding to the peak of CRP, while death occurred within the interval in which the maximum peak of NLR was observed.
The correlation between NLR and P/F was not statistically significant in ICU-admitted patients (ρ = -0.03 p = 0.7203,
Figure 3a), but resulted significant in deceased patients (ρ = -0.26 p = 0.0025;
Figure 3b). An opposite phenomenon was observed in the correlation between CRP and P/F ratio in the two subgroups identifying the outcome (ρ = -0.24, p = 0.0129; ρ = -0.13 p = 0.865) (
Figure 3c,d).
Furthermore, we analysed these relationships using multiple linear regression analysis, being P/F ratio the dependent variable, and then using univariate and multivariate logistic regression analyses with single outcomes as dependent variable (
Table 2 and
Table 3). All multivariate regression models were adjusted for age, sex, and comorbidities.
In the multiple linear regression model (
Table 2) adjusted for age, sex, and comorbidities, NLR and CRP significantly predicted P/F in both whole population and in survivors, while only CRP predicted P/F in ICU-admitted patients; conversely, only NLR predicted P/F in deceased patients.
In Cox proportional hazard regression, NLR predicted mortality independently of CRP, which did not result statistically significant, and other confounders. ICU admission was significantly associated with both biomarkers, although CRP showed a higher HR, a more significant p value, and a more restricted confidence interval (
Table 3).
Based on our results, showing that NLR and CRP are differently associated with P/F in deceased and ICU-admitted patients, we wondered whether the relationship between CRP and P/F could be mediated by another variable that could partly explains why CRP has a greater weight than NLR in influencing P/F in ICU patients, with no significant association with P/F in deceased patients.
Observing the trend of NLR over time, we used as mediating variables the two cell populations from which it is calculated, namely absolute neutrophil and lymphocyte counts. We first demonstrated by repeated measures ANOVA that the mean levels of neutrophils and lymphocytes differed significantly in the various time intervals, comparing patients grouped for outcome (
Figure 4).
Finally, mediation analysis was made to examine the mediating effect of neutrophils, lymphocytes, and age on the relationship between CRP and P/F. Age, ANC, CRP and Lymphocytes significantly and directly influenced P/F, while the influence of CRP on P/F was also mediated by ANC, with no mediating effect of lymphocytes and age (
Figure 5).
In
Table 4 are shown details of the significant indirect effect of CRP on P/F, mediated by neutrophils (b= -0.035, p=0.001).
Furthermore, the direct effect of CRP on P/F in presence of the neutrophils as mediator was stronger (83.7 %; b = -2.849 p < 0.001), as compared to indirect effect mediated by neutrophils (16.3%; b-0.035; p=0.001). Hence, these results suggest that neutrophils partially mediated the relationship between CRP and P/F.
4. Discussion
To the best of our knowledge this is the first study evaluating the pathways of main immuno-inflammatory circulating biomarkers, as well as their direction, in COVID-19 patients. Of note, we also demonstrated the relationships between two markers we have already described as prognostic factors in COVID-19 patients[
14], namely NLR and CRP, with P/F, that is in turn a marker of respiratory failure secondary to respiratory distress .
Type 1 respiratory failure in COVID-19 patients is characterized by high compliance and a lack of apparent respiratory distress, at least in the early phase, also known as silent hypoxemia[
5]. This could lead to an underestimation of the disease’s severity, and could seriously impact on decision-making, so delaying ICU admission. Consequently, emergency departments are facing a sneaky clinical scenario, sometimes leading to underestimate disease severity.
Recently Sinatti et al.[
15] highlighted the potential of PaO2/FiO2 in predicting pneumonia progression towards ARDS in COVID-19 patients.
Of note, the acute respiratory failure in COVID-19 is characterized by a specific pathological substrate: the
primum movens is endothelial damage, as a consequence of the hyperinflammatory state, with a massive cytokine and immune cell storm. In this direction, ARDS is the first result of V/Q mismatch due to a vascular damage[
16]. In this respect, recent studies have highlighted the link between inflammation, arterial stiffening, and cardiovascular events[
17]. In this direction, our group [
18] recently demonstrated the stiffening of elastic arteries following COVID-19 infection, followed by partial regression in survivors, in contrast to the higher brachial-ankle pulse wave velocity (PWV) shown in deceased.
Overall, mechanisms underlying the close connection between respiratory failure and systemic inflammation are still under debate. We aimed to investigate this link in order to identify the time-course and dynamics of release of inflammatory substances and activation of immune system during hospitalization, that should be considered an alarm bell.
Laboratory abnormalities have been linked to adverse outcomes in COVID-19 patients. In our previous study[
14], we assessed the predictive value of inflammatory biomarkers, such as NLR and CRP, on mortality and severe COVID-19 disease by examining their trends during hospitalization. In our previous study we demonstrated the different CRP and NLR behaviour, as related to mortality and ICU admission.
In the present study we attempted to understand whether these two biomarkers are differently involved in determining outcome and to identify the relationships, if any, between both CRP and NLR with P/F ratio. We demonstrated for the first time the influence of CRP and NLR on P/F, so giving ground to the concept that a derangement in host immune-inflammatory response can influence disease severity.
CRP is a non-specific acute phase protein, induced by IL-6 in the liver, and a sensitive biomarker of inflammation, infection and tissue damage. The increase of CRP levels generally starts as early as 4-8 hours after the initiation of the inflammatory process and peaks at 48 hours, although the duration of this peak is variable, and directly proportional to the stimulus persistence, with a half-life of about 19 hours[
19]. The higher CRP levels observed in COVID-19 patients requiring intensive treatment is the consequence of systemic hyper-inflammatory state occurring in severe COVID-19 cases[
20]. Despite its low specificity, CRP is an helpful marker in a lot of acute conditions characterized by inflammation/infection[
12], and its elevation was shown to be associated with a bad outcome in cardiovascular diseases[
21,
22]. It is crucial to emphasize, in this respect, the functional role of CRP in the inflammatory process, aimed at recruiting complement components, with positive feedback on inflammation, especially involving endothelium[
23,
24,
25]. In COVID-19 patients CRP has been already described as bad outcome predictor, together with the cytokines associated with its expression (IL-6, IL-10)[
26,
27].
In the present study, the prominent peak in NLR (
Figure 2A) should be finalized to face the disease. However, the early activation of the innate immune response (increase in neutrophils) and the subsequent fall in lymphocytes (adaptive immune response) characterizes the host response to COVID-19 infection[
28].
Lowery et al. [
29] described the kinetics of innate immune system, demonstrating that COVID-19 pathophysiology depends on the so-called cytokine storm, with production of cytokines and chemokines (TNF, IL-6, CXCL10, CCL2, CCL5 and IFN-II), and/or the lack of early IFN-I and IFN-III expression. The consequent prolonged activation of the innate immune system and a continuous suppression of lymphocytes release into circulation are main characteristics of COVID-19 disease. An elevated absolute neutrophil count (ANC) has been described as negative predictor of outcome in COVID-19 patients[
30,
31], so reflecting the hyper-activation of the innate immune response caused by both virus-triggered or cytokine-dependent mechanisms[
13]. In particular, neutrophils are involved in platelet activation, over-production of inflammatory cytokines, epithelial and endothelial cells damage, especially through a process named NETsosis, that is a molecular mechanism leading to formation of neutrophil extracellular traps (NETs)[
32]. As a matter of fact, an elevated expression of NETs could be associated with death and disease progression in COVID-19 patients[
33].
By contrast, adaptative immune response induces a reduction of lymphocyte absolute count in COVID-19, as a consequence of extended TNF-α-induced apoptosis, peripheral consumption, direct ACE-2-cytopathic effect, or through the interaction with CD147[
29,
34,
35,
36]. Peripheral lymphopenia in COVID-19 patients has been largely described in recent literature assessing its predictive value for disease severity and mortality[
37,
38]. In addition, neutrophilia itself leads to a suppression of lymphocytes, through a cytotoxic indirect effect[
39,
40].
NLR represents the balance between innate and immune response. It readily increases as a consequence of physiological and pathophysiological response to acute stress[
29]. NLR may be considered a marker of subclinical inflammation, with higher values in acute exacerbations, such as CAP, COPD, sepsis, cancer and many cardiac diseases[
41,
42,
43,
44,
45]. It is strictly linked to the immune system derangement and could be used as a predictor of disease severity and mortality, especially in conditions characterized by systemic inflammation, such as COVID-19[
14,
46,
47].
Considering that respiratory failure in COVID-19 is characterized by a strong inflammatory involvement, we sought to identify a potential link between the host response to viral load and severity degree of COVID-19 disease. Some studies have already identified an association between CRP and respiratory failure. Poggiali et al.[
48] and Herold et al.[
49], in two single-center studies, demonstrated that higher CRP levels are associated with lung function worsening, but data on their time-course throughout hospitalization and their correlation with outcome were lacking. On the other hand, Mueller et al.[
50], in a single-center study, showed the temporal relationship between CRP and P/F, without any data referred to outcome.
In a recent study we emphasized that NLR independently predicted mortality and ICU admission in COVID-19 patients, being this prediction abolished after adjustment for P/F[
14].
Here we further tested the pathways linking neutrophils, lymphocytes, age and CRP to P/F. In consideration of the differences in the time-course of NLR and CRP between ICU admitted and deceased during hospitalization (
Figure 2A,B), we investigated the relationships of these biomarkers with P/F, according to outcome, and found a significant inverse correlation between NLR and P/F in deceased (
Figure 3B) and a significant inverse correlation between CRP and P/F in the subgroup of ICU-admitted (
Figure 3C).
Furthermore, we corroborated the robustness of these correlations by a multiple linear regression analysis, demonstrating that NLR and CRP significantly predicted P/F in the whole population, while only CRP predicted P/F in ICU-admitted patients, whereas only NLR predicted P/F in deceased patients (
Table 2). These findings would imply that inflammation plays a key role in the worsening of respiratory failure, but also that a dysfunction of immune system, as shown by neutrophilia/lymphopenia, would be closely associated with decease.
Clinical implications of our findings look interesting. For each increase of 1 unit of standard deviation of NLR, the risk of mortality increased by 77% in the whole sample. As to ICU admission, NLR and CRP have both a specific weight, especially CRP, whose increment of 1 unit of standard deviation corresponded to a 70% increase in the risk of ICU admission.
Finally, we investigated the pathways and their direction influencing P/F. In this way, we demonstrated that mean levels of neutrophils and lymphocytes, whose absolute counts determine NLR, have a significantly different trend during hospitalization in ICU-admitted and deceased (
Figure 4). Besutti et al. [
51] have recently reported that persistent lung abnormalities had a link to inflammatory burden, especially to CRP, reflecting the intensity and the duration of inflammatory reaction. In keeping with Besutti’s data [
51], we considered the potential influence of neutrophils and lymphocytes in the relationship between CRP, a marker of inflammation, and P/F, a marker of lung function.
We observed that, in the worsening of lung function, CRP had a large direct effect (83.7%) (Tab. 4) on P/F: the higher was CRP at hospital admission, the lower resulted P/F (
Figure 5) (path coefficient: - 0.19; p<0.05). Neutrophils significantly mediated the deleterious effect of inflammation on lung function. In fact, neutrophils count indirectly influenced (16.3%) the relationship between CRP and P/F (
Table 4,
Figure 5). Taken together, these data support the notion that neutrophilia would potentiate the impact of CRP on P/F, therefore paving the way to the need for ICU admission.
To the best of our knowledge, our study is therefore the first to demonstrate some of the pathways, and their direction, involved in the pathogenetic chain of respiratory failure in COVID-19 patients. Of note, we carried out this retrospective survey involving two different centers recruiting a selected cohort of hospitalized COVID-19 patients free of confounders affecting lymphocyte and neutrophil counts and/or serum CRP values.
Author Contributions
Conceptualization, L.M. and M.R. (Matteo Regolo); methodology, M.R. (Matteo Regolo), A.S. and L.M.; software, A.S.; validation, L.M., A.S. and M.R. (Matteo Regolo); formal analysis, A.S. and M.R. (Matteo Regolo); investigation, M.R. (Matteo Regolo), F.G. and D.F.; resources, L.M.; data curation, M.R. (Matteo Regolo), I.I., F.C., F.G., D.F., P.A., A.L., M.V., B.S., M.C., G.N., M.R. (Mario Russo), N.S. and M.M.; writing—original draft preparation, M.R. (Matteo Regolo) and L.M.; writing—review and editing, L.M. and M.R. (Matteo Regolo); visualization, M.R. (Matteo Regolo); supervision, L.M., P.N. and B.S.; project administration, L.M.; funding acquisition, L.M. All authors have read and agreed to the published version of the manuscript.