Abstract
COVID-19 has affected more than half a billion people worldwide, with more than 6.3 million deaths, but the pathophysiological mechanisms involved in lethal cases and the host determinants that determine the different clinical outcomes are still unclear. In this study, we assessed lung autopsies of 47 COVID-19 patients and examined the inflammatory profiles, viral loads, and inflammasome activation. Additionally, we correlated these factors with the patient’s clinical and histopathological conditions. Robust inflammasome activation, mediated by macrophages and endothelial cells, was detected in the lungs of lethal cases of SARS-CoV-2. An analysis of gene expression allowed for the classification of COVID-19 patients into two different clusters. Cluster 1 died with higher viral loads and exhibited a reduced inflammatory profile than Cluster 2. Illness time, mechanical ventilation time, pulmonary fibrosis, respiratory functions, histopathological status, thrombosis, viral loads and inflammasome activation significantly differed between the two clusters. Our data demonstrated two distinct profiles in lethal cases of COVID-19, thus indicating that the balance of viral replication and inflammasome-mediated pulmonary inflammation led to different clinical outcomes. We provide important information to understand clinical variations in severe COVID-19, a process that is critical for decisions between immune-mediated or antiviral-mediated therapies for the treatment of critical cases of COVID-19.
Introduction
After the COVID-19 pandemic, coronaviruses became frequent viral agents of acute respiratory distress syndrome (ARDS). Between 2019 and 2021, there were more than 255,000,000 confirmed cases and 5,127,000 deaths caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide (WHO, 2020). Severe ARDS is characterized by major pulmonary involvement, respiratory distress, systemic thrombosis, and death (Batah and Fabro, 2021; Bösmüller et al., 2021; Polak et al., 2020). The activation of the NLRP3 inflammasome has been described in response to SARS-CoV-2 infections (Cama et al., 2021; Eisfeld et al., 2021; Junqueira et al., 2022; Lucas, 2020; Rodrigues et al., 2021; Sefik et al., 2022), and this process contributes to excessive inflammation and poor clinical outcome. The recruitment of immune cells to the lungs of infected individuals culminates in the excessive release of cytokines that cause structural damage to the lungs (Kommoss et al., 2020; Xu et al., 2020). Moreover, the major reported lung pathologies of severe COVID-19 include diffuse alveolar damage (DAD), acute fibrinoid organizing pneumonia (AFOP), and chronic interstitial pneumonia (Batah and Fabro, 2021; Bösmüller et al., 2021; Polak et al., 2020; Sauter et al., 2020). NLRP3 inflammasome activation also occurs in the lung parenchyma of these patients (Rodrigues et al., 2021; Toldo et al., 2021). However, it is still unknown how inflammasome activation promotes the exacerbated inflammatory process during SARS-CoV-2 infection and how it impacts viral replication and the development of important clinical characteristics associated with death in COVID-19 patients. In the present study, postmortem lung tissue samples from 47 fatal COVID-19 patients were examined for inflammasome activation and gene expression; in addition, these factors were correlated with the clinical conditions of the patients to understand the molecular mechanisms underlying the pathological processes that lead to the death of these patients. We demonstrated the presence of robust inflammasome activation in the lungs of patients infected with SARS-CoV-2 and defined the specific cell types that contribute most to inflammasome activation. An analysis of gene expression in pulmonary tissues allowed for the classification of COVID-19 patients into two different clusters: one cluster with patients who died with higher viral loads and reduced inflammatory profiles, which opposes the data from the other cluster. Our data establish that the magnitude of inflammasome activation contributes to the pathology of COVID-19 by balancing a thrombotic versus a fibrotic process that impacts the clinical outcome of the disease and patient death.
Results
Lethal cases of COVID-19 develop severe acute respiratory distress syndrome (ARDS) and higher inflammasome activation
We evaluated inflammasome activation in the lung parenchyma of 47 patients who died from SARS-CoV-2 infection from April to July 2020. These patients were infected with the ancestral strain of SARS-CoV-2 before the development of the COVID-19 variants of concern and before the use of vaccinations. As noninfected controls, we evaluated inflammasome activation in biopsies of the benign area of the lungs from patients who died due to lung adenocarcinoma (referred as uninfected controls). The data in Table 1 show the demographic and clinical characteristics of the 47 patients. As reported in the early COVID-19 cases (that occurred in 2020), we found a high frequency of comorbidities, such as hypertension, obesity, diabetes, heart and lung disease. Furthermore, the patients exhibited altered values of CRP, D-dimers, LDH, creatinine, urea, AST, ALT, PT (INR), and blood glucose levels. The patients had an average illness time of 18 days, and the majority of the patients (76.59%) used MV, with an average of 12.38 days of use of mechanical ventilation. The patients had considerably altered PaO2 and PaO2/FiO2 characterizing a moderate to severe ARDS condition and consistent with the respiratory statuses. The histopathological analyses of COVID-19 patients demonstrated a high fibrotic phase of DAD, OP, and pneumonitis (Table 1). To analyze the histological features of patients’ lungs, we measured the lung parenchyma area in the histological sections to assess loss of airspace. In support of the data described above, we observed that COVID-19 patients had a higher parenchyma area (loss of airspace), suggesting that SARS-CoV-2 triggers robust inflammatory infiltration to the lung in severe cases of the disease (Supplementary Fig. 1A-C). These data are in agreement with histopathological findings showing an intense inflammatory process in patients who died from SARS-CoV-2 (Batah and Fabro, 2021; Hariri et al., 2021). To assess inflammasome activation in the lungs of COVID-19 patients, the lung parenchyma area and cell count were scanned in histological sections by using multiphoton microscopy. The parenchyma area was calculated by using ImageJ software, and inflammasome activation was microscopically scored by the presence of characteristic NLRP3 or ASC puncta/specks (Hauenstein et al., 2015). The parenchyma area was used to normalize all of the counts between the patients. Inflammasome activation was quantified by counting ASC and NLRP3 puncta in the lung tissue of these patients. We observed that patients with SARS-CoV-2 infection had robust inflammasome activation, as shown by the abundant presence of ASC and NLRP3 puncta in the lungs (Supplementary Fig. 1D-G). Representative images of ASC and NLRP3 puncta in the patients’ lungs are shown (Supplementary Fig. 1E, G). Importantly, we observed ASC colocalization in nearly all of the scored NLRP3 puncta, thus confirming that these puncta structures that were abundantly found in the patient’s lungs are indeed the NLRP3/ASC inflammasome (Supplementary Fig. 1H, I). Collectively, these data indicate that COVID-19 patients have high inflammasome activation and that these patients evolved to a fibrotic phase of DAD, OP, and pneumonitis.
Macrophages and endothelial cells contribute to inflammasome activation in lethal cases of COVID-19
To comprehensively investigate the cell types operating during SARS-CoV-2 infection, we aimed to assess which cell types promote inflammasome activation in the lung of COVID-19 patients. We found abundant ASC and NLRP3 puncta in macrophages (CD64+), endothelial cells (CD34+), type 1 pneumocytes (PDPN+), and type 2 pneumocytes (SFTPC+) (Figure 1A-H). We scored ASC and NLRP3 puncta in CD64+, CD34+, PDPN+, and SFTPC+ and compared them with uninfected controls. We observed that in the lungs of COVID-19 patients, both ASC and NLRP3 puncta were significantly higher in macrophages (Figure 1I, M) and endothelial cells (Figure 1J, N) than in controls. In contrast, in the assessment of type 1 and type 2 pneumocytes, we found that ASC and NLRP3 puncta did not differ from uninfected controls, suggesting that these cells contribute less to overall inflammasome activation (Figure 1K, L, O, P). Analyzing the total numbers of macrophages, endothelial cells, type 1 pneumocytes, and type 2 pneumocytes in the lungs of patients who died from SARS-CoV-2 infection, we did not find differences in total counts of these cell types comparing COVID-19 and uninfected controls (Supplementary Fig. 2A-D). Together, these data demonstrate the main contribution of macrophages and endothelial cells but not pneumocytes, to inflammasome activation during SARS-CoV-2 infection.
Subsequently, we performed immunohistochemistry analyses of the lung parenchyma to assess the expression of inflammasome components and observed variable expression of NLRP3, ASC, cleaved GSDMD, IL-1β, and caspase-1 in COVID-19 patients (Supplementary Fig. 3A-E). To assess which lung parenchyma cells expressed these inflammasome components, we used immunohistochemistry to evaluate the expression of ASC, NLRP3, caspase-1, cleaved GSDMD, and IL-1β on macrophages, endothelial cells, and type 2 pneumocytes. We also stained these tissues with anti-Spike of SARS-CoV-2 and found that virus-associated endothelial cells, macrophages, and type 2 pneumocytes expressed ASC, NLRP3, Caspase-1, and IL-1β (Supplementary Fig. 3F-H).
Disease progression in fatal cases of COVID-19 occurs with decreasing viral load and increasing inflammasome activation
To assess whether inflammasome activation was associated with specific patient clinical conditions, we performed Pearson’s correlations in COVID-19 patients. Even though all of these patients died, we observed a strong negative correlation between viral load and the time of the disease (symptoms onset to death) (Figure 2A, B). In addition, we observed a positive correlation between the amount of NLRP3 and ASC puncta with the time of disease (Figure 2C, D), thus suggesting that, in general, whereas the viral load is reduced, inflammasome activation increases during hospitalization in these lethal cases of COVID-19. We also observed significant negative correlations between NLRP3 puncta and PaO2/FiO2 (Figure 2E), suggesting that inflammasome activation is related to an overall worsening pulmonary function. We detected no statistically significant correlations between NLRP3 or ASC puncta and viral load in the tissues (Supplementary Fig. 4).
Inflammasome activation and pulmonary viral loads define two distinct clinical outcomes in lethal cases of COVID-19
The imbalance of inflammatory and anti-inflammatory processes leads to an excessive release of cytokines into the systemic circulation with potentially deleterious consequences, including systemic inflammatory response syndrome (SIRS), circulatory shock, multiorgan dysfunction syndrome (MODS), and death (Sinha et al., 2020). Several studies have reported the association of SARS-CoV-2 with hyperinflammatory syndrome related to disease severity (Chen et al., 2021; Huang et al., 2005; Qin et al., 2020), and inflammasome activation may be associated with Cytokine Release Syndrome (CRS) (Cui and Zhang, 2020; Dolinay et al., 2012; Olajide et al., 2021). Thus, to correlate inflammasome activation with the hyperinflammatory profile observed in lethal cases of COVID-19, we analyzed the expression of genes related to the inflammatory process and inflammasome activation. We evaluated gene expression in the lungs of COVID-19 patients and observed that several inflammatory cytokine genes positively correlate with genes involved in inflammasome activation (Supplementary Fig. 5A). We did not observe a significant difference in the expression of these genes when we compared COVID-19 patients with uninfected controls (Supplementary Fig. 5B-U), possibly due to the large dispersion of data observed in COVID-19 patients.
Due to this dispersion and high variation in gene expression detected in COVID-19 patients, we performed an unsupervised heatmap constructed with data from relative gene expression and revealed the formation of two clusters in COVID-19 patients (Figure 3). Cluster 1 was characterized by a higher viral load and lower expression of inflammasome and inflammatory genes. In contrast to Cluster 1, Cluster 2 was comprised of patients who died with a lower viral load and increased expression of inflammasome and inflammatory genes (Figure 3A). To gain insights into key differences between these two clusters, we analyzed PaO2/FiO2 and A-a O2 gradient kinetics in these patients. We observed that patients in Cluster 2, who had overall increased inflammation, had worsening pulmonary function compared to those in Cluster 1 (Figure 3B, C). Basic demographic information stratified by cluster is provided in Table 2. Strikingly, patients who belonged to Cluster 2 had a long illness time (Figure 3D), worse pulmonary function (as indicated by the PaO2/FiO2 and A-a O2 gradients) (Figure 3E, F), greater inflammasome activation, as measured by NLRP3 puncta formation (Figure 3G), and increased area of pulmonary parenchyma, indicating loss or airways (Figure 3H). Images of the lung parenchyma of three representative patients from each cluster are shown (Figure 3I-J). Importantly, we found that patients belonging to Cluster 2 had an overall lower viral load than patients from Cluster 1; this was quantified both by RT-PCR for N2 and E gene expression (Figure 4A, B), and by expression of Spike protein (Figure 4C-E). Next, we assessed radiological analyses of 15 patients from Cluster 1 and 31 patients from Cluster 2 by comparing the first and last chest X-ray images (CXR) and observed overall worsening pulmonary conditions in patients from Cluster 2 (Supplementary Fig. 6A). Representative images of CXR of two patients from Cluster 1, with Mild/Moderate pneumonia (stable conditions) (Supplementary Fig. 6B) and two patients from Cluster 2 with severe pneumonia (worsening) are shown (Supplementary Fig. 6C). Table 3 summarizes the general imaging evaluation findings for the patients in Cluster 1 and Cluster 2. Taken together, our data suggest the existence of two distinct groups of patients who succumbed to COVID-19, one group with lower viral loads, higher inflammation, and worse pulmonary conditions than the other group.
Differential regulation of fibrinolysis-related genes and growth factors dictate the induction of fibrosis or disseminated intravascular coagulation in fatal cases of COVID-19
Histopathological analyses of the lungs from lethal cases of COVID-19 belonging to Clusters 1 and 2 allowed us to detect an increased development of fibrosis in patients from Cluster 2 (Figure 5A, and Table 2). By contrast, patients belonging to Cluster 1 exhibited a markedly increased disseminated intravascular coagulation (Figure 5B, and Table 2), according to the parameters recommended by the International Society of Thrombosis and Hemostasis (ISTH) (Taylor et al., 2001). In support of our analyses demonstrating increased disseminated intravascular coagulation, we found reduced platelet counts in patients from Cluster 1 (Figure 5C, and Table 2), supporting that vascular dysfunctions impact clinical outcomes in patients belonging to Cluster 1.
To further investigate the biological pathways involved in the induction of fibrosis in Cluster 2 patients, we assessed the expression of fibrosis-related genes in the patient’s lungs. We observed an increased expression of genes related to fibrinolysis, inflammatory cytokines, nitric oxide-induced, extracellular matrix deposition, and extracellular matrix remodeling in patients belonging to Cluster 2 (Figure 5D). A volcano plot indicates genes differentially expressed in samples from Cluster 2 compared to Cluster 1, including Plg (Figure 5E), Plat (Figure 5F), Plau (Figure 5G); Ccr2 (Figure 5H), Ccl3 (Figure 5I); Cdkna1 (Figure 5J); Serpina1 (Figure 5K), Timp2 (Figure 5L) and Fgf1 (Figure 5M).
To further illustrate the differences between fibrosis and thrombosis in patients from these different clusters, we performed Masson-Goldner staining to assess erythrocytes and clot formation (bright red) and collagen deposition (green) in Patients’ lungs. By assessing three patients from each cluster, we found an increased collagen deposition in Cluster 2 patients, suggesting an increased fibrotic process (Figure 5N). By contrast, samples from Cluster 1 show increased clot formation, suggesting increased thrombotic processes (Figure 5O). Together, these data support the hypothesis that fatal cases of COVID-19 progress via induction of disseminated intravascular coagulation, culminating in a thrombotic process or via a fibrotic process, which is aggravated by the inflammasome-induced exacerbated inflammation.
Discussion
COVID-19 is significantly lethal in nonvaccinated individuals, and although inflammation and cytokine storm are associated with poor clinical outcomes, the mechanisms underlying dysregulated inflammatory processes are unknown. The revelation that exacerbated inflammasome activation contributes to COVID-19 pathology (Junqueira et al., 2022; Rodrigues et al., 2021; Sefik et al., 2022) advanced the understanding of disease pathology, but a significant proportion of the lethal cases progresses with lower inflammasome activation. Our analysis of 47 fatal COVID-19 cases allowed for the classification of the patients into two groups. Cluster 2 showed a remarkably high inflammasome activation and hyperexpression of inflammatory genes with increased pulmonary fibrosis and worsened respiratory functions. By contrast, patients belonging to Cluster 1 died faster, with higher viral loads, reduced inflammatory process, and increased disseminated intracellular coagulation. Our data reveal two distinct profiles in lethal cases of COVID-19, thus indicating that the balance of viral replication and inflammasome-mediated pulmonary inflammation led to different clinical outcomes.
Patients belonging to Cluster 2, died of poor respiratory functions and increased fibrosis induced by the excessive inflammatory process. Inflammatory cytokines released upon inflammasome activation have been previously linked to the development of pulmonary fibrosis (Cho et al., 2020; Hussain et al., 2014; Lv et al., 2018; Meng et al., 2015; Meng et al., 2019; Sun et al., 2016). In addition, cytokines such as IL-1β, IL-18, and IL-1α have been described as triggering the activation of fibroblasts and stimulating the synthesis and accumulation of type I collagen, TIMP, collagenase, and PGE2 (Hussain et al., 2014; Postlethwaite et al., 1988), which may contribute to the triggering of the deleterious effects of inflammasomes in patients’ lungs. Interestingly, patients from Cluster 2 died with low viral loads and higher inflammasome activation; in many patients from this cluster, viral loads were not detected, suggesting that SARS-CoV-2 per see may not be required for continuous inflammasome activation in the lungs of these patients. How NLRP3 inflammasome is activated in the tissues is still unknown, it is possible that Damage-Associated Molecular Patterns (DAMPs) participate in this process. Alternatively, it is possible that viral proteins produced before viral elimination or undetectable levels of SARS-CoV-2 present in the tissues are sufficient to activate NLRP3 inflammasome. Future studies will be required to clarify these questions. We also tested whether the mechanical ventilation used in COVID-19 patients would interfere with inflammasome activation. This hypothesis is supported by data indicating that mechanical ventilation-induced hyperoxia can induce potassium efflux through the P2X7 receptor, thus leading to inflammasome activation and the secretion of proinflammatory cytokines (Dolinay et al., 2012; Jones et al., 2014; Kolliputi et al., 2010; Kuipers et al., 2012; Wu et al., 2013; Zhang et al., 2014). This does not appear to be the case in our study, as when we separated the COVID-19 patients into two groups according to the use or non-use of mechanical ventilation, we detected no differences in histopathological analyses and inflammasome activation (Table 4). Moreover, we confirmed that patients who underwent mechanical ventilation had a longer illness time, higher hypertension, less CRP and albumin, higher amounts of urea, and worsened pulmonary function (according to the PaO2/FiO2 and A-an O2 gradients). Thus, our data do not support the hypothesis that mechanical ventilation is directly associated with inflammasome activation.
The mechanisms underlying dysregulated inflammatory processes and the cell types involved in inflammasome activation in COVID-19 are largely unknown. Macrophages are possibly the most effective cell types that trigger inflammasome activation (Broz and Dixit, 2016); according to this information, our data showed that macrophages in the lungs of COVID-19 patients are effectively infected and strongly induce inflammasome activation. These data are in agreement with previously published articles indicating high levels of viral RNA in lung monocytes and macrophages of COVID-19 patients (Delorey et al., 2021; Pontelli et al., 2022), as well as pronounced inflammasome activation in macrophages from COVID-19 patients (Junqueira et al., 2022; Rodrigues et al., 2021; Sefik et al., 2022). In addition to inflammasome activation in macrophages, the analysis of inflammasome activation in different pulmonary cell types indicates that endothelial cells exposed to SARS-CoV-2 express inflammasome proteins and contain active NLRP3/ASC puncta, thus demonstrating inflammasome activation. Despite reports indicating that endothelial cells are not productively infected by SARS-CoV-2 (Schimmel et al., 2021), precursors of hematopoietic and endothelial cells stimulated with SARS-CoV-2 spike protein increase the expression of AIM2, NLRP1, NLRP3, IL1Β, and ASC and trigger caspase-1 activation (Kucia et al., 2021; Ratajczak et al., 2021). In addition, inflammasome activation in endothelial cells has been previously reported (Paul et al., 2021; Xiang et al., 2011; Xu et al., 2013; Yang et al., 2016). These observations are consistent with our data indicating the expression of inflammasome components and inflammasome activation in CD34+ cells in the lungs of COVID-19 patients. Although these cells were stained positive for the SARS-CoV-2 spike protein, it is unknown whether infection and viral replication are required for inflammasome activation in these cells. It is possible that exosomes, released from COVID-19-infected cells, trigger NLRP3 inflammasome in endothelial cells (Sur et al., 2022).
Alternatively, It is possible that spike stimulation is sufficient for inflammasome activation, which is a feature that would explain reports indicating damage and dysfunction of these cells during COVID-19 and the characterization of COVID-19 as an endothelial disease (Hottz et al., 2020; Libby and Lüscher, 2020; Liu et al., 2021; Nuovo et al., 2021; Varga et al., 2020; Ward et al., 2021). Importantly, vascular endothelium is actively involved in the regulation of inflammation and thrombus formation. This is particularly important in COVID-19 because the interplay between inflammasome activation in macrophages and the endothelium can lead to pyroptotic macrophages releasing tissue factor (TF), which is an essential initiator of coagulation cascades that are frequently observed in severe cases of COVID-19 (Campos et al., 2021; Wu et al., 2019; Zhang et al., 2021). These data support the association of inflammasomes with the coagulopathy that is observed in COVID-19.
In this study, we demonstrated specific cell types involved in inflammasome activation in patients’ lungs and identified two distinct profiles in lethal cases of COVID-19. The revealed balance of viral-induced intravascular coagulation versus inflammasome-mediated pulmonary inflammation contributes to our understanding of disease pathophysiology and may contribute to decisions between immune-mediated or antiviral-mediated therapies for the treatment of critical cases of COVID-19.
Methods
Samples and Study Approval
Minimally invasive autopsy was performed on 47 patients diagnosed with SARS-CoV-2 at the Hospital das Clínicas da Faculdade de Medicina de Ribeirão Preto da Universidade de São Paulo, Brazil (Ribeirão Preto, SP, Brazil) from April to July, 2020, by Serviço de Patologia (SERPAT). Minimally invasive autopsy was done at bedside through post-mort surgical lung biopsy by a matching 14-gauge cutting needle (Magnum Needles, Bard) and a biopsy gun (Magnum, Bard). Moreover, a 3 cm incision on the more affected side of the chest between the fourth and fifth ribs were also used to provide extra lung tissue. All tissue samples were embedded in paraffin and fixed in formalin (Formalin-Fixed Paraffin-Embedded, FFPE). Lung tissue from biopsy of Lung Adenocarcinoma patients (n=5) were obtained from SERPAT. The project was approved by the Research Ethics Committee of FMRP/USP under protocol n° 4,089.567.
Immunohistochemistry (IHC)
Tissue section in paraffin blocks were tested by immunohistochemistry using antibodies for the detection of NLRP3 (clone D2P5E; 1:3,000; Cell Signaling), ASC (1:2,000; Adipogen AL177), SFTPC (1:200, ThermoFisher PA5-71680), CD68 (1:200, Dako M0814), CD34 (1:500, Zeta Z2063ML), PDPN (1:200, ThermoFisher 14-9381-82), Anti-cleaved N-terminal GSDMD (1:200, Abcam ab215203), Caspase-1 (1:200, Abcam ab207802), IL-1β (1:200, Abcam ab2105) for the in situ detection of these inflammasome proteins. The Sequential Immunoperoxidase Labeling and Erasing (SIMPLE) technique was used to evaluate all markers in the same tissue, as previously published (Glass et al., 2009). Briefly, after incubation with primary antibody (overnight at 4°C), slides were incubated with immune peroxidase polymer anti-mouse visualization system (SPD-125, Spring Bioscience, 345 Biogen) and with chromogen-substrate AEC peroxidase system kit (SK-4200, 346 Vector Laboratories, Burlingame, CA). After high resolution scanning by a VS120 Olympus microscope, the coverslips were removed in PBS and the slides were dehydrated in an ethanol gradient to 95% ethanol. The slides were incubated in a series of ethanol to erase the AEC marking. Afterwards the slides were rehydrated and the antibodies were removed with an incubation for 2 min in a solution of 0.15 mM 351 KMnO4/0.01 M H2SO4, followed immediately by a wash in distilled water. Tissues were then remarked.
Immunofluorescence
The slides were incubated with the primary antibodies, rabbit anti-human NLRP3 mAb (clone D2P5E; 1:300; Cell Signaling), rabbit anti-human ASC polyclonal antibody (1:200; Adipogen AL177), SFTPC (1:200, ThermoFisher PA5-71680), CD64 (1:200, Biolegend 139304), CD34 (1:500, Zeta Z2063ML), PDPN (1:200, ThermoFisher 14-9381-82), overnight at 4°C and with the secondary antibodies Goat anti-mouse Alexa fluor-647 (Invitrogen) or Goat anti-rabbit Alexa fluor-594 (Invitrogen). Images were acquired by the Axio Observer system combined with the LSM 780 confocal device microscope at 63x magnification (Carl Zeiss).
Histological Evaluation
Paraffin-embedded lung tissues sections (3 μm) were stained by standard hematoxylin and eosin (H&E). Morphological lung injury patterns were evaluated by specialized pulmonary pathologists (ATF) blinded to clinical history. They were classified as absent or present with its extent of lung injury area by 5% cut-offs in Fibrosis, Organizing Pneumonia (OP), Acute Fibrinous and Organizing Pneumonia (AFOP), Diffuse Alveolar Damage (DAD), Cellular Pneumonitis, Thrombus Formation.
RNA extraction and Real-Time Polymerase Chain Reaction for inflammatory genes
Total RNA from fresh lung tissue of SARS-CoV-2 patients and controls was obtained using Trizol reagent, and purification was performed according to the manufacturer’s instructions. The RNA was quantified by spectrophotometry in a NanoDrop 2000c spectrophotometer. The concentration was adjusted to 1 μg/µL, and the RNA was stored at −70 ◦C until reverse transcription. For inflammatory genes qPCR the total RNA was transcribed into complementary DNA (cDNA) using a High-Capacity cDNA Reverse Transcription kit (without an inhibitor) according to the protocol provided by the manufacturer (Thermo Fisher, Carlsbad, CA, USA). The reaction was prepared in a final volume of 20.0 µL containing 4.2 µL of H2O, 2.0 µL of buffer, 2.0 µL of random primers, 0.8 µL of dNTP Mix (100 mM), 1.0 µL of reverse transcriptase (RT) enzyme and 1 µL of RNA (1 ug/µL). The solution was then placed into a thermocycler with the following program: 25 ◦C for 10 min, 37 ◦C for 120 min and 85 ◦C for 5 min. The real-time PCR was performed in 96-well plates using Sybr Green reagents (Applied Biosystems, Waltham, MA, USA) and a Quant studio real-time PCR system (Applied Biosystems, Foster City, CA, USA). The real time-RT-PCR was carried out in a total volume of 20 μl on a 96-well MicroAmp Fast Optical plate (Applied Biosystems). Each well contained 10 μl SYBR Green qPCR Master Mix (Thermofisher), 1 μl of each primer (Supplementary Table 1), 2 μl cDNA (20 ng) and 7 μl RNase free water using the following protocol: initial denaturation at 95°C for 10 min, 40 cycles of denaturation at 95°C for 15 s followed by annealing/extension at 60°C for 60 s. Each PCR was followed by a dissociation curve analysis between 60-95°C. The Ct values were analyzed by the comparative Ct (ΔΔCt) method and normalized to the endogenous control GAPDH. Fold difference was calculated as 2−ΔΔCt
Real-Time Polymerase Chain Reaction for Viral RNA
Detection and quantification of SARS-CoV-2 genes was performed with primer-probe sets for 2019-nCoV_N2 and gene E, according to US Centers for Disease Control and Prevention (Lu et al., 2020) and Charité group protocols (Corman et al., 2020). The genes evaluated (N2, E, and RNase-P housekeeping gene) were tested by one-step real-time RT-PCR using total nucleic acids extracted with TRIzol (Invitrogen). All real-time PCR assays were done on a Quant studio real-time PCR system (Applied Biosystems, Foster City, CA, USA). A total of 70 ng of RNA was used for genome amplification, adding specific primers (20 µM), and probe (5 µM), and with TaqPath 1-Step quantitative RT-PCR Master Mix (Applied Biosystems), with the following parameters: 25°C for 2 min, 50°C for 15 min, and 95°C for 2 min, followed by 45 cycles of 94°C for 5 s and 60°C for 30 s. Primers used were the following: N2 forward: 5′-TTACAAACATTGGCCGCAAA-3′, N2 reverse: 5′-GCGCGACATTCCGAAGAA-3′; N2 probe: 5′-FAM-ACAATTTGCCCCCAGCGCTTCAG-BHQ1-3′ (Lu et al., 2020); E forward: 5′-ACAGGTACGTTAATAGTTAATAGCGT-3′, E reverse: 5′-ATATTGCAGCAGTACGCACACA-3′; E probe: 5′-AM-ACACTAGCCATCCTTACTGCGCTTCG-BHQ-1-3′ (Corman et al., 2020); RNase-P forward: 5′-AGATTTGGACCTGCGAGCG-3′, RNase-P reverse: 5′-GAGCGGCTGTCTCCACAAGT-3′; and RNase-P probe: 5′-FAM-TTCTGACCTGAAGGCTCTGCGCG-BHQ-1-3′ (Lu et al., 2020). A plasmid of N2 protein and E protein was used for a standard curve construction for viral load quantification.
Chest computed tomography images acquisition and evaluation
Chest x-radiography (CXR) and computed tomography (CT) exams were performed as part of the routine clinical evaluation. Chest radiographies were performed in conventional equipment, mainly in the anteroposterior incidence. CT images were performed in multidetector scanners (Brilliance CT Big Bore 16 - Philips, Holland, or Aquilion Prime 160 - Toshiba, Japan), using similar protocols for the acquisition of high-resolution images of the lungs (Raghu et al., 2018). Patients were scanned in the supine position without the administration of intravenous contrast media. Typical acquisition parameters were: 120 kVp tube voltage, 100–140 ref mAs (Koenigkam-Santos et al., 2022), 0.3–0.7 s gantry rotation time, reconstruction matrix size of 512×512, slice thickness, and increment of 1.0 mm, using standard (soft) and hard kernel filters. Imaging exams were independently evaluated by two thoracic radiologists (MKS and DTW), blinded to clinical data, laboratory, and pathology results as described (Koenigkam-Santos et al., 2022). Divergences were solved by consensus. All CXRs available were classified for the presence and grade of viral pneumonia (Revel et al., 2020). Pulmonary disease evolution on imaging was evaluated considering all exams from the initial to the last image before death. Tomographic images were evaluated similarly to CXR images.
Statistical Analysis
For puncta quantification, all histological sections were viewed on a 63x objective for digitalizing random images using the LSM 780 system in the Axio Observer microscope, covering an area of about ∼1.7 mm2 of lung parenchyma analyzed per case. Manual counting of puncta and cells was blinded and performed using the acquired images. Morphometric analyzes were performed as described (Weibel, 1963). The quantification of expression by immunohistochemistry was performed by calculating the percentage of marked area, the scanned images were opened in the ImageJ software, using the IHQ Toolbox plugin, which consists of a semi-automatic color selection tool that selects the pixels positive for immunohistochemical marking, differentiating them from background and H&E marking, after selecting the positive pixels, the images were transformed into 8-bits and the area occupied by these tones was calculated. The distribution of the gene expression and biochemical marker and puncta count data was evaluated using the Shapiro–Wilk test; the data were analyzed using non-parametric Kruskal-Wallis, Mann-Whitney test and Spearman correlation. The data as violin plot graphs show median and quartiles. Statistical analyzes were performed using the GraphPad PRISM 5.0 program, with p<0.05 being considered statistically significant. Heatmaps were constructed using the heatmap.2 function in the R program (Project for Statistical Computing, version 3.4.1), the hierarchical clustering method used a correlation distance measure with ward.D2 and Canberra analysis.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Funding
Fundação de Amparo à Pesquisa do Estado de Sao Paulo, FAPESP grants 2013/08216-2, 2019/11342-6 and 2020/04964-8. Conselho Nacional de Desenvolvimento Científico e Tecnológico, CNPq grant 303021/2020-9. Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, CAPES grant 88887.507253/2020-00.
Competing interests
Authors declare that they have no competing interests.
Supplementary Figures and Supplementary Figure Legends
Acknowledgments
We would like to thank Maira Nakamura, Amanda Zuin and Dr. Roberta Sales for technical support.
Footnotes
The data obtained from samples of Influenza patients have been removed from this version of the manuscript. In addition, this version contains additional data from samples of COVID-19 patients (e.g., Figure 5). Supplemental files updated.