Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Persistent pulmonary pathology after COVID-19 is associated with high viral load, weak antibody response, and high levels of matrix metalloproteinase-9

Abstract

The association between pulmonary sequelae and markers of disease severity, as well as pro-fibrotic mediators, were studied in 108 patients 3 months after hospital admission for COVID-19. The COPD assessment test (CAT-score), spirometry, diffusion capacity of the lungs (DLCO), and chest-CT were performed at 23 Norwegian hospitals included in the NOR-SOLIDARITY trial, an open-labelled, randomised clinical trial, investigating the efficacy of remdesivir and hydroxychloroquine (HCQ). Thirty-eight percent had a CAT-score ≥ 10. DLCO was below the lower limit of normal in 29.6%. Ground-glass opacities were present in 39.8% on chest-CT, parenchymal bands were found in 41.7%. At admission, low pO2/FiO2 ratio, ICU treatment, high viral load, and low antibody levels, were predictors of a poorer pulmonary outcome after 3 months. High levels of matrix metalloproteinase (MMP)-9 during hospitalisation and at 3 months were associated with persistent CT-findings. Except for a negative effect of remdesivir on CAT-score, we found no effect of remdesivir or HCQ on long-term pulmonary outcomes. Three months after hospital admission for COVID-19, a high prevalence of respiratory symptoms, reduced DLCO, and persistent CT-findings was observed. Low pO2/FiO2 ratio, ICU-admission, high viral load, low antibody levels, and high levels of MMP-9 were associated with a worse pulmonary outcome.

Introduction

The abrupt spread of SARS-CoV-2, leading to the COVID-19 pandemic, has had a devastating impact on the entire world’s population1. With the lungs as the primary target, knowledge about the prevalence and causes of long-term pulmonary sequelae in COVID-19 survivors is paramount.

We recently reported that 50% of COVID-19 survivors struggle with dyspnoea 3 months following hospital admission, with 25% presenting with reduced diffusion capacity of the lungs for carbon monoxide (DLCO), and 20% having signs of early pulmonary fibrosis on chest computed tomography (CT)2. Similar results have been reported by others3. With various drugs emerging as possible therapeutic candidates against SARS-CoV-2, long-term pulmonary function following treatment and potential targets for preventing pulmonary sequelae warrant elucidation4.

The NOR-SOLIDARITY trial is an independent add-on study to the WHO SOLIDARITY trial, evaluating the effects of remdesivir and hydroxychloroquine (HCQ) compared to standard of care (SoC) in hospitalised COVID-19 patients, assessed by in-hospital mortality, the need for mechanical ventilation and the length of hospital stays5. To explore end-organ deterioration following the infection, the NOR-SOLIDARITY trial included biobanking during hospital stay and at a 3-month follow-up visit6.

Studies of patients hospitalised for SARS, caused by another beta coronavirus (SARS-CoV) with high similarity to SARS-CoV-2, have highlighted the role of an overwhelming host response in mediating lung injury and subsequent development of pulmonary fibrosis7. In other forms of pulmonary fibrosis, several mediators are under investigation as markers and potential mediators of disease severity and prognosis, such as various matrix metalloproteinases (MMP), surfactant protein D (SP-D), and vascular endothelial growth factor (VEGF) A, as a key regulator of angiogenesis8,9,10. Some of these mediators have also been reported to be elevated in the acute phase of COVID-1911,12. However, the association between these regulators of extracellular matrix remodelling and the development of persistent lung parenchyma pathology and reduced pulmonary function following COVID-19, has not been fully elucidated.

We aimed to investigate pulmonary function through COPD Assessment Test Score (CAT-score), spirometry, DLCO, and chest CT 3 months after hospital admission for COVID-19, and whether there was a difference in these variables between participants treated with remdesivir, HCQ, or local SoC alone. As explorative outcomes, we evaluated associations between clinical characteristics during hospital admission, including viral load in oropharynx and the presence of antibodies against SARS-CoV-2, with the pulmonary outcomes. In addition, we examined if plasma levels of factors suspected to be involved in development of pulmonary fibrosis (i.e., MMP-9, SP-D, and VEGF-A), obtained during hospital admission and after 3 months, were related to patient-reported outcome measures, as well as pulmonary function and persistent findings on chest CT.

Materials and methods

Design

The NOR-SOLIDARITY trial is an adaptive, multicentred, randomised, open-labelled clinical trial based on the WHO trial SOLIDARITY, where the target was to evaluate the efficacy and safety of remdesivir and HCQ versus their respective SoC group in hospital-admitted COVID-19 patients. A 3-month follow-up visit was organized to evaluate lung function by spirometry, DLCO, and chest CT. Peripheral blood was collected both during hospital stay and at the 3-month follow up. Biobanking material from oropharynx was collected during hospital admission. Recent work by the SOLIDARITY and NOR-SOLIDARITY consortia describes in detail the design of the trial, the recruitment of participants, the randomisation process, interventions, and secondary outcomes5,6. The study was approved by the Committee for Medical Research Ethics Region South East Norway (REK 118684) and registered on ClinicalTrials.gov on the 25.03.2020 (NCT04321616). Monitors from the Clinical Trial Unit at Oslo University Hospital monitored the conduct of the trial. Good clinical practice was followed throughout the study. All participants gave informed consent prior to inclusion, either themselves or by a legally authorised representative.

