Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Fraction of plasma exomeres and low-density lipoprotein cholesterol as a predictor of fatal outcome of COVID-19

  • Tatiana Usenko ,

    Roles Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft

    tatiana.s.usenko@gmail.com

    Affiliations Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia, Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia, Kurchatov Genome Center—PNPI, Saint-Petersburg, Russia

  • Valentina Miroshnikova,

    Roles Visualization, Writing – review & editing

    Affiliations Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia, Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Anastasia Bezrukova,

    Roles Formal analysis, Investigation

    Affiliation Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia

  • Katerina Basharova,

    Roles Investigation, Visualization

    Affiliation Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia

  • Sergey Landa,

    Roles Investigation

    Affiliations Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia, Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Zoia Korobova,

    Roles Investigation

    Affiliation Saint Petersburg Pasteur Institute, Saint-Petersburg, Russia

  • Natalia Liubimova,

    Roles Investigation

    Affiliation Saint Petersburg Pasteur Institute, Saint-Petersburg, Russia

  • Ivan Vlasov,

    Roles Investigation

    Affiliation Institute of Molecular Genetics, National Research Center "Kurchatov Institute", Moscow, Russia

  • Mikhael Nikolaev,

    Roles Investigation

    Affiliations Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia, Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Artem Izyumchenko,

    Roles Investigation

    Affiliation Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia

  • Elena Gavrilova,

    Roles Investigation

    Affiliation Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Irina Shlyk,

    Roles Investigation

    Affiliation Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Elena Chernitskaya,

    Roles Investigation

    Affiliation Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Yurii Kovalchuk,

    Roles Investigation

    Affiliation Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  • Petr Slominsky,

    Roles Writing – review & editing

    Affiliation Institute of Molecular Genetics, National Research Center "Kurchatov Institute", Moscow, Russia

  • Areg Totolian,

    Roles Writing – review & editing

    Affiliation Saint Petersburg Pasteur Institute, Saint-Petersburg, Russia

  • Yurii Polushin,

    Roles Writing – review & editing

    Affiliation Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia

  •  [ ... ],
  • Sofya Pchelina

    Roles Conceptualization, Data curation, Project administration, Supervision, Writing – review & editing

    Affiliations Petersburg Nuclear Physics Institute Named by B.P. Konstantinov of National Research Centre «Kurchatov Institute», Gatchina, Russia, Pavlov First Saint-Petersburg State Medical University, Saint-Petersburg, Russia, Kurchatov Genome Center—PNPI, Saint-Petersburg, Russia

  • [ view all ]
  • [ view less ]

Abstract

Transcriptomic analysis conducted by us previously revealed upregulation of genes involved in low-density lipoprotein particle receptor (LDLR) activity pathway in lethal COVID-19 caused by SARS-CoV-2 virus (severe acute respiratory syndrome coronavirus 2). Last data suggested the possible role of extracellular vesicles in COVID-19 pathogenesis. The aim of the present study was to retrospectively evaluate parameters of cholesterol metabolism and newly identified EVs, exomeres, as possible predictors of fatal outcome of COVID-19 patients infected by the Alpha and the Delta variants of SARS-CoV-2 virus. Blood from 67 patients with severe COVID-19 were collected at the time of admission to the intensive care unit (ICU) and 7 days after admission to the ICU. After 30 days patients were divided into two subgroups according to outcome—34 non-survivors and 33 survivors. This study demonstrated that plasma low- and high-density lipoprotein cholesterol levels (LDL-C and HDL-C) were decreased in non-survivors compared to controls at the time of admission to the ICU. The conjoint fraction of exomeres and LDL particles measured by dynamic light scattering (DLS) was decreased in non-survivors infected by the Alpha and the Delta variants compared to survivors at the time of admission to the ICU. We first showed that reduction of exomeres fraction may be critical in fatal outcome of COVID-19.

Introduction

Coronavirus infection (COVID-19) is caused by SARS-CoV-2 virus (severe acute respiratory syndrome coronavirus 2) and significant number of patients with COVID-19 develops critical illness with fatal disease outcome [1]. Despite of the numerous studies devoted to the study of the pathogenesis of COVID-19, there are still many questions, in particular, what factors may determine the lethal outcome. Last data revealed the link between cholesterol metabolism and susceptibility and severity of COVID-19 [2, 3]. In particular, serum lipid levels decrease in patients with COVID-19 and low level of serum low-density lipoprotein cholesterol (LDL-C) have been proposed to be a predictor of poor disease prognosis [2, 4]. In our previous study transcriptome analysis of peripheral blood mononuclear cells (PBMCs) of patients with COVID-19 upon admission to the ICU revealed an activation of the low-density lipoprotein particle receptor (LDLR) activity pathway (GO:0005041) in deceased patients when compared with survivors [5]. Our results point out this hyperactivation of LDLR activity pathway could be linked to cholesterol metabolism disturbances expressed in changes of lipid spectrum observed in another studies. According to the latest point of view extracellular vesicles (EVs) could be considered as alternative cholesterol vehicles in blood plasma [6, 7]. We and others demonstrated the variety in the size of EVs in non-survivors compared to patients with positive outcome, that is in accordance with other studies showing the heterogeneous nature of EVs from COVID-19 patients [810]. EVs include large microvesicles (200–1000 nm), small exosomes (50–200 nm), and newly identified exomeres (<50 nm) that was shown to be very different from exosomes in their lipid composition [11, 12]. We suppose that exomeres may play a role in COVID-19 pathogenesis as an essential role of exomeres in cholesterol transport was shown [13]. To further analyze the role of cholesterol metabolism in poor prognosis of COVID-19 we evaluate retrospectively the parameters of cholesterol metabolism (plasma lipids concentrations and exomere fraction, the expression levels of LDLR activity pathway genes selected through our previous transcriptome analysis) as well as cytokine profile as possible predictors of fatal outcome in patients with severe COVID-19 infected by the Alpha and the Delta variants of SARS-CoV-2 virus.

Materials and methods

Ethics statement

The study involving human participants was approved by the Ethics Committee of the Pavlov First State Medical University of St. Petersburg (Russia). A formal written consent form was provided to all included subjects to read and sign prior to the study at the time of admission to intensive care (ICU).

Subjects

The current study included 67 patients with severe COVID-19 infected with the Alpha and the Delta variants of SARS-CoV-2 who were admitted to the to the intensive care unit (ICU) of Pavlov First Saint-Petersburg State Medical University during the Alpha variant predominant period (between 1st November 2020 to 25th February 2021) and the Delta variant predominant period (between 1st July 2021 and 31th August 2021), respectively. Biological samples (plasma and PBMCs) were collected at the time to admission to the ICU and one week (7 days) after admission to the ICU. Depending on the disease outcome during the time of observation (30 days) patients were divided into subgroups depending of outcome: 33 survivors and 34 non-survivors. Also, all 67 COVID-10 patients were separated in accordance by variant of COVID-19: 28 COVID-19 patients infected by the Alpha variant of SARS-CoV-2: 13 non-survivors and 15 survivors and 39 COVID-19 patients infected by the Delta variant of SARS-CoV-2: 21 non-survivors and 18 survivors. The inclusion criteria were Russian ethnicity, age between 40 and 80 years with the absence of chronic comorbidities such as cancer, cerebrovascular diseases, heart failure, or renal failure but we did not exclude obesity or hypertension, as the common features in patients among patients with severed COVID-19. All patients admitted to the intensive care unit and met the inclusion criteria were included in the study. No one received the lipid-lowering therapy before and during the study. The control group included 40 healthy individuals who had never been infected with SARS-CoV-2 that was additionally confirmed by negative PCR and IgG antibody tests. The characteristics of studies groups are presented in Table 1.

Blood plasma collection

Blood plasma was isolated from peripheral venous blood collected in vacuum containers with EDTA by centrifugation for 20 min, 3000 g. Centrifugation was carried out no later than 20 min after blood taking. Plasma samples were stored at -80°C until the study.

Plasma lipid profile

Total cholesterol (TC), high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), triglycerides (TG) in plasma were determined by OLYMPUS AU400 apparatus.

Plasma cytokine profile

The concentration of cytokines in blood plasma was assessed by multiplex analysis using a MagPix analyzer (Luminex Corporation, USA) using a kit for simultaneous detection of 25 analytes (IL-17F, GM-CSF, IFNγ, IL-10, CCL20/MIP3ɑ, IL-12(p70), IL-13, IL-15, IL-17A, IL-22, IL-9, IL-1β, IL-33, IL-2, IL-21, IL-4, IL-23, IL -5, Il-6, IL-17E/IL-25, IL-27, IL-31, TNFɑ, TNFβ, IL-28A) (Milliplex MAP Human Th17 Magnetic Premix 25 Plex Kit (HT17MG-14K-PX-25, Merk-Millipore, USA) according to the manufacturer’s conditions. The results were processed using xPONENT 4.2 software.

The method of dynamic light scattering (DLS) for detection of exomeres

Relative levels of exomeres and LDL particles were assessed in plasma by means of a combination of immunoadsorption and dynamic light scattering (DLS) methods as described earlier [13, 14]. This approach allows to distinguish particles by size (using hydrodynamic radius) including after preliminary depletion of unnecessary particles by immunoprecipitation with appropriate antibodies to clean particles of interest. Plasma was processed separately with CD9, HSP90 and ApoB-100 anti-bodies with subsequent separation of corresponding particles via centrifugation of 1,000 g for 10 min. Supernatants and initial plasma were analyzed by DLS. Principle of the method is represented in S1 Fig. Conjoint fraction of exomeres and LDLs particles is determined by analysis of corresponding peak from measurement of unprocessed plasma (S1A Fig). Relative level of LDL particles is determined as a mean from analysis of resulting peak after CD9+ and HSP90+ particles (exosomes and exomeres) were removed (S1B and S1C Fig). Relative level of exomeres was determined by analysis of resulting peak after ApoB-100 containing particles were removed (S1D Fig). The measurements were carried out using a laser correlation spectrometer DLS (INTOX MED LLC, St. Petersburg, Russia) with a heterogeneous measurement scheme [15]. Mathematical processing of the obtained data was carried out using the algorithm [15] using the QELSspec (version 3.4) software package, Gatchina, Russia.

Isolation of peripheral blood mononuclear cells

PBMCs were isolated from 8 mL EDTA-anticoagulated venous blood by a Ficoll-Paque gradient method (Ficoll-Paque PLUS, GE Healthcare, Chicago, IL, USA) [16]. After centrifugation, PBMCs were collected from the interface and washed twice with PBS (pH 7.4) to remove the platelet-rich plasma fraction. The PBMCs cell pellets were aliquoted and immediately frozen at −80°C.

RNA extraction, cDNA synthesis and quantitative real-time PCR in PBMCs

Total RNA was extracted from PBMCs using GeneJET RNA Purification Kit (K0731, Thermo Fisher Scientific, USA) and complementary DNA (cDNA) was synthesized by Revert Aid First cDNA Synthesis kit (K1622, Thermo Fisher Scientific, USA. Primers and probes were designed by Primer3 program (https://bioinfo.ut.ee/primer3-0.4.0/) (S1 Table). The expression level of cholesterol metabolism genes (LDLR, PPARG, LRP6, CD36, ANXA2, STAB1) was assessed with TaqMan probes by real-time quantitative PCR on CFX96 (Bio-Rad, USA). mRNA levels of the studied genes were normalized to mRNA levels of RPLP0 (Ribosomal Protein Lateral Stalk Subunit P0) and ACTB (beta actin). Each sample was tested in triplicate. The value of the expression level relative to the calibrator was determined by the formula: 2- ΔΔC [17].

Droplet digital PCR (ddPCR) assay

Nucleic acid copy counts of STAB1 were determined by means of the RPLP0 gene as housekeeping droplets using droplet digital PCR technology (Bio-Rad, USA) following the manufacturer´s protocol and using ddPCR EvaGreen Supermix (Bio-Rad, USA). The thermal cycling conditions for STAB1 was standard as describe in the manufacturer´s protocol of ddPCR EvaGreen Supermix (Bio-Rad, USA) with annealing temperature 57°C. The ddPCR reaction with primers listed above was performed in a T100 Thermal Cycler. Measurement of positive droplets per μl sample was performed on a QX200 ddPCR Droplet Reader (Bio-Rad, USA). Based on the droplet counts and according to Poisson distribution, absolute nucleic acid copy counts were calculated utilizing the software QuantaSoft (Bio-Rad, USA).

Statistical analysis

Conformity of findings to normal distribution was tested using the Shapiro–Wilk test. Comparison of variation between the studied groups was performed by Kruskal-Wallace test. Comparison of variation between two groups was conducted using the nonparametric Mann-Whitney u-test. Correlation analysis was conducted with Spearman coefficient. Linear regression analysis was used to test for association between outcome of COVID-19 and studied parameters, adjusting for age, sex, obesity and hypertension. Significance was established at p<0.05. Statistical analysis was performed using R software (version 4.1.2). Clinical data expressed as ± standard deviation (SD), experimental data as the median (min–max).

Results

This study included the analysis of plasma lipid profile measured by standard technics and DLS. All patients admitted to the ICU were characterized with respiratory failure and lung injury with CT-SKAN scores equal to 3–4. The demographic, clinical and laboratory characteristics of the studied populations are shown in Table 1. COVID-19 patients were characterized with lymphopenia and the elevated values of neutrophils, C-reactive protein (CRP) and D-dimer that is in accordance with previous data [18]. There were no differences between survivors and non-survivors except more pronounced lymphopenia in non-survivors (Table 1). Also, patients with COVID-19 were characterized by increased secretion of proinflammatory cytokines CCL20/MIP3ɑ, IL-10, IL-15, IL-27 in non-survivors than in controls (S1 File).

Plasma lipid profile for all patients with various outcomes of COVID-19 (survivors / non-survivors) at the time and 7 days after admission to the ICU and in the control group

Lipid profile was estimated in plasma of all patients with COVID-19 with different outcomes of COVID-19 (survivors / non-survivors) included in the current study at the time and 7 days after admission to the ICU and in the control group by standard technics (Table 2).

thumbnail
Table 2. Plasma lipid profile for patients with COVID-19 and controls.

https://doi.org/10.1371/journal.pone.0278083.t002

At the time of admission to the ICU we showed decreased TC and LDL levels in non-survivors compared to survivors and control group (TC: p = 0.014, p<0.0001; LDL: p = 0.011, p<0.0001; respectively) and decreased TC and LDL levels in survivors compared to the control group (p = 0.012, p = 0.028, respectively). HDL level was decreased in non-survivors and survivors when comparing to the control group (p<0.0001, p = 0.00064, respectively). No statistically significant differences in TG concentration between studied groups were found (p>0.05) at the time of admission to the ICU. 7 days after of admission to the ICU, TC, LDL and HDL levels were decreased in non-survivors compared to survivors and controls (TC: p = 0.00034, p<0.0001; LDL: p = 0.0002, p<0.0001; HDL: p = 0.018, p<0.0001; respectively). At the same time decreased level of HDL was observed in survivors compared to the control group 7 days after of admission to the ICU (p = 0.0018). TG level was increased in non-survivors and survivors compared to the control group after 7 days of admission to the ICU (p = 0.0039, p = 0.00086, respectively). TC and TG levels were increased in survivors (p = 0.018, p = 0.0054, respectively) and TG in non-survivors at 7 days after of admission to the ICU compared to the time of admission to the ICU (p = 0.042).

Also, plasma lipid profiles were compared in patients with COVID-19 depending on the variant of the SARS-CoV-2 virus (the Alpha, the Delta) (Table 2). TC, LDL and HDL cholesterol in the plasma were decreased in non-survived patients with COVID-19 infected by the Alpha variant compared to controls (p = 0.00028, p = 0.002, p = 0.00023. respectively) and TC and LDL had a tendency to decrease compared with survived patients with COVID-19 (p = 0.054, p = 0.087, respectively) at the time to admission to the ICU. Also, we demonstrated a decrease of plasma HDL cholesterol in survived patients with COVID-19 when compared with controls (p = 0.0038). No statistically significant differences in TG concentration between studied groups were found (p>0.05). In order to assess the dynamic change in plasma lipid profile of survivors and non-survivors during severe COVID-19 we additionally evaluated lipid levels in 7 days after admission to the ICU (Table 2). 7 days after admission to the ICU TC and LDL concentration were decreased in non-survivors compared with controls (p = 0.025, 0.011, respectively) and survivors (p = 0.018, p = 0.011, respectively). HDL concentration was decreased and TG concentration had tendency to increase in non-survivors and survivors compared with controls (p = 0.00054, p = 0.071; p = 0.00082, p = 0.061, respectively). A pairwise analysis comparing plasma lipid profile change in non-survivors and survivors between two observation points (the time and 7 days after admission to the ICU) did not reveal significant differences(p>0.05).

With regard to changes in the lipid profile in patients infected by the Delta variant, we also showed a decreased level of LDL and HDL cholesterol in the plasma of non-survived patients with COVID-19 when compared with controls (p = 0.0013, p = 0.018, p = 0.0021, respectively) and tendency to decrease of LDL, HDL and TC levels in survived patients with COVID-19 (p = 0.062, p = 0.067, p = 0.063, respectively). No statistically significant differences in TG concentration between studied groups were found (p>0.05).

7 days after admission HDL and LDL cholesterol as well as TC levels were decreased in non-survivors compared with survivors (p = 0.0043, p = 0.0067, p = 0.0035 respectively) and the control group. (p = 0.001, p = 0.0011, p = 0.00088, respectively) (Table 2). A pairwise analysis comparing plasma lipid profile change in non-survivors and survivors between two observation points was performed. TC concentration was increased in survivors in 7 days after admission to the ICU when compared with the time of admission (p = 0.037). At the same time non-survivors were characterized by continuing decrease of HDL and LDL cholesterol levels in 7 days after admission to the ICU when compared with the time of admission (p = 0.049, p = 0.0014, respectively). Plasma TG concentration did not differ between the studied groups neither at the time of admission nor 7 days after (p>0.05).

Overall, deceased patients with COVID-19, regardless of the variant with which they were infected (the Alpha, the Delta), were characterized by decreased LDL level at the time of admission. However, the Delta variant was associated with continuous decrease in LDL and HDL levels in deceased patients 7 days after admission in ICU. Thus, more pronounce cholesterol metabolism disturbances in patients affected with the Delta variant could be supposed. Therefore, our further study was focused on the evaluation of LDL-related parameters (exomeres, analysis of gene expression involved in LDLR activity pathway revealed in our previous study) in patients with COVID-19 with the Delta variant [5].

Relative levels of plasma exomeres and LDL particles detected by means of DLS for patients with various outcomes of COVID-19 (survivors / non-survivors) at the time and 7 days after admission to the ICU and in the control group

Exomeres are a class of extracellular vesicles enriched in cholesterol and comparable in size with LDLs. Our approach allowed to distinguish particles by size (using hydrodynamic radius). Relative levels of exomeres and LDLs as well as conjoint fraction of exomeres and LDL particles were assessed by combination of immunoadsorption and DLS. We analyzed conjoint level of particles with hydrodynamic radius of 20 nm (exomers and LDLs), next we analyzed the same peak after immunoadsorption of exomeres using antibodies against CD9 and HSP90 (LDL particles remain after cleaning) and antibodies against apolipoproteins B100 (exomeres remain after cleaning). It should be noted that relative LDLs level determined by this method as relatively real quantity of particles is highly positively correlated with those measured by standard technique.

Level of plasma exomeres was decreased in non-survivors infected by the Delta variant compared with survivors at time to admission to the ICU (p = 0.011), still there were not any differences between the studied groups in 7 days after admission (Fig 1, S2 Table). Level of LDL particles was decreased in non-survivors at the time of admission to the UCI compared with controls (p = 0.039). 7 days after to admission to the ICU level of LDL particles was lower in non-survivors compared with survivors (p = 0.0049) and controls (p = 0.0095). Interestingly, if one takes into account conjoint fraction of exomeres and LDL particles (ExoM_LDL) significant reduction of this parameter in non-survivors when compared with survivors and the control group is traced both at the time of admission to the ICU (p<0.0001, p<0.0001, respectively) and 7 days after (p = 0.00043, p = 0.0089, respectively). Similar changes in the level of ExoM_LDL were observed on the first day of admission to the ICU but not in 7 days in the group of COVID-19 patients infected with the Alpha variant (S2b Fig, S2 Table).

thumbnail
Fig 1.

Relative level of plasma exomeres in patients with the Delta variant and controls with concentration of TC bound by LDL: a. ExoM (admission to the ICU), b. ExoM_LDL (admission to the ICU), c. LDL (admission to the ICU), d. ExoM (7 days after admission to the ICU), e. ExoM_LDL (7 days after admission to the ICU), f. LDL (7 days after admission to the ICU).

https://doi.org/10.1371/journal.pone.0278083.g001

Survival predictive value of studied parameters by ROC analysis

ROC analysis was conducted to determine the threshold value of lipid concentrations (TC, LDL, HDL) measured by standard technique between survived and non-survived all patients with COVID-19, patients infected by the Alpha and the Delta variants separately. Significant threshold value was not found (data are not presented).

Also, ROC was performed to determine the threshold value of relative level of exomeres, LDL particles and conjoint fraction of exomeres and LDLs measured by DLS in patients with the Delta variants. Among all studied parameters reflecting lipid spectrum conjoint fraction of exomeres and LDL particles was the best to distinguish non-survivors and survivors at two estimated points (the time and 7 days after to admission to the ICU) (Fig 2). We have identified a threshold value for ExoM_LDL at the time of admission to the ICU as 0.272 (AUC = 0.903 (CI95%:0.793–1.00), specificity = 0.809, sensitivity = 0.889, accuracy = 0.846, p<0.0001) (Fig 2A). The threshold value for ExoM_LDL 7 days after the time of admission to the ICU was 0.284 (AUC = 0.909 (CI 95%: 0.773–1.000), specificity = 0.75, sensitivity = 1.00, accuracy = 0.80, p = 0.0007) (Fig 2B). Similar results have been demonstrated for patients with the Alpha variant. We have identified a threshold value for ExoM_LDL at the time of admission to the ICU as 0.296 (AUC = 0.796 (CI95%:0.616–0.976), specificity = 0.692, sensitivity = 0.800, accuracy = 0.750, p = 0.011) (S3H Fig).

thumbnail
Fig 2.

ROC analysis for analytes in blood plasma of patients with the Delta variant: A. ExoM_LDL (admission to the ICU), B. ExoM_LDL (7 days after admission to the ICU).

https://doi.org/10.1371/journal.pone.0278083.g002

Correlation analysis between lipid profile measured by DLS and standard technics

The matrices of correlation pairs of lipid profile and the ratio of the contributions to scattering parameters for exomeres (ExoM, ExoM_LDL, LDL) in the three above-mentioned groups (non-survivors with COVID-19, survivors with COVID-19, controls) in both time points for all COVID-19 patients and separately for COVID-19 patients infected by the Alpha and the Delta variants are presented in S3 and S4 Tables. We revealed a positive correlation between LDL measured by DLS and standard technique, namely, between LDL (rplU) with plasma LDL in survivors and non-survivors in both time points (at time to admission to the ICU and 7 days after admission to the ICU for all COVID-19 patients and separately for COVID-19 patients infected by the Alpha and the Delta variants (S3 and S4 Tables). Interesting, fraction ExoM_LDL was strong positive correlated with TC and LDL in survivors and non-survivors in all COVID-19 patients and separately for COVID-19 patients infected by the Alpha variant and only in survivors in group of COVID-19 patients infected by the Delta variant at time to admission to the ICU and in non-survivors in entire cohort of COVID-19 patients and separately for COVID-19 patients infected by the Alpha variant 7 days after to the ICU) (S3 and S4 Tables).

Prediction of fatal outcome of COVID-19 using multivariable regression analysis

Multivariable regression analysis was performed to identify contribution of the parameters of lipid profile and exomeres to fatal outcome of all COVID-19 patients, patients infected by the Alpha and the Delta separately adjusted for hypertension, obesity, age and sex. We found the association between the higher relative level of conjoint fraction of exomeres and LDL particles (ExoM_LDL) and lower risk of fatal outcome of COVID-19 at the time of admission to the ICU (OR = 1.574e+21, CI95%:1.515e+09–1.635e+33, p = 0.0005; OR = 577.776, CI95%:4.473–74621.069, p = 0.0185, OR = 1.075, CI95%:1.032–1.120, p = 0.0015, respectively) (S5 Table).

Validation of differential gene expression received in our transcriptomic data in PBMCs in the studied patients with COVID-19 with different outcome (survivor, non-survivor) infected by the Delta variant of SARS-CoV-2 at the time and 7 days after admission to the ICU and in the control group

We evaluated the expression levels of genes involved in cholesterol metabolism (PPARG, LDLR, LRP6, CD36, STAB1, ANXA2) in PBMCs of patients with various outcomes of COVID-19 (convalescence / death) infected by the Delta variant as variant associated with pronounce decrease of LDL level in non-survivors over time, as well as in controls (S6 Table). Studied genes were selected through our previous transcriptomic analysis that revealed upregulation of genes involved in LDLR activity pathway. At the time of admission to the ICU, increased PPARG and CD36 gene expression levels were found both in non-survivors and survivors when compared with controls (PPARG: p = 0.00024, p = 0.00048, respectively; CD36: p = 0.0027, p<0.0001, respectively), as well as ANXA2 gene expression level was higher only in survivors (p = 0.011) (S6 Table). After 7 days in the ICU increased STAB1, CD36, LDLR and simultaneously decreased ANXA2 gene expression were detected in both non-survivors and survivors compared with the control group (p<0.05). Also, decreased expression levels of PPARG and ANXA2 after 7 days treatment compared with the time of admission were revealed both in non-survivors and survivors (PPARG: p = 0.016, p = 0.016, respectively; ANXA2: p = 0.00098, p = 0.0029, respectively). LRP6 mRNA level was decreased in survivors when compared with the control group in 7 days after admission to the ICU (p = 0.0089). Additionally in survivors LRP6 decreased 7 days after admission to the ICU when compared with the time to admission the ICU (p = 0.02).

It should be noted that only for the STAB1 gene differential expression between non-survivors and survivors was demonstrated (p = 0.017) and this was specific only for 7 days after admission to the ICU. In non-survivors STAB1 expression level was increased 7 days after admission to the ICU compared with the time of the admission (p = 0.0078). ddPCR was used for confirmation of this result. mRNA copy counts were determined for the STAB1 gene in non-survived patients with COVID-19, survived patients with COVID-19 in 7 days after admission to the ICU and in control individuals. Increased expression of STAB1 of non-survivors compared with survivors and controls was revealed (p = 0.024, p<0.0001, respectively). So, the result obtained by RT-PCR and demonstrating distinguishing differences between expression level of the STAB1 gene in non-survivors and survivors was supported by ddPCR (Fig 3). However, ddPCR did not reveal the differences in expression level of the STAB1 gene between controls and survived patients with COVID-19 (p>0.05) in 7 days after admission to the ICU.

thumbnail
Fig 3. STAB1 expression level in PBMCs of COVID-19 patients infected by the Delta variant with different outcomes and controls.

https://doi.org/10.1371/journal.pone.0278083.g003

Discussion

Here, we for the first time carried out simultaneous assessment of the complex of parameters that determine cholesterol metabolism in order to find relevant biomarkers of poor prognosis and fatal outcome of severe COVID-19. Plasma LDL and HDL levels were decreased in non-survivors infected by the Alpha or the Delta variants compared to controls at the time of admission to the ICU. However, the Delta variant has been associated with a more pronounce decrease in LDL levels in non-surviving patients during the time spent in the ICU. The conjoint fraction of exomeres and LDL particles measured by dynamic light scattering (DLS) was decreased in non-survivors infected by the Alpha and the Delta variants compared to survivors and controls at the time of admission to the ICU that may be considered as potential predictor of fatal outcome.

Over the past two years, link between cholesterol metabolism and severity of COVID-19 was shown. As, it has also been demonstrated that COVID-19 is accompanied by dyslipidemia [2, 3]. In particular, serum lipid levels decrease after the acute onset of COVID-19 and continue to decrease until the patient recovers [3]. Serum LDL levels have been found to be a predictor of poor disease prognosis of COVID-19 [4, 19]. Patients with COVID-19 with a fatal outcome were characterized by irreversible decrease of LDL reduced by up to 60% compared with the level at hospitalization [19]. However, this study operated on a small cohort of patients with COVID-19 which included 21 patients with COVID-19 with concomitant chronic pathologies, among them only 4 cases of death. Aparisi and colleagues also demonstrated the association of low LDL serum level with 30-days all-cause mortality in patients with COVID-19 [20]. This study included representative sample number (654 patients, 149 cases of death), but expected the patient group was heterogeneous with apparent prevalence of diabetes mellitus, dyslipidemia, as well as other comorbidities such as chronic kidney disease, ischemic heart disease, which were maximally excluded in our study. Our data generally confirmed these findings: LDL level was decreased 7 days after admission to the ICU in non-survivors compared to survivors independent of SARS-CoV-2 variants. It should be noting that proinflammatory cytokines induced by viral infection may modulate lipid metabolism [21]. Our study found increased secretion of pro-inflammatory cytokines such as CCL20/MIP3ɑ, IL-10, IL-15, IL-27 in non-survivors than in controls that is consistent with previously published data (Supplementary materials) [22]. Also, our finding is that increased expression levels of LDL receptors genes (LDLR, CD36, STAB1) were found in non-survivors and survivors compared to controls at time to admission to the ICU. LDLR encodes canonical LDL receptor, CD36 and STAB1 encode important scavenger receptors for modified, primarily oxidized LDL (oxLDL). It was shown earlier that COVID-19 severity is associated with oxidative stress and products of lipid peroxidation accumulate during disease [23]. Thus, elevated CD36 and STAB1 gene expression may indicate increased oxidized lipoproteins uptake during severe COVID-19 and this in turn partially could explain lipid profile abnormalities observed in our study as well as in others. Survived patients also have increased STAB1 gene expression compared with controls but significantly lower than in non-survivors 7 day after that supported the results of our previous transcriptome analysis [5]. In particular plasma TC, LDL and HDL cholesterol levels were markedly reduced in non-survivors when compared with the control group at the time of admission to the ICU and 7 days after. And while lipid concentrations tended to go up in survivors during treatment in the ICU they were continuously reduced in non-survivors in case of the Delta variant of SARS-CoV-2.

Increased activity of LDL receptors genes in our study reinforces data on enhanced lipid accumulation within cellular lipid droplets observed in SARS-CoV-2 infected cells, as well as in the lungs of deceased patients with COVID-19 [24]. Cholesterol is essential for coronavirus entry, membrane fusion, and pathological syncytia formation [25]. Interaction with ACE2 is facilitated by binding cholesterol to receptor-binding domain (RBD) of SARS-CoV-2 spike protein [26]. At the same time mutations in the RBD were linked to stability of S protein, increased ACE2 binding affinity and contagiousness of the Delta variant [27] therefore this variant could more actively influence cholesterol metabolism and could lead to more pronounced decrease of lipid concentrations.

Additionally, to lipid profile we analyzed cholesterol containing particles in plasma primarily LDLs and extracellular vesicles by alternative method based on immunoadsorption combined with DLS. Our previous study linked production of exomeres to cellular cholesterol synthesis [13]. We and others believe now that exomeres are responsible for essential part of plasma cholesterol [13, 28]. Here we first assessed the сonjoint level of particles with hydrodynamic radius of 20 nm (exomers and LDLs) as well as exomeres and LDLs fractions separately. Relative LDL level determined using this approach generally repeated LDL profiles in the studied groups.

At the same time, we first observed marked decrease in total proportion of exomeres and LDL particles in non-survivors independently of the Alpha or the Delta variants when compared with both survivors and the control group even at the time of admission to the ICU. The role of conjoint fraction of exomeres and LDL particles as predictor of fatal outcome of COVID-19 was also supported via multivariable regression analysis adjusted for hypertension, obesity, age and sex. Various opinions exist concerning specificity of decreased LDL level and COVID19 outcome. Different points of view suggests that a significant decrease in LDL cholesterol level is also a feature of mortality during sepsis [29] or, alternatively, ICU patients receiving membrane oxygenation but not lethality [30]. According to our data conjoin fraction of exomeres and LDL particles may be more accurate than simple estimation of plasma LDL cholesterol concentration at the time of admission to the ICU as a relevant predictor of fatal outcome. However, the nature of the decrease of these classes of particles quantity is unknown. Previous study showed that exomeres are enriched in ecto-domain fragments of ACE2 receptor and thus can bind to SARS-CoV-2 S protein S1 subunit [31]. Possible application of EVs and exomeres containing fully post-translationally modified ACE2 that can be used for an attenuation of SARS-CoV-2 infection ius being discussed [6, 31, 32]. The exact role of exomeres in SARS-CoV-2 infection needs further investigations.

Conclusion

The current study confirmed that the severity of changes in lipid profile may correlate with the outcome of COVID-19. Overall, deceased patients with COVID-19, regardless of the variant with which they were infected (the Alpha, the Delta), were characterized by decreased LDL level at the time of admission. Decreased conjoint fraction of exomeres and LDL particles was identified as predictor of a fatal outcome of COVID-19. The Delta variant was associated with negative dynamics of HDL and LDL plasma levels during the disease and upregulation of LDLR pathway genes. Further studies are needed to analyze if this observation could be referred to other critical illnesses or is specific to COVID-19 infection.

Supporting information

S1 Fig.

Principle of particle count estimation by dynamic light scattering (DLS) combined with immunoadsorption: A. Analysis of unprocessed plasma, B. Analysis of plasma after CD9 positive particles (exosomes, exomeres) depletion, C. Analysis of plasma after HSP90 positive particles (exomeres) depletion, D. Analysis of plasma after apolipoprotein B100 positive particles (LDLs, VLDLs) depletion. Abbreviations on the figure: ApoB100—apolipoprotein B100, Dh–hydrodynamic diameter, HSA–human serum albumin, LDL–low density lipoproteins, VLDL–very low-density lipoproteins.

https://doi.org/10.1371/journal.pone.0278083.s001

(TIF)

S2 Fig.

Relative level of plasma exomeres with concentration of TC bound by LDL for patients with COVID-19 infected by Alpha variant: A. ExoM (admission to the ICU), B. ExoM_LDL (admission to the ICU), C. LDL (admission to the ICU), D. ExoM (7 days after admission to the ICU), E. ExoM_LDL (7 days after admission to the ICU), F. LDL (7 days after admission to the ICU).

https://doi.org/10.1371/journal.pone.0278083.s002

(TIF)

S3 Fig.

ROC analysis for in blood plasma for all patients with COVID-19 (A-F): A. ExoM (admission to the ICU), B. ExoM_LDL (admission to the ICU), C. LDL (admission to the ICU), D. ExoM (7 days after admission to the ICU), E. ExoM_LDL (7 days after admission to the ICU), F. LDL (7 days after admission to the ICU); for patients with COVID-19 infected by the Alpha variant (G-L): G. ExoM (admission to the ICU), H. ExoM_LDL (admission to the ICU), I. LDL (admission to the ICU), J. ExoM (7 days after admission to the ICU), K. ExoM_LDL (7 days after admission to the ICU), L. LDL (7 days after admission to the ICU); for patients with COVID-19 infected by the Delta variant (M-Q): M. ExoM (admission to the ICU), N. LDL (admission to the ICU), O. LDL (admission to the ICU), P. ExoM (7 days after admission to the ICU), Q. ExoM_LDL (7 days after admission to the ICU), R. LDL (7 days after admission to the ICU).

https://doi.org/10.1371/journal.pone.0278083.s003

(TIF)

S1 File. Plasma cytokine profile in the studied patients with COVID-19 with different outcome (survivor, non-survivor) infected by the Delta variant of SARS-CoV-2 at the time to admission to the ICU and in the control group.

https://doi.org/10.1371/journal.pone.0278083.s004

(PDF)

S1 Table. Primers and probe sequences for estimation expression levels of the studied genes.

https://doi.org/10.1371/journal.pone.0278083.s005

(DOCX)

S2 Table. Relative level of plasma exomeres for patients with COVID-19 with concentration of TC bound by LDL.

https://doi.org/10.1371/journal.pone.0278083.s006

(DOCX)

S3 Table. Correlation analysis between lipid profile parameters and SIC exomeres in patients with COVID-19 at the time of admission to ICU and in control group.

https://doi.org/10.1371/journal.pone.0278083.s007

(DOCX)

S4 Table. Correlation analysis between lipid profile parameters and SIC exomeres in patients with COVID-19 7 days after admission to the ICU.

https://doi.org/10.1371/journal.pone.0278083.s008

(DOCX)

S5 Table. The association between studied parameters of lipid parameters measured by DLS and standard technics for patients with COVID-19.

https://doi.org/10.1371/journal.pone.0278083.s009

(DOCX)

S6 Table. mRNA levels of the studied genes in PBMCs of patients with severe COVID-19 with Delta variant and in the control group.

https://doi.org/10.1371/journal.pone.0278083.s010

(DOCX)

References

  1. 1. Berekaa MM. Insights into the COVID-19 pandemic: origin, pathogenesis, diagnosis, and therapeutic interventions. Front Biosci (Elite Ed). 2021;13: 117–139. pmid:33048778
  2. 2. Wei X, Zeng W, Su J, Wan H, Yu X, Cao X, et al. Hypolipidemia is associated with the severity of COVID-19. J Clin Lipidol. 2020;14: 297–304. pmid:32430154
  3. 3. Sorokin AV, Karathanasis SK, Yang ZH, Freeman L, Kotani K, Remaley AT. COVID-19-associated dyslipidemia: implications for mechanism of impaired resolution and novel therapeutic approaches. FASEB J Off Publ Fed Am Soc Exp Biol. 2020;34: 9843–9853. pmid:32588493
  4. 4. Hu X, Chen D, Wu L, He G, Ye W. Declined serum high density lipoprotein cholesterol is associated with the severity of COVID-19 infection. Clin Chim Acta. 2020;510: 105–110. pmid:32653486
  5. 5. Vlasov I, Panteleeva A, Usenko T, Nikolaev M, Izumchenko A, Gavrilova E, et al. Transcriptomic profiles reveal downregulation of low-density lipoprotein particle receptor pathway activity in patients surviving severe COVID-19. Cells. 2021;10: 3495. pmid:34944005
  6. 6. Pfrieger FW, Vitale N. Cholesterol and the journey of extracellular vesicles. J Lipid Res. 2018;59: 2255–2261. pmid:29678958
  7. 7. Wang W, Zhu N, Yan T, Shi YN, Chen J, Zhang CJ. The crosstalk: exosomes and lipid metabolism. Cell Commun Signal. 2020;18: 119. pmid:32746850
  8. 8. Sirotkina OV, Ulitina AS, Kulabukhova DG, Nikolaev MA, Izyumchenko AD, Garaeva LA, et al. Correlation of laboratory markers of hemo-static system activation with concentration and size of plasma extracellular microparticles in patients with COVID-19. The Scientific Notes of the Pavlov University. 2022;29: 28–36. (In Russ.)
  9. 9. Krishnamachary B, Cook C, Spikes L, Chalise P, Dhillon NK. Extracellular vesicle‐mediated endothelial apoptosis and EV‐associated proteins correlate with COVID‐19 disease severity. J Extracell Vesicles. 2021;10: e12117. pmid:34262673
  10. 10. Sun B, Tang N, Peluso MJ, Iyer NS, Torres L, Donatelli JL, et al. Characterization and biomarker analyses of post-COVID-19 complications and neurological manifestations. Cells. 2021;10: 386. pmid:33668514
  11. 11. Anand S, Samuel M, Mathivanan S. Exomeres: a new member of extracellular vesicles family. Subcell Biochem. 2021;97: 89–97. pmid:33779915
  12. 12. Zhang Q, Jeppesen DK, Higginbotham JN, Graves-Deal R, Trinh VQ, Ramirez MA, et al. Supermeres are functional extracellular nanoparticles replete with disease biomarkers and therapeutic targets. Nat Cell Biol. 2021;23: 1240–1254. pmid:34887515
  13. 13. Landa S, Verlov N, Fedorova N, Filatov M, Pantina R, Burdakov V, et al. Extracellular particles as carriers of cholesterol not associated with lipoproteins. Membranes (Basel). 2022;12: 618. pmid:35736324
  14. 14. Landa SB, Korabliov PV, Semenova EV, Filatov MV. Peculiarities of the formation and subsequent removal of the circulating immune complexes from the bloodstream during the process of digestion. F1000Research. 2018;7: 618. pmid:30079242
  15. 15. Lebedev AD, Ivanova MA, Lomakin AV, Noskin VA. Heterodyne quasi-elastic light-scattering instrument for biomedical diagnostics. Appl Opt. 1997;36: 7518–7522. pmid:18264264
  16. 16. Kanof ME, Smith PD, Zola H. Isolation of whole mononuclear cells from peripheral blood and cord blood. Curr Protoc Immunol. 2001;7: 7.1.1–7.1.7. pmid:18432823
  17. 17. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25: 402–408.
  18. 18. Ucciferri C, Caiazzo L, Di Nicola M, Borrelli P, Pontolillo M, Auricchio A, et al. Parameters associated with diagnosis of COVID-19 in emergency department. Immunity, inflammation and disease. 2021;9: 851–861. pmid:33960733
  19. 19. Fan J, Wang H, Ye G, Cao X, Xu X, Tan W, et al. Letter to the Editor: low-density lipoprotein is a potential predictor of poor prognosis in patients with coronavirus disease 2019. Metabolism. 2020;107: 154243. pmid:32320740
  20. 20. Aparisi Á, Iglesias-Echeverría C, Ybarra-Falcón C, Cusácovich I, Uribarri A, García-Gómez M, et al. Low-density lipoprotein cholesterol levels are associated with poor clinical outcomes in COVID-19. Nutr Metab Cardiovasc Dis. 2021;31: 2619–2627. pmid:34353699
  21. 21. Vu CN, Ruiz-Esponda R, Yang E, Chang E, Gillard B, Pownall HJ, et al. Altered relationship of plasma triglycerides to HDL cholesterol in patients with HIV/HAART-associated dyslipidemia: further evidence for a unique form of metabolic syndrome in HIV patients. Metabolism. 2013;62: 1014–1020. pmid:23522788
  22. 22. Monserrat J, Gómez-Lahoz A, Ortega MA, Sanz J, Muñoz B, Arévalo-Serrano J, et al. Role of innate and adaptive cytokines in the survival of COVID-19 patients. International journal of molecular sciences. 2022;23(18): 10344. pmid:36142255
  23. 23. Martín-Fernández M, Aller R, Heredia-Rodríguez M, Gómez-Sánchez E, Martínez-Paz P, Gonzalo-Benito H, et al. Lipid peroxidation as a hallmark of severity in COVID-19 patients. Redox Biol. 2021;48: 102181. pmid:34768063
  24. 24. Nardacci R, Colavita F, Castilletti C, Lapa D, Matusali G, Meschi S, et al. Evidences for Lipid Involvement in SARS-CoV-2 Cytopathogenesis. Cell Death Dis. 2021;12: 263. pmid:33712574
  25. 25. Dai J, Wang H, Liao Y, Tan L, Sun Y, Song C, et al. Coronavirus Infection and Cholesterol Metabolism. Frontiers in immunology. 2022;13: 791267. pmid:35529872
  26. 26. Wei C, Wan L, Yan Q, Wang X, Zhang J, Yang X, et al. HDL-scavenger receptor B type 1 facilitates SARS-CoV-2 entry. Nat Metab. 2020;2: 1391–1400. pmid:33244168
  27. 27. Goher SS, Ali F, Amin M. The Delta Variant Mutations in the Receptor Binding Domain of SARS-CoV-2 Show Enhanced Electrostatic Interactions with the ACE2. Medicine in drug discovery. 2021; 100114. pmid:34901826
  28. 28. Zhang H, Freitas D, Kim HS, Fabijanic K, Li Z, Chen H, et al. Identification of distinct nanoparticles and subsets of extracellular vesicles by asymmetric flow field-flow fractionation. Nat Cell Biol. 2018;20: 332–343. pmid:29459780
  29. 29. Tanaka S, De Tymowski C, Stern J, Bouzid D, Zappella N, Snauwaert A, et al. Relationship between liver dysfunction, lipoprotein concentration and mortality during sepsis. PloS one. 2022;17: e0272352. pmid:35994439
  30. 30. Tanaka S, De Tymowski C, Zappella N, Snauwaert A, Robert T, Lortat-Jacob B, et al. Lipoprotein concentration in patients requiring extracorporeal membrane oxygenation. Scientific reports. 2021;11: 17225. pmid:34446802
  31. 31. Zhang Q, Jeppesen DK, Higginbotham JN, Franklin JL, Crowe JE Jr, Coffey RJ. Angiotensin-converting enzyme 2-containing small extracellular vesicles and exomeres bind the severe acute respiratory syndrome coronavirus 2 spike protein. Gastroenterology. 2021;160: 958–961.e3. pmid:33022277
  32. 32. El-Shennawy L, Hoffmann AD, Dashzeveg NK, McAndrews KM, Mehl PJ, Cornish D, et al. Circulating ACE2-expressing extracellular vesicles block broad strains of SARS-CoV-2. Nat Commun. 2022;13: 405. pmid:35058437