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

Xuan, Do Thi Minh; Yeh, I-Jeng; Wu, Chung-Che; Su, Che-Yu; Liu, Hsin-Liang; Chiao, Chung-Chieh; Ku, Su-Chi; Jiang, Jia-Zhen; Sun, Zhengda; Ta, Hoang Dang Khoa; Anuraga, Gangga; Wang, Chih-Yang; Yen, Meng-Chi

doi:https://doi.org/10.1155/2022/3883822

Journal of Immunology Research

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 3883822 | https://doi.org/10.1155/2022/3883822

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

Do Thi Minh Xuan,¹I-Jeng Yeh,^2,3Chung-Che Wu,^4,5Che-Yu Su,²Hsin-Liang Liu,²Chung-Chieh Chiao,⁶Su-Chi Ku,^1,7Jia-Zhen Jiang,⁸Zhengda Sun,⁹Hoang Dang Khoa Ta,⁶Gangga Anuraga,^6,10and Chih-Yang Wang^1,6,11 et al.

Academic Editor: Baohui Xu

Received29 Jun 2022

Accepted21 Jul 2022

Published01 Sept 2022

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.

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 1

Schematic illustration of the study design. High-throughput data of monkeypox-infected 7-hours of Macaca mulatta kidney epithelial (MK2) cells and monkeypox-infected 6-hours of human immortal epithelial cancer (HeLa) cells were both acquired from the GEO database. The top 5% of differentially expressed genes in the two monkeypox-infected models were determined through a Venn diagram analysis. Results of pathway analyses and functional interpretations were analyzed using DAVID, GO, STRING, and MetaCore.

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 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 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 -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.

(a)

(b)

Figure 2

Analytical results of Gene Ontology (GO) enrichment accompanied by heatmap visualization in two separate monkeypox-infected models. (a) Comparison of gene expression patterns between the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cells and the mock control. Cluster 1 (c1) contained the most significant highly expressed genes and pathways in the monkeypox-infected group, whereas c2 exerted an opposite trend. The GO analysis suggested that inflammatory and chemokine-related peptidoglycan pathways were significantly correlated with the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model. (b) Comparison of gene expression patterns between the monkeypox-infected human immortal epithelial cancer (HeLa) cell model and the mock control. Cluster 1 (c1) contained the most significantly highly expressed genes and pathways in the monkeypox-infected group, whereas c2 exerted an opposite trend. The GO analysis suggested that nucleosome assembly and cytoskeleton-related pathways were found to be significantly correlated with the human immortal epithelial cancer (HeLa) cell model.

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).

(a)

(b)

(c)

(d)

(e)

Figure 3

Overexpressed pathways shared by the two in vitro monkeypox-infected models. (a) The top 5% of overexpressed genes shared by the GSE21001 and GSE36854 datasets, filtered by a Venn diagram, represent the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model and the human immortal epithelial cancer (HeLa) cell model. (b) The most enriched biological pathways corresponding to genes from the shortlist of upregulated genes shared by the two models, in order of decreasing log ( values). (c) The most enriched metabolic pathways as components of the Endogenous Metabolic Networks corresponding to genes from the shortlist of upregulated genes shared by the two models, in order of decreasing log ( values). (d) Related pathways and network analyses by MetaCore confirmed the vital role of the “immune response_IL-1 signaling” pathway in both monkeypox-infected models. (e) Endogenous Metabolic Networks analyses by MetaCore confirmed the vital role of the “N-acyl-sphingosine phosphate pathway” in both monkeypox-infected models.

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 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)).

(a)

(b)

(c)

Figure 4

Overview of protein-protein-interacting networks inferred from overexpressed genes in the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model and the human immortal epithelial cancer (HeLa) cell model. (a) Similar transcriptome patterns obtained from the GSE21001 and GSE36854 datasets were filtered by a Venn diagram. Shortlist of the top 1% of upregulated genes shared by the two monkeypox-infected models. (b) Distribution of 24 coregulated genes inferred from the above shortlist, presented using box plot visualizations. (c) STRING analyses revealed the protein-protein-interacting network of the 24 co-regulated genes, separated by clusters in different colors. Colored nodes represent 24 target genes as input, while gray nodes represent the corresponding proteins.

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 5

Major characteristics of Gene Ontology (GO) and cellular component annotations are inferred from highly expressed genes shared by the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model and the human immortal epithelial cancer (HeLa) cell model. (a) Top highly enriched GO cellular components in the MK2 model. (b) Top highly enriched GO cellular components in the HeLa model. (c) Conclusion of enriched maps discovered by bioinformatic analysis of upregulated genes from the MK2 and HeLa models. The data revealed that interleukins (ILs), G protein-coupled receptors (GPCRs), heat shock proteins (HSPs), Toll-like receptors (TLRs), and metabolic-related pathways were the elevated expression of both monkeypox-infected monkey and human cell lines. Besides, cluster of differentiation 40 (CD40), plasmin, and histamine serve as major regulators in the monkeypox-infected monkey cell line model, while interferons (IFNs), macrophages, and neutrophil-related signaling dominate the monkeypox-infected human cell line model.

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.

Data Availability

The datasets used and analyzed during the current study are available from the corresponding authors on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

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.

Acknowledgments

This research was supported by grants from the Ministry of Science and Technology (MOST) of Taiwan (MOST-109-2320-B-038-009-MY2 and MOST-107-2320-B-037-011-MY3), from Taipei Medical University Hospital (109TMU-TMUH-19 and 111TMU-TMUH-08), and from Kaohsiung Medical University Hospital (KMUH110-0T06, KMUH109-9R82 and KMUH110-0M75). This work was financially supported by the Higher Education Sprout Project of the Ministry of Education (MOE) in Taiwan. The authors acknowledge the online platform for data analysis and visualization (http://www.bioinformatics.com.cn/), and we also thank the statistical/computational/technical support of the Clinical Data Center, Office of Data Science, Taipei Medical University, Taiwan.

Supplementary Materials

Figure S1: MetaCore-generated pathways regulated by top upregulated genes in the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cell model. Figure S2: MetaCore-generated Endogenous Metabolic Network among top upregulated genes in the monkeypox-infected Macaca mulatta kidney epithelial (MK2) cells, performed on the GSE21001 dataset. Figure S3: MetaCore-generated pathways related to top upregulated genes in the human immortal epithelial cancer (HeLa) cell model. Figure S4: MetaCore-generated Endogenous Metabolic Network among top upregulated genes in the human immortal epithelial cancer (HeLa) cell model, performed on the GSE36854 dataset. Figure S5: Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathways in Formula A analysis in different models. Figure S6: heatmap displaying specific differentially regulated genes within the monkeypox-infected human immortal epithelial cancer (HeLa) cell model (GSE24125) versus the mock-infected group. Table S1: top 50 upregulated pathways in the MetaCore analysis of potential maps. Comparison between monkeypox-infected Macaca mulatta kidney epithelial (MK2) cells and human epithelial (HeLa) cells from the GSE21001 and GSE36854 datasets. Table S2: top 50 upregulated pathways in the MetaCore analysis of potential Endogenous Metabolic Networks. Comparison between monkeypox-infected Macaca mulatta kidney epithelial (MK2) cells and human epithelial (HeLa) cells from the GSE21001 and GSE36854 datasets. Table S3: top 50 upregulated pathways in the MetaCore analysis of potential maps. Comparison between monkeypox-infected and mock-infectedMacaca mulatta kidney epithelial (MK2) cells from the GSE21001 dataset. Table S4: top 50 upregulated pathways in the MetaCore analysis of potential Endogenous Metabolic Networks. Comparison between monkeypox-infected and mock-infectedMacaca mulatta kidney epithelial (MK2) cells from the GSE21001 dataset. Table S5: top 50 upregulated pathways in the MetaCore analysis of potential maps. Comparison between monkeypox-infected and mock-infected human epithelial (HeLa) cells from the GSE36854 dataset. Table S6: top 50 upregulated pathways in the MetaCore analysis of potential Endogenous Metabolic Networks. Comparison between monkeypox-infected and mock-infected human epithelial (HeLa) cells from the GSE36854 dataset. (Supplementary Materials)

References

E. Tambo and A. M. Al-Nazawi, “Combating the global spread of poverty-related monkeypox outbreaks and beyond,” Infectious Diseases of Poverty, vol. 11, no. 1, p. 80, 2022.
View at: Publisher Site | Google Scholar
S. Essbauer, M. Pfeffer, and H. Meyer, “Zoonotic poxviruses,” Veterinary Microbiology, vol. 140, no. 3-4, pp. 229–236, 2010.
View at: Publisher Site | Google Scholar
S. N. Shchelkunov, A. V. Totmenin, I. V. Babkin et al., “Human monkeypox and smallpox viruses: genomic comparison,” Human monkeypox and smallpox viruses: genomic comparison, vol. 509, no. 1, pp. 66–70, 2001.
View at: Publisher Site | Google Scholar
S. Parker, D. A. Schultz, H. Meyer, and R. M. Buller, “Smallpox and Monkeypox Viruses☆,” in Reference Module in Biomedical Sciences, Elsevier, 2014.
View at: Google Scholar
D. B. Di Giulio and P. B. Eckburg, “Human monkeypox: an emerging zoonosis,” The Lancet Infectious Diseases, vol. 4, no. 1, pp. 15–25, 2004.
View at: Publisher Site | Google Scholar
M. Kozlov, “Monkeypox goes global: why scientists are on alert,” Nature, vol. 606, no. 7912, pp. 15-16, 2022.
View at: Publisher Site | Google Scholar
J. G. Breman, R. Kalisa, M. V. Steniowski, E. Zanotto, A. I. Gromyko, and I. Arita, “Human monkeypox, 1970-79,” Bulletin of the World Health Organization, vol. 58, no. 2, pp. 165–182, 1980.
View at: Google Scholar
H. Meyer, M. Perrichot, M. Stemmler et al., “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, 2002.
View at: Publisher Site | Google Scholar
P. Y. Nguyen, W. S. Ajisegiri, V. Costantino, A. A. Chughtai, and 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, 2021.
View at: Publisher Site | Google Scholar
R. Vivancos, C. Anderson, P. Blomquist et al., “Community transmission of monkeypox in the United Kingdom, April to May 2022,” Eurosurveillance, vol. 27, no. 22, 2022.
View at: Publisher Site | Google Scholar
B. Titanji, B. Tegomoh, S. Nematollahi, M. Konomos, and P. A. Kulkarni, “Monkeypox-a contemporary review for healthcare professionals,” in Open Forum Infectious Diseases, vol. 9, no. 7, 2022.
View at: Publisher Site | Google Scholar
N. Kumar, A. Acharya, H. E. Gendelman, and S. N. Byrareddy, “The 2022 outbreak and the pathobiology of the monkeypox virus,” Journal of Autoimmunity, vol. 131, p. 102855, 2022.
View at: Publisher Site | Google Scholar
J. Heskin, A. Belfield, C. Milne et al., “Transmission of monkeypox virus through sexual contact - a novel route of infection,” Journal of Infection, vol. 85, no. 3, pp. 334–363, 2022.
View at: Publisher Site | Google Scholar
E. Petersen, A. Zumla, D. Hui et al., “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, 2022.
View at: Publisher Site | Google Scholar
W H Organization, “Surveillance, case investigation and contact tracing for monkeypox: interim guidance, 22 May 2022,” World Health Organization, 2022.
View at: Google Scholar
P. Dashraath, K. Nielsen-Saines, C. Mattar, D. Musso, P. Tambyah, and D. Baud, “Guidelines for pregnant individuals with monkeypox virus exposure,” The Lancet, vol. 400, no. 10345, pp. 21-22, 2022.
View at: Publisher Site | Google Scholar
J. R. Kugelman, S. C. Johnston, P. M. Mulembakani et al., “Genomic variability of monkeypox virus among humans, Democratic Republic of the Congo,” Emerging Infectious Diseases, vol. 20, no. 2, pp. 232–239, 2014.
View at: Publisher Site | Google Scholar
U. A. Awan, S. Riasat, W. Naeem, S. Kamran, A. A. Khattak, and S. Khan, “Monkeypox: a new threat at our doorstep!,” The Journal of Infection, vol. 85, no. 2, pp. e47–e48, 2022.
View at: Publisher Site | Google Scholar
A. Alkhalil, R. Hammamieh, J. Hardick, M. A. Ichou, M. Jett, and S. Ibrahim, “Gene expression profiling of monkeypox virus-infected cells reveals novel interfaces for host-virus interactions,” Virology Journal, vol. 7, no. 1, 2010.
View at: Publisher Site | Google Scholar
D. Bourquain, P. W. Dabrowski, and 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, 2013.
View at: Publisher Site | Google Scholar
C.-Y. Wang, Y.-C. Chang, Y.-L. Kuo et al., “Mutation of the PTCH1 gene predicts recurrence of breast cancer,” Scientific reports, vol. 9, no. 1, pp. 1–14, 2019.
View at: Google Scholar
D. L. Cooke, D. B. McCoy, V. V. Halbach et al., “Endovascular biopsy: in vivo cerebral aneurysm endothelial cell sampling and gene expression analysis,” Translational Stroke Research, vol. 9, no. 1, pp. 20–33, 2018.
View at: Publisher Site | Google Scholar
T. Y. Weng, H. F. Wu, C. Y. Li et al., “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, p. 8216, 2018.
View at: Publisher Site | Google Scholar
P. Shahi, C. Y. Wang, D. A. Lawson et al., “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, 2017.
View at: Publisher Site | Google Scholar
N. N. Phan, C. Y. Wang, and 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, 2014.
View at: Publisher Site | Google Scholar
Y. H. Hung, Y. S. Chan, Y. S. Chang et al., “Fatty acid metabolic enzyme acyl-CoA thioesterase 8 promotes the development of hepatocellular carcinoma,” Oncology Reports, vol. 31, no. 6, pp. 2797–2803, 2014.
View at: Publisher Site | Google Scholar
W. H. Da, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
W. H. Da, B. T. Sherman, and R. A. Lempicki, “Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists,” Nucleic Acids Research, vol. 37, no. 1, pp. 1–13, 2009.
View at: Publisher Site | Google Scholar
A. Subramanian, P. Tamayo, V. K. Mootha et al., “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, 2005.
View at: Publisher Site | Google Scholar
C. Y. Wang, C. C. Chiao, N. N. Phan et al., “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.
View at: Google Scholar
T. J. Kao, C. C. Wu, N. N. Phan et al., “Prognoses and genomic analyses of proteasome 26S subunit, ATPase (PSMC) family genes in clinical breast cancer,” Aging (Albany NY), vol. 13, no. 14, p. 17970, 2021.
View at: Publisher Site | Google Scholar
C. Y. Wang, Y. J. Chao, Y. L. Chen et al., “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, 2021.
View at: Publisher Site | Google Scholar
H. L. Liu, I. J. Yeh, N. N. Phan et al., “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, 2020.
View at: Publisher Site | Google Scholar
Y. H. Wu, I. J. Yeh, N. N. Phan et al., “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, 2021.
View at: Publisher Site | Google Scholar
M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene Ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, K. C. Nastou et al., “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, 2021.
View at: Publisher Site | Google Scholar
Y. Wang, Y. Luo, Y. Yao et al., “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, 2020.
View at: Publisher Site | Google Scholar
Y. Wang, C. Liu, Y. Chen et al., “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, 2022.
View at: Publisher Site | Google Scholar
R. Rajeev, R. Prathiviraj, G. S. Kiran, and J. Selvin, “Zoonotic evolution and implications of microbiome in viral transmission and infection,” Virus Research, vol. 290, article 198175, 2020.
View at: Publisher Site | Google Scholar
C. M. Chao, C. C. Lai, and W. L. Yu, “COVID-19 associated mucormycosis - an emerging threat,” Journal of Microbiology, Immunology, and Infection, vol. 55, no. 2, pp. 183–190, 2022.
View at: Publisher Site | Google Scholar
C. K. Fan, C. W. Liao, and Y. C. Cheng, “Factors affecting disease manifestation of toxocarosis in humans: genetics and environment,” Veterinary Parasitology, vol. 193, no. 4, pp. 342–352, 2013.
View at: Publisher Site | Google Scholar
G. H. Bai, S. C. Lin, Y. H. Hsu, and S. Y. Chen, “The human virome: viral metagenomics, relations with human diseases, and therapeutic applications,” Viruses, vol. 14, no. 2, p. 278, 2022.
View at: Publisher Site | Google Scholar
G. Y. Lan, Y. J. Lee, J. C. Wu et al., “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, 2022.
View at: Publisher Site | Google Scholar
P. F. Mwale, C. H. Lee, L. T. Lin et al., “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, p. 492, 2020.
View at: Publisher Site | Google Scholar
Z. W. Ye, S. Yuan, K. S. Yuen, S. Y. Fung, C. P. Chan, and D. Y. Jin, “Zoonotic origins of human coronaviruses,” International Journal of Biological Sciences, vol. 16, no. 10, pp. 1686–1697, 2020.
View at: Publisher Site | Google Scholar
J. Zahradník, S. Marciano, M. Shemesh et al., “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, 2021.
View at: Publisher Site | Google Scholar
B. Xu, G. Li, J. Guo et al., “Angiotensin-converting enzyme 2, coronavirus disease 2019, and abdominal aortic aneurysms,” Journal of Vascular Surgery, vol. 74, no. 5, pp. 1740–1751, 2021.
View at: Publisher Site | Google Scholar
C. H. Liu, C. H. Lu, and 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, 2022.
View at: Publisher Site | Google Scholar
H. G. Lim, S. H. Hsiao, Y. C. Fann, and 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, p. 686, 2022.
View at: Publisher Site | Google Scholar
W. L. Yu, H. S. Toh, C. T. Liao, and 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, 2021.
View at: Publisher Site | Google Scholar
T. N. Poly, M. M. Islam, Y. J. Li et al., “Application of artificial intelligence for screening COVID-19 patients using digital images: meta-analysis,” JMIR Medical Informatics, vol. 9, no. 4, article e21394, 2021.
View at: Publisher Site | Google Scholar
C. S. Yu, S. S. Chang, T. H. Chang et al., “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, 2021.
View at: Publisher Site | Google Scholar
H. Zhang, T. Chen, Y. Zhang et al., “Crucial genes in aortic dissection identified by weighted gene coexpression network analysis,” Journal of Immunology Research, vol. 2022, Article ID 7585149, 18 pages, 2022.
View at: Publisher Site | Google Scholar
Y. Li, Y. Wang, Y. Yao et al., “Systematic study of the immune components after ischemic stroke using CyTOF techniques,” Journal of Immunology Research, vol. 2020, Article ID 9132410, 2020.
View at: Publisher Site | Google Scholar
J. Ruby, H. Bluethmann, and 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, 1997.
View at: Publisher Site | Google Scholar
M. H. Orzalli, A. Smith, K. A. Jurado, A. Iwasaki, J. A. Garlick, and 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, 2018.
View at: Publisher Site | Google Scholar
G. Wang, L. Jiang, J. Wang et al., “The G protein-coupled receptor FFAR2 promotes internalization during influenza a virus entry,” Journal of Virology, vol. 94, no. 2, 2020.
View at: Publisher Site | Google Scholar
N. Tsutsumi, Q. Qu, M. Mavri et al., “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, 2021.
View at: Publisher Site | Google Scholar
H. Fausther-Bovendo, X. Qiu, S. McCorrister et al., “Ebola virus infection induces autoimmunity against dsDNA and HSP60,” Scientific Reports, vol. 7, no. 1, pp. 42147–42147, 2017.
View at: Publisher Site | Google Scholar
S. G. Park and G. Jung, “Human hepatitis B virus polymerase interacts with the molecular chaperonin Hsp60,” Journal of Virology, vol. 75, no. 15, pp. 6962–6968, 2001.
View at: Publisher Site | Google Scholar
S. Pujhari, M. Brustolin, V. M. Macias et al., “Heat shock protein 70 (Hsp70) mediates Zika virus entry, replication, and egress from host cells,” Emerg Microbes Infect, vol. 8, no. 1, pp. 8–16, 2019.
View at: Publisher Site | Google Scholar
Y. Jang, M. Jin, and 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, 2018.
View at: Publisher Site | Google Scholar
T. Silva, C. Jeewandara, L. Gomes et al., “Urinary leukotrienes and histamine in patients with varying severity of acute dengue,” PLoS One, vol. 16, no. 2, article e0245926, 2021.
View at: Publisher Site | Google Scholar
J. A. Tree, G. Hall, G. Pearson et al., “Sequence of pathogenic events in cynomolgus macaques infected with aerosolized monkeypox virus,” Journal of Virology, vol. 89, no. 8, pp. 4335–4344, 2015.
View at: Publisher Site | Google Scholar
J. Maelfait, K. Roose, P. Bogaert et al., “A20 (Tnfaip3) deficiency in myeloid cells protects against influenza A virus infection,” PLoS Pathogens, vol. 8, no. 3, article e1002570, 2012.
View at: Publisher Site | Google Scholar
S. S. D. de Azevedo, M. Ribeiro-Alves, F. H. Côrtes et al., “Increased expression of CDKN1A/p21 in HIV-1 controllers is correlated with upregulation of ZC3H12A/MCPIP1,” Retrovirology, vol. 17, no. 1, p. 18, 2020.
View at: Publisher Site | Google Scholar
D. S. Chandrashekar, M. Athar, U. Manne, and S. Varambally, “Comparative transcriptome analyses reveal genes associated with SARS-CoV-2 infection of human lung epithelial cells,” Scientific Reports, vol. 11, no. 1, p. 16212, 2021.
View at: Publisher Site | Google Scholar
L. C. Cheng, T. J. Kao, N. N. Phan et al., “Novel signaling pathways regulate SARS-CoV and SARS-CoV-2 infectious disease,” Medicine (Baltimore), vol. 100, no. 7, article e24321, 2021.
View at: Publisher Site | Google Scholar
Y. H. Wu, I. J. Yeh, N. N. Phan et al., “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, 2021.
View at: Publisher Site | Google Scholar
E. Brunetta, M. Folci, B. Bottazzi et al., “Macrophage expression and prognostic significance of the long pentraxin PTX3 in COVID-19,” Nature Immunology, vol. 22, no. 1, pp. 19–24, 2021.
View at: Publisher Site | Google Scholar
B. Dahal, S. C. Lin, B. D. Carey et al., “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, 2020.
View at: Publisher Site | Google Scholar
R. Vlahos and S. Selemidis, “NADPH oxidases as novel pharmacologic targets against influenza A virus infection,” Molecular Pharmacology, vol. 86, no. 6, pp. 747–759, 2014.
View at: Publisher Site | Google Scholar
Z. Saud, M. D. Hitchings, and 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, p. 17758, 2021.
View at: Publisher Site | Google Scholar
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, 2015.
View at: Publisher Site | Google Scholar
E. L. McGrath, S. L. Rossi, J. Gao et al., “Differential responses of human fetal brain neural stem cells to Zika virus infection,” Stem Cell Reports, vol. 8, no. 3, pp. 715–727, 2017.
View at: Publisher Site | Google Scholar
R. Brahma, S. Gurumayum, L. D. Naorem, M. Muthaiyan, J. Gopal, and 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, 2018.
View at: Publisher Site | Google Scholar
K. H. Rubins, L. E. Hensley, D. A. Relman, and 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, 2011.
View at: Publisher Site | Google Scholar
Z. Jezek, M. Szczeniowski, K. M. Paluku, M. Mutombo, and B. Grab, “Human monkeypox: confusion with chickenpox,” Acta Tropica, vol. 45, no. 4, pp. 297–307, 1988.
View at: Google Scholar

Copyright

Copyright © 2022 Do Thi Minh Xuan et al. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

1625

Downloads

1582

Citations