Outcomes

The outcomes of this sub-study were CAT-score, spirometry, DLCO, and chest CT-findings 3 months following hospital admission for COVID-19, and the effect of remdesivir and HCQ versus their respective SoC on these variables. Explorative outcomes were to correlate clinical characteristics during hospital admission with the pulmonary outcomes. These included pO2/FiO2 (P/F)-ratio, respiratory rate (RR), admission to an intensive care unit (ICU), the receipt of mechanical ventilation, as well as plasma levels of C-reactive protein (CRP), lactate dehydrogenase (LDH), D-dimer, ferritin, viral load in oropharynx, and antibodies against the receptor binding domain (RBD) of SARS-CoV-213. Secondary explorative outcomes were to correlate the levels of MMP-9, SP-D, and VEGF-A during hospital admission and at 3 months following discharge for COVID-19, with the pulmonary outcomes.

Participants

Twenty-three Norwegian hospitals participated in the NOR-SOLIDARITY trial. Adult patients (≥ 18 years) admitted to hospital with a PCR-confirmed SARS-CoV-2 infection from March 28th until October 5th 2020 were eligible for participation. The exclusion criteria are described in detail in a recently published report6. For the purpose of this report, participants having completed the 3-month follow-up visit by the 15th of October 2020 were included (n = 108). Participants were randomised to receive remdesivir, HCQ or SoC. During the trial period, the availability of the drugs varied, implying that some participants served as SoC control subjects for one drug, while others served as SoC control subjects for both remdesivir and HCQ. This might give the impression that there are more than 108 participants included. An overview of missing parameters (i.e., CT scans, blood analyses, and pulmonary function tests) are presented in Supplementary Table S1.

Intervention

In the NOR-SOLIDARITY trial, the participants were allocated to one of three treatment arms; SoC, remdesivir (200 mg day 1, then 100 mg daily up to nine days), or HCQ (800 mg twice daily day 1, then 400 mg twice daily up to nine days). All trial regiments were discontinued at discharge. Unintended effects are described in a previous report6.

Patient-reported outcome measures

CAT-score is an eight-item questionnaire designed to assess and quantify the impact of COPD symptoms on health status on a scale from 0 (low) to 40 (very high)14. The eight items cover respiratory symptoms (cough, mucus, tightness in chest) and impaired quality of life (limitations in daily activities, sleep impairment, loss of energy). Scores ≥ 10, indicating a moderate to very high symptom impact on health status, were used as a cut-off15. Participants responded to the questionnaire at the 3-month follow-up visit.

Pulmonary function tests

Pulmonary function tests consisted of spirometry (Jaeger, Höechberg, Germany and Carefusion, Yorba Linda, CA, USA) and DLCO, and were performed according to the European Respiratory Society’s and the American Thoracic Society’s guidelines, thoroughly described in detail in a previously published report2,16,17,18,19.

Acquisition and review of chest CT images

Low-dose, thin-section chest CT images were obtained in supine and prone position during breath-holding in deep inspiration. Supplementary expiratory scans were obtained in 18 (17.6%) participants to differentiate between air trapping and hypoperfusion in areas with decreased attenuation on inspiratory scans. The CT protocol and method of interpretation according to the Fleischner Society, are described in detail in a previously published report2,20,21. For the purpose of this report, we assessed the prevalence of any ground-glass opacities (GGO), GGO ≥ 10% in at least one of the four lung zones, or any mosaic pattern (classified as potentially reversible changes and grouped together), any consolidations, reticular pattern, parenchymal bands, interlobular septal thickening, or any bronchiectasis (classified as irreversible changes and grouped together).

Routine laboratory analyses and measurements of SARS-COV-2 viral load and antibodies

Routine blood samples (CRP, LDH, ferritin, D-dimer, white blood cell counts, the number of neutrophils, lymphocytes, and platelets, and haemoglobin levels) were analysed by the routine laboratory at each participating centre. Viral load in oropharynx was analysed by RT-PCR, and SARS-CoV-2 receptor binding domain (RBD) IgG antibodies were detected by multiplexed bead-based flow cytometric assay, referred to as microsphere affinity proteomics as previously described6.

Measurement of MMP-9, tissue inhibitor of metalloproteinase (TIMP) 2, SP-D, AND VEGF-A

Per protocol, biobanking was performed during hospital admission and at the 3-month follow up visit. Plasma was collected in sterile EDTA containing tubes, put on melting ice and centrifuged at 2000×g for 20 min to obtain platelet-poor plasma, before storage at − 80 °C. The samples were thawed only once. Plasma levels of MMP-9, TIMP-2, SP-D, and VEGF-A were measured in duplicate by enzyme immunoassays (EIA) using commercially available reagents (Cat#DY911, DY971, DY1920, DY293B; RnD systems, Stillwater, MN, USA) in a 384-well format. Absorption was read at 450 nm with wavelength correction set to 540 nm using an EIA plate reader (BioTek). Intra- and inter-assay coefficients of variation were < 10% for all EIAs. Levels of MMP-9, TIMP-2, SP-D, and VEGF-A in the plasma of 24 healthy blood donors (15 (62.5%) males, mean age 44 years, one (4.2%) former smoker, one (4.2%) taking statins), were also measured to establish a point of reference.

Statistics

The primary outcome of the study design has already been published, and analyses in this report are to be considered as secondary and explorative analyses5,6. Demographic data were described by median (25th and 75th percentile (IQR)) and percentages. CAT-score, pulmonary function tests, and CT-findings at the 3-months follow-up are presented as median values and IQR. We assessed the difference in pulmonary function, defined by FEV1 and DLCO in percent of predicted (DLCO%) as continuous variables, CAT-score, and CT-findings between the intervention groups and their respective SoC-group 3 months following treatment for COVID-19 by 95% confidence intervals from the two-sample t-distribution, or p-value from the 2 × 2 Chi-squared distribution. To evaluate participants’ symptom load to pulmonary function, we associated CAT-score ≥ 10 with DLCO%, and the presence of any reversible, and any irreversible CT-findings using logistic regression.

Associations between clinical characteristics during hospital stay and pulmonary outcomes after 3 months, were performed by logistic and linear regression, and adjusted by age, sex, treatment, and history of smoking. As an explorative outcome we wanted to study whether these clinical characteristics could predict pulmonary dysfunction after 3 months, defined by CAT-score ≥ 10, DLCO%, and the presence of either GGO in more than ≥ 10% of at least one of four lung zones as a measure of possible reversible changes, and irreversible changes on chest CT.

Further explorative outcomes were associations between the maximum plasma levels of MMP-9, TIMP-2, SP-D, and VEGF-A during hospital stay and at the 3-month follow-up visit, with DLCO%, and the predefined changes on chest CT. Comparisons of pro-fibrotic mediators between the healthy control group, the maximum levels during admission and at the 3-months follow-up, were performed with Mann–Whitney U test and Wilcoxon paired test. Statistical analyses were performed using the statistical software R (version 4.0.2, 2020.06.22, R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/) and SPSS (IBM Corp. released 2019. IBM SPSS Statistics for Windows, Version 26.0. Armonk, NY: IBM Corp). Figures were produced using Prism 9 (GraphPad Software, LCC, version 9.1.0., San Diego, California, USA).

Results

Participant flow

Hospitalised COVID-19 patients were included from 23 Norwegian hospitals from the 28th of March 2020. For the purpose of this paper, the data collection from the 3-month follow up visit was ended on the 15th of October, 2020. As demonstrated in Fig. 1, 108 participants completed the 3-month follow-up visit, with a median time of 93.5 (IQR 88, 106) days from hospital admission. Missing parameters are presented in Supplementary Table S1.

Figure 1
figure 1

Flow chart of the participant flow. *Other: Emigration, progression of cancer diseases, general health deterioration, prolonged corrected QT-time during treatment with HCQ; Excluded from the full analysis set. HCQ hydroxychloroquine.

Demographics

Demographic and clinical characteristics from the hospital admission for all participants and treatment groups with their respective SoC-group are presented in Table 1. Overall, the median age was 58 years (range 27.4, 89.8), 41 (38%) were female, few were current smokers (1.9%), and the median body mass index (BMI) was 27.6 (IQR 25.3, 30.7). Of comorbidities, hypertension (31.5%), diabetes mellitus (15.7%), and obesity (25.9%) were the most common. Overall, 40 (37.0%) presented with respiratory failure at inclusion, defined as P/F-ratio ≤ 40 kPa. The median duration of symptoms prior to hospital admission was eight days (IQR 5, 10).

Table 1 Characteristics of population at baseline.

CAT-score, pulmonary function tests, and chest ct-findings 3 months following treatment for COVID-19

Table 2 presents CAT-score, pulmonary function tests, and all registered CT-findings 3 months following hospital admission for the total COVID-19 population, as well as the distribution between the intervention groups and their respective SoC groups. As demonstrated in Table 3, there was a higher prevalence of self-reported symptoms in the remdesivir group than in the corresponding SoC group. Otherwise, there were no significant differences between the intervention groups and their respective SoC groups with regard to pulmonary function tests or CT-findings 3 months following inclusion. Forty-one (38.0%) participants had a CAT-score ≥ 10, indicating a moderate to high symptom burden 3 months following hospitalisation (Table 2), with dyspnoea on exertion and loss of energy being the most prominent symptoms (Fig. 2).

Table 2 Presentation of lung function and CT findings 3 months after treatment for COVID-19.
Table 3 Effect of medication on patient reported outcomes, lung function and CT-findings 3 months after intervention.
Figure 2
figure 2

The CAT-score results for the overall population 3 months following hospital admission for COVID-19. CAT-score is an eight-item questionnaire designed to assess and quantify the impact of COPD symptoms on health status. The items comprise symptoms of chough, phlegm, feeling of chest tightness, breathlessness, limitations in activities, confidence, sleep problems, and energy. Each item ranging from 0 = no symptoms, to 5 = high symptom burden. CAT-score COPD Assessment Test score, COPD chronic obstructive pulmonary disease.

For the pulmonary function tests, the majority of participants had values within the limits of normal. For FVC and FEV1, 13 (12%) and 15 (13.9%) of the overall population had values below LLN, respectively. Thirty-two (29.6%) had DLCO below LLN.

The most dominant finding on chest CT was parenchymal bands, a potential sign of pulmonary fibrosis, present in 45 (41.7%) participants. GGO, suggesting sustained inflammation, were also a prevalent finding in the total study population; 43 participants (39.8%). GGO present in ≥ 10% in one of the four lung zones, were found in 19 (17.6%) participants in the overall population.

CAT-score associated with DLCO% and chest CT-findings

No association between pulmonary-related symptomatology at 3 months, determined by CAT-score, and pulmonary function tests and pathology on chest CT at the same point, was observed (Supplementary Table S2).

Association between clinical characteristics during hospital admission and pulmonary outcomes

Table 4 presents explorative associations between demographic, clinical and biochemical characteristics at baseline and the pulmonary outcomes after 3 months, adjusted for age, sex, treatment, and a history of smoking.

Table 4 Association of respiratory symptoms, DLCO, ground glass opacities, and irreversible CT-findings, with clinical and demographic characteristics during hospital admission.

CAT-score ≥ 10 demonstrated a significant positive correlation with a history of smoking. Three-months’ values of DLCO%, the most frequently affected pulmonary function test variable, were positively correlated with P/F-ratio, negatively correlated with admission to ICU, negatively correlated with viral load at baseline, and positively correlated with the presence of IgG antibodies against SARS-CoV-2 RBD at baseline. Both the presence of GGO in ≥ 10% in at least one of the lung zones and the presence of irreversible CT-findings at 3 months were negatively correlated with female sex and positively correlated with age. Notably, we found no association between markers of inflammation during hospital admission (i.e., CRP, LDH, D dimer, and ferritin) and pulmonary outcomes at 3 months (Table 4).

Association of pro-fibrotic mediators with DLCO% and chest CT-findings

The maximum levels of both MMP-9 and VEGF-A during the hospital stay were significantly higher among participants, than in the healthy control group (Fig. 3). All mediators showed a significant decrease from the acute phase to the 3-month follow-up. At the 3-month follow-up visit, only VEGF-A levels were significantly higher than in the healthy control subjects. In this explorative part of the study, a significant correlation between maximum levels of MMP-9 during hospital stay and the presence of mosaic pattern on chest CT, as well as parenchymal bands, was observed (Table 5). The levels of TIMP-2, a known inhibitor of MMP-9 activity, were not increased during the study period, compared to healthy control subjects (Fig. 3).

Figure 3
figure 3

Group comparison of pro-fibrotic markers between 24 healthy control subjects, and all participants’ maximum plasma levels during hospital stay and at the 3-months follow-up evaluation, performed with Mann–Whitney U test and Wilcoxon paired test. Controls, 24 healthy control subjects. Maximum, maximum plasma levels during the hospital stay. 3 months, plasma levels at the 3-month follow-up. MMP-9 matrix metalloproteinase 9, TIMP-2 tissue inhibitor of metalloproteinases 2, SP-D surfactant protein D, VEGF-A vascular endothelial growth factor. Median MMP-9 level healthy control subjects: 2.70 ng/mL, median maximum level MMP-9: 3.83 ng/mL, median MMP-9 level at 3 months: 2.93 ng/mL. Median TIMP-2 level healthy control subjects: 7.78 ng/mL, median maximum level TIMP-2: 7.85 ng/mL, median TIMP-2 level at 3 months: 8.76 ng/mL. Median SP-D level healthy control subjects: 20.9 ng/mL, median maximum level SP-D 16.4, median SP-D level at 3 months: 14.1 ng/mL. Median VEGF-A level healthy control subjects: 125 pg/mL, median maximum VEGF-A level: 277 pg/mL, median VEGF-A level at 3 months: 216 pg/mL. *p < 0.05, **p < 0.01, ***p < 0.001.

Table 5 Association of pro-fibrotic markers during hospital admission and after 3 months, with CT findings 3 months after hospital admission for COVID-19.

Discussion

In this sub-study of the NOR-SOLIDARITY trial, nearly 40% of the participants reported persistent airway symptoms, and nearly 30% had DLCO below LLN 3 months following hospital admission and treatment for COVID-19. Findings on chest CT were also common, with 40% presenting with GGO and a similar percentage with parenchymal bands 3 months after hospital admission. We did not find significant differences between participants randomised to remdesivir and HCQ, and their respective SoC, with regard to pulmonary function and CT-findings, except for a somewhat higher symptom burden in the remdesivir group compared to the SoC group. Of the clinical characteristics during hospitalisation, a low P/F-ratio, the need for ICU admission, a high viral load, and low levels of IgG antibodies to SARS-CoV-2 seemed to predict reduced DLCO% at 3 months. Finally, high plasma levels of MMP-9, both during hospital stay and 3 months after hospitalisation, were associated with signs of both reversible and irreversible changes on chest CT.

The majority of our study participants were men, and the most common comorbidities were hypertension, diabetes mellitus, and obesity, which is consistent with findings in other studies22. Former smoking was relatively common, but very few were current smokers. The participants had a duration of symptoms of 8 days prior to admission, and nearly 40% tested positive for SARS-CoV-2 IgG antibodies at admission, reflecting the natural course of COVID-19 infection described by others23. Nearly 40% of the participants presented with respiratory failure, defined by P/F-ratio < 40 kPa, and the majority had moderately elevated inflammatory markers, which is in accordance with studies from other countries24,25,26. Although previously described as a marker of COVID-19 severity, lymphopenia was not a common trait in our population at baseline27. Kåsine et al. recently characterised the leukocyte subpopulations in our participants, demonstrating neutrophilia as more dominant than lymphopenia during the first week of hospitalisation in our cohort, along with an association between high neutrophil counts and disease severity28.

In accordance with the WHO SOLIDARITY study and the NOR-SOLIDARITY sub study demonstrating no effect on in-hospital mortality, the need for ICU treatment, or P/F ratio, we show that treatment with remdesivir and HCQ for COVID-19 had no significant effect on pulmonary function test results 3 months following discharge for COVID-19, compared to SoC5,6. However, participants treated with remdesivir reported a higher CAT-score 3 months after hospital admission for COVID-19 than those receiving SoC.

We have previously reported decreased pulmonary function and persistent CT pathology in COVID-19 patients admitted to hospital, with results comparable to those presented in the current study2. However, by extending the analyses, we show a relationship between disease severity during hospital admission, as indicated by both low P/F-ratio and ICU admission, and impaired pulmonary function 3 months after discharge, as shown by reduced DLCO% in the fully adjusted model. The positive correlation between P/F-ratio during hospitalisation and DLCO% at the follow-up indicates that a reduced diffusion capacity of the lungs may persist at least 3 months after hospital admission.

Currently, it is not known how viral load affects the long-term pulmonary outcome following COVID-19. Another study found no difference in viral load between mild and severe disease in the acute phase29. Here we present a significant negative association between viral load in the oropharynx at hospital admission and DLCO% 3 months following hospitalisation. We also found a positive correlation between the presence of IgG antibodies to SARS-CoV-2 during the acute phase, and preserved diffusion capacity. Taken together, a low viral load combined with a robust and early antibody response against SARS-CoV-2 seem to protect against long-term impairment of pulmonary function. Finally, although COVID-19 seems to be characterized by sustained inflammation, plasma levels of inflammatory markers such as CRP and ferritin during hospital admission did not seem to predict a poorer pulmonary outcome at 3 months26.

In the present study, we also showed that nearly 40% of the participants had a CAT-score ≥ 10, which indicates a medium to high symptom burden after 3 months. The most prevalent symptoms were dyspnoea on exertion and loss of energy. Smokers are known to have a more severe course of COVID-19, and we found a higher symptom burden among participants with a history of smoking30. Interestingly, we found no association between CAT-score, and the pulmonary outcomes DLCO%, and chest CT-findings in these participants.

As an explorative outcome, we investigated the association between known pro-fibrotic mediators and the set of pulmonary outcomes. MMP-9 and VEGF-A were significantly elevated during the acute phase in these COVID-19 patients compared with healthy control subjects. From a hypothesis-generating point of view, it was interesting to see that MMP-9 levels, both during the acute phase and after 3-months, were significantly associated with persistent CT changes. Investigations have found MMP-9 in several cell types of the lung, such as alveolar macrophages, airway epithelium, and in neutrophils recruited to the lung during inflammatory processes31,32, but the major cellular sources of circulating MMP-9 in COVID-19 patients are still not evident. MMP-9 is involved in the complex process of extracellular matrix remodelling within the lung, and its activity is regulated by different TIMPs, including TIMP-233. Increased levels of MMP-9, combined with our finding of TIMP-2 levels within normal limits throughout the observation period, suggest enhanced MMP-9 activity in these COVID-19 patients. Based on the role of MMP-9 in the inflammatory process of acute lung injury, as well as its role in the development of pulmonary fibrosis34, MMP-9 has been suggested as a therapeutic target in the treatment of interstitial and acute inflammatory lung diseases35,36. Our findings may indicate that a similar approach could be of interest in COVID-19-induced lung injury, as well.

Our cohort had significantly elevated plasma levels of VEGF-A compared to healthy control subjects, both during hospital admission and 3 months after hospitalisation. Alveolar epithelial cells are claimed to be a major source of VEGF-A in humans, and VEGF-A has for long been a mediator of interest in understanding the pathophysiology in acute lung injury and ARDS, as well as in interstitial lung diseases, like idiopathic pulmonary fibrosis (IPF)37. This is reflected by the fact that an existing anti-fibrotic agent for treatment of IPF is targeting VEGF-A38. Even though we found no association between the increased VEGF-A levels and persistent findings on chest CT, it was interesting to observe that the level of VEGF-A was still elevated 3 months after hospitalisation. Other studies have found similar results, with significantly elevated VEGF-A levels in plasma of moderate to severe COVID-19 patients39,40. These findings might indicate that VEGF-A plays a role in the development of COVID-19 sequelae.

Our findings are limited by a small sample size, both regarding the evaluation of treatment effects and the interpretation of sub-analyses. Although we did not find significant effects of remdesivir and HCQ, except for a negative effect of remdesivir on CAT score, we cannot firmly exclude any treatment effects. Moreover, the participants’ pulmonary status prior to inclusion was not known, and established chronic lung disease prior to COVID-19 might have influenced the results. Due to a lack of additional plasma aliquots, we were not able to pursue our findings on matrix regulating mediators by analysing MMP-9 activity, other MMPs, and additional TIMPs. Moreover, the explorative part of this study presents data of descriptive nature and we lack data on the cellular sources of MMP-9 in these COVID-19 patients. Finally, associations do not necessarily reflect a causal relationship.

Conclusion

In this sub study of the NOR-SOLIDARITY trial, a high prevalence of respiratory symptoms, diffusion capacity impairment, and persistent chest CT-findings were observed 3 months after hospital admission for COVID-19. CAT-score seemed to be negatively influenced by treatment with remdesivir, while no difference in pulmonary function tests and CT-findings between remdesivir and HCQ, and their respective SoC group was found. A high viral load, low P/F-ratio, and admission to ICU seemed to predict reduced diffusion capacity in the lungs at 3 months, whereas high antibody levels at baseline had a positive effect on pulmonary function. The association between the pro-fibrotic mediator MMP-9 and persistent CT-findings highlights the need for more knowledge about the development of pulmonary fibrosis following hospitalisation for COVID-19.

Data availability

De-identified patient level data for participants in this study, statistical analysis plan and statistical coding can be made available after the approval of the institutional review board. Requests should be made to the corresponding author.

References

  1. Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273. https://doi.org/10.1038/s41586-020-2012-7 (2020).

    CAS  Article  PubMed  PubMed Central  ADS  Google Scholar 

  2. Lerum, T. V. et al. Dyspnoea, lung function and CT findings 3 months after hospital admission for COVID-19. Eur. Respir. J. https://doi.org/10.1183/13993003.03448-2020 (2020).

    Article  Google Scholar 

  3. Guler, S. A. et al. Pulmonary function and radiological features four months after COVID-19: First results from the national prospective observational Swiss COVID-19 lung study. Eur. Respir. J. 57, 2003690. https://doi.org/10.1183/13993003.03690-2020 (2021).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. Salvi, R. & Patankar, P. Emerging pharmacotherapies for COVID-19. Biomed. Pharmacother. 128, 110267. https://doi.org/10.1016/j.biopha.2020.110267 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. WHO Solidarity Trial Consortium. Repurposed antiviral drugs for covid-19—Interim WHO solidarity trial results. N. Engl. J. Med. https://doi.org/10.1056/NEJMoa2023184 (2020).

    Article  PubMed Central  Google Scholar 

  6. Barratt-Due, A. et al. Evaluation of the effects of remdesivir and hydroxychloroquine on viral clearance in COVID-19: A randomized trial. Ann. Intern. Med. https://doi.org/10.7326/M21-0653 (2021).

    Article  PubMed  PubMed Central  Google Scholar 

  7. Venkataraman, T. & Frieman, M. B. The role of epidermal growth factor receptor (EGFR) signaling in SARS coronavirus-induced pulmonary fibrosis. Antiviral Res. 143, 142–150. https://doi.org/10.1016/j.antiviral.2017.03.022 (2017).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. Davey, A., McAuley, D. F. & O’Kane, C. M. Matrix metalloproteinases in acute lung injury: Mediators of injury and drivers of repair. Eur. Respir. J. 38, 959–970. https://doi.org/10.1183/09031936.00032111 (2011).

    CAS  Article  PubMed  Google Scholar 

  9. Ikeda, K. et al. Serum surfactant protein D predicts the outcome of patients with idiopathic pulmonary fibrosis treated with pirfenidone. Respir. Med. 131, 184–191. https://doi.org/10.1016/j.rmed.2017.08.021 (2017).

    Article  PubMed  Google Scholar 

  10. Ando, M. et al. Significance of serum vascular endothelial growth factor level in patients with idiopathic pulmonary fibrosis. Lung 188, 247–252. https://doi.org/10.1007/s00408-009-9223-x (2010).

    CAS  Article  PubMed  Google Scholar 

  11. Ueland, T. et al. Distinct and early increase in circulating MMP-9 in COVID-19 patients with respiratory failure. J. Infect. 81, e41–e43. https://doi.org/10.1016/j.jinf.2020.06.061 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. Kerget, B. et al. Are serum interleukin 6 and surfactant protein D levels associated with the clinical course of COVID-19? Lung 198, 777–784. https://doi.org/10.1007/s00408-020-00393-8 (2020).

    CAS  Article  PubMed  Google Scholar 

  13. Brown, S. M. et al. Nonlinear imputation of PaO2/FIO2 from SpO2/FIO2 among mechanically ventilated patients in the ICU: A prospective. Observational study. Crit. Care Med. 45, 1317–1324. https://doi.org/10.1097/CCM.0000000000002514 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  14. Jones, P. W. et al. Development and first validation of the COPD assessment test. Eur. Respir. J. 34, 648–654. https://doi.org/10.1183/09031936.00102509 (2009).

    CAS  Article  PubMed  Google Scholar 

  15. Ghobadi, H., Ahari, S. S., Kameli, A. & Lari, S. M. The relationship between COPD assessment test (CAT) scores and severity of airflow obstruction in stable COPD patients. Tanaffos 11, 22–26 (2012).

    PubMed  PubMed Central  Google Scholar 

  16. Graham, B. L. et al. 2017 ERS/ATS standards for single-breath carbon monoxide uptake in the lung. Eur. Respir. J. https://doi.org/10.1183/13993003.00016-2016 (2017).

    Article  PubMed  Google Scholar 

  17. Graham, B. L. et al. Standardization of spirometry 2019 update. An Official American Thoracic Society and European Respiratory Society Technical Statement. Am. J. Respir. Crit. Care Med. 200, e70–e88. https://doi.org/10.1164/rccm.201908-1590ST (2019).

    Article  PubMed  PubMed Central  Google Scholar 

  18. Quanjer, P. H. et al. Multi-ethnic reference values for spirometry for the 3–95-yr age range: The global lung function 2012 equations. Eur. Respir. J. 40, 1324–1343. https://doi.org/10.1183/09031936.00080312 (2012).

    Article  PubMed  PubMed Central  Google Scholar 

  19. Stanojevic, S. et al. Official ERS technical standards: Global lung function initiative reference values for the carbon monoxide transfer factor for Caucasians. Eur. Respir. J. https://doi.org/10.1183/13993003.00010-2017 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  20. Hansell, D. M. et al. Fleischner society: Glossary of terms for thoracic imaging. Radiology 246, 697–722. https://doi.org/10.1148/radiol.2462070712 (2008).

    Article  PubMed  Google Scholar 

  21. Soyseth, V. et al. Diagnosis of biopsy verified usual interstitial pneumonia by computed tomography. Respir. Med. 109, 897–903. https://doi.org/10.1016/j.rmed.2015.05.008 (2015).

    Article  PubMed  Google Scholar 

  22. Richardson, S. et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City area. JAMA 323, 2052–2059. https://doi.org/10.1001/jama.2020.6775 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. Zhao, J. et al. Antibody responses to SARS-CoV-2 in patients with novel coronavirus disease 2019. Clin. Infect. Dis. 71, 2027–2034. https://doi.org/10.1093/cid/ciaa344 (2020).

    CAS  Article  PubMed  Google Scholar 

  24. Zhou, F. et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: A retrospective cohort study. Lancet 395, 1054–1062. https://doi.org/10.1016/S0140-6736(20)30566-3 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. Cummings, M. J. et al. Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: A prospective cohort study. Lancet 395, 1763–1770. https://doi.org/10.1016/S0140-6736(20)31189-2 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. Tay, M. Z., Poh, C. M., Renia, L., MacAry, P. A. & Ng, L. F. P. The trinity of COVID-19: Immunity, inflammation and intervention. Nat. Rev. Immunol. 20, 363–374. https://doi.org/10.1038/s41577-020-0311-8 (2020).

    CAS  Article  PubMed  Google Scholar 

  27. Tan, L. et al. Lymphopenia predicts disease severity of COVID-19: A descriptive and predictive study. Signal Transduct. Target Ther. 5, 33. https://doi.org/10.1038/s41392-020-0148-4 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Kasine, T. et al. Neutrophil count predicts clinical outcome in hospitalized covid-19 patients: Results from the NOR-Solidarity trial. J. Intern. Med. https://doi.org/10.1111/joim.13377 (2021).

    Article  PubMed  PubMed Central  Google Scholar 

  29. Yilmaz, A. et al. Upper respiratory tract levels of severe acute respiratory syndrome coronavirus 2 RNA and duration of viral RNA shedding do not differ between patients with mild and severe/critical coronavirus disease 2019. J. Infect. Dis. 223, 15–18. https://doi.org/10.1093/infdis/jiaa632 (2021).

    CAS  Article  PubMed  Google Scholar 

  30. Reddy, R. K. et al. The effect of smoking on COVID-19 severity: A systematic review and meta-analysis. J. Med. Virol. 93, 1045–1056. https://doi.org/10.1002/jmv.26389 (2021).

    CAS  Article  PubMed  Google Scholar 

  31. Vermeer, P. D. et al. MMP9 modulates tight junction integrity and cell viability in human airway epithelia. Am. J. Physiol. Lung Cell Mol. Physiol. 296, 751–762. https://doi.org/10.1152/ajplung.90578.2008 (2009).

    CAS  Article  Google Scholar 

  32. Wong, S., Birrell, M. A. & Belvisi, M. G. Matrix Metalloproteinases in Tissue Remodelling and Inflammation 1–20 (Springer, 2008).

    Book  Google Scholar 

  33. Cabral-Pacheco, G. A. et al. The roles of matrix metalloproteinases and their inhibitors in human diseases. Int. J. Mol. Sci. https://doi.org/10.3390/ijms21249739 (2020).

    Article  PubMed  PubMed Central  Google Scholar 

  34. Elkington, P. T. & Friedland, J. S. Matrix metalloproteinases in destructive pulmonary pathology. Thorax 61, 259–266. https://doi.org/10.1136/thx.2005.051979 (2006).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. Espindola, M. S. et al. Differential responses to targeting matrix metalloproteinase 9 in idiopathic pulmonary fibrosis. Am. J. Respir. Crit. Care Med. 203, 458–470. https://doi.org/10.1164/rccm.201910-1977OC (2021).

    CAS  Article  PubMed  Google Scholar 

  36. Solun, B. & Shoenfeld, Y. Inhibition of metalloproteinases in therapy for severe lung injury due to COVID-19. Med. Drug Discov. 7, 100052. https://doi.org/10.1016/j.medidd.2020.100052 (2020).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. Medford, A. R. & Millar, A. B. Vascular endothelial growth factor (VEGF) in acute lung injury (ALI) and acute respiratory distress syndrome (ARDS): Paradox or paradigm? Thorax 61, 621–626. https://doi.org/10.1136/thx.2005.040204 (2006).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. Grimminger, F., Gunther, A. & Vancheri, C. The role of tyrosine kinases in the pathogenesis of idiopathic pulmonary fibrosis. Eur. Respir. J. 45, 1426–1433. https://doi.org/10.1183/09031936.00149614 (2015).

    CAS  Article  PubMed  Google Scholar 

  39. Syed, F. et al. Excessive matrix metalloproteinase-1 and hyperactivation of endothelial cells occurred in COVID-19 patients and were associated with the severity of COVID-19. J. Infect. Dis. 224, 60–69. https://doi.org/10.1093/infdis/jiab167 (2021).

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. Kristensen, M. K., Plovsing, R. R., Berg, R. M. G., Krogh-Madsen, R. & Ronit, A. Cell adhesion molecules and vascular endothelial growth factor at the systemic and alveolar level in coronavirus disease 2019 acute respiratory distress syndrome. J. Infect. Dis. 224, 1101–1103. https://doi.org/10.1093/infdis/jiab347 (2021).

    CAS  Article  PubMed  Google Scholar 

Download references

Funding

Gunnar Einvik reports Grants from Boehringer Ingelheim in relation to previous work, unrelated to the submitted report. The National Program of Clinical Therapy Research in the Specialist Health Services, Norway (Grant no. 2020201), contributed with funding to the NOR-SOLIDARITY trial. The program had no role in the collection, analysis, or interpretation of the data in this sub-study.

Author information

Affiliations

Authors

Consortia

Contributions

P.A., A.B.D., A.M.D.R., I.C.O., T.K., K.N.H., and M.T. were responsible for the management, coordination, research activity planning, and execution of the NOR-SOLIDARITY trial. The 3-month follow-up protocol for pulmonary function and chest CT imaging was designed by T.M.A., P.A., A.B.D., G.E., T.V.L., and O.H.S. in collaboration. Patient data collection was performed by T.V.L., B.S., B.B., B.R.K., M.H., A.J., R.H., S.A., A.B.K., R.E., and the NOR-SOLIDARITY Consortium. T.M.A., J.R.R. and C.M. interpreted and analysed the CT images. T.B.D. was responsible for the biobank. F.L.J. was responsible for the antibody analyses, and F.M. for the viral analyses. Biomarkers were analysed by T.U. Statistical analyses were performed by N.N.M., I.C.O., and T.V.L. The first draft of the article was written by T.V.L., P.A., A.B.D., and O.H.S. All authors contributed with critical review of the manuscript and approval of the final version.

Corresponding author

Correspondence to Tøri Vigeland Lerum.

Ethics declarations

Competing interests

Gunnar Einvik reports grants from Boehringer Ingelheim in relation to previous work, unrelated to the submitted report. The other authors declare no competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

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/.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lerum, T.V., Maltzahn, N.N., Aukrust, P. et al. Persistent pulmonary pathology after COVID-19 is associated with high viral load, weak antibody response, and high levels of matrix metalloproteinase-9. Sci Rep 11, 23205 (2021). https://doi.org/10.1038/s41598-021-02547-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41598-021-02547-x

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing