Abstract

Role of microRNA (miRNA) has been highlighted in pathogen-host interactions recently. To identify cellular miRNAs involved in the host response to enterovirus 71 (EV71) infection, we performed a comprehensive miRNA profiling in EV71-infected Hep2 cells through deep sequencing. 64 miRNAs were found whose expression levels changed for more than 2-fold in response to EV71 infection. Gene ontology analysis revealed that many of these mRNAs play roles in neurological process, immune response, and cell death pathways, which are known to be associated with the extreme virulence of EV71. To our knowledge, this is the first study on host miRNAs expression alteration response to EV71 infection. Our findings supported the hypothesis that certain miRNAs might be essential in the host-pathogen interactions.

1. Introduction

Hand, foot, and mouth disease (HFMD), a common febrile illness in children, is usually caused by human enteroviruses. Enterovirus 71 (EV71) and coxsackievirus A16 (CA16) are two major causative agents of HFMD. EV71 and CA16 infections manifesting as lesions on the skin and oral mucosa are clinically similar, but EV71 infection is more frequently associated with serious neurological diseases such as aseptic meningitis, encephalitis, and acute flaccid paralysis and fatalities [14] while the CA16-associated HFMD has a milder outcome [5]. More than 500,000 HFMD cases caused by EV71 were reported nationwide, including 176 fatal cases in China since March 2008 [6]. Elucidating the cellular events following EV71 infection will facilitate the development of strategies to prevent and treat this virus. However, the molecular mechanisms of the host response to EV71 infection are not completely understood.

microRNAs (miRNAs) have emerged as key regulators in many biological processes, from development to defense, at almost all organismal levels through mRNA degradation or translational repression of their targets [7, 8]. Recently, their role has been highlighted in pathogen-host interactions. Not only the miRNA encoded by viral genomes but host encoded miRNAs have been found participating in host-virus interactions. They can stimulate as well as suppress viral infections. For example, liver-specic miR-122 is an indispensable factor in supporting hepatitis C virus (HCV) replication [9], whereas, miR-125b and miR-223 directly target human immunodeficiency virus-1 (HIV-1) mRNA, thereby attenuating viral gene expression in resting CD4+ T cells [10]. It has been showed that Epstein-Barr virus (EBV) encodes five miRNAs in its large DNA genome. These miRNAs can potentially regulate several genes encoding proteins involved in apoptosis, cell proliferation, signal transduction, transcription regulation, and immune response [11]. In contrast, the global changes, miRNAs expression during EV71 infection have not yet to be extensively elucidated. To determine which cellular microRNAs play a role in the host response to enterovirus infection, in this study, we performed a comprehensive miRNA proling in EV71-infected Hep2 cells through deep sequencing.

2. Materials and Methods

2.1. Cell Culture and Virus Infection

Human epidermoid carcinoma (Hep2) cells were grown in RPMI 1640 medium (Invitrogen) supplemented with 10% fetal calf serum (FCS). When the Hep2 cells were grown to 70% confluence in 25 cm2 flasks, they were infected with EV71 at a multiplicity of infection (m.o.i.) of 0.03 50% tissue culture infectious doses (TCID50) and maintained after infection at in RPMI 1640 medium with 2% FCS.

2.2. Isolation of RNA and miRNA

Hep2 cells were infected with EV71 as described above. At 6 h, 24 h, 48 h, 72 h, and 96 h post infection, total RNA or miRNA was extracted from cells using the mirVana kit according to the manufacturer’s protocol (Ambion). RNA was also extracted from noninfected control cells at the same times. The quality and the concentration of the RNA samples were monitored by gel electrophoresis and absorbance at A260/280 ratio.

2.3. Library Construction, SOLiD Sequencing and Analysis

The miRNA sequencing library construction followed standard procedure of SOLiD small RNA expression kit (Applied Biosystem). All SOLiD run parameters followed standard Applied Biosystems protocols. Different barcodes were introduced to two samples in the polymerase chain reaction of library construction, and all the samples were sequenced in a single sequencing run.

2.4. Data Analysis

SOLiD data were first analyzed by SOLiD System Small RNA Analysis Pipeline Tool (RNA2MAP). The miRBase sequences (Sanger) of human being were downloaded from miRBase (http://www.mirbase.org/).The number of bases to use when generating initial seeds locations was 18 with a tolerance of 3 mismatches. After extension step, at most 6 mismatches were allowed in full length mapping.

Potential conserved target genes of differentially expressed miRNAs were firstly predicted by targetscan (http://www.targetscan.org/) [1214]. Target genes of some miRNAs, such as mir-1972, mir-1974, mir-1975, mir-1979, and mir-764, could not be predicted in targetscan database, we further predicted those miRNAs targets using DIANA-microT v3.0 [15, 16]. In brief, each differentially expressed miRNA was submitted to targetscan individually and all of its targets predicted in targetscan or microT v3.0 were used for the following gene ontology (GO) analysis (http://www.babelomics.bioinfo.cipf.es/). All targets of induced and repressed miRNAs were submitted to FatiGO program [17]. Functional category enrichment based on the GO terms was evaluated on the targets of these differentially expressed miRNAs.

2.5. Confirmation of Differentially Expressed miRNAs by Real-time Quantitative RT-PCR

To confirm the expression of miRNAs by deep sequencing approach, stem-loop quantitative RT-PCR (qRT-PCR) was performed. In brief, cDNA was synthesized from total RNA by using AMV reverse transcriptase (TaKaRa). The 20  l reactions were incubated for 15 min at C, 30 min at C and 5 min at C, and then held at C. Subsequently, real-time quantification was performed using an Applied Biosystems 7500 Sequence Detection system (Applied Biosystems). The 20  l PCR reactions included 1  l RT-PCR product, 10  l Premix Ex Taq (TaKaRa), and 1  l SYBR green (Invitrogen). The reactions were incubated in a 96-well optical plate at C for 10 min, followed by 40 cycles of C for 15 s and C for 1 min. All reactions were run in triplicate. After reaction, the threshold cycle value (CT) data were determined using default threshold settings, and the mean CT was determined from the duplicate PCRs. Human small nuclear RNA U6 was used for normalization. The expression levels of miRNAs were measured in terms of CT and normalized to U6 using [18].

2.6. Statistical Analysis

We firstly used the Z test to determine the statistical significance of the differences between the two libraries [19]. This approach is to look at the number of copies of a specic miRNA per cell as a fraction or proportion of the total number of miRNA molecules in that cell. The same proportion of specic tags should be present in the miRNA library of all sequenced tags. In this test a false discovery rate less than 5% was selected. miRNAs were considered significantly altered only when they fullled three criteria: (1) mean fold change 2 or 0.5, (2) having at least 10 copies by SOLiD sequencing, and (3) Z score 1.96 or −1.96.

3. Results

3.1. Replication of EV71 in Hep2 Cells

Cytopathic effect (CPE) was first observed at 24 h post infection and progressed to moderate and severe CPE at 72 and 96 h, respectively. miRNA extracted from cells after 72 h infection was used for SOLiD deep sequencing.

3.2. Deep Sequencing

A total of 411,419 and 479,414 filtered high quality reads were obtained from infected and control cells by deep sequencing, respectively. After filtrating reads contaminated by rRNA, tRNA, snRNA, and snoRNA, 35,272 and 78,143 reads were obtained from infected and non-infected cells, respectively. Out of these reads, 11,188 of these high quality reads were exact matches while other 68,041 reads were loose matches to known human miRNAs (Table 1). Loose matches were defined by sequence reads that aligned with human miRNA consensus sequence with 1–4 mismatches. These may represent sequencing errors (when occurring in low copy numbers), mutations, and/or RNA editing events.

The size distribution of sequence reads showed that the majority of miRNAs was 18–25 for both libraries (>90%), with 22 nt small RNA being the most abundant (Figure 1), which is within the typical size range of human miRNAs.

After reads were compared with an miRBase database (release 14.0), 569 miRNAs were detected in EV71-infected cells while 540 miRNAs were detected in non-infected control cells.

3.3. Aberrant miRNAs Expression in EV71-Infected Cells

On the basis of differentially expressed miRNA, we found 64 that miRNAs were differentially expressed between infected and non-infected cells. More miRNAs (42 out of 64 miRNAs) were upregulated than down-regulated during EV71 infection in Hep2 cells (Tables 2 and 3).

3.4. Confirmation of Differentially Expressed miRNAs

Quantitative RT-PCR assays were used to confirm the expression pattern of differentially expressed miRNAs in Hep2 cells. There was general consistency between quantitative RT-PCR assay and deep sequence analysis in four miRNAs (miR-1246, miR-1237, miR-30a, and miR-222) in terms of directions of regulation and significance. Specifically, there was a 1.92-fold upregulation (7.35-fold in deep sequencing analysis) in miR-1246, 1.58-fold upregulation (2.58-fold in deep sequencing analysis) in mir-1237, ( ), 2.45-fold down-regulation (2.94-fold in deep sequencing analysis) in miR-30a, and 1.42-fold down-regulation (2.02-fold in deep sequencing analysis) in miR-222.

3.5. Gene Ontology Analysis

Targets were predicted for all identified differentially expressed miRNA families. In total 5765 unique target genes were predicted for 64 of the differentially expressed miRNAs ( see Table S1 in Supplementray Material available online at doi:115/2010/425939). On the basis of the biological functions described by FatiGO program (http://www.babelomics.bioinfo.cipf.es/), these target genes can be grouped into 72 categories (S2). The top 30 Gene Ontology terms are shown in Figure 2. The majority of targets fall into the category of metabolic process, regulation of biological process, and cell communication indicating intense biological change in Hep2 cells after EV71 infection. Several other groups contain genes regulating death (255 target genes including 239 apoptosis-related target genes), neurological process (211 target genes), and immune response (120 target genes).

Hundreds of altered transcripts in response to EV71 infection were found in two transcriptomic studies [20, 21]. To determine whether miRNAs might be modulators of mRNAs that were differentially expressed, we investigated whether those differentially expressed mRNAs are enriched for predicted targets and interrogated their inversely correlated targets for functional associations. To our surprise, targets for differentially expressed miRNAs or miRNA families in the present study account for 12.1% (19 out of 157) of the differentially expressed mRNAs in one study [20] whereas they account for 22.4% (13 out of 58) of the transcripts that were significantly altered in another study (Tables 4 and 5) [21].

4. Discussion

Many technologies have been developed for miRNA profiling, including real-time quantitative RT-PCR [22], northern blotting [23, 24], and microarray analyses based on either direct hybridization or hybridization coupled with enzymatic extension [25, 26]. These methods have been used successfully in a variety of studies. But they still have some technical limitations. For example, first two methods were not high-throughput while microarray method needs large amounts of starting materials. Recent progress in high-throughput sequencing technologies allows deep sequencing of large libraries of short RNAs [2730]. The longest reads are obtained by the 454 technology, which currently gives reads of about 400 base pairs (bp). However, this technology yields much less reads than other techniques (about 400 000 per sample) [27]. The Solexa (Illumina) and SOLiD platform (Applied Biosystems) generates shorter reads (up to 35 bp) but yields 1–3 million reads per sample [28, 29]. The other high-throughput technique, such as massively parallel sequencing (MPSS), gives more reads than Solexa but the reads are even shorter, only 17 bp [30]. In the present study, we used SOLiD platform to explore differential miRNA proling of Hep2 cells response to EV71 infection. We obtained about 0.4 million reads in both samples, well below the capacity of SOLiD. One possibility is that we used cell line as infection model but not tissue. Numbers of small RNAs were not expressed in Hep2 cells. Even now, we found that there was general consistency between quantitative RT-PCR assay and deep sequence analysis. So, SOLiD deep sequencing successfully revealed miRNA proling in EV71-infected and control Hep2 cells.

Based on a comprehensive examination of miRNA expression from EV71-infected and non-infected control Hep2 cells, we identified 64 miRNAs that were differentially expressed, with most of them (65.6%) upregulated in EV71-infected cells. The altered patterns of cellular miRNAs we observed for EV71-infected cells are similar in some respects to changes seen for cells infected with Hepatitis C virus (HCV) [31] but contrast with those with Epstein-Barr virus (EBV) [32] and human cytomegalovirus (HCMV) in which more miRNAs were dow-regulated in response to virus infection [33]. Most of the miRNAs affected by EV71 are different from those affected by HCV, EBV, and HCMV. We observed miR-636 and miR-584 expressed in only infected or non-infected cells. Though their functions during virus infection have not been explored, their predicted target genes have been identified to be involved in virus entry, replication and propagation. For example, reticulon 3 (RTN3), one predicted target gene of miR-636, can bind the 2C protein of enterovirus 71 and is required for viral replication, [34]. Poliovirus receptor-related 1 (herpesvirus entry mediator C) (PVRL1, also known as nectin-1), another target gene of miR-636, can serve as receptor for herpes simplex virus and pseudorabies virus entry [35]. Abhydrolase domain containing 2 (ABHD2), one predicted target gene of miR-584, is essential for Hepatitis B virus propagation [36]. Thus, it will be important to investigate the mechanisms of regulation of miRNA levels during virus infection, which could be at the stages of transcription, maturation, and/or degradation.

The host response to viral infection represents complex orchestration of divergent pathways deigned to eliminate the virus and protect the host. Viruses impact on many aspects of the host cell’s biology and function [20]. As seen from gene ontology analysis, the top 3 gene ontology terms related to metabolic process which indicated EV71 infection have enormous effect on Hep2 cells metabolism. Intriguingly, neurological process, apoptosis, and immune response related GO terms in biological process were also enriched from the predicted targets. Considering the fact that EV71 can cause lethal encephalitis or myocarditis, both apoptosis and immune response contribute to EV71 pathogenesis. Myocarditis, for example, represents an intricate interplay between virus and patient responses, in which both direct viral injury and immunopathologic damage caused by innocent-bystander phenomena affect the disease course [37]. Further study of the functions of those underlying miRNAs related to neurological process, apoptosis, and immune response will help to elucidate the molecular mechanisms of EV71 pathogenesis.

Our study has some limitations that will need to be addressed in future studies. We did not assess the roles in infection of miRNAs whose expression were altered after infection. As microRNAs predominately function as repressors of target gene expression, we indeed found numbers of targets in an another transcriptomic research whose expression was inversely correlated with the expression of dysregulated miRNAs in the present study. Although all assays were executed when cells demonstrated similar CPE, the virus may interact with miRNA regulatory pathways differently in different cell types in which human neural SF268 cells and rhabdomyosarcoma cells were used in those studies, respectively. Even then, these findings suggest that expression modification of host miRNAs during EV71 infection could be related to a number of cellular physiological processes that eventually control the cell fate.

In summary, in this study, we identified the miRNAs involved in the host response to EV71 infection using deep sequencing technology. Our findings provide a deeper understanding of the mechanisms underlying EV71 infection. Once the role of these miRNAs in the regulation of host-EV71 interaction has been determined, it will improve the protection and treatment strategies in enterovirus infection.

Acknowledgments

The first two authors contributed equally to this study. The study was supported by the grant from Important National Science & Technology Specific Projects during the eleventh five-year plan period (2008ZX10004-002), the National Natural Science Foundation of China (30901285), the Natural Science Foundation of Jiangsu Province (SBK200922783), and the Open Research Fund of State Key Laboratory of Bioelectronics Southeast University.

Supplementary Materials

Supplementary materials contains “Gene Ontology Analysis”. Targets were predicted for all identified differentially expressed miRNA families. In total 5765 unique target genes were predicted for 64 of the differentially expressed miRNAs.

  1. Supplementary Table 1
  2. Supplementary Table 2