Molecular pathology of acute respiratory distress syndrome, mechanical ventilation and abnormal coagulation in severe COVID-19

Antoine Soulé, William Ma, Katelyn Yixiu Liu, Catherine Allard, Salman Qureshi, Karine Tremblay, Amin Emad, Simon Rousseau
doi: https://doi.org/10.1101/2023.03.09.23286797

Abstract

Systemic inflammation in critically ill patients can lead to serious consequences such as acute respiratory distress syndrome (ARDS), a condition characterized by the presence of lung inflammation, edema, and impaired gas exchange, associated with poor survival. Understanding molecular pathobiology is essential to improve critical care of these patients. To this end, we use multimodal profiles of SARS-CoV-2 infected hospitalized participants to the Biobanque Québécoise de la COVID-19 (BQC-19) to characterize endophenotypes associated with different degrees of disease severity. Proteomic, metabolomic, and genomic characterization supported a role for neutrophil-associated procoagulant activity in severe COVID-19 ARDS that is inversely correlated with sphinghosine-1 phosphate plasma levels. Fibroblast Growth Factor Receptor (FGFR) and SH2-containing transforming protein 4 (SHC4) signaling were identified as molecular features associated with endophenotype 6 (EP6). Mechanical ventilation in EP6 was associated with alterations in lipoprotein metabolism. These findings help define the molecular mechanisms related to specific severe outcomes, that can be used to identify early unfavorable clinical trajectories and treatable traits to improve the survival of critically ill patients.

Introduction

Acute respiratory distress syndrome (ARDS), a condition characterized by the presence of lung inflammation, edema, and impaired gas exchange, is observed in severe coronavirus disease 2019 (COVID-19) and associated with poor survival. Mechanical ventilation is an important critical care management strategy for severe COVID-19 (1). However, during the first waves of the pandemic, almost half of patients undergoing invasive mechanical ventilation died based on case fatality rate reports (2). Defining the molecular mechanisms related to specific severe outcomes, such as ARDS or acute respiratory failure that require invasive mechanical ventilation, is important to identify treatable traits and improve the survival of critically ill patients.

Endophenotypes are subgroups which are inapparent to traditional methods but share a common set of factors that dictate their response in a manner that is distinct from other sub-groups (3). Endophenotypes have been described in airway diseases such as asthma (4), COPD, (5) and chronic rhinosinusitis (6). Recent studies have identified endophenotypes associated with COVID-19 severity using different large data sets, such as bulk transcriptomic from blood (7), bulk transcriptomic from neutrophil-enriched fractions (8), or proteomic in our accompanying manuscript (9).

In this manuscript, we report multimodal characterization of one of the six endophenotypes (EPs) identified in our accompanying manuscript (9), which we found to be associated with the worst clinical trajectory during hospitalization following SARS-CoV-2 infection (including ARDS). This analysis reveals known and novel molecular effectors of immuno-thrombosis, ARDS and mechanical ventilation associated with COVID-19.

Results

Serprocidins are enriched in an endophenotype associated with COVID-19 severity and ARDS

To gain insight into the molecular mechanisms underlying the severity in COVID-19, we investigated the aptamers enriched within an endophenotype (EP6) (Figure 1 and Table S1) that we previously discovered using the circulating proteome of SARS-CoV-2 infected and hospitalized participants to BQC19 (n = 731) and found to be associated with severe outcome and ARDS (9). These aptamers correspond to measurements obtained using a multiplex SOMAmer affinity array. We used two-sided Mann–Whitney U (MWU) tests and the p-values were corrected for the number of aptamers and EPs using Benjamini–Hochberg (BH) false discovery rate (FDR). Consistent with the elevated neutrophil numbers in EP6, the most enriched aptamer was Cathepsin G (FDR = 4.10E-58), a serine protease localized in neutrophil azurophil granules. Notably, two other serprocidins, neutrophil elastase and proteinase-3, were also significantly enriched in EP6 (FDR = 1.04E-10 and FDR = 8.27E-7, respectively, Figure 1 and Table S1). Moreover, another NETosis-associated protein, Olfactomedin 4 (OLFM4) (10), was greatly enriched in EP6 (top 5%; FDR = 2.61E-46, Figure 1 and Table S1). This suggests an important role for neutrophil degranulation and/or NETosis in severe COVID-19.

Figure 1:
Figure 1:

Distribution of aptamer expressions in EP6 versus other EPs.

The violin plots show the distribution of aptamer expressions for patients in EP6 (magenta) versus EP1-5 (teal). Two-sided Mann–Whitney U (MWU) tests followed by Benjamini–Hochberg method for adjusting the p-values were used to calculate false discovery rates.

Pathway enrichment analysis identified FGFR-signaling in severe COVID-19 acute respiratory distress syndrome

To capture a more comprehensive landscape of the signaling networks active in EP6, we performed pathway enrichment analysis using the Knowledge Engine for Genomics (KnowEnG) (11). Since aptamers were used to identify the EPs, it is expected that many of them would have significant associations with the EPs. As a result, we selected a very strict threshold of FDR < 1E-20 (two-sided MWU test, BH FDR) to identify top aptamers associated with EP6 (9). We then used the gene set characterization pipeline of KnowEnG with Reactome pathway collection (12). This analysis showed that EP6 is characterized by pathways associated with interleukins and cytokine signaling in the immune system. In addition, fibroblast growth factor receptor (FGFR) signaling was identified in EP6, suggesting this pathway as a potential driver of severe pathology (Table 1 and S1).

View this table:
Table 1:

Reactome pathways associated with EP6, based on expression of aptamers. KnowEnG analytical platform was used. The p-values were calculated using a one-sided Fisher’s exact test and were corrected for multiple tests using Benjamini–Hochberg method. Only signaling pathways with FDR < 5E-4 are shown in this table (see Table S1 for the full list).

EP6 metabolomic signature was associated with pro-coagulant activities and liver dysfunction

To further characterize the pathobiology of severe COVID-19 disease, metabolomic profiling of plasma samples was performed in parallel to the aptamer analysis. The results yielded data on 1,435 metabolites, of which 576 were found significantly altered in EP6 (two-sided MWU FDR < 0.01) (Table S2). Characterization of the plasma samples supported the distinction among the different EPs at the levels of metabolomic sub-pathways and individual metabolites (Figure 2 and Table S2). The two top sub-pathways (FDR < 0.01) “Methionine, Cysteine, SAM and Taurine Metabolism” and “phosphatidinylethanolamine (PE)”, are known to interact (13, 14) and are associated with pro-coagulant activities. Phospholipids-containing microparticles from platelet activation contribute to Tissue Factor activation and pro-thrombinase activity (15), linking PE to coagulation and D-dimers level, which are elevated in EP6 (9). PE methylation acts as a major consumer of S-Adenosylmethionine (SAM) leading to the synthesis of S-Adenosylhomocysteine (SAH) and cystathionine (14), upstream of 2-hydroxybutyrate and 2-aminobutyrate in a model of hepatoxicity (16). SAH, cystathionine, 2-hydroxybutyrate and 2-aminobutyrate are all significantly enriched in EP6 (Figure 2B and 3A). To identify relationships that may shed light on factors influencing clinical laboratory results that define EP6, we also performed Spearman’s rank correlation analyses between 21 blood variables and metabolites (Table S3). SAH is positively correlated with creatinine, while cystathionine is positively correlated with creatinine and negatively correlated with albumin (Figures 4, S1-S2), congruent with the data from the model of hepatoxicity (16), linking abnormal coagulation and liver damage in severe COVID-19. We also observed kynurenine to be enriched in EP6 but depleted in EP1 (Figures 3 and 5). Kynurenine levels were also found to be positively correlated with D-Dimers levels (Figure 4 and S1-S2), suggesting a role of this metabolite of the tryptophan pathway in abnormal coagulation in severe COVID-19. Interestingly, TNFSF18 (GITR), a T-cell receptor that can activate indoleamine 2,3-dioxygenase (IDO) to increase kynurenine synthesis (17) in allergy, is enriched in EP6 and depleted in EP1 (Figures 3 and 5).

Figure 2:
Figure 2:

Metabolite characteristics of endophenotypes (EPs).

A) The heatmap shows the enrichment (one-sided Fisher’s exact test) of EPs in different metabolite sub-pathways. See Table S2 for the full list. B) The heatmap shows the over-expression (red) and under-expression (blue) of metabolites in different EPs (two-sided Mann– Whitney U test). Row names show metabolites followed by the sub-pathway to which they belong in parentheses. Only top 15 metabolites (based on false discovery rate for EP6) for which a definite name and sub-pathway was available are shown. Full list is provided in Table S2.

Figure 3:
Figure 3:

Distribution of metabolite levels in EP6 versus other EPs.

The violin plots show the distribution of metabolite levels for patients in EP6 (magenta) versus EP1-5 (teal). Two-sided Mann–Whitney U (MWU) tests followed by Benjamini–Hochberg method for adjusting the p-values were used to calculate false discovery rates.

Figure 4:
Figure 4:

Correlation between different variables in EP6.

Each circle in the scatter plots correspond to a patient in EP6. The circles are colored based on the density of points surrounding them for clearer visualization. SCC corresponds to Spearman’s rank correlation coefficient. A) Scatter plots reflecting the correlation between metabolites (y-axis) and blood variables (x-axis). B) Scatter plot corresponding to correlation of 6-oxopiperidine-2-carboxylate and SHC4.

Figure 5:
Figure 5:

Distribution of metabolite levels in EP1 versus other EPs.

The violin plots show the distribution of metabolite levels for patients in EP1 (red) versus EP2-6 (blue). Two-sided Mann–Whitney U (MWU) tests followed by Benjamini–Hochberg method for adjusting the p-values were used to calculate false discovery rates.

Low levels of sphingosine-1 phosphate are associated with COVID-19 ARDS

Strikingly, sphingosine-1 phosphate was enriched in EP1 and depleted in EP6 (Figures 3 and 5), confirming the previously reported negative association of this metabolite with COVID-19 ARDS (18). Accordingly, the aptamer detecting neutral ceramidase, an enzyme converting ceramides into sphingosine, is enriched in EP1 (FDR = 2.65E-10, Table S1). While changes in protein abundance as detected by aptamers may not necessarily reflect changes in enzymatic activity, it is interesting to note that dihydroceramides and ceramides are depleted in EP1 (Figure 5, Table S1). Conversely, dihydroceramides and ceramides are significantly enriched in EP6 (Figure 3, Table S1). Moreover, one of the top aptamers found enriched in EP6 (FDR = 2.10E-56, Figure 1) is serine palmitoyltransferase 2 (SPTLC2), an enzyme that is responsible for the rate-limiting step in de novo sphingolipid biosynthesis (19, 20).

Mechanical ventilation in EP6 is associated with lipoprotein metabolism

As noted in the introduction, mechanical ventilation in COVID-19 is associated with a high mortality rate. Within EP6, 55 out of the 118 individuals received mechanical ventilation. Within all EP6 members that share the overall similar aptamer signature, we identified a subset of aptamers and metabolites that were significantly different between those with or without mechanical ventilation (Figure 6A and Table S4). The second most enriched aptamer (Table 2 and S4), angiopoietin-like 3 (ANGPTL3), corresponds to a liver-secreted regulator of lipoprotein metabolism (21), which may reflect alterations observed in fatty acid metabolites (Table S4). There were no significant differences between the two groups in terms of sex, age, or body mass index (BMI) (Figure 6B-6C and Table S4). Finally, we investigated whether molecular information (proteomic and/or metabolomic) was associated with mechanical ventilation duration. We identified alpha-L-iduronidase (IDUA) to be inversely correlated with duration of ventilation in the patients with mechanical ventilation in EP6 (Spearman’s rank correlation = −0.60, FDR = 2.57E-2, Figure 6D and Table S4). Overall, these results point to an association between altered lipoprotein metabolism and mechanical ventilation in EP6.

