# Identification of microRNAs Involved in the Host Response to Enterovirus 71 Infection by a Deep Sequencing Approach

Authors: Lunbiao Cui; Xiling  Guo; Yuhua  Qi; Xian  Qi; Yiyue  Ge; Zhiyang  Shi; Tao  Wu; Jun  Shan; Yunfeng  Shan; Zheng  Zhu; Hua  Wang
Journal: Journal of Biomedicine and Biotechnology
 (2010)

---

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


---

## Body

## 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 [1–4] 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-speciﬁc 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 proﬁling 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 37°C 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/) [12–14]. 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 16°C, 30 min at 42°C and 5 min at 85°C, and then held at 4°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 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°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 2-ΔΔCT [18].


### 2.6. Statistical Analysis

We firstly used theZ 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 speciﬁc miRNA per cell as a fraction or proportion of the total number of miRNA molecules in that cell. The same proportion of speciﬁc 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 fulﬁlled three criteria: (1) mean fold change>2or <0.5, (2) having at least 10 copies by SOLiD sequencing, and (3) Z score>1.96 or <−1.96.


## 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 37°C 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/) [12–14]. 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 16°C, 30 min at 42°C and 5 min at 85°C, and then held at 4°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 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°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 2-ΔΔCT [18].


## 2.6. Statistical Analysis

We firstly used theZ 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 speciﬁc miRNA per cell as a fraction or proportion of the total number of miRNA molecules in that cell. The same proportion of speciﬁc 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 fulﬁlled three criteria: (1) mean fold change>2or <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 (Table1). 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.Table 1
Number of reads of miRNAs from EV71 infected and non-infected Hep2 cells.
Infected cellsControl cellsHigh quality/both adapter35,27278,143Exact match to known human miRNAs2,7318,457Loose match to known human miRNAs20,38947,652The 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 (Figure1), which is within the typical size range of human miRNAs.Figure 1
Size distribution of sequenced short RNAs.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 (Tables2 and 3).Table 2
Upexpressed miRNAs in infected cells compared with non-infected cells.
miRNAReads in infectedReads in non-infectedZ scoreRatio infected/non-infected (Normalized)hsa-mir-6361105.05958439#hsa-mir-6193054123.5944353417.31099765hsa-mir-1302-83875426.4642829116.67720702hsa-mir-1302-2156022853.0515521615.92193571hsa-mir-1302-6,73875726.2815161115.79945928hsa-mir-1302-33905826.3497765615.64741958hsa-mir-1302-13665925.2411671614.4356116hsa-mir-1302-53055022.9806197814.19501807hsa-mir-1302-43005022.7356495313.96231285hsa-mir-12731623715.8090752710.18871479hsa-mir-12903487.1948481519.889971606hsa-mir-12681745.087180489.889971606hsa-mir-11781133.9592633538.532524522hsa-mir-124638212122.482476167.346561309hsa-mir-197234126.4883324846.593314404hsa-mir-1285-1963510.794916496.382771591hsa-mir-5661353.9019006126.05033557hsa-mir-12261043.3763390145.817630356hsa-mir-6351253.6450089945.584925142hsa-mir-6271674.1327980285.318976326hsa-mir-518c1363.6474438065.041946309hsa-mir-1185-220104.3738045874.654104285hsa-mir-1289-11473.6592809774.654104285hsa-mir-5391163.1160086464.266262261hsa-mir-124718103.950663744.188693856hsa-mir-7641273.1470563353.989232244hsa-mir-197933020116.155625063.820533368hsa-mir-12721172.8737054273.656796224hsa-mir-32428184.5586161363.619858888hsa-mir-197582954824.50785213.520303332hsa-mir-74428194.4142396933.429339999hsa-mir-54119163.1266306692.763374419hsa-mir-125413112.5762832222.750152532hsa-mir-42142384.3792709312.572005hsa-mir-123711102.2313814742.559757357hsa-mir-320c-174775.0811965782.236387773hsa-mir-61516172.3091015492.190166722hsa-mir-72014152.1400610262.171915333hsa-mir-1401842087.2992473532.058546126hsa-mir-4511882147.316425942.044326181hsa-mir-213445399231.356789482.008190038Table 3
Down-expressed miRNAs in infected cells compared with non-infected cells.
miRNAReads in infectedReads in non-infectedZ scoreRatio infected/non-infected (Normalized)hsa-mir-584010-2.073060hsa-mir-221470-4.62610.132974hsa-let-7e23199-6.405150.268956hsa-mir-1180543-2.968150.270587hsa-mir-1259646-2.912930.303529hsa-mir-338128905-12.42480.32913hsa-mir-30a3042083-18.60080.339618hsa-mir-19b-14872952-20.56040.383901hsa-mir-19b-24912946-20.39310.387842hsa-mir-545741-2.338450.397302hsa-mir-433635-2.154070.398923hsa-mir-58293531-8.272170.407563hsa-mir-26a-1129735-9.724250.408421hsa-mir-26a-2132746-9.736020.411757hsa-mir-27b123677-9.080340.422788hsa-mir-45231158-4.096660.456574hsa-mir-197438188-4.345080.470362hsa-mir-27a2801382-11.80870.471472hsa-mir-26b35172-4.128630.473528hsa-mir-30e86419-6.395210.477629hsa-mir-22250236-4.642150.49302hsa-mir-66048225-4.49920.496438


### 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, (P<.05), 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).Figure 2
Predicted targets gene ontology terms in the biological process category.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].Table 4
Reverse regulatory association of up-regulated miRNAs and their predicted mRNA targets.
Up-regulated miRNATarget genesFunctionsReferencehsa-mir-1289-1Interferon regulatory factor 2IFN response[20]hsa-mir-140,hsa-mir-320cv-yes-1 Yamaguchi sarcoma viral oncogene homolog 1Oncogene[20]hsa-mir-421v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog[21]hsa-mir-1226Interleukin 10 receptor, betaImmune response[20]hsa-mir-1290Nuclear receptor coactivator 2Transcriptional regulators[20]hsa-mir-320c,hsa-mir-1290Transcription factor AP-2 beta[20]hsa-mir-320cpre-B-cell leukemia homeobox 3[21]hsa-mir-636,hsa-miR-1979,hsa-mir-1302,hsa-mir-518c,hsa-mir-1226,hsa-mir-1290,hsa-mir-421,hsa-mir-21,hsa-mir-140ribosomal protein S6 kinaseKinases and phosphatases[20]hsa-mir-1289-1,hsa-mir-1272,hsa-mir-421Splicing factor, arginine/serine-rich 1RNA synthesis and modification[20]hsa-mir-1272,hsa-mir-140,hsa-mir-636Neuro-oncological ventral antigen 1Neuron specific[20]hsa-mir-539DEAD (Asp-Glu-Ala-Asp) box polypeptide 5[21]hsa-mir-1254,hsa-mir-627,hsa-mir-320c,hsa-mir-421,hsa-mir-21,hsa-mir-140cyclin-dependent kinase 6cell cycle[21]hsa-mir-619protocadherin 9cell adhesion[21]hsa-mir-627RER1 retention in endoplasmic reticulum 1 homolog (S. cerevisiae)trafficking and targeting proteins[21]hsa-mir-140sorting nexin 2[21]hsa-mir-1289-1xenotropic and polytropic retrovirus receptorCellular receptors[21]hsa-mir-21interleukin 12A (natural killer cell stimulatory factor 1, cytotoxic lymphocyte maturation factor 1, p35)cell signalling, extracellular communication[21]hsa-mir-1178,hsa-mir-140transforming growth factor, alpha[21]hsa-mir-1178dual specificity phosphatase 6intracellular trasducers/effector/modulators[21]Table 5
Reverse regulatory association of down-regulated miRNAs and their predicted mRNA targets.
Down-regulated miRNATarget genesFunctionsReferencehsa-mir-27Cytochrome P450,superfamily I polypeptide 1Mitochondrial function[20]hsa-mir-26a,hsa-mir-27BCL2-antagonist/killer 1Apoptosis[20]hsa-mir-26aProgrammed cell death 10[20]hsa-mir-582Amyloid beta (A4) precursor protein[20]hsa-mir-19BCL2-associated athanogene 5[20]hsa-mir-1180,hsa-mir-19,hsa-mir-582,hsa-mir-26a,hsa-mir-30a,hsa-let-7e,hsa-mir-545calcium channel, voltage-dependent, L type, alpha 1C subunitMembrane transporters[20]hsa-mir-433,hsa-let-7eCathepsin CProtein degradation[20]hsa-mir-26a,hsa-mir-545Muscle RAS oncogene homologOncogene[20]hsa-mir-19RAB5B, member RAS oncogene family[20]hsa-mir-19neurotrophic tyrosine kinase, receptor, type 2Kinases and phosphatases[20]hsa-mir-19Fas-activated serine/threonine kinase[20]hsa-mir-582,hsa-let-7eadaptor-related protein complex 1, sigma 1 subunittrafficking and targeting proteins[21]hsa-mir-27actin, alpha 2, smooth muscle, aortacytoskeleton/motility[21]


