Identification of common differentially expressed genes (cDEGs) between COVID-19 and VAE
As regards the COVID-19, 814 differentially expressed genes were identified from GSE147507 dataset, including 419 down-regulated genes and 395 up-regulated genes. As for the VAE, 477 differentially expressed genes were identified from GSE30326, including 474 down-regulated genes and 3 up-regulated genes. Eighty-eight cDEGs were identified in both datasets (Fig. 1).
Enrichment analysis of the GO and KEGG
The current study analyzes GO terms including biological process (BP), cellular component (CC) and molecular functions (MF) for 88 cDEGs, as well as KEGG pathway. Gene ontology revealed that in terms of BP, these cDEGs were mainly involved in mitotic nuclear division, defense response to virus and organelle fission. From the perspective of CC, it mainly involves in spindle and condensed chromosome. As for MF, cDEGs were involved in binding and motor activity, histone kinase activity and chemokine activity. KEGG Enrichment analysis on the 88 cDEGs showed that the top 5 signal pathways which were significantly enriched were in Influenza A, NOD-like receptor signaling pathway, Cell cycle, Viral protein interaction with cytokine and cytokine receptor and Oocyte meiosis (Tab. 1). Visualization of GO and KEGG enrichment analysis see Fig. 2
Tab. 1 Top 5 GO_KEGG terms and their corresponding P values for cGEGs
ONTOLOGY
|
ID
|
Description
|
P. adjust
|
BP
|
GO:0140014
|
mitotic nuclear division
|
6.04e-18
|
|
GO:0051607
|
defense response to virus
|
1.56e-17
|
|
GO:0009615
|
response to virus
|
1.79e-17
|
|
GO:0048285
|
organelle fission
|
7.76e-15
|
|
GO:0000280
|
nuclear division
|
1.36e-14
|
CC
|
GO:0005819
|
spindle
|
1.54e-09
|
|
GO:0005876
|
spindle microtubule
|
4.22e-09
|
|
GO:0072686
|
mitotic spindle
|
4.22e-09
|
|
GO:0000779
|
condensed chromosome, centromeric region
|
1.01e-08
|
|
GO:0000775
|
chromosome, centromeric region
|
2.01e-08
|
MF
|
GO:0008017
|
microtubule binding
|
9.70e-05
|
|
GO:0035173
|
histone kinase activity
|
8.00e-04
|
|
GO:0042379
|
chemokine receptor binding
|
8.00e-04
|
|
GO:0003774
|
motor activity
|
8.00e-04
|
|
GO:0008009
|
chemokine activity
|
8.00e-04
|
KEGG
|
hsa05164
|
Influenza A
|
5.66e-05
|
|
hsa04621
|
NOD-like receptor signaling pathway
|
3.31e-04
|
|
hsa04110
|
Cell cycle
|
7.31e-04
|
|
hsa04061
|
Viral protein interaction with cytokine and cytokine receptor
|
0.003
|
|
hsa04114
|
Oocyte meiosis
|
0.003
|
Construction of Protein-protein interactions (PPIs) network
The analyzed network holds 87 nodes and 638 edges (Fig. 3). Average node degree: 14.7. Expected number of edges: 31. Avg. local clustering coefficient: 0.632. It is suggested that these proteins have more interactions among themselves than what would be expected for a random set of proteins of similar size, drawn from the genome.
Detection of hub genes and module analysis
According to topology analysis, 10 cDEGs degree values are the highest, and these genes are considered as Hub genes (Tab. 3). The proteins encoded by the hub genes have rich interactions with other proteins, including 52 nodes and 618 edges (Fig. 4). Module analysis showed that TRIM38, GBP3, STAT1, CASP5, C19orf66 and other genes were high-density modules, including 41 spots and 126 edges. In addition, except for STAT1, which is both hub gene and high-density module, other hub genes are closely related to high-density module (Fig. 5).
Tab. 3 Topological result exploration for top ten hub genes
Hub
gene
|
Degree
|
EPC
|
BottleNeck
|
EcCentricity
|
Closeness
|
Radiality
|
Betweenness
|
Stress
|
Clustering
Coefficient
|
STAT1
|
43
|
24.3
|
6
|
0.3
|
52.16
|
4.20
|
421.65
|
2256
|
0.58
|
MX1
|
42
|
25.3
|
4
|
0.3
|
51.33
|
4.25
|
158.83
|
1406
|
0.63
|
IFIH1
|
42
|
24.3
|
1
|
0.3
|
51.16
|
4.24
|
206.80
|
1576
|
0.62
|
ISG15
|
40
|
24.2
|
3
|
0.3
|
50
|
4.19
|
124.12
|
1206
|
0.66
|
DDX58
|
40
|
23.9
|
5
|
0.3
|
50.33
|
4.22
|
274.77
|
1708
|
0.65
|
RSAD2
|
40
|
24.4
|
2
|
0.3
|
50.17
|
4.20
|
105.99
|
1110
|
0.67
|
IFIT3
|
39
|
25.0
|
2
|
0.3
|
49.33
|
4.17
|
83.57
|
982
|
0.70
|
IRF7
|
38
|
24.3
|
1
|
0.3
|
49.17
|
4.18
|
98.44
|
958
|
0.68
|
CXCL10
|
38
|
22.6
|
6
|
0.3
|
50
|
4.25
|
639.64
|
2948
|
0.60
|
IFIT1
|
37
|
24.2
|
4
|
0.3
|
48.33
|
4.14
|
47.54
|
748
|
0.76
|
Transcriptional factor regulatory network analysis of hub genes
By NetworkAnalyst analysis, 123 TFs regulating eight hub genes expression were identified. The main hub genes regulated by TF have ISG15, RSAD2, IFIT3. NF-κb1, SP1 and RALA are the most widely acting transcription factors that regulate the ISG15, RSAD2, IRF7, IFIT1 and IFIT3 genes (Fig. 6).
TF-miRNA co-regulatory network
By NetworkAnalyst analysis, 308 miRNAs that regulated the expression level of hub genes were identified, and the main co-regulated hub genes were STAT1, CXCL10, IDG15, DDX58, IFIH1, IFIT3 et al. (Fig. 7). The frequency of co-regulation of miR-186, miR-135a, miR-362-5p, miR-409-3p and miR-496 with TF was higher.
Screening of potential therapeutic drugs
Enrichr Platform analysis found that suloctidil HL60 UP, 3'-Azido-3' -Deoxythymidine CTD 00007047, estradiol CTD 00005920 and acetaminophen CTD 00005295 have been identified as potential therapeutic targets for both COVID-19 and VAE, which target genes contain hub genes (Tab. 4). According to hub Genes, TCM formulas with therapeutic effects were searched in Pubmed database, and the results showed that YuPingFengSan and QingFeiPaiDu Decoction had therapeutic effects (Tab. 5).
Tab. 4 Potential drug components for treatment of both COVID-19 and VAE
Drug
|
Adj. P value
|
cDEGs
|
suloctidil HL60 UP
|
6.85E-57
|
SAMD9;IFIT5;IFI6;UBE2L6;ZC3HAV1;IFIT1;IFI44L;
IFIT3;IFIT2;OASL;IFIH1;HERC5;CCL8;TNFSF10;
TRIM21;TRIM22;HERC6;RSAD2;GCH1;DDX58;STAT1;PHF11;MX1;IFI44;EIF2AK2;ISG15;CXCL10;CXCL11;AIM2;OAS2;OAS3;TFEC;IRF7;C19ORF66;CD69;XAF1;SP140L;IRF9
|
3'-Azido-3'-deoxythymidine CTD 00007047
|
2.15E-23
|
SAMD9L;IFIT5;IFI6;DDX60L;IFIT1;IFI44L;IFIT2;IFIH1;CASP5;LAMP3;C3AR1;EPSTI1;TNFSF10;LGALS9;TRIM22;GBP4;HERC6;RSAD2;STAT1;MX1;ISG15;PARP9;OAS2;OAS3;IRF7;C19ORF66;PIK3AP1
|
estradiol CTD 00005920
|
6.29E-08
|
SAMD9L;SAMD9;IFIT5;IFI6;UBE2L6;DDX60L;IFIT1;TARP;IFI44L;IFIT3;IFIT2;OASL;IFIH1;HERC5;CASP5;LAMP3;TNFSF10;STAT4;SLC39A8;TRIM21;TNS1;HERC6;RSAD2;GCH1;DDX58;STAT1;SPHK1;ANKRD22;P2RY14;MX1;CYBB;IFI44;EIF2AK2;ISG15;NMI;PARP9;CXCL10;CXCL11;OAS2;OAS3;IRF7;CMPK2;XAF1;HAMP;IRF9
|
acetaminophen CTD 00005295
|
1.74E-07
|
IFI6;TMEM51;UBE2L6;DDX60L;IFIT1;IFIT3;IFIT2;IFIH1;PSTPIP2;LAMP3;C3AR1;TNFSF10;NBN;TRIM21;TRIM22;GBP4;TNS1;GBP3;RSAD2;DYNLT1;DDX58;STAT1;SPHK1;PHF11;MX1;IFI44;EIF2AK2;NMI;RNASE2;PARP9;PATL1;CXCL10;XRN1;OAS3;IRF7;CMPK2;C19ORF66;CD69;TRIM38;XAF1;SP140L;HAMP;IRF9
|
Tetradioxin CTD 00006848
|
1.70E-06
|
SAMD9;IFIT5;IFI6;IFIT1;TARP;IFI44L;IFIT3;OASL;IFIH1;HERC5;LAMP3;TNFSF10;TRIM21;GBP4;KCNJ2;TNS1;HERC6;GBP3;GBP5;RSAD2;DYNLT1;DDX58;STAT1;MX1;IFI44;ISG15;NMI;PARP9;CXCL10;CXCL11;OAS2;LPAR6;OAS3;IRF7;CD69;HAMP;IRF9;IDO1;C1QC
|
cyclosporin A CTD 00007121
|
6.70E-04
|
SAMD9L;SAMD9;IFIT5;IFI6;UBE2L6;ZC3HAV1;IFIT1;IFI44L;IFIT3;IFIT2;OASL;LILRA5;HERC5;PSTPIP2;LAMP3;TNFSF10;SLC39A8;TRIM22;TNS1;HERC6;GCH1;DDX58;STAT1;MX1;IFI44;ISG15;NMI;PARP9;F5;CXCL10;OAS2;LPAR6;OAS3;CARD16;CD69;TRIM38;PIK3AP1;XAF1;IRF9
|
Tab. 5 Potential TCM formulas for treatment of both COVID-19 and VAE
TCM Formula
|
Composition
|
TCM Theory
|
Hub gene
|
Yu Ping Feng San
|
Fructus Forsythiae, Flos Lonicerae, Radix Platycodonis, Herba Menthae, Herba Lophatheri, Radix Glycyrrhizae, Herba Schizonepetae, Fermented soybean, Fructus arctii, and Rhizoma Phragmitis
|
Disperses wind-heat, clears heat, and relieves toxicity
|
ISG15
|
Qing Fei Pai Du decoction
|
Herba Ephedrae, Radix Glycyrrhizae, Semen Armeniacae Amarum, Gypsum Fibrosum, Ramulus Cinnamomi, Rhizoma Alismatis, Polyporus Umbellatus, Rhizoma
Atractylodis Macrocephalae, Poria, Radix Bupleuri, Radix Scutellariae, Rhizome Pinelliae Preparata, Rhizoma Zingiberis Recens, Radix Asteris, Flos
Farfarae, Rhizoma elamcandae, Herba Asari, Rhizoma Dioscoreae, Fructus Aurantii Immaturus, Pericarpium Citri Reticulatae, Herba ogostemonis
|
Dispelling cold,
clears damp and
heat
|
IFIT1
IFIT3
IFIH1
|