Hand, foot, and mouth disease (HFMD) is a common endemic childhood disease worldwide, particularly in Asia. Generally, it is triggered by two major causative agents: enterovirus 71 (EV71) and coxsackievirus A16 (CA16) (Kaminska et al. 2013; Koh et al. 2016), followed by coxsackievirus A6 (CA6) and coxsackie virus A10 (CA10) (Wang et al. 2018). This disease majorly affects children under five years old (Omana-Capeda et al. 2016). The common symptoms include fever, rashes on the volar regions of the hands and feet, herpangina, and difficulties in eating and drinking. Severe circumstances can exhibit potentially fatal complications involving the nervous system, such as brain stem encephalitis, or cardiopulmonary systems, such as pulmonary edema. Severe cases may show drastic progression and die from complicated symptoms.
There have been limited tools for effective diagnosis of HFMD. Although enterovirus infection triggers various pathological responses, the exact mechanisms underlying HFMD development remain largely unknown. For example, the regulatory roles played by EV71/CA16 infection towards endothelial cells and neural systems remain unclear. Previous studies have explored the disordered signals in different aspects, such as inflammatory profiles in cytokine expression (Teo et al. 2018; Linghua et al. 2019), long non-coding RNA (lncRNA) profiles (Meng et al. 2017), and immune cell changes (Wang et al. 2014). It is conceivable and supported that patients with HFMD may exhibit changes in expression profiles of microRNAs and mRNAs, mostly derived from blood samples (Cui et al. 2011; Hu et al. 2016; Yee et al. 2016; Zhu et al. 2016; Song et al. 2017; Hu et al. 2018; Jia et al. 2018; Li et al. 2018; Mi et al. 2018; Song et al. 2018). In addition, a comprehensive understanding of HFMD-related expression profiles and identification key markers may provide substantial diagnostic and prognostic values. To dig the clinical significance of micro-RNA and mRNA changes under HFMD virus infection, we conducted this bioinformatics analysis to clarify crucial differentially expressed genes (DEGs), functional and pathway enrichment, and the relative regulatory network, using four published datasets. This analysis may provide deeper insight into the mechanism of HFMD pathological development and novel strategies to prevent HFMD outbreaks.
Differentially expressed miRNAs after HFMD virus infection in 16HBE cells based on the GSE85829 dataset.
(A) Left: EV71 infection-induced 25 DE-miRNAs; Right: CA16 infection-induced 13 DE-miRNAs.
(B) The Venn diagram of 7 common DE-miRNAs between EV71 and CA16 infection.
Differentially expressed miRNAs after HFMD virus infection in the HUVEC cell line based on the GSE94551 dataset.
(A) Left: 39 DE-miRNAs were observed after EV71 infection; Right: 99 DE-miRNAs were found after CA16 infection.
(B) The Venn diagram of 7 common DE-miRNAs between EV71 and CA16 infection.
(C) The Venn diagram of common DE-miRNAs among four groups of HFMD virus infection.
MicroRNAs in serum exosomes from HFMD patients were analyzed based on the GSE52780 dataset.
(A) Heatmap of 258 DE-miRNAs was found, including 85 up-regulated and 173 down-regulated ones.
(B) The top ten up-regulated and down-regulated DE-miRNAs.
GO and KEGG enrichment of key miRNAs and DE-mRNAs. Based on the above five key miRNAs, GO functional, and KEGG pathway enrichment was performed using the miRPath v. 3 database.
(A) GO enrichment based on the intersection targets of the five key miRNAs (only four miRNAs were included in the miRPath database).
(B) KEGG enrichment based on the intersection targets of the four key miRNAs.
(C) GO enrichment verification through the union of key miRNA related GO terms, presented as a heatmap.
(D) Verification of KEGG pathway through the union of key miRNA related KEGG pathways presented as a heatmap.
(E) The enriched GO functions of up-regulated DEGs in SH-SY5Y cells infected by EV71 from the GSE45589 dataset.
(F) The enriched GO functions of down-regulated DEGs in SH-SY5Y cells infected by EV71 from the GSE45589 dataset.
(G) Enriched KEGG pathways in SH-SY5Y cells infected by EV71 from the GSE45589 dataset.
Common differential mRNAs and PPI network.
(A) Using the TarBase V8 tool, we acquired 1,520 potential targets (mRNA) from the 5 key DE-miRNAs, among which 11 DEGs were also included by the 159 DE-mRNAs in the GSE45589 dataset: MACF1, MARS, SF3B3, SMARCC1, BRMS1L, SMC1A, SPHK2, LIG1, CSF3, CYR61, and FGFR1OP.
(B) GO functional analysis showed three terms might be enriched according to these DEGs.
(C) These common DEGs showed a PPI network mainly connected by SMC1A, SMARCC1, SF3B3, LIG1, and BRMS1L.
(D) The PPI network constructed by the 159 DE-mRNAs in the GSE45589 dataset (the isolated nodes were removed).
In this study, we used five datasets to identify key RNA members in HFMD development. We were interested in five key miRNAs, 11 mRNAs, and several important GO and KEGG enrichment after filtering progressively. Our results might provide some theoretical perspectives about HFMD development and a potential strategy in its early warning.
At the miRNA level, several potentially useful markers have been proposed in clinical diagnosis. A survey in Singapore reported a 6-miRNA scoring model which predicts HFMD with an overall accuracy of 85.11% in the training set and 92.86% in the blinded test set, and circulating salivary miRNA hsa-miR-221 (downregulated in that work) was regarded as a highly validated marker (Mi et al. 2018). Song et al. have applied rhesus monkey peripheral blood mononuclear cells to search DE-miRNAs, and they identified 13 novel DE-miRNAs with 2501 targets (Song et al. 2018). Zhu et al. (2016) performed the microarray examination and noticed 27 DE-miRNAs (15 up-regulated and 12 downregulated) associated with CA16 and EV71 infection. There were some other important findings of specific miRNAs. MiR-1303 has been known to promote CNS lesions following CA16 infections by targeting MMP9 (Song et al. 2018). EV71 can evade the immune surveillance system to proliferate by activating miR-21 (Feng et al. 2017), antagonize the antiviral activity of host STAT3 and IL-6R through miR-124 (Chang et al. 2017), and induce autophagy by regulating miR-30a to promote viral replication (Fu et al. 2015). So far, there have been very few direct reports about the relationship between the five key miRNAs and HFMD. Only one study had surveyed the miRNA expression profile in the exosome of HFMD patients (Jia et al. 2014), and it revealed that the expression level of three miRNAs (miR-671-5p, miR-16-5p, and miR-150-3p) was significantly different between mild HFMD, extremely severe SHFMD, and the healthy controls. We also noticed that miR-671-5p was among the key miRNAs in HFMD.
The PPI network implied that five targets, SMC1A, SMARCC1, SF3B3, LIG1, and BRMS1L, might play the most crucial roles during HFMD progression. However, none of them has been paid enough attention to date, and they are worth more concerns in further researches. Taken different datasets together, we found three common enriched GO terms between miRNA-derived prediction and mRNA-derived analysis: biosynthetic process, cytosol, and nucleoplasm; and a common KEGG pathway, cell cycle, was noticed (Fig. 4). These functions and pathways suggest that HFMD viruses strongly drive the host proliferation. This fact could be also be supported by the GO functional enrichment constructed based on the 11 key mRNAs (Fig. 5B), which exhibited that positive regulation of cell proliferation was the most enriched functional term. This process may contribute to virus amplification and also be a homeostasis response to fight against virus invasion, particularly for epithelial cells. However, the definite mechanism needs more evidence to unravel.
However, this study has some limitations. First, when we probed the key roles, common DEGs were screened between EV71 and CA16 infection. However, despite belonging to the same genus, Enterovirus, these two viruses do not necessarily have similar effects. For example, Chinese scholars identified that miR-4516 presented down-regulation in EV71 infection and upregulation in CA16 infection, and it was an important regulator of intercellular junctions by targeting PVRL1 (Hu et al. 2016). Liu et al. had analyzed miRNA profiles and acquired diverse outcomes induced by EV71 and CA16 infection (Song et al. 2017). The inconsistency was also shown in Fig. 2B. There were two common up-regulated miRNAs, two common down-regulated miRNAs, and three inconsistent ones (miR-502-5p, miR-503-5p, and miR-542-3p). Besides, the direct regulatory relationship between five key miRNAs and 11 key mRNAs had not been validated in the present study. Our further efforts would focus on the construction of diagnostic and prognostic models based on the large real-world sample using these miRNA and mRNA factors.
Our results shed light on the potentially crucial roles of the five miRNAs, 11 coding genes, and several functions and pathways in HFMD. A combination of these roles may benefit the early diagnosis and treatment of HFMD.