Comparison of Transcriptomic Signatures between Monkeypox-Infected Monkey and Human Cell Lines

Abstract

Monkeypox virus (MPV) is a smallpox-like virus belonging to the genus Orthopoxvirus of the family Poxviridae. Unlike smallpox with no animal reservoir identified and patients suffering from milder symptoms with less mortality, several animals were confirmed to serve as natural hosts of MPV. The reemergence of a recently reported monkeypox epidemic outbreak in nonendemic countries has raised concerns about a global outburst. Since the underlying mechanism of animal-to-human transmission remains largely unknown, comprehensive analyses to discover principal differences in gene signatures during disease progression have become ever more critical. In this study, two MPV-infected in vitro models, including human immortal epithelial cancer (HeLa) cells and rhesus monkey (Macaca mulatta) kidney epithelial (MK2) cells, were chosen as the two subjects to identify alterations in gene expression profiles, together with co-regulated genes and pathways that are affected during monkeypox disease progression. Using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and MetaCore analyses, we discovered that elevated expression of genes associated with interleukins (ILs), G protein-coupled receptors (GPCRs), heat shock proteins (HSPs), Toll-like receptors (TLRs), and metabolic-related pathways play major roles in disease progression of both monkeypox-infected monkey MK2 and human HeLa cell lines. Interestingly, our analytical results also revealed that a cluster of differentiation 40 (CD40), plasmin, and histamine served as major regulators in the monkeypox-infected monkey MK2 cell line model, while interferons (IFNs), macrophages, and neutrophil-related signaling pathways dominated the monkeypox-infected human HeLa cell line model. Among immune pathways of interest, apart from traditional monkeypox-regulated signaling pathways such as nuclear factor- (NF-κB), mitogen-activated protein kinases (MAPKs), and tumor necrosis factors (TNFs), we also identified highly significantly expressed genes in both monkey and human models that played pivotal roles during the progression of monkeypox infection, including CXCL1, TNFAIP3, BIRC3, IL6, CCL2, ZC3H12A, IL11, CSF2, LIF, PTX3, IER3, EGR1, ADORA2A, and DUOX1, together with several epigenetic regulators, such as histone cluster family gene members, HIST1H3D, HIST1H2BJ, etc. These findings might contribute to specific underlying mechanisms related to the pathophysiology and provide suggestions regarding modes of transmission, post-infectious sequelae, and vaccine development for monkeypox in the future.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Monkeypox virus (MPV/MPXV) is commonly known as one of the rare systemic infections caused by one specie of the genus Orthopoxvirus in the family Poxviridae [1, 2]. Other than the smallpox virus (variola virus (VARV)), the best-known species of this genus, monkeypox virus, together with cowpox virus (CPXV), horsepox virus (HSPV), raccoonpox virus (RCNV), camelpox virus (CMLV), mousepox virus (also referred to as Ectromelia virus (ECTV)), and Alaskapox (AKPV) virus, are species believed to infect healthy humans through different modes of zoonotic transmission [3–5]. After the eradication of smallpox, a result of universal vaccination, almost no natural outbreaks related to Orthopoxvirus have been recorded until recently, after a significant number of confirmed monkeypox cases were reported worldwide [6]. The first discovery of this virus in monkeys in the lab (in 1958), followed by the first human infection detected in Africa (in 1970), meant that monkeypox was no longer endemic in animals as previously assumed; at least 400 human cases scattered across the western and central parts of the African continent have been diagnosed with similar clinical characteristics to smallpox but causing milder symptoms and being less lethal [7–9]. As of May 2022, multiple cases of MPV have been documented in at least 11 regions outside of Africa, including the European Zone, Britain, Middle East, Australia, the North American continent (including the United States, Canada, and South America), resulting in unprecedented outbreaks in those developed countries [8–10]. Further, public health investigations revealed the approximately 20-fold increase in incidences that resulted from the discontinuation of smallpox vaccination since 1980, while the probable source of infection could be either an animal reservoir (mainly rodents) or cross-infection [11, 12]. Transmission through person-to-person contact is claimed to primarily occur through close and prolonged contact with fluids from skin lesions, respiratory droplets, fomites of infected people, and even during sexual intercourse among populations of men who have sex with men (MSM) [13, 14].

To date, definite diagnoses of monkeypox are largely based on either a polymerase chain reaction (PCR) of skin lesion specimens or immunohistochemical (IHC) staining of biopsies taken from those infected, as misdiagnoses of other blister-like diseases such as smallpox and chickenpox may occur in cases where only a visual skin examination and/or dermoscopy are solely employed [15, 16]. Although the genomes of several variants of the MPV have been fully sequenced, there are limited reports on the underlying mechanisms of viral transmission and pathogenic pathways of this disease [17, 18]. Moreover, proven treatment and preventive vaccines specifically designed for monkeypox have not yet been formally approved by the agencies, thus clinical care only focuses on symptomatic treatment. Therefore, developing a research design that focuses on the underlying molecular mechanisms during disease progression would greatly contribute to improving the situation as it provides necessary background knowledge for subsequent studies, such as diagnostic test development, epidemiological tracking, designs of vaccines and antiviral drugs, and other related adjuvant treatments.

In this study, we attempted to provide a concise and informative overview of principal differences in gene expression profiles between two MPV-infected in vitro models, including human immortal epithelial cancer (HeLa) cells and rhesus monkey (Macaca mulatta) kidney epithelial (MK2) cells, together with co-regulated genes and pathways that are affected during disease progression. The Gene Expression Omnibus (GEO) database and bioinformatics approaches were employed to perform analyses of available samples related to these two models (Figure 1). We also revealed several critical regulatory downstream networks and investigated the functions of potential genes to determine whether they can serve as clinical biomarkers for MPV. These findings may contribute to an understanding of specific underlying mechanisms related to the pathophysiology and provide suggestions regarding modes of transmission, post-infection sequelae, and vaccine development for monkeypox in the future.

[figure(s) omitted; refer to PDF]

2. Materials and Methods

2.1. Bioinformatics Application for Data Acquisition and Processing

In this study, we investigated an in vitro system using gene expression data of two separate subjects, including human immortal epithelial (HeLa) cells and monkey Macaca mulatta kidney epithelial (MK2) cells to characterize the effects of physical monkeypox viral infection. Relevant microarray gene expression datasets from the GEO (http://www.ncbi.nlm.nih.gov/geo)—an international public functional genomics data repository affiliated with the National Center for Biotechnology Information (NCBI; https://www.ncbi.nlm.nih.gov/)—were retrieved for subsequent analyses. In particular, MPV-induced alterations in transcriptomic profiles of 7-hour monkeypox-infected Macaca mulatta kidney epithelial (MK2) cells were obtained from the GSE21001 dataset [19], whereas similar information of 6-hour monkeypox-infected human immortal epithelial (HeLa) cells was acquired from the GSE36854 dataset [20]. The following analyses of both studies are based on comparison of each subject with its corresponding mock-infected counterpart as control.

Before being mapped onto Ensembl features using packages attached to the biomaRt tool (vers. 2.26.1), the CLC Genomics Workbench vers. 10.1 (CLC bio, Aarhus, Denmark) was utilized to process and standardize the raw data [21–23]. Signals were subsequently presented as a clustered heatmap based on messenger (m)RNA expression patterns using pheatmap (vers. 1.0.12) in the R environment [24–26], as well as the online platform available at http://www.bioinformatics.com.cn/srplot. For subsequent discoveries of enriched functionally related gene groups, the Database for Annotation, Visualization, and Integrated Discovery (DAVID), which is maintained by the Laboratory of Human Retrovirology and Immunoinformatic (LHRI) [27–29], was chosen as a powerful tool by the various functions it offers, especially in term of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). As previously described [30–32], the top 5% of upregulated genes derived from monkeypox-infected groups were further processed as input data for the GO analyses [33–35], while a $p$ value of <0.05 was set as the significance level. After biological networks and related pathways were constructed, a cut-off value of 0.05 was applied to select significantly enriched pathways or groups of annotated genes.

2.2. Pathway-Based Analyses and Network Enrichment Analyses

MetaCore (Enrichment Analysis Workflow and analysis network; GeneGo, St. Joseph, MN, USA) was designed to identify biological processes related to gene microarray data. The top 5% of upregulated genes with substantial differences in transcriptome levels were processed as input to the MetaCore program to compare average levels of gene expressions in the two models. Signal transduction pathways were examined. A statistically significant difference was indicated by a $p$ value of <0.05.

2.3. Construction of a Protein-Protein-Interacting (PPI) Network using the STRING Analysis

The STRING database (vers. 11.0) (https://string-db.org/) serves as a search engine and a protein information resource that provides a considerable number of proteins and known interaction activities upon which differentially expressed genes (DEGs) were examined [36–38]. The Cytoscape stringApp was also employed in this work to construct a PPI-network, while protein cluster classification was built using the $k$-means clustering algorithm.

3. Results

3.1. GO Analysis of Monkeypox-Infected Models

For the monkey cell line model, transcriptomic data of the 7-hour monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell retrieved from the GSE21001 dataset were matched with a mock treatment group for a comparative analysis of DEGs. GO analytical results associated with monkeypox infection are depicted in Figure 2(a). Pathways assigned to cluster 1 (c1) contained the most significant highly expressed genes and pathways in the monkeypox-infected group, whereas cluster 2 (c2) exerted an opposite trend. The GO analysis suggested that c1 was related to inflammatory responses, chemokine activities, immune system processes, responses to peptidoglycans, cellular responses to lipopolysaccharide, nucleosome assembly, chromatin assembly or disassembly, responses to ionizing radiation, and responses to temperature stimuli. Pathways within c2 were associated with kinase binding, thymus development, responses to endoplasmic reticular stress, transfer (t)RNA aminoacylation for protein translation, and responses to amino acid stimulation.

[figure(s) omitted; refer to PDF]

For the 6-hour monkeypox-infected human cell line model, transcriptomic data of monkeypox-infected human HeLa cells were retrieved from the GSE36854 dataset. Pathways assigned to cluster 1 (c1) were those most highly correlated with monkeypox infection, whereas cluster 2 (c2) exerted an opposite trend. The GO analysis suggested that pathways within c1 were related to growth factor activities, leukocyte activation, responses to organic cyclic compounds, inflammatory responses, regulation of the extracellular signal-regulated kinase 1 (ERK1) and the extracellular signal-regulated kinase 2 (ERK2) cascade, chemokine activities, and regulation of nitric oxide biosynthetic processes. Pathways in c2 were related to nucleosome assembly, protein heterodimerization activities, telomere maintenance, extracellular regions, the actin cytoskeleton, responses to protozoans, and coagulation.

3.2. Pathway Analysis of the Two Monkeypox-Infected Models

Venn diagrams were employed to filter upregulated genes in common between the two models mentioned above. The top 5% of messenger (m)RNAs within the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model and the human immortal epithelial cancer (HeLa) counterpart were highly expressed compared to those within the mock control (Figure 3(a)). In total, 149 overlapping genes that were highly expressed in both monkey and human models were subsequently input to the MetaCore tool for a map analysis. Several standard maps were found to be related to monkeypox infection, including the “immune response_IL-1 signaling pathway,” “NF-κB pathway in multiple myeloma,” “immune response_CD40 signaling in B cells,” “TNF-alpha-induced inflammation,” “signaling in normal and asthmatic airway epithelium,” “immune response_plasmin signaling,” “glomerular injury in lupus nephritis,” “immune response_IL-33 signaling pathway,” “macrophage and dendritic cell phenotype shift in cancer,” “renal tubulointerstitial injury in lupus nephritis,” “signal transduction_nonapoptotic FasR (CD95) signaling,” and “inflammatory mechanisms of pancreatic cancerogenesis” (Figure 3(b), Supplementary Table A1). Meanwhile, Endogenous Metabolic Networks corresponding to genes from the shortlist of upregulated genes shared by the two models revealed that pathways were associated with phospholipid biosynthesis, including the “Ceramide pathway,” “N-acyl-sphingosine phosphate pathway,” “2-arachidonoyl-glycerol 3-phosphocholine pathway,” “phosphatidylethanolamine pathway,” “(L)-valine pathways and transport,” “carbohydrate metabolism_TCA and tricarboxylic acid transport,” and “vitamin, mediator and cofactor metabolism_alpha-tocotrienol,” which are metabolic signalling pathways dominated by these genes within both aforementioned models (Figures 3(c)–3(e)), Supplementary Table A2).

[figure(s) omitted; refer to PDF]

To further understand biological processes in the individual models, we analyzed pathways related to monkeypox infection in each case. For the 7-hour monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model, the top 5% of upregulated genes obtained compared to the mock treatment group were input into the MetaCore platform for a comprehensive map analysis. The “immune response_IL-1 signaling pathway,” “immune response_IL-33 signaling pathway,” “signal transduction_nonapoptotic FasR (CD95) signaling,” “immune response_histamine H1 receptor signaling in immune response,” “IL-1 signaling in melanoma,” “inflammatory mechanisms of pancreatic cancerogenesis,” “immune response_HSP60 and HSP70/TLR signaling pathway,” “G protein-coupled receptors signaling in lung cancer,” “TNF-alpha-induced inflammatory signaling in normal and asthmatic airway epithelium,” and “resistance of pancreatic cancer cells to death receptor signaling” were among the most significant maps detected (Supplementary Figure A1, Table A3). Meanwhile, the endogenous Metabolic Networks corresponding to genes from the shortlist of upregulated genes within the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cells revealed that pathways associated with phospholipid biosynthesis and amino acid synthesis included the “phosphatidylethanolamine pathway,” “N-acyl-sphingosine phosphate pathway,” “lyso-phosphatidylserine pathway,” “1-alkyl-glycerol_3-phosphoethanolamine pathway,” “1-oleoyl-glycerol_3-phosphate pathway,” “2-oleoyl-glycerol_3-phosphate pathway,” and “aminoacid metabolism_asparagine, and aspartic acid metabolism and transport,” which are metabolic signaling pathways dominated by these genes within this model (Supplementary Figure A2, Table A4).

For the human immortal epithelial cancer (HeLa) cell model, the top 5% of upregulated genes compared to the mock treatment groups were input into the MetaCore platform for a comprehensive map analysis. The “immune response_IL-1 signaling pathway,” “inflammatory mechanisms of pancreatic cancerogenesis,” “G protein-coupled receptors signaling in lung cancer,” “immune response_plasmin signaling,” “immune response_CD40 signaling in dendritic cells, monocytes, and macrophages,” “interleukins-induced inflammatory response in asthmatic airway fibroblasts,” “TNF-alpha-induced inflammatory signaling in normal and asthmatic airway epithelium,” “NF-κB-, AP-1-, and MAPKs-mediated proinflammatory cytokine production by eosinophils in asthma,” “immune response_CD40 signaling in B cells,” and “immune response_IL-17 signaling pathways” were among the most significant maps detected (Supplementary Figure A3, Table A5). Through the Metabolic Networks (Endogenous) corresponding to genes from the shortlist of upregulated genes shared by the monkeypox-infected human immortal epithelial cancer (HeLa) cells, data revealed that pathways associated with phospholipid biosynthesis included the “phosphatidylcholine pathway,” “1-docosahexaenoyl-glycerol_3-phosphocholine pathway,” “lysophosphatidic acid pathway,” “1-icosatrienoyl-sn-glycero-3-phosphocholine pathway,” “2-arachidonoyl-glycerol_3-phosphocholine pathway,” and “GalNAcbeta1-3Gal pathway,” which are metabolic signaling pathways dominated by these genes within this model (Supplementary Figure A4, Table A6).

In addition to the pathways and network analysis for leading-edge genes generated by MetaCore as mentioned before, we also used the KEGG database to validate these data. The top enriched KEGG pathways were analyzed for common genes shared by the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model and the human immortal epithelial cancer (HeLa) cell model, and pathways were ranked by $p$ values (Supplementary Figure A5).

3.3. Analytical Results of Co-regulated Interactions and Cellular Component Annotations in the Monkeypox-Infected Models

To further explore the most significance of co-upregulated genes from both models, we used a Venn diagram to filter the 24 genes (among the top 1% of genes) in common, which were highly expressed in both models (Figure 4(a)). Upregulated genes in monkeypox-infected cells included CXCL1, TNFAIP3, BIRC3, IL-6, CCL2, ZC3H12A, IL-11, CSF2, LIF, PTX3, IER3, EGR1, ADORA2A, DUOX1, and histone clusters family members, such as HIST1H3D, HIST1H2BJ, HIST1H2AK, HIST1H2AD, HIST1H2AC, HIST1H1B, HIST2H2AB, HIST1H2BM, HIST1H2BH, and HIST1H2BB (Figure 4(b)). These 24 genes were further imported into the STRING platform to construct PPI networks, which may play pivotal roles in manifestations after monkeypox infection (Figure 4(c)).

[figure(s) omitted; refer to PDF]

Upregulated genes shared by these two models were also imported into the GO platform for cellular component annotations under monkeypox infection circumstances. In the case of the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model, these genes mostly functioned in apical plasma membranes, neuronal cell bodies, the cytosol, nucleoplasm cytoplasm, extracellular exosomes, an integral component of presynaptic membranes, postsynaptic membranes, cytoplasmic exosomes, and the perikaryon (Figure 5(a)). The same analytical flowchart was applied to the monkeypox-infected human immortal epithelial cancer (HeLa) cell model to identify cellular components under monkeypox infection circumstances. Genes with the greatest significance functioned in the extracellular space, receptor complexes, integral component of membranes, caveolae, external side of plasma membranes, extracellular matrices, voltage-gated calcium channel complexes, presynaptic active zones, glutamatergic synapses, and dendritic shafts (Figure 5(b)). From such observations gained from the two in vitro monkeypox-infected models, we concluded that certain cellular mechanisms of monkeypox infection are almost identical but slightly different when comparing the two species (Figure 5(c)).

[figure(s) omitted; refer to PDF]

4. Discussion

In most cases, viral infections only occur when a virus that causes disease successfully invades a host, beginning by penetrating the antiviral defenses, inserting its genetic materials, taking control of host cells’ machinery to rapidly create a huge number of viral replications, and ultimately becoming transmissible between individuals. In the context of the world having gone through many epidemic zoonotic outbreaks that originated from animals [39–41], including the human immunodeficiency virus (HIV; which jumped from chimpanzees), influenza (from birds and pigs), bovine spongiform encephalopathy (BSE; from cows), Ebola (from bats), and at least three zoonotic coronavirus diseases in chronological order, including the severe acute respiratory syndrome coronavirus (SARS-CoV), the Middle Eastern respiratory syndrome coronavirus (MERS-CoV), and the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Therefore, gaining a comprehensive understanding of the mechanisms of infection in both human and non-human hosts is urgently needed [42–45].

In attempts to prevent further outbreaks and deaths from these diseases, worldwide efforts to discover rapid early-detection methods, effective antiviral drugs, and preventive vaccines have been increasing over the years [46, 47]. Other than conventional approaches, novel and advantageous technologies, such as high-throughput sequencing, searchable database repositories of gene expression data, and assistance from web-based or computer-based bioinformatics tools, have been employed to elucidate underlying molecular mechanisms during viral infection and explore how different host cells respond via performing analyses of alterations in gene expressions [48–54].

In the case of monkeypox infection, according to previous studies, a considerable number of immunomodulation-related cytokines, including interleukin-1 (IL-1), tumor necrosis factors (TNFs), and interferons (IFNs), were reported to be involved in disease progression. The ILs superfamily (including IL-1, IL-18, and IL-33) and TNF signaling-mediated immune responses are important for antiviral activity [55, 56]; however, the roles of these signaling pathways in monkeypox-infected somatic cells remain unclear. Interestingly, in case of other zoonotic viral diseases mentioned before, numerous cellular and molecular immunomodulatory pathways were also proven to play crucial roles. In particular, the G protein-coupled receptors (GPCRs) not only helps stimulate internalization during infection with the influenza A virus [57] but also takes part in Epstein-Barr virus- (EBV-) mediated immunosuppression and oncogenesis [58]. Studies on the molecular chaperonin HSP60 confirmed its interactions with Ebola virus infection [59] and human hepatitis B virus polymerase [60], while HSP70 plays a crucial role in mediating the progression of Zika virus in infected host cells [61]. Histamine, on the other hand, contributes to severe pneumonia in pigs infected with the H1N1 influenza virus [62], whereas levels of histamine and leukotrienes presented in acute dengue patients are closely related to disease severity [63].

Our research findings were consistent with previous studies on monkeypox infection as mentioned above; the GO and pathway analyses revealed similar enriched pathways for both monkeypox-infected monkey and human models, including IL-1, IL-33, IL-17, IL-18, NF-κB, MAPKs, and TNF-R2 signaling (Figure 3). However, interestingly, we also discovered relationships between monkeypox and GPCRs, HSP60/70, histamine, plasmin, and histone cluster-related signaling (Supplementary Figures A1, A2) that have rarely been reported before.

Regarding 24 common genes that were highly expressed in both in vitro models (Figure 4), when compared to the previous literature, similar expression patterns were also observed in case of CXCL1, TNFAIP3, BIRC3, IL6, CCL2, ZC3H12A, IL11, CSF2, LIF, PTX3, IER3, EGR1, ADORA2A, DUOX1, and several histone cluster family members such as HIST1H3D and HIST1H2BJ. According to a recent study, there were dramatic increases in serum concentrations of IL1-RA, IL-6, and IFN-γ in cynomolgus macaques after being infected with aerosolized MPV [64]. In myeloid cells, TNFAIP3 deficiency helps protect host cells against invasion by the influenza A virus [65]. Both ZC3H12A/MCPIP1 and CDKN1A/p21 may contribute to the synergistic effect on replication control and pharmacological manipulation of HIV-1 [66]. Infection with SARS-CoV-2 and several influenza viruses can lead to alterations in CSF2 and IL-11 expressions [67–69]. PTX3, the expression of which was reported to be upregulated by TNF-α during acute lung injury, was recently found to serve as a reliable prognostic indicator in predicting short-term mortality from SARS-CoV-2 [70]. Apoptosis and impacted viral replication were found to be induced by infection with the Venezuelan equine encephalitis virus; however, knockdown of early growth response 1 (EGR1) significantly inhibited these signaling pathways [71]. Two members of the NADPH oxidase (NOX) family, including the DUOX1 and DUOX2 proteins, are potential therapeutic targets against infectious diseases caused by the influenza A virus [72]. Previous research demonstrated that DNA viruses exploit host cellular epigenetic processes to their advantage during infection [73, 74]. Most interestingly, several histone cluster family members were reported to be involved in differential responses of human fetal brain neural stem cells to Zika virus infection [75]. For example, three members of the histone cluster 1 H2B (HIST1H2B) family, namely HIST1H2BB, HIST1H2BK, and HIST1H2BO, were among those hub genes and pathways that contributed to Zika virus infection [76]. However, only a limited number of studies have described such relationships in terms of monkeypox. Interestingly, results of our bioinformatics analysis revealed that expression level of certain histone cluster family members, such as HIST1H3D, HIST1H2BJ, HIST1H2AK, HIST1H2AD, HIST1H2AC, HIST1H1B, HIST2H2AB, HIST1H2BM, HIST1H2BH, and HIST1H2BB, were elevated in both monkeypox-infected cell models of monkey and human. These data suggested that some of these epigenetic regulators may also play important roles in monkeypox infection (Figure 4). Cross-validation performed in the 6-hour monkeypox-infected human immortal epithelial cancer (HeLa) cell model (GSE24125) [77] revealed similar patterns of elevated expressions of the aforementioned genes (Supplementary Figure A6). As the most challenging aspect of diagnosing monkeypox infection is distinguishing monkeypox vesicular rash eruption from those of other diseases [78], our study provides a prospect of the histone cluster family as potential regulators correlated with monkeypox infection. Since the current study only focused on signaling pathways regulated by DEGs caused by monkeypox infection, further investigation into the underlying pathogenesis is required for consolidated guidelines for monkeypox treatment in the future.

5. Conclusion

In the current study, we focused on investigating genetic signatures and their related pathways in two different monkeypox-infected models. Our research findings revealed the critical roles of multiple genes and their regulatory pathways in both monkeypox-infected models, which not only directly in line with previous literature but also provided novel functions as well as highly expressed genes related to MPV infection that had been less well revealed in the past. In the context of widespread epidemics of viral infectious diseases that are prevalent today, this work may contribute to bridging the traditional gap between bench research and clinical applications. However, with attempts to validate and examine the potentials of these genes as novel antiviral therapies, further expanding research on disease basics may provide more-concrete perspectives.

Authors’ Contributions

T.M.X.D., I.J.Y., C.Y.W., and M.C.Y. designed the studies; C.C.W., J.Z.J., C.Y.S., H.L.L., and Z.S. collected the data; T.M.X.D., C.C.W., C.Y.S., and H.L.L. performed the bioinformatics analysis; H.D.K.T., G.A., and S.C.K. provided biostatistics analyses; T.M.X.D., I.J.Y., C.C.C., and C.Y.W. wrote the initial manuscript draft; M.C.Y. edited the manuscript, and all authors read and approved the final manuscript. Do Thi Minh Xuan and I-Jeng Yeh contributed equally to this work.

References

[1] E. Tambo, A. M. Al-Nazawi, "Combating the global spread of poverty-related monkeypox outbreaks and beyond," Infectious Diseases of Poverty, vol. 11 no. 1,DOI: 10.1186/s40249-022-01004-9, 2022.

[2] S. Essbauer, M. Pfeffer, H. Meyer, "Zoonotic poxviruses," Veterinary Microbiology, vol. 140 no. 3-4, pp. 229-236, DOI: 10.1016/j.vetmic.2009.08.026, 2010.

[3] S. N. Shchelkunov, A. V. Totmenin, I. V. Babkin, P. F. Safronov, O. I. Ryazankina, N. A. Petrov, V. V. Gutorov, E. A. Uvarova, M. V. Mikheev, J. R. Sisler, J. J. Esposito, P. B. Jahrling, B. Moss, L. S. Sandakhchiev, "Human monkeypox and smallpox viruses: genomic comparison," Human monkeypox and smallpox viruses: genomic comparison, vol. 509 no. 1, pp. 66-70, DOI: 10.1016/S0014-5793(01)03144-1, 2001.

[4] S. Parker, D. A. Schultz, H. Meyer, R. M. Buller, "Smallpox and Monkeypox Viruses☆," in Reference Module in Biomedical Sciences, 2014.

[5] D. B. Di Giulio, P. B. Eckburg, "Human monkeypox: an emerging zoonosis," The Lancet Infectious Diseases, vol. 4 no. 1, pp. 15-25, DOI: 10.1016/S1473-3099(03)00856-9, 2004.

[6] M. Kozlov, "Monkeypox goes global: why scientists are on alert," Nature, vol. 606 no. 7912, pp. 15-16, DOI: 10.1038/d41586-022-01421-8, 2022.

[7] J. G. Breman, R. Kalisa, M. V. Steniowski, E. Zanotto, A. I. Gromyko, I. Arita, "Human monkeypox, 1970-79," Bulletin of the World Health Organization, vol. 58 no. 2, pp. 165-182, 1980.

[8] H. Meyer, M. Perrichot, M. Stemmler, P. Emmerich, H. Schmitz, F. Varaine, R. Shungu, F. Tshioko, P. Formenty, "Outbreaks of disease suspected of being due to human monkeypox virus infection in the Democratic Republic of Congo in 2001," Journal of Clinical Microbiology, vol. 40 no. 8, pp. 2919-2921, DOI: 10.1128/JCM.40.8.2919-2921.2002, 2002.

[9] P. Y. Nguyen, W. S. Ajisegiri, V. Costantino, A. A. Chughtai, C. R. MacIntyre, "Reemergence of human monkeypox and declining population immunity in the context of urbanization, Nigeria, 2017-2020," Emerging Infectious Diseases, vol. 27 no. 4,DOI: 10.3201/eid2704.203569, 2021.

[10] R. Vivancos, C. Anderson, P. Blomquist, S. Balasegaram, A. Bell, L. Bishop, C. S. Brown, Y. Chow, O. Edeghere, I. Florence, S. Logan, P. Manley, W. Crowe, A. McAuley, A. G. Shankar, B. Mora-Peris, K. Paranthaman, M. Prochazka, C. Ryan, D. Simons, R. Vipond, C. Byers, N. A. Watkins, W. Welfare, E. Whittaker, C. Dewsnap, A. Wilson, Y. Young, M. Chand, S. Riley, S. Hopkins, "Community transmission of monkeypox in the United Kingdom, April to May 2022," Eurosurveillance, vol. 27 no. 22,DOI: 10.2807/1560-7917.ES.2022.27.22.2200422, 2022.

[11] B. Titanji, B. Tegomoh, S. Nematollahi, M. Konomos, P. A. Kulkarni, "Monkeypox-a contemporary review for healthcare professionals," in Open Forum Infectious Diseases, vol. 9 no. 7,DOI: 10.1093/ofid/ofac310, 2022.

[12] N. Kumar, A. Acharya, H. E. Gendelman, S. N. Byrareddy, "The 2022 outbreak and the pathobiology of the monkeypox virus," Journal of Autoimmunity, vol. 131,DOI: 10.1016/j.jaut.2022.102855, 2022.

[13] J. Heskin, A. Belfield, C. Milne, N. Brown, Y. Walters, C. Scott, M. Bracchi, L. S. Moore, N. Mughal, T. Rampling, A. Winston, M. Nelson, S. Duncan, R. Jones, D. A. Price, B. Mora-Peris, "Transmission of monkeypox virus through sexual contact - a novel route of infection," Journal of Infection, vol. 85 no. 3, pp. 334-363, DOI: 10.1016/j.jinf.2022.05.028, 2022.

[14] E. Petersen, A. Zumla, D. Hui, L. Blumberg, S. Valdoleiros, L. Amao, F. Ntoumi, D. Asogun, L. Simonsen, N. Haider, T. Traore, N. Kapata, O. Dar, J. Nachega, A. Abbara, A. al Balushi, R. Kock, M. Maeurer, S. S. Lee, D. R. Lucey, G. Ippolito, M. P. G. Koopmans, "Vaccination for monkeypox prevention in persons with high-risk sexual behaviours to control on-going outbreak of monkeypox virus clade 3," International Journal of Infectious Diseases, vol. 122, pp. 569-571, DOI: 10.1016/j.ijid.2022.06.047, 2022.

[15] W H Organization, "Surveillance, case investigation and contact tracing for monkeypox: interim guidance, 22 May 2022," World Health Organization, 2022.

[16] P. Dashraath, K. Nielsen-Saines, C. Mattar, D. Musso, P. Tambyah, D. Baud, "Guidelines for pregnant individuals with monkeypox virus exposure," The Lancet, vol. 400 no. 10345, pp. 21-22, DOI: 10.1016/S0140-6736(22)01063-7, 2022.

[17] J. R. Kugelman, S. C. Johnston, P. M. Mulembakani, N. Kisalu, M. S. Lee, G. Koroleva, S. E. McCarthy, M. C. Gestole, N. D. Wolfe, J. N. Fair, B. S. Schneider, L. L. Wright, J. Huggins, C. A. Whitehouse, E. O. Wemakoy, J. J. Muyembe-Tamfum, L. E. Hensley, G. F. Palacios, A. W. Rimoin, "Genomic variability of monkeypox virus among humans, Democratic Republic of the Congo," Emerging Infectious Diseases, vol. 20 no. 2, pp. 232-239, DOI: 10.3201/eid2002.130118, 2014.

[18] U. A. Awan, S. Riasat, W. Naeem, S. Kamran, A. A. Khattak, S. Khan, "Monkeypox: a new threat at our doorstep!," The Journal of Infection, vol. 85 no. 2, pp. e47-e48, DOI: 10.1016/j.jinf.2022.05.027, 2022.

[19] A. Alkhalil, R. Hammamieh, J. Hardick, M. A. Ichou, M. Jett, S. Ibrahim, "Gene expression profiling of monkeypox virus-infected cells reveals novel interfaces for host-virus interactions," Virology Journal, vol. 7 no. 1,DOI: 10.1186/1743-422X-7-173, 2010.

[20] D. Bourquain, P. W. Dabrowski, A. Nitsche, "Comparison of host cell gene expression in cowpox, monkeypox or vaccinia virus-infected cells reveals virus-specific regulation of immune response genes," Virology Journal, vol. 10 no. 1,DOI: 10.1186/1743-422X-10-61, 2013.

[21] C.-Y. Wang, Y.-C. Chang, Y.-L. Kuo, K.-T. Lee, P.-S. Chen, C. H. A. Cheung, C.-P. Chang, N. N. Phan, M.-R. Shen, H.-P. Hsu, "Mutation of the PTCH1 gene predicts recurrence of breast cancer," Scientific reports, vol. 9 no. 1, 2019.

[22] D. L. Cooke, D. B. McCoy, V. V. Halbach, S. W. Hetts, M. R. Amans, C. F. Dowd, R. T. Higashida, D. Lawson, J. Nelson, C. Y. Wang, H. Kim, Z. Werb, C. McCulloch, T. Hashimoto, H. Su, Z. Sun, "Endovascular biopsy: in vivo cerebral aneurysm endothelial cell sampling and gene expression analysis," Translational Stroke Research, vol. 9 no. 1, pp. 20-33, DOI: 10.1007/s12975-017-0560-4, 2018.

[23] T. Y. Weng, H. F. Wu, C. Y. Li, Y. H. Hung, Y. W. Chang, Y. L. Chen, H. P. Hsu, Y. H. Chen, C. Y. Wang, J. Y. Chang, M. D. Lai, "Homoharringtonine induced immune alteration for an efficient anti-tumor response in mouse models of non-small cell lung adenocarcinoma expressing Kras mutation," Scientific Reports, vol. 8 no. 1,DOI: 10.1038/s41598-018-26454-w, 2018.

[24] P. Shahi, C. Y. Wang, D. A. Lawson, E. M. Slorach, A. Lu, Y. Yu, M. D. Lai, H. Gonzalez Velozo, Z. Werb, "ZNF503/Zpo2drives aggressive breast cancer progression by down-regulation of GATA3 expression," Proceedings of the National Academy of Sciences of the United States of America, vol. 114 no. 12, pp. 3169-3174, DOI: 10.1073/pnas.1701690114, 2017.

[25] N. N. Phan, C. Y. Wang, Y. C. Lin, "The novel regulations of MEF2A, CAMKK2, CALM3, and TNNI3 in ventricular hypertrophy induced by arsenic exposure in rats," Toxicology, vol. 324, pp. 123-135, DOI: 10.1016/j.tox.2014.07.010, 2014.

[26] Y. H. Hung, Y. S. Chan, Y. S. Chang, K. T. Lee, H. P. Hsu, M. C. Yen, W. C. Chen, C. Y. Wang, M. D. Lai, "Fatty acid metabolic enzyme acyl-CoA thioesterase 8 promotes the development of hepatocellular carcinoma," Oncology Reports, vol. 31 no. 6, pp. 2797-2803, DOI: 10.3892/or.2014.3155, 2014.

[27] W. H. Da, B. T. Sherman, R. A. Lempicki, "Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources," Nature Protocols, vol. 4 no. 1, pp. 44-57, DOI: 10.1038/nprot.2008.211, 2009.

[28] W. H. Da, B. T. Sherman, R. A. Lempicki, "Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists," Nucleic Acids Research, vol. 37 no. 1,DOI: 10.1093/nar/gkn923, 2009.

[29] A. Subramanian, P. Tamayo, V. K. Mootha, S. Mukherjee, B. L. Ebert, M. A. Gillette, A. Paulovich, S. L. Pomeroy, T. R. Golub, E. S. Lander, J. P. Mesirov, "Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles," Proceedings of the National Academy of Sciences of the United States of America, vol. 102 no. 43, pp. 15545-15550, DOI: 10.1073/pnas.0506580102, 2005.

[30] C. Y. Wang, C. C. Chiao, N. N. Phan, C. Y. Li, Z. D. Sun, J. Z. Jiang, J. H. Hung, Y. L. Chen, M. C. Yen, T. Y. Weng, W. C. Chen, H. P. Hsu, M. D. Lai, "Gene signatures and potential therapeutic targets of amino acid metabolism in estrogen receptor-positive breast cancer," American Journal of Cancer Research, vol. 10 no. 1, pp. 95-113, 2020.

[31] T. J. Kao, C. C. Wu, N. N. Phan, Y. H. Liu, H. D. K. Ta, G. Anuraga, Y. F. Wu, K. H. Lee, J. Y. Chuang, C. Y. Wang, "Prognoses and genomic analyses of proteasome 26S subunit, ATPase (PSMC) family genes in clinical breast cancer," Aging (Albany NY), vol. 13 no. 14,DOI: 10.18632/aging.203345, 2021.

[32] C. Y. Wang, Y. J. Chao, Y. L. Chen, T. W. Wang, N. N. Phan, H. P. Hsu, Y. S. Shan, M. D. Lai, "Upregulation of peroxisome proliferator-activated receptor- α and the lipid metabolism pathway promotes carcinogenesis of ampullary cancer," International Journal of Medical Sciences, vol. 18 no. 1, pp. 256-269, DOI: 10.7150/ijms.48123, 2021.

[33] H. L. Liu, I. J. Yeh, N. N. Phan, Y. H. Wu, M. C. Yen, J. H. Hung, C. C. Chiao, C. F. Chen, Z. Sun, J. Z. Jiang, H. P. Hsu, C. Y. Wang, M. D. Lai, "Gene signatures of SARS-CoV/SARS-CoV-2-infected ferret lungs in short- and long-term models," Infection, Genetics and Evolution, vol. 85, article 104438,DOI: 10.1016/j.meegid.2020.104438, 2020.

[34] Y. H. Wu, I. J. Yeh, N. N. Phan, M. C. Yen, J. H. Hung, C. C. Chiao, C. F. Chen, Z. Sun, H. P. Hsu, C. Y. Wang, M. D. Lai, "Gene signatures and potential therapeutic targets of Middle East respiratory syndrome coronavirus (MERS-CoV)-infected human lung adenocarcinoma epithelial cells," Journal of Microbiology, Immunology, and Infection, vol. 54 no. 5, pp. 845-857, DOI: 10.1016/j.jmii.2021.03.007, 2021.

[35] M. Ashburner, C. A. Ball, J. A. Blake, D. Botstein, H. Butler, J. M. Cherry, A. P. Davis, K. Dolinski, S. S. Dwight, J. T. Eppig, M. A. Harris, D. P. Hill, L. Issel-Tarver, A. Kasarskis, S. Lewis, J. C. Matese, J. E. Richardson, M. Ringwald, G. M. Rubin, G. Sherlock, "Gene Ontology: tool for the unification of biology," Nature Genetics, vol. 25 no. 1, pp. 25-29, DOI: 10.1038/75556, 2000.

[36] D. Szklarczyk, A. L. Gable, K. C. Nastou, D. Lyon, R. Kirsch, S. Pyysalo, N. T. Doncheva, M. Legeay, T. Fang, P. Bork, L. J. Jensen, C. von Mering, "The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets," Nucleic Acids Research, vol. 49 no. D1, pp. D605-d612, DOI: 10.1093/nar/gkaa1074, 2021.

[37] Y. Wang, Y. Luo, Y. Yao, Y. Ji, L. Feng, F. Du, X. Zheng, T. Tao, X. Zhai, Y. Li, P. Han, B. Xu, H. Zhao, "Silencing the lncRNA Maclpil in pro-inflammatory macrophages attenuates acute experimental ischemic stroke via LCP1 in mice," Journal of Cerebral Blood Flow and Metabolism, vol. 40 no. 4, pp. 747-759, DOI: 10.1177/0271678X19836118, 2020.

[38] Y. Wang, C. Liu, Y. Chen, T. Chen, T. Han, L. Xue, B. Xu, "Systemically silencing long non-coding RNAs Maclpil with short interfering RNA nanoparticles alleviates experimental ischemic stroke by promoting macrophage apoptosis and anti-inflammatory activation," Front Cardiovasc Med, vol. 9, article 876087,DOI: 10.3389/fcvm.2022.876087, 2022.

[39] R. Rajeev, R. Prathiviraj, G. S. Kiran, J. Selvin, "Zoonotic evolution and implications of microbiome in viral transmission and infection," Virus Research, vol. 290, article 198175,DOI: 10.1016/j.virusres.2020.198175, 2020.

[40] C. M. Chao, C. C. Lai, W. L. Yu, "COVID-19 associated mucormycosis - an emerging threat," Journal of Microbiology, Immunology, and Infection, vol. 55 no. 2, pp. 183-190, DOI: 10.1016/j.jmii.2021.12.007, 2022.

[41] C. K. Fan, C. W. Liao, Y. C. Cheng, "Factors affecting disease manifestation of toxocarosis in humans: genetics and environment," Veterinary Parasitology, vol. 193 no. 4, pp. 342-352, DOI: 10.1016/j.vetpar.2012.12.030, 2013.

[42] G. H. Bai, S. C. Lin, Y. H. Hsu, S. Y. Chen, "The human virome: viral metagenomics, relations with human diseases, and therapeutic applications," Viruses, vol. 14 no. 2,DOI: 10.3390/v14020278, 2022.

[43] G. Y. Lan, Y. J. Lee, J. C. Wu, H. Y. Lai, H. Y. Liu, H. C. Chuang, K. Li-Chun Hsieh, "Serial quantitative chest computed tomography imaging as prognosticators of coronavirus disease 2019 pneumonia," Journal of the Formosan Medical Association, vol. 121 no. 3, pp. 718-722, DOI: 10.1016/j.jfma.2021.07.021, 2022.

[44] P. F. Mwale, C. H. Lee, L. T. Lin, S. J. Leu, Y. J. Huang, L. C. Chiang, Y. C. Mao, Y. Y. Yang, "Expression, purification, and characterization of anti-zika virus envelope protein: polyclonal and chicken-derived single chain variable fragment antibodies," International Journal of Molecular Sciences, vol. 21 no. 2,DOI: 10.3390/ijms21020492, 2020.

[45] Z. W. Ye, S. Yuan, K. S. Yuen, S. Y. Fung, C. P. Chan, D. Y. Jin, "Zoonotic origins of human coronaviruses," International Journal of Biological Sciences, vol. 16 no. 10, pp. 1686-1697, DOI: 10.7150/ijbs.45472, 2020.

[46] J. Zahradník, S. Marciano, M. Shemesh, E. Zoler, D. Harari, J. Chiaravalli, B. Meyer, Y. Rudich, C. Li, I. Marton, O. Dym, N. Elad, M. G. Lewis, H. Andersen, M. Gagne, R. A. Seder, D. C. Douek, G. Schreiber, "SARS-CoV-2 variant prediction and antiviral drug design are enabled by RBD in vitro evolution," Nature Microbiology, vol. 6 no. 9, pp. 1188-1198, DOI: 10.1038/s41564-021-00954-4, 2021.

[47] B. Xu, G. Li, J. Guo, T. Ikezoe, K. Kasirajan, S. Zhao, R. L. Dalman, "Angiotensin-converting enzyme 2, coronavirus disease 2019, and abdominal aortic aneurysms," Journal of Vascular Surgery, vol. 74 no. 5, pp. 1740-1751, DOI: 10.1016/j.jvs.2021.01.051, 2021.

[48] C. H. Liu, C. H. Lu, L. T. Lin, "Pandemic strategies with computational and structural biology against COVID-19: a retrospective," Computational and Structural Biotechnology Journal, vol. 20, pp. 187-192, DOI: 10.1016/j.csbj.2021.11.040, 2022.

[49] H. G. Lim, S. H. Hsiao, Y. C. Fann, Y. G. Lee, "Robust mutation profiling of SARS-CoV-2 variants from multiple raw illumina sequencing data with cloud workflow," Genes (Basel), vol. 13 no. 4,DOI: 10.3390/genes13040686, 2022.

[50] W. L. Yu, H. S. Toh, C. T. Liao, W. T. Chang, "A double-edged sword-cardiovascular concerns of potential anti-COVID-19 drugs," Cardiovascular Drugs and Therapy, vol. 35 no. 2, pp. 205-214, DOI: 10.1007/s10557-020-07024-7, 2021.

[51] T. N. Poly, M. M. Islam, Y. J. Li, B. Alsinglawi, M. H. Hsu, W. S. Jian, H. C. Yang, "Application of artificial intelligence for screening COVID-19 patients using digital images: meta-analysis," JMIR Medical Informatics, vol. 9 no. 4, article e21394,DOI: 10.2196/21394, 2021.

[52] C. S. Yu, S. S. Chang, T. H. Chang, J. L. Wu, Y. J. Lin, H. F. Chien, R. J. Chen, "A COVID-19 pandemic artificial intelligence-based system with deep learning forecasting and automatic statistical data acquisition: development and implementation study," Journal of Medical Internet Research, vol. 23 no. 5, article e27806,DOI: 10.2196/27806, 2021.

[53] H. Zhang, T. Chen, Y. Zhang, J. Lin, W. Zhao, Y. Shi, H. Lau, Y. Zhang, M. Yang, C. Xu, L. Tang, B. Xu, J. Jiang, X. Chen, "Crucial genes in aortic dissection identified by weighted gene coexpression network analysis," Journal of Immunology Research, vol. 2022,DOI: 10.1155/2022/7585149, 2022.

[54] Y. Li, Y. Wang, Y. Yao, B. B. Griffiths, L. Feng, T. Tao, F. Wang, B. Xu, C. M. Stary, H. Zhao, "Systematic study of the immune components after ischemic stroke using CyTOF techniques," Journal of Immunology Research, vol. 2020,DOI: 10.1155/2020/9132410, 2020.

[55] J. Ruby, H. Bluethmann, J. J. Peschon, "Antiviral activity of tumor necrosis factor (TNF) is mediated via p55 and p75 TNF receptors," The Journal of Experimental Medicine, vol. 186 no. 9, pp. 1591-1596, DOI: 10.1084/jem.186.9.1591, 1997.

[56] M. H. Orzalli, A. Smith, K. A. Jurado, A. Iwasaki, J. A. Garlick, J. C. Kagan, "An antiviral branch of the IL-1 signaling pathway restricts immune-evasive virus replication," Molecular Cell, vol. 71 no. 5, pp. 825-840.e6, DOI: 10.1016/j.molcel.2018.07.009, 2018.

[57] G. Wang, L. Jiang, J. Wang, J. Zhang, F. Kong, Q. Li, Y. Yan, S. Huang, Y. Zhao, L. Liang, J. Li, N. Sun, Y. Hu, W. Shi, G. Deng, P. Chen, L. Liu, X. Zeng, G. Tian, Z. Bu, H. Chen, C. Li, "The G protein-coupled receptor FFAR2 promotes internalization during influenza a virus entry," Journal of Virology, vol. 94 no. 2,DOI: 10.1128/JVI.01707-19, 2020.

[58] N. Tsutsumi, Q. Qu, M. Mavri, M. S. Baggesen, S. Maeda, D. Waghray, C. Berg, B. K. Kobilka, M. M. Rosenkilde, G. Skiniotis, K. C. Garcia, "Structural basis for the constitutive activity and immunomodulatory properties of the Epstein-Barr virus-encoded G protein-coupled receptor BILF1," Immunity, vol. 54 no. 7, pp. 1405-1416.e7, DOI: 10.1016/j.immuni.2021.06.001, 2021.

[59] H. Fausther-Bovendo, X. Qiu, S. McCorrister, G. Westmacott, P. Sandstrom, C. Castilletti, A. Di Caro, G. Ippolito, G. P. Kobinger, "Ebola virus infection induces autoimmunity against dsDNA and HSP60," Scientific Reports, vol. 7 no. 1, pp. 42147-42147, DOI: 10.1038/srep42147, 2017.

[60] S. G. Park, G. Jung, "Human hepatitis B virus polymerase interacts with the molecular chaperonin Hsp60," Journal of Virology, vol. 75 no. 15, pp. 6962-6968, DOI: 10.1128/JVI.75.15.6962-6968.2001, 2001.

[61] S. Pujhari, M. Brustolin, V. M. Macias, R. H. Nissly, M. Nomura, S. V. Kuchipudi, J. L. Rasgon, "Heat shock protein 70 (Hsp70) mediates Zika virus entry, replication, and egress from host cells," Emerg Microbes Infect, vol. 8 no. 1,DOI: 10.1080/22221751.2018.1557988, 2019.

[62] Y. Jang, M. Jin, S. H. Seo, "Histamine contributes to severe pneumonia in pigs infected with 2009 pandemic H1N1 influenza virus," Archives of Virology, vol. 163 no. 11, pp. 3015-3022, DOI: 10.1007/s00705-018-3973-9, 2018.

[63] T. Silva, C. Jeewandara, L. Gomes, C. Gangani, S. D. Mahapatuna, T. Pathmanathan, A. Wijewickrama, G. S. Ogg, G. N. Malavige, "Urinary leukotrienes and histamine in patients with varying severity of acute dengue," PLoS One, vol. 16 no. 2, article e0245926,DOI: 10.1371/journal.pone.0245926, 2021.

[64] J. A. Tree, G. Hall, G. Pearson, E. Rayner, V. A. Graham, K. Steeds, K. R. Bewley, G. J. Hatch, M. Dennis, I. Taylor, A. D. Roberts, S. G. Funnell, J. Vipond, "Sequence of pathogenic events in cynomolgus macaques infected with aerosolized monkeypox virus," Journal of Virology, vol. 89 no. 8, pp. 4335-4344, DOI: 10.1128/JVI.03029-14, 2015.

[65] J. Maelfait, K. Roose, P. Bogaert, M. Sze, X. Saelens, M. Pasparakis, I. Carpentier, G. van Loo, R. Beyaert, "A20 (Tnfaip3) deficiency in myeloid cells protects against influenza A virus infection," PLoS Pathogens, vol. 8 no. 3, article e1002570,DOI: 10.1371/journal.ppat.1002570, 2012.

[66] S. S. D. de Azevedo, M. Ribeiro-Alves, F. H. Côrtes, E. Delatorre, L. Spangenberg, H. Naya, L. N. Seito, B. Hoagland, B. Grinsztejn, V. G. Veloso, M. G. Morgado, T. M. L. Souza, G. Bello, "Increased expression of CDKN1A/p21 in HIV-1 controllers is correlated with upregulation of ZC3H12A/MCPIP1," Retrovirology, vol. 17 no. 1,DOI: 10.1186/s12977-020-00522-4, 2020.

[67] D. S. Chandrashekar, M. Athar, U. Manne, S. Varambally, "Comparative transcriptome analyses reveal genes associated with SARS-CoV-2 infection of human lung epithelial cells," Scientific Reports, vol. 11 no. 1,DOI: 10.1038/s41598-021-95733-w, 2021.

[68] L. C. Cheng, T. J. Kao, N. N. Phan, C. C. Chiao, M. C. Yen, C. F. Chen, J. H. Hung, J. Z. Jiang, Z. Sun, C. Y. Wang, H. P. Hsu, "Novel signaling pathways regulate SARS-CoV and SARS-CoV-2 infectious disease," Medicine (Baltimore), vol. 100 no. 7, article e24321,DOI: 10.1097/MD.0000000000024321, 2021.

[69] Y. H. Wu, I. J. Yeh, N. N. Phan, M. C. Yen, H. L. Liu, C. Y. Wang, H. P. Hsu, "Severe acute respiratory syndrome coronavirus (SARS-CoV)-2 infection induces dysregulation of immunity: in silico gene expression analysis," International Journal of Medical Sciences, vol. 18 no. 5, pp. 1143-1152, DOI: 10.7150/ijms.52256, 2021.

[70] E. Brunetta, M. Folci, B. Bottazzi, M. De Santis, G. Gritti, A. Protti, S. N. Mapelli, S. Bonovas, D. Piovani, R. Leone, I. My, V. Zanon, G. Spata, M. Bacci, D. Supino, S. Carnevale, M. Sironi, S. Davoudian, C. Peano, F. Landi, F. Di Marco, F. Raimondi, A. Gianatti, C. Angelini, A. Rambaldi, C. Garlanda, M. Ciccarelli, M. Cecconi, A. Mantovani, "Macrophage expression and prognostic significance of the long pentraxin PTX3 in COVID-19," Nature Immunology, vol. 22 no. 1, pp. 19-24, DOI: 10.1038/s41590-020-00832-x, 2021.

[71] B. Dahal, S. C. Lin, B. D. Carey, J. L. Jacobs, J. D. Dinman, M. L. van Hoek, A. A. Adams, K. Kehn-Hall, "EGR1 upregulation following Venezuelan equine encephalitis virus infection is regulated by ERK and PERK pathways contributing to cell death," Virology, vol. 539, pp. 121-128, DOI: 10.1016/j.virol.2019.10.016, 2020.

[72] R. Vlahos, S. Selemidis, "NADPH oxidases as novel pharmacologic targets against influenza A virus infection," Molecular Pharmacology, vol. 86 no. 6, pp. 747-759, DOI: 10.1124/mol.114.095216, 2014.

[73] Z. Saud, M. D. Hitchings, T. M. Butt, "Nanopore sequencing and de novo assembly of a misidentified camelpox vaccine reveals putative epigenetic modifications and alternate protein signal peptides," Scientific Reports, vol. 11 no. 1,DOI: 10.1038/s41598-021-97158-x, 2021.

[74] D. M. Knipe, "Nuclear sensing of viral DNA, epigenetic regulation of herpes simplex virus infection, and innate immunity," Virology, vol. 479-480, pp. 153-159, DOI: 10.1016/j.virol.2015.02.009, 2015.

[75] E. L. McGrath, S. L. Rossi, J. Gao, S. G. Widen, A. C. Grant, T. J. Dunn, S. R. Azar, C. M. Roundy, Y. Xiong, D. J. Prusak, B. D. Loucas, T. G. Wood, Y. Yu, I. Fernández-Salas, S. C. Weaver, N. Vasilakis, P. Wu, "Differential responses of human fetal brain neural stem cells to Zika virus infection," Stem Cell Reports, vol. 8 no. 3, pp. 715-727, DOI: 10.1016/j.stemcr.2017.01.008, 2017.

[76] R. Brahma, S. Gurumayum, L. D. Naorem, M. Muthaiyan, J. Gopal, A. Venkatesan, "Identification of hub genes and pathways in Zika virus infection using RNA-Seq data: a network-based computational approach," Viral Immunology, vol. 31 no. 4, pp. 321-332, DOI: 10.1089/vim.2017.0116, 2018.

[77] K. H. Rubins, L. E. Hensley, D. A. Relman, P. O. Brown, "Stunned silence: gene expression programs in human cells infected with monkeypox or vaccinia virus," PLoS One, vol. 6 no. 1, article e15615,DOI: 10.1371/journal.pone.0015615, 2011.

[78] Z. Jezek, M. Szczeniowski, K. M. Paluku, M. Mutombo, B. Grab, "Human monkeypox: confusion with chickenpox," Acta Tropica, vol. 45 no. 4, pp. 297-307, 1988.