- 1Laboratório de Investigação em Saúde Global e Doenças Negligenciadas, Instituto Gonçalo Moniz, Salvador, Bahia, Brazil
- 2Faculdade de Biomedicina, Universidade Católica do Salvador, Salvador, Bahia, Brazil
- 3Laboratório de Pesquisa Clínica e Translacional, Instituto Gonçalo Moniz, Salvador, Bahia, Brazil
- 4Faculdade de Farmácia, Universidade Federal da Bahia, Salvador, Bahia, Brazil
- 5Departamento de Bioquímica e Biofísica, Instituto de Ciências da Saúde, Universidade Federal da Bahia, Salvador, Bahia, Brazil
- 6Programa de Pós-graduação em Imunologia, Instituto de Ciências da Saúde, Universidade Federal da Bahia, Salvador, Bahia, Brazil
- 7Departamento de Hematologia, Hospital Universitário Professor Edgard Santos, Salvador, Bahia, Brazil
Introduction: Coronavirus disease 2019 (COVID-19) has been widely associated with intense systemic inflammation, endothelial injury, and a high incidence of thrombotic complications, which together contribute to disease severity and poor clinical outcomes. While endothelial dysfunction, dysregulated cytokine production, and oxidative stress are recognized features of severe COVID-19, the direct impact of circulating factors from infected individuals on endothelial cell behavior remains insufficiently characterized. Here, we examined how serum from patients with severe COVID-19 and from convalescent individuals modulates endothelial activation, inflammatory responses, and oxidative stress using human umbilical vein endothelial cells as an in vitro model.
Methods: Venous blood samples were collected from individuals with severe COVID-19 (n = 13), convalescent patients (n = 11), and healthy volunteers (n = 7) during the initial phase of the COVID-19 pandemic. Human umbilical vein endothelial cells (HUVEC) were maintained in culture and exposed to 15% serum from each study group after a period of serum deprivation. The expression of genes associated with endothelial activation, thrombosis, inflammation, and oxidative stress was analyzed by quantitative real-time PCR at defined time points. In addition, the endothelial secretory profile was evaluated in cell culture supernatants using multiplex bead-based immunoassays. Statistical analyses were performed using one-way ANOVA followed by appropriate post hoc tests, receiver operating characteristic (ROC) curve analysis to assess the discriminatory capacity of biomarkers, and multivariate linear regression to identify factors associated with disease severity.
Results and discussion: We investigated the role of the endothelium in modulating the cytokine storm in severe COVID-19. HUVEC were stimulated with serum from patients with severe COVID-19, convalescent individuals, and healthy volunteers. Stimulation with serum from severe cases induces significant increases in VCAM1, F3, PROCR, IL6, IL12A, NFE2L2, HMOX1, GPX1, and GSR expression within 60 minutes. Antioxidant genes SOD1 and CAT were upregulated later, after 120 minutes. HUVEC stimulated with severe COVID-19 sera showed increased levels of sICAM-1, sVCAM-1, P-selectin, sE-selectin, PECAM-1, tissue factor, thrombomodulin, and a broad range of cytokines and growth factors, such as IL-1α, IL-1Ra, IL-5, IL-6, IL-10, IL-12(p40), IL-18, IL-27, TNF-α, TGF-α, FGF-2, G-CSF, M-CSF, FLT-3L, fractalkine, eotaxin, MIG, IP-10, MIP-1β, MDC, GROa and PDGF-AB/BB. In contrast, convalescent sera induced fewer markers, specifically IL-12(p40), IL-18, FGF-2, MIP-1β, MDC, GROa, and PDGF-AB/BB, while HV sera induced significant increases in IL-12(p40), IL-27, TNF-α, VEGF, MDC, eotaxin, and GROa. ROC curve analysis revealed that P-selectin and MIP-1β levels clearly distinguish severe cases from HV. When comparing severe and convalescent groups, we observed increases in IL-27, TGF-α, sVCAM-1, IL-1α, and G-CSF levels. Furthermore, Multivariable logistic regression analysis associated disease severity with decreased IL-10 and increased MIP-1β, sICAM-1, and P-selectin.
Conclusion: These findings suggest that HUVEC serves as a promising biological sensor for detecting inflammatory responses in COVID-19 patients and shows the crucial role of the endothelium in sustaining the cytokine storm that contributes to patient severity and mortality.
1 Introduction
Coronavirus disease 2019 (COVID-19), declared a pandemic by the World Health Organization, is caused by the SARS-CoV-2 virus and ranges from mild respiratory symptoms to severe acute respiratory distress syndrome, often resulting in high morbidity and mortality among hospitalized patients. Previous studies have indicated that approximately 15% of infected patients face a heightened risk of thromboembolic events (1, 2).
COVID-19 is characterized as an inflammatory disease marked by hypoxia, pulmonary edema, and cytokine storm, with complications intensifying in the presence of underlying conditions such as diabetes, hypertension, and atherosclerosis. These conditions can exacerbate the disease, frequently leading to multiple organ failure (3, 4).
Patients exhibiting mild COVID-19 symptoms, such as fever and cough, often do not require hospitalization unless they experience dyspnea or hypoxemia (oxygen saturation ≤94% or need for oxygenation or ventilatory support), which are indicative of severe disease states (5, 6).
SARS-CoV-2 may trigger an exacerbated immune response known as a cytokine storm, characterized by a significant rise in pro-inflammatory cytokines, such as interleukin (IL)-1 and IL-6, chemokines, and a T helper (Th) 1 response (7–9). This response can activate monocytes, induce tissue factor expression, and release other procoagulant substances, along with platelet activation (10–12). Associated endothelial dysfunction includes glycocalyx degradation and loss of the endothelium’s anti-inflammatory, antithrombotic, and permeability-regulating functions (13–15).
Endothelial damage is closely linked to increased organ failure risks in COVID-19 patients, with studies demonstrating how the virus triggers endothelial cell dysregulation, enhances pro-inflammatory cytokine production, oxidative stress, and expression of adhesion molecules such as ICAM1, VCAM1, and E-selectin (16–19). These changes lead to increased immune cell recruitment and endothelial cell activation, intensifying the inflammatory response (20–24).
Excessive activation of the endothelium and platelets, alongside an exacerbated immune response, often results from an unbalanced reaction between the infectious agent and the host, leading to oxidative system activation and contributing to COVID-19’s pathogenesis (25).
Oxidative stress is a critical factor in COVID-19’s pathogenesis, occurring when there is an imbalance between the production of reactive oxygen species (ROS) and the body’s antioxidant capacity. This imbalance can cause cellular and tissue damage, triggering inflammatory and apoptotic processes (26, 27). Key molecules, such as nuclear factor erythroid 2-related factor 2 (NFE2L2 or NRF2), glutathione, catalase, hemopexin, and superoxide dismutase, regulate this process and are vital for managing oxidative stress in severe pneumonia cases (28–32).
The interaction of endothelial cells with SARS-CoV-2 in vivo, particularly how they react to the virus initially present in lung epithelial cells, is still not fully elucidated. This study used human umbilical cord vein endothelial cells (HUVEC) as a model to explore their potential involvement in cytokine storm, oxidative stress, and endothelial dysfunction in patients with COVID-19.
2 Methods
2.1 Collection of samples of individuals with COVID-19 and healthy volunteers
Venous blood samples from patients with severe COVID-19 (Severe, n = 13), convalescents (Conv., n = 11), and healthy volunteers (HV, n = 7) were obtained after written consent from all study participants or their legal guardians. The protocol was conducted according to the 1975 Helsinki statement and Brazilian ethical guidelines (466-CNS-2012) and was approved by the institutional review board of the Fiocruz/Bahia Foundation (protocol number 41130920.3.0000.0040). Samples from patients with severe disease, convalescent individuals, and HV were obtained from the Intensive Care Unit (ICU) of the Salvador/BA Suburban Hospital, Irecê/BA Hospital, and Instituto Gonçalo Moniz (IGM), Salvador/BA, respectively. Severely ill patients were defined as individuals with a confirmed COVID-19 diagnosis who exhibited a severe clinical profile and required ICU admission. Convalescent individuals were defined as those who had COVID-19, were treated, and were discharged, while HV were individuals who never had COVID-19. All individuals included in this study had not received COVID-19 vaccines for, and all samples were collected during the initial period of the pandemic.
2.2 Cell culture
Human umbilical vein endothelial cells (HUVEC) were kindly donated by Ph.D. Ana Moretti and M.D. Heraldo Possolo de Souza, from the Faculdade de Medicina da USP/São Paulo, were cultivated in cell culture vials with a 25 cm2 surface area (Jet Biofil, Guangzhou, China) containing 5 mL RPMI 1640 (Gibco, New York, NY, USA) supplemented with 10% bovine fetal serum (10% SFB) heat inactivated (Gibco, New York, NY, USA) and antibiotics (100 U/mL penicillin and 10 mg/mL streptomycin) (Sigma Aldrich, St. Louis, MO, USA). Experiments were performed using cells at passages 1 to 3, and cultures were maintained in a humidified atmosphere at 37 °C with 5% CO2.
2.3 HUVEC stimulation
For the experiments, 2.5 × 105 cells/well (0.38 mL) were cultivated in 24-well plates (Costa, Corning, NY, USA) for 24 h under the above conditions. After confluence, the HUVEC were deprived of RPMI 1640 (SFB 2%) for 12 h and stimulated with medium RPMI 1640 (SFB 2%) containing 15% serum from patients with severe COVID-19 (n = 13), convalescents (n = 11), and HV (n = 7) for different periods of time, depending on the test.
2.4 Gene expression assays
RNA was extracted from HUVEC after 30 min, 60 min, and 120 min of stimulation with severe, convalescent patient, and HV serum using the RNeasy Mini-Kit (Qiagen, USA), according to the manufacturer’s specifications. The concentration and purity of the extracted RNA were determined based on spectrophotometric readings using a Nanodrop 2000 (Thermo Fisher Scientific, Rockford, IL, USA) with Absorbance A260/280. Reverse cDNA synthesis (500 ng) was performed using the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific, Rockford, IL, USA), according to the manufacturer’s specifications. The synthesized cDNAs were then diluted to a concentration of 2.5 ng/µL. Real-time PCR (qRT-PCR) was performed using an ABI Prism 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). To do this, mixtures containing SYBR® Green PCR Master Mix plus target gene (8 μL) and cDNA (2 μL) were used. The target genes evaluated in the races were: tissue factor (F3), protein 1 of vascular cellular adhesion (VCAM1), intercellular adhesion molecule 1 (ICAM1), endothelin 1 (EDN1), endothelial protein receptor (PROCR), interleukin 6 (IL6), interleukin 12 alpha (IL12A), alpha tumor necrosis factor (TNFA), nuclear factor erythroid 2 related to factor 2 (NFE2L2), heme oxygenase 1 (HMOX1), superoxide dismutase 1 (SOD1), catalase (CAT), glutathione peroxidase 1 (GPX1), glutathione S-reductase (GSR). Glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH) and Tubulin alpha-1C chain gene (TUBA1C) were used as endogenous controls for the PCR reactions. All evaluated target primers were used at a concentration of 250 pM.
2.5 Dosage of HUVEC supernatants stimulated with the sera
The cells were stimulated with sera from patients with severe, convalescent patients, and controls after 24 h of deprivation. Subsequently, the supernatants were collected and stored at −80 °C for subsequent analysis. The dosages were performed using the Milliplex® Human Cytokine/Chemokine/Growth Factor Panel to 48 Plex Kit-Hcyta-60K-PX48 (Merck/Sigma Aldrich, USA) and Milliplex® MAP Human Cardiovascular Disease Magnetic Bead Panel 2 (Merck/Sigma Aldrich, USA) according to the manufacturer’s specifications.
2.6 Statistical analysis
All experiments were performed in triplicates. GraphPad Prism v.6.0 (San Diego, CA, USA) and IBM SPSS Statistics v.21.0 (Armonk, NY, USA) were used for the statistical analyses. Statistical significance was set at p <0.05. One-way ANOVA followed by Tukey’s post-hoc test was applied to test the variance among multiple groups. Receiver operatoring characteristics (ROC) curve analysis was performed to evaluate the sensitivity and specificity of the biomarker concentrations in the discrimination of health conditions (severe vs. convalescent or severe vs. HV). Moreover, Multivariate Linear Regression (MLR) analyses were performed to investigate the weight of each variable, which presented a significant difference in association analyses (severe vs. HV) when combined.
3 Results
3.1 Serum from patients with severe COVID-19 induced significant expression of genes associated with endothelial dysfunction, thrombosis, inflammation, and antioxidant response in HUVECs
Gene expression analysis targeting markers associated with endothelial dysfunction, thrombosis, and inflammation revealed that serum from patients with severe COVID-19 significantly upregulated the expression of VCAM1 (1.43 ± 0.31-fold; p <0.01), F3 (1.59 ± 0.31-fold; p < 0.0001), PROCR (3.32 ± 0.65-fold; p < 0.0001), IL6 (1.72 ± 0.49-fold; p <0.01), and IL12A (2.26 ± 0.60-fold; p < 0.0001) in HUVECs after 60 min of stimulation (Figures 1, 2). We highlighted the most pronounced increase in the PROCR gene. In contrast, serum from convalescent patients and HV did not induce these genes. When we analyzed the effect of sera from patients with severe COVID-19 on the TNFA gene, we noticed that there was a significant reduction in the expression of this gene (0.55 ± 0.31-fold; p <0.05) after 30 min of incubation, followed by an increase (1.15 ± 0.29-fold) equivalent to the negative control (medium; 1.09 ± 0.14-fold) after 60 min, and a new significant decrease (0.58 ± 0.38-fold; p <0.05) after 120 min (Figure 1C). Additionally, stimulation with these sera showed no significant changes in EDN1 expression (Figure 2D).
Figure 1. Gene expression profile related to inflammation in HUVEC stimulated with serum from severe COVID-19 patients, convalescents, and healthy volunteers (HV). HUVEC were stimulated at intervals of 30 min, 60 min, and 120 min with serum from distinct groups: severe (n = 12), convalescents (n = 11), HV (n = 7), as well as culture medium serving as the basal control (n = 6). The data presented are mean ± standard deviation. Gene expression values were determined by relative quantification using the equation: Fold-change = 2−Δ(ΔCT), where ΔCT = CTtarget − CTendogenous and Δ (ΔCT) = ΔCTtreatment − ΔCTcontrol (medium). Each gene’s expression was normalized, expressed as a fold increase over the control. Kruskal–Wallis test, p < 0.05; Dunn’s post-test: serum vs medium, *p < 0.05; **p < 0.01; ****p < 0.0001. The existing differences between the different groups were represented by a bar, where #p < 0.05; ##p < 0.01; ###p < 0.001; and ####p < 0.0001. IL6, interleukin 6 gene (A); IL12A, interleukin 12 alpha subunit gene (B); and TNFA, tumor necrosis factor alpha gene (C).
Figure 2. Gene expression profile related to endothelial dysfunction and thrombosis in HUVEC stimulated with serum from severe COVID-19 patients, convalescents, and healthy volunteers (HV). HUVEC were stimulated for 30 min, 60 min, and 120 min with serum from different groups: severe (n = 12), convalescent (n = 11), HV (n = 7), and culture medium (n = 6). Each value corresponds to the mean ± standard deviation. Expression values were determined by relative quantification using the following expression: Fold-change = 2−Δ(ΔCT), where ΔCT = CTtarget − CTendogenous and Δ(ΔCT) = ΔCTtreatment − ΔCTcontrol (medium). Data are normalized as expression folded over control for each gene. Kruskal–Wallis test, p < 0.05; Dunn’s post-test: sera vs medium, *p < 0.05; **p < 0.01; ****p < 0.0001. The existing differences between the different groups were represented by a bar, where #p < 0.05; ##p < 0.01; ###p < 0.001; and ####p < 0.0001. VCAM1, vascular cell adhesion molecule 1 gene (A); ICAM1, intercellular adhesion molecule 1 gene (B); F3, tissue factor gene (C); EDN1, endothelin gene (D); and PROCR, endothelial protein C receptor (E).
Concerning the genes involved in the cellular antioxidant response, serum from patients with severe COVID-19 prompted increases in NFE2L2 (1.49 ± 0.38-fold; p < 0.01), HMOX1 (1.35 ± 0.30-fold; p < 0.01), GPX1 (1.44 ± 0.17-fold; p < 0.001), and GSR (1.58 ± 0.19-fold; p < 0.01) after 60 min of stimulation (Figure 3). The expression of SOD1 and CAT occurred later (after 120 min) post-stimulation. Serum from severe cases resulted in increased SOD1 (1.85 ± 0.35-fold; p <0.001) and CAT (1.78 ± 0.19-fold; p <0.0001) levels, while serum from convalescent patients induced SOD1 and CAT increases of 1.97 ± 0.32-fold (p <0.001) and 1.83 ± 0.35-fold (p <0.01), respectively (Figure 3).
Figure 3. Gene expression profile related to antioxidant response in HUVEC stimulated with serum from patients with severe COVID-19, convalescents and healthy volunteers (HV). HUVEC were stimulated for 30 min, 60 min, and 120 min with serum from severe (n = 12), convalescent (n = 11), HV (n = 7), as well as with culture medium (n=6) as a control. Each value corresponds to the mean ± standard deviation. Expression values were determined by relative quantification using the following expression: Fold-change = 2−Δ(ΔCT), where ΔCT = CTtarget − CTendogenous and Δ(ΔCT) = ΔCTtreatment − ΔCTcontrol (medium). Data are normalized as expression folded over control for each gene. Kruskal–Wallis test, p < 0.05; Dunn’s post-test: sera vs medium, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. The existing differences between the different sera were represented by a bar, where #p < 0.05; ##p < 0.01; ###p < 0.001; and ####p < 0.0001. NFE2L2, factor-related erythroid nuclear factor 2 gene (A); HMOX1, heme oxygenase-1 gene (B); SOD1, superoxide dismutase 1 gene (C); CAT, catalase gene (D); GPX1, glutathione peroxidase1 gene (E); and GSR, glutathione S-reductase gene (F).
3.2 Serum from patients with severe COVID-19 and convalescents induce the production of markers associated with endothelial dysfunction and thrombosis in HUVEC
HUVEC supernatants stimulated with serum from patients with severe COVID-19 exhibited elevated levels of sICAM-1, sVCAM-1, P-selectin, sE-selectin, PECAM-1, tissue factor, and thrombomodulin compared to those stimulated with control medium. When contrasted with serum from HV, serum from patients with severe COVID-19 significantly increased the production of sICAM-1 in HUVEC supernatants (Figure 4). Serum from convalescent patients specifically enhanced thrombomodulin secretion in HUVEC (Figure 4G).
Figure 4. Profile of endothelial dysfunction and thrombosis markers in HUVEC supernatants stimulated with serum from severe COVID-19 patients, convalescents and healthy volunteers (HV). HUVEC were stimulated with serum from severe (n = 13), convalescent (n = 11), HV (n = 7), as well as with medium alone (n = 6) as a control. and culture medium (n = 6). This stimulation occurred over 24 h following serum deprivation. Each value corresponds to the mean ± standard deviation. Kruskal–Wallis test, p < 0.05; Dunn’s post-test: sera vs medium, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. The existing differences between the different sera were represented by a bar, where #p < 0.05; ##p < 0.01. sICAM-1, soluble intercellular adhesion molecule 1 (A); sVCAM-1, soluble vascular cell adhesion molecule 1 (B); P-selectin (C); sE-Selectin, soluble E-Selectin (D); PECAM-1, platelet endothelial cell adhesion molecule (E); tissue factor (F); and thrombomodulin (G).
3.3 Serum from patients with severe COVID-19 and convalescents stimulate the production of pro-inflammatory cytokines and mediators in HUVEC
HUVEC supernatants, when stimulated with serum from patients with severe COVID-19 compared to control medium, showed a marked increase in a range of cytokines including IL-1α, IL-1Ra, IL-5, IL-6, IL-10, IL-12(p40), IL-18, IL-27, TNF-α, TGF-α, FGF-2, G-CSF, M-CSF, FLT-3L, fractalkine, MIG, eotaxin, IP-10, MIP-1β, MDC, and GROa (Figures 5, 6).
Figure 5. Cytokine profile in HUVEC supernatants stimulated with serum from severe COVID-19, convalescents and healthy volunteers (HV). HUVEC were stimulated with serum from severe (n = 13), convalescent (n = 11), and HV (n = 7) along with a culture medium (n = 6). This stimulation was conducted for 24 h following serum deprivation. Each value corresponds to the mean ± standard deviation. Kruskal–Wallis test, p < 0.05; Dunn’s post-test: sera vs medium, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. The existing differences between the different sera were represented by a bar, where #p < 0.05; ##p < 0.01. IL-1α, interleukin 1 alpha (A); IL-1Ra, interleukin 1 alpha receptor antagonist (B); IL-3, interleukin 3 (C); IL-5, interleukin 5 (D); IL-6, interleukin 6 (E); IL-10, interleukin 10 (F); IL-12(p40), interleukin 12 p40 subunit (G); IL-18, interleukin 18 (H); IL-27, interleukin 27 (I); and TNF-α, tumor necrosis factor alpha (J).
Figure 6. Quantification of growth factors and chemokines in HUVEC supernatants stimulated with serum from severe COVID-19 patients, convalescents and healthy volunteers (HV). HUVEC were stimulated with serum from severe (n = 13), convalescent (n = 11), and HV (n = 7), as well as culture medium (n = 6). This stimulation occurred over 24 h following serum deprivation. Each value corresponds to the mean ± standard deviation. Kruskal–Wallis test, p < 0.05; Dunn’s post-test: sera vs medium, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. The existing differences between the different sera were represented by a bar, where ##p < 0.01. TGF-α, Transforming growth factor alpha (A); VEGF-β, vascular endothelial growth factor B (B); FGF-2 (or bFGF), basic fibroblast growth factor (C); G-CSF, granulocyte colony stimulating factor (D); M-CSF, monocyte colony stimulating factor (E); FLT-3L, FMS-related tyrosine kinase 3 ligand (F); fractalkine (CX3CL1) (G); MIG, interferon gamma-induced monokine (H); eotaxin (I); IP-10, interferon gamma-induced protein 10 (J); MIP-1β, macrophage inflammatory protein-1 beta (K); MDC, macrophage-derived chemokine (L); GROa (or CXCL1), growth-regulated alpha protein (M); and PDGF-AB/BB, platelet-derived growth factor-AB/BB (N).
In contrast, HUVEC supernatants stimulated with serum from convalescent patients compared to the control medium, exhibited an increase in a more limited set of cytokines, specifically IL-12(p40), IL-18, FGF-2, eotaxin, MIP-1β, MDC, and GROa (Figures 5, 6). Furthermore, stimulation with convalescent serum led to a significant reduction in IL-3 levels (Figure 5).
When comparing HUVEC supernatants stimulated with serum from HV to the control medium, there was a significant increase in the levels of IL-12(p40), IL-18, IL-27, TNF-α, VEGF, eotaxin, MDC, GROa, and PDGF-AB/BB (Figures 5, 6). The results for the other cytokines, growth factors, and chemokines that did not show significant changes are detailed in Supplementary Figure 1.
Furthermore, serum from patients with severe COVID-19 significantly stimulated the production of IL-1ra and IL-10 compared to HV. Elevated levels of MIP-1b levels were also observed in HUVEC supernatants after stimulation with serum from severe COVID-19 cases compared to HV. Moreover, serum from convalescent patients demonstrated reduced production of TGF-α and G-CSF compared with that from patients with severe COVID-19. Notably, convalescent serum showed significant decreases in sICAM-1, sVCAM-1, P-selectin, IL-1α, IL-1ra, IL-10, IL-27, TNF-α, TGF-α, and G-CSF compared to serum from severe COVID-19 patients.
3.4 Assessing biomarker efficacy in distinguishing health conditions
ROC curve analysis was conducted to assess the sensitivity and specificity of biomarker concentrations in discriminating health conditions. Comparing the severe group and HV, the analyses demonstrated statistically significant curves for IL-1Ra (AUC = 1.000; cut off >18.50 pg/mL), IL-10 (AUC = 0.956; cut off >12.00 pg/mL), sICAM-1 (AUC = 1.000; cut off >16,700 ng/mL), MIP-1β (AUC = 0.885; cut off >18.50 pg/mL), and P-selectin (AUC = 0.940; cut off >1,025 ng/mL). Comparisons between the severe and convalescent groups also showed significant curves for IL-1α (AUC = 0.839; cut off >16.25 pg/mL), IL-1Ra (AUC = 0.982; cut off >21.50 pg/mL), IL-10 (AUC = 0.920; cut off >11.50 pg/mL), IL-27 (AUC = 0.913; cut off >21.75 pg/mL), TGF-α (AUC = 0.909; cut off >75.50 pg/mL), G-CSF (AUC = 0.811; cut off >18.00 pg/mL), sICAM-1 (AUC = 1.000; cut off >14,350 ng/mL), and sVCAM-1 (AUC = 0.857; cut off >19,675 ng/mL) (Figure 7; Supplementary Table 1).
Figure 7. ROC curve analysis of biomarkers in HUVEC stimulated by serum from patients with different health conditions (COVID-19). The graphs show the sensitivity and specificity of markers that have statistical significance of patients with severe COVID-19 vs. HV and severe COVID-19 vs. HV. AUC, Area under the curve; CI, confidence interval. IL-1α, interleukin 1 alpha (A); IL-1Ra, interleukin 1 alpha receptor antagonist (B); IL-10, interleukin 10 (C); IL-27, interleukin 27 (D); TGF-α, transforming growth factor alpha (E); G-CSF, granulocyte colony stimulating factor (F); MIP-1β, macrophage inflammatory protein-1 beta (G); sICAM-1, soluble intercellular adhesion molecule 1 (H); sVCAM-1, soluble vascular cell adhesion molecule 1 (I); and P-selectin (J).
Table 1 presents the MLR data. The proposed model significantly predicted health conditions (F (5, 19) = 8.947, p = 0.001). The adjusted R2 of 0.507 indicates that the model accounts for 50.7% of the variance in health conditions. Notably, reduced IL-10 levels were significantly associated with severe disease (β = −1.130), whereas increased levels of MIP-1β, sICAM-1, and P-selectin were strongly linked to severe disease (β = 1.145, 0.724, and 0.464, respectively).
Table 1. Multivariate linear regression model comparing HUVECs stimulated by serum from patients with different health conditions (COVID-19).
4 Discussion
The findings indicate that HUVEC exhibit distinct responses when stimulated with serum from patients with severe COVID-19, convalescents, and HV, demonstrating the endothelium’s critical role in the body’s reaction to COVID-19. Notably, serum from patients with severe COVID-19 induced significant upregulation of genes associated with endothelial dysfunction, thrombosis, inflammation, and antioxidant responses in HUVECs. This suggests robust endothelial activation and dysfunction in severe COVID-19 cases, which could contribute to the increased thromboembolic and inflammatory complications observed in these patients.
Initially, we assessed the expression of IL6, IL12A, and TNFA, which are key pro-inflammatory genes implicated in the cytokine storm associated with COVID-19. Our results revealed that serum from patients with severe COVID-19 significantly upregulated IL6 and IL12A expression in HUVEC within 60 min of stimulation.
Although there was no significant increase in TNFA expression at this early stage, we noted a marked decrease in its expression 120 min after stimulation with serum from both severe and convalescent patients compared to that from HV. Additionally, substantial reductions were observed in the expression of IL6 and IL12A at later time points. This rapid decline may be attributed to the swift synthesis and subsequent conversion of mRNA into protein products, illustrating that the serum from patients induces a distinctly different response in HUVEC than that induced by serum from HV. The increase in IL-6, IL-12, and TNF-α levels in HUVEC supernatants further supports our hypothesis regarding the dynamic regulation of inflammatory mediators in response to COVID-19.
Serum from patients with severe COVID-19 significantly upregulated genes involved in the cellular antioxidant response in HUVEC, including NFE2L2, HMOX1, SOD1, CAT, GPX1, and GSR. The later expression of SOD1 and CAT after 120 min could indicate a delayed but essential antioxidant defense mechanism against oxidative stress induced by severe inflammation. NFE2L2 (also known as NRF2) plays a pivotal role in the cytoprotective response to oxidative stress, although the mechanisms by which it combats SARS-CoV-2 infection remain poorly understood (33).
Genetic variation in NFE2L2 and KEAP1 may influence acute lung injury by modifying the activation of this protective pathway, potentially accounting for varying susceptibilities to different clinical stages of SARS-CoV-2 infection (32). This highlights the critical importance of maintaining NFE2L2 integrity to trigger antioxidant enzyme induction, providing a cytoprotective response against pro-oxidant insults from the virus.
Catalase is widely distributed in epithelial cells and leukocytes, among others, and protects them against oxidative damage, regulates cytokine production, and slows the replication of SARS-CoV-2 (34). Mo et al. (35) showed that enhancing SOD production could prevent plasma cell apoptosis upon re-exposure and dramatically prolong SARS-CoV-specific antibody production. Another study showed that serum concentrations of SOD and GPx were significantly higher in patients with COVID-19 than in controls, indicating an activated antioxidant response (36). Conversely, significantly lower GSR and elevated levels of IL-10 in COVID-19 patients suggest deficiencies in antioxidant capabilities and immune function, respectively. Our data indicated a distinct response in HUVEC, with an increase in GSR expression. Supporting this finding, Naghashpour et al. (37) showed that serum from patients with severe COVID-19 strongly induced GSR expression, underscoring the potential role of enhanced antioxidant responses in severe disease stages.
The elevated expression of genes such as VCAM1, F3, PROCR, IL6, and IL12A within 60 min of stimulation indicates an acute and rapid endothelial response to the severe inflammatory state induced by COVID-19 infection. The lack of gene induction by serum from convalescent patients and healthy volunteers underscores the unique pathological environment presented by severe COVID-19, which appears to persist beyond the acute phase of disease.
The production of soluble markers, such as sICAM-1, sVCAM-1, P-selectin, sE-selectin, PECAM-1, tissue factor, and thrombomodulin, further illustrates extensive endothelial activation and the potential for increased vascular permeability and thrombosis. This is particularly notable in severe cases compared to convalescent patients, where only thrombomodulin showed a significant increase.
Vascular cell adhesion molecule-1 (VCAM-1) and ICAM-1 levels were found to be elevated in the blood of patients with mild COVID-19 and dramatically elevated in severe cases, decreasing during the convalescence phase (38). This trend may explain the increased expression of the corresponding genes in the present study. Pro-inflammatory cytokines, such as IL-6, IL-1β, and TNF-α, are known to induce F3 gene expression in both immune and non-immune cells, potentially contributing to the thromboembolic events observed in COVID-19 patients, particularly in severe cases (39). Furthermore, a meta-analysis showed that tissue factor levels were positively associated with senescence and hypercoagulation gene signatures in patients with COVID-19. This study underscores the role of tissue factors as a critical link between inflammation, thrombosis, and senescence in respiratory viral infection (40).
The distinct cytokine profile, characterized by elevated levels of pro-inflammatory cytokines in HUVECs supernatants stimulated by serum from severe COVID-19 patients, suggests a hyperactive inflammatory response that could exacerbate endothelial injury and dysfunction. This hyperactivity contrasts with the more restrained cytokine production observed in HUVECs stimulated with serum from convalescent patients, indicating a gradual resolution of the inflammatory response.
HUVEC substantially express the PROCR gene, which encodes the endothelial protein C receptor (EPCR). EPCR is associated with anticoagulant and anti-inflammatory roles, regulating hemostasis and inflammation (41). Some studies have reported that mutations in the PROCR gene are associated with the development of thrombosis (42, 43). EPCR is an important controller of the protein C anticoagulant pathway, binding to protein C and potentiating its activation by the thrombomodulin-thrombin complex. EPCR binds to protein C and activated protein C (APC) with high affinity.
A recent study by Won et al. (44) showed that pulmonary endothelial cells from patients who died of severe COVID-19 had increased expression of von Willebrand procoagulant factor (VWF) and decreased expression of the anticoagulant thrombomodulin and EPCR. Although we did not measure EPCR in the culture supernatant, the remarkable increase in PROCR gene expression suggests that critically ill patients respond by inducing the antithrombotic system to compensate for the thrombogenic condition in COVID-19.
We then assessed the cytokine levels in HUVEC supernatants stimulated with sera. Serum from patients with severe COVID-19 induced increased levels of IL-1α, IL-1Ra, IL-5, IL-6, IL-10, IL-12(p40), IL-18, IL-27, and TNF-α.
Previous reports have shown that the pathogenesis of COVID-19 is supported by a potent inflammatory response involving various mediators, such as IL-6 and IL-10. These pleiotropic cytokines are produced at inflammation sites by different cell types, including lymphocytes, endothelial cells, and epithelial cells, and are released into the circulation during sepsis and acute organ injury (45). Furthermore, IL-6, IL-12, and TNF-α are cytokines classically known to have pro-inflammatory properties and are elevated in COVID-19 (46). Studies carried out with sera from patients with various degrees of severity of COVID-19 showed that patients had higher levels of these markers (47, 48). Our results corroborate these findings, suggesting that endothelial cells may contribute to the production of these cytokines in response to stimulation of sera from patients infected with SARS-CoV-2 in severe conditions. Previous studies have shown that endothelial cells can produce cytokines after infection (49, 50).
Although none of the sera showed a significant increase in TNFA expression, the results of assaying HUVEC supernatants showed that sera from patients with severe COVID-19 significantly induced the production of TNF-α. Based on this information, it is possible to suggest that the expression of the TNFA gene may have occurred shortly after the stimuli, with subsequent production of the cytokine.
Shao et al. (51) showed that Human Mammary Epithelial Cells (HMEC) have trained immunity in response to subsequent MERS-CoV-2 infections, which contributes to upregulation of pro-inflammatory cytokines genes, such as TNFA, IL6, CSF1, CSF3 and IL32 (trained immunity marker), and indicates that trained immunity in infected HMEC can generate cytokine storms.
Elevated levels of IL-1α were associated with lung injury in patients with severe COVID-19 (52). In contrast, the IL-1Ra is a member of the IL-1 cytokine family and is secreted by various cell types, including epithelial cells, immune cells, hepatocytes, keratinocytes, and fibroblasts. It acts as a natural inhibitor of the proinflammatory effect of IL-1β (53) and is significantly associated with the severity and progression of COVID-19 (24).
Han et al. (45) showed that high levels of cytokine storms are directly associated with disease severity and are associated with IL-6, suggesting that these markers can be used as predictors for the rapid diagnosis of patients at higher risk of disease severity. In contrast, IL-10 is a Th2-type cytokine that exhibits anti-inflammatory action against cytokines such as IL-6, IL-1β, IFN-γ, and TNF-α in different cell types, preventing the maturation of dendritic cells by blocking the production of IL-12 (54).
IL-18 is involved in the activation of Th1 and natural killer cells, being able to induce IFN-γ in the presence of IL-12 or IL-15. IL-18 acts on CD4, CD8, and NK cells to induce IFN-γ production (55). Furthermore, macrophages can produce IFN-γ when activated by IL-18 and IL-12 (56).
IL-27 has a dual role, being able to stimulate the differentiation of Th1 cells that express IFN-γ, as well as IL-10-producing T cells (57). IL-27 acts on several cell types, including T and B cells, macrophages, dendritic cells, NK cells, and non-hematopoietic cells (58). The expansion of Th1 cells occurs through the transcription-mediated T-bet signaling pathway (STAT) 1 (59). In contrast, IL-27 also acts by suppressing immune responses through the inhibition of Th17 cell development and induction of IL-10 production in a STAT1- and STAT3-dependent manner (60, 61). Circulating IL-27 levels are positively correlated with severity and worse prognosis in patients with severe COVID-19 (62). Xu et al. (63) showed in a recent study that serum IL-27 levels are positively associated with severity and poor prognosis in patients with pneumonia, indicating that IL-27 may be involved in the pathophysiological process of the disease. Therefore, these data suggest that IL-27 may be involved in the pathogenesis of inflammatory lung diseases.
In individuals infected with SARS-CoV-2, Th1 and Th2 cytokines are elevated, Wolf et al. (64) reported a consistently elevated expression of type 1 interferons, Th1 and Th2 interleukins, and chemokines early in SARS-CoV-2 infection. The observed increased IL-5 levels corroborate the findings of Montazersaheb et al. (65) who reported a correlation between the severity of COVID-19 and the reduction in T cells. This leads to the overexpression of Th2-derived cytokines such as IL-4, IL-5, and IL-13 (65, 66).
Previous studies have shown that FGF2 is involved in the strong apoptotic response induced by MERS-CoV, which is fundamental for the development of increased replication by the lytic cycle in the lungs (67). Furthermore, FGF2 plays a role in intussusceptive angiogenesis and induces angiogenesis due to the marked inflammatory process and hypoxia that occur in the disease (68). Meini et al. (69) suggest that intussusceptive angiogenesis may have the function of restoring a functional vascular plexus that results from extensive endothelialitis and alveolar capillary microthrombosis in COVID-19.
Fractalkine (CX3CL1) is a chemokine produced by the activated endothelium in response to stimuli with pro-inflammatory cytokines, such as TNF-α and IL-1β (70, 71), and acts as an important regulator of endothelial function and dysfunction, including thrombosis (72). CX3CL1 may be upregulated in the endothelium during SARS-CoV-2 infection and contribute to the perpetuation of a prothrombotic cycle (72).
MIG (CXCL9) and IP-10 (CXCL10) are involved in Th1, CD8, and NK cell migration and early upregulation after SARS-CoV-2 infection (73). Higher MIG levels are detected in mildly symptomatic and severe patients than in asymptomatic patients (74). IP-10 is produced by activated bronchial and alveolar epithelial cells and is implicated in apoptosis and the development of lymphopenia and neutrophilic infiltration, and its levels increase with disease severity or death (74).
Studies have shown that patients who died with COVID-19 had elevated plasma levels of CX3CL1, CXCL9, and CXCL10, among others, compared to patients with severe and/or mild COVID-19 (73).
The Flt-3L ligand is a hematopoietic cytokine that stimulates the proliferation and differentiation of several blood cell progenitors and is an important growth factor that stimulates the growth of dendritic cells (75).
In summary, serum from patients with severe COVID-19 significantly induced the expression of genes associated with endothelial dysfunction, thrombosis, inflammation, and oxidative stress. Conversely, serum from convalescents exhibited an expression profile similar to that of healthy volunteers, with no notable increases in target gene expression at the evaluated times, except for a discrete yet significant induction of the tissue factor gene.
Additionally, HUVEC supernatants stimulated with serum from severe cases showed a substantially higher presence of cytokines, chemokines, and inflammation-related growth factors than those from convalescents and healthy volunteers.
Specifically, P-selectin and MIP-1β levels were elevated in HUVEC supernatants stimulated with serum from patients with severe COVID-19 compared to those from HV. Comparing the severe and convalescent groups, we observed an observed increase in IL-27, TGF-α, sVCAM-1, IL-1α, and G-CSF. Analysis of the ROC curve results revealed common markers between patients with severe COVID-19 and both HV and convalescent groups, including IL-1Ra, IL-10, and sICAM-1, as diagnostic tools to distinguish between health conditions with high sensitivity and specificity.
Moreover, MLR analysis reinforced the significance of IL-10, MIP-1 β, s-ICAM-1, and P-selectin as indicators of disease severity. Decreased IL-10 and increased MIP-1β, sICAM-1, and P-selectin levels were associated with greater disease severity, emphasizing the critical importance of these markers in understanding and managing COVID-19.
These findings suggest that these biomarkers can be useful for distinguishing the health status of critically ill and convalescent patients, underscoring the significant role of the endothelium in the pathophysiology of COVID-19. This role persists even in convalescent individuals, as SARS-CoV-2 infection can lead to long-term complications. However, additional research is required to substantiate this hypothesis.
5 Conclusion
Based on our findings, HUVEC demonstrates significant potential as biological sensors in response to serum stimuli from patients with COVID-19. These cells contribute to the production of cytokines, chemokines, and growth factors, which are crucial in the inflammatory response, oxidative stress, and procoagulant state observed in response to inflammatory mediators present in the serum of patients with the disease. This highlights the critical role of the endothelium in the pathophysiology of COVID-19, underscoring its importance in both disease progression and potential therapeutic targeting.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by the institutional review board of the Fiocruz/Bahia Foundation (protocol number 41130920.3.0000.0040). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
SS: Data curation, Writing – review & editing, Formal analysis, Investigation, Methodology, Writing – original draft. SY: Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. CA: Investigation, Methodology, Writing – review & editing. JS: Investigation, Methodology, Writing – review & editing. HS: Investigation, Methodology, Writing – review & editing. CB: Writing – review & editing, Resources, Visualization. TP: Visualization, Writing – review & editing. VB: Visualization, Writing – review & editing. VF: Visualization, Writing – review & editing. IL: Visualization, Writing – review & editing, Resources. MG: Visualization, Writing – review & editing, Conceptualization, Data curation, Funding acquisition, Project administration, Resources, Supervision, Validation.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by grants (N° 48116285170491) from Programa Fiocruz de Fomento à Inovação (INOVA Fiocruz). The sponsors of this study are public or nonprofit organizations that support science in general.
Acknowledgments
The authors would like to thank the Programa Fiocruz de Fomento à Inovação (INOVA Fiocruz) and the participants enrolled in the study.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1689772/full#supplementary-material
References
1. Levy JH and Iba T. COVID-19: Thrombosis, thromboinflammation, and anticoagulation considerations. Int J Lab Hematol. (2021) 43 Suppl 1:29–35. doi: 10.1111/ijlh.13500
2. Yin S, Huang M, Li D, and Tang N. Difference of coagulation features between severe pneumonia induced by SARS-CoV2 and non-SARS-CoV2. J Thromb Thrombolysis. (2021) 51:1107–10. doi: 10.1007/s11239-020-02105-8
3. Berlin DA, Gulick RM, and Martinez FJ. Severe covid-19. N Engl J Med. (2020) 383:2451–60. doi: 10.1056/NEJMcp2009575
4. Matthay MA, Aldrich JM, and Gotts JE. Treatment for severe acute respiratory distress syndrome from COVID-19. Lancet Respir Med. (2020) 8:433–4. doi: 10.1016/S2213-2600(20)30127-2
5. Gandhi RT, Lynch JB, and Del Rio C. Mild or moderate Covid-19. N Engl J Med. (2020) 383:1757–66. doi: 10.1056/NEJMcp2009249
6. Gulick RM, Pau AK, Daar E, Evans L, Gandhi RT, Tebas P, et al. National institutes of health COVID-19 treatment guidelines panel: perspectives and lessons learned. Ann Intern Med. (2024) 177:1547–57. doi: 10.7326/ANNALS-24-00464
7. Buja LM, Wolf DA, Zhao B, Akkanti B, Mcdonald M, Lelenwa L, et al. The emerging spectrum of cardiopulmonary pathology of the coronavirus disease 2019 (COVID-19): Report of 3 autopsies from Houston, Texas, and review of autopsy findings from other United States cities. Cardiovasc Pathol. (2020) 48:107233. doi: 10.1016/j.carpath.2020.107233
8. Fajgenbaum DC and June CH. Cytokine Storm. N Eng J Med. (2020) 383:2255–73. doi: 10.1056/NEJMra2026131
9. Hu B and Huang S. The cytokine storm and COVID-19. J Med Virol. (2021) 93:250–6. doi: 10.1002/jmv.26232
10. Hottz ED, Azevedo-Quintanilha IG, Palhinha L, Teixeira L, Barreto EA, Pão CRR, et al. Platelet activation and platelet-monocyte aggregate formation trigger tissue factor expression in patients with severe COVID-19. Blood. (2020) 136:1330–41. doi: 10.1182/blood.2020007252
11. Wang J, Pendurthi UR, Yi G, and Rao LVM. SARS-CoV-2 infection induces the activation of tissue factor–mediated coagulation via activation of acid sphingomyelinase. Blood. (2021) 138:344–9. doi: 10.1182/blood.2021010685
12. Mackman N. Tissue factor and COVID-19 associated thrombosis. Arteriosc Thromb Vasc Biol. (2024) 44:523–9. doi: 10.1161/ATVBAHA.123.320144
13. Yamaoka-Tojo M. Vascular endothelial glycocalyx damage in COVID-19. Int J Mol Sci. (2020) 21. doi: 10.3390/ijms21249712
14. Xu SW, Ilyas I, and Weng JP. Endothelial dysfunction in COVID-19: an overview of evidence, biomarkers, mechanisms and potential therapies. Acta Pharmacol Sin. (2023) 44:695–709. doi: 10.1038/s41401-022-00998-0
15. Tay EA, Vijayakumar V, Morales RF, Lee ES, and Teo A. Protecting the endothelial glycocalyx in COVID-19. PLoS Pathog. (2024) 20:e1012203. doi: 10.1371/journal.ppat.1012203
16. Chernyak BV, Popova EN, Prikhodko AS, Grebenchikov OA, Zinovkina LA, and Zinovkin RA. COVID-19 and oxidative stress. Biochem (Mosc). (2020) 85:1543–53. doi: 10.1134/S0006297920120068
17. Dupont A, Rauch A, Staessens S, Moussa M, Rosa M, Corseaux D, et al. Vascular endothelial damage in the pathogenesis of organ injury in severe COVID-19. Arterioscler Thromb Vasc Biol. (2021) 41:1760–73. doi: 10.1161/ATVBAHA.120.315595
18. Qian Y, Lei T, Patel Parth S, Lee Chi H, Monaghan-Nichols P, Xin H-B, et al. Direct activation of endothelial cells by SARS-CoV-2 nucleocapsid protein is blocked by simvastatin. J Virol. (2021) 95:10.1128/jvi.01396–01321. doi: 10.1128/JVI.01396-21
19. Bruni F, Charitos P, Lampart M, Moser S, Siegemund M, Bingisser R, et al. Complement and endothelial cell activation in COVID-19 patients compared to controls with suspected SARS-CoV-2 infection: A prospective cohort study. Front Immunol. (2022) 13. doi: 10.3389/fimmu.2022.941742
20. Copin MC, Parmentier E, Duburcq T, Poissy J, and Mathieu D. Time to consider histologic pattern of lung injury to treat critically ill patients with COVID-19 infection. Intensive Care Med. (2020) 46:1124–6. doi: 10.1007/s00134-020-06057-8
21. Jin Y, Yang H, Ji W, Wu W, and Chen S. Virology, epidemiology, pathogenesis, and control of COVID-19. Viruses. (2020) 12:. doi: 10.3390/v12040372
22. Pons S, Fodil S, Azoulay E, and Zafrani L. The vascular endothelium: the cornerstone of organ dysfunction in severe SARS-CoV-2 infection. Crit Care. (2020) 24:353. doi: 10.1186/s13054-020-03062-7
23. Sardu C, Gambardella J, Morelli MB, Wang X, Marfella R, and Santulli G. Hypertension, thrombosis, kidney failure, and diabetes: is COVID-19 an endothelial disease? A comprehensive evaluation of clinical and basic evidence. J Clin Med. (2020) 9:. doi: 10.20944/preprints202004.0204.v1
24. Yang Y, Shen C, Li J, Yuan J, Yang M, Wang F, et al. Exuberant elevation of IP-10, MCP-3 and IL-1ra during SARS-CoV-2 infection is associated with disease severity and fatal outcome. medRxiv. (2020). doi: 10.1101/2020.03.02.20029975
25. Blanco-Melo D, Nilsson-Payant BE, Liu WC, Uhl S, Hoagland D, Møller R, et al. Imbalanced host response to SARS-CoV-2 drives development of COVID-19. Cell. (2020) 181:1036–1045.e1039. doi: 10.1016/j.cell.2020.04.026
26. Jian Z, Guo H, Liu H, Cui H, Fang J, Zuo Z, et al. Oxidative stress, apoptosis and inflammatory responses involved in copper-induced pulmonary toxicity in mice. Aging (Albany NY). (2020) 12:16867–86. doi: 10.18632/aging.103585
27. Forcados GE, Muhammad A, Oladipo OO, Makama S, and Meseko CA. Metabolic implications of oxidative stress and inflammatory process in SARS-CoV-2 pathogenesis: therapeutic potential of natural antioxidants. Front Cell Infect Microbiol. (2021) 11:654813. doi: 10.3389/fcimb.2021.654813
28. Reichard JF, Motz GT, and Puga A. Heme oxygenase-1 induction by NRF2 requires inactivation of the transcriptional repressor BACH1. Nucleic Acids Res. (2007) 35:7074–86. doi: 10.1093/nar/gkm638
29. Baird L and Yamamoto M. The molecular mechanisms regulating the KEAP1-NRF2 pathway. Mol Cell Biol. (2020) 40. doi: 10.1128/MCB.00099-20
30. Zuo Y, Yalavarthi S, Shi H, Gockman K, Zuo M, Madison JA, et al. Neutrophil extracellular traps in COVID-19. JCI Insight. (2020) 5. doi: 10.1172/jci.insight.138999
31. Gayashani Sandamalika WM, Kwon H, Lim C, Yang H, and Lee J. The possible role of catalase in innate immunity and diminution of cellular oxidative stress: Insights into its molecular characteristics, antioxidant activity, DNA protection, and transcriptional regulation in response to immune stimuli in yellowtail clownfish (Amphiprion clarkii). Fish Shellfish Immunol. (2021) 113:106–17. doi: 10.1016/j.fsi.2021.03.022
32. Soto ME, Fuentevilla-Álvarez G, Palacios-Chavarría A, Vázquez RRV, Herrera-Bello H, Moreno-Castañeda L, et al. Impact on the clinical evolution of patients with COVID-19 pneumonia and the participation of the NFE2L2/KEAP1 polymorphisms in regulating SARS-CoV-2 infection. Int J Mol Sci. (2023) 24:415. doi: 10.3390/ijms24010415
33. Kensler TW, Wakabayashi N, and Biswal S. Cell survival responses to environmental stresses via the Keap1-Nrf2-ARE pathway. Annu Rev Pharmacol Toxicol. (2007) 47:89–116. doi: 10.1146/annurev.pharmtox.46.120604.141046
34. Qin M, Cao Z, Wen J, Yu Q, Liu C, Wang F, et al. An antioxidant enzyme therapeutic for COVID-19. Adv Mater. (2020) 32:2004901. doi: 10.1002/adma.202004901
35. Mo L-H, Luo X-Q, Ma K, Shao J-B, Zhang G-H, Liu Z-G, et al. Superoxide dismutase prevents SARS-coV-2-induced plasma cell apoptosis and stabilizes specific antibody induction. Oxid Med Cell Longevity. (2022) 2022:5397733. doi: 10.1155/2022/5397733
36. Golabi S, Ghasemi S, Adelipour M, Bagheri R, Suzuki K, Wong A, et al. Oxidative stress and inflammatory status in COVID-19 outpatients: A health center-based analytical cross-sectional study. Antioxidants. (2022) 11:606. doi: 10.3390/antiox11040606
37. Naghashpour M, Ghiassian H, and Mobarak S. Profiling serum levels of glutathione reductase and interleukin-10 in positive and negative-PCR COVID-19 outpatients: A comparative study from southwestern Iran. J Med Virol. (2022) 94:1457–64. doi: 10.1002/jmv.27464
38. Tong M, Jiang Y, Xia D, Xiong Y, Zheng Q, Chen F, et al. Elevated expression of serum endothelial cell adhesion molecules in COVID-19 patients. J Infect Dis. (2020) 222:894–8. doi: 10.1093/infdis/jiaa349
39. Subramaniam S, Kothari H, and Bosmann M. Tissue factor in COVID-19-associated coagulopathy. Thromb Res. (2022) 220:35–47. doi: 10.1016/j.thromres.2022.09.025
40. Nguyen D, Jeon H-M, and Lee J. Tissue factor links inflammation, thrombosis, and senescence in COVID-19. Sci Rep. (2022) 12:19842. doi: 10.1038/s41598-022-23950-y
41. Esmon CT. Structure and functions of the endothelial cell protein C receptor. Crit Care Med. (2004) 32:S298–301. doi: 10.1097/01.CCM.0000126128.64614.81
42. Simmonds RE and Lane DA. Structural and functional implications of the intron/exon organization of the human endothelial cell protein C/activated protein C receptor (EPCR) gene: comparison with the structure of CD1/major histocompatibility complex alpha1 and alpha2 domains. Blood. (1999) 94:632–41. doi: 10.1182/blood.V94.2.632
43. Biguzzi E, Merati G, Liaw PC, Bucciarelli P, Oganesyan N, Qu D, et al. A 23bp insertion in the endothelial protein C receptor (EPCR) gene impairs EPCR function. Thromb Haemost. (2001) 86:945–8. doi: 10.7326/ANNALS-24-00464
44. Won T, Wood MK, Hughes DM, Talor MV, Ma Z, Schneider J, et al. Endothelial thrombomodulin downregulation caused by hypoxia contributes to severe infiltration and coagulopathy in COVID-19 patient lungs. EBioMedicine. (2022) 75:103812. doi: 10.1016/j.ebiom.2022.103812
45. Han H, Ma Q, Li C, Liu R, Zhao L, Wang W, et al. Profiling serum cytokines in COVID-19 patients reveals IL-6 and IL-10 are disease severity predictors. Emerg Microbes Infect. (2020) 9:1123–30. doi: 10.1080/22221751.2020.1770129
46. Hasanvand A. COVID-19 and the role of cytokines in this disease. Inflammopharmacology. (2022) 30:789–98. doi: 10.1007/s10787-022-00992-2
47. Del Valle DM, Kim-Schulze S, Huang HH, Beckmann ND, and Nirenberg S. An inflammatory cytokine signature predicts COVID-19 severity and survival. Nat Med. (2020) 26:1636–43. doi: 10.1038/s41591-020-1051-9
48. Tjan LH, Furukawa K, Nagano T, Kiriu T, Nishimura M, Arii J, et al. Early differences in cytokine production by severity of coronavirus disease 2019. J Infect Dis. (2021) 223:1145–9. doi: 10.1093/infdis/jiab005
49. Gadient RA and Patterson PH. Leukemia inhibitory factor, interleukin 6, and other cytokines using the GP130 transducing receptor: roles in inflammation and injury. Stem Cells. (1999) 17:127–37. doi: 10.1002/stem.170127
50. Robles JP, Zamora M, Adan-Castro E, and Siqueiros-Marquez L. The spike protein of SARS-CoV-2 induces endothelial inflammation through integrin α5β1 and NF-κB signaling. J Biol Chem. (2022) 298:101695. doi: 10.1016/j.jbc.2022.101695
51. Shao Y, Saredy J, Xu K, Sun Y, Saaoud F, Drummer CT, et al. Endothelial immunity trained by coronavirus infections, DAMP stimulations and regulated by anti-oxidant NRF2 may contribute to inflammations, myelopoiesis, COVID-19 cytokine storms and thromboembolism. Front Immunol. (2021) 12:653110. doi: 10.3389/fimmu.2021.653110
52. Liu Y, Yang Y, Zhang C, Huang F, Wang F, Yuan J, et al. Clinical and biochemical indexes from 2019-nCoV infected patients linked to viral loads and lung injury. Sci China Life Sci. (2020) 63:364–74. doi: 10.1007/s11427-020-1643-8
53. Arend WP, Malyak M, Guthridge CJ, and Gabay C. Interleukin-1 receptor antagonist: role in biology. Annu Rev Immunol. (1998) 16:27–55. doi: 10.1146/annurev.immunol.16.1.27
54. Costela-Ruiz VJ, Illescas-Montes R, Puerta-Puerta JM, Ruiz C, and Melguizo-Rodríguez L. SARS-CoV-2 infection: The role of cytokines in COVID-19 disease. Cytokine Growth Factor Rev. (2020) 54:62–75. doi: 10.1016/j.cytogfr.2020.06.001
55. Kaplanski G. Interleukin-18: Biological properties and role in disease pathogenesis. Immunol Rev. (2018) 281:138–53. doi: 10.1111/imr.12616
56. Munder M, Mallo M, Eichmann K, and Modolell M. Murine macrophages secrete interferon gamma upon combined stimulation with interleukin (IL)-12 and IL-18: A novel pathway of autocrine macrophage activation. J Exp Med. (1998) 187:2103–8. doi: 10.1084/jem.187.12.2103
57. Yoshida H and Hunter CA. The immunobiology of interleukin-27. Annu Rev Immunol. (2015) 33:417–43. doi: 10.1146/annurev-immunol-032414-112134
58. Iwasaki Y, Fujio K, Okamura T, and Yamamoto K. Interleukin-27 in T cell immunity. Int J Mol Sci. (2015) 16:2851–63. doi: 10.3390/ijms16022851
59. Takeda A, Hamano S, Yamanaka A, Hanada T, Ishibashi T, Mak TW, et al. Cutting edge: role of IL-27/WSX-1 signaling for induction of T-bet through activation of STAT1 during initial Th1 commitment. J Immunol. (2003) 170:4886–90. doi: 10.4049/jimmunol.170.10.4886
60. Kastelein RA, Hunter CA, and Cua DJ. Discovery and biology of IL-23 and IL-27: related but functionally distinct regulators of inflammation. Annu Rev Immunol. (2007) 25:221–42. doi: 10.1146/annurev.immunol.22.012703.104758
61. Stumhofer JS, Silver JS, Laurence A, Porrett PM, Harris TH, Turka LA, et al. Interleukins 27 and 6 induce STAT3-mediated T cell production of interleukin 10. Nat Immunol. (2007) 8:1363–71. doi: 10.1038/ni1537
62. Angioni R, Sánchez-Rodríguez R, Munari F, Bertoldi N, Arcidiacono D, Cavinato S, et al. Age-severity matched cytokine profiling reveals specific signatures in Covid-19 patients. Cell Death Dis. (2020) 11:957. doi: 10.1038/s41419-020-03151-z
63. Xu Z, Wang XM, Cao P, Zhang C, Feng CM, Zheng L, et al. Serum IL-27 predicts the severity and prognosis in patients with community-acquired pneumonia: a prospective cohort study. Int J Med Sci. (2022) 19:74–81. doi: 10.7150/ijms.67028
64. Wolf A, Khimani F, Yoon B, Gerhart C, Endsley D, Ray AK, et al. The mechanistic basis linking cytokine storm to thrombosis in COVID-19. Thromb Update. (2022) 8:100110. doi: 10.1016/j.tru.2022.100110
65. Montazersaheb S, Hosseiniyan Khatibi SM, Hejazi MS, Tarhriz V, Farjami A, Ghasemian Sorbeni F, et al. COVID-19 infection: an overview on cytokine storm and related interventions. Virol J. (2022) 19:92. doi: 10.1186/s12985-022-01814-1
66. Mehta P, Mcauley DF, Brown M, Sanchez E, Tattersall RS, and Manson JJ. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet. (2020) 395:1033–4. doi: 10.1016/S0140-6736(20)30628-0
67. Yeung ML, Yao Y, Jia L, Chan JF, Chan KH, Cheung KF, et al. MERS coronavirus induces apoptosis in kidney and lung by upregulating Smad7 and FGF2. Nat Microbiol. (2016) 1:16004. doi: 10.1038/nmicrobiol.2016.4
68. Presta M, Andrés G, Leali D, Dell’era P, and Ronca R. Inflammatory cells and chemokines sustain FGF2-induced angiogenesis. Eur Cytokine Netw. (2009) 20:39–50.
69. Meini S, Giani T, and Tascini C. Intussusceptive angiogenesis in Covid-19: hypothesis on the significance and focus on the possible role of FGF2. Mol Biol Rep. (2020) 47:8301–4. doi: 10.1007/s11033-020-05831-7
70. Ahn SY, Cho CH, Park KG, Lee HJ, Lee S, Park SK, et al. Tumor necrosis factor-alpha induces fractalkine expression preferentially in arterial endothelial cells and mithramycin A suppresses TNF-alpha-induced fractalkine expression. Am J Pathol. (2004) 164:1663–72. doi: 10.1016/S0002-9440(10)63725-X
71. Imaizumi T, Yoshida H, and Satoh K. Regulation of CX3CL1/fractalkine expression in endothelial cells. J Atheroscler Thromb. (2004) 11:15–21. doi: 10.5551/jat.11.15
72. Rivas-Fuentes S, Valdés VJ, Espinosa B, Gorocica-Rosete P, and Salgado-Aguayo A. Could SARS-CoV-2 blocking of ACE2 in endothelial cells result in upregulation of CX3CL1, promoting thrombosis in COVID-19 patients? Med Hypotheses. (2021) 151:110570. doi: 10.1016/j.mehy.2021.110570
73. Hue S, Beldi-Ferchiou A, Bendib I, Surenaud M, Fourati S, Frapard T, et al. Uncontrolled innate and impaired adaptive immune responses in patients with COVID-19 acute respiratory distress syndrome. Am J Respir Crit Care Med. (2020) 202:1509–19. doi: 10.1164/rccm.202005-1885OC
74. Khalil BA, Elemam NM, and Maghazachi AA. Chemokines and chemokine receptors during COVID-19 infection. Comput Struct Biotechnol J. (2021) 19:976–88. doi: 10.1016/j.csbj.2021.01.034
Keywords: COVID-19 severity, endothelial dysfunction, HUVEC, cytokine storm, thromboinflammation, oxidative stress
Citation: Santana SS, Yahouédéhou SCMA, Adanho CSA, Silva JR, Santos HM, Barbosa CG, Pitanga TN, Borges VM, Fortuna V, Lyra IM and Goncalves MS (2026) Assessing the endothelium’s role in COVID-19 severity using the HUVEC model. Front. Immunol. 16:1689772. doi: 10.3389/fimmu.2025.1689772
Received: 20 August 2025; Accepted: 08 December 2025; Revised: 25 November 2025;
Published: 08 January 2026.
Edited by:
Hugo Caire Castro-Faria-Neto, Oswaldo Cruz Foundation (Fiocruz), BrazilReviewed by:
Dan Predescu, Rush University, United StatesShetty Ravi Dyavar, Adicet Bio, Inc, United States
Copyright © 2026 Santana, Yahouédéhou, Adanho, Silva, Santos, Barbosa, Pitanga, Borges, Fortuna, Lyra and Goncalves. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marilda Souza Goncalves, bWFyaWxkYS5nb25jYWx2ZXNAZmlvY3J1ei5icg==
†These authors have contributed equally to this study
Hayna Malta Santos3