key: cord-0694165-wapgp9k8
authors: Besutti, Giulia; Giorgi Rossi, Paolo; Ottone, Marta; Spaggiari, Lucia; Canovi, Simone; Monelli, Filippo; Bonelli, Efrem; Fasano, Tommaso; Sverzellati, Nicola; Caruso, Andrea; Facciolongo, Nicola; Ghidoni, Giulia; Simonazzi, Anna; Iori, Mauro; Nitrosi, Andrea; Fugazzaro, Stefania; Costi, Stefania; Croci, Stefania; Teopompi, Elisabetta; Gallina, Annalisa; Massari, Marco; Dolci, Giovanni; Sampaolesi, Fabio; Pattacini, Pierpaolo; Salvarani, Carlo
title: Inflammatory burden and persistent CT lung abnormalities in COVID-19 patients
date: 2022-03-11
journal: Sci Rep
DOI: 10.1038/s41598-022-08026-1
sha: 646abeb4437741478a8c93baa559d330836315de
doc_id: 694165
cord_uid: wapgp9k8

Inflammatory burden is associated with COVID-19 severity and outcomes. Residual computed tomography (CT) lung abnormalities have been reported after COVID-19. The aim was to evaluate the association between inflammatory burden during COVID-19 and residual lung CT abnormalities collected on follow-up CT scans performed 2–3 and 6–7 months after COVID-19, in severe COVID-19 pneumonia survivors. C-reactive protein (CRP) curves describing inflammatory burden during the clinical course were built, and CRP peaks, velocities of increase, and integrals were calculated. Other putative determinants were age, sex, mechanical ventilation, lowest PaO2/FiO2 ratio, D-dimer peak, and length of hospital stay (LOS). Of the 259 included patients (median age 65 years; 30.5% females), 202 (78%) and 100 (38.6%) had residual, predominantly non-fibrotic, abnormalities at 2–3 and 6–7 months, respectively. In age- and sex-adjusted models, best CRP predictors for residual abnormalities were CRP peak (odds ratio [OR] for one standard deviation [SD] increase = 1.79; 95% confidence interval [CI] = 1.23–2.62) at 2–3 months and CRP integral (OR for one SD increase = 2.24; 95%CI = 1.53–3.28) at 6–7 months. Hence, inflammation is associated with short- and medium-term lung damage in COVID-19. Other severity measures, including mechanical ventilation and LOS, but not D-dimer, were mediators of the relationship between CRP and residual abnormalities.

.

Early studies on radiological follow-up of severe COVID-19 pneumonia showed residual abnormalities at follow-up CT scan, with a high prevalence especially in the short-and medium-term [4] [5] [6] [7] [8] . The clinical significance of these abnormal findings is unclear and the degree of persistency of residual pulmonary abnormalities will remain uncertain until more longer-term follow-up data are available [6] [7] [8] [9] .

In a study including 114 survivors of severe COVID-19 with six-month follow-up CT scan, residual groundglass opacities (GGO) or interstitial thickening, and CT fibrotic-like features (e.g., traction bronchiectasis, parenchymal bands, and/or honeycombing) were reported in 27%, and 35% of cases, respectively 5 . The presence of residual fibrotic-like abnormalities on CT was associated with older age, indices of higher disease severity, and mechanical ventilation 5 .

During the disease course, pulmonary and multiorgan damage is highly dependent on immune response and inflammation. In fact, inflammatory markers have been linked with poorer disease outcomes [10] [11] [12] . C-reactive protein (CRP) is the most widely used indicator of acute-phase response in clinical practice. CRP circulating concentrations increase rapidly and markedly in response to the onset of inflammation, and its levels decline with a half-life of approximately 20 h after resolution of the underlying inflammatory stimulus 13 . The rapid kinetics of circulating CRP explains why both CRP peak and the velocity of increase have been previously used in the definition of hyperinflammation in COVID-19 patients 12 .

The aim of this study was to evaluate the association between the degree of systemic inflammatory response during COVID-19 disease course and the presence of residual CT abnormalities at 2-3 and 6-7 months after COVID-19 pneumonia.

Setting. In the Reggio Emilia province (Northern Italy, 532,000 inhabitants) the first case of SARS-CoV-2 infection was diagnosed on 27 February 2020. The first pandemic wave lasted in Italy until May 2020. As of May 15, 2020, there had been 4863 RT-PCR-confirmed COVID-19 cases in the province. Details on hospital care in the province are reported in the Supplementary Materials. Study design and population. This was a single-centre retrospective cohort study based on routinely collected data, approved by the Area Vasta Emilia Nord (AVEN) Ethics Committee on 28 July 2020 (protocol number 855/2020/OSS/AUSLRE). The study was conducted in accordance with relevant guidelines and regulations, particularly with the Declaration of Helsinki. All consecutive COVID-19 patients who underwent a follow-up CT scan 2-3 months after COVID-19 pneumonia diagnosis were enrolled. In our institution, a routine 2-3 month follow-up CT scan was proposed to all COVID-19 survivors who had been hospitalised during the disease course for severe pneumonia with the following clinical-radiological features: respiratory failure during hospital stay (history of invasive or non-invasive mechanical ventilation and/or tocilizumab administration) or total extent of disease > 40% at baseline CT. A second CT scan at 6-7 months was proposed if residual CT abnormalities were present at 2-3 months. Given the retrospective nature of the cohort building, the need for informed consent was waived by the Area Vasta Emilia Nord Ethics Committee, after all reasonable efforts have been made to contact the patients to collect consent. Data collection. The COVID-19 Surveillance Registry, coordinated by the National Institute of Health and implemented in each Local Health Authority, was used to retrieve data on the dates of symptom onset, diagnosis, and hospitalisation. Patients' medical records were reviewed to collect data on comorbidities, COVID-19 treatment, and mechanical ventilation (invasive and non-invasive). Blood tests at Emergency Department presentation, and all the available levels of CRP, D-dimer, and PaO2/FiO2 measured during hospitalization were retrieved from the laboratory information system. Detailed measurement methods are reported in the Supplementary Materials. The visually estimated extent of disease was retrieved from the structured reports of baseline CT scans 14 .

Putative determinants. To study the relationship between inflammatory reaction and persistent lung abnormalities, three descriptors of CRP curves during disease course were considered: the peak (i.e., the maximum observed CRP value), the integral (i.e., the area under the curve of CRP levels), and the velocity of increase (i.e., the peak value divided by the number of days from symptom onset to the day on which the peak was registered). Details of the methods used to compute the CRP descriptors are reported in the Supplementary Materials and Supplementary Fig. 1 . Other known determinants of COVID-19 prognosis and residual lung abnormalities were considered as potential confounding factors and included in the adjusted models: age, sex, and history of chronic pulmonary conditions preceding COVID-19 pneumonia (current or previous smoking history and/ or concurrent respiratory diseases such as asthma, chronic obstructive pulmonary disease, and so on). Other www.nature.com/scientificreports/ indicators and markers of disease severity were considered as potential mediators of the inflammatory effect: mechanical ventilation, lowest PaO2/FiO2 ratio and D-dimer peak during the disease course, and length of hospital stay (LOS). Other known prognostic factors for COVID-19 which are intrinsically linked to the severity of inflammation as upstream causes or long-term effects of chronic inflammation (i.e., cardiovascular diseases, type 2 diabetes, and obesity) were not used to adjust the models because they probably lie on the same causal pathway linking inflammation to COVID-19 outcomes.

Materials. Two radiologists retrospectively reviewed follow-up CT scans at 2-3 and 6-7 months by consensus reading, and a gestalt assessment (i.e. global visual evaluation of the dominant CT pattern 15 ), comparing serial CT scans was jointly performed by the two radiologists. Three main CT evolution patterns were systematically assessed for both short-and medium-term follow-up CT scans: no residual or trivial CT abnormalities (e.g., subtle GGO and reticulations occupying < 5% lung parenchyma); residual CT abnormalities suggestive of organizing pneumonia (OP) sequelae without overt fibrotic features (e.g., GGO, perilobular opacities, parenchymal bands, bronchial dilatation, with no predominant subpleural distribution); residual CT abnormalities with fibrotic features suggestive of non-specific or usual interstitial pneumonia (GGO and/or reticulation with traction bronchiectasis and/or honeycombing). Moreover, the presence of CT findings typical of post-ventilatory damage (cicatricial emphysema, parenchymal bands, and traction bronchiectasis in the anterior portions of the upper lobes) was collected. CT extension of residual abnormalities at 2-3 and 6-7 months was collected as the visually-estimated percentage of lung parenchyma involved by residual abnormalities, and categorized in 4 groups: < 20%, 20-39%, 40-59%, and ≥ 60%. Early resolution was defined as no residual CT abnormality at 2-3 months, late resolution was considered for residual CT abnormalities (both non-fibrotic and fibrotic) at 2-3 months and no residual CT abnormalities at 6-7 months, while persistent CT abnormalities were defined when residual CT abnormalities (both non-fibrotic and fibrotic) were still visible at 6-7 months. Data analysis. We report continuous variables as median (interquartile range-IQR) and categorical variables as numbers and percentages. P-values are reported as continuous measures, and no significance threshold was set. Spearman' correlations between the three descriptors of CRP are reported.

In regression analysis, CRP peak, the logarithmic form of CRP velocity of increase, CRP integral, PaO2/ FiO2, D-dimer, and LOS were standardized. We built logistic models adjusted for age and sex to evaluate the association between CRP curve descriptors and the presence of CT lung abnormalities at 2-3 and 6-7 months. The performances of models with different descriptors were compared using Pseudo R 2 , Akaike information criterion (AIC), and area under ROC curve (AUROC), and the models with the best values of these three statistics were used to identify the CRP descriptors for subsequent analyses. We also present models adjusted by age, sex and including one by one the other COVID-19 severity indicators, i.e. mechanical ventilation and standardized D-dimer peak, lowest PaO2/FiO2, and LOS. Sensitivity analyses were also conducted on the selected models excluding patients with post-ventilatory lung CT abnormalities and patients with reported active smoking habit to exclude that the associations were only linked to these conditions. Linear regression models adjusted for age, sex only and for age, sex, and baseline CT extension, were used to test the association of the selected CRP descriptors and the extension of residual abnormalities at 2-3 and 6-7 months.

To evaluate the role of other measures of disease severity, we performed a formal mediation analysis 16,17 using mechanical ventilation (invasive and non-invasive), and the standardized variables D-dimer peak, lowest PaO2/FiO2, and LOS as putative mediators of the CRP descriptors that were selected to be the best predictors of persistent lung abnormalities; models were always adjusted for age and sex. Details on mediation analysis are reported in the Supplementary Materials.

Given the different recommendations on mechanical ventilation in older patients, we tested for an interaction between age and PaO2/FiO2 or mechanical ventilation in their effect on outcomes.

Mediation analysis was performed using "medflex" package in R 3.4.0; the other analyses were conducted using Stata 13.0 SE (Stata Corporation, Texas, TX).

February -15 May 2020), 1513 residents were hospitalised due to COVID-19. Of these, 361 patients died. Among the 1152 survivors, 194 patients experienced a respiratory failure during hospitalisation, and 157 other patients had a baseline total CT extent of ≥ 40% (Fig. 1) . Of the patients who had an indication for follow-up CT scan, 259 patients (median age 65 years, 30.5% females) underwent a follow-up CT scan at 2-3 months and were included in the study; 241 of them had a second follow-up CT scan at 6-7 months.

Clinical characteristics of the included patients are reported in Table 1 and Supplementary Table 1 . CRP was measured on average every 3 days (79 h), for a total average of 7.4 times during the clinical course of disease. All CRP descriptors were higher in patients with persistent abnormal findings ( Table 1 ). The distribution of CRP curve descriptors in the whole population are described in Supplementary Fig. 2 . CRP peak was strongly correlated with CRP velocity of increase and CRP integral, while the correlation between CRP integral and velocity of increase was weaker (Supplementary Table 2 6% of the initial sample) still had residual abnormalities at 6-7 months (Fig. 2) . Most of the patients showed non-fibrotic CT abnormalities suggestive of OP sequelae ( Fig. 3 and Supplementary Fig. 3 ). Only 9/100 patients showed fibrotic abnormalities at 6-7 months (Fig. 4) , while 7 had CT findings typical of ventilation-induced damage (Fig. 5 ).

CRP predictors of lung CT abnormalities at 2-3 and 6-7 months. By building different models adjusted for age and sex, we evaluated different combinations of standardized CRP curve descriptors in the prediction of persistent lung CT abnormalities. Other considered putative determinants, such as asthma or COPD and/or active smoking habit, were also tested as covariates, but they provided no added value to the models. The most accurate model for CT abnormalities at 2-3 months included CRP peak only (AUROC = 0.752, R 2 = 0.138), while for residual abnormalities at 6-7 months the most accurate model included CRP integral only (AUROC = 0.748, R 2 = 0.141) (Supplementary Table 3 ). In particular, a CRP peak increase of one SD was associated with an 80% increase in the risk of having residual abnormalities at 2-3 months (OR = 1.79; 95% CI = 1.23-2.62), while one SD increase in CRP integral was associated with a two-fold increase in the risk of having residual Table 4 ). In linear regression models, CRP peak and CRP integral were associated with CT extension of residual abnormalities at 2-3 and 6-7 months, respectively, both when adjusting for age and sex only and when adjusting also for baseline CT extension (Supplementary Table 5 ).

