ORIGINAL RESEARCH article

Front. Genet., 25 May 2022

Sec. Computational Genomics

Volume 13 - 2022 | https://doi.org/10.3389/fgene.2022.889348

Screening the Potential Biomarkers of COVID-19-Related Thrombosis Through Bioinformatics Analysis

  • 1. Department of Emergency, First Medical Center of Chinese PLA General Hospital, Beijing, China

  • 2. Department of Nephrology, First Medical Center of Chinese PLA General Hospital, Beijing, China

Article metrics

View details

4

Citations

3,3k

Views

1,2k

Downloads

Abstract

A high proportion of critically ill patients with coronavirus disease 2019 (COVID-19) experience thrombosis, and there is a strong correlation between anticoagulant therapy and the COVID-19 survival rate, indicating that common COVID-19 and thrombosis targets have potential therapeutic value for severe COVID-19.Gene expression profiling data were downloaded from Gene Expression Omnibus (GEO), and common differentially expressed genes (co-DEGs) were identified. The potential biological functions of these co-DEGs were explored by functional enrichment analysis, and protein–protein interaction (PPI) networks were constructed to elucidate the molecular mechanisms of the co-DEGs. Finally, hub genes in the co-DEG network were identified, and correlation analysis was performed.We identified 8320 upregulated genes and 7651 downregulated genes from blood samples of COVID-19 patients and 368 upregulated genes and 240 downregulated genes from blood samples of thrombosis patients. The enriched cellular component terms were mainly related to cytosolic ribosomes and ribosomal subunits. The enriched molecular function terms were mainly related to structural constituents of ribosomes and electron transfer activity. Construction of the PPI network and identification of hub genes ultimately confirmed that RPS7, IGF1R, DICER1, ERH, MCTS1, and TNPO1 were jointly upregulated hub genes, and FLNA and PXN were jointly downregulated hub genes.The identification of novel potential biomarkers provides new options for treating COVID-19-related thrombosis and reducing the rate of severe COVID-19.

Introduction

The coronavirus disease 2019 (COVID-19) pandemic has had a significant impact on global public health, social interactions and the economy (Han et al., 2020; Global Burden of Disease 2020 Health Financing Collaborator Network, 2021; Lal et al., 2021; Mcintyre et al., 2022). As of February 2022, nearly 400 million people have been diagnosed with COVID-19 worldwide, and more than 5.7 million people have died. COVID-19 can induce different complications and affect different organs (Drake et al., 2021; Wang et al., 2021), mainly the respiratory tract, heart, gastrointestinal tract (Serban et al., 2021; Bulte et al., 2022), central nervous system (Sriwastava et al., 2021), liver, and kidney (Gross et al., 2020). Elderly individuals and patients with underlying diseases are prone to severe cases, and their mortality rates are significantly higher than those of other populations (Ghasemiyeh et al., 2021). Patients who die from COVID-19 infection in the late stage of disease progression can develop acute respiratory distress syndrome (ARDS) (Zangrillo et al., 2022), cytokine storm syndrome (CSS) (Yongzhi, 2021), acute kidney injury (AKI) (Ronco et al., 2020), septic shock (Wei et al., 2021), multiple organ failure (MOF) (Devaux et al., 2020) and cardiac arrest (Shao et al., 2020), resulting in death. Especially in patients with severe COVID-19, cardiac manifestations such as arrhythmia, myocardial infarction, acute heart failure and cardiogenic shock are common and are associated with increased mortality (Bangalore et al., 2020; Lala et al., 2020). Researchers have proposed various hypotheses about the direct and indirect mechanisms of COVID-19-related heart damage, but the exact mechanism is still unclear (Chung et al., 2021). The earlier theory that COVID-19 heart damage is due to acute viral myocarditis has been disproved by the lack of cases meeting the histological criteria for myocarditis (Halushka and Vander Heide, 2021). Subsequently, further studies have suggested that coagulation abnormalities, including thrombotic microangiopathy and complications of venous and arterial thromboembolism, are common in critically ill patients with COVID-19 (Levi et al., 2020). The pathological results suggest that the most common pathological cause of myocardial cell necrosis is the occurrence of microthrombi (Pellegrini et al., 2021). Moreover, a series of studies have also found that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection promotes hypercoagulability, and there is a strong correlation between the use of anticoagulants and COVID-19 survival (Nadkarni et al., 2020; Fröhlich et al., 2021; Lawler et al., 2021). A series of studies have confirmed that COVID-19 is associated with venous thromboembolic events, including pulmonary embolism, lower extremity deep venous thrombosis and cerebral venous sinus thrombosis, and there is a correlation between thrombosis and mortality in patients with COVID-19 (Rosovsky et al., 2020). Some studies have even found a 51.5-fold increased risk of in-hospital mortality in patients with D-2mer ≥2.0 μg/ml (Zhang L. et al., 2020), but the actual prevalence and characteristics of venous thromboembolism (VTE) remain unclear (Sastry et al., 2022). The incidence of VTE reported by different studies varies significantly (ranging from 0.7% to 57%) due to differences in disease severity and diagnostic imaging techniques used (Cui et al., 2020; Klok et al., 2020; Lodigiani et al., 2020; Suh et al., 2021). Even with anticoagulant prophylaxis, the incidence of VTE in COVID-19 patients remains high and is much higher than that in non-COVID patients receiving thromboprophylaxis (Lodigiani et al., 2020). Cytokine storms were initially believed to be the main factor for COVID-19 thrombosis (Lowenstein and Solomon, 2020), antiphospholipid antibodies were later detected in COVID-19 (Zhang Y. et al., 2020; Zuo et al., 2020), and the presence of prothrombotic antibodies also constituted one of the risk factors. Studies have also suggested that elevated levels of neutrophil extracellular traps are closely related to thrombotic events (Zuo et al., 2021), but their overall role in COVID-19-related thrombosis remains unclear (Bonaventura et al., 2021). At present, most reports on thrombosis in patients with COVID-19 are case reports (Achim et al., 2022; Bleiziffer et al., 2022; Jeckelmann et al., 2022; Russell et al., 2022; Takasu et al., 2022) and some observational studies (Maier et al., 2020). Even after strict ethical review, only a small number of large studies have confirmed the relationship between the use of anticoagulants and the survival rate of COVID-19 (Lawler et al., 2021). There is a good correlation. Therefore, it is still necessary to obtain a theoretical basis for the mechanism of COVID-19 thrombosis. In summary, exploring targets that can be used to treat both COVID-19 and thrombosis has theoretical guiding significance for the assessment and treatment of COVID-19 patients. However, there are various difficulties in clinical research on COVID-19, including ethical requirements and limited patient contact. The use of public databases solves the above problems. Further analysis of relevant datasets from public databases we can be used to identify and develop reliable biomarkers and treatment strategies for thrombosis in patients with COVID-19 and thus to better understand the mechanism of thrombosis in patients with COVID-19 and identify the corresponding targets. Such studies will promote follow-up drug research and development. We identified 8320 upregulated genes and 7651 downregulated genes from blood samples of COVID-19 patients and 368 upregulated genes and 240 downregulated genes from blood samples of thrombosis patients. We identified hub genes that may contribute to the evaluation and treatment of COVID-19-related thrombosis.

Materials and Methods

Data Collection and Processing

The GEO database was searched with “COVID-19” and “thrombosis” as keywords, and we searched using keywords and restricted the screening conditions to “Homo sapiens” and “Expression profiling by array” and limited the sample to peripheral blood. After reading the original studies on the data sets that met the screening conditions individually, we focused on the experimental design and methodology. The data set that met the principle of randomized control and did not violate ethical guidelines was taken as the research object. After screening each dataset, GSE19151 and GSE164805 were ultimately selected for further research. The GSE19151 dataset included 70 adults who had used warfarin one or more times for venous thromboembolism and 63 healthy controls. This dataset excluded patients with antiphospholipid syndrome and cancer and was generated using the GPL571 [HG-U133A_2] Affymetrix Human Genome U133A 2.0 Array. The GSE164805 dataset included 10 patients diagnosed with COVID-19 and 5 healthy controls and was generated using the GPL26963 Agilent-085982 Arraystar Human LncRNA V5 microarray. R software (version 4.0.2, http://r-project.org/) was used to download the two datasets and relevant information. The raw data were processed through the “affy” (Gautier et al., 2004) package, and the RMA algorithm was used to perform background correction and data normalization. The normalized results were visualized using boxplots (Supplementary Figures S1, S2).

Identification of DEGs

The “limma” package of R software was used to analyze the disease group and healthy control group of the two datasets, and the differentially expressed genes (DEGs) between the two groups were obtained. To identify the most meaningful DEGs, the following limits were used to determine which DEGs would be retained: (1) | log2 (fold-change) | > 0.5 and (2) adjusted p < 0.05. In addition, adjusted p < 0.05 was considered to indicate a significant difference (Ritchie et al., 2015; Shankar et al., 2020).

GO and KEGG Pathway Enrichment Analysis of the DEGs

Gene Ontology (GO) functional analysis (including cellular component [CC], biological process [BP], and molecular function [MF] analysis) is a powerful bioinformatics tool for classifying gene expression and associated properties (Ashburner et al., 2000). Kyoto Encylopedia of Genes and Genomes (KEGG) pathway analysis was used to determine which cellular pathways may be related to the DEGs (Altermann and Klaenhammer, 2005). The “clusterProfiler" (Yu et al., 2012) package of R software was used for enrichment analysis of the DEGs, the “org.Hs.eg.db” package was used for ID conversion, and the “Goplot” package was used to calculate Z scores (Walter et al., 2015). The results are visually presented as chord graphs and tables.

Identification of Co-DEGs and Construction of the PPI Network

The co-DEGs of the two datasets were analyzed using a Venn diagram, and the results were visualized using the “ggplot2” (Rosenthal et al., 2019) package of R software. The co-upregulated genes and co-downregulated genes were subsequently uploaded to the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, http://string-db.org/) database for protein–protein interaction (PPI) network construction, and 0.4 was used as the minimum required interaction score for analysis. The results were input into Cytoscape (Shannon et al., 2003) to construct and visualize the molecular interaction network.

Identification of Hub Genes

Hub genes are the DEGs that show the highest association with other network components, and they are important nodes in the PPI network. The Cyto-Hubba (Chin et al., 2014) plug-in in Cytoscape software was used to determine the hub genes in the PPI network. To reduce the selection bias caused by different algorithms, the MCC, DMNC, MNC, and EPC algorithms were used to determine the hub genes, and the final results were confirmed and visualized with a Venn diagram.

Functional Enrichment Analysis of Hub Genes

The selected hub genes were analyzed by the PANTHER(Mi et al., 2019) classification system (http://pantherdb.org/), and the basic classification of each gene was obtained. Then, the Human Protein Atlas (Uhlen et al., 2019) (HPA) (https://www.proteinatlas.org/) was used to explore the expression profile of each hub gene in human immune cells. Subsequently, the Deeply Integrated human Single-Cell Omics (DISCO) database (Li et al., 2022) (https://www.immunesinglecell.org/) was used to clarify the expression level of the hub genes in all blood cells of COVID-19 patients. Finally, the “clusterProfiler” (Yu et al., 2012) package of R software was used for the enrichment analysis of the hub genes, and the “ggplot2” package of R software was used to visually present the results of the enrichment analysis. The results are displayed in a histogram and table.

Results

Identification of DEGs

The “limma” package of R software was used to identify the DEGs of the two datasets. In the GSE19151 dataset, included 70 adults who had used warfarin one or more times for venous thromboembolism and 63 healthy controls, the number of genes after data processing was 21144, including 368 upregulated genes that met the log2 (FC) > 0.5 and adjusted p < 0.05 criteria and 240 downregulated genes that met the log2 (FC) < −0.5 and adjusted p < 0.05 criteria. In the GSE164805 dataset, included 10 patients diagnosed with COVID-19 and 5 healthy controls, the number of genes after data processing was 60494, including 8320 upregulated genes meeting the log2 (FC) > 0.5 and adjusted p < 0.05 criteria and 7651 downregulate genes meeting the log2 (FC) < −0.5 and adjusted p < 0.05 criteria. A volcano plot was generated using the “ggplot2” package to visualize the DEGs (Figure 1).

FIGURE 1

FIGURE 1

DEGs in the two datasets. (A) Volcano plot of GSE19151; Red represents upregulated genes, blue represents downregulated genes, and gray represents no significantly expressed genes. (B) Volcano plot of GSE164805. The criteria for statistically significant difference of DEGs was adjusted |log2(FC)|>0.5 and p.adj<0.05 in expression.

GO and KEGG Pathway Enrichment Analyses of DEGs

The main enriched BP terms of the DEGs in the venous thrombosis dataset were SRP-dependent cotranslational protein targeting to membrane (GO:0006614) and cotranslational protein targeting to membrane (GO:0006613). The main enriched CC terms were cytosolic ribosome (GO:0022626) and ribosomal subunit (GO:0044391), and the main enriched MF terms were structural constituent of ribosome (GO:0003735) and electron transfer activity (GO:0009055). The KEGG pathway analysis mainly showed enrichment of the terms ribosome (hsa03010) and cardiac muscle contraction (hsa04260). For the DEGs identified from the COVID-19 dataset, the enriched BP terms mainly included nuclear-transcribed mRNA catabolic process, nonsense-mediated decay (GO:0000184) and translational initiation (GO:0006413), the enriched CC terms mainly included cytosolic ribosome (GO:0022626) and ribosomal subunit (GO:0044391), and the enriched MF term was structural constituent of ribosome (GO:0003735). The KEGG pathway enrichment analysis mainly identified the terms leishmaniasis (hsa05140) and graft-versus-host disease as enriched (hsa05332). For both datasets, the CC and MF terms were similar. All the results are presented in chord diagrams (Figure 2) and tables (Tables 1, 2).

FIGURE 2

FIGURE 2

GO terms and KEGG pathway enrichment. (A) Details of GO terms and KEGG pathway enrichment in GSE19151; The left side of the figure indicated that the genes were sorted according to logFC. Red represented up-regulation, blue represented down-regulation, and color depth represented change multiples. The right side of the figure indicates the order of enrichment analysis term according to the change of strength. There is a link between genes and term, indicating that they have a subordinate relationship. (B) Details of GO terms and KEGG pathway enrichment in GSE164805.

TABLE 1

OntologyIDDescriptionGeneRatioBgRatiop valuep.adjustq value
BPGO:0006614SRP-dependent cotranslational protein targeting to membrane12/99105/186703.24e-135.61e-104.88e-10
BPGO:0006613cotranslational protein targeting to membrane12/99109/186705.11e-135.61e-104.88e-10
CCGO:0022626cytosolic ribosome14/102112/197176.19e-161.70e-131.52e-13
CCGO:0044391ribosomal subunit15/102190/197175.51e-147.55e-126.76e-12
MFGO:0003735structural constituent of ribosome15/98202/176973.43e-131.02e-109.46e-11
MFGO:0009055electron transfer activity7/98114/176973.26e-064.83e-044.50e-04
KEGGhsa03010Ribosome14/56158/80761.93e-123.16e-102.88e-10
KEGGhsa04260Cardiac muscle contraction6/5687/80762.77e-050.0020.002

Details of GO terms and KEGG pathway enrichment in GSE19151.

TABLE 2

OntologyIDDescriptionGeneRatioBgRatiopvaluep.adjustq value
BPGO:0000184nuclear-transcribed mRNA catabolic process, nonsense-mediated decay30/1653120/186701.28e-077.54e-047.17e-04
BPGO:0006413translational initiation39/1653193/186708.02e-070.0020.002
CCGO:0022626cytosolic ribosome31/1734112/197175.12e-093.58e-063.31e-06
CCGO:0044391ribosomal subunit42/1734190/197171.80e-086.29e-065.83e-06
MFGO:0003735structural constituent of ribosome37/1638202/176974.24e-050.0450.044
KEGGhsa05140Leishmaniasis24/77477/80761.00e-073.15e-052.69e-05
KEGGhsa05332Graft-versus-host disease14/77442/80762.01e-050.0030.003

Details of GO terms and KEGG pathway enrichment in GSE164805.

PPI Network Construction and Hub Gene Identification

The “VennDiagram” package of R software was used to extract the co-DEGs of the two datasets (Figure 3C). A total of 103 co-upregulated genes and 106 co-downregulated genes were identified. The co-upregulated genes and co-downregulated genes were uploaded to the STRING database to construct a PPI network (Figures 3A,B), and the results were imported into Cytoscape software for further analysis. The PPI network of the co-upregulated genes included 59 nodes and 57 edges, and the PPI network of the co-downregulated genes included 82 nodes and 126 edges (Figure 4). Next, the Cyto-Hubba plugin was used to identify the hub genes. To reduce the selection bias caused by different algorithms, the MCC, DMNC, MNC, and EPC algorithms were used to extract the hub genes, and the “VennDiagram” package was used to merge the hub genes extracted by the four algorithms (Figure 5). Ultimately, RPS7, IGF1R, DICER1, ERH, MCTS1, and TNPO1 were identified as jointly upregulated hub genes, and FLNA and PXN were identified as jointly downregulated hub genes.

FIGURE 3

FIGURE 3

PPI network and Venn diagram. (A) A PPI network were constructed by co-upregulated genes. (B) A PPI network were constructed by co-downregulated genes. (C) The Venn diagram showed the co-DEGs of the two datasets.

FIGURE 4

FIGURE 4

PPI network of DEGs. (A) The PPI network of the co-upregulated genes included 59 nodes and 57 edges. (B) The PPI network of the co-downregulated genes included 82 nodes and 126 edges.

FIGURE 5

FIGURE 5

The hub genes were merged by the MCC, DMNC, MNC, and EPC algorithms. (A) RPS7, IGF1R, DICER1, ERH, MCTS1, and TNPO1 were jointly upregulated hub genes. (B) FLNA and PXN were identified as jointly downregulated hub genes.

Analysis of Hub Genes

The PANTHER (Mi et al., 2019) classification system was used to analyze the selected hub genes, and the basic classification of each gene was obtained. The results are listed in a table (Table 3). Complex interactions between innate immunity, platelet activation and coagulation lead to immune thrombosis. Immune thrombosis is an important host defense mechanism that limits the systematic transmission of pathogens through blood flow. However, abnormal activation of immune thrombosis can also lead to myocardial infarction, stroke and venous thromboembolism (Stark and Massberg, 2021). To explore the expression of the hub genes in human immune cells, the HPA database was used to analyze each hub gene (Figure 6). HPA RNA-seq data are displayed in different colors based on cell type, including B cells, T cells, natural killer (NK) cells, monocytes, granulocytes, dendritic cells and total peripheral blood mononuclear cells (PBMCs). The DISCO (Li et al., 2022) database was used annotate the single-cell sequencing data of blood cells from patients with COVID-19. The expression level of each hub gene in each cell is shown. The deeper the color is, the higher the expression of the gene in the cell (Figure 7).

TABLE 3

IDGene nameGene IDPANTHER Family/SubfamilyPANTHER Protein Class
MCTS1Malignant T-cell-amplified sequence 1HUMAN|HGNC = 23357|UniProtKB = Q9ULC4MALIGNANT T-CELL-AMPLIFIED SEQUENCE 1 (PTHR22798:SF12)translation initiation factor
DICER1Endoribonuclease DicerHUMAN|HGNC = 17098|UniProtKB = Q9UPY3ENDORIBONUCLEASE DICER (PTHR14950:SF37)
ERHEnhancer of rudimentary homologHUMAN|HGNC = 3447|UniProtKB = P84090ENHANCER OF RUDIMENTARY HOMOLOG (PTHR12373:SF9)DNA-binding transcription factor
TNPO1Transportin-1HUMAN|HGNC = 6401|UniProtKB = Q92973TRANSPORTIN-1 (PTHR10527:SF21)transporter
RPS740S ribosomal protein S7HUMAN|HGNC = 10440|UniProtKB = P6208140S RIBOSOMAL PROTEIN S7 (PTHR11278:SF5)ribosomal protein
IGF1RInsulin-like growth factor 1 receptorHUMAN|HGNC = 5465|UniProtKB = P08069INSULIN-LIKE GROWTH FACTOR 1 RECEPTOR (PTHR24416:SF106)transmembrane signal receptor
PXNPaxillinHUMAN|HGNC = 9718|UniProtKB = P49023PAXILLIN (PTHR24216:SF11)actin or actin-binding cytoskeletal protein
FLNAFilamin-AHUMAN|HGNC = 3754|UniProtKB = P21333FILAMIN-A (PTHR38537:SF5)

The PANTHER classification system of Hub genes.

FIGURE 6

FIGURE 6

Blood cell type expression overview shows RNA-seq data from internally generated Human Protein Atlas (HPA) data. Color-coding is based on blood cell type lineages. The resulting transcript expression values calculated as nTPM, resulting from the internal normalization pipeline for 18 immune cell types and total PBMC. (A) DICER1 (B) ERH (C) RPS7 (D) IGF1R (E) MCTS1 (F) TNP01 (G) PXN (H) FLNA.

FIGURE 7

FIGURE 7

The DISCO database was used annotate the single-cell sequencing data of blood cells from patients with COVID-19. The deeper the color is, the higher the expression of the gene in the cell. (A) DICER1 (B) ERH (C) RPS7 (D) IGF1R (E) MCTS1 (F) TNP01 (G) PXN (H) FLNA.

Functional Enrichment Analysis of Hub Genes

Functional enrichment analysis of the six coupregulated hub genes (RPS7, IGF1R, DICER1, ERH, MCTS1, and TNPO1) was performed. The results showed that the main enriched BP terms were the ribonucleoprotein complex assembly (GO: 0022618), viral gene expression (GO: 0019080) and translational initiation (GO: 0006413). The main enriched CC terms were the cytosolic part (GO: 0044445) and ribosomal subunit (GO: 0044391), and the main enriched MF term was peptide binding (GO: 0042277). The KEGG pathway analysis mainly showed enrichment of the terms adherens junction (hsa04520), melanoma (hsa05218) and longevity regulating pathway - multiple species (hsa04213). All results are displayed in a bubble diagram (Figure 8) and table (Table 4).

FIGURE 8

FIGURE 8

Functional enrichment analysis of the six co-upregulated hub genes. (A) Shows the results of biological process terms enriched by BP analysis; (B) Shows the results of biological process terms enriched by CC analysis; (C) Shows the results of biological process terms enriched by MF analysis; (D) Shows the enriched pathway by KEGG analysis. The coloured dots represent the p-value for that term, with red representing greater significance. The size of the dots represents the number of involved genes.

TABLE 4

OntologyIDDescriptionGeneRatioBgRatiop valuep.adjustq value
BPGO:0035148tube formation2/6148/186709.17e-040.0630.029
BPGO:0019080viral gene expression2/6191/186700.0020.0630.029
BPGO:0006413translational initiation2/6193/186700.0020.0630.029
BPGO:0022618ribonucleoprotein complex assembly2/6277/186700.0030.0630.029
BPGO:0075522IRES-dependent viral translational initiation1/610/186700.0030.0630.029
CCGO:0022627cytosolic small ribosomal subunit2/644/197177.26e-050.0027.64e-04
CCGO:0015935small ribosomal subunit2/672/197171.95e-040.0020.001
CCGO:0022626cytosolic ribosome2/6112/197174.73e-040.0040.002
CCGO:0044391ribosomal subunit2/6190/197170.0010.0080.004
CCGO:0044445cytosolic part2/6247/197170.0020.0110.005
MFGO:0035197siRNA binding1/611/176970.0040.0290.005
MFGO:0043560insulin receptor substrate binding1/611/176970.0040.0290.005
MFGO:0042277peptide binding2/6295/176970.0040.0290.005
MFGO:0031994insulin-like growth factor I binding1/612/176970.0040.0290.005
MFGO:0016888endodeoxyribonuclease activity, producing 5′-phosphomonoesters1/616/176970.0050.0290.005
KEGGhsa04913Ovarian steroidogenesis1/351/80760.0190.0860.045
KEGGhsa04730Long-term depression1/360/80760.0220.0860.045
KEGGhsa04213Longevity regulating pathway - multiple species1/362/80760.0230.0860.045
KEGGhsa04520Adherens junction1/371/80760.0260.0860.045
KEGGhsa05218Melanoma1/372/80760.0270.0860.045

Details of GO terms and KEGG pathway enrichment in co-upregulated Hub genes.

Verification of Hub Genes

To explore whether the hub genes of COVID-19 and thrombosis are common targets caused by viral induction, we included a third-party data set for corresponding verification. Previous studies have confirmed that viral diseases induce platelet and neutrophil recruitment to the microvascular system to protect host cells from viral infection (Jenne et al., 2013). However, recent research has also shown that compared with that in patients with COVID-19 pneumonia, immunothrombotic vessel occlusion in the lungs is less prominent in patients with influenza pneumonia (Ackermann et al., 2020; Nicolai et al., 2021). We searched the GEO database with “SARS” as the key word, searched using keywords, restricted the screening conditions to “Homo sapiens” and “Expression profiling by array” and limited the sample to peripheral blood. After reading the original studies on the data sets that met the screening conditions individually, we focused on the experimental design and methodology. The data set that met the principle of randomized control and did not violate ethical guidelines was taken as the research object. Finally, GSE5972 was selected for further research. The GSE5972 data set included 50 SARS patients, regardless of age, sex or past medical history. Positive PCR or serum results confirmed that each patient was infected with SARS-CoV. Ten healthy volunteers, five men and five women, were recruited. The chip platform was GPL4387 UHNMAC Homo sapiens 19K Hu19Kv8. Using the aforementioned methodology, we performed background correction and data normalization on the GSE5972 data set and then used the “limma” package of R software to analyze the disease group and the healthy control group to identify DEGs. The inclusion criteria of DEGs were the same as those for the above two data sets. The DEGs of the GSE5972 and GSE19151 data sets were combined and analyzed, the co-DEGs were screened out, and a PPI network was constructed. Then, the hub genes were extracted by Cytoscape software. The final results showed that the upregulated hub genes were IFI44L, TNFSF10, RSAD2, ETFA and MDH1 and that the downregulated hub genes were PXN, MYH9, SPTAN1, TLN1, CXCR4, SUPT6H and ENO1. It can be seen from the results that the hub genes of COVID-19 and thrombosis are not the same as the hub genes of SARS-CoV infection and thrombosis, which also explains why there are many studies on the mechanism of virus-induced thrombosis (Boilard et al., 2014; Tripathi et al., 2014; Koupenova et al., 2019) but the specific mechanism of COVID-19 thrombosis is still unclear (Farkouh et al., 2022; Thomas and Scully, 2022; Vincent et al., 2022). In addition, there have been an increasing number of recent reports of thrombosis after COVID-19 vaccination (De Michele et al., 2022), all of which indicate the existence of characteristic factors of thrombosis in COVID-19, which deserves further exploration.

Discussion

Thrombosis and the inflammatory response were once thought to be two independent physiological processes, but in recent years, studies have gradually found that there is a clear correlation between them. Activation of the coagulation system is part of the host response to pathogens and primarily serves to limit their spread in the bloodstream, which is achieved through interactions between innate immune cells and platelets, termed immunothrombosis (Engelmann and Massberg, 2013). Due to the COVID-19 pandemic, this normal physiological response has received more attention (Stark and Massberg, 2021). Histological analysis of pulmonary vessels in COVID-19 patients showed extensive thrombosis with microangiopathy, and the incidence of pulmonary capillary microthrombosis in COVID-19 patients was 9 times higher than that in influenza patients; the difference was statistically significant (Ackermann et al., 2020). Furthermore, data show that immune thrombosis is also a key indicator of the severity of COVID-19 (Nicolai et al., 2020). Preliminary studies in Wuhan, China, showed that most patients with COVID-19 who died of disease had coagulation disorders, characterized by prolonged prothrombin time and partial prothrombin kinase time, similar to the characteristics observed in patients with disseminated intravascular coagulation (Zhou et al., 2020). When pathogens invade the bloodstream and pose a lethal threat to the host, the body tries to control the infection by all available means, and viral diseases induce the recruitment of platelets and neutrophils to the microvasculature to protect host cells from viral infection (Jenne et al., 2013). The main cellular drivers of this process are platelets and innate immune cells, including neutrophils, monocytes and macrophages. In general, the response to viral disease should be similar; however, immune thrombotic vaso-occlusion was less prominent in patients with influenza-related pneumonia than in patients with COVID-19 (Nicolai et al., 2021). This finding illustrates that the immunothrombosis in COVID-19 is not just a general feature of viral infections but may also be a characteristic factor, suggesting that it is worthy of further exploration to reduce the disease severity and mortality of patients with COVID-19. Bioinformatics analysis enables us to understand the molecular mechanisms of disease initiation and progression, providing a new and effective method for identifying potential biomarkers and therapeutic targets for assessing and treating COVID-19 thrombosis. In this study, six co-upregulated genes, RPS7, IGF1R, DICER1, ERH, MCTS1, and TNPO1, and two co-downregulated genes, FLNA and PXN, were identified from COVID-19 and thrombosis datasets. Ribosomal protein S7 (rps7) is an essential component of the 40s subunit of the ribosome and is essential in the translation process and assembly of the ribosome (Fatica and Tollervey, 2002). Studies have shown that RPS7 regulates cell proliferation, the cell cycle and apoptosis through the PI3K/AKT and MAPK signaling pathways (Wang et al., 2013). However, whether RPS7 is specifically overexpressed in patients with COVID-19 with a high risk of thrombosis remains to be determined. Insulin-like growth factor 1 receptor (IGF1R) is a widely expressed membrane-bound tyrosine kinase receptor that can recognize two main ligands, IGF1 and IGF2. IGF1R expression on lymphocytes and macrophages can cause the proliferation of target cells and the production and activation of cytokines. Phosphorylation of IGF1R exacerbates inflammation, and its overexpression increases cytokine levels during influenza infection (Li et al., 2019). IGF1R has also recently been defined by some studies as a novel biomarker for predicting mortality in severe COVID-19 patients (Fraser et al., 2020). Dicer 1, ribonuclease III (DICER1) is a highly conserved RNase III ribonuclease that plays a key role in biological processes related to microRNAs (miRNAs). Recent studies have shown that DICER1 may maintain homeostasis during inflammation (Jiang et al., 2019). If the inflammatory response is out of balance, the crosstalk between inflammation and thrombosis can have serious consequences because not only does inflammation lead to thrombosis, but thrombosis in turn directly triggers inflammation. Thrombin cleaves fibrinogen and activates the cytokine IL-1α, providing a direct link between coagulation and inflammation. ERH mRNA splicing and mitosis factor (ERH) is a small protein that is highly conserved in eukaryotes. Though it was discovered nearly 30 years ago, its biological function is still largely unknown. It is considered to play a role in transcriptional regulation and the cell cycle. Although a series of studies have found that ERH has some characteristics related to RNA metabolism, the molecular functions of the ERH protein family still need to be further explored (Graille and Rougemaille, 2020). Malignant T cell-amplified sequence 1 (MCTS1) is an that was oncogene originally found in human T-cell lymphoma. MCTS1 can affect the cell cycle by shortening the G1/S conversion time, which is achieved by an imbalance between G1 cyclin activators and inhibitors. It has been confirmed that the expression of MCTS1 increases after DNA damage (Huang et al., 2021). In addition, MCTS1 also acts as a translation enhancer. Translation control has been proven to be an important step in regulating gene expression and subsequently cell growth, differentiation and necrosis (Bohlen et al., 2020). These functions enable MCTS1 to play multiple roles in the occurrence and development of malignant tumors, but its mechanisms in the process of inflammation and thrombosis have not been reported and are worth further exploration. Transportin 1 (TNPO1) is a protein involved in protein transduction into the nucleus, and its main target is RNA binding proteins. Studies have shown that TNPO1 is related to the influenza A virus core shell (Larson et al., 2019; Miyake et al., 2019). Moreover, its expression level can significantly predict ARDS mortality (Liao et al., 2021). This finding indicates that TNPO1 plays a key role in the inflammatory response and the occurrence of viral diseases, which show similar pathway activation to COVID-19 and thrombosis. During the COVID-19 pandemic, the most worrisome complications were ARDS and cardiovascular events, such as venous thromboembolism, myocardial infarction and stroke (Bangalore et al., 2020; Zhang Y. et al., 2020; Poissy et al., 2020; Stefanini et al., 2020). This theoretical evidence provides the basis for studying the mechanism of TNPO1 in COVID-19. In addition, two co-downregulated genes, FLNA and PXN, were obtained through bioinformatics analysis. Filamin A (FLNA) is a large actin-binding cytoskeletal protein. It is important for cell movement because it stabilizes the actin network and integrates it with the cell membrane. There is increasing evidence that FLNA is involved in the pathogenesis of cardiovascular and respiratory diseases. The interaction of FLNA with transcription factors and cell signaling molecules determines the function of vascular cells (Bandaru et al., 2021). Paxillin (PXN) is a focal adhesion protein consisting of an LD motif and a conserved cysteine-rich domain. PXN is thought to play a role in active muscle protein uptake and cell adhesion sites. Although PXN itself lacks enzymatic activity, it contains a variety of domains that can bind signals and structural proteins (Tsai et al., 2017). Therefore, PXN is considered to be a cytoskeletal protein that can effectively transmit signals (Stiegler et al., 2012). Inducers of these two co-downregulated genes have theoretical therapeutic value for COVID-19 patients with thrombosis. Based on the correlation between thrombosis and the severity of COVID-19 (Nicolai et al., 2020), further research may even identify a specific mechanism to reduce the rate of severe COVID-19, thereby reducing the threat of COVID-19 to the elderly and people with underlying diseases.

A limitation of our study is that these identified biomarkers were only verified at a theoretical level. Further experimental studies and clinical trials should be performed for accurate validation.

Conclusion

Based on bioinformatics analysis and previous studies, this study identified hub genes (RPS7, IGF1R, DICER1, ERH, MCTS1, TNPO1, FLNA, and PXN) that may contribute to the evaluation and treatment of COVID-19 and thrombosis. These findings provide new ideas to treat COVID-19 thrombosis and reduce the rate of severe COVID-19, and further clinical studies are needed to explore these possibilities.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE19151, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE164805.

Author contributions

PQ and MH performed bioinformatics and were major contributors in writing the manuscript. TL guided the writing of the manuscript, and critically reviewed the manuscript. All the authors read and approved the final manuscript.

Funding

This research was supported by the Major Project of Military Logistics Scientific Research (AWS17J004), National Key Research and Development Plan for Science and Technology for the Winter Olympics of the Ministry of Science and Technology of China (2019YFF0302300), National Natural Science Foundation of China (82000631), Beijing Natural Science Foundation (7222169), Young Elite Scientist Sponsorship Program by CAST 2020QNRC001 (to MJ.H.) and Military Medical Youth Special Project of PLA General Hospital (QNF19035).

Acknowledgments

The authors thank all staff in the Department of Emergency and Nephrology for their support during the study.

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.

Publisher’s note

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

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.889348/full#supplementary-material

Supplementary Figure S1

Box plots of GSE19151. The dataset included 70 adults who had used warfarin one or more times for venous thromboembolism and 63 healthy controls.

Supplementary Figure S2

Box plots of GSE164805. The dataset included 10 patients diagnosed with COVID-19 and 5 healthy controls.

References

  • 1

    AchimA.KákonyiK.JambrikZ.RuzsaZ. (2022). A Case Report of COVID-19-Associated Acute Hand Ischaemia in a Young Professional Volleyball Player. Eur. Heart J. Case Rep.6, ytac099. 10.1093/ehjcr/ytac099

  • 2

    AckermannM.VerledenS. E.KuehnelM.HaverichA.WelteT.LaengerF.et al (2020). Pulmonary Vascular Endothelialitis, Thrombosis, and Angiogenesis in Covid-19. N. Engl. J. Med.383, 120128. 10.1056/nejmoa2015432

  • 3

    AltermannE.KlaenhammerT. R. (2005). PathwayVoyager: Pathway Mapping Using the Kyoto Encyclopedia of Genes and Genomes (KEGG) Database. BMC Genomics6, 60. 10.1186/1471-2164-6-60

  • 4

    AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene Ontology: Tool for the Unification of Biology. The Gene Ontology Consortium. Nat. Genet.25, 2529. 10.1038/75556

  • 5

    BandaruS.AlaC.ZhouA.-X.AkyürekL. M. (2021). Filamin A Regulates Cardiovascular Remodeling. Ijms22 (12), 6555. 10.3390/ijms22126555

  • 6

    BangaloreS.SharmaA.SlotwinerA.YatskarL.HarariR.ShahB.et al (2020). ST-Segment Elevation in Patients with Covid-19 - A Case Series. N. Engl. J. Med.382, 24782480. 10.1056/nejmc2009020

  • 7

    BleizifferS.ZabelR.GummertJ. F.DeutschM. A. (2022). Large Thrombus in Transit in a COVID-19 Patient. Interact. Cardiovasc Thorac. Surg.ivac077. [Epub ahead of print].10.1093/icvts/ivac077

  • 8

    BohlenJ.HarbrechtL.BlancoS.Clemm Von HohenbergK.FenzlK.KramerG.et al (2020). DENR Promotes Translation Reinitiation via Ribosome Recycling to Drive Expression of Oncogenes Including ATF4. Nat. Commun.11, 4676. 10.1038/s41467-020-18452-2

  • 9

    BoilardE.ParéG.RousseauM.CloutierN.DubucI.LévesqueT.et al (2014). Influenza Virus H1N1 Activates Platelets through FcγRIIA Signaling and Thrombin Generation. Blood123, 28542863. 10.1182/blood-2013-07-515536

  • 10

    BonaventuraA.VecchiéA.DagnaL.MartinodK.DixonD. L.Van TassellB. W.et al (2021). Endothelial Dysfunction and Immunothrombosis as Key Pathogenic Mechanisms in COVID-19. Nat. Rev. Immunol.21, 319329. 10.1038/s41577-021-00536-9

  • 11

    BulteJ. P.PostmaN.BeukemaM.InbergB.StegemanA. G.Van Der HoevenH. (2022). COVID 19 and the Risk of Gastro-Intestinal Perforation: A Case Series and Literature Review. J. Crit. Care67, 100103. 10.1016/j.jcrc.2021.10.020

  • 12

    ChinC. H.ChenS. H.WuH. H.HoC. W.KoM. T.LinC. Y. (2014). cytoHubba: Identifying Hub Objects and Sub-networks from Complex Interactome. BMC Syst. Biol.8 (Suppl. 4), S11. 10.1186/1752-0509-8-S4-S11

  • 13

    ChungM. K.ZidarD. A.BristowM. R.CameronS. J.ChanT.HardingC. V.et al (2021). COVID-19 and Cardiovascular Disease: From Bench to Bedside. Circ. Res.128, 12141236. 10.1161/circresaha.121.317997

  • 14

    CuiS.ChenS.LiX.LiuS.WangF. (2020). Prevalence of Venous Thromboembolism in Patients with Severe Novel Coronavirus Pneumonia. J. Thromb. Haemost.18, 14211424. 10.1111/jth.14830

  • 15

    De micheleM.KahanJ.BertoI.SchiavoO. G.IacobucciM.ToniD.et al (2022). Cerebrovascular Complications of COVID-19 and COVID-19 Vaccination. Circ. Res.130, 11871203. 10.1161/circresaha.122.319954

  • 16

    DevauxC. A.RolainJ.-M.RaoultD. (2020). ACE2 Receptor Polymorphism: Susceptibility to SARS-CoV-2, Hypertension, Multi-Organ Failure, and COVID-19 Disease Outcome. J. Microbiol. Immunol. Infect.53, 425435. 10.1016/j.jmii.2020.04.015

  • 17

    DrakeT. M.RiadA. M.FairfieldC. J.EganC.KnightS. R.PiusR.et al (2021). Characterisation of In-Hospital Complications Associated with COVID-19 Using the ISARIC WHO Clinical Characterisation Protocol UK: a Prospective, Multicentre Cohort Study. Lancet398, 223237. 10.1016/S0140-6736(21)00799-6

  • 18

    EngelmannB.MassbergS. (2013). Thrombosis as an Intravascular Effector of Innate Immunity. Nat. Rev. Immunol.13, 3445. 10.1038/nri3345

  • 19

    FarkouhM. E.StoneG. W.LalaA.BagiellaE.MorenoP. R.NadkarniG. N.et al (2022). Anticoagulation in Patients with COVID-19: JACC Review Topic of the Week. J. Am. Coll. Cardiol.79, 917928. 10.1016/j.jacc.2021.12.023

  • 20

    FaticaA.TollerveyD. (2002). Making Ribosomes. Curr. Opin. Cell Biol.14, 313318. 10.1016/s0955-0674(02)00336-8

  • 21

    FraserD. D.CepinskasG.PattersonE. K.SlessarevM.MartinC.DaleyM.et al (2020). Novel Outcome Biomarkers Identified with Targeted Proteomic Analyses of Plasma from Critically Ill Coronavirus Disease 2019 Patients. Crit. Care Explor2, e0189. 10.1097/cce.0000000000000189

  • 22

    FröhlichG. M.JeschkeE.EichlerU.ThieleH.AlhaririL.ReinthalerM.et al (2021). Impact of Oral Anticoagulation on Clinical Outcomes of COVID-19: a Nationwide Cohort Study of Hospitalized Patients in Germany. Clin. Res. Cardiol.110, 10411050. 10.1007/s00392-020-01783-x

  • 23

    GautierL.CopeL.BolstadB. M.IrizarryR. A. (2004). affy--analysis of Affymetrix GeneChip Data at the Probe Level. Bioinformatics20, 307315. 10.1093/bioinformatics/btg405

  • 24

    Global Burden of Disease 2020 Health Financing Collaborator Network (2021). Tracking Development Assistance for Health and for COVID-19: a Review of Development Assistance, Government, Out-Of-Pocket, and Other Private Spending on Health for 204 Countries and Territories, 1990-2050. Lancet398, 13171343. 10.1016/S0140-6736(21)01258-7

  • 25

    GhasemiyehP.Mohammadi-SamaniS.FirouzabadiN.DehshahriA.VazinA. (2021). A Focused Review on Technologies, Mechanisms, Safety, and Efficacy of Available COVID-19 Vaccines. Int. Immunopharmacol.100, 108162. 10.1016/j.intimp.2021.108162

  • 26

    GrailleM.RougemailleM. (2020). ERH Proteins: Connecting RNA Processing to Tumorigenesis?Curr. Genet.66, 689692. 10.1007/s00294-020-01065-z

  • 27

    GrossO.MoererO.WeberM.HuberT. B.ScheithauerS. (2020). COVID-19-associated Nephritis: Early Warning for Disease Severity and Complications?Lancet395, e87e88. 10.1016/s0140-6736(20)31041-2

  • 28

    HalushkaM. K.Vander heideR. S. (2021). Myocarditis Is Rare in COVID-19 Autopsies: Cardiovascular Findings across 277 Postmortem Examinations. Cardiovasc. Pathol.50, 107300. 10.1016/j.carpath.2020.107300

  • 29

    HanE.TanM. M. J.TurkE.SridharD.LeungG. M.ShibuyaK.et al (2020). Lessons Learnt from Easing COVID-19 Restrictions: an Analysis of Countries and Regions in Asia Pacific and Europe. Lancet396, 15251534. 10.1016/s0140-6736(20)32007-9

  • 30

    HuangZ.SuQ.LiW.RenH.HuangH.WangA. (2021). MCTS1 Promotes Invasion and Metastasis of Oral Cancer by Modifying the EMT Process. Ann. Transl. Med.9, 997. 10.21037/atm-21-2361

  • 31

    JeckelmannC.DjokicB.DuchatelleV.GirodG. (2022). Case Report: Mechanical Mitral Prosthetic Valve Thrombosis in the Context of COVID-19 Despite Effective Anticoagulation. Eur. Heart J. Case Rep.6, ytac006. 10.1093/ehjcr/ytac006

  • 32

    JenneC. N.WongC. H. Y.ZempF. J.McdonaldB.RahmanM. M.ForsythP. A.et al (2013). Neutrophils Recruited to Sites of Infection Protect from Virus Challenge by Releasing Neutrophil Extracellular Traps. Cell Host Microbe13, 169180. 10.1016/j.chom.2013.01.005

  • 33

    JiangC.XuJ.ZhuW.CaiY.WangS.GuoY.et al (2019). Abnormal Expression of DICER1 Leads to Dysregulation of Inflammatory Effectors in Human Synoviocytes. Mediat. Inflamm.2019, 6768504. 10.1155/2019/6768504

  • 34

    KlokF. A.KruipM. J. H. A.Van Der MeerN. J. M.ArbousM. S.GommersD. A. M. P. J.KantK. M.et al (2020). Incidence of Thrombotic Complications in Critically Ill ICU Patients with COVID-19. Thrombosis Res.191, 145147. 10.1016/j.thromres.2020.04.013

  • 35

    KoupenovaM.CorkreyH. A.VitsevaO.ManniG.PangC. J.ClancyL.et al (2019). The Role of Platelets in Mediating a Response to Human Influenza Infection. Nat. Commun.10, 1780. 10.1038/s41467-019-09607-x

  • 36

    LalA.EronduN. A.HeymannD. L.GitahiG.YatesR. (2021). Fragmented Health Systems in COVID-19: Rectifying the Misalignment between Global Health Security and Universal Health Coverage. Lancet397, 6167. 10.1016/s0140-6736(20)32228-5

  • 37

    LalaA.JohnsonK. W.JanuzziJ. L.RussakA. J.ParanjpeI.RichterF.et al (2020). Prevalence and Impact of Myocardial Injury in Patients Hospitalized with COVID-19 Infection. J. Am. Coll. Cardiol.76, 533546. 10.1016/j.jacc.2020.06.007

  • 38

    LarsonG. P.TranV.S.CaìY.HigginsC. A.SmithD. M.et al (2019). EPS8 Facilitates Uncoating of Influenza A Virus. Cell Rep.29, 21752183. 10.1016/j.celrep.2019.10.064

  • 39

    LawlerP. R.GoligherE. C.BergerJ. S.LawlerP. R.GoligherE. C.BergerJ. S.et al (2021). Therapeutic Anticoagulation with Heparin in Noncritically Ill Patients with Covid-19. N. Engl. J. Med.385, 790802. 10.1056/NEJMoa2105911

  • 40

    LeviM.ThachilJ.IbaT.LevyJ. H. (2020). Coagulation Abnormalities and Thrombosis in Patients with COVID-19. Lancet Haematol.7, e438e440. 10.1016/s2352-3026(20)30145-9

  • 41

    LiG.ZhouL.ZhangC.ShiY.DongD.BaiM.et al (2019). Insulin-Like Growth Factor 1 Regulates Acute Inflammatory Lung Injury Mediated by Influenza Virus Infection. Front. Microbiol.10, 2541. 10.3389/fmicb.2019.02541

  • 42

    LiM.ZhangX.AngK. S.LingJ.SethiR.LeeN. Y. S.et al (2022). DISCO: a Database of Deeply Integrated Human Single-Cell Omics Data. Nucleic Acids Res.50, D596d602. 10.1093/nar/gkab1020

  • 43

    LiaoS. Y.CasanovaN. G.BimeC.CampS. M.LynnH.GarciaJ. G. N. (2021). Identification of Early and Intermediate Biomarkers for ARDS Mortality by Multi-Omic Approaches. Sci. Rep.11, 18874. 10.1038/s41598-021-98053-1

  • 44

    LodigianiC.IapichinoG.CarenzoL.CecconiM.FerrazziP.SebastianT.et al (2020). Venous and Arterial Thromboembolic Complications in COVID-19 Patients Admitted to an Academic Hospital in Milan, Italy. Thrombosis Res.191, 914. 10.1016/j.thromres.2020.04.024

  • 45

    LowensteinC. J.SolomonS. D. (2020). Severe COVID-19 Is a Microvascular Disease. Circulation142, 16091611. 10.1161/circulationaha.120.050354

  • 46

    MaierC. L.TruongA. D.AuldS. C.PollyD. M.TanksleyC.-L.DuncanA. (2020). COVID-19-associated Hyperviscosity: a Link between Inflammation and Thrombophilia?Lancet395, 17581759. 10.1016/s0140-6736(20)31209-5

  • 47

    McintyreP. B.AggarwalR.JaniI.JawadJ.KochharS.MacdonaldN.et al (2022). COVID-19 Vaccine Strategies Must Focus on Severe Disease and Global Equity. Lancet399, 406410. 10.1016/s0140-6736(21)02835-x

  • 48

    MiH.MuruganujanA.HuangX.EbertD.MillsC.GuoX.et al (2019). Protocol Update for Large-Scale Genome and Gene Function Analysis with the PANTHER Classification System (v.14.0). Nat. Protoc.14, 703721. 10.1038/s41596-019-0128-8

  • 49

    MiyakeY.KeuschJ. J.DecampsL.Ho-XuanH.IketaniS.GutH.et al (2019). Influenza Virus Uses Transportin 1 for vRNP Debundling during Cell Entry. Nat. Microbiol.4, 578586. 10.1038/s41564-018-0332-2

  • 50

    NadkarniG. N.LalaA.BagiellaE.ChangH. L.MorenoP. R.PujadasE.et al (2020). Anticoagulation, Bleeding, Mortality, and Pathology in Hospitalized Patients with COVID-19. J. Am. Coll. Cardiol.76, 18151826. 10.1016/j.jacc.2020.08.041

  • 51

    NicolaiL.LeunigA.BrambsS.KaiserR.JoppichM.HoffknechtM. L.et al (2021). Vascular Neutrophilic Inflammation and Immunothrombosis Distinguish Severe COVID‐19 from Influenza Pneumonia. J. Thromb. Haemost.19, 574581. 10.1111/jth.15179

  • 52

    NicolaiL.LeunigA.BrambsS.KaiserR.WeinbergerT.WeigandM.et al (2020). Immunothrombotic Dysregulation in COVID-19 Pneumonia Is Associated with Respiratory Failure and Coagulopathy. Circulation142, 11761189. 10.1161/circulationaha.120.048488

  • 53

    PellegriniD.KawakamiR.GuagliumiG.SakamotoA.KawaiK.GianattiA.et al (2021). Microthrombi as a Major Cause of Cardiac Injury in COVID-19: A Pathologic Study. Circulation143, 10311042. 10.1161/circulationaha.120.051828

  • 54

    PoissyJ.GoutayJ.CaplanM.ParmentierE.DuburcqT.LassalleF.et al (2020). Pulmonary Embolism in Patients with COVID-19: Awareness of an Increased Prevalence. Circulation142, 184186. 10.1161/circulationaha.120.047430

  • 55

    RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic acids Res.43, e47. 10.1093/nar/gkv007

  • 56

    RoncoC.ReisT.Husain-syedF. (2020). Management of Acute Kidney Injury in Patients with COVID-19. Lancet Respir. Med.8, 738742. 10.1016/s2213-2600(20)30229-0

  • 57

    RosenthalR.CadieuxE. L.CadieuxE. L.SalgadoR.BakirM. A.MooreD. A.et al (2019). Neoantigen-directed Immune Escape in Lung Cancer Evolution. Nature567, 479485. 10.1038/s41586-019-1032-7

  • 58

    RosovskyR. P.SanfilippoK. M.WangT. F.RajanS. K.ShahS.MartinK. A.et al (2020). Anticoagulation Practice Patterns in COVID‐19: A Global Survey. Res. Pract. Thromb. Haemost.4, 969983. 10.1002/rth2.12414

  • 59

    RussellJ.WagonerM.DupontJ.MyersD.MuthuK.ThotakuraS. (2022). Left Ventricular Thrombus of Unknown Etiology in a Patient with COVID-19 Disease with No Significant Medical History. Cardiovasc. Revasc. Med. S15538389(22)00121-X. 10.1016/j.carrev.2022.03.010

  • 60

    SastryS.CuomoF.MuthusamyJ. (2022). COVID-19 and Thrombosis: The Role of Hemodynamics. Thrombosis Res.212, 5157. 10.1016/j.thromres.2022.02.016

  • 61

    SerbanD.TribusL. C.VanceaG.StoianA. P.DascaluA. M.SuceveanuA. I.et al (2021). Acute Mesenteric Ischemia in COVID-19 Patients. Jcm11 (1), 200. 10.3390/jcm11010200

  • 62

    ShankarR.LeimanisM. L.NewburyP. A.LiuK.XingJ.NedveckD.et al (2020). Gene Expression Signatures Identify Paediatric Patients with Multiple Organ Dysfunction Who Require Advanced Life Support in the Intensive Care Unit. EBioMedicine62, 103122. 10.1016/j.ebiom.2020.103122

  • 63

    ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res.13, 24982504. 10.1101/gr.1239303

  • 64

    ShaoF.XuS.MaX.XuZ.LyuJ.NgM.et al (2020). In-hospital Cardiac Arrest Outcomes Among Patients with COVID-19 Pneumonia in Wuhan, China. Resuscitation151, 1823. 10.1016/j.resuscitation.2020.04.005

  • 65

    SriwastavaS.TandonM.PoduryS.PrasadA.WenS.GuthrieG.et al (2021). COVID-19 and Neuroinflammation: a Literature Review of Relevant Neuroimaging and CSF Markers in Central Nervous System Inflammatory Disorders from SARS-COV2. J. Neurol.268, 44484478. 10.1007/s00415-021-10611-9

  • 66

    StarkK.MassbergS. (2021). Interplay between Inflammation and Thrombosis in Cardiovascular Pathology. Nat. Rev. Cardiol.18, 666682. 10.1038/s41569-021-00552-1

  • 67

    StefaniniG. G.MontorfanoM.TrabattoniD.AndreiniD.FerranteG.AnconaM.et al (2020). ST-elevation Myocardial Infarction in Patients with COVID-19: Clinical and Angiographic Outcomes. Circulation141, 21132116. 10.1161/circulationaha.120.047525

  • 68

    StieglerA. L.DraheimK. M.LiX.ChayenN. E.CalderwoodD. A.BoggonT. J. (2012). Structural Basis for Paxillin Binding and Focal Adhesion Targeting of β-Parvin. J. Biol. Chem.287, 3256632577. 10.1074/jbc.m112.367342

  • 69

    SuhY. J.HongH.OhanaM.BompardF.RevelM.-P.ValleC.et al (2021). Pulmonary Embolism and Deep Vein Thrombosis in COVID-19: A Systematic Review and Meta-Analysis. Radiology298, E70e80. 10.1148/radiol.2020203557

  • 70

    TakasuS.AriizumiM.MatsumotoS.NakagawaH.IwadateK. (2022). Cerebral Venous Sinus Thrombosis Associated with COVID-19: an Autopsy Case Report. Forensic Sci. Med. Pathol.18, 8085. 10.1007/s12024-022-00458-5

  • 71

    ThomasM.ScullyM. (2022). Clinical Features of Thrombosis and Bleeding in COVID-19. Blood. 10.1182/blood.2021012247

  • 72

    TripathiS.VermaA.KimE.-J.WhiteM. R.HartshornK. L. (2014). LL-37 Modulates Human Neutrophil Responses to Influenza A Virus. J. Leukoc. Biol.96, 931938. 10.1189/jlb.4a1113-604rr

  • 73

    TsaiW.-C.YuT.-Y.LinL.-P.LinM.-S.TsaiT.-T.PangJ.-H. S. (2017). Platelet Rich Plasma Promotes Skeletal Muscle Cell Migration in Association with Up-Regulation of FAK, Paxillin, and F-Actin Formation. J. Orthop. Res.35, 25062512. 10.1002/jor.23547

  • 74

    UhlenM.KarlssonM. J.ZhongW.TebaniA.PouC.MikesJ.et al (2019). A Genome-wide Transcriptomic Analysis of Protein-Coding Genes in Human Blood Cells. Science366 (6472), eaax9198. 10.1126/science.aax9198

  • 75

    VincentJ.-L.LeviM.HuntB. J. (2022). Prevention and Management of Thrombosis in Hospitalised Patients with COVID-19 Pneumonia. Lancet Respir. Med.10, 214220. 10.1016/s2213-2600(21)00455-0

  • 76

    WalterW.Sánchez-CaboF.RicoteM. (2015). GOplot: an R Package for Visually Combining Expression Data with Functional Analysis: Fig. 1. Bioinformatics31, 29122914. 10.1093/bioinformatics/btv300

  • 77

    WangJ.ZhuK.XueY.WenG.TaoL. (2021). Research Progress in the Treatment of Complications and Sequelae of COVID-19. Front. Med.8, 757605. 10.3389/fmed.2021.757605

  • 78

    WangZ.HouJ.LuL.QiZ.SunJ.GaoW.et al (2013). Small Ribosomal Protein Subunit S7 Suppresses Ovarian Tumorigenesis through Regulation of the PI3K/AKT and MAPK Pathways. PLoS One8, e79117. 10.1371/journal.pone.0079117

  • 79

    WeiH.LiangG.VeraR. M. (2021). Dantrolene Repurposed to Treat Sepsis or Septic Shock and COVID-19 Patients. Eur. Rev. Med. Pharmacol. Sci.25, 31363144. 10.26355/eurrev_202104_25569

  • 80

    YongzhiX. (2021). COVID-19-associated Cytokine Storm Syndrome and Diagnostic Principles: an Old and New Issue. Emerg. Microbes Infect.10, 266276. 10.1080/22221751.2021.1884503

  • 81

    YuG.WangL.-G.HanY.HeQ.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS A J. Integr. Biol.16, 284287. 10.1089/omi.2011.0118

  • 82

    ZangrilloA.BellettiA.PalumboD.CalviM. R.GuzzoF.FominskiyE. V.et al (2022). One-Year Multidisciplinary Follow-Up of Patients with COVID-19 Requiring Invasive Mechanical Ventilation. J. Cardiothorac. Vasc. Anesth.36, 13541363. 10.1053/j.jvca.2021.11.032

  • 83

    ZhangL.YanX.FanQ.LiuH.LiuX.LiuZ.et al (2020). D‐dimer Levels on Admission to Predict In‐hospital Mortality in Patients with Covid‐19. J. Thromb. Haemost.18, 13241329. 10.1111/jth.14859

  • 84

    ZhangY.XiaoM.ZhangS.XiaP.CaoW.JiangW.et al (2020). Coagulopathy and Antiphospholipid Antibodies in Patients with Covid-19. N. Engl. J. Med.382, e38. 10.1056/nejmc2007575

  • 85

    ZhouF.YuT.DuR.FanG.LiuY.LiuZ.et al (2020). Clinical Course and Risk Factors for Mortality of Adult Inpatients with COVID-19 in Wuhan, China: a Retrospective Cohort Study. Lancet395, 10541062. 10.1016/s0140-6736(20)30566-3

  • 86

    ZuoY.EstesS. K.AliR. A.GandhiA. A.YalavarthiS.ShiH.et al (2020). Prothrombotic Autoantibodies in Serum from Patients Hospitalized with COVID-19. Sci. Transl. Med.12 (570), eabd3876. 10.1126/scitranslmed.abd3876

  • 87

    ZuoY.ZuoM.YalavarthiS.GockmanK.MadisonJ. A.ShiH.et al (2021). Neutrophil Extracellular Traps and Thrombosis in COVID-19. J. Thromb. Thrombolysis51, 446453. 10.1007/s11239-020-02324-z

Summary

Keywords

COVID-19, thrombosis, therapy, bioinformatic analysis, differentially expressed genes

Citation

Qi P, Huang M and Li T (2022) Screening the Potential Biomarkers of COVID-19-Related Thrombosis Through Bioinformatics Analysis. Front. Genet. 13:889348. doi: 10.3389/fgene.2022.889348

Received

07 March 2022

Accepted

09 May 2022

Published

25 May 2022

Volume

13 - 2022

Edited by

Manal S. Fawzy, Faculty of Medicine, Suez Canal University, Egypt

Reviewed by

Aaron N. Johnson, Washington University in St. Louis, United States

Chien-Ling Huang, Hong Kong Polytechnic University, Hong Kong SAR, China

Updates

Copyright

*Correspondence: Tanshi Li,

†These authors have contributed equally to this work and share first authorship

This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics

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