ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 18 August 2022

Sec. Atherosclerosis and Vascular Medicine

Volume 9 - 2022 | https://doi.org/10.3389/fcvm.2022.907665

Bioinformatics approach to identify the influences of SARS-COV2 infections on atherosclerosis

  • Department of Neurology, The First Affiliated Hospital of Harbin Medical University, Harbin, China

Article metrics

View details

10

Citations

3,3k

Views

1,6k

Downloads

Abstract

Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been a global pandemic since early 2020. Understanding the relationship between various systemic disease and COVID-19 through disease ontology (DO) analysis, an approach based on disease similarity studies, has found that COVID-19 is most strongly associated with atherosclerosis. The study provides new insights for the common pathogenesis of COVID-19 and atherosclerosis by looking for common transcriptional features. Two datasets (GSE152418 and GSE100927) were downloaded from GEO database to search for common differentially expressed genes (DEGs) and shared pathways. A total of 34 DEGs were identified. Among them, ten hub genes with high degrees of connectivity were picked out, namely C1QA, C1QB, C1QC, CD163, SIGLEC1, APOE, MS4A4A, VSIG4, CCR1 and STAB1. This study suggests the critical role played by Complement and coagulation cascades in COVID-19 and atherosclerosis. Our findings underscore the importance of C1q in the pathogenesis of COVID-19 and atherosclerosis. Activation of the complement system can lead to endothelial dysfunction. The DEGs identified in this study provide new biomarkers and potential therapeutic targets for the prevention of atherosclerosis.

Introduction

Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been a global pandemic since early 2020. According to the WHO, until March 9, 2022, the number of confirmed cases worldwide was 448,313,293, including 6,011,482 deaths (1). In response to the COVID-19 pandemic, a global effort is in progress to develop a vaccine against SARS-CoV-2. Vaccination can help control SARS-CoV-2 outbreaks by preventing infection, reducing disease severity, and blocking transmission (2). COVID-19 is typically characterized by upper respiratory symptoms, including fever, cough, and fatigue, and it is often accompanied by pulmonary infection (3). In addition to typical symptoms, some patients have serious cardiovascular damage, or even the first symptoms. (4) Except for the traditional established risk factors for atherosclerosis, such as age, smoking, hyperlipidemia, and hypertension, viral infection has been supposed to be a potential implication in atherosclerosis (5). SARS-CoV-2 binds to ACE2 to gain intracellular entry, leading to endothelial dysfunction (6). SARS-CoV-2 also promotes the accumulation of perivascular adipose tissue (7). These may exacerbate the underlying pathology of cardiovascular disease, leading to accelerated progression of atherosclerosis.

The purpose of this study was to explore the pathophysiological association between SARS-CoV-2 and atherosclerosis, and to better understand the underlying mechanisms, so as to facilitate early detection and prevention of atherosclerosis. Two gene expression datasets (GSE152418 and GSE100927) were downloaded from Gene Expression Omnibus (GEO) database. We used bioinformatics and enrichment analysis to determine the common DEGs and their functions for COVID-19 and atherosclerosis. In addition, protein protein interaction (PPI) networks were established to reveal hub genes. These data can better understand the potential link between the two diseases and provide evidence for therapeutic targets.

Materials and methods

Microarray data

The GSE152418 and GSE100927 gene expression profile were obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo) [Illumina NovaSeq 6000 (Homo sapiens)] platform was used for the GSE152418 dataset where samples were got from seventeen COVID-19 patients, and seventeen healthy people. On the contrary, for the GSE100927 dataset, GPL17077 [Agilent-039494 SurePrint G3 Human GE v2 8x60K Microarray 039381 (Probe Name version)] platform was adopted where samples were collected from sixty-nine atherosclerotic patients and thirty-five control subjects (Table 1).

Table 1

Disease nameDataset IDSubjectsGEO platformNumber of samples(control/disease)
COVID-19GSE152418Peripheral blood mononuclear cellGPL2467617/17
AtherosclerosisGSE100927Carotid, femoral and infra-popliteal arteriesGPL1707735/69

Basic information of the two microarray databases derived from the GEO database.

GEO, Gene Expression Omnibus; COVID-19, Coron a Virus Disease 2019.

Disease ontology (DO) analysis

Firstly, the “edgeR” package was applied to screen Differentially Expressed Genes (DEGs) from the GSE152418 dataset and a adjusted P less than 0.05, and |log2 Fold change (FC)| more than or equal to 1 was set as a cut-off point for selecting DEGs. The “DOSE” (8). and “ClusterProfiler” (9). packages were then used for DO analysis to study the disease mechanism by looking for disease correlation. DO analysis is a method based on the study of disease similarity, and it plays a vital role in understanding the pathogenesis of complex diseases, the early prevention and diagnosis of major diseases, new drug development, and drug safety evaluation.

Acquisition of common genes

The LIMMA package was used to detect the DEGs between atherosclerotic patients and healthy control from the GSE100927 dataset, and the adjusted P-value and |log2FC| were calculated. Genes that met the cutoff criteria, adjusted P < 0.05 and |log2FC| more than or equal to 1.0, were considered as DEGs. Then the common genes of the GSE152418 and GSE100927 sets were identified by using the Venn diagram webtool (bioinformatics.psb.ugent.be/webtools/Venn/).

Enrichment analysis of common genes

To further analyze biological processes of common DEGs, GO annotation analysis and KEGG pathway enrichment analysis were carried out through the Database for Annotation, Visualization and Integrated Discovery [DAVID (2021 Update), https://david.ncifcrf.gov/]. P-Value < 0.05 was used as the enrichment screening condition.

Construction PPI network and selection hub genes

The PPI network was predicted using Search Tool for the Retrieval of Interacting Genes (STRING, version 11.5, http://string-db.org/) online database. The PPI pairs were extracted with a interaction score more than or equal to 0.15, and then the PPI network was visualized by Cytoscape software (www.cytoscape.org/). Here, we used Degree to evaluate and select hub genes.

Results

DO analysis

Based on the cut-off criteria of adjusted P < 0.05 and |log2FC| more than or equal to 1, a total of 2080 DEGs were identified from GSE152418, including 1905 upregulated genes and 175 downregulated genes. p.adjust < 0.05 and gene counts more than or equal to 20 were used as the DO screening condition. Figure 1 shows the top ten most significantly enriched diseases, coronary artery disease, atherosclerosis, arteriosclerotic cardiovascular disease, arteriosclerosis, myocardial infarction, congestive heart failure, acute myocardial infarction, pulmonary hypertension, focal epilepsy and temporal lobe epilepsy. The enrichment results of other diseases by DO analysis are shown in Table 2.

Figure 1

Figure 1

Disease Ontology (DO) analysis of the DEGs from the GSE152418 dataset. The size of the circle represents the number of genes involved, and the abscissa represents the frequency of the genes involved in the term total genes.

Table 2

DO IDDescriptionCountP-ValueP-Adjust
DOID:3393Coronary artery disease856.97E−095.83E−06
DOID:5844Myocardial infarction712.17E−089.07E−06
DOID:2234Focal epilepsy241.52E−074.24E−05
DOID:6000Congestive heart failure592.20E−074.59E−05
DOID:3328Temporal lobe epilepsy213.18E−075.32E−05
DOID:1936Atherosclerosis793.97E−075.38E−05
DOID:2348Arteriosclerotic cardiovascular disease794.50E−075.38E−05
DOID:6432Pulmonary hypertension279.09E−079.50E−05
DOID:9408Acute myocardial infarction301.34E−060.000124898
DOID:2349Arteriosclerosis791.70E−060.000142212
DOID:5679Retinal disease777.92E−060.000601666
DOID:1168Familial hyperlipidemia261.20E−050.000834386
DOID:3146Lipid metabolism disorder281.34E−050.00085861
DOID:1793Pancreatic cancer691.94E−050.00115929
DOID:850Lung disease982.75E−050.001532539
DOID:4450Renal cell carcinoma723.34E−050.001746893
DOID:8466Retinal degeneration583.98E−050.001958846
DOID:6364Migraine234.23E−050.00196426
DOID:3324Mood disorder455.17E−050.00227529
DOID:0080000Muscular disease825.54E−050.002316532
DOID:263Kidney cancer868.04E−050.003053437
DOID:3459Breast carcinoma779.26E−050.003364892
DOID:0060037Developmental disorder of mental health750.0001136020.003794572
DOID:1826Epilepsy syndrome460.0001189770.003794572
DOID:4451Renal carcinoma760.0001225520.003794572
DOID:1686Glaucoma290.000145540.004217764
DOID:2355Anemia530.00015870.004416723
DOID:2742Auditory system disease270.0001872630.00470284
DOID:936Brain disease870.0001922210.00470284
DOID:15Reproductive system disease760.0002058830.00470284
DOID:120Female reproductive organ cancer870.0002077720.00470284
DOID:74Hematopoietic system disease900.000208140.00470284
DOID:3996Urinary system cancer940.0002172880.00473576
DOID:4074Pancreas adenocarcinoma380.0002349280.00473576
DOID:0060040Pervasive developmental disorder450.0002396730.00473576
DOID:0060116Sensory system cancer340.0002413850.00473576
DOID:2174Ocular cancer340.0002413850.00473576
DOID:0060041Autism spectrum disorder430.0002549150.00473576
DOID:12849Autistic disorder430.0002549150.00473576
DOID:18Urinary system disease900.0002808640.005097229
DOID:3083Chronic obstructive pulmonary disease480.0002865670.005097229
DOID:423Myopathy770.000330990.005647089
DOID:66Muscle tissue disease770.000330990.005647089
DOID:374Nutrition disease670.0003990540.006672186
DOID:0050700Cardiomyopathy350.0004493180.007240401
DOID:6713Cerebrovascular disease290.0004590210.007240401
DOID:654Overnutrition640.0004965620.007574304
DOID:4905Pancreatic carcinoma500.0004983090.007574304
DOID:557Kidney disease860.0005223920.007798559
DOID:4645Retinal cancer290.0006208660.008650729
DOID:9970Obesity620.0006634980.008946528
DOID:1115Sarcoma420.0007112310.009147531
DOID:229Female reproductive system disease420.0007112310.009147531
DOID:2320Obstructive lung disease610.0007344420.009302929
DOID:10534Stomach cancer550.0009504680.011515815
DOID:768Retinoblastoma280.0010339770.012138426
DOID:771Retinal cell cancer280.0010339770.012138426
DOID:3312Bipolar disorder350.0010916380.012501503
DOID:9352Type 2 diabetes mellitus450.0011107940.012548966
DOID:5041Esophageal cancer330.0011581250.012909231
DOID:3770Pulmonary fibrosis290.0012613540.013519123
DOID:403Mouth disease400.0014846730.015136423
DOID:26Pancreas disease370.0016054150.015977706
DOID:633Myositis210.0016482370.016210895
DOID:0060085Organ system benign neoplasm520.001679460.01632591
DOID:4607Biliary tract cancer380.0018077650.017371169
DOID:0060084Cell type benign neoplasm850.0020610290.019359779
DOID:657Adenoma640.0021097740.019597453
DOID:1074Kidney failure340.0021454240.019709612
DOID:48Male reproductive system disease300.0022006060.019996808
DOID:865Vasculitis280.0023255290.020904752
DOID:8398Osteoarthritis390.0024735370.021998689
DOID:0060100Musculoskeletal system cancer790.0025082470.022023954
DOID:299Adenocarcinoma340.0026767630.022177269
DOID:5223Infertility420.0026952570.022177269
DOID:3082Interstitial lung disease360.002719930.022177269
DOID:1575Rheumatic disease390.0027323670.022177269
DOID:418Systemic scleroderma390.0027323670.022177269
DOID:419Scleroderma390.0027323670.022177269
DOID:2394Ovarian cancer590.0027808590.022331756
DOID:201Connective tissue cancer680.0028949520.022331756
DOID:10952Nephritis320.0029116760.022331756
DOID:3620Central nervous system cancer280.0029881880.022505631
DOID:0080015Physical disorder300.0031473880.023285098
DOID:4960Bone marrow cancer610.0035375220.025494557
DOID:0070004Myeloma600.0039059170.026264796
DOID:2621Autonomic nervous system neoplasm680.0040493390.026264796
DOID:769Neuroblastoma680.0040493390.026264796
DOID:1091Tooth disease340.0040842390.026264796
DOID:10825Essential hypertension270.0047081320.028812356
DOID:289Endometriosis210.004832720.029276479
DOID:854Collagen disease400.0052163310.03128014
DOID:1107Esophageal carcinoma270.0052900250.031364972
DOID:0050737Autosomal recessive disease610.0053655840.031588932
DOID:0060036Intrinsic cardiomyopathy290.005444470.03160817
DOID:127Leiomyoma220.0058284570.032922905
DOID:37Skin disease630.006560530.035847078
DOID:1192Peripheral nervous system neoplasm700.0066798090.036261818
DOID:552Pneumonia240.0070579170.037823198
DOID:4766Embryoma630.0074507640.039423029
DOID:3388Periodontal disease290.0083003140.04322301
DOID:16Integumentary system disease690.0087585120.044109131
DOID:0060038Specific developmental disorder420.0090062640.044816886
DOID:12930Dilated cardiomyopathy210.0093801290.04640111
DOID:230Lateral sclerosis240.0099877270.047987012

Significantly enriched DO terms of DEGs.

DEG, Differentially Expressed Gene; DO, Disease Ontology; ID, Identity Document.

Identification of common DEGs

From GSE100927, 418 DEGs including 295 upregulated genes and 123 downregulated genes were identified. We analyzed the intersection of the DEG profiles using Venn (Figure 2). Ultimately, 34 DEGs were significantly differentially expressed in two datasets, of which 33 were significantly upregulated genes and 1 was downregulated gene.

Figure 2

Figure 2

Venn diagram. (A) Upregulated common DEGs of the GSE152418 and GSE100927 datasets. (B) Downregulated common DEGs of the GSE152418 and GSE100927 datasets.

Gene ontology and pathway enrichment analysis

GO and KEGG pathway analyses for DEGs were performed using the DAVID. The biological processes of DEGs were primarily associated with synapse disassembly, complement activation and innate immune response. For the cell component, the DEGs were enriched in extracellular region, blood microparticle, hemoglobin complex, collagen trimer, and so on. Molecular functions analysis showed that the DEGs were significantly enriched in oxygen transporter activity, oxygen binding, scavenger receptor activity, voltagE–gated potassium channel activity involved in atrial cardiac muscle cell action potential repolarization, phosphatidylcholinE–sterol O-acyltransferase activator activity, haptoglobin binding, organic acid binding and heme binding (Table 3). In addition, the KEGG pathway analysis showed that the DEGs were significantly enriched in Complement and coagulation cascades, Pertussis, Coronavirus disease—COVID-19, Staphylococcus aureus infection, Chagas disease, Systemic lupus erythematosus and Alcoholic liver disease (Table 4).

Table 3

GO IDDescriptionCountP-ValueGenes
Biological process
GO:0098883Synapse disassembly36.04E−05C1QB, C1QA, C1QC
GO:0006958Complement activation, classical pathway58.67E−05C1QB, C1QA, IGLL5, IGLL1, C1QC
GO:0045087Innate immune response72.32E−04C1QB, C1QA, IGLL5, VNN1, IGLL1, C1QC, OASL
GO:0006898Receptor-mediated endocytosis40.001825851CD163, STAB1, HBA2, APOE
GO:0006954Inflammatory response50.002915761CCR1, VNN1, STAB1, SPP1, SIGLEC1
GO:0098869Cellular oxidant detoxification30.00591978HBA2, HBD, APOE
GO:0042159Lipoprotein catabolic process20.007472267APOE, CTSD
GO:0098914Membrane repolarization during atrial cardiac muscle cell action potential20.007472267KCNJ5, KCNA5
GO:0006956Complement activation30.008884327C1QB, C1QA, C1QC
GO:0034447Very-low-density lipoprotein particle clearance20.008960238APOC1, APOE
GO:0034382Chylomicron remnant clearance20.010446055APOC1, APOE
GO:0030449Regulation of complement activation30.012175792C1QB, C1QA, C1QC
GO:0010873Positive regulation of cholesterol esterification20.013411239APOC1, APOE
GO:0033700Phospholipid efflux20.017842935APOC1, APOE
GO:0044267Cellular protein metabolic process30.021136392MMP1, SPP1, APOE
GO:0015671Oxygen transport20.022255408HBA2, HBD
GO:0015909Long-chain fatty acid transport20.025186416FABP5, APOE
GO:0034375High-density lipoprotein particle remodeling20.026648738APOC1, APOE
GO:0042157Lipoprotein metabolic process20.032476871APOC1, APOE
GO:0033344Cholesterol efflux20.036825845APOC1, APOE
GO:0045671Negative regulation of osteoclast differentiation20.039714668MAFB, LILRB4
GO:0032703Negative regulation of interleukin-2 production20.041155941VSIG4, LILRB4
GO:0042744Hydrogen peroxide catabolic process20.041155941HBA2, HBD
GO:0010033Response to organic substance20.042595125AQP9, KCNA5
GO:0007267Cell-cell signaling30.045927718CCR1, C1QA, STAB1
GO:0042742Defense response to bacterium30.046288292IGLL5, IGLL1, STAB1
Cellular component
GO:0005576Extracellular region172.40E−08C1QB, C1QA, CD163, CD163L1, MMP1, HBA2, VNN1, FNDC1, FABP5, IGLL1, APOC1, SPP1, PLBD1, SIGLEC1, APOE, CTSD, C1QC
GO:0072562Blood microparticle57.68E−05C1QB, HBA2, HBD, APOE, C1QC
GO:0005833Hemoglobin complex32.19E−04HBA2, HBD
GO:0005581Collagen trimer44.28E−04C1QB, C1QA, MMP1, C1QC
GO:0009897External side of plasma membrane66.43E−04CCR1, KCNJ5, IGLL5, CD163, CD163L1, IGLL1
GO:0005602Complement component C1 complex20.003171533C1QB, C1QA
GO:0098794Postsynapse30.012921039C1QB, C1QA, C1QC
GO:0031838Haptoglobin-hemoglobin complex20.017323254HBA2, HBD
GO:0042627Chylomicron20.021997095APOC1, APOE
GO:0071682Endocytic vesicle lumen20.028195396HBA2, APOE
GO:0016021Integral component of membrane150.028427769PTCRA, CCR1, KCNJ5, CD163, CD163L1, AQP9, KCNA5, HBD, LILRB4, MS4A4A, VNN1, SLCO2B1, STAB1, SIGLEC1, VSIG4
GO:0034361Very-low-density lipoprotein particle20.03281913APOC1, APOE
GO:0045202Synapse40.041387143C1QB, C1QA, FABP5, C1QC
GO:0034364High-density lipoprotein particle20.042002749APOC1, APOE
Molecular function
GO:0005344Oxygen transporter activity33.12E−04HBA2, HBD
GO:0019825Oxygen binding30.001605913HBA2, HBD
GO:0005044Scavenger receptor activity30.002957949CD163, CD163L1, STAB1
GO:0086089Voltage–gated potassium channel activity involved in atrial cardiac muscle cell action potential repolarization20.006590464KCNJ5, KCNA5
GO:0060228Phosphatidylcholine–sterol O-acyltransferase activator activity20.009869914APOC1, APOE
GO:0031720Haptoglobin binding20.016397412HBA2, HBD
GO:0043177Organic acid binding20.018022768HBA2, HBD
GO:0020037Heme binding30.025666884HBA2, HBD

Significantly enriched GO terms of DEGs.

DEG, Differentially Expressed Gene; GO, Gene Ontology; ID, Identity Document.

Table 4

KEGG IDDescriptionCountP-ValueGenes
hsa04610Complement and coagulation cascades40.001112855C1QB, C1QA, VSIG4, C1QC
hsa05133Pertussis30.014757214C1QB, C1QA, C1QC
hsa05171Coronavirus disease—COVID-1940.018374812C1QB, C1QA, MMP1, C1QC
hsa05150Staphylococcus aureus infection30.022928086C1QB, C1QA, C1QC
hsa05142Chagas disease30.02567277C1QB, C1QA, C1QC
hsa05322Systemic lupus erythematosus30.043532485C1QB, C1QA, C1QC
hsa04936Alcoholic liver disease30.047059029C1QB, C1QA, C1QC

Significantly enriched KEGG terms of DEGs.

KEGG, Kyoto Encyclopedia of Genes and Genomes; ID, Identity Document; DEG, Differentially Expressed Gene.

PPI network construction and hub gene identification

Using STRING tools, we predicted protein interactions among DEGs. The PPI network presented in Figure 3 consists of 34 nodes and 209 edges. Based on the PPI network, we identified 10 genes with the highest connectivity degree (Table 5). The results showed that C1QA was the most outstanding gene with connectivity degree = 24, followed by C1QB (degree = 23), C1QC (degree = 22), CD163 (degree = 22), SIGLEC1 (degree = 21), APOE (degree = 19), MS4A4A (degree = 19), VSIG4 (degree = 18), CCR1 (degree = 18), STAB1 (degree = 18).

Figure 3

Figure 3

Protein-protein interaction (PPI) network of common DEGs among SRAS-CoV-2 and atherosclerosis. In the figure, the circle nodes represent DEGs and edges represent interactions between nodes.

Table 5

Gene symbolGene descriptionDegree
C1QAComplement C1q A chain24
C1QBComplement C1q B chain23
C1QCComplement C1q C chain22
CD163CD163 molecule22
SIGLEC1Sialic acid binding Ig like lectin 121
APOEApolipoprotein E19
MS4A4AMembrane spanning 4-domains A4A19
VSIG4V-set and immunoglobulin domain containing 418
CCR1C-C motif chemokine receptor 118
STAB1Stabilin 118

Top ten hub genes with higher degree of connectivity.

Discussion

Some diseases, thought to be unrelated, share the same biological processes (10). We conducted the DO analysis on the GSE152418 dataset to find the similarity between diseases and COVID-19, and found that COVID-19 was most significantly associated with atherosclerosis among various diseases. Our results suggest that COVID-19 will lead to faster atherosclerosis. Then, we took the intersection of two datasets, GSE152418 and GSE100927, to identify common genes between COVID-19 and atherosclerosis. After obtaining 34 common genes, the GO, pathway, PPI networks were further analyzed.

GO enrichment analysis showed that C1QA, C1QB, C1QC were significantly enriched in synapse disassembly, complement activation, and innate immune response. Complement 1q (C1q) is composed of six subunits, which form a molecule containing 18 polypeptide chains, while C1qA, C1qB, and C1qC genes encode three types of polypeptide chains, A, B, and C of the subunit of C1q, respectively (11). C1q is an important recognition molecule to initiate the classical pathway involved in the complement activation and function, playing a major role in the connection between innate and specific immunity. (12, 13). After identifying the complement binding site on the antibody Fc segment of the IgM or IgG immune complex, the complement cascade will be activated to clear the antigen-antibody complexes (14). Complement proteins specifically locate apoptotic, immature or weak developing synapses in the central nervous system (15). The number of those apoptotic markers in the synapse is equal to the localization of C1q, which promotes synaptic pruning (16). A study found that of 281 patients diagnosed with COVID-19, 21.1% had dementia and 8.9% had mild cognitive impairment (MCI) (17). Moreover, high activation of C1q leads to a large number of synaptic loss which is associated with the development of Alzheimer's disease (18). Then, does the activation of complement system C1q cause cognitive impairment in COVID-19 patients?

KEGG enrichment analysis is the best way to reflect the changes of pathways in organisms. Those results indicate that complement and coagulation cascades change most significantly in atherosclerosis and COVID-19. Macor et al. found positive lung C1q staining which suggests that the classical pathway is important for complement activation which may be triggered by IgG, antibodies widely distributed in patients' lungs (19). In atherosclerosis plaques, C1q activates the classical complement pathway by recognizing oxidized low-density lipoprotein auto-antibodies or directly binding modified lipoprotein and cholesterol crystal (20). Endothelial dysfunction, an important mechanism for the formation and development of atherosclerosis, can be caused by the activation of the complement system can lead to (20). Gao et al. demonstrated that subsequent endothelial dysfunction persisted in COVID-19 survivors even 327 days after diagnosis (6). The activated fragments generated after the activation of the complement system may be closely related to the coagulation and fibrinolytic system and inflammation in COVID-19 patients, so additional studies on the changes in the number of fragments and tissue distribution are needed.

The 10 hub genes selected by PPI were C1QA, C1QB, C1QC, CD163, SIGLEC1, APOE, MS4A4A, VSIG4, CCR1, and STAB1. The C1QA, C1QB, and C1QC genes had the highest degree in the PPI networks. Then v-set and immunoglobulin domain containing 4 (VSIG4) is the receptor of complement component 3 fragments C3b and iC3b, which activates macrophage immunity through C3b/iC3b binding (21). VSIG4 may be involved in lung injury through induction of phagocytosis (22). VSIG4 activate macrophages, through induction of chemokines, promote the migration of inflammatory cells to the lesion area, and participate in the pathogenesis of arteriosclerosis (23). Increased expressions of C1QA, C1QB, C1QC, and VSIG4 all relate to enhanced complement system. CD163, a scavenger receptor, is a major component of inflammation and the immune response. Among plasmacytoid dendritic cells, type I interferon is induced with the appearance of CD163+ SIGLEC1+ macrophages with increased angiotensin converting enzyme 2 (ACE2) levels (24). Macrophages are highly enriched in the lungs of macaques at peak viremia and harbor the SARS-CoV-2 virus while also expressing an interferon-driven innate antiviral gene signature (25). CD163(+) macrophages promote angiogenesis, vascular permeability and inflammation in atherosclerosis via the CD163/HIF1α/VEGF-A pathway. The increased expression of CD163 was revealed in ruptured coronary plaques (26). There are three APOE isoforms, namely APOE epsilon2 (APOE2), APOE epsilon3 (APOE3) and APOE epsilon4 (APOE4) located on chromosome 19q13.2 (27). APOE can function as an endogenous, concentration-dependent pulmonary danger signal that primes and activates the NLPR3 inflammasome in bronchoalveolar lavage fluid macrophages from asthmatic subjects to secrete IL-1β (28). A recent study in the UK Biobank Cohort, APOE4 has been shown to associate with increased susceptibility to SARS-CoV-2 infection and COVID-19 mortality (29). APOE is a therapeutic target for statins that inhibit inflammation in patients with atherosclerotic vascular disease.

Statins possess antiviral, immunomodulatory, antithrombotic, and anti-inflammatory properties, which may improve short- and long-term outcomes in COVID-19 patients.

STAB1 encodes an unusual type of multifunctional scavenger receptor that causes increased lipid uptake and transient lipid depletion in virus-infected areas and is associated with poor prognosis for COVID-19 (30). STAB1 expression may contribute to foam cell formation, monocyte adhesion/migration, and regulation of inflammation in atherosclerotic lesions (31). Lectins such as sialic acid-binding Ig-like lectin 1 (SIGLEC1/CD169) mediate the attachment of viruses to Antigen-presenting cells (APCs) (32). SIGLEC1 expression is induced on APCs upon IFN-α or LPS exposure and increased in myeloid cells of COVID-19 patients (33) Inhibition of Siglec-1 prevents monocytes from adhering to vascular endothelial cells in the early stage of atherosclerosis, and reduces lipid phagocytosis and chemokine secretion of macrophages, alleviating the inflammatory response of established fat streaking lesions (34). CCR1 is critical mediators of monocyte/macrophage polarization and tissue infiltration, which are pathogenic hallmarks of severe COVID-19 (35). The use of monocyte CCR1 in arterial recruitment is due in part to activated chemokines of platelet deposition, which is important in the early stages of atherosclerosis (36). MS4A4A is a novel M2 macrophage cell surface marker, which is essential for dectin-1-dependent activation of NK cell-mediated anti-metastatic properties (37). Silva-Gomes et al. found MS4A4A was expressed by MΦs or alveolar MΦs in COVID-19 bronchoalveolar lavage fluid (38).

Through DO analysis, we also found several neurological disorders associated with COVID-19, such as focal epilepsy, temporal lobe epilepsy, migraine, epilepsy syndrome, neuroblastoma, and lateral sclerosis. There have been a large number of reported cases of these conditions, with a seizure prevalence ranging from 0 to 26% in COVID-19 patients (39, 40). Moreover, seizures may be related to cerebrovascular disease and central nervous system infection. Vascular endothelial injury leads to hypercoagulability and microembolism, resulting in reduced cortical blood flow accompanied by hypoxia. Vascular endothelial dysfunction can lead to changes in the nervous system, resulting in neurological sequelae (41). The Atherosclerosis Risk in Communities (ARIC) study also revealed that migraine patients were more susceptible to retinopathy (retinal hemorrhage, macular oedema, retinal microvascular abnormalities, venous bleeding, etc.) than non-migraine patients, and retinopathy was more strongly associated with migraine in people without a history of diabetes or hypertension (42). Interestingly, we also discovered DEGs enrichment in retinopathy. Besides, previous animal-based experimental studies of the coronavirus infection reported retinal diseases such as retinal vasculitis and retinal degeneration. Moreover, blood-retinal barrier breakdown revealed the possibility of immune-privileged site infectivity by SARS-CoV-2 (43). We believe that SARS-CoV-2 causes vascular injury and may lead to retinal degeneration. Results also revealed different types of cancer, such as pancreatic cancer, kidney cancer, breast carcinoma, stomach cancer, esophageal cancer, and ovarian cancer. In patients with COVID19, severe illness and mortality are closely related to cancer. SARS CoV 2 may promote tumor progression and stimulate metabolic switching in tumor cells to initiate tumor metabolic modes with higher production efficiency, such as glycolysis, for facilitating the replication of SARS CoV 2 (44). Meanwhile, we also established that muscular disease, such as myositis, is associated with COVID-19. Previous studies have demonstrated that patients with dermatomyositis have three immunogenic linear epitopes with a high degree of sequence identity to the SARS-CoV-2 protein, so potential exposure to the coronavirus family may lead to the development of dermatomyositis (45). Effective Janus kinase (JAK) inhibitors for dermatomyositis, including tofacitinib, ruxolitinib, and baricitinib, may provide new directions for COVID-19 treatment.

In conclusion, the study provides new insights for the common pathogenesis of COVID-19 and atherosclerosis by looking for common transcriptional features. The DEGs identified by bioinformatics data analysis, including C1QA, C1QB, C1QC, CD163, SIGLEC1, APOE, MS4A4A, VSIG4, CCR1, and STAB1, may be therapeutic targets for the atherosclerosis caused by COVID-19. However, more wet lab-based studies are required to validate the impact of COVID-19 severity on atherosclerosis. Studies on the long-term effects of SARS-CoV-2 infection, the effect of persistent endothelial dysfunction on atherosclerosis, and the role of preventive therapy are also needed.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

Author contributions

JZ performed the data analyses and wrote the manuscript. LZ helped perform the analysis with constructive discussions. Both authors approved the final version of the manuscript.

Acknowledgments

We acknowledge GEO database for providing their platforms and contributors for uploading their meaningful datasets.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    WHO COVID-19 Dashboard (2020). Geneva: World Health Organization. Available online: https://covid19.who.int/ (accessed March 9, 2022).

  • 2.

    CallawayE. The race for coronavirus vaccines: a graphical guide. Nature. (2020) 580:5767. 10.1038/d41586-020-01221-y

  • 3.

    ZhengYXuHYangMZengYChenHLiuRet al. Epidemiological characteristics and clinical features of 32 critical and 67 noncritical cases of COVID-19 in Chengdu. J Clin Virol. (2020) 127:104366. 10.1016/j.jcv.2020.104366

  • 4.

    HuangCWangYLiXRenLZhaoJHuYet al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet (London, England). (2020) 395:497506. 10.1016/S0140-6736(20)30183-5

  • 5.

    GaoYPZhouWHuangPNLiuHYBiXJZhuYet al. Persistent Endothelial dysfunction in coronavirus diseasE−2019 survivors late after recovery. Front. Med. (2022) 9:809033. 10.3389/fmed.2022.809033

  • 6.

    RobsonA. Preventing cardiac damage in patients with COVID-19. Nat Rev Cardiol. (2021) 18:387. 10.1038/s41569-021-00550-3

  • 7.

    MesterABenedekIRatNTolescuCPolexaSABenedekTet al. Imaging cardiovascular inflammation in the COVID-19 era. Diagnostics (Basel, Switzerland). (2021) 11:1114. 10.3390/diagnostics11061114

  • 8.

    YuGWangLGYanGRHeQY. DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis. Bioinformatics (Oxford, England). (2015) 31:6089. 10.1093/bioinformatics/btu684

  • 9.

    WuTHuEXuSChenMGuoPDaiZet al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (New York. NY). (2021) 2:100141. 10.1016/j.xinn.2021.100141

  • 10.

    ButteAJKohaneIS. Creation and implications of a phenomE–genome network. Nat Biotechnol. (2006) 24:5562. 10.1038/nbt1150

  • 11.

    LuJHTehBKWangLDWangYNTanYSLaiMCet al. The classical and regulatory functions of C1q in immunity and autoimmunity. Cell Mol Immunol. (2008) 5:921. 10.1038/cmi.2008.2

  • 12.

    TennerAJ. C1q interactions with cell surface receptors. Behring Institute Mitteilungen. (1989) 84:2209.

  • 13.

    MortensenSASanderBJensenRKPedersenJSGolasMMJenseniusJCet al. Structure and activation of C1, the complex initiating the classical pathway of the complement cascade. Proc Natl Acad Sci USA. (2017) 114:98691. 10.1073/pnas.1616998114

  • 14.

    SjöbergAPTrouwLABlomAM. Complement activation and inhibition: a delicate balance. Trends Immunol. (2009) 30:8390. 10.1016/j.it.2008.11.003

  • 15.

    CornellJSalinasSHuangHYZhouM. Microglia regulation of synaptic plasticity and learning and memory. Neural Regener Res. (2022) 17:70516. 10.4103/1673-5374.322423

  • 16.

    GyörffyBAKunJTörökGBulyákiÉBorhegyiZGulyássyPet al. Local apoptotic-like mechanisms underlie complement-mediated synaptic pruning. Proc Natl Acad Sci USA. (2018) 115:63038. 10.1073/pnas.1722613115

  • 17.

    Martín-JiménezPMuñoz-GarcíaMISeoaneDRoca-RodríguezLGarcía-ReyneALaluezaAet al. Cognitive impairment is a common comorbidity in deceased covid-19 patients: a hospital-based retrospective cohort study. J Alzheimer's Dis. (2020) 78:136772. 10.3233/JAD-200937

  • 18.

    StephanAHBarresBAStevensB. The complement system: an unexpected role in synaptic pruning during development and disease. Annu Rev Neurosci. (2012) 35:36989. 10.1146/annurev-neuro-061010-113810

  • 19.

    MacorPDuriguttoPMangognaABussaniRDe MasoLD'ErricoPLet al. MultiplE–organ complement deposition on vascular endothelium in COVID-19 patients. Biomedicines. (2021) 9:1003. 10.3390/biomedicines9081003

  • 20.

    SamstadEONiyonzimaNNymoSAuneMHRyanLBakkeSSet al. Cholesterol crystals induce complement-dependent inflammasome activation and cytokine release. J Immunol (Baltimore, Md. : 1950). (2014) 192:283745. 10.4049/jimmunol.1302484

  • 21.

    HertleEStehouwerCDvan GreevenbroekMM. The complement system in human cardiometabolic disease. Mol Immunol. (2014) 61:13548. 10.1016/j.molimm.2014.06.031

  • 22.

    HelmyKYKatschkeKJGorganiNNKljavinNMElliottJMDiehlLet al. CRIg: a macrophage complement receptor required for phagocytosis of circulating pathogens. Cell. (2006) 124:91527. 10.1016/j.cell.2005.12.039

  • 23.

    ZhaoCMoJZhengXWuZLiQFengJet al. Identification of an alveolar macrophagE–related core gene set in acute respiratory distress syndrome. J Inflamm Res. (2021) 14:235361. 10.2147/JIR.S306136

  • 24.

    LeeMYKimWJKangYJJungYMKangYMSukKet al. Z39Ig is expressed on macrophages and may mediate inflammatory reactions in arthritis and atherosclerosis. J Leukoc Biol. (2006) 80:9228. 10.1189/jlb.0306160

  • 25.

    SinghDKAladyevaEDasSSinghBEsaulovaESwainAet al. Myeloid cell interferon responses correlate with clearance of SARS-CoV-2. Nat Commun. (2022) 13:679. 10.1038/s41467-022-28315-7

  • 26.

    GuoLAkahoriHHarariESmithSLPolavarapuRKarmaliVet al. CD163+ macrophages promote angiogenesis and vascular permeability accompanied by inflammation in atherosclerosis. J Clin Invest. (2018) 128:110624. 10.1172/JCI93025

  • 27.

    WeisgraberKHRallSCMahleyRW. Human E apoprotein heterogeneity. CysteinE–arginine interchanges in the amino acid sequence of the apo-E isoforms. J Biol Chem. (1981) 256:907783. 10.1016/S0021-9258(19)52510-8

  • 28.

    GordonEMYaoXXuHKarkowskyWKalerMKalchiem-DekelOet al. Apolipoprotein E is a concentration-dependent pulmonary danger signal that activates the NLRP3 inflammasome and IL-1β secretion by bronchoalveolar fluid macrophages from asthmatic subjects. J Allerg Clin Immunol. (2019) 144:426441.e3. 10.1016/j.jaci.2019.02.027

  • 29.

    KuoCLPillingLCAtkinsJLMasoliJDelgadoJKuchelGAet al. ApoE e4e4 genotype and mortality with COVID-19 in UK Biobank. J Gerontol Ser A, Biol Sci Med Sci. (2020) 75:18013. 10.1093/gerona/glaa169

  • 30.

    VlasovIPanteleevaAUsenkoTNikolaevMIzumchenkoAGavrilovaEet al. Transcriptomic profiles reveal downregulation of low-density lipoprotein particle receptor pathway activity in patients surviving severe COVID-19. Cells. (2021) 10:3495. 10.3390/cells10123495

  • 31.

    BrochériouIMaoucheSDurandHBraunersreutherVLe NaourGGratchevAet al. Antagonistic regulation of macrophage phenotype by M-CSF and GM-CSF: implication in atherosclerosis. Atherosclerosis. (2011) 214:31624. 10.1016/j.atherosclerosis.2010.11.023

  • 32.

    Perez-ZsoltDMuñoz-BasagoitiJRodonJElosua-BayesMRaïch-ReguéDRiscoCet al. SARS-CoV-2 interaction with Siglec-1 mediates trans-infection by dendritic cells. Cell Mol Immunol. (2021) 18:26768. 10.1038/s41423-021-00794-6

  • 33.

    BedinASMakinsonAPicotMCMennechetFMalergueFPisoniAet al. Monocyte CD169 expression as a biomarker in the early diagnosis of coronavirus disease 2019. J Infect Dis. (2021) 223:5627. 10.1093/infdis/jiaa724

  • 34.

    XiongYSWuALMuDYuJZengPSunYet al. Inhibition of siglec-1 by lentivirus mediated small interfering RNA attenuates atherogenesis in apoE–deficient mice. Clin Immunol (Orlando, Fla). (2017) 174:3240. 10.1016/j.clim.2016.11.005

  • 35.

    StikkerBStikGHendriksRWStadhoudersR. Severe COVID-19 associated variants linked to chemokine receptor gene control in monocytes and macrophages. bioRxiv : the preprint server for biology. (2021) 01, 22.427813. 10.1101/2021.01.22.427813

  • 36.

    DrechslerMMegensRTvan ZandvoortMWeberCSoehnleinO. Hyperlipidemia-triggered neutrophilia promotes early atherosclerosis. Circulation. (2010) 122:183745. 10.1161/CIRCULATIONAHA.110.961714

  • 37.

    HuoQLiZChenSWangJLiJXieNet al. VWCE as a potential biomarker associated with immune infiltrates in breast cancer. Cancer Cell Int. (2021) 21:272. 10.1186/s12935-021-01955-3

  • 38.

    Silva-GomesRMapelliSNBoutetMAMattiolaISironiMGrizziFet al. Differential expression and regulation of MS4A family members in myeloid cells in physiological and pathological conditions. J leukocyte Biol. (2021) 111:81736. 10.1002/JLB.2A0421-200R

  • 39.

    PellinenJHolmesMG. Evaluation and treatment of seizures and epilepsy during the COVID-19 pandemic. Curr Neurol Neurosci Rep. (2022) 17. 10.1007/s11910-022-01174-x

  • 40.

    DanounOAZillgittAHillCZutshiDHarrisDOsmanGet al. Outcomes of seizures, status epilepticus, and EEG findings in critically ill patient with COVID-19. Epilepsy Behav: EandB. (2021) 118:107923. 10.1016/j.yebeh.2021.107923

  • 41.

    QinYWuJChenTLiJZhangGWuDet al. Long-term microstructure and cerebral blood flow changes in patients recovered from COVID-19 without neurological manifestations. J Clin Invest. (2021) 131:e147329. 10.1172/JCI147329

  • 42.

    RoseKMWongTYCarsonAPCouperDJKleinRSharrettARet al. Migraine and retinal microvascular abnormalities: the atherosclerosis risk in communities study. Neurology. (2007) 68:1694700. 10.1212/01.wnl.0000261916.42871.05

  • 43.

    ZhouXZhouYNAliALiangCYeZChenXet al. Case report: a rE–positive case of SARS-CoV-2 Associated With Glaucoma. Front Immunol. (2021) 12:701295. 10.3389/fimmu.2021.701295

  • 44.

    LiYSRenHCCaoJH. Correlation of SARS-CoV-2 to cancer: Carcinogenic or anticancer? (Review). Int J Oncol. (2022) 60:42. 10.3892/ijo.2022.5332

  • 45.

    MegremisSWalkerTHeXOllierWChinoyHHampsonLet al. Antibodies against immunogenic epitopes with high sequence identity to SARS-CoV-2 in patients with autoimmune dermatomyositis. Ann Rheum Dis. (2020) 79:13836. 10.1136/annrheumdis-2020-217522

Summary

Keywords

COVID-19, atherosclerosis, C1q, SARS-CoV-2, immune

Citation

Zhang J and Zhang L (2022) Bioinformatics approach to identify the influences of SARS-COV2 infections on atherosclerosis. Front. Cardiovasc. Med. 9:907665. doi: 10.3389/fcvm.2022.907665

Received

30 March 2022

Accepted

11 July 2022

Published

18 August 2022

Volume

9 - 2022

Edited by

Mark Slevin, Manchester Metropolitan University, United Kingdom

Reviewed by

Sakir Ahmed, Kalinga Institute of Medical Sciences (KIMS), India; Erkan Cüre, Bagcilar Medilife Hospital, Turkey

Updates

Copyright

*Correspondence: Liming Zhang

This article was submitted to Atherosclerosis and Vascular Medicine, a section of the journal Frontiers in Cardiovascular Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics