1. Introduction
Monkeypox is a zoonotic disease caused by the monkeypox virus (MPXV), belonging to the Orthopoxvirus genus, which includes double-stranded DNA viruses [Reference Isidro, Borges and Pinto1, Reference Faye, Pratt and Faye2]. The MPXV has been reported in Western and Central African countries, the United States, and Europe [Reference Adler, Gould and Hine3]. MPXV penetrates the body through direct contact via the skin, airways, or mucous membranes such as nose, eyes, or mouth from an infected person, through respiratory droplets of face-to-face contact and through indirect contact via contaminated materials [Reference Murphy4]. The symptoms of MPXV infection are similar to those of a milder form of smallpox [Reference Jayswal and Kakadiya5]. The distinction is that lymphadenopathy is caused by MPXV infection rather than smallpox [Reference Kumar, Acharya, Gendelman and Byrareddy6]. The MPXV has split into two distinct clades: West Africa and the Congo Basin [Reference Likos, Sammons and Olson7]. The diseases caused by the two MPXV clades differ in clinical and epidemiological characteristics. The Congo Basin strain has a 10% mortality rate [Reference Doshi, Guagliardo and Doty8]. West Africa has a fatality rate of about 1%, with patients with the human immunodeficiency virus (HIV) coinfection having a higher mortality rate [Reference Ray and Maunsell9].
In comparison to the West African clade, the Congo Basin monkeypox clade selectively suppresses host responses such as apoptosis and growth factor responses [Reference Kindrachuk, Arsenault and Kusalik10]. By repressing cognate T cell activation, MPXV avoids antiviral CD4+ and CD8+ T cell responses [Reference Hammarlund, Dasgupta, Pinilla, Norori, Fruh and Slifka11]. Realegeno et al. have shown that the high-throughput genetic screen of MPXV infection in haploid cell is involved in Golgi-associated retrograde protein complex [Reference Realegeno, Puschnik and Kumar12]. Monkeypox requires heat shock factor 1 (HSF1) for infection [Reference Filone, Caballero and Dower13]. MPXV infection resulted in significant increases in the number of natural killer cells and cell proliferation [Reference Song, Josleyn and Janosko14].
The outbreak of coronavirus disease (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread globally in December 8, 2020, according to the World Health Organization [Reference Bhimraj, Morgan and Shumaker15]. COVID-19 has affected in various diseases such as inflammatory response and nervous system symptoms such as ischemic stroke and delineation [Reference Nordvig, Fong and Willey16, Reference Ramos-Casals, Brito-Zeron and Mariette17]. Molecular and immunological diagnostic techniques have been developed to identify SARS-CoV-2 infection [Reference Rotondo, Martini and Maritati18]. Currently, the most reliable method for diagnosing SARS-CoV-2 is real-time PCR-based assays for viral RNA detection [Reference Rotondo, Martini and Maritati18, Reference Kevadiya, Machhi and Herskovitz19]. Furthermore, more advanced molecular approaches for detecting SARS-CoV-2 RNA are being developed, including next-generation sequencing (NGS), droplet-digital PCR (ddPCR), and clustered regularly interspaced short palindromic repeats (CRISPRs) [Reference Rotondo, Martini and Maritati18, Reference Hou, Zeng and Yang20–Reference Chen, Kang and Luo22]. Recently, the number of monkeypox cases has increased in 2022 after the COVID-19 outbreak; however, the natural reservoir and function of monkeypox are unknown.
Several studies have reported the possible relationship between monkeypox and COVID-19. One connection is that patients with SARS-CoV-2 infection are more prone to coinfect with MPXV due to the reduced amount of circulating immune cells [Reference Roushdy23]. Another link between the SARS-CoV-2 and MPXV outbreaks is a lack of immunity in younger generations to the smallpox virus. Smallpox is in the same family as monkeypox, and the discontinuation of its vaccine contributes to increased monkeypox cases as well as spread of SARS-CoV-2 by lowering immunity among local residents [Reference Roushdy23, Reference Haider, Guitian and Simons24]. It is still unclear whether the new monkeypox outbreak is a distinct phenomenon or if it is exacerbated by the COVID-19 pandemic [Reference Farahat, Abdelaal and Shah25]. Because of the current small patient sample size, the exact cause-effect relationship between SARS-CoV-2 and MPXV coinfection cannot yet be established [Reference Farahat, Abdelaal and Shah25, Reference El-Qushayri, Reda and Shah26].
In this study, we aimed to investigate the possible molecular function of monkeypox by performing bioinformatic analyses. In addition, by comparing our data of monkeypox infection with published transcriptome of SARS-CoV-2 infection in other cell lines, we found that the function of monkey pox and COVID-19 infection is similar in terms of cytokine signaling in the immune system and TNF signaling. Therefore, the potential functional link of monkeypox and COVID-19 is described in this bioinformatic study.
2. Materials and Methods
2.1. The Collection of Databases and the Identification of DEGs
To create a dataset of monkeypox infection, DEGs of monkeypox infection were obtained from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) database. Two GEO datasets were collected to obtain the DEGs related to monkeypox, including GSE21001 and GSE36854, and the datasets were analyzed using the GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/) web tool. GSE21001 was the database for GeneChip rhesus macaque genome microarrays on a 3- and 7-hour postinfection (hpi) time point after an MPXV infection on Macaca mulatta kidney epithelial cells (MK2) [Reference Alkhalil, Hammamieh, Hardick, Ichou, Jett and Ibrahim27]. GSE36854 was the database for the gene expression profile at 6 hpi after infection with MPXV strain MSF#6 on HeLa cells [Reference Bourquain, Dabrowski and Nitsche28]. The cutoff criteria were obtained for GSE21001 and GSE36854 using an adjusted P value <0.05.
2.2. Gene Ontology (GO) and KEGG Pathway Enrichment Analyses of DEGs of Monkeypox
To examine the enrichment analysis of 212 common differentially expressed genes (DEGs) from the GSE21001 and GSE36854 datasets of MPXV infection, KEGG pathway and GO analyses were performed using the web tool Shiny Go (ShinyGO 0.76 (https://sdstate.edu)) [Reference Ge, Jung and Yao29, Reference Kanehisa, Furumichi, Sato, Ishiguro-Watanabe and Tanabe30]. For the enrichment biological pathways of the comparison between monkeypox and COVID-19, Metascape/COVID analysis (https://metascape.org/COVID) was carried out [Reference Zhou, Zhou and Pache31].
2.3. Identify Hub Genes of MPXV Infection
The cytoHubba plugins of Cytoscape were used to identify the top 20 nodes (hub genes), ranked using the maximal clique centrality (MCC) algorithm of cytoHubba. The Metascape tool [Reference Zhou, Zhou and Pache31] was to identify the Molecular Complex Detection (MCODE) components for functional gene clusters.
2.4. DisGeNET, PaGenBase, and TRRUST for Enrichment Analyses
DisGeNET, PaGenBase, and TRRUST were used to analyze the monkeypox infection, and comparison between monkeypox and COVID-19 was performed through Metascape (https://metascape.org). DisGeNET is a platform for the genetic underpinning of human diseases [Reference Pinero, Queralt-Rosinach and Bravo32, Reference Pinero, Bravo and Queralt-Rosinach33]. PaGenBase is a database for obtaining information on pattern genes under various tissue and time specific conditions [Reference Pan, Hu and Shi34]. TRRUST is a platform for transcriptional regulatory networks [Reference Han, Shim and Shin35].
3. Results
3.1. Design
To determine the function of the MPXV infection, the public datasets GSE36854 and GSE21001 were used. Figure 1(a) shows a flowchart of this study. GEO2R was used to obtain DEGs from the GEO database. Adjusted P values <0.05 were considered as the criteria for the DEGs. Figure 1(b) shows the intersections of the two groups of GSE36854 and GSE21001, indicating the common DEGs.
3.2. Screening of Hub Genes and Important Modules of MPXV Infection
We constructed the protein-protein interaction networks for 211 DEGs of GSE36854 and GSE21001 of the MPXV infection using the Cytoscape and cytoHubba plugins. Twenty hub genes from 211 DEGs of GSE36854 and GSE21001, including HIST1H2AE, HIST2H2A B, HST1H2BM, HIST2H2AC, HIST1H3D, HIST1H2BH, HIST1H2AC, HIST1H2BJ, FOS, CSF2, RELB, IL-10, TNFRSF1B, CXCL8, ICAM1, CCL2, TNF, IL-6, TLR2, and TNFAIP3, were obtained (Figure 2(a)).
To obtain the functional clusters of MPXV infection, a module analysis was carried out using the Metascape tool. In the MCODE analysis by Metascape, 211 DEGs of GSE36854 and GSE21001 of MPXV infection were grouped into five modules. The pathways in module 1 were mainly associated with systemic lupus erythematosus, histone deacetylases (HDACs), and alcoholism. Interleukin -10, -4, -13 signaling, and overview of proinflammatory and profibrotic mediators were involved in module 2. The MCODE3 module included TNFR2 noncanonical NF-kB pathway, TNF receptor superfamily (TNFSF) members mediating noncanonical NF-kB pathway, and NF-kappa B signaling pathway. The MCODE4 module included EGFR downregulation, signaling by EGFR, and branch elongation of an epithelium. Finally, the MCODE5 module included protein ubiquitination and protein modification by small protein conjugation (Figure 2(b)).
3.3. KEGG and GO Analyses of MPXV Infection
To determine the functional role of MPXV infection, the KEGG and GO analyses of 212 DEGs of GSE36854 and GSE21001 were carried out with the online tool SHINEY GO (https://bioinformatics.sdstate.edu/go/). The network of GO BP of 211 DEGs from the MPXV infection using Shiny Go was mainly involved in the regulation of multicellular organism development, regulation of phosphate metabolic process, regulation of cell differentiation, and positive regulation (Figure 3(a)). As shown in Figure 3(b), KEGG analysis of 212 DEGs from GSE36854 and GSE21001 showed involvement in coronavirus disease, cytokine-cytokine receptor interaction, inflammatory bowel disease, Chagas disease, lipid and atherosclerosis, T cell receptor signaling pathway, C-type lectin receptor signaling pathway, rheumatoid arthritis, and malaria.
3.4. Predicted Drug-Target Relationship Using 212 DEGs of Monkeypox Infection
We used the drug matador from the SHINEY GO online tool to predict drug-target relationships using 212 DEGs of GSE36854 and GSE21001 for MPXV infection. The top ten drugs related to monkeypox include troglitazone, fluticasone propionate, indomethacin, atorvastatin, nitroglycerin, cerivastatin, budesonide, 3-morpholinosydnonimine, quinapril, and valsartan (Figure 3(c)).
3.5. DisGNET, PaGenBase, and TRRUST of 212 DEGs from MPXV Infection
DisGeNET database disease enrichment analysis by using the Metascape tool demonstrated that the hub genes of 212 DEGs from GSE36854 and GSE21001 were associated with proliferative vitreoretinopathy, pneumonitis, juvenile arthritis, anemia, airway disease, and so on (Figure 4(a)). PaGenBase was used to perform a tissue characteristic enrichment analysis with the Metascape tool and showed that the hub genes of 212 DEGs from GSE36854 and GSE21001 were mainly enriched in the lung (Figure 4(b)). The TRRUST database analysis by using the Metascape tool demonstrated that RELA, SP1, and NF-κB1 are core transcription factors regulating the hub genes of 212 DEGs of MPXV infection (Figure 4(c)).
3.6. The Functional Analysis of MPXV and SARS-CoV-2 Infections
KEGG analysis of 212 DEGs from GSE36854 and GSE21001 for MPXV infection revealed COVID-19 relationship with COVID-19. Therefore, to determine the relationship between MPXV and SARS-CoV-2 infections, the Coronascape COVID database (https://metascape.org/COVID) was utilized. The biological function of coexpressed genes between SARS-CoV-2 infection and 212 DEGs from GSE36854 and GSE21001 for MPXV infection included GO1902532 negative regulation of intracellular signal transduction, GO0010942 positive regulation of cell death, hsa05202 transcriptional misregulation in cancer, GO0045596 negative regulation of cell differentiation, GO0043408 regulation of the MPAK cascade, hsa05417 NF-kappa B signaling pathway, GO 0009617 response to bacterium, and hsa04668 TNF signaling pathway (Figure 5(a)).
The network was visualized using Cytoscape5, wherein each node represents an enriched term and was colored first by its cluster ID and MCODE algorithm to see densely connected network components of genes coexpressed between SARS-CoV-2 and MPXV infections, as shown in Figure 5(b). The PPI interactions of the coexpressed genes were subjected to pathway and process enrichment analyses. The network of enriched term by clusters is involved in cytokine signaling in the immune system, TNF signaling, response to growth factor, regulation of MAPK cascade, positive regulation of cellular components, COVID-19, and so on (Figure 5(b)).
3.7. Key Transcription Factors and Associated Human Disease by Comparing between MPXV and SARS-CoV-2 Infections in the Cells
The DisGeNET discovery platform (https://www.disgenet.org) provided the Coronascape, a COVID-19 database (https://metascape.org/COVID), to determine human diseases associated with COVID-19 and monkeypox. We used 212 DEGs from GSE36854 and GSE21001 to compare the MPXV and SARS-CoV-2 infections. DisGeNET database disease enrichment analysis of our data (212 DEGs from GSE36854 and GSE21001 for MPXV infection) and COVID-19 database demonstrated that the hub genes were associated with lung disease, reperfusion injury, anoxia, dermatitis, pneumonitis, inflammatory dermatosis, and others (Figure 6(a)). Analysis of the transcription factors using the TRRUST database from our data (212 DEGs from GSE36854 and GSE21001) for MPXV infection and COVID-19 database by Metascape identified JUN, REL, STAT3, NFKB1, RELA, SP1, and so on (Figure 6(b)).
4. Discussion
The outbreak of monkeypox has spread across several countries in 2022 after the COVID-19 pandemic. However, there is no information of the function of monkeypox and its relation with COVID-19. In the present study, the function of MPXV infection was identified by bioinformatics analyses. In addition, key pathways between SARS-CoV-2 and monkeypox infections were suggested using the bioinformatics analyses tool.
Although SARS-CoV-2 and MPXV are both enveloped viruses that replicate in the cytoplasm and require host machinery to replicate and make new viral proteins, there are distinctions between the two viruses [Reference Roushdy23]. MPXV is a double-stranded DNA virus, whereas SARS-CoV-2 is a single-stranded RNA virus. Within the last two years, SARS-CoV-2 has several different clades being mutated from its alpha form to beta, gamma, delta, and omicron clades [Reference Roushdy23, Reference Ding, Jiang, Xiong and Lei36]. MPXV has two clades such the Congo Basin and Nigerian clades [Reference Kumar, Acharya, Gendelman and Byrareddy6]. MPXV is transmitted primarily through very close physical contact from human-human transmission and long-term contact with an infected animal, whereas COVID-19 illness is caused by droplet, direct contact, and airborne transmission of SARS-CoV-2 [Reference Roushdy23]. The SARS-CoV-2 pandemic spreads in the form of tides and waves, whereas MPXV spreads more linearly [Reference Roushdy23]. MPXV enters cells via macropinocytosis, whereas SARS-CoV-2 penetrates the human body by attaching its spikes to endothelial cells angiotensin converting enzyme 2 receptors (ACE-2) [Reference Roushdy23].
In this study, the hub genes of 212 DEGs of MPXV infection by cytoHubba analysis were involved in cytokine-related genes such IL-6, IL-10, TNF, CCL2, ICAM1, CXCL8, and HDAC-related genes. In addition, the MCODE analysis also associated with systemic lupus erythematosus, HDACs, and IL‐10 and ‐4 signaling. HDACs and histone acetyl transferases play major roles in cell survival, growth, and proliferation as well as the regulation of gene transcription [Reference Wang, Fu, Mani, Wadler, Senderowicz and Pestell37]. IL-10 plays an important role as an anti-inflammatory cytokine by preventing pathogenic infection and excessive immune response [Reference Ni, Zhang and Yang38, Reference Verma, Balakrishnan and Sharma39], and IL-4 regulates the immune responses [Reference Nelms, Keegan, Zamorano, Ryan and Paul40]. Systemic lupus erythematosus is well known as a chronic autoimmune disease that produces autoantibodies and affects multiple organs [Reference Pan, Lu, Wang, Xu and Yang41].
In the present study, KEGG analysis of 212 DEGs of GSE36854 and GSE21001 for MPXV infection showed its involvement in COVID-19, cytokine-cytokine receptor interaction, and TNF signaling. In addition, the biological function of coexpressed genes between SARS-CoV-2 and MPXV infections in the cell lines by Metascape/COVID-19 analysis was enriched in cytokine signaling in the immune system, TNF signaling, response to bacterium, response to growth factor, and regulation of the MAPK cascade. There is evidence that cytokine modulation is associated with the severity of the monkeypox disease in humans [Reference Johnston, Johnson and Stonier42]. MPXV infection in Hela cells involves genes for proinflammatory cytokine and leukocyte chemotaxis [Reference Bourquain, Dabrowski and Nitsche28] and that in MK2 cells enhanced the IL-8 gene [Reference Alkhalil, Hammamieh, Hardick, Ichou, Jett and Ibrahim27, Reference Singh, Khan and Dinner43]. Many studies have suggested that SARS-CoV-2 is associated with a dysregulated inflammatory response, such as cytokine storm and inflammatory molecules [Reference Mulchandani, Lyngdoh and Kakkar44, Reference Blanco-Melo, Nilsson-Payant and Liu45]. Therefore, monkeypox and COVID-19 might be associated with the inflammatory immune response.
A previous study demonstrated that transcription factors play an important role in the regulation of transcription, metabolism, and immune response [Reference Singh, Khan and Dinner43, Reference Desvergne, Michalik and Wahli46], and they are closely associated with diseases [Reference Lee and Young47]. In this study, JUN, REL, NFKB1, and STAT3 are mainly involved transcriptional factor according to the TRRUST analysis when comparing MPXV and SARS-CoV-2 infections. There is evidence that JUN, REL, NFKB1, and STAT3 are involved in inflammation [Reference Cartwright, Perkins and L Wilson48–Reference Fullard, Wilson and Oakley51].
5. Conclusions
In conclusion, the hub genes based on the cytoHubba and MCODE analyses of MPXV infection included cytokine and HDAC-related genes. KEGG analysis of MPXV infection showed involvement of MPXV in COVID-19, cytokine-cytokine receptor interaction, and TNF signaling pathway. Our study data suggest that the function of monkeypox is similar to that COVID-19 in terms of cytokine signaling in the immune system and TNF signaling. Our study provides pathological insight of monkeypox and molecular connections between COVID-19 and monkeypox. Our study has some limitations due to the lack of experimental verification of computational data. It will inevitably need in-depth research in the future.
Data Availability
Data used in this study are available at the NCBI (National Center for Biotechnology Information Gene Expression Omnibus), https://www.ncbi.nlm.nih.gov/geo/.
Conflicts of Interest
The author declares that there are no conflicts of interest.
Acknowledgments
This study was supported by the National Research Foundation of Korea, funded by the Ministry of Education (grant no. NRF-2021R1I1A1A01052609).