To test whether other indicators of disease severity not directly linked to hyperinflammation may have an independent role in predicting persistent lung CT abnormalities, we added other disease severity measures to the previously selected models (Table 2 ). An interaction between mechanical ventilation, PaO2/FiO2, and age in their effect on the outcomes was excluded. Only a slight increase in the models' performance was observed when adding LOS or the lowest Pa/FiO2, while when adding mechanical ventilation, the AUC of the models improved more noticeably (from 0.752 to 0.824 for residual abnormalities at 2-3 months and from 0.748 to 0.780 for residual abnormalities at 6-7 months). Model performance was not directly comparable when adding D-dimer peak due to the presence of a remarkable number of missing values (n = 63). However, the increase in R 2 suggests that D-dimer peak also adds a predictive value. www.nature.com/scientificreports/ Adding other measures of disease severity partially reduced the association between CRP measures and the outcomes. This effect was absent or scarce for D-dimer and more apparent for the other measures of disease severity, except for the lowest PaO2/FiO2 in the model for 6-7-month abnormalities. The decrease in the association with the outcome is stronger for CRP peak when adding mechanical ventilation and for CRP integral when adding LOS.

The mediation analysis confirmed that mechanical ventilation and LOS were mediators of the effect of CRP descriptors on persistent CT abnormalities (Fig. 6) . For example, 42% of the effect of CRP peak on short-term abnormalities and 23% of the effect of CRP integral on medium-term abnormalities was due to mechanical ventilation. A similar effect was found for the lowest PaO2/FiO2 ratio when considering abnormalities at 2-3 months. Instead, when considering abnormalities persisting at 6-7 months, PaO2/FiO2 was not associated with the outcome and therefore it could not be a mediator of CRP effect. Finally, the mediation effect of D-dimer peak on outcomes was negligible and could have been due to random fluctuations (the proportion of exposure effect due to the mediator on CT lung abnormalities at 2-3 months and at 6-7 months was 16% and 15%, respectively). www.nature.com/scientificreports/ Figure 6 . Path diagrams representing the mediation effect of mechanical ventilation, lowest PaO2/FiO2, D-dimer peak, and length of hospital stay in the relationship between CRP peak and CT lung abnormalities at 2-3 months and between CRP integral and CT lung abnormalities at 6-7 months. θ1 is the estimated effect (Ln OR) of CRP descriptors on mechanical ventilation adjusting for confounders. θ2 is the estimated effect (Ln OR) of mediators on persistent lung abnormalities adjusting for confounders. β is the correlation coefficient of CRP descriptors on mediators. Confounders include age and sex. In the tables we report the decomposition of the effect of CRP descriptors on persistent lung abnormalities. OR, odds ratio; TE, total effect; NDE, natural direct effect; NIE, natural indirect effect; NIE/TE, % mediation proportion. 

In a population of COVID-19 patients who were hospitalised for severe pneumonia, lung CT abnormalities were still visible in almost 80% patients after 2-3 months from disease onset, and in almost 40% after 6-7 months. CRP as a measure of systemic inflammatory response predicted residual CT abnormalities. The best predictors were CRP peak for residual abnormalities at 2-3 months and CRP integral for persistent abnormalities at 6-7 months. Other measures of disease severity, such as LOS, the lowest PaO2/FiO2, and the need for mechanical ventilation were mediators in the relationship between CRP and persistent CT lung abnormalities. While D-dimer was associated with persistent lung abnormalities, it was not a mediator of the inflammation effect.

Our results are consistent with prior follow-up studies, which reported late resolution or residual abnormalities at 6 months in up to 72% of subjects who had had COVID-19 pneumonia [4] [5] [6] [7] . Most patients with residual CT abnormalities did not have overt CT features of lung fibrosis. Furthermore, more than 40% of patients who showed residual abnormalities at 2-3-month achieved complete resolution at 6-7 months. These findings corroborate the hypothesis that residual CT abnormalities are the expression of slow recovery in most cases, despite data on longer-term follow-up being needed to confirm that persistent fibrotic changes are not common. In fact, non-fibrotic abnormalities could be the expression of immature fibrosis (i.e. "fibroblastic changes" due to diffuse alveolar damage) remodelling 9 .

In the few available studies on the topic, several measures of disease severity predicted the presence of persistent lung abnormalities in COVID-19 patients 4,5 , but the role of inflammatory burden throughout the disease course has never been investigated.

Baseline CRP has been widely studied as a predictor of COVID-19 severity in the acute phase of the disease [10] [11] [12] . Some studies have also evaluated CRP trajectories during hospitalisation, showing that rapidly rising CRP levels may predict respiratory failure and that CRP peak is higher in patients who require intubation or die [12] [13] [14] [15] [16] [17] [18] . We show that persistent lung abnormalities are also linked to inflammatory burden, with most important predictors being CRP peak for persistent damage in a subacute phase (2-3 months), and CRP integral, reflecting both the intensity and the duration of the inflammatory reaction, for persistent abnormalities at 6-7 months.

Given the known association between CRP and lung damage in the acute phase 19 , it is not surprising that much of the effect of inflammatory burden on persistent lung abnormalities is mediated by proxies of lung failure severity, such as the need for mechanical ventilation and, at least for residual abnormalities at 2-3 months, the lowest PaO2/FiO2 as a measure of hypoxemia. The association between mechanical ventilation and other disease severity measures, including CRP, should be considered cautiously in our cohort, since the presence of mechanical ventilation was one of the sufficient, but not necessary inclusion criteria. Surprisingly, we did not observe an association between the lowest PaO2/FiO2 and 6-7 month abnormalities, therefore we did not find any mediation effect for this measure of disease severity.

Higher D-dimer levels have been widely associated with COVID-19 severity and mortality [20] [21] [22] , as the infection induces a coagulopathy with secondary hyperfibrinolysis in severe cases 23 . In our study, D-dimer peak was associated with persistent lung abnormalities, adding predictive value to the models without subtracting from the effect of CRP. Accordingly, its role as a mediator of the effect of inflammation on persistent abnormalities was modest, if any. Apparently, even if inflammation and coagulopathy are linked in COVID-19 patients 24-26 , CRP and D-dimer effects on persistent lung damage are independent, acting through different causal pathways.

This study has some limitations. First, only COVID-19 survivors who experienced severe pneumonia (extension > 40% at hospital admission and/or respiratory failure during hospitalisation) were routinely evaluated with follow-up CT scans and included in the study. Consequently, selection, survivorship and colliding bias cannot be excluded, and some of the observed associations between the severity measures themselves and between severity measures and the outcomes may be under or over-estimated. Mediation analysis should reduce the impact of colliding bias; while it must be clear that our findings are only meaningful for survivors of COVID-19 pneumonia identified as severe by baseline imaging or respiratory failure occurrence. Second, CRP curves were not always complete, and extrapolations were done. Many values were missing for D-dimer. CRP curve may be influenced by factors occurring during the disease course, such as bacterial superinfection, anti-inflammatory drugs, and post-ventilatory lung damage, but also by pre-existing pathological conditions. In our study, we can only observe the association between the inflammatory burden during COVID-19 and residual lung abnormalities, whatever was the mechanism underlying the inflammation, i.e. the direct response to the virus or rather an inflammatory response whose intensity can be modulated by therapies or by other pre-existing conditions. In fact, these factors cannot be considered confounders of our association, but factors probably lying on the same causal pathway, that starts from the existence of prognostic factors favoring an intense inflammatory response when SARS-CoV-2 infection occurs, and then includes the immune response following infection and the therapies we administered in response to infection and disease progression. Furthermore, we had no information on other potential confounders, such as pre-existing interstitial lung abnormalities. Nevertheless, our model building strategy was guided by the attempt to describe causal pathways, and our data suggest that therapy-linked variables, i.e. mechanical ventilation and LOS, are on the same causal pathway linking inflammation to persistent lung abnormalities. Our findings suggest that future studies trying to understand the etiopathogenesis of COVID-19 should carefully consider which variables should be placed in multivariable analyses. Otherwise, the risk is of hiding the causal links with upstream determinants, such as hyperinflammation markers, and of introducing effect mediators. Moreover, our study focused on persistent morphological abnormalities evaluated by means of CT scan, but the evaluation of the impact of residual CT abnormalities on symptoms and lung function tests is ongoing and is beyond the scope of the present study. Finally, our attempt to classify the pattern of CT abnormalities may be questionable. In fact, the non-fibrotic CT pattern may include features (e.g., GGO), that require more time to become overtly fibrotic on CT 9 . Also, the use of quantitative extension of residual changes does not necessarily improve the evaluation of CT findings evolution, which is variable and may include diffuse but www.nature.com/scientificreports/ very mild abnormalities (e.g. subtle GGO). However, the association of CRP with residual abnormalities was confirmed also by using CT visual extension as outcome measure.

In conclusion, higher inflammatory burden is associated with short-and medium-term lung damage. The peak value of CRP is the most important determinant of damage at 2-3 months, while both the duration and intensity of the inflammatory response become important for more long-lasting damage (6-7 months). Other measures of disease severity seem to be mediators of the inflammatory effect, except for D-dimer peak, confirming that the two pathways of inflammation and coagulopathy are partially independent, although both contribute to lung damage.

Participant data that underlie the results reported in this manuscript will be shared after de-identification, beginning 6 months and ending at least 7 years after article publication, to researchers who provide a methodologically sound proposal with objectives consistent with those of the original study. Proposals and data access requests should be directed to the Area Vasta Emilia Nord (AVEN) Ethics Committee at CEReggioemilia@ausl.re.it as well as to the Authors at the Epidemiology Unit of AUSL-IRCCS di Reggio Emilia at info.epi@ausl.re.it, who are the data guardians. To gain access, data requestors will need to sign a data access agreement. 

COVID-19) Dashboard

An interactive web-based dashboard to track COVID-19 in real time

Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a singlecentered, retrospective, observational study

Radiological and functional lung sequelae of COVID-19: a systematic review and meta-analysis

Six-month Follow-up Chest CT Findings after Severe COVID-19 Pneumonia

3-month, 6-month, 9-month, and 12-month respiratory outcomes in patients following COVID-19-related hospitalisation: a prospective study

Postacute Sequelae of COVID-19 Pneumonia: 6-month Chest CT Follow-up

Eight months follow-up study on pulmonary function, lung radiographic, and related physiological characteristics in COVID-19 survivors

Interstitial lung disease after COVID-19 infection: a catalog of uncertainties

Inflammatory and hematologic markers as predictors of severe outcomes in COVID-19 infection: A systematic review and meta-analysis

The role of cytokine profile and lymphocyte subsets in the severity of coronavirus disease 2019 (COVID-19): a systematic review and meta-analysis

COVID-19-associated hyperinflammation and escalation of patient care: a retrospective longitudinal cohort study

How C-reactive protein structural isoforms with distinctive bioactivities affect disease progression

Accuracy of CT in a cohort of symptomatic patients with suspected COVID-19 pneumonia during the outbreak peak in Italy

Gestalt theory: implications for radiology education

Odds ratios for mediation analysis for a dichotomous outcome

Assessing natural direct and indirect effects for a continuous exposure and a dichotomous outcome

Inflammatory biomarker trends predict respiratory decline in COVID-19 patients

The association between clinical laboratory data and chest CT findings explains disease severity in a large Italian cohort of COVID-19 patients

D-dimer is associated with severity of coronavirus disease 2019: a pooled analysis

Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia

Immunity, endothelial injury and complement-induced coagulopathy in COVID-19

Coagulopathy of coronavirus disease 2019

D-dimer correlateswith proinflammatory cytokine levels and outcomes in critically ill patients

Pathological findings of COVID-19 associated with acute respiratory distress syndrome

Pathogenic human coronavirus infections: causes and consequences of cytokine storm and immunopathology

were responsible for the study conceptualization

We thank Jacqueline Costa for the English language editing.

This work was supported by the Italian Ministry of Health, grant no. COVID-2020-12371808. The funder had no role in the definition of the study design, data collection and analysis, interpretation of the results or preparation of the manuscript.

PGR conducted negotiations with Roche, Hologic and Becton Dickinson, as the PI of an independent study funded by the Italian Ministry of Health, to obtain reagents at reduced price or for free. All other Authors have no conflict of interests do declare.

Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-08026-1.Correspondence and requests for materials should be addressed to G.B.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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:// creat iveco mmons. org/ licen ses/ by/4. 0/.