The role of microRNAs in metabolic interactions between viruses and their hosts

Polygonum cillinerve polysaccharide inhibits transmissible gastroenteritis virus by regulating microRNA-181

Transmissible gastroenteritis virus (TGEV) is an important pathogen capable of altering the expression profile of cellular miRNA. In this study, the potential of Polygonum cillinerve polysaccharide (PCP) to treat TGEV-infected piglets was evaluated through in vivo experiments. High-throughput sequencing technology was employed to identify 9 up-regulated and 17 down-regulated miRNAs during PCP-mediated inhibition of TGEV infection in PK15 cells. Additionally, miR-181 was found to be associated with target genes of key proteins in the apoptosis pathway. PK15 cells were treated with various concentrations of PCP following transfection with miR-181 mimic or inhibitor. Real-time PCR assessed the impact on TGEV replication, while electron microscopy (TEM) and Hoechst fluorescence staining evaluated cellular functionality. Western blot analysis was utilized to assess the expression of key signaling factors-cytochrome C (cyt C), caspase 9, and P53-in the apoptotic signaling pathway. The results showed that compared with the control group, 250 μg/mL PCP significantly inhibited TGEV gRNA replication and gene N expression (P < 0.01). Microscopic examination revealed uniform cell morphology and fewer floating cells in PCP-treated groups (250 and 125 μg/mL). TEM analysis showed no typical virus structure in the 250 μg/mL PCP group, and apoptosis staining indicated a significant reduction in apoptotic cells at this concentration. Furthermore, PCP may inhibit TGEV-induced apoptosis via the Caspase-dependent mitochondrial pathway following miR-181 transfection. These findings provide a theoretical basis for further exploration into the mechanism of PCP's anti-TGEV properties.

Cytomegalovirus microRNAs level determination in kidney recipients post transplantation

BACKGROUND

Human cytomegalovirus (CMV) can establish a latent infection with periodic or sporadic reactivation after the first infection happens. Primary and recurrent infection, results in different problems in patients with impaired or immature immune systems, such as kidney transplant recipients (KTRs). MicroRNAs (miRNAs, miRs) are important regulatory molecules in the outcome of CMV-infected KTRs. Therefore, in this study the expression level of CMV miRNAs were evaluated in active vs. latent CMV infected KTRs.

METHODS

Expression of viral miRNAs were studied in 61 KTRs which were divided into 30 active CMV and 31 latent CMV infected individuals. In order to study the expression level of selected miRNAs, SYBR Green Real-time PCR technique was exploited. Also, mature miRNAs expression level that were produced from one precursor, studied both in active and latent situations.

RESULTS

Among studied miRNAs' expression level, CMV miR-UL112-3p/5p, -UL22A-3p/5p, -US25-1-5p, -US25-2-3p/5p, -UL36-3p/5p and -UL70-3p showed significant increase in active CMV infected KTRs in comparison to latent ones. The ROC curve analysis results for miR-UL112-3p, -UL22A-3p, -US25-2-3p, -UL36-3p and -UL70-3p showed significant difference between two studied patient groups.

CONCLUSION

This study revealed an extremely high expression level in CMV miR-UL112-3p/5p, -UL22A-3p/5p, -US25-1-5p, -US25-2-3p/5p, -UL36-3p/5p and -UL70-3p in active CMV infected KTRs in comparison to latent ones. Further studies might help in finding the capability of miRNAs to differentiate active from latent stage of CMV infection in KTRs.

Extracellular Vesicles Regulated by Viruses and Antiviral Strategies

Extracellular vesicles (EVs), consisting of exosomes, micro-vesicles, and other vesicles, mainly originate from the multi-vesicular body (MVB) pathway or plasma membrane. EVs are increasingly recognized as a tool to mediate the intercellular communication and are closely related to human health. Viral infection is associated with various diseases, including respiratory diseases, neurological diseases, and cancers. Accumulating studies have shown that viruses could modulate their infection ability and pathogenicity through regulating the component and function of EVs. Non-coding RNA (ncRNA) molecules are often targets of viruses and also serve as the main functional cargo of virus-related EVs, which have an important role in the epigenetic regulation of target cells. In this review, we summarize the research progress of EVs under the regulation of viruses, highlighting the content alteration and function of virus-regulated EVs, emphasizing their isolation methods in the context of virus infection, and potential antiviral strategies based on their use. This review would promote the understanding of the viral pathogenesis and the development of antiviral research.

MicroRNA Mimics or Inhibitors as Antiviral Therapeutic Approaches Against COVID-19

Coronaviruses, such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) responsible for the coronavirus disease 2019 (COVID-19) pandemic, present a significant threat to human health by inflicting a wide variety of health complications and even death. While conventional therapeutics often involve administering small molecules to fight viral infections, small non-coding RNA sequences, known as microRNAs (miRNAs/miR-), may present a novel antiviral strategy. We can take advantage of their ability to modulate host-virus interactions through mediating RNA degradation or translational inhibition. Investigations into miRNA and SARS-CoV-2 interactions can reveal novel therapeutic approaches against this virus. The viral genomes of SARS-CoV-2, severe acute respiratory syndrome coronavirus (SARS-CoV), and Middle East respiratory syndrome coronavirus (MERS-CoV) were searched using the Nucleotide Basic Local Alignment Search Tool (BLASTn) for highly similar sequences, to identify potential binding sites for miRNAs hypothesized to play a role in SARS-CoV-2 infection. miRNAs that target angiotensin-converting enzyme 2 (ACE2), the receptor used by SARS-CoV-2 and SARS-CoV for host cell entry, were also predicted. Several relevant miRNAs were identified, and their potential roles in regulating SARS-CoV-2 infections were further assessed. Current treatment options for SARS-CoV-2 are limited and have not generated sufficient evidence on safety and efficacy for treating COVID-19. Therefore, by investigating the interactions between miRNAs and SARS-CoV-2, miRNA-based antiviral therapies, including miRNA mimics and inhibitors, may be developed as an alternative strategy to fight COVID-19.

Recent Progress on Relevant microRNAs in Autism Spectrum Disorders

Autism spectrum disorder (ASD) is a neurodevelopmental disorder whose pathogenesis is unclear and is affected by both genetic and environmental factors. The microRNAs (miRNAs) are a kind of single-stranded non-coding RNA with 20-22 nucleotides, which normally inhibit their target mRNAs at a post-transcriptional level. miRNAs are involved in almost all biological processes and are closely related to ASD and many other diseases. In this review, we summarize relevant miRNAs in ASD, and analyze dysregulated miRNAs in brain tissues and body fluids of ASD patients, which may contribute to the pathogenesis and diagnosis of ASD.

Understanding the cross-talk between host and virus in poultry from the perspectives of microRNA

In poultry, viral infections (e.g., Marek's disease virus, avian leukosis virus, influenza A virus, and so on) can cause devastating mortality and economic losses. Because viruses are solely dependent on host cells to propagate, they alter the host intracellular microenvironment. Thus, understanding the virus-host interaction is important for antiviral immunity and drug development in the poultry industry. MicroRNAs are crucial posttranscriptional regulators of gene expression in a wide spectrum of biological processes, including viral infection. Recently, microRNAs have been identified as key players in virus-host interactions. In this review, we will discuss the intricacies involved in the virus-host cross-talk mediated by host and viral microRNAs in poultry (i.e., chicken and ducks), as well as recent trends and challenges in this field. These findings may provide some insights into the rapidly developing area of research regarding viral pathogenesis and antiviral immunity in poultry production.

miR-142a-3p promotes the proliferation of porcine hemagglutinating encephalomyelitis virus by targeting Rab3a

Porcine hemagglutinating encephalomyelitis virus (PHEV) is a typical neurotropic coronavirus that mainly invades the central nervous system (CNS) in piglets and causes vomiting and wasting disease. Emerging evidence suggests that PHEV alters microRNA (miRNA) expression profiles, and miRNA has also been postulated to be involved in its pathogenesis, but the mechanisms underlying this process have not been fully explored. In this study, we found that PHEV infection upregulates miR-142a-3p RNA expression in N2a cells and in the CNS of mice. Downregulation of miR-142a-3p by an miRNA inhibitor led to a significant repression of viral proliferation, implying that it acts as a positive regulator of PHEV proliferation. Using a dual-luciferase reporter assay, miR-142a-3p was found to bind directly bound to the 3’ untranslated region (3’UTR) of Rab3a mRNA and downregulate its expression. Knockdown of Rab3a expression by transfection with an miR-142a-3p mimic or Rab3a siRNA significantly increased PHEV replication in N2a cells. Conversely, the use of an miR-142a-3p inhibitor or overexpression of Rab3a resulted in a marked restriction of viral production at both the mRNA and protein level. Our data demonstrate that miR-142a-3p promotes PHEV proliferation by directly targeting Rab3a mRNA, and this provides new insights into the mechanisms of PHEV-related pathogenesis and virus-host interactions.

A conserved miRNA-183 cluster regulates the innate antiviral response

Upon immune recognition of viruses, the mammalian innate immune response activates a complex signal transduction network to combat infection. This activation requires phosphorylation of key transcription factors regulating IFN production and signaling, including IFN regulatory factor 3 (IRF3) and STAT1. The mechanisms regulating these STAT1 and IRF3 phosphorylation events remain unclear. Here, using human and mouse cell lines along with gene microarrays, quantitative RT-PCR, viral infection and plaque assays, and reporter gene assays, we demonstrate that a microRNA cluster conserved among bilaterian animals, encoding miR-96, miR-182, and miR-183, regulates IFN signaling. In particular, we observed that the miR-183 cluster promotes IFN production and signaling, mediated by enhancing IRF3 and STAT1 phosphorylation. We also found that the miR-183 cluster activates the IFN pathway and inhibits vesicular stomatitis virus infection by directly targeting several negative regulators of IRF3 and STAT1 activities, including protein phosphatase 2A (PPP2CA) and tripartite motif-containing 27 (TRIM27). Overall, our work reveals an important role of the evolutionarily conserved miR-183 cluster in the regulation of mammalian innate immunity.

Identification of potential microRNA markers related to Crimean-Congo hemorrhagic fever disease

Crimean-Congo hemorrhagic fever (CCHF) is a tick-borne disease caused by the arbovirus Crimean-Congo hemorrhagic fever virus (CCHFV). The CCHFV has a single-stranded RNA genome of negative sense. MicroRNAs (miRNAs) are key players in virus-host interactions and viral pathogenesis. We investigated the miRNA gene expression profiles in patients with CCHF using microarray for the first time in the world. Microarray analysis was performed using mirBase Ver 21 (Agilent Technologies, Santa Clara, CA). All statistical analyses were performed across the case-control, fatal-control, and fatal-nonfatal case groups using Genespring (Ver 3.0). Fifteen miRNAs were statistical significant in patients with CCHF compared with the controls (5 were upregulated, 10 were downregulated). Seventy-five and sixty-six miRNAs are in fatal compared with control and nonfatal case, respectively (fold change ([FC] ≥50) were statistically significant. In this study, the target genes of important miRNAs were identified and Gene Ontology analyses were performed across all groups. As a result of this study, we propose that the detection of miRNAs in patients with CCHF will allow the determination of therapeutic targets in diseases. CCHF is an important public health problem that can often be fatal. In this study, we investigated miRNA expression in case-control, fatal-control, and fatal-nonfatal case groups. Significant miRNAs associated with fatality were detected in CCHF. This study will serve as a source of data for the development of an antagomir-based therapy against CCHF using miRNAs in the future.

miREM: an expectation-maximization approach for prioritizing miRNAs associated with gene-set

Background

The knowledge of miRNAs regulating the expression of sets of mRNAs has led to novel insights into numerous and diverse cellular mechanisms. While a single miRNA may regulate many genes, one gene can be regulated by multiple miRNAs, presenting a complex relationship to model for accurate predictions.

Results

Here, we introduce miREM, a program that couples an expectation-maximization (EM) algorithm to the common approach of hypergeometric probability (HP), which improves the prediction and prioritization of miRNAs from gene-sets of interest. miREM has been made available through a web-server (https://bioinfo-csi.nus.edu.sg/mirem2/) that can be accessed through an intuitive graphical user interface. The program incorporates a large compendium of human/mouse miRNA-target prediction databases to enhance prediction. Users may upload their genes of interest in various formats as an input and select whether to consider non-conserved miRNAs, amongst filtering options. Results are reported in a rich graphical interface that allows users to: (i) prioritize predicted miRNAs through a scatterplot of HP p-values and EM scores; (ii) visualize the predicted miRNAs and corresponding genes through a heatmap; and (iii) identify and filter homologous or duplicated predictions by clustering them according to their seed sequences.

Conclusion

We tested miREM using RNAseq datasets from two single “spiked” knock-in miRNA experiments and two double knock-out miRNA experiments. miREM predicted these manipulated miRNAs as having high EM scores from the gene set signatures (i.e. top predictions for single knock-in and double knock-out miRNA experiments). Finally, we have demonstrated that miREM predictions are either similar or better than results provided by existing programs.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2292-1) contains supplementary material, which is available to authorized users.

piRNA Profiling of Dengue Virus Type 2-Infected Asian Tiger Mosquito and Midgut Tissues

The Asian tiger mosquito, Aedes albopictus, is a competent vector for the majority of arboviruses. The mosquito innate immune response is a primary determinant for arthropod-borne virus transmission, and the midgut is the first barrier to pathogen transmission. Mosquito antiviral immunity is primarily mediated by the small interfering RNA pathway. However, the roles that the P-element induced wimpy testis (PIWI)-interacting RNA (piRNA) pathway play in antiviral immunity in Ae. albopictus and its midgut still need further exploration. This study aimed to explore the profiles of both viral-derived and host-originated piRNAs in the whole body and midgut infected with Dengue virus 2 (DENV-2) in Ae. albopictus, and to elucidate gene expression profile differences of the PIWI protein family between adult females and their midguts. A deep sequencing-based method was used to identify and analyze small non-coding RNAs, especially the piRNA profiles in DENV-2-infected Ae. albopictus and its midgut. The top-ranked, differentially-expressed piRNAs were further validated using Stem-loop qRT-PCR. Bioinformatics analyses and reverse-transcription PCR (RT-PCR) methods were used to detect PIWI protein family members, and their expression profiles. DENV-2 derived piRNAs (vpiRNA, 24–30 nts) were observed in both infected Ae. albopictus and its midgut; however, only vpiRNA in the whole-body library had a weak preference for adenine at position 10 (10A) in the sense molecules as a feature of secondary piRNA. These vpiRNAs were not equally distributed, instead they were derived from a few specific regions of the genome, especially several hot spots, and displayed an obvious positive strand bias. We refer to the differentially expressed host piRNAs after DENV infection as virus-induced host endogenous piRNAs (vepiRNAs). However, we found that vepiRNAs were abundant in mosquito whole-body tissue, but deficient in the midgut. A total of eleven PIWI family genes were identified in Ae. albopictus; however, only AalPiwi5–7 and AalAgo3(1–2) were readily detected in the midgut. The characteristics of piRNAs in DENV-2-infected Ae. albopictus adult females were similar to those previously described for flavivirus infections but were not observed in the midgut. The reduced levels of vepiRNAs and incomplete expression of PIWI pathway genes in midgut samples from DENV-2-infected Ae. albopictus suggests that viral regulation of host piRNAs may not be an important factor in the midgut.

The expression profile of human peripheral blood mononuclear cell miRNA is altered by antibody-dependent enhancement of infection with dengue virus serotype 3

Background

Antibody-dependent enhancement (ADE) of dengue virus (DENV) infection has been identified as the main risk factor for severe dengue disease, although the underlying mechanisms leading to severe dengue fever remain unclear. MicroRNAs (miRNAs) participate in numerous pathological and biological processes, including host responses to viral infections.

Method

Here, we aimed to investigate the differences in miRNA expression patterns in human peripheral blood mononuclear cells (PBMCs) infected with DENV-3 and DENV-3-ADE at various time points employing high-throughput sequencing.

Results

According to miRNAs high-throughput sequencing, a total of 50 known miRNAs exhibited significant differences. GO (Gene Ontology) and pathway analysis of the predicted targets showed enrichment in the regulation of transcription, including multicellular organismal development, DNA-dependent transcription, negative regulation of cell differentiation and transcription. Afterwards, regulatory networks of miRNA predicted targets, miRNA transcription factors, miRNA pathways and miRNA GOs were formulated to expose the complex regulatory mechanisms of miRNAs during the infection phase. Finally, we analyzed hierarchical GO categories of the predicted targets involved in the MAPK signaling pathway, the cGMP-PKG signaling pathway, the cAMP signaling pathway, the endocytosis effect, and our analyses indicated that innate and adaptive immunity following DENV-3 and DENV-3-ADE infections may be signally distinct.

Conclusion

Our results demonstrate a novel describing miRNA expression profiles in human PBMCs with DENV-3 and DENV-3-ADE infections using high-throughput sequencing. Our findings could provide a beneficial basis for further studies on the regulatory roles of miRNAs relevant to the different immune responses caused by DENV-3 and DENV-3-ADE infections.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-0963-1) contains supplementary material, which is available to authorized users.

An interferon-independent lncRNA promotes viral replication by modulating cellular metabolism

Viruses regulate host metabolic networks to improve their survival. The molecules that are responsive to viral infection and regulate such metabolic changes are hardly known, but are essential for understanding viral infection. Here we identify a long noncoding RNA (lncRNA) that is induced by multiple viruses, but not by type I interferon (IFN-I), and facilitates viral replication in mouse and human cells. In vivo deficiency of lncRNA-ACOD1 (a lncRNA identified by its nearest coding gene Acod1, aconitate decarboxylase 1) significantly attenuates viral infection through IFN-I-IRF3 (interferon regulatory factor 3)-independent pathways. Cytoplasmic lncRNA-ACOD1 directly binds the metabolic enzyme glutamic-oxaloacetic transaminase (GOT2) near the substrate niche, enhancing its catalytic activity. Recombinant GOT2 protein and its metabolites could rescue viral replication upon lncRNA-ACOD1 deficiency and increase lethality. This work reveals a feedback mechanism of virus-induced lncRNA-mediated metabolic promotion of viral infection and a potential target for developing broad-acting antiviral therapeutics.

Expression of miR-155 associated with Toll-like receptors 3, 7, and 9 transcription in the olfactory bulbs of cattle naturally infected with BHV5

Bovine herpesvirus 5 (BHV5) infection of young cattle is frequently associated with fatal neurological disease and, as such, represents an attractive model for studying the pathogenesis of viral-induced meningoencephalitis. Following replication in the nasal mucosa, BHV5 invades the central nervous system (CNS) mainly through the olfactory pathway. The innate immune response triggered by the host face to virus replication through the olfactory route is poorly understood. Recently, an upregulation of conserved pathogen-associated molecular pattern, as Toll-like receptors (TLRs), has been demonstrated in the CNS of BHV5 experimentally infected cows. A new perspective to understand host-pathogen interactions has emerged elucidating microRNAs (miRNAs) network that interact with innate immune response during neurotropic viral infections. In this study, we demonstrated a link between the expression of TLRs 3, 7, and 9 and miR-155 transcription in the olfactory bulbs (OB) of 16 cows suffering from acute BHV5-induced neurological disease. The OBs were analyzed for viral antigens and genome, miR-155 and TLR 3, 7, and 9 expression considering three major regions: olfactory receptor neurons (ORNs), glomerular layer (GL), and mitral cell layer (ML). BHV5 antigens and viral genomes, corresponding to glycol-C gene, were detected in all OBs regions by fluorescent antibody assay (FA) and PCR, respectively. TLR 3, 7, and 9 transcripts were upregulated in ORNs and ML, yet only ORN layers revealed a positive correlation between TLR3 and miR-155 transcription. In ML, miR-155 correlated positively with all TLRs studied. Herein, our results evidence miR-155 transcription in BHV5 infected OB tissue associated to TLRs expression specifically ORNs which may be a new window for further studies.

Viral Ubiquitin Ligase Stimulates Selective Host MicroRNA Expression by Targeting ZEB Transcriptional Repressors

Infection with herpes simplex virus-1 (HSV-1) brings numerous changes in cellular gene expression. Levels of most host mRNAs are reduced, limiting synthesis of host proteins, especially those involved in antiviral defenses. The impact of HSV-1 on host microRNAs (miRNAs), an extensive network of short non-coding RNAs that regulate mRNA stability/translation, remains largely unexplored. Here we show that transcription of the miR-183 cluster (miR-183, miR-96, and miR-182) is selectively induced by HSV-1 during productive infection of primary fibroblasts and neurons. ICP0, a viral E3 ubiquitin ligase expressed as an immediate-early protein, is both necessary and sufficient for this induction. Nuclear exclusion of ICP0 or removal of the RING (really interesting new gene) finger domain that is required for E3 ligase activity prevents induction. ICP0 promotes the degradation of numerous host proteins and for the most part, the downstream consequences are unknown. Induction of the miR-183 cluster can be mimicked by depletion of host transcriptional repressors zinc finger E-box binding homeobox 1 (ZEB1)/-crystallin enhancer binding factor 1 (δEF1) and zinc finger E-box binding homeobox 2 (ZEB2)/Smad-interacting protein 1 (SIP1), which we establish as new substrates for ICP0-mediated degradation. Thus, HSV-1 selectively stimulates expression of the miR-183 cluster by ICP0-mediated degradation of ZEB transcriptional repressors.

MicroRNA-125a-5p plays a role as a tumor suppressor in lung carcinoma cells by directly targeting STAT3

Increasing evidence supports that the dysregulation of microRNA expression plays an important role in the process of tumor occurrence and development. Studies have found that mir-125a-5p expression was downregulated in a variety of tumors, but the effects and mechanism of mir-125a-5p in lung cancer are still unclear. The aim of this study is to detect the expression of mir-125a-5p in lung cancer tissues and lung cancer cell lines and to explore the effects of mir-125a-5p on the biological characteristics of lung cancer cells; thus, this study aims to provide new methods and new strategies for the treatment of lung cancer. The result from quantitative reverse transcription polymerase chain reaction showed that the expression of miR-125a-5p was significantly lower in lung cancer tissues and lung cancer cell lines (95-D, A549, HCC827, and NCI-H1299) than that in normal tissue adjacent to lung cancer or normal human bronchial epithelial cells. In order to explore the function and mechanism of mir-125a-5p in lung cancer cells, miR-125a-5p mimic or mir-125a-5p inhibitor was transfected into A549 cells. Mir-125a-5p displayed an obvious upregulation in A549 cells transfected with miR-125a-5p and an obvious downregulation in A549 cells transfected with mir-125a-5p inhibitor compared to that in A549 cells transfected with control miRNA. 3-(4,5-Dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide, BrdU staining, flow cytometry, and Transwell assay showed that the upregulation of miR-125a-5p could significantly decrease the cell viability, proliferation, and invasion of lung cancer cells and increase apoptosis of lung cancer cells. The downregulation of miR-125a-5p provided very contrasting results. Computational algorithms predicted that the STAT3 is a target of miR-125a-5p. Here, we validated that miR-125a-5p could directly bind to the 3'-untranslated region of STAT3, and miR-125a-5p overexpression could significantly inhibit the protein expression of STAT3. These results suggested that mir-125a-5p can regulate the expression of STAT3 in lung cancer cells. To further verify whether mir-125a-5p can play a biological role through regulating STAT3, 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide, flow cytometry, and Transwell analysis demonstrated that overexpression of STAT3 can reverse the cells' biological effects induced by mir-125a-5p overexpression. Mir-125a-5p downregulated in lung cancer tissue and cell lines can negatively regulate STAT3 protein expression. Taken together, mir-125a-5p inhibited the proliferation and invasion of lung cancer cells and facilitated lung cancer cell apoptosis through suppressing STAT3. Enhancing the expression of miR-125a-5p is expected to benefit the therapy for the patients with lung cancer.

Different microRNA alterations contribute to diverse outcomes following EV71 and CA16 infections: Insights from high-throughput sequencing in rhesus monkey peripheral blood mononuclear cells

Enterovirus 71 (EV71) and Coxsackievirus A16 (CA16) are the predominant pathogens of hand, foot, and mouth disease (HFMD). Although these viruses exhibit genetic homology, the clinical manifestations caused by the two viruses have some discrepancies. In addition, the underlying mechanisms leading to these differences remain unclear. microRNAs (miRNAs) participate in numerous biological or pathological processes, including host responses to viral infections. Here, we focused on differences in miRNA expression patterns in rhesus monkey peripheral blood mononuclear cells (PBMCs) infected with EV71 and CA16 at various time points using high-throughput sequencing. The results demonstrated that 106 known and 13 novel miRNAs exhibited significant differences, and 32 key miRNAs among them for target prediction presented opposite trends in the EV71- and CA16-infected samples. GO and pathway analysis of the predicted targets showed enrichment in 14 biological processes, 10 molecular functions, 8 cellular components and 104 pathways. Subsequently, regulatory networks of miRNA-transcription factors, miRNA-predicted targets, miRNA-GOs and miRNA-pathways were constructed to reveal the complex regulatory mechanisms of miRNAs during the infection phase. Ultimately, we analysed hierarchical GO categories of the predicted targets involved in immune system processes, which indicated that the innate and adaptive immunity following EV71 and CA16 infections may be remarkably distinct. In conclusion, this report is the first describing miRNA expression profiles in PBMCs with EV71 and CA16 infections using high-throughput sequencing. Our findings could provide a valuable basis for further studies on the regulatory roles of miRNAs related to the different immune responses caused by EV71 and CA16 infections.