High-Throughput Sequencing of MicroRNAs in Adenovirus Type 3 Infected Human Laryngeal Epithelial Cells

Adenovirus infection can cause various illnesses depending on the infecting serotype, such as gastroenteritis, conjunctivitis, cystitis, and rash illness, but the infection mechanism is still unknown. MicroRNAs (miRNA) have been reported to play essential roles in cell proliferation, cell differentiation, and pathogenesis of human diseases including viral infections. We analyzed the miRNA expression profiles from adenovirus type 3 (AD3) infected Human laryngeal epithelial (Hep2) cells using a SOLiD deep sequencing. 492 precursor miRNAs were identified in the AD3 infected Hep2 cells, and 540 precursor miRNAs were identified in the control. A total of 44 miRNAs demonstrated high expression and 36 miRNAs showed lower expression in the AD3 infected cells than control. The biogenesis of miRNAs has been analyzed, and some of the SOLiD results were confirmed by Quantitative PCR analysis. The present studies may provide a useful clue for the biological function research into AD3 infection.


Introduction
Adenovirus (AD) belongs to a family of nonenveloped icosahedral viruses containing a double-stranded DNA genome, Adenoviridae, and genus Mastadenovirus. There are 53 human AD serotypes divided into seven subgroups (A to G) based on immunologic, biologic, and biochemical characteristics [1]. Among these, viruses in groups B1, C, and E cause respiratory infections, group B2 viruses infect the kidney and urinary tract; group F viruses cause gastroenteritis; and several group D serotypes are associated with epidemic keratoconjunctivitis [1]. Most people recover from AD infections by themselves, but people with immunodeficiency and small children are at higher risk of fatal adenovirus infection. Studies showed that AD can evade host immune responses, such as inhibition of interferon functions by RNA and E1A, and inhibition of intrinsic cellular apoptosis in infected cells [2]. The interaction between the virus and its host is still not clear. Recently, several studies evaluated human adenoviruses which play an important role in the study of tumorigenesis, including oncogenic transformation of cells, immune evasion, angiogenesis, and metastasis [3,4]. There are no antiviral drugs to treat adenoviral infections, so treatment is largely directed at the symptoms. The discovery of the AD yielded further evidence that viruses can sometimes interfere with each other during the growth and replication within a host animal. Other researchers found that this interference can be mediated by interferon produced by the host animal. So interferon has become prominent in the treatment of a variety of cancers and infectious diseases, such as hepatitis C [5]. Recently, several studies demonstrated that RNA interference (RNAi) can also be used for treatment with viral infection in vitro, such as hepatitis B virus, coxsackievirus B3, and Coxsackievirus A16 [6][7][8]. These findings indicate that those specific siRNAs targeted viruses are not as effective as expected.
RNAi can also be triggered by microRNAs (miRNAs), which are short ∼22 nucleotide RNA sequences that bind complementary sequences in the 3 UTR of multiple target Journal of Biomedicine and Biotechnology mRNAs, playing an important role in the development, proliferation, differentiation and apoptosis of organisms. Growing evidences have indicated a strong association between microRNAs and various viral infections. For example, miR-125b, mirR-150, miR-382, and miR-223 have been reported to contribute to the maintenance of human immunodeficiency virus-1 (HIV-1) latency in resting CD4 + T cells [9]. Epstein-Barr virus and Kaposi sarcoma herpes viruses (KSHVs) have been reported to encode miRNAs [10], but the functions of most of them are not known. miRNAs have been reported to be used as a tool for gene specific therapeutics against viral infections [11]. The potential use of miRNAs as novel noninvasive biomarkers for diagnosis has been suggested in other studies [12][13][14].
In this study, we used a strategy of initial screening by Applied Biosystems' next generation sequencing system which is Sequencing by Oligonucleotide Ligation and Detection (SOLiD) and validation by quantitative RT-PCR (qRT-PCR), to analyze the different miRNA expression profiling in AD3 infected Hep2 cells.

Cell Culture and Virus Infection.
Human laryngeal epithelial (Hep2) cells were cultivated in complete RPMI medium 1640 (Invitrogen) supplemented with 10% fetal calf serum (FCS) at 37 • C in an atmosphere of 5% CO 2 . When the Hep2 cells were grown in 25-cm 2 flasks reached to 70% confluence, they were infected with AD3 at a multiplicity of infection (m.o.i.) of 0.4 under 50% tissue culture infectious doses (TCID 50 ), and maintained after infection at 37 • C in RPMI 1640 medium with 2% FCS.

RNA Isolation.
Hep2 cells were infected with AD3 as described above. At 6 hours, 24 hours, 48 hours, and 72 hours post infection (p.i.), RNA samples were prepared from AD3 infected Hep2 cells and controlled Hep2 cells using a mirVana miRNA Isolation Kit (Ambion): 10 2 -10 7 cultured cells were washed in 1% PBS for three times, washed cells were mixed with 600 μl Lysis/Binding Solution and 1/10 volume of miRNA Homogenate Additive thoroughly by vortexing and incubated on ice for 10 minutes An equal volume of acid/phenol/chloroform (Ambion) was then added to the each aliquot. The solution was centrifuged for 5 minutes at 10,000 g at room temperature. For miRNA isolation, the extraction was mixed thoroughly with 1/3 volume of 100% ethanol and passed through a column (Ambion). The flow-through was contained the small RNAs. The filtrate was then mixed thoroughly with 2/3 volume 100% ethanol and passed through the second column (Ambion). The column was washed following the protocol, and miRNAs were eluted in 100 μl of elution buffer (95 • C). For total RNA isolation, the extraction was mixed thoroughly with 1.25 volume of 100% ethanol and passed through a column (Ambion). The column was washed following the protocol and total RNA were eluted in 100 μl of elution buffer (95 • C). A total of 1 μl of the eluate was quantified by Nanodrop.

SOLiD Sequencing and Sequence Analysis.
MiRNA samples (100 ng) isolated from AD3 infected Hep2 cells and control Hep2 cells were processed into sequencing libraries using the Small RNA Expression Kit (Applied Biosystems). Briefly, RNA was ligated overnight with the adapters from the kit, reverse transcribed, RNAse H-treated, and PCRamplified before size selection on Agarose gels to contain 16-61 nt of inserted sequences. Libraries were amplified onto beads using emulsion PCR, deposited on slides, and sequenced using the SOLiD v 2 sequencing system (Applied Biosystems) [15] at the State Key Lab of Bioelectronics Laboratory, Southeast University of China. Finally, SOLiD data was first analyzed by SOLiD System Small RNA Analysis Pipeline Tool (RNA2MAP). Firstly, the parameters set for alignments were the maximum length (18 nt) and the tolerance of 3 mismatches when generating initial seeds locations. In the extension step, less than 6 mismatches were allowed in full length mapping. The miRBase sequences (Sanger) of human being was download from miRBase (http://www.mirbase.org/). Acceptable sequences were compared to the miRBase database (release 14.0). To increase signal above noise, we conservatively selected only those alignments corresponding to beads sampled a minimum of 10 times in any of the libraries.

Quantitative RT-PCR (qRT-PCR).
To confirm the expression of miRNAs by deep sequencing approach, stemloop quantitative RT-PCR (qRT-PCR) was performed. Briefly, cDNA was synthesized from total RNA by using AMV reverse transcriptase (TaKaRa). The 20 μL reactions were incubated for 15 minutes at 16 • C, 30 minutes at 42 • C, and 5 minutes at 85 • C, and then held at 4 • C. Subsequently, 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 minutes, followed by 40 cycles of 95 • C for 15 seconds, and 60 • C for 1 minute. Q-PCR assay was performed in triplicate using an ABI 7500 machine (Applied Biosystems). The expression levels of candidate miRNAs were measured in terms of threshold cycle value (C T ) and normalized to Human small nuclear RNA U6. The ratio of two groups of miRNAs was calculated by using the equation

Statistical Analysis.
We firstly used the Z-test to determine the statistical significance of the differences between the two libraries. 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 (p) of specic tags should be present in the miRNA library of all sequenced tags [17] Journal of Biomedicine and Biotechnology 3 So the Z-value for a specific miRNA could be calculated as following: The proportion p 0 , the expected proportion when the null hypothesis is true, is calculated as p 0 = (n test + n control )/(N test + N control ). The Z statistic is approximately normally distributed and can be compared with critical Z value for the two-sided significance level α [17]. A Z-test was used to deal with the raw data by considering both the copies and fold change of miRNAs. Candidate miRNAs should be full three criteria: (1) mean fold change > 2, (2) having at least 10 copies by SOLiD sequencing, and (3) Z score > 1.96 or < −1.96 for P < .05.

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

SOLiD Sequencing of miRNAs from Hep2 Cells Infected by AD3 and Control Hep2 Cells.
To search for the miRNAs expression in AD3 infected Hep2 cells and the control, we prepared these two miRNA samples for sequencing by the SOLiD system. 358,297 effective reads were obtained in AD3 infected cells and 479,414 effective reads in control cells. After filtrating reads contaminated by rRNAetc (rRNA, tRNA, snRNA, and snoRNA), repeat fractions, and degraded fractions, the remaining effective reads were mapped to the human precursor libraries. 37,800 and 78,143 reads were obtained from infected and control cells, respectively. In these reads, the most abundant size class was 22 nt in both of AD3 infected Hep2 cells and the control. The percentage was 29.3% and 34.6% respectively. Following was 23 nt (27.7% in AD3 infected Hep2 cells and 26.3% in the control) ( Figure 1). Compared to the miRBase (14.0), 492 precursor, miRNA out of 734 known precursor miRNAs were identified among the AD3 infected Hep2 cells, and 540 precursor miRNAs were identified among the control Hep2 cells. All of the known miRNAs have been characterized in the Sanger miRBase sequence database. Figure 2 showed the distribution of genomic loci in the chromosome. All of the human chromosome were located in the precursor miRNAs. Chromosomes 19 and X contain more precursor miRNAs (12%, and 11%, resp.). 314 miRNAs (51%) were located in introns of protein-coding genes, 231 miRNAs (38%) were located in intergenic regions, while only 43 miRNAs (7%) were located in exons of noncoding RNAs or UTR of protein coding genes. Interestedly, 23 miRNAs (4%) were located in either an exon or an intron depending on alternative splicing of the host transcript (Table 1).

Expression Profiles of miRNAs in AD3 Infected
Cells. To study the differential expression profile of known miRNAs     between AD3 infected Hep2 cells and the control, both the Z-test and fold-change analysis were made based on the sequencing results. There are 80 precursor miRNAs to meet the three criteria described above, and could be selected as candidates. Among these candidates, 44 precursor miRNAs were shown upregulated expression in AD3 infected cells (   (Table 3). Among this there are 11 miRNAs that have the higher Z values (over 10). Table 3 showed the ratio of AD3 infected cells/control cells (expression level for precursor miRNAs). Among this miR-181b-2 and miR-1274b showed lowest change in AD3 infected cells (0.19 fold) compared with control cells.

Q-PCR Detection for miRNAs Profiles in AD3 Infected
Cells. To confirm the differential expression of miRNAs in AD3 infected Hep2 cells, we performed some miRNAs tests with stem-loop Q-PCR assays in total RNA isolated from AD3 infected Hep2 cells and the control, respectively. The expression levels of candidate miRNAs were measured in terms of threshold cycle value (C T ) and normalized to Human small nuclear RNA U6. The ratio of miRNAs in two groups was calculated by using the equation 2 −ΔΔCT . The results showed that there was general consistency between Q-PCR assay and SOLiD sequencing analysis. The seven upregulated miRNAs in SOLiD results showed increased ratios measured in Q-PCR assays (Figure 3). For example, miR-1974, miR-1975, and miR-7 showed increasing values of 2.51, 4.87 and 2.44 respectively measured in Q-PCR assay, and increasing values of 2.95, 7.47 and 6.34 were obtained in SOLiD sequencing ( Table 2). The seven downregulated  Table 4 exhibited increased expression in AD3 infected cells at 6 hours post infection. However, at 48 hours p.i. and 72 hours p.i., miR-27b, miR-125b, and miR-181b showed downregulated expression in AD3 infected cells comparing controls (

Discussion
As second-generation sequencing platforms have matured, all three major high-throughput sequencing systems-Illumina's Genome Analyzer (Solexa), Life Technologies' ABI SOLiD, and Roche's 454 GS FLX have been used to study miRNA expression profiles associated with disease such as cancer, viral and metabolic diseases. Compared with the microarray technology, deep sequencing can generate millions of small RNA sequence reads from a given sample, giving possible to discover the novel miRNAs or study the function of virus encoded miRNAs by using bioinformatics and molecular biology tools. miRNAs have recently been recognized as major regulators of gene expression in various viral infections. Virus-encoded miRNAs seem to evolve rapidly and regulate both the viral life cycle and the interaction between viruses and their hosts. In the present study, we have screened out some miRNA candidates that have obvious differential expression from AD3 infected Hep2 cells by SOLiD deep sequencing system. This study identified 492 precursor miRNAs in the AD3 infected cells and more precursor miRNAs in control cells by mapping the human miRBase database. Out of these almost 38% were located in intergenic regions, which usually contain their own miRNA gene promoter and regulatory units [18]. However, as much as 51% miRNAs were located in introns, and only 7% were located in exons of noncoding RNAs or UTR of protein coding genes. These usually show a concurrent transcription and regulation expression profile originating from a common promoter with their host genes [19][20][21]. Some pre-miRNAs existing in drosphila and c elegants have been reported to be spliced from the introns in which they reside without having to undergo the microprocessor machinery. These miRNAs are called mirtrons [22,23]. Since our study identified 51% miRNAs were located in introns, we think more mirtons could be discovered through bioinformatics and molecular biology tools.
The interaction between the AD and its host is still unknown. One of the mechanisms is that the virus must interfere with the host cell antiviral defense mechanisms to maximize escape and spread of the progeny virus during a virus infection [24]. Previous studies determine the changes in the host cell gene expression profiles upon infection with AD using DNA microarray technique [24,25]. For example, many transcriptions factor E2F-dependent host cell genes showed different expression during AD infection. 45% of the upregulated genes and 25% of the downregulated genes contained potential E2F binding sites [24]. It is showen that E2F-dependent host cell genes are subjected to a selective regulation. However, the mechanism of AD infection after transcription is still unknown. miRNA can regulate gene expression through translational repression or mRNA cleavage. Previous research has reported that the expression of E2F1 is negatively regulated by two c-mycregulated miRNAs, miR-17-5p, and miR-20a [26]. In our study, upregulated miR-17 may target with transcription factor E2F1 [26], playing important role in cell differentiation, proliferation, and apoptosis. Out of these, miR-17 has been reported to be overexpression in B-cell lymphoma samples [27]. The high miR-17 cluster miRNAs suppress cell death [27]. This suggests that during a virus infection, cellular apoptosis must delay long enough and the virus needs to establish optimal conditions for replication to ensure efficient production of progeny virus. Our study may provide a useful clue for the biological function research into AD3 infected host cells. Further investigation needs to clarify the roles of identified miRNAs in the pathogenesis of AD3 infected host cells.

Conclusion
SOLiD sequencing provides a useful method for identification of the miRNAs profiles in AD3 infected Hep2 cells. 492 precursor miRNAs were identified in the AD3 infected Hep2 cells, and 540 precursor miRNAs were identified in the control. A total of 44 miRNAs demonstrated high expression and 36 miRNAs showed lower expression in the AD3 infected cells than the control. Further studies are required to identify the miRNA target genes and the functions of the miRNAs in the complex molecular network regulation during the virus infection host cells using bioinformatics tools.