Skip to main content

COVID-19 induces more pronounced extracellular matrix deposition than other causes of ARDS

Abstract

Background

Lung fibrosis is a major concern in severe COVID-19 patients undergoing mechanical ventilation (MV). Lung fibrosis frequency in post-COVID syndrome is highly variable and even if the risk is proportionally small, many patients could be affected. However, there is still no data on lung extracellular matrix (ECM) composition in severe COVID-19 and whether it is different from other aetiologies of ARDS.

Methods

We have quantified different ECM elements and TGF-β expression in lung tissue of 28 fatal COVID-19 cases and compared to 27 patients that died of other causes of ARDS, divided according to MV duration (up to six days or seven days or more). In COVID-19 cases, ECM elements were correlated with lung transcriptomics and cytokines profile.

Results

We observed that COVID-19 cases presented significant increased deposition of collagen, fibronectin, versican, and TGF-β, and decreased decorin density when compared to non-COVID-19 cases of similar MV duration. TGF-β was precociously increased in COVID-19 patients with MV duration up to six days. Lung collagen was higher in women with COVID-19, with a transition of upregulated genes related to fibrillogenesis to collagen production and ECM disassembly along the MV course.

Conclusions

Fatal COVID-19 is associated with an early TGF-β expression lung environment after the MV onset, followed by a disordered ECM assembly. This uncontrolled process resulted in a prominent collagen deposition when compared to other causes of ARDS. Our data provides pathological substrates to better understand the high prevalence of pulmonary abnormalities in patients surviving COVID-19.

Introduction

The main cause of death in severe cases of COVID-19 is respiratory failure due to diffuse alveolar damage (DAD), the histological hallmark of acute respiratory distress syndrome (ARDS) [1].

The extension of pulmonary fibroproliferation in patients with early ARDS of different causes predicts increased mortality, ventilator dependency, and susceptibility to multiple organ failure [2]. The mortality rate is higher in patients with increased levels of TGF-β and procollagen type III in lung lavage fluids [3]. The pathological fibroproliferation associated with the lack of reabsorption of the provisional extracellular matrix (ECM), decreased fibrin clearance, matrix metalloproteinases (MMP), and tissue inhibitor of metalloproteinase (TIMP) imbalance, lead to residual fibrosis with pulmonary dysfunction in ARDS survivors [4]. Patients who die of ARDS show pulmonary fibrosis, even when they die in relatively early stages [3].

In severe COVID-19, the development of lung fibrosis is also a major concern. Wendish et al. [5] showed that SARS-CoV-2 triggered profibrotic macrophage responses and pronounced fibroproliferative ARDS. The profibrotic cytokine TGF-β seems to be involved in the early development of the disease, regulating the immune response, antibody production, and tissue remodelling [6]. An autopsy study conducted by Li et al. [7], found that 43% of fatal cases of COVID-19 showed fibrosing DAD.

Lung fibrosis in survivors of COVID-19 has a highly variable prevalence. A study from the United Kingdom showed that 7% of the patients presented lung fibrosis four months after hospital discharge [8], whereas a Chinese study showed that 60% of the patients presented lung fibrosis 3–5 months after the disease onset [1]. Given that millions of people were infected (and re-infected) with SARS-CoV-2, even if the long-term risk is proportionally small, there still could be a significant number of patients developing lung fibrosis.

The ECM has a broad range of functions, from the cellular scaffold to cell adhesion, migration, and survival being involved in several important biological processes. In the lungs, the ECM composition is especially important to its homeostasis, injury-repair response, and mechanical properties [9]. Collagen, elastic fibres, fibronectin and proteoglycans are major components of the ECM in the lungs, having different roles in maintaining structural and cell homeostasis [10]. In our previous studies, we have observed an increased fibroblast proliferation [11] and an early predominance of intermediate DAD, with loose aggregates of fibroblasts admixed with scattered inflammatory cells, and collagen deposition shortly after the onset of MV [12] in the lungs of patients with fatal COVID-19. However, there is still no data on ECM composition of lung remodelling in severe COVID-19 and whether it is different from other aetiologies of ARDS.

Therefore, in this study, we have analysed different ECM elements and TGF-β expression in lung tissue of fatal COVID-19 cases compared to other causes of ARDS. Furthermore, we analysed the transcriptomics and cytokines profile, in COVID-19 cases.

Materials and methods

Ethics

The National Research Ethics Commission (CONEP) approved this study (CAAE #30364720.0.0000.0068). Autopsies were performed after written informed consent by the next-of-kin.

Study Population

We included all the autopsies of adult individuals performed between March and May of 2020 due to COVID-19 in Brazil (n = 28) when no vaccines were available. All COVID-19 patients had confirmation of SARS-CoV-2 infection by a positive RT-PCR result on the naso/oropharyngeal swab and/or lung tissue. During this period, there were 1898 admissions due to COVID-19 in the hospitals linked to the Medical School, with 659 deaths.

From May 2020 until July 2021, all autopsies were suspended in Brazil, and we had the approval to perform minimally invasive autopsy (MIA) for research purposes. MIA/US protocol in COVID-19 was described by Duarte-Neto et al. [13]. Lung tissues were collected during the MIA and the post-mortem interval ranged from 3 to 23 h and 51 min (mean of 22 h and 23 min). Lung tissue was sampled from eight regions, in a combination of the upper and lower chest (lateral and medial, 4 sites) and right and left lungs (2 sides). In each sampling site, six samples were collected. Immediately after sampling, pulmonary samples either were immersed in 10% formalin solution for 24 h or immediately snap frozen at -80 °C. For the immunohistochemistry analysis, two representative formalin-fixed paraffin-embedded biopsies (3–6 samples per slide) presenting DAD pattern and no histological signs of secondary bacterial/fungus infection were selected.

As controls, we selected retrospectively post-mortem lung tissue of 27 cases with criteria for non-COVID-19 ARDS diagnosis. Clinical diagnosis of ARDS was defined according to the Berlin definition [14] in addition to histological findings of diffuse alveolar damage (DAD) and the absence of chronic lung diseases. The cause of non-COVID-19 DAD in these patients was, in most of them, pulmonary or extrapulmonary sepsis. Lung tissues were collected during the conventional autopsy between 2004 and 2010. The post-mortem interval ranged from 4 h to 50 min to 21 h and 50 min (mean of 11 h and 18 min).

All patients were submitted to MV and both groups were further divided by MV duration: up to 6 days (L7) and 7 days or more (H7).

xMAP Luminex Cytokine Multiplex Assay of COVID-19 cases

Samples of frozen lung tissue from one post-mortem biopsy site per case of the COVID-19 patients were prepared using Bio-plex Cell Lysis Kit (Bio-rad Laboratories, Hercules, CA, USA), following the manufacturer’s instructions. There was no availability of frozen tissue from the non-COVID-19 cases. Cytokines and chemokines were assessed using the following kits: Bio-plex TGF-β, Bio-plex Human Chemokine (Bio-rad, Hercules, California), and Milliplex MAP Human Cytokine/Chemokine (Millipore Corp., Billerica, MA) according to the kit-specific protocols. Measurements were made using a Magpix analytical test instrument by Luminex (Austin, TX). XPONENT 4.2 software and Milliplex Analyst 5.1 software (Millipore Corp., Billerica, MA, USA) were used for data analysis. To standardize the results, the total proteins of the lung samples were quantified with the TP2 kit (Roche Diagnostics GmbH, Mannheim, Germany) using COBAS C11 equipment (Roche Instrument Center, Rotkreuz, Switzerland). Data are expressed in pg/g and were correlated with clinical and histological parameters.

Extracellular matrix analysis

Lung tissue slides were stained with Sirius Red (collagen), and Weigert’s resorcin-fuchsin with oxidation (elastic fibres). Lung tissue was immunostained using anti-Fibronectin, anti-Versican, anti-Decorin, and anti-TGF-β antibodies. Table S1 (Additional file 1) shows the commercial sources and dilutions of the antibodies. The ECM elements were quantified in the lung parenchyma, excluding airways, large blood vessels (only capillaries were not excluded), and pleura. The proportional area, positive stained area per lung parenchyma area (µm²/µm²) were measured in 15 high-power fields using the Image-Pro Plus 4.1 software (Media Cybernetics, Silver Spring, MD, USA). We selected a range of colors encompassing the positive staining for each marker, avoiding any background in the slide. There was no statistical difference in the analysed lung tissue area of each marker between the COVID-19 and non-COVID-19 groups (data not shown).

Transcriptomics

For differential expression analysis, twelve COVID-19 cases were divided into 2 groups according to the MV duration: L7 (n = 5) and H7 (n = 7). We included control cases (n = 4) of individuals that died due to non-pulmonary causes. Characteristics of the control individuals are presented in Table S2 (Additional file 1). Data have been deposited in NCBI’s Gene Expression Omnibus (GSE205099).

Details about RNA extraction and sequencing are presented in supplementary material (Additional file 1). Data, read filtering and gene expression estimation were done as described in Erjefält et al. [11]. A preliminary analysis detected 23 genes related to heart tissue in two samples (Additional file 2 - Table S3) that were removed from further analysis. Heart tissue might contaminate lung biopsy during MIAS puncture trajectory. Due to the small amount of tissue recovery, the frozen sections were submitted to RNA extraction, without histological analysis. DESeq2 [15] was used to identify the differentially expressed genes (DEGs) between (1) L7 vs. control; (2) H7 vs. control; and (3) L7 vs. H7. For further details on the DEGs analysis, see the supplementary material (Additional file 1).

Real-time PCR

Total RNA was extracted using Trizol (Thermo Fisher Scientific, Waltham, MA, USA) followed by a DNAse treatment (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s protocol. Reverse-transcribed cDNA samples were synthesized using High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific, Waltham, MA, USA) also according to the manufacturer’s instructions. Real-time PCR reactions were conducted using QuantStudio 12 Instrument K Flex Real-Time PCR System (Applied Biosystems, Waltham, MA, USA) in the presence of Fast SYBR Green Master Mix (Thermo Fisher Scientific, Vilnius, Lithuania) and specific primers for the Collagen type I alpha 1 chain (COL1A1). The relative quantification was calculated after normalization against the levels of the reference gene ribosomal protein S18 (RPS18). Primer sequences are described in Table S4 (Additional file 1).

Statistical analysis

SPSS 21 software (SPSS Inc/IBM Chicago, USA) was used for the statistical analyses. Data distribution was assessed by Shapiro-Wilk normality test. Depending on the data distribution, T-student or Mann-Whitney test were used to compare two groups and Kruskal–Wallis or one-way analysis of variance (ANOVA) test, followed by the Bonferroni or Tukey posthoc test to compare four groups. Statistical difference was assumed at the 5% significance level.

Factorial analysis of covariance (ANCOVA) was performed to explore influence factors on each dependent variable. We included COVID-19 and sex as factors and age, MV duration, BMI, and PaO2/FiO2 ratio as covariates. The covariates were added as main effects only and no interaction terms involving covariates were included. In addition, we calculated the partial eta squared (η2) as measures of effect size. Furthermore, we performed the Spearman correlation test between variables; coefficients (r) were considered statistically significant at p < 0.05.

Results

Demographics and clinical features

Lung samples from 28 COVID-19 patients (11 F/17 M), mean age of 56.9 years (range 32–88), were selected for this study. Part of this population has been previously presented [11, 12]. Table 1 and S5 (Additional file 1) show COVID-19 patients’ demographics, comorbidities, and clinical characteristics. There were 10 patients in the L7 group and 18 in the H7 group (demographics are shown in Table S7 - Additional file 1).

Table 1 Demographic and clinical characteristics of COVID-19 and non-COVID-19 patients

The 27 non-COVID-19 DAD patients (15 F/12 M), mean age of 44.6 years (range 18–81), had non-pulmonary (n = 14) or pulmonary ARDS (n = 13). The most common underlying conditions for respiratory failure were bronchopneumonia (n = 14), followed by liver diseases (n = 11), extrapulmonary sepsis (n = 9), gastrointestinal diseases (n = 7) and bleeding (n = 5), cardiovascular diseases (n = 5), neurologic diseases (n = 5), and renal diseases (n = 4). Table 1 and S6 (Additional file 1) show non-COVID-19 patients’ demographics, comorbidities, and clinical characteristics. There were 17 patients in the L7 group and 10 in the H7 group (demographics are shown in Table S7 - Additional file 1).

There was no statistical difference in the MV and hospitalization duration between the COVID-19 and non-COVID-19 groups.

ECM composition: COVID-19 DAD x non-COVID-19 DAD

The collagen density was increased in the COVID-19 cases compared to the non-COVID-19 (p < 0.0001) (Fig. 1). There was no statistically significant difference in elastic fibre density between the groups (Fig. 1). The fibronectin (p < 0.0001) and the versican (p = 0.044) densities were increased in the COVID-19 cases compared to the non-COVID-19, whereas the decorin density was decreased in the COVID-19 cases compared to the non-COVID-19 (p = 0.033) (Figs. 1 and 2). In addition, the expression of TGF-β was higher in the COVID-19 than in the non-COVID-19 cases (p = 0.003) (Fig. 2).

Fig. 1
figure 1

Graphical representation of the lung tissue extracellular matrix composition: Comparison between the COVID-19 and non-COVID-19 cases. A: Collagen (µm²/ µm²). Photomicrographs of non-COVID-19 (B) and COVID-19 (C) stained with Sirius red. D: Elastic Fibres (µm²/ µm²). Photomicrographs of non-COVID-19 (E) and COVID-19 (F) stained with Weigert’s resorcin-fuchsin with oxidation. G: Fibronectin (µm²/ µm²). Photomicrographs of non-COVID-19 (H) and COVID-19 (I) immunostainned with anti-fibronectin antibody. Bars show the mean and each point represents one patient. Scale bar = 100 μm

Fig. 2
figure 2

Graphical representation of the lung tissue extracellular matrix composition: Comparison between the COVID-19 and non-COVID-19 cases. A: Decorin (µm²/ µm²). Photomicrographs of non-COVID-19 (B) and COVID-19 (C) immunostainned with anti-decorin antibody. D: Versican (µm²/ µm²). Photomicrographs of non-COVID-19 (E) and COVID-19 (F) immunostainned with anti-versican antibody. G: TGF-beta (µm²/ µm²). Photomicrographs of non-COVID-19 (H) and COVID-19 (I) immunostainned with anti-TGF-beta antibody. Bars show the mean and each point represents one patient. Scale bar = 100 μm

The influence of MV duration, age, BMI, and PaO 2/FiO2 ratio on the ECM composition and TGF-β

Table 2 shows the results of the ANCOVA analysis.

Table 2 Factorial ANCOVA analysis

All assessed ECM components, collagen (p < 0.0001; η2 = 0.38), fibronectin (p < 0.0001; η2 = 0.34), versican (p = 0.013; η2 = 0.15), decorin (p = 0.028; η2 = 0.12), and TGF-β (p = 0.018; η2 = 0.14), had their expression influenced by COVID-19, except the elastic fibres (p = 0.271; η2 = 0.03).

The only component that was influenced by sex was collagen (p = 0.004; η2 = 0.19). The COVID-19 patients of the female sex presented higher density of collagen in the lung tissue (p = 0.004) (Fig. 3A). Between the COVID-19 patients of the female and male sex, we did not observe a difference in age, MV duration, BMI, or PaO2/FiO2 ratio. There was no difference in the collagen between males and females in the non-COVID-19 group. The characteristics of males and females are presented in Table S8 (Additional file 1).

Fig. 3
figure 3

Graphical representation of the lung tissue extracellular matrix composition. A: Collagen density – COVID-19 Males x COVID-19 Females (µm²/ µm²). B: TGF-β expression divided by the duration of the mechanical ventilation (µm²/ µm²). C: Elastic fibre density divided by the duration of the mechanical ventilation (µm²/ µm²). D: Decorin density – Adult COVID-19 x Elderly COVID-19 (> 60 years of age) (µm²/ µm²). L7: up to 6 days of mechanical ventilation; H7: more than 7 days of mechanical ventilation. *p < 0.05 compared to the COVID-19 L7 group. Bars show the mean and each point represents one patient

The components influenced by the MV duration were the TGF-β (p = 0.044, η2 = 0.11) and elastic fibres (p = 0.011, η2 = 0.16). When the groups are divided by the MV duration, the COVID-19 L7 group had higher TGF-β density than the non-COVID-19 L7 and the non-COVID-19 H7 groups (p = 0.003 and p = 0.001, respectively). In addition, the COVID-19 L7 has higher density of TGF-β than the COVID-19 H7 group (p = 0.043) (Fig. 3B). Considering all patients, regardless of COVID-19, there was a decrease in elastic fibres in the H7 cases when compared to L7 (p = 0.008) (Fig. 3C).

There was a positive correlation between age and decorin density (r = 0.411; p = 0.022). Elderly COVID-19 cases (> 60 years of age) had higher decorin than the adult COVID-19 cases (p = 0.033) (Fig. 3D).

Correlations between the ECM components, cytokines/chemokines, and clinical variables in COVID-19 cases

Data showing the mean concentration levels of the cytokines/chemokines are shown in Table S9 (Additional file 1). Figure 4 shows the correlation plot.

Fig. 4
figure 4

Correlation plot showing the correlations between histological analysis of the lung tissue extracellular matrix elements, TGF-beta protein expression, cytokines and clinical data of COVID-19 patients. The specific correlation coefficients and p-values of the statistically significant correlations are presented in the main text. BMI: body mass index. P/F ratio: PaO2/FiO2 ratio. MV: mechanical ventilation. Symptoms onset – Death: Time from symptom onset to death. Symptoms onset - Hospitalization: Time from symptom onset to hospitalization. TGF-beta: Transforming growth factor-beta. MMP: metalloproteinase. IFN: interferon. TNF: Tumor necrosis factor. IL: Interleukin. IP-10: Interferon-gamma inducible Protein 10 kDa. MCP: Monocyte chemoattractant protein. MIP: Macrophage inflammatory protein. CCL: Chemokine (C-C motif) ligand. MIG: Monokine induced by interferon-gamma. CXCL: chemokine (C-X-C motif) ligand. TARC: Thymus and activation-regulated chemokine

The collagen density was positively correlated with the elastic fibre density (r = 0.528; p = 0.005) and TIMP2 (r = 0.435; p = 0.034).

The elastic fibre density was positively correlated with the MMP2 (r = 0.483; p = 0.02) and MMP10 (r = 0.489; p = 0.011). Moreover, it was negatively correlated with the MV duration (r=-0.483; p = 0.011) and with the time elapsed between hospital admission and death (r=-0.439; p = 0.022).

Versican density was positively correlated with TGF-β3 (r = 0.423; p = 0.045), and the decorin density was positively correlated with the TGF-β1 levels (r = 0.428; p = 0.042) and age (r = 0.411; p = 0.022). Decorin was also negatively correlated with the interval between the onset of symptoms and hospital admission (r=-0.414; p = 0.032). TGF-β density was negatively correlated with IFN-α2 (r=-0.604, p = 0.022).

Increased expression of collagen-related genes and fibrillogenesis in severe COVID-19 patients

When compared to the control cases, the COVID-19 L7 cases showed 70 DEGs (54 up- and 16 down-regulated) and the COVID-19 H7 showed 208 DEGs (139 up- and 69 down-regulated). The COVID-19 L7 compared to the COVID-19 H7 showed 46 DEGs (34 up- and 12 down-regulated).

DEG analysis revealed increased expression of genes related to collagen assembly and organization in the L7 group compared to the control, such as tolloid-like 2, ADAMTS8, and collagen III. Group L7 had up-regulated genes related to complement activation, humoral immune response and B cell activation, phagocytosis, activation of the immune response, and mitosis (Table S10 - Additional file 3).

The DEG analysis in the H7 group compared to the control showed increased expression of genes related to fibrillogenesis and ECM degradation, such as biglycan, collagen types I, III, and V, CTHRC1, MMP-8, MMP-9, and cathepsin G. The up-regulated genes were related to phagocytosis, complement activation, immune response, B cell activation and immunoglobulin production, mitosis, neutrophil-mediated immunity, collagen synthesis, ECM disassembly, and regulation of homotypic cell-cell adhesion (Table S11 - Additional file 4).

When we compared COVID-19 L7 to the COVID-19 H7 group, the downregulated genes were related to response to interleukin-1, leukocyte migration and neutrophil chemotaxis, humoral immune response, cell adhesion molecule production, and negative regulation of proteolysis (Table S12 - Additional file 5).

Although the gene expression of COL1A1 assessed by real-time PCR was increased in the COVID-19 H7 (4.36 ± 2.9; n = 7) and compared to the control (1.01 ± 0.21; n = 3), we did not observe the statistical difference (p = 0.143), probably due to the low number of samples analysed.

Discussion

In this study, we have shown that COVID-19 cases present increased deposition of collagen, fibronectin, versican, and TGF-β in the lungs when compared to non-COVID-19 ARDS cases of similar ventilation duration. ECM deposition occurred early after MV onset in both groups, and TGF-β was precociously increased in COVID-19. There was a transition of upregulated genes related to fibrillogenesis to collagen production and ECM disassembly along the MV course, with lung collagen was higher in women with COVID-19. To our knowledge, this is the first study to address the ECM composition and its influence factors in COVID-19.

Our data reassure the central role of TGF-β in the remodelling taking place in severe COVID-19. TGF-β can be secreted by different cell populations of lung inflammatory and structural cells [6]. Xu et al. [16] showed that TGF-β expression was increased in alveolar epithelial cells after SARS-CoV-2 infection, which could block apoptosis of the infected host cells, promoting tissue fibrosis [17]. Vaz de Paula et al. [18] showed increased lung TGF-β expression in fatal COVID-19 when compared to H1N1 patients.

There are some mechanisms to help explain the early and excessive TGF-β increase in fatal COVID-19 cases, as observed in this study. In severe cases, SARS-CoV-2 induces poor interferon production [19], which could explain the increased TGF-β expression, since IFN-γ specifically inhibits an early step of TGF-β-induced activation of Smad3 [20]. In our data, we also observed a negative correlation between TGF-β density and IFN-α2. Accordingly, Hu et al. [21] found a negative correlation between the lung fibrotic volume and two circulating interferons, IFN-gamma and IFN-α2, at the discharge of COVID-19 patients.

We observed a decrease in the decorin density in COVID-19 lungs. Decorin regulates fibrillogenesis acting as an anti-fibrotic agent, partially via binding and neutralizing the TGF-β [22] and it may also bind to collagen fibrils aiding in the stabilization of collagen [23]. We observed a positive correlation between decorin density and TGF-β1 levels. TGF-β1 can downregulate the decorin synthesis in fibroblasts and, in turn, decorin can inhibit TGF-β1 [24]. We also observed an increased biglycan gene expression in COVID-19 cases with more than 7 days of VM. Unlike decorin, biglycan may be responsible for the sequestration of TGF-β to cell-surface receptors, promoting ECM formation [25].

We observed an important increase in fibronectin in the COVID-19 cases. Xu et al. [16] showed that the SARS-CoV-2 infection increases the fibronectin gene expression in alveolar epithelial cells. In situations of lung remodelling, fibronectin is one of the first molecules to be deposited in the extracellular space, and it regulates assembly and stability of several ECM components [26]. Furthermore, the binding of fibronectin with latent TGF-β binding proteins and later activation by proteases and integrins are major regulators of TGF-β bioavailability and signalling [26]. Fibronectin levels were increased in the plasma of COVID-19 patients compared to non-COVID-19 patients [27] and even higher in fatal COVID-19 [28].

The versican density was also increased in COVID-19. Versican plays an important role in immunity and inflammation influencing key events such as cell adhesion, proliferation, migration, and ECM remodelling [29]. Versican mRNA can be differentially spliced generating five isoforms (V0-V4), which differently affect cell behaviour [29, 30]. We observed a positive correlation between versican and TGF-β3. TFG-β3 induces the expression of the variants V0 and V1 in fibroblasts [31], which may perpetuate tissue remodelling by increasing proliferation and resistance to apoptosis of fibroblasts [32]. Moreover, we observed the increased expression of several genes linked to mitosis including Ki-67, as per our previous study [11], indicating that COVID-19 drives the lungs to a highly fibroproliferative state.

There was an increased deposition of collagen in COVID-19 lungs, with an early increased expression of genes linked to ECM orchestration and organization, such as Tolloid-like 2, ADAMTS8, and collagen III. After seven days of VM, there was increased gene expression of collagen I, III, V, and CTHRC-1. Interestingly, a lung cell atlas study identified a subpopulation of CTHRC-1 + fibroblasts, present exclusively within fibrotic areas expressing collagen I, III, and V [33]. Increased gene expression of collagens was accompanied by the upregulation of MMP-8, MMP-9, and cathepsin G. In addition, there was a positive correlation between collagen density and TIMP-2 and a higher collagen density in females compared to males. Although sex hormones are likely to have some influence over fibrillogenesis, their role is poorly understood. Although it is not unanimous, many studies already showed that the female sex is associated with an increased risk for long COVID-19 [34,35,36,37].

The increased ECM deposition and protease expression represent a dynamic process of deposition and degradation that may lead to destabilization. Furthermore, endogenous molecules, such as versican and biglycan, can interact with toll-like receptors, inducing the production of proinflammatory cytokines [30]. Interestingly, we observed a negative correlation between the elastic fibres and VM duration, and a positive correlation with MMP-2 and MMP-10. It has been described that the products resulting from elastic fibre hydrolysis can induce leukocyte activation, inhibit elastase activity, and stimulates fibroblast proliferation [38].

Taken together, our data suggest that in severe cases of COVID-19, an early expression of TGF-β environment and ECM-rich synthesis and degradation takes place, leading to a more pronounced lung remodelling than other causes of ARDS. Our data complement the work of Mothes et al. [39], which showed the role of a non-resolving local immune response that extends beyond active viral infection and perpetuates tissue remodelling.

Our study has some limitations, such as the low number of the studied cases, mainly in the transcriptomic analysis and PCR, due to the difficulties in obtaining COVID-19 autopsy tissue during the pandemics, including ones with good RNA quality. Another limitation was the lack of frozen tissue from the non-COVID-19 cases which limits our interpretation regarding the differences between COVID-19 and non-COVID-19. We do not have cases with very prolonged ventilation periods, as not infrequently observed in some COVID-19 cases along the pandemic course, as these cases were collected very early during the pandemic in Brazil. Since all patients had a fatal outcome, we do not know how long ECM abnormalities in the lung would persist. A study performed on COVID-19 survivors hospitalized at our Medical School in Sao Paulo during the same period showed that 30% had persistent dyspnoea and 32% had reduced forced vital capacity [37].

In summary, our data show that severe COVID-19 is associated with an early expression of TGF-β lung environment after MV onset. This particular environment led to a disordered ECM assembly with upregulation of genes related to ECM fibrillogenesis and disassembly. This uncontrolled process resulted in a prominent collagen deposition when compared to other causes of ARDS. Our data provides pathological substrates to better understand the high prevalence of pulmonary abnormalities in patients surviving COVID-19.

Data Availability

Data on transcriptome reads have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE205099 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE205099). Further data are available from the corresponding author upon reasonable request.

References

  1. Li X, Shen C, Wang L, Majumder S, Zhang D, Deen MJ, et al. Pulmonary fibrosis and its related factors in discharged patients with new coronavirus pneumonia: a cohort study. Respir Res. 2021;22(1):203. https://doi.org/10.1186/s12931-021-01798-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Ichikado K, Muranaka H, Gushima Y, Kotani T, Nader HM, Fujimoto K, et al. Fibroproliferative changes on high-resolution CT in the acute respiratory distress syndrome predict mortality and ventilator dependency: a prospective observational cohort study. BMJ Open. 2012;2(2):e000545. https://doi.org/10.1136/bmjopen-2011-000545.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Cabrera-Benitez NE, Laffey JG, Parotto M, Spieth PM, Villar J, Zhang H, et al. Mechanical ventilation-associated lung fibrosis in acute respiratory distress syndrome: a significant contributor to poor outcome. Anesthesiology. 2014;121(1):189–98. https://doi.org/10.1097/ALN.0000000000000264.

    Article  PubMed  Google Scholar 

  4. Burnham EL, Janssen WJ, Riches DW, Moss M, Downey GP. The fibroproliferative response in acute respiratory distress syndrome: mechanisms and clinical significance. Eur Respir J. 2014;43(1):276–85. https://doi.org/10.1183/09031936.00196412.

    Article  PubMed  Google Scholar 

  5. Wendisch D, Dietrich O, Mari T, von Stillfried S, Ibarra IL, Mittermaier M, et al. SARS-CoV-2 infection triggers profibrotic macrophage responses and lung fibrosis. Cell. 2021;184(26):6243–6261e27. https://doi.org/10.1016/j.cell.2021.11.033.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Ferreira-Gomes M, Kruglov A, Durek P, et al. SARS-CoV-2 in severe COVID-19 induces a TGF-β-dominated chronic immune response that does not target itself. Nat Commun. 2021;12(1):1961. https://doi.org/10.1038/s41467-021-22210-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Li Y, Wu J, Wang S, Heinrich F, Tizian C, Heinz GA, et al. Progression to fibrosing diffuse alveolar damage in a series of 30 minimally invasive autopsies with COVID-19 pneumonia in Wuhan, China. Histopathology. 2021;78(4):542–55. https://doi.org/10.1111/his.14249.

    Article  PubMed  Google Scholar 

  8. Robey RC, Kemp K, Hayton P, Mudawi D, Wang R, Greaves M, et al. Pulmonary sequelae at 4 months after COVID-19 infection: a single-centre experience of a COVID Follow-Up service. Adv Ther. 2021;38(8):4505–19. https://doi.org/10.1007/s12325-021-01833-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Zhou Y, Horowitz JC, Naba A, Ambalavanan N, Atabai K, Balestrini J, Bitterman PB, Corley RA, Ding BS, Engler AJ, Hansen KC, Hagood JS, Kheradmand F, Lin QS, Neptune E, Niklason L, Ortiz LA, Parks WC, Tschumperlin DJ, White ES, Chapman HA, Thannickal VJ. Extracellular matrix in lung development, homeostasis and disease. Matrix Biol. 2018;73:77–104. https://doi.org/10.1016/j.matbio.2018.03.005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Annoni R, Lanças T, Yukimatsu Tanigawa R, de Medeiros Matsushita M, de Morais Fernezlian S, Bruno A, Fernando Ferraz da Silva L, Roughley PJ, Battaglia S, Dolhnikoff M, Hiemstra PS, Sterk PJ, Rabe KF, Mauad T. Extracellular matrix composition in COPD. Eur Respir J. 2012;40(6):1362–73. https://doi.org/10.1183/09031936.00192611.

    Article  PubMed  Google Scholar 

  11. Erjefält JS, de Souza Xavier Costa N, Jönsson J, Cozzolino O, Dantas KC, Clausson CM et al. Diffuse alveolar damage patterns reflect the immunological and molecular heterogeneity in fatal COVID-19. EBioMedicine 2022; 83: 104229. https://doi.org/10.1016/j.ebiom.2022.104229.

  12. Mauad T, Duarte-Neto AN, da Silva LFF, de Oliveira EP, de Brito JM, do, Nascimento ECT et al. Tracking the time course of pathological patterns of lung injury in severe COVID-19. Respir Res 2021; 22 (1): 32. https://doi.org/10.1186/s12931-021-01628-9.

  13. Duarte-Neto AN, Monteiro RAA, da Silva LFF, Malheiros DMAC, de Oliveira EP, Theodoro-Filho J, et al. Pulmonary and systemic involvement in COVID-19 patients assessed with ultrasound-guided minimally invasive autopsy. Histopathology. 2020;77(2):186–97. https://doi.org/10.1111/his.14160.

    Article  PubMed  PubMed Central  Google Scholar 

  14. ARDS Definition Task Force, Ranieri VM, Rubenfeld GD, Thompson BT, Ferguson ND, Caldwell E, et al. Acute respiratory distress syndrome: the Berlin definition. JAMA. 2012;307(23):2526–33. https://doi.org/10.1001/jama.2012.5669.

    Article  CAS  Google Scholar 

  15. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Xu J, Xu X, Jiang L, Dua K, Hansbro PM, Liu G. SARS-CoV-2 induces transcriptional signatures in human lung epithelial cells that promote lung fibrosis. Respir Res. 2020;21(1):182. https://doi.org/10.1186/s12931-020-01445-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hemmat N, Asadzadeh Z, Ahangar NK, Alemohammad H, Najafzadeh B, Derakhshani A, et al. The roles of signaling pathways in SARS-CoV-2 infection; lessons learned from SARS-CoV and MERS-CoV. Arch Virol. 2021;166(3):675–96. https://doi.org/10.1007/s00705-021-04958-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Vaz de Paula CB, Nagashima S, Liberalesso V, Collete M, da Silva FPG, Oricil AGG, et al. COVID-19: immunohistochemical analysis of TGF-β signaling pathways in Pulmonary Fibrosis. Int J Mol Sci. 2021;23(1):168. https://doi.org/10.3390/ijms23010168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Znaidia M, Demeret C, van der Werf S, Komarova AV. Characterization of SARS-CoV-2 evasion: Interferon Pathway and Therapeutic Options. Viruses. 2022;14(6):1247. https://doi.org/10.3390/v14061247.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Ulloa L, Doody J, Massagué J. Inhibition of transforming growth factor-beta/SMAD signalling by the interferon-gamma/STAT pathway. Nature. 1999;397:710–3. https://doi.org/10.1038/17826.

    Article  CAS  PubMed  Google Scholar 

  21. Hu ZJ, Xu J, Yin JM, Li L, Hou W, Zhang LL, et al. Lower circulating Interferon-Gamma is a risk factor for lung fibrosis in COVID-19 patients. Front Immunol. 2020;11:585647. https://doi.org/10.3389/fimmu.2020.585647.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Goetsch KP, Niesler CU. The extracellular matrix regulates the effect of decorin and transforming growth factor β-2 (TGF-β2) on myoblast migration. Biochem Biophys Res Commun. 2016;479(2):351–7. https://doi.org/10.1016/j.bbrc.2016.09.079.

    Article  CAS  PubMed  Google Scholar 

  23. Keene DR, San Antonio JD, Mayne R, McQuillan DJ, Sarris G, Santoro SA, Iozzo RV. Decorin binds near the C terminus of type I collagen. J Biol Chem. 2000;275(29):21801–4. https://doi.org/10.1074/jbc.C000278200.

    Article  CAS  PubMed  Google Scholar 

  24. Yamaguchi Y, Mann DM, Ruoslahti E. Negative regulation of transforming growth factor-beta by the proteoglycan decorin. Nature. 1990;346(6281):281–4. https://doi.org/10.1038/346281a0.

    Article  CAS  PubMed  Google Scholar 

  25. Okada H. Tolloid-like proteinases orchestrate extracellular matrix formation. J Am Soc Nephrol. 2011;22(4):588–9. https://doi.org/10.1681/ASN.2011020131.

    Article  PubMed  Google Scholar 

  26. Todorovic V, Rifkin DB. LTBPs, more than just an escort service. J Cell Biochem. 2012;113(2):410–8. https://doi.org/10.1002/jcb.23385.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Overmyer KA, Shishkova E, Miller IJ, Balnis J, Bernstein MN, Peters-Clarke TM, et al. Large-scale multi-omic analysis of COVID-19 severity. Cell Syst. 2021;12(1):23–40e7. https://doi.org/10.1016/j.cels.2020.10.003.

    Article  CAS  PubMed  Google Scholar 

  28. Lemańska-Perek A, Krzyżanowska-Gołąb D, Dragan B, Tyszko M, Adamik B. Fibronectin as a marker of Disease Severity in critically ill COVID-19 patients. Cells. 2022;11(9):1566. https://doi.org/10.3390/cells11091566.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Wight TN, Kang I, Evanko SP, Harten IA, Chang MY, Pearce OMT, et al. Versican-A critical Extracellular Matrix Regulator of immunity and inflammation. Front Immunol. 2020;11:512. https://doi.org/10.3389/fimmu.2020.00512.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Andersson-Sjöland A, Hallgren O, Rolandsson S, Weitoft M, Tykesson E, Larsson-Callerfelt AK, et al. Versican in inflammation and tissue remodeling: the impact on lung disorders. Glycobiology. 2015;25(3):243–51. https://doi.org/10.1093/glycob/cwu120.

    Article  CAS  PubMed  Google Scholar 

  31. Norian JM, Malik M, Parker CY, Joseph D, Leppert PC, Segars JH, et al. Transforming growth factor beta 3 regulates the versican variants in the extracellular matrix-rich uterine leiomyomas. Reprod Sci. 2009;16(12):1153–64. https://doi.org/10.1177/1933719109343310.

    Article  CAS  PubMed  Google Scholar 

  32. Sheng W, Wang G, Wang Y, Liang J, Wen J, Zheng PS, et al. The roles of versican V1 and V2 isoforms in cell proliferation and apoptosis. Mol Biol Cell. 2005;16(3):1330–40. https://doi.org/10.1091/mbc.e04-04-0295.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Tsukui T, Sun KH, Wetter JB, Wilson-Kanamori JR, Hazelwood LA, Henderson NC, et al. Collagen-producing lung cell atlas identifies multiple subsets with distinct localization and relevance to fibrosis. Nat Commun. 2020;11(1):1920. https://doi.org/10.1038/s41467-020-15647-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Al-Husinat L, Nusir M, Al-Gharaibeh H, Alomari AA, Smadi MM, Battaglini D, et al. Post-COVID-19 syndrome symptoms after mild and moderate SARS-CoV-2 infection. Front Med (Lausanne). 2022;9:1017257. https://doi.org/10.3389/fmed.2022.1017257.

    Article  PubMed  Google Scholar 

  35. Bai F, Tomasoni D, Falcinella C, Barbanotti D, Castoldi R, Mulè G, et al. Female gender is associated with long COVID syndrome: a prospective cohort study. Clin Microbiol Infect. 2022;28(4):611e9–16. https://doi.org/10.1016/j.cmi.2021.11.002.

    Article  CAS  Google Scholar 

  36. Chen C, Haupert SR, Zimmermann L, Shi X, Fritsche LG, Mukherjee B, et al. Global prevalence of Post-Coronavirus Disease 2019 (COVID-19) Condition or Long COVID: a Meta-analysis and systematic review. J Infect Dis. 2022;226(9):1593–607. https://doi.org/10.1093/infdis/jiac136.

    Article  CAS  PubMed  Google Scholar 

  37. Ferreira JC, Moreira TCL, de Araújo AL, Imamura M, Damiano RF, Garcia ML, et al. Clinical, sociodemographic and environmental factors impact post-COVID-19 syndrome. J Glob Health. 2022;12:05029. https://doi.org/10.7189/jogh.12.05029.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Maquart FX, Pasco S, Ramont L, Hornebeck W, Monboisse JC. An introduction to matrikines: extracellular matrix-derived peptides which regulate cell activity. Implication in tumor invasion. Crit Rev Oncol Hematol. 2004;49(3):199–202. https://doi.org/10.1016/j.critrevonc.2003.06.007.

    Article  PubMed  Google Scholar 

  39. Mothes R, Pascual-Reguant A, Koehler R, Liebeskind J, Liebheit A, Bauherr S, et al. Distinct tissue niches direct lung immunopathology via CCL18 and CCL21 in severe COVID-19. Nat Commun. 2023;14(1):791. https://doi.org/10.1038/s41467-023-36333-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Kely Cristina Soares Bispo, Reginaldo Silva do Nascimento, Thabata Larissa Luciano Ferreira Leite, and Catia Sales de Moura for technical support. The authors would like to thank Conselho Nacional de Desenvolvimento Científico e Tecnológico, Bill and Melinda Gates Foundation, Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, and Fundação de Amparo à Pesquisa do Estado de São Paulo for financial support.

Funding

This study was funded by: Conselho Nacional de Desenvolvimento Científico e Tecnológico [401825/2020-5 and 304277/2019-3]; Bill and Melinda Gates Foundation [INV 002396]; Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, HC Convida [HC-02.18/2020]; Fundação de Amparo à Pesquisa do Estado de São Paulo [2020/10281-0]. The funders had no role in the study design, data collection, analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: TM. Methodology: NSXC, LA, JCS, JRRP, TM. Investigation: NSXC, GRJ, ECTN, JMB, LA, CSF, JSM, JCS, JRRP, RVP, MS, GSC, JDCCP, RAAM, ANDN. Visualization: NSXC, JMB, JSM. Funding acquisition: PHNS, LFFS, MD, TM.Project administration: LFFS, MD, TM, PHNS. Supervision: TM, MD, PHNS. Writing – original draft: NSXC, JCS, TM. Writing – review & editing: NSXC, MD, TM. All authors approved the final version for publication.

Corresponding author

Correspondence to Thais Mauad.

Ethics declarations

Competing interests

Authors declare that they have no competing interests.

Ethics approval and consent to participate

The National Research Ethics Commission (CONEP) approved this study (CAAE #30364720.0.0000.0068). Autopsies were performed after written informed consent by the next-of-kin.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

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

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

de Souza Xavier Costa, N., Ribeiro Júnior, G., do Nascimento, E.C.T. et al. COVID-19 induces more pronounced extracellular matrix deposition than other causes of ARDS. Respir Res 24, 281 (2023). https://doi.org/10.1186/s12931-023-02555-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12931-023-02555-7

Keywords