View this table:
Table 2:

Aptamers associated with mechanical ventilation in EP6.

In this table, p-values are obtained using a two-sided Mann–Whitney U test and is corrected for multiple tests using Benjamini–Hochberg FDR. A value of +1 (−1) in the “sign” column shows that the value of aptamer was higher (lower) in the samples of 55 patients that received mechanical ventilation in EP6, compared to the rest of EP6 members. Only 14 targets with FDR < 1.2E-02 are shown. The full list is provided in Table S4.

Figure 6:
Figure 6:

Mechanical ventilation in EP6.

A) The violin plots show the distribution of aptamer expressions for patients in EP6 that received mechanical ventilation (MV+) versus those that did not (MV-). Two-sided Mann–Whitney U (MWU) tests followed by Benjamini–Hochberg method for adjusting p-values were used to calculate false discovery rates. B) Distribution of age and BMI of MV+ and MV-patients in EP6. P-values are calculated using two-sided MWU tests. C) The status of sex for MV+ and MV-patients in EP6. P-value is calculated using a Fisher’s exact test. D) The scatter plot shows the correlation between expression of IDUA aptamer and number of days that patients received mechanical ventilation (n = 50).

SHC4 genotype and protein expression levels are associated with higher odds of belonging to EP6

To obtain independent support for the involvement of specific proteins identified using aptamers, we leveraged genome-wide association study (GWAS) genotyping data corresponding to these patients. We identified 25 single nucleotide variations (SNVs) distributed in 13 annotated genes that were below a p-value threshold of 1E-4 and distinguished EP6 from the rest (Table 3). We then investigated each of the SNVs to which we could assign a gene and an aptamer, to assess whether their protein product in circulation was differentially regulated by the genotype (Table 4). We discovered two genes, SHC4 (encoding SHC adaptor protein 4) and CACNA2D3 (encoding calcium voltage-gated channel auxiliary subunit alpha2 delta3) for which there was a significant association between genotype and protein expression levels (p-values < 0.05). While CACNA2D3 may have mild impact on EP6 membership (odds ratio = 0.61, Table 4), SHC4 was one of the top-enriched aptamers (position 32 out of 4,985, Figure 1) for inclusion in EP6, with an odds ratio of 11.98 for the protein product and 2.00 for the alternate allele, respectively. Finally, to gain further insight into the potential role of SHC4 in COVID-19 disease severity, we investigated the metabolites associated with SHC4 in EP6. The only significantly correlated metabolite was 6-oxopiperidine-2-carboxylate (Spearman’s rank correlation = 0.40, FDR = 0.012, Figures 4, S1-S2), where the p-values were corrected for the number of metabolites using BH FDR. Therefore, the GWAS analysis revealed that the signaling adaptor protein SHC4 may play an important mechanistic role in contributing to severe disease pathology.

View this table:
Table 3:

SNVs differentiating EP6 against all other endophenotype clusters.

View this table:
Table 4:

Association between genotypes and aptamer expression levels

Discussion

In this study, we leveraged the molecular features associated with endophenotypes discovered using unsupervised hierarchical clustering of circulating proteome of SARS-CoV-2 positive hospitalized participants to BQC19 to improve our understanding of pathobiology. This approach enabled the identification of the following molecular features: 1) the level of Sphingosine-1 phosphate is highest in EP1 and lowest in EP6, reflecting an inverse association with time to admission to ICU and/or death; 2) interleukins and FGFR signaling are molecular pathways associated with severe COVID-19; 3) genetic and environmental determinants of SHC4 expression are associated with inclusion in EP6; 4) the combined proteomic and metabolomic signature of EP6 supports a pathologic role for neutrophil-associated procoagulant activity in COVID-19 ARDS; and 5) mechanical ventilation in EP6 is associated with altered lipoprotein metabolism.

COVID-19 ARDS molecular pathology: a counterbalance between ceramides and sphingosine is associated with critical illness

The datasets used in this study carry rich molecular information on mechanisms of disease for COVID-19. EP6 is characterized by its enrichment in ARDS (50% vs 9% in EP1), acute kidney injury (50%) and liver dysfunction (22%) (9). Low levels of Sphingosine 1-phosphate (S1P) are associated with ARDS and have been shown to be associated with greater ICU admission and decreased survival in COVID-19 (18), and is supported by this study where S1P levels were found depleted in EP6 but enriched in EP1 (the endophenotype with the best prognostic). The study also supports a shunting away from sphingosine towards more pro-inflammatory ceramides that are associated with metabolic disorders (22), suggesting a counterbalance between ceramides and sphingosine, where the former is associated with poorer outcomes during critical illness, whereas higher levels of the latter is associated with more favorable outcomes in ARDS.

Molecular mechanisms of mechanical ventilation in EP6: a potential role for ANGPTL3

ANGPTL3 is a liver-secreted enzyme that can inhibit lipoprotein and endothelial lipase (23). Loss-of-function of ANGPTL3 is associated with hypolipidemia (24). Therefore, increased circulating levels of ANGPTL3 would be expected to raise plasma levels of triglycerides and phospholipids which is consistent with the observed metabolomic profiles (Table S4). Whether these alterations are sufficient to increase the likelihood of mechanical ventilation remains speculative. Nevertheless, it is interesting to note that administration of medium- and long-chain triglycerides can aggravate gas exchange abnormalities in ARDS (25) and may warrant the preclinical investigation of ANGTPL3 antagonists like evinacumab or other functionally similar therapeutic agents (26) in ameliorating outcomes in ARDS. As for the duration of mechanical ventilation, we discovered an inverse correlation with IDUA, a lysosomal enzyme required for the degradation of heparan sulfate and dermatan sulfate (27). Our current data cannot distinguish between the presence of this enzyme as a biomarker or as an active mediator that influences the course of COVID-19-associated acute respiratory failure.

Multimodal profiling of EP6 supports a role for neutrophil-derived procoagulant activity-mediated organ damage in COVID-19

COVID-19 ARDS has been found to be more frequently associated with thrombotic complications than non-COVID-19 ARDS (28, 29). In addition, histopathologic evidence of pulmonary microthrombi is frequently observed in autopsies of individuals deceased with COVID-19 (30). Early reports suggested a role for neutrophil-mediated release of NETs contributing to endothelial dysfunction as a mechanism of microthrombosis in COVID-19 ARDS (31-42). The neutrophil-derived serprocidins cathepsin G, neutrophil elastase, and proteinase-3, are significantly enriched in EP6. In addition to their roles in bacterial killing, serprocidins couple innate immunity to coagulation, where they promote coagulation by enhancing tissue factor and factor XII-dependent coagulation (43). Serprocidins are also important regulators of NETosis (44). In pediatric sepsis, another NETosis-associated protein olfactomedin 4 (OLFM4) (10), enriched in EP6, is associated with increased odds of having greater organ failure and death (45). The procoagulant environment associated with COVID-19-associated organ damage is further supported by the alteration in metabolic profiles in EP6. Phosphatidinylethanolamines (PE) become exposed at the surface of cell membranes upon exposure to stress, inflammation, and cell death (46, 47). In a Syrian hamster model, infection with SARS-CoV-2 markedly increased circulating PE expression in the animals that were fed a high salt, high fat diet, demonstrating the interaction between infection and metabolic disorder (48). Furthermore, the membranes of azurophilic granules, which contain serprocidins, are enriched in PE (49). Accordingly, EP6 is associated with decreased platelet counts, increased D-Dimers, international normalized ratio (INR) and activated partial thromboplastin time (aPTT) (9), hallmarks of abnormally activated coagulation pathways that are associated with serious and often lethal complication of sepsis (50). Accordingly, we report an association between kynurenine and D-dimer, pointing to an involvement of tryptophan metabolism in abnormal coagulation. Interestingly, TNFSF18, a T-cell receptor that promotes leukocyte adhesion to endothelial cells (51), is enriched in EP6. A SNV found in the proximity of TNFSF18 was found associated with severity of COVID19 in the western Indian population (52). Unfortunately, this SNV was absent from the array used in BQC19 preventing the confirmation of that observation, but it is nevertheless suggestive of a link between TNFSF18, kynurenine, and abnormal coagulation in severe COVID-19. Collectively, these findings point to immuno-thrombosis activity in EP6 in response to SARS-CoV-2 infection that is associated with abnormal coagulation. Recent results from the Antithrombotic Therapy to Ameliorate Complications of COVID-19 (ATTACC) study showed anticoagulation to be efficacious in moderate but not severe COVID-19 disease (53, 54), revealing a need for alternative management approaches in critically ill patients.

Novel molecular markers of COVID-19 pathology: FGFR and SHC4

Two novel molecular factors associated with COVID-19 ARDS that were identified by our study are FGFR and SHC4. Circulating levels of the pro-angiogenic FGF-2 has been associated with COVID-19 severity and creatinine levels in a study of 208 SARS-CoV-2 positive participants (55). It is also noteworthy that the use of nintedanib, an inhibitor of FGFR, vascular endothelial growth factor receptor (VEGFR), and platelet-derived growth factor receptor (PDGF-R) that is approved for use in interstitial lung disease, improved pulmonary inflammation and helped three middle-aged obese COVID-19 patients with markedly impaired lung function wean off mechanical ventilation (56). While the adaptor protein SHC4 has not been experimentally demonstrated to modulate FGFR signaling, a 12-gene biomarker signature associated with melanoma contains FGFR2, FGFR3 and SHC4 (57). It is attractive to speculate that SHC4 may act downstream of FGFR or other associated growth factor receptors, favoring heightened procoagulant activity associated with COVID-19 ARDS. In view of the limited knowledge of this understudied member of the SHC family, we looked at the metabolites associated with SHC4 to gain insights into its possible functions. Interestingly, 6-oxopiperidine-2-carboxylate (the only metabolite with a significant correlation with SHC4 in EP6) was found to be negatively associated with glomerular filtration rate in a GWAS study of kidney disease and hypertension in African Americans (58). Acute kidney injury is a frequent complication of acute liver failure (59), and liver dysfunction is associated with abnormal coagulation (60). Thus, the overall molecular information coming from the multi-modal analysis of EP6 points to a triad between liver function, kidney function, and hemostasis that becomes dysfunctional following ARDS-associated inflammation driven by SARS-CoV-2 infection. Whether the presence of SHC4 in circulation is a marker of dysfunction of this triad, or an active mediator, remains to be determined. Moreover, the identity of the cells expressing SHC4 leading to its presence in the circulation is not known. Taken together, our identification of FGFR and SHC4 signaling pathways distinguishing EP6 from other endophenotypes provides a rationale for investigation of their potential utility as biomarkers of severe disease activity and/or the use of antagonists of those pathways to treat severe COVID-19 disease manifestations.

Methods

Datasets and preprocessing

We obtained data corresponding to clinical, genomics, metabolomics, and proteomics (circulating proteome) of n = 731 hospitalized and SARS-CoV-2 positive patients (based on qRT-PCR) from the Biobanque Québécoise de la COVID-19 (BQC19; www.quebeccovidbiobank.ca). The circulating proteome measurements obtained using a multiplex SOMAmer affinity array (SomaLogic, 4,985 aptamers) between April 2020 and April 2021 were preprocessed and used to identify six endophenotypes (EP1-EP6). The details of the cohort, circulating proteome data processing, and identification of endophenotypes are provided in the accompanying manuscript (9).

In this study, we used the batch-normalized, missing-values-imputed, and log-transformed version of the BQC19 metabolomic data (1,435 metabolites). BQC19 GWAS imputation data were generated by Tomoko Nakanishi at Brent Richards lab, Jewish General Hospital and McGill University. Detailed codes used for generating the data can be found in: https://github.com/richardslab/BQC19_genotype_pipeline.

Analyses of the GWAS dataset

For the GWAS analyses, annotation of SNVs were done using the biomaRt package (61) from R (62), and all analyses were done using R version 4.1.3. Quality control steps were derived in majority from a 2017 QC tutorial article (63). At the beginning, we had 867,450 markers and 2,429 samples. We imported Plink format data into R using the “read_plink” function from genio R package (64). We removed 103,592 non ACGT bi-allelic markers. We calculated the predicted sex by looking at the rate of homozygote markers on chromosome 23. We removed 3,588 markers with call rates < 98%, 448,932 monomorphic markers and markers with MAF < 0.05, and 28,092 markers with Hardy–Weinberg equilibrium < 1E-6 (calculated by the “HWE.exact” function from the genetics package (65) from R. For the EP6 cluster group analyses, we additionally removed 1,747 samples that were not in the cluster analysis, 8 samples with a sex discrepancy (based on predicted sex calculated earlier and reported sex), and 3 samples with a heterozygosity rate > 3 standard deviations. We finally removed a pair of samples which had approximately the same genome, likely due to an error of manipulation. As we could not know which one was correct, we removed both. We also found 2 pairs of individuals with a pi-hat of ∼0.5 (meaning first degree relatives), so we kept the one sample per pair with the higher call rate. All other pairs of individuals had a pi-hat < 0.21 that is judged acceptable considering our population. Pi-hats were calculated with the “snpgdsIBDMoM” function from SNPRelate R package (66). At the end of quality controls, 283,246 markers on 655 samples had been used to perform association analyses.

To perform the principal component analyses (PCA), we took a subsample of independent markers (pruning) with a maximum sliding window of 500,000 base pairs and a linkage disequilibrium (LD) threshold of 0.2 using the “snpgdsLDpruning” function from the SNPRelate R package. We ran the PCA with the “snpgdsPCA” function from SNPRelate package from R. The first 2 principal components (PCs) were considered significant.

For the GWAS analyses of the 283,246 remaining markers between EP6 cluster compared to all others, we used a logistic regression model, with the dichotomous variable indicating if the participants belong to the EP6 cluster as the outcome variable. We used the additive model for markers as an independent variable, and we adjusted the models with the first two PCs. Odds ratio and p-values were calculated on each model. Quantile–quantile (Q–Q) plots have been performed as quality control of the models; p-values were plotted using “qqplot.pvalues” function from gaston R package (67) (data not shown). We compared the aptamers’ normalized level of expression between the three groups of genotypes for each studied gene by performing standard analysis of variance (ANOVA) analyses followed by Tukey post hoc tests (referred to in this study as pQTL analysis). Since aptamers tested are limited compared to SNVs, we fixed significance p-value threshold below 1E-4 to report more SNVs instead of the more common 1E-5 suggestive threshold. Finally, to identify if EP6 may be characterized by the aptamers and the SNVs, we performed multiple logistic regression analysis including aptamers expression values, SNV genotypes (additive model), and the two principal components. Odds ratios are reported with 95% confidence intervals.

Metabolomic pathway characterization of EP6

The 1,435 metabolites measured were organized into 122 sub-pathways in the original dataset (denoted as “SUB_PATHWAYS”). We first identified metabolites whose values were significantly higher or lower in EP6 compared to other EPs (two-sided MWU test, FDR<0.01). Then, we used these metabolites to perform pathway enrichment analysis (one-sided Fisher’s exact test) based on 122 pathways. The resulting p-values were then corrected for multiple tests using Benjamini– Hochberg FDR.

Statistics

Several non-parametric tests, including the MWU test, Fisher’s exact test, and Spearman’s rank correlation, were used in this study. Benjamini–Hochberg FDR was used to adjust the p-values for multiple tests.

Study approval

The study was approved by the Institutional Ethics Review Board of the “Centre intégré universitaire de santé et de services sociaux du Saguenay-Lac-Saint-Jean” (CIUSSS-SLSJ) affiliated to the Université de Sherbrooke [protocol #2021-369, 2021-014 CMDO – COVID19].

Data Availability

Input data corresponding to the cohort can be obtained form BQC19 (www.quebeccovidbiobank.ca). The data generated as a result of the analyses are provided as tables and supplementary tables.

https://www.quebeccovidbiobank.ca

Conflict of Interests

The authors have declared that no conflict of interest exists.

Acknowledgements

This work was made possible through open sharing of data and samples from the Biobanque Québécoise de la COVID-19, funded by the Fonds de recherche du Québec - Santé, Génome Québec, the Public Health Agency of Canada and, as of March 2022, the Ministère de la Santé et des Services Sociaux du Québec. We thank all participants to BQC19 for their contribution. This study was supported by the Fonds de recherche du Québec - Santé (FRQS)-Cardiometabolic Health, Diabetes and Obesity Research Network (CMDO)-Initiative. This work was also supported by Natural Sciences and Engineering Research Council of Canada (NSERC) grant RGPIN-2019-04460 (AE).

References

  1. 1.↵
    Sanders JM, Monogue ML, Jodlowski TZ, Cutrell JB. Pharmacologic Treatments for Coronavirus Disease 2019 (COVID-19): A Review. JAMA 2020; 323: 1824–1836.
    OpenUrlPubMed
  2. 2.↵
    Lim ZJ, Subramaniam A, Ponnapa Reddy M, Blecher G, Kadam U, Afroz A, et al. Case Fatality Rates for Patients with COVID-19 Requiring Invasive Mechanical Ventilation. A Meta-analysis. Am J Respir Crit Care Med 2021; 203: 54–66.
    OpenUrlCrossRef
  3. 3.↵
    Berrettini WH. Genetic bases for endophenotypes in psychiatric disorders. Dialogues Clin Neurosci 2005; 7: 95–101.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Froidure A, Mouthuy J, Durham SR, Chanez P, Sibille Y, Pilette C. Asthma phenotypes and IgE responses. Eur Respir J 2016; 47: 304–19.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    Barnes PJ. Inflammatory endotypes in COPD. Allergy 2019; 74: 1249–1256.
    OpenUrl
  6. 6.↵
    Koennecke M, Klimek L, Mullol J, Gevaert P, Wollenberg B. Subtyping of polyposis nasi: phenotypes, endotypes and comorbidities. Allergo J Int 2018; 27: 56–65.
    OpenUrl
  7. 7.↵
    Baghela A, An A, Zhang P, Acton E, Gauthier J, Brunet-Ratnasingham E, et al. Predicting severity in COVID-19 disease using sepsis blood gene expression signatures. Sci Rep 2023; 13: 1247.
    OpenUrl
  8. 8.↵
    LaSalle TJ, Gonye ALK, Freeman SS, Kaplonek P, Gushterova I, Kays KR, et al. Longitudinal characterization of circulating neutrophils uncovers phenotypes associated with severity in hospitalized COVID-19 patients. Cell Rep Med 2022; 3: 100779.
    OpenUrl
  9. 9.↵
    Ma W, Soulé A, Tremblay K, Rousseau S, Emad A. A circulating proteome-informed prognostic model of COVID-19 disease activity that relies on routinely available clinical laboratories. medRxiv 2022: 2022.11.02.22281834.
  10. 10.↵
    Welin A, Amirbeagi F, Christenson K, Björkman L, Björnsdottir H, Forsman H, et al. The Human Neutrophil Subsets Defined by the Presence or Absence of OLFM4 Both Transmigrate into Tissue In Vivo and Give Rise to Distinct NETs In Vitro. PLOS ONE 2013; 8: e69575.
    OpenUrlCrossRefPubMed
  11. 11.↵
    Blatti III C, Emad A, Berry MJ, Gatzke L, Epstein M, Lanier D, et al. Knowledge-guided analysis of “omics” data using the KnowEnG cloud platform. PLoS biology 2020; 18: e3000583.
    OpenUrlCrossRefPubMed
  12. 12.↵
    Gillespie M, Jassal B, Stephan R, Milacic M, Rothfels K, Senff-Ribeiro A, et al. The reactome pathway knowledgebase 2022. Nucleic Acids Research 2021; 50: D687–D692.
    OpenUrlCrossRef
  13. 13.↵
    Blachier F, Andriamihaja M, Blais A. Sulfur-Containing Amino Acids and Lipid Metabolism. The Journal of Nutrition 2020; 150: 2524S–2531S.
    OpenUrlCrossRef
  14. 14.↵
    Ye C, Sutter BM, Wang Y, Kuang Z, Tu BP. A metabolic function for phospholipid and histone methylation. Molecular cell 2017; 66: 180-193. e8.
    OpenUrlCrossRefPubMed
  15. 15.↵
    Tavoosi N, Davis-Harrison RL, Pogorelov TV, Ohkubo YZ, Arcario MJ, Clay MC, et al. Molecular Determinants of Phospholipid Synergy in Blood Clotting. Journal of Biological Chemistry 2011; 286: 23247–23253.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    Parman T, Bunin DI, Ng HH, McDunn JE, Wulff JE, Wang A, et al. Toxicogenomics and metabolomics of pentamethylchromanol (PMCol)-induced hepatotoxicity. Toxicological Sciences 2011; 124: 487–501.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    Grohmann U, Volpi C, Fallarino F, Bozza S, Bianchi R, Vacca C, et al. Reverse signaling through GITR ligand enables dexamethasone to activate IDO in allergy. Nature Medicine 2007; 13: 579–586.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Marfia G, Navone S, Guarnaccia L, Campanella R, Mondoni M, Locatelli M, et al. Decreased serum level of sphingosine-1-phosphate: a novel predictor of clinical severity in COVID-EMBO molecular medicine 2021; 13: e13424.
    OpenUrl
  19. 19.↵
    Lone MA, Hülsmeier AJ, Saied EM, Karsai G, Arenz C, von Eckardstein A, et al. Subunit composition of the mammalian serine-palmitoyltransferase defines the spectrum of straight and methyl-branched long-chain bases. Proceedings of the National Academy of Sciences 2020; 117: 15591–15598.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    Han G, Gupta SD, Gable K, Niranjanakumari S, Moitra P, Eichler F, et al. Identification of small subunits of mammalian serine palmitoyltransferase that confer distinct acyl-CoA substrate specificities. Proceedings of the National Academy of Sciences 2009; 106: 8186–8191.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    Chen PY, Gao WY, Liou JW, Lin CY, Wu MJ, Yen JH. Angiopoietin-Like Protein 3 (ANGPTL3) Modulates Lipoprotein Metabolism and Dyslipidemia. Int J Mol Sci 2021; 22.
  22. 22.↵
    Lachkar F, Ferré P, Foufelle F, Papaioannou A. Dihydroceramides: their emerging physiological roles and functions in cancer and metabolic diseases. American Journal of Physiology-Endocrinology and Metabolism 2021; 320: E122–E130.
    OpenUrl
  23. 23.↵
    Koishi R, Ando Y, Ono M, Shimamura M, Yasumo H, Fujiwara T, et al. Angptl3 regulates lipid metabolism in mice. Nature Genetics 2002; 30: 151–157.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    Pisciotta L, Favari E, Magnolo L, Simonelli S, Adorni MP, Sallo R, et al. Characterization of Three Kindreds With Familial Combined Hypolipidemia Caused by Loss-of-Function Mutations of ANGPTL3. Circulation: Cardiovascular Genetics 2012; 5: 42–50.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    Lekka ME, Liokatis S, Nathanail C, Galani V, Nakos G. The Impact of Intravenous Fat Emulsion Administration in Acute Lung Injury. American Journal of Respiratory and Critical Care Medicine 2004; 169: 638–644.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    Ahmad Z, Banerjee P, Hamon S, Chan K-C, Bouzelmat A, Sasiela WJ, et al. Inhibition of Angiopoietin-Like Protein 3 With a Monoclonal Antibody Reduces Triglycerides in Hypertriglyceridemia. Circulation 2019; 140: 470–486.
    OpenUrl
  27. 27.↵
    Bie H, Yin J, He X, Kermode AR, Goddard-Borger ED, Withers SG, et al. Insights into mucopolysaccharidosis I from the structure and action of α-L-iduronidase. Nature Chemical Biology 2013; 9: 739–745.
    OpenUrl
  28. 28.↵
    Helms J, Tacquard C, Severac F, Leonard-Lorant I, Ohana M, Delabranche X, et al. High risk of thrombosis in patients with severe SARS-CoV-2 infection: a multicenter prospective cohort study. Intensive Care Medicine 2020; 46: 1089–1098.
    OpenUrlCrossRefPubMed
  29. 29.↵
    Asakura H, Ogawa H. COVID-19-associated coagulopathy and disseminated intravascular coagulation. International Journal of Hematology 2021; 113: 45–57.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Fahmy OH, Daas FM, Salunkhe V, Petrey JL, Cosar EF, Ramirez J, et al. Is Microthrombosis the Main Pathology in Coronavirus Disease 2019 Severity?-A Systematic Review of the Postmortem Pathologic Findings. Crit Care Explor 2021; 3: e0427.
    OpenUrlPubMed
  31. 31.↵
    Barnes BJ, Adrover JM, Baxter-Stoltzfus A, Borczuk A, Cools-Lartigue J, Crawford JM, et al. Targeting potential drivers of COVID-19: Neutrophil extracellular traps. Journal of Experimental Medicine 2020; 217.
  32. 32.
    Bonaventura A, Vecchié A, Dagna L, Martinod K, Dixon DL, Van Tassell BW, et al. Endothelial dysfunction and immunothrombosis as key pathogenic mechanisms in COVID-19 Nature Reviews Immunology 2021; 21: 319–329.
    OpenUrl
  33. 33.
    Ding J, Hostallero DE, El Khili MR, Fonseca GJ, Milette S, Noorah N, et al. A network-informed analysis of SARS-CoV-2 and hemophagocytic lymphohistiocytosis genes’ interactions points to Neutrophil extracellular traps as mediators of thrombosis in COVID-19 PLoS Computational Biology 2021; 17: e1008810.
    OpenUrl
  34. 34.
    Blasco A, Coronado M-J, Hernández-Terciado F, Martín P, Royuela A, Ramil E, et al. Assessment of neutrophil extracellular traps in coronary thrombus of a case series of patients with COVID-19 and myocardial infarction. JAMA cardiology 2021; 6: 469–474.
    OpenUrl
  35. 35.
    Desilles JP, Solo Nomenjanahary M, Consoli A, Ollivier V, Faille D, Bourrienne MC, et al. Impact of COVID-19 on thrombus composition and response to thrombolysis: Insights from a monocentric cohort population of COVID-19 patients with acute ischemic stroke. Journal of Thrombosis and Haemostasis 2022; 20: 919–928.
    OpenUrl
  36. 36.
    Englert H, Rangaswamy C, Deppermann C, Sperhake J-P, Krisp C, Schreier D, et al. Defective NET clearance contributes to sustained FXII activation in COVID-19-associated pulmonary thrombo-inflammation. EBioMedicine 2021; 67: 103382.
    OpenUrl
  37. 37.
    Leppkes M, Knopf J, Naschberger E, Lindemann A, Singh J, Herrmann I, et al. Vascular occlusion by neutrophil extracellular traps in COVID-19. EBioMedicine 2020; 58: 102925.
    OpenUrl
  38. 38.
    Middleton EA, He X-Y, Denorme F, Campbell RA, Ng D, Salvatore SP, et al. Neutrophil extracellular traps contribute to immunothrombosis in COVID-19 acute respiratory distress syndrome. Blood 2020; 136: 1169–1179.
    OpenUrlCrossRefPubMed
  39. 39.
    Obermayer A, Jakob L-M, Haslbauer JD, Matter MS, Tzankov A, Stoiber W. Neutrophil extracellular traps in fatal COVID-19-associated lung injury. Disease markers 2021; 2021.
  40. 40.
    Ouwendijk WJ, Raadsen MP, Van Kampen JJ, Verdijk RM, Von Der Thusen JH, Guo L, et al. High levels of neutrophil extracellular traps persist in the lower respiratory tract of critically ill patients with coronavirus disease 2019. The Journal of infectious diseases 2021; 223: 1512–1521.
    OpenUrl
  41. 41.
    Petito E, Falcinelli E, Paliani U, Cesari E, Vaudo G, Sebastiano M, et al. Association of neutrophil activation, more than platelet activation, with thrombotic complications in coronavirus disease 2019. The Journal of infectious diseases 2021; 223: 933–944.
    OpenUrl
  42. 42.↵
    Skendros P, Mitsios A, Chrysanthopoulou A, Mastellos DC, Metallidis S, Rafailidis P, et al. Complement and tissue factor–enriched neutrophil extracellular traps are key drivers in COVID-19 immunothrombosis. The Journal of clinical investigation 2020; 130: 6151–6157.
    OpenUrl
  43. 43.↵
    Massberg S, Grahl L, von Bruehl M-L, Manukyan D, Pfeiler S, Goosmann C, et al. Reciprocal coupling of coagulation and innate immunity via neutrophil serine proteases. Nature Medicine 2010; 16: 887–896.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    Papayannopoulos V, Metzler KD, Hakkim A, Zychlinsky A. Neutrophil elastase and myeloperoxidase regulate the formation of neutrophil extracellular traps. Journal of Cell Biology 2010; 191: 677–691.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    Alder MN, Opoka AM, Lahni P, Hildeman DA, Wong HR. Olfactomedin-4 Is a Candidate Marker for a Pathogenic Neutrophil Subset in Septic Shock. Critical Care Medicine 2017; 45: e426–e432.
    OpenUrl
  46. 46.↵
    Stafford JH, Thorpe PE. Increased exposure of phosphatidylethanolamine on the surface of tumor vascular endothelium. Neoplasia 2011; 13: 299–IN2.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    Ran S, Downes A, Thorpe PE. Increased exposure of anionic phospholipids on the surface of tumor blood vessels. Cancer research 2002; 62: 6132–6140.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    Port JR, Adney DR, Schwarz B, Schulz JE, Sturdevant DE, Smith BJ, et al. High-fat high-sugar diet-induced changes in the lipid metabolism are associated with mildly increased COVID-19 severity and delayed recovery in the Syrian hamster. Viruses 2021; 13: 2506.
    OpenUrl
  49. 49.↵
    MacDonald J, Sprecher H. Distribution of arachidonic acid in choline-and ethanolamine-containing phosphoglycerides in subfractionated human neutrophils. Journal of Biological Chemistry 1989; 264: 17718–17726.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    Papageorgiou C, Jourdi G, Adjambri E, Walborn A, Patel P, Fareed J, et al. Disseminated intravascular coagulation: an update on pathogenesis, diagnosis, and therapeutic strategies. Clinical and Applied Thrombosis/Hemostasis 2018; 24: 8S–28S.
    OpenUrlCrossRef
  51. 51.↵
    Lacal PM, Petrillo MG, Ruffini F, Muzi A, Bianchini R, Ronchetti S, et al. Glucocorticoid-induced tumor necrosis factor receptor family-related ligand triggering upregulates vascular cell adhesion molecule-1 and intercellular adhesion molecule-1 and promotes leukocyte adhesion. J Pharmacol Exp Ther 2013; 347: 164–72.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    Pandit R, Singh I, Ansari A, Raval J, Patel Z, Dixit R, et al. First report on genome wide association study in western Indian population reveals host genetic factors for COVID-19 severity and outcome. Genomics 2022; 114: 110399.
    OpenUrl
  53. 53.↵
    Goligher EC, Bradbury CA, McVerry BJ, Lawler PR, Berger JS, Gong MN, et al. Therapeutic Anticoagulation with Heparin in Critically Ill Patients with Covid-19. N Engl J Med 2021; 385: 777–789.
    OpenUrlCrossRefPubMed
  54. 54.↵
    Lawler PR, Goligher EC, Berger JS, Neal MD, McVerry BJ, Nicolau JC, et al. Therapeutic Anticoagulation with Heparin in Noncritically Ill Patients with Covid-19. N Engl J Med 2021; 385: 790–802.
    OpenUrlCrossRefPubMed
  55. 55.↵
    Smadja DM, Philippe A, Bory O, Gendron N, Beauvais A, Gruest M, et al. Placental growth factor level in plasma predicts COVID-19 severity and in-hospital mortality. Journal of Thrombosis and Haemostasis 2021; 19: 1823–1830.
    OpenUrl
  56. 56.↵
    Bussolari C, Palumbo D, Fominsky E, Nardelli P, De Lorenzo R, Vitali G, et al. Case Report: Nintedaninb May Accelerate Lung Recovery in Critical Coronavirus Disease 2019. Frontiers in Medicine 2021; 8: 766486.
    OpenUrl
  57. 57.↵
    Liu W, Peng Y, Tobin DJ. A new 12-gene diagnostic biomarker signature of melanoma revealed by integrated microarray analysis. PeerJ 2013; 1: e49.
    OpenUrlCrossRef
  58. 58.↵
    Luo S, Feofanova EV, Tin A, Tung S, Rhee EP, Coresh J, et al. Genome-wide association study of serum metabolites in the African American Study of Kidney Disease and Hypertension. Kidney International 2021; 100: 430–439.
    OpenUrl
  59. 59.↵
    Agarwal B, Gatt A, Riddell A, Wright G, Chowdary P, Jalan R, et al. Hemostasis in patients with acute kidney injury secondary to acute liver failure. Kidney International 2013; 84: 158–163.
    OpenUrlCrossRefPubMed
  60. 60.↵
    Williams R, Schalm SW, O’Grady JG. Acute liver failure: redefining the syndromes. The Lancet 1993; 342: 273–275.
    OpenUrl
  61. 61.↵
    Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nature Protocols 2009; 4: 1184–1191.
    OpenUrl
  62. 62.↵
    R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2022. Available from: https://www.R-project.org/.
  63. 63.↵
    Marees AT, de Kluiver H, Stringer S, Vorspan F, Curis E, Marie-Claire C, et al. A tutorial on conducting genome-wide association studies: Quality control and statistical analysis. Int J Methods Psychiatr Res 2018; 27: e1608.
    OpenUrlCrossRefPubMed
  64. 64.↵
    Ochoa A. genio: Genetics Input/Output Functions. R package version 1.1.1; 2022. Available from: https://CRAN.R-project.org/package=genio.
  65. 65.↵
    Warnes G, Gorjanc G, Leisch F, Man M. genetics: Population Genetics. R package version 1.3.8.1.3; 2021. Available from: https://CRAN.R-project.org/package=genetics.
  66. 66.↵
    Zheng X, Levine D, Shen J, Gogarten SM, Laurie C, Weir BS. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics 2012; 28: 3326–3328.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    Perdry H, Dandine-Roulland C. gaston: Genetic Data Handling (QC, GRM, LD, P.A. & Linear Mixed Models. R package version 1.5.7; 2020. Available from: https://CRAN.R-project.org/package=gaston.
Back to top
Posted March 10, 2023.
Download PDF

Supplementary Material

Data/Code
Email
Molecular pathology of acute respiratory distress syndrome, mechanical ventilation and abnormal coagulation in severe COVID-19
Antoine Soulé, William Ma, Katelyn Yixiu Liu, Catherine Allard, Salman Qureshi, Karine Tremblay, Amin Emad, Simon Rousseau
medRxiv 2023.03.09.23286797; doi: https://doi.org/10.1101/2023.03.09.23286797
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area