key: cord-1016262-hzkj6943
authors: Yu, Xiaosi; Xu, Xianqun; Wu, Tianpeng; Huang, Wenjie; Xu, Chen; Xie, Wen; Long, Xinghua
title: APOA1 Level is Negatively Correlated with the Severity of COVID-19
date: 2022-01-18
journal: Int J Gen Med
DOI: 10.2147/ijgm.s332956
sha: a91b297a41cc2262022a664ab3f67fe62943822e
doc_id: 1016262
cord_uid: hzkj6943

OBJECTIVE: We used public data to analyze the proteomics, metabolomics and transcriptomics characteristics of COVID-19 patients to identify potential therapeutic targets. More importantly, we also collected clinical data for verification to make the analysis results more reliable. METHODS: Download the serum proteomics and metabolomics data of COVID-19 patients and describe their changes in serum proteins and metabolites, and use bioinformatics analysis methods to identify potential biomarkers and therapeutic targets. Finally, clinical data and experimental data of cell infection were combined for verification. RESULTS: It was found that the serum apolipoprotein A1 (APOA1) protein level in COVID-19 patients was down-regulated (log2FC = −0.39, false discovery rate (FDR) < 0.001), and the degree of reduction in the severe group was more significant (kruskal-test p = 2.5e-05). What is more, APOA1 was not only expressed lower in male patients (Wilcox-test p = 0.012), but also negatively correlated with C-reactive protein (CRP, r = −0.37, p = 0.019). The experiment data from SARS-CoV-2 infected cells further showed that the protein and transcript level of APOA1 gradually decreased as the infection time increased, and the transcription level (log(2)FC = −8.3, FDR = 0.0015) was more down-regulated than protein level (log2FC = −0.95, FDR = 0.0014). More importantly, the collected clinical data also confirmed that APOA1 was down-regulated in COVID-19 patients (kruskal-test p = 0.001), and APOA1 levels are negatively correlated with IL6 (r = −0.396, p = 2.22e-07), D-dimers (DD, r = −0.262, p = 8.19e-04), prothrombin time (PT, r = −0.464, p = 6.68e-10) and thrombin time (TT, r = −0.279, p = 3.46e-04). CONCLUSION: The degree of down-regulation of APOA1 is positively correlated with the severity of COVID-19, and the expression level of APOA1 is negatively correlated with CRP, IL6, DD, PT, TT, and positively correlated with HD and LDL. This indicates that APOA1 may be a key molecule in tandem acute inflammatory response, coagulation abnormalities and cholesterol metabolism disorder in COVID-19, and could be a potential therapeutic target.

Coronavirus disease 2019 (COVID-19) is a rapidly spreading disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). It was first discovered in Wuhan, China, and has rapidly spread and become a major global disease. 1, 2 According to the World Health Organization (WHO, https://www.who.int/emergencies/diseases/ novel-coronavirus-2019), as of October 30, 2020, there were 44,351,506 confirmed cases in 219 countries, including 1,171,255 deaths, and the number of infections is still rising rapidly. While there are many clinical trials underway for COVID-19, there is no specific vaccine or effective treatment.

One of the known pathogenesis of COVID-19 is an excessive inflammatory response, 3 a sharp increase in proinflammatory cytokines, known as "cytokine storms", leading to respiratory distress and coagulation disorders in patients. The patient developed elevated D-dimer, deep vein thrombosis. 4 Cardiovascular injury has been reported in 20-30% of COVID19 patients. 5 And hyperlipidemia and diabetes impair immune responses and lead to persistent chronic inflammation, which increases the risk of cardiovascular disease (CVD). [6] [7] [8] The application of proteomics and metabolomics research technology has greatly promoted the progress of biological research and is of great significance for the identification of biomarkers. [9] [10] [11] In this study, we performed a series of statistical analyses based on the serum proteomics and metabolomics data of COVID-19 patients to determine potential therapeutic targets and find diagnostic biomarker. And verify it in the proteome and transcriptome data of the SARS-CoV-2 infected cell experiment. Finally, combined with the clinical data of the hospital, we successfully verified that the expression of APOA1 in the patient's serum was significantly downregulated, and that the degree of decline was positively correlated with the severity of the disease. This indicates that APOA1 is likely to be a potential therapeutic target for COVID-19.

Serum proteomic and metabolomic data for COVID-19 was downloaded from Shen et al. 11 Only 41 COVID-19 patients and 18 healthy controls with both proteomic and metabolomic data were included in the follow-up analysis. Urine proteomic data for COVID-19 were downloaded from Li et al, 12 data from 6 patients and 32 healthy control samples were included. Proteomic data and transcriptomic data from SARS-CoV-2 infected cell experiments were downloaded from Bojkova et al. 13 

The "Metaboanalyst" website (https://www.metaboanalyst. ca/) 14 was used for data standardization. We removed the feature of missing > 50% values, and replaced the missing values with 1/5 of the minimum positive value of each variable. Finally, we perform log transformations to standardize the data for subsequent analysis.

Gene Set Enrichment Analysis (GSEA) GSEA (https://www.gsea-msigdb.org/gsea/msigdb) 15 uses gene expression matrices and gene function annotation files to compare differences in gene function enrichment between the two groups. In order to explore the potential molecular mechanism of COVID-19 infection, GSEA was used to compare the differences in the KEGG signaling pathway between COVID-19 patients and healthy people. Nominal p-value < 0.05 was considered statistically significant.

"Limma" package 16 was used for differential analysis, | log2FC| > 0.25 and FDR < 0.05 were considered as the cutoff values for differentially expressed proteins (DEPs) screening. Volcano plot and heatmap were performed through the corresponding package of R (version 3.6.3). In addition, function enrichment analysis was carried out by the Metascape website (https://metascape.org/gp/index. html). 17 Weighted Gene Co-Expression Network Analysis (WGCNA) WGCNA 18 can identify highly synergistically altered gene sets and identify potential biomarker genes or therapeutic targets based on their internal connectivity and association with phenotypes. Therefore, the "WGCNA" package was used to put the features with similar expression patterns into the module by average linkage clustering to identify gene modules related to clinical traits.

The intersection of the significant module and the differentially expressed proteins were used for subsequent analysis. The PPI network was retrieved from STRING (https://string-db.org/) 19 website and visualized by Cytoscape software (version 3.6.0). 20 The "CytoHubba" 21 plug-in, which can use 12 different algorithms to calculate the importance of each point, were used to find important features. Besides, the plug-in of "MCODE" 22 was used to find clusters based on topology to locate densely connected regions. SPSS (Version 24.0) was used to calculate the poisson correlation. |Correlation coefficient (r)| > 0.2 and p < 0.05 were selected. 

We retrospectively collected blood test data of inpatients in Zhongnan Hospital of Wuhan University from January 20 to March 20. First, we screened patients with confirmed COVID-19. Then we screened out the data of 160 COVID-19 patients who have tested APOA1, APOB, HDL, LDL, APTT, TT, PT, DD and IL6 at the same time, of which 6 patients were diagnosed with severe pneumonia. At the same time, we also collected 31 healthy control of ApoA1, APOB, HDL and LDL blood test data during routine health examination in January 2020.

This study was reviewed and approved by the Ethical Committee of Zhongnan Hospital of Wuhan University. This was a retrospective study, written informed consent was waived by the Ethics Commission of the designated hospital for emerging infectious diseases. The data was kept confidential and complies with the Declaration of Helsinki.

The serum proteome data of Shen et al were analyzed by GSEA to explore the changes in the KEGG signaling pathway between COVID-19 patients and healthy controls. Among them, 7 pathways were significantly different. The two KEGG signaling pathways, complement and coagulation cascade and systemic lupus erythematosus, were significantly up-regulated in COVID-19 patients ( Figure 1A and B). Pathways in cancer, ECM receptor interaction, focal adhesion, chemokine and TGF-β signaling pathway were significantly down-regulated in healthy people ( Figure 1C -G). This indicates that these signaling pathways may be related to the pathogenesis of COVID-19.

Analysis of the differential expression of serum proteins in COVID-19 patients showed that 67 proteins were upregulated and 44 proteins were down-regulated ( Figure 1H and I). Enrichment analysis showed that the top three significantly enriched terms in up-regulated proteins were "acute inflammatory response", "acute-phase response" and "negative regulation of proteolysis" ( Figure 1J ). While in down-regulated proteins were "platelet degranulation", "extracellular structure organization" and "high-density lipoprotein particle remodeling" ( Figure 1K ). These were consistent with the known excessive inflammatory response and coagulation dysfunction in patients with COVID-19.

WGCNA was used to put the proteins with similar expression patterns into modules and analyzed the correlation between each module and disease severity, gender, age and BMI. It was found that turquoise module was significantly positively correlated with disease severity (r = 0.87, p = 2e-18) and had the highest connectivity (r = 0.94, p =1.8e-115) (Figure 2A-D) . The turquoise module contains 245 proteins, of which 116 proteins with p < 0.05 were considered as significant proteins. The enrichment analysis showed that the top three biological functions of significant enrichment were "acute inflammatory response", "platelet degranulation", and "regulation of insulin-like growth factor (IGF) transport and uptake by insulin-like growth factor Bi" ( Figure 2E ). These results further confirmed the important role of acute inflammation and platelet degranulation in the pathogenesis of COVID-19.

A total of 73 proteins with 111 DEPs and 116 significant proteins intersections were used for subsequent analysis ( Figure 3A ). The protein interaction network was constructed through the "String" website and introduced into "Cytoscape" for visualization. "CytoHubba" plug-in was used to screen the key genes in the PPI network. It was found that APOA1 have the most significant among various algorithms. Figure 3B showed the top 5 notable proteins in degree, betweenness, closeness, EcCentricity and maximal clique centrality (MCC). Through the "MCODE" plug-in, we obtained 3 modules. APOA1 was present in MCODE1, which contains 22 proteins, and enrichment analysis showed that it was significantly related to "platelet degranulation" and "acute phase response" (Figure 3C and D). We further compared the expression of APOA1 protein between different groups, and found that the expression level of severe group was the lowest (kruskal-test p = 2.5e-05, Figure 3E ). Moreover, the APOA1 protein level of men was lower than that of women (Wilcox-test p = 0.012, Figure 3F Figure 3G ).

In order to further understand the changes of related metabolites, a difference analysis was carried out on the metabolome data. 224 metabolites were differentially expressed (DEMs, |log2FC| > 0.25, FDR < 0.05), of which 32 were up-regulated and 192 were down-regulated ( Figure 3H ). Subsequently, we analyzed the correlation between the DEMs and APOA1. 25 DEMs were significantly positively correlated with APOA1, and 3 were significantly negatively correlated (|r| > 0.2, p < 0.05, Figure 3I ). Among them, the top three with the highest correlation are docosahexaenoylcholine, 1-(1-alkenyl-palmitoyl) -2-linoleoyl-GPC (P-16:0/18:2) and 1-(1-alkenyl-stearoyl) -2-linoleoyl-GPE (P-18:0/18:2). These metabolites were closely related to cholesterol metabolism. 

In the urine proteomic data of COVID-19 patients, APOA1 was down-regulated as well (Wilcox-test p = 0.00058), which further confirms the analysis results of plasma proteomics ( Figure 4A ).

The analysis of the proteome and transcriptome data of the SARS-CoV-2 infection cell experiment found that with the increase of infection time, the protein and transcription levels of APOA1 decreased significantly ( Figure 4B and C). In order to further understand the relevant transcriptional regulation mechanism, the differential expression analysis was carried out on the transcription data of transfection 24h, and 99 genes were differentially expressed (DEGs, |log2FC| >2, FDR < 0.05, Figure 4D ). Next, the correlation between DEGs and APOA1 protein expression was analyzed (|r| > 0.5, p < 0.05). It was found that these DEGs were significantly correlated with APOA1 protein expression, of which 77 down-regulated genes were significantly positively correlated, and 22 up-regulated genes were significantly negatively correlated ( Figure 4E ). Enrichment analysis showed that "cholesterol homeostasis", "lysosome" and "signaling by WNT" were significantly enriched in 77 genes positively correlated with APOA1 protein expression ( Figure 4F ). Meanwhile, "sig regulation of the actin cytosketon by rho gtpases", "platelet activation" and "salmonella infection" were significantly enriched in 22 negatively related genes ( Figure 4G ).

In order to clinically verify the analysis results, blood test data were collected from inpatients confirmed with 

As we all known excessive inflammation, coagulation dysfunction, and elevated IL6 and D-dimers were the most common pathological features in patients with COVID-19. This means that APOA1 was closely related to the pathogenesis of COVID-19 and was likely to be a potential therapeutic target.

At present, COVID-19 has caused millions of deaths. Although there are a large number of clinical trials in progress, there is still a lack of effective treatment methods. 23 What is exciting is that many breakthrough results have been achieved in the research of COVID-19, 24,25 but it is undeniable that we have not fully grasped the pathogenic mechanism of this new virus. 2, 26 In this study, we comprehensively analyzed the changes in the serum proteome and metabolome of COVID-19 to explore the underlying pathogenesis and therapeutic targets. Moreover, we also used SARS-CoV-2 infected cell experimental data and COVID-19 patient's clinical examination data for experimental and clinical verification. It was found that the APOA1 protein level and transcription level of patients were down-regulated, and the transcription level was down-regulated more significantly. What's more, the more severe the disease, the lower the expression level of APOA1.

We used bioinformatics to analyze the serum proteome and metabolome data of COVID-19 patients from Shen et al. GSEA analysis shown that the two KEGG signaling pathways, "complement and coagulation cascade" and "systemic lupus erythematosus", were significantly enriched in COVID-19 patients. The enrichment analysis of DEPs showed that the up-regulated proteins were mainly related to "acute inflammatory response", and the down-regulated proteins were mainly related to "platelet degranulation". In addition, the subsequent gene enrichment results were highly consistent with these. According to reports, invasive inflammation caused by virus replication was the main pathogenic mechanism of SARS-CoV infection. 27 In COVID-19, patients' inflammatory markers also increased significantly, 28 and even cytokine storm appeared. 29 Moreover, COVID-19 patients have also been reported to be at high risk of thrombotic complications, especially in critically ill patients. 4, 30 These indicated that acute inflammation and venous thromboembolism were both important pathogenic mechanisms of SRSA-Cov-2, which are consistent with our analysis results.

Further analysis showed that APOA1 was not only differentially expressed in COVID-19 patients, but also negatively related to the severity of the disease, which was a critical protein in COVID-19. The serum protein data of Shen et al 11 shown that APOA1 was downregulated in COVID-19 patients, and down-regulated more significantly in critically ill patients and men, and its expression level was negatively correlated with CRP. Analysis of APOA1 related metabolites found that these metabolites were closely related to cholesterol metabolism. Furthermore, the clinical data we collected further showed that APOA1 was down-regulated more significantly in severe patients, and its expression level was positively correlated with HDL and LDL, and negatively correlated with IL6, DD, PT, and TT. CRP reflected the patient's acute inflammatory response, and IL6 was a key factor in the cytokine storm caused by excessive inflammation. IL6 and CRP were important indicators to assess the severity of patients and predict the outcome. 31 DD, PT, and TT reflected the coagulation and fibrinolysis status. Abnormal coagulation was also a common complication of COVID-19 patients, and it was positively related to the risk of death. 32 The most common coagulation abnormalities in COVID-19 patients were elevated DD and PT. 33 Low levels of HDL and LDL have been reported to increase the risk of COVID-19 developing into a serious event. [34] [35] [36] Studies have suggested that targeting lipid rafts could be a potential therapy for COVID-19. 37 Excitingly, our analysis results were completely consistent with the reported articles, which further confirmed the correctness of our analysis results. APOA1 was closely related to acute inflammation, cytokine storm, coagulation abnormalities and cholesterol metabolism. It concatenated nonspecific symptoms of COVID-19 and was likely to be a key molecule in unraveling the pathogenesis of COVID-19.

APOA1 is a major protein component of HDL, which promotes the extranodation of cholesterol from tissues to the liver, and is involved in inflammatory response and immune regulation. 38 APOA1 can regulate platelet activity and was used in the treatment of atherosclerosis. 39, 40 Moreover, it has been found that APOA1 has antiviral activity 41 which not only shows host defense function in Yersinia enterocolitis, 42 but also targets lysosomal membrane to promote the dissolution of trypanosomes. 43 However, the inhibitory effect of APOA1 on SARS-CoV-2 has not been realized.

In summary, using bioinformatics analysis and combined with sarS-COV-2 infected cell experiments and 696 clinical data, we found that the expression level of APOA was negatively correlated with the severity of COVID-19 patients, and significantly correlated with acute inflammation, coagulation abnormalities, and cholesterol metabolism disorder caused by SARS-COV-2 infection. APOA1 is likely to be a key molecule in decoding the pathogenesis of SARS-COV-2 and a potential target for the treatment of COVID-19. In the current study, the potential link between APOA1 and COVID-19 pathogenesis was clearly proposed for the first time. However, the research is limited to sample size for data analysis and theoretical verification, and further study was needed. The research can provide new biomarker and target for diagnosis and treatment of severe COVID-19.

COVID-19, coronavirus disease 2019; SARS-CoV-2, severe acute respiratory syndrome coronavirus-2; APOA1, apolipoprotein A1; WGCNA, weighted gene coexpression network analysis; FDR, false discovery rate; GSEA, gene set enrichment analysis; PPI, protein-protein interaction networks; R, correlation coefficient; DEPs, differentially expressed proteins; DEMs, differentially expressed metabolites; DEGs, differentially expressed genes; CRP, C-reactive protein; DD, D-dimers; PT, prothrombin time; TT, thrombin time; HDL, high density lipoprotein; LDL, low density lipoprotein.

All data used are contained within the article.

A novel coronavirus outbreak of global health concern

Covid-19: WHO declares pandemic because of "alarming levels" of spread, severity, and inaction

Correlation analysis between disease severity and inflammation-related parameters in patients with COVID-19 pneumonia

Incidence of venous thromboembolism in hospitalized patients with COVID-19

COVID-19 and the Heart

Disease tolerance and immunity in host protection against infection

Inflammation, immunity, and infection in atherothrombosis: JACC review topic of the week

Benefits of flu vaccination for persons with diabetes mellitus: a review

Proteomics and metabolomics in kidney disease, including insights into etiology, treatment, and prevention

Identification of bioactive metabolites using activity metabolomics

Proteomic and metabolomic characterization of COVID-19 patient sera

Urine proteome of COVID-19 patients

Proteomics of SARS-CoV-2-infected host cells reveals therapy targets

Using metaboAnalyst 4.0 for comprehensive and integrative metabolomics data analysis

Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles

limma powers differential expression analyses for RNA-sequencing and microarray studies

Metascape provides a biologist-oriented resource for the analysis of systems-level datasets

WGCNA: an R package for weighted correlation network analysis

STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

Cytoscape: a software environment for integrated models of biomolecular interaction networks

cytoHubba: identifying hub objects and sub-networks from complex interactome

An automated method for finding molecular complexes in large protein interaction networks

Emergence of Novel Coronavirus 2019-nCoV: need for rapid vaccine and biologics development

Remdesivir as a possible therapeutic option for the COVID-19

Breakthrough: chloroquine phosphate has shown apparent efficacy in treatment of COVID-19 associated pneumonia in clinical studies

Potential rapid diagnostics, vaccine and therapeutics for 2019 Novel Coronavirus (2019-nCoV): a systematic review

Plasma inflammatory cytokines and chemokines in severe acute respiratory syndrome

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

The cytokine storm and COVID-19

Confirmation of the high cumulative incidence of thrombotic complications in critically ill ICU patients with COVID-19: an updated analysis

Prognostic value of interleukin-6, C-reactive protein, and procalcitonin in patients with COVID-19

Coagulation abnormalities and thrombosis in patients with COVID-19

ISTH interim guidance on recognition and management of coagulopathy in COVID-19

The potential role of dyslipidemia in COVID-19 severity: an umbrella review of systematic reviews

Hypolipidemia is associated with the severity of COVID-19

Declined serum high density lipoprotein cholesterol is associated with the severity of COVID-19 infection

Targeting lipid rafts-a potential therapy for COVID-19

Apolipoprotein A-I (ApoA-I), immunity, inflammation and cancer

Evaluation of potential antiplatelet effects of CSL112 (Apolipoprotein A-I [Human]) in patients with atherosclerosis: results from a phase 2a study

Lipoproteins as biomarkers and therapeutic targets in the setting of acute coronary syndrome

HDL in innate and adaptive immunity

Apolipoprotein A-I exerts bactericidal activity against Yersinia enterocolitica serotype O:3

Apolipoprotein L-I promotes trypanosome lysis by forming pores in lysosomal membranes

This work was supported by grant from National Natural Science Foundation of China (81272372) and by Zhongnan Hospital of Wuhan University Science, Technology and Innovation Seed Fund (znpy2016033). The Project-sponsored by SRF for ROCS, SEM.

All authors declare no conflicts of interest in relation to this work.

The International Journal of General Medicine is an international, peer-reviewed open-access journal that focuses on general and internal medicine, pathogenesis, epidemiology, diagnosis, monitoring and treatment protocols. The journal is characterized by the rapid reporting of reviews, original research and clinical studies across all disease areas. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/ testimonials.php to read real quotes from published authors.