## 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 (Table1). 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.Table 1
Number of reads of miRNAs from EV71 infected and non-infected Hep2 cells.
Infected cellsControl cellsHigh quality/both adapter35,27278,143Exact match to known human miRNAs2,7318,457Loose match to known human miRNAs20,38947,652The 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 (Figure1), which is within the typical size range of human miRNAs.Figure 1
Size distribution of sequenced short RNAs.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 (Tables2 and 3).Table 2
Upexpressed miRNAs in infected cells compared with non-infected cells.
miRNAReads in infectedReads in non-infectedZ scoreRatio infected/non-infected (Normalized)hsa-mir-6361105.05958439#hsa-mir-6193054123.5944353417.31099765hsa-mir-1302-83875426.4642829116.67720702hsa-mir-1302-2156022853.0515521615.92193571hsa-mir-1302-6,73875726.2815161115.79945928hsa-mir-1302-33905826.3497765615.64741958hsa-mir-1302-13665925.2411671614.4356116hsa-mir-1302-53055022.9806197814.19501807hsa-mir-1302-43005022.7356495313.96231285hsa-mir-12731623715.8090752710.18871479hsa-mir-12903487.1948481519.889971606hsa-mir-12681745.087180489.889971606hsa-mir-11781133.9592633538.532524522hsa-mir-124638212122.482476167.346561309hsa-mir-197234126.4883324846.593314404hsa-mir-1285-1963510.794916496.382771591hsa-mir-5661353.9019006126.05033557hsa-mir-12261043.3763390145.817630356hsa-mir-6351253.6450089945.584925142hsa-mir-6271674.1327980285.318976326hsa-mir-518c1363.6474438065.041946309hsa-mir-1185-220104.3738045874.654104285hsa-mir-1289-11473.6592809774.654104285hsa-mir-5391163.1160086464.266262261hsa-mir-124718103.950663744.188693856hsa-mir-7641273.1470563353.989232244hsa-mir-197933020116.155625063.820533368hsa-mir-12721172.8737054273.656796224hsa-mir-32428184.5586161363.619858888hsa-mir-197582954824.50785213.520303332hsa-mir-74428194.4142396933.429339999hsa-mir-54119163.1266306692.763374419hsa-mir-125413112.5762832222.750152532hsa-mir-42142384.3792709312.572005hsa-mir-123711102.2313814742.559757357hsa-mir-320c-174775.0811965782.236387773hsa-mir-61516172.3091015492.190166722hsa-mir-72014152.1400610262.171915333hsa-mir-1401842087.2992473532.058546126hsa-mir-4511882147.316425942.044326181hsa-mir-213445399231.356789482.008190038Table 3
Down-expressed miRNAs in infected cells compared with non-infected cells.
miRNAReads in infectedReads in non-infectedZ scoreRatio infected/non-infected (Normalized)hsa-mir-584010-2.073060hsa-mir-221470-4.62610.132974hsa-let-7e23199-6.405150.268956hsa-mir-1180543-2.968150.270587hsa-mir-1259646-2.912930.303529hsa-mir-338128905-12.42480.32913hsa-mir-30a3042083-18.60080.339618hsa-mir-19b-14872952-20.56040.383901hsa-mir-19b-24912946-20.39310.387842hsa-mir-545741-2.338450.397302hsa-mir-433635-2.154070.398923hsa-mir-58293531-8.272170.407563hsa-mir-26a-1129735-9.724250.408421hsa-mir-26a-2132746-9.736020.411757hsa-mir-27b123677-9.080340.422788hsa-mir-45231158-4.096660.456574hsa-mir-197438188-4.345080.470362hsa-mir-27a2801382-11.80870.471472hsa-mir-26b35172-4.128630.473528hsa-mir-30e86419-6.395210.477629hsa-mir-22250236-4.642150.49302hsa-mir-66048225-4.49920.496438


## 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, (P<.05), 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).Figure 2
Predicted targets gene ontology terms in the biological process category.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].Table 4
Reverse regulatory association of up-regulated miRNAs and their predicted mRNA targets.
Up-regulated miRNATarget genesFunctionsReferencehsa-mir-1289-1Interferon regulatory factor 2IFN response[20]hsa-mir-140,hsa-mir-320cv-yes-1 Yamaguchi sarcoma viral oncogene homolog 1Oncogene[20]hsa-mir-421v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog[21]hsa-mir-1226Interleukin 10 receptor, betaImmune response[20]hsa-mir-1290Nuclear receptor coactivator 2Transcriptional regulators[20]hsa-mir-320c,hsa-mir-1290Transcription factor AP-2 beta[20]hsa-mir-320cpre-B-cell leukemia homeobox 3[21]hsa-mir-636,hsa-miR-1979,hsa-mir-1302,hsa-mir-518c,hsa-mir-1226,hsa-mir-1290,hsa-mir-421,hsa-mir-21,hsa-mir-140ribosomal protein S6 kinaseKinases and phosphatases[20]hsa-mir-1289-1,hsa-mir-1272,hsa-mir-421Splicing factor, arginine/serine-rich 1RNA synthesis and modification[20]hsa-mir-1272,hsa-mir-140,hsa-mir-636Neuro-oncological ventral antigen 1Neuron specific[20]hsa-mir-539DEAD (Asp-Glu-Ala-Asp) box polypeptide 5[21]hsa-mir-1254,hsa-mir-627,hsa-mir-320c,hsa-mir-421,hsa-mir-21,hsa-mir-140cyclin-dependent kinase 6cell cycle[21]hsa-mir-619protocadherin 9cell adhesion[21]hsa-mir-627RER1 retention in endoplasmic reticulum 1 homolog (S. cerevisiae)trafficking and targeting proteins[21]hsa-mir-140sorting nexin 2[21]hsa-mir-1289-1xenotropic and polytropic retrovirus receptorCellular receptors[21]hsa-mir-21interleukin 12A (natural killer cell stimulatory factor 1, cytotoxic lymphocyte maturation factor 1, p35)cell signalling, extracellular communication[21]hsa-mir-1178,hsa-mir-140transforming growth factor, alpha[21]hsa-mir-1178dual specificity phosphatase 6intracellular trasducers/effector/modulators[21]Table 5
Reverse regulatory association of down-regulated miRNAs and their predicted mRNA targets.
Down-regulated miRNATarget genesFunctionsReferencehsa-mir-27Cytochrome P450,superfamily I polypeptide 1Mitochondrial function[20]hsa-mir-26a,hsa-mir-27BCL2-antagonist/killer 1Apoptosis[20]hsa-mir-26aProgrammed cell death 10[20]hsa-mir-582Amyloid beta (A4) precursor protein[20]hsa-mir-19BCL2-associated athanogene 5[20]hsa-mir-1180,hsa-mir-19,hsa-mir-582,hsa-mir-26a,hsa-mir-30a,hsa-let-7e,hsa-mir-545calcium channel, voltage-dependent, L type, alpha 1C subunitMembrane transporters[20]hsa-mir-433,hsa-let-7eCathepsin CProtein degradation[20]hsa-mir-26a,hsa-mir-545Muscle RAS oncogene homologOncogene[20]hsa-mir-19RAB5B, member RAS oncogene family[20]hsa-mir-19neurotrophic tyrosine kinase, receptor, type 2Kinases and phosphatases[20]hsa-mir-19Fas-activated serine/threonine kinase[20]hsa-mir-582,hsa-let-7eadaptor-related protein complex 1, sigma 1 subunittrafficking and targeting proteins[21]hsa-mir-27actin, alpha 2, smooth muscle, aortacytoskeleton/motility[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 [27–30]. 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 proﬁling 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 proﬁling 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.



---

*Source: 425939-2010-06-10.xml*