Functional similarity analysis of human virus-encoded miRNAs
Journal of Clinical Bioinformatics volume 1, Article number: 15 (2011)
miRNAs are a class of small RNAs that regulate gene expression via RNA silencing machinery. Some viruses also encode miRNAs, contributing to the complex virus-host interactions. A better understanding of viral miRNA functions would be useful in designing new preventive strategies for treating diseases induced by viruses. To meet the challenge for how viruses module host gene expression by their encoded miRNAs, we measured the functional similarities among human viral miRNAs by using a method we reported previously. Higher order functions regulated by viral miRNAs were also identified by KEGG pathway analysis on their targets. Our study demonstrated the biological processes involved in virus-host interactions via viral miRNAs. Phylogenetic analysis suggested that viral miRNAs have distinct evolution rates compared with their corresponding genome.
miRNAs, about 22 nucleotides in length, constitute a large family of non-coding RNAs that regulate gene expression posttranscriptionally, leading their target mRNAs to direct destructive cleavage or translational repression by base pairing with the 3' untranslated regions (3' UTRs). miRNA-mediated regulation plays crucial roles in a wide spectrum of biological processes, including proliferation , apoptosis , development , immune system regulation , and oncogenesis .
Recent discoveries on viral miRNAs, mostly in herpesvirus family , threw lights on a new level of cross-talk between virus and host in viral infections and pathogenesis . Viral miRNAs have been reported to participate in immune evasion by directly down-regulating host immune defence genes, and even to cooperate with viral proteins to target the same process . The combination of protein-mediated and miRNA-mediated regulations forms an intricate strategy for viruses to resist host defence system and thus increase the opportunities of their survival.
The research on viral miRNAs is still far from exhausted, with many unknown miRNA functions yet to be discovered. miRNA identification using computational tools is the most widely used method. In contrast to most eukaryotic miRNAs, virus-encoded miRNAs do not have homologs in other viral genomes or in the genome of the human host , and thus are difficult to be identified using existing miRNA gene prediction tools. Cloning and sequencing small RNA libraries to identify and characterize miRNAs is the basic method for miRNA discovery, since computationally predicted miRNAs should also be confirmed by experimental methods. Reverse ligation-mediated RT-PCR  is a widely used method in the identification of mature miRNAs and has been used to detect maturely processed MuHV-4 miRNAs . Experimental validation is still a barrier in miRNA identification, especially in host cells infected by viruses. Currently, only a small fraction of viral miRNAs has been identified, and the functions of most of these viral miRNAs remain unknown. To bridge the gap in understanding the targets regulated by these virus-encoded miRNAs, we used computational method to predict host targets of viral miRNAs and measured their functional similarities to reveal the interspecies cross-talk between virus and host by viral miRNAs.
Materials and methods
Host target gene prediction of viral miRNAs
In order to determine how viruses reshape the physiological states of human cells by their encoded miRNAs, we first predicted host genes targeted by viral miRNAs. We collected viral miRNAs encoded by BK polyomavirus (BKV), Epstein-Barr virus (EBV), human cytomegalovirus (HCMV), human immunodeficiency virus 1 (HIV1), human herpesvirus 1 (HSV1), human herpesvirus 2 (HSV2), and Kaposi's sarcoma-associated herpesvirus (KSHV). Viral miRNA sequences were retrieved from miRBase  release 16 (Sep 2010). We extracted 3' UTR sequences in a single FASTA format file from human genome (version 18) that was downloaded from UCSC .
Host target genes of virus-encoded miRNAs were predicted by the algorithm of Probability of Interaction by Target Accessibility (PITA) that computes the difference between the free energy gained from the formation of the miRNA-mRNA duplex and the energetic cost of unpairing the mRNA to make it accessible to the miRNA . We chose PITA for viral miRNA target prediction because it had been demonstrated to reach high accuracy, and more importantly, it takes advantage of the target accessibility but not conservation information to reduce false positive. Conservation information, which was used by most of other methods, is not suitable for predicting target genes of the less evolutionarily conserved viral miRNAs . We used a flank of 3 upstream and 15 downstream nucleotides when performing prediction, since miRNA-mRNA interaction requires unpairing of bases flanking the targets. To reduce false positive, the prediction results were narrowed down by the criteria of 7-8 bases seed length, with no G:U wobble or loops, no mismatch, and ∆∆G < -20 kcal/mol.
Functional similarity measurement of viral miRNAs
We have previously proposed a method for systematic study of functional similarities among miRNAs by using their target genes Gene Ontology (GO) semantic similarities . As described in our previous study, the functional similarity of human miRNAs, obtained by our method, showed positive correlation with expression similarity, and the clustering results derived from the functional similarity were coherent with biological knowledge in many aspects including disease association, genome conservation, and the cross-talk between hosts and viruses . The method is reliable to calculate functional similarities and sensible to cluster miRNAs, and thus can be used to predict novel miRNA functions.
Here, we applied our method to measure functional similarities among viral miRNAs. As suggested in our previous study , the measurement was fundamentally based on host target genes of viral miRNAs. Biological process ontology was used to annotate target genes, and Wang's method  was used to calculate semantic similarity. Semantic similarity calculation was implemented by our in-house developed R package GOSemSim .
Similarity scores were then analyzed by R package pvclust , which used multi-scale bootstrap re-sampling to evaluate the uncertainty of cluster analysis. The agglomerative method, average linkage, was used, and 10,000 bootstrap replications were run. All clusters were extracted with approximately unbiased (AU) p-value > 0.95, meaning that the hypothesis with "the cluster does not exist" is rejected with significance level of 0.05.
GO enrichment analysis of significant clusters
The common biological processes regulated by these significant miRNA clusters were evaluated by GOstats  with p < 0.001. GOstats using hypergeometric model to assess whether the number of selected genes associated with the GO term is larger than expected. This method had been used to predict the functions of miRNAs  and can be used to provide biological insights of viral strategies.
KEGG enrichment analysis of genes targeted by viral miRNAs
In order to uncover higher order functions of how viruses transform cellular states by their encoded miRNAs, we adopted KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis to identify pathways regulated by viral miRNAs to provide biological insights. KEGG pathway is a collection of manually drawn pathway maps representing molecular interactions and reaction networks, and has been widely used for biological interpretation of higher level systemic functions . KEGG enrichment analysis is calculated by R package SubpathwayMiner , which implements hypergeometric test to measure p-value for evaluating enrichment significance of pathways. SubpathwayMiner also provides the FDR-corrected q-values to reduce the false positive discovery rate .
Comparing viral miRNA regulated pathways
Significant KEGG pathways regulated by different viruses were compared and visualized using our in-house developed R package clusterProfiler http://bioconductor.org/packages/2.8/bioc/html/clusterProfiler.html. ClusterProfiler, which was implemented based on R and its plotting system ggplot2 , is released under the Artistic-2.0 license within Bioconductor project . ClusterProfiler was designed to provide statistical analysis of GO and KEGG and visualization tools for comparing functional profiles among gene clusters. More details on the use of clusterProfiler are available in the package vignette.
We built phylogenetic trees of human viruses based on the functions their miRNAs encoded. Phylogenetic trees were constructed by R package phangorn  using the popular neighbour-joining (NJ) method. For validating our phylogenetic analysis, we compared our results with phylogenetic trees obtained from whole genome sequence alignment. Complete genome sequences of viruses were obtained from NCBI nucleotide database. Multiple sequence alignment and phylogenetic tree construction were done by ClustalX (version 2.0.12)  using NJ algorithm. Robinson-Foulds (RF) metric , the most widely used method in comparing phylogenetic trees, was adopted to compute the topological distance between phylogenetic trees. RF rate was obtained by normalizing the RF distances by the number of total edges for representing the relationship among trees . RF rate measures the dissimilarity between two trees.
Results and Discussion
The functional similarity matrix of the pairwise viral miRNAs was then analyzed by R package pvclust to assess the uncertainty of clustering result . We obtained 3 clusters with AU p-value > 0.95. These 3 clusters contain 2 (ebv-miR-BART20-5p and hsv2-miR-H6*), 7 (hcmv-miR-UL70-3p, hiv1-miR-H1, hsv1-miR-H1, hsv1-miR-H6-5p, hsv2-miR-H10, hsv2-miR-H22 and kshv-miR-K12-12), and 3 (ebv-miR-BART17-5p, hcmv-miR-UL148D and hsv1-miR-H6-3p) miRNAs as illustrated in Figure 2.
GO enrichment analysis was performed across these three significant clusters to discover their biological themes. As a result, Cluster 1 suggests the down-regulation of xylosyltransferase activity, involved in O-glycan processing. O-glycans had been described to play roles in cell polarity , which involves in the formation of immunological synapse , indicating that viruses prevent the formation of immunological synapse by inhibiting the xylosyltransferase activity. Cluster 2 represses a wide range of binding activities, including protein binding, DNA binding, receptor binding, and enzyme binding. Especially, the inhibition of MHC protein binding and CD40 receptor binding suggests that viruses use miRNAs to interfere the activation of antigen presenting cells. This may be the strategy for viruses to extend the life of the infected cells and to establish a favourable environment for their replication. Cluster 3 down-regulates transcription factor activity to favour viral latency. EBV BART miRNAs were expressed in latent infection . Hsv1-miR-H6-3p had been reported to promote latency by inhibiting the expression of HSV-1-encoded transcription factor, ICP4, that is required for the expression of most HSV-1 genes during productive infection [30, 31]. It has been reported that viruses encode proteins to interfere with transcription factors, and that miRNAs are more versatile to reshape the cellular status to escape host immune system and to hijack cellular machinery for their replication [32, 33].
The average similarity among 29 viral miRNAs is only 0.434 and most of the miRNAs cannot be clustered with AU p-value > 0.95, indicating that a majority of these viral miRNAs have distinct functions, with the versatilities and flexibilities of viral regulations.
Viral infection generally results in dramatic alterations in cellular mRNA expression. We thus further identified cellular pathways perturbed by viral miRNAs using KEGG enrichment analysis to gain a higher level perception. The statistically and significantly enriched pathways perturbed by different viruses were then compared and illustrated in Figure 3.
As shown in Figure 3, different viruses have distinct strategies to reshape cellular status. It seems that viral miRNAs were designed to against many important pathways to favour their pathogenesis. KSHV-encoded miRNAs had been described to directly down regulate a major regulator of cell adhesion, THBS1 , that is involved in the recruitment of monocytes and T cells to the sites of infection . Down regulation of THBS1 by KSHV miRNAs may aid KSHV-infected cells in avoiding detection by the host immune system . HIV1-encoded miRNAs play critical roles in oncogenic transformation , and three miRNAs encoded by EBV are crucial for efficient B cell transformation . These biological findings are consistent with our analyses. In addition, many pathways in our analyses have not been reported yet, and thus can serve as putative functions played by viral miRNAs for further investigations.
Reconstructing the tree of virus phylogeny is still the cardinal challenges in biology. Here we used the similarity index by functions that viral miRNAs encoded to rebuild the phylogenetic tree. We then compared our tree with phylogenetic tree obtained by genome alignment as shown in Figure 4. Although the tree based on genome alignment included biases like horizontal gene transfer (HGT) , genome alignment is still the de facto standard for phylogenetic tree construction.
We evaluated the similarity between these two trees. The topological distance between them was calculated by RF metric to be 8, and the corresponding RF rate is 0.727, and thus the similarity between the two trees is 0.273. Surprisingly, viral miRNAs have distinct evolution rates compared to their corresponding genome based on our functional analysis. We thus measured the evolutionary distance among viruses by their encoded miRNA sequences. RF distance between phylogenetic trees obtained from genome sequences and miRNA sequences is 6, and the corresponding RF rate is 0.545, and thus the similarity between the two trees is 0.455.
Viral miRNAs have different properties compared with viral proteins, such as small and non-immunogenic, and thus they may serve as ideal tools to interpolate cellular environment in the ways that benefit virus replication. This would mean an evolutionary reward for rapid adjustment to the host and environmental statuses. Viral miRNAs do not share a high level of homology even within the members of the same family [6, 39]. Phylogenetic analysis of all previously known virus miRNA genes showed that most of the known viral miRNAs have long distant relationships and could be classified into specific miRNA families . These findings are consistent with our phylogenetic analysis, suggesting that viral miRNAs may evolve more rapidly than their genome. Especially, the functions of viral miRNAs evolve even more rapidly than their sequences.
Obviously, miRNAs are ideal for the tight space constraints characteristic of viral genomes and the evolution of a miRNA down-regulating a new target gene can presumably be achieved more easily than the evolution of a new protein . It must be pointed out, however, our current method only provides a perception of viral miRNA perspective and may contain some biases, as it did not consider the fact that the activation of viral miRNAs depends on the viral life cycle in various latent or at lytic stages, and the specific infected cell types.
Intimately connected with various kinds of diseases, viruses pose a crucial health problem on host. Host cellular expression profiles altered by virus-encoded miRNAs form a new regulatory layer. Though studies into pathogenesis by viral miRNAs are still in its infancy, the interspecies regulation at the miRNA level fuels the spark of the investigation into the repertoire of virus-host interactions. Here, we applied our method to assess the functional similarity among viral miRNAs. Our analyses showed that viral miRNAs have diverse functions. We then summarized cellular pathways regulated by viral miRNAs by the GO and KEGG enrichment analyses. Phylogenetic trees were reconstructed to reveal the evolutionary distance at the perspective of viral miRNAs.
Experimental validation of computational results is still a challenge, a hindrance towards understanding the functions of viral miRNAs. We believe that the integration of bioinformatics with microarray and proteomic data would be a promising way to elucidate the whole picture of virus-host interaction mediated by viral miRNAs. In addition, the identification of roles played by viral miRNAs in pathogenesis would help in designing new preventive and therapeutic approaches. This has also been described as new therapeutics to correct the aberrant activity of miRNA-mRNA interaction by using anti-miRNA oligonucleotides (AMOs) . We hope that this work can provide a better understanding of basic biological processes involved in latency and oncogenic transformation mediated by viral miRNAs.
Manni I, Artuso S, Careccia S, Rizzo MG, Baserga R, Piaggio G, Sacchi A: The microRNA miR-92 increases proliferation of myeloid cells and by targeting p63 modulates the abundance of its isoforms. The FASEB Journal. 2009, 23: 3957-3966. 10.1096/fj.09-131847.
Wang Y, Lee CGL: MicroRNA and cancer--focus on apoptosis. J Cell Mol Med. 2009, 13: 12-23.
Wienholds E, Koudijs MJ, van Eeden FJM, Cuppen E, Plasterk RHA: The microRNA-producing enzyme Dicer1 is essential for zebrafish development. Nat Genet. 2003, 35: 217-218. 10.1038/ng1251.
Lu LF, Liston A: MicroRNA in the immune system, microRNA as an immune system. Immunology. 2009, 127: 291-298. 10.1111/j.1365-2567.2009.03092.x.
Cho W: OncomiRs: the discovery and progress of microRNAs in cancers. Molecular Cancer. 2007, 6: 60-10.1186/1476-4598-6-60.
Pfeffer S, Sewer A, Lagos-Quintana M, Sheridan R, Sander C, Grasser FA, van Dyk LF, Ho CK, Shuman S, Chien M, Russo JJ, Ju J, Randall G, Lindenbach BD, Rice CM, Simon V, Ho DD, Zavolan M, Tuschl T: Identification of microRNAs of the herpesvirus family. Nat Meth. 2005, 2: 269-276. 10.1038/nmeth746.
Veksler-Lublinsky I, Shemer-Avni Y, Kedem K, Ziv-Ukelson M: Gene bi-targeting by viral and human miRNAs. BMC Bioinformatics. 2010, 11: 249-10.1186/1471-2105-11-249.
Stern-Ginossar N, Elefant N, Zimmermann A, Wolf DG, Saleh N, Biton M, Horwitz E, Prokocimer Z, Prichard M, Hahn G, Goldman-Wohl D, Greenfield C, Yagel S, Hengel H, Altuvia Y, Margalit H, Mandelboim O: Host Immune System Gene Targeting by a Viral miRNA. Science. 2007, 317: 376-381. 10.1126/science.1140956.
Zhu JY, Strehle M, Frohn A, Kremmer E, Hofig KP, Meister G, Adler H: Identification and Analysis of Expression of Novel MicroRNAs of Murine Gammaherpesvirus 68. J Virol. 2010, 84: 10266-10275. 10.1128/JVI.01119-10.
Diebel KW, Smith AL, van Dyk LF: Mature and functional viral miRNAs transcribed from novel RNA polymerase III promoters. RNA. 2010, 16: 170-185. 10.1261/rna.1873910.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-144. 10.1093/nar/gkj112.
Pheasant M, Meyer L, Hsu F, Hinrichs AS, Harte RA, Giardine B, Fujita P, Diekhans M, Dreszer T, Clawson H, Kuhn RM, Barber GP, Haussler D, Kent WJ, Karolchik D, Zweig AS, Wang T, Smith KE, Rosenbloom KR, Rhead B, Raney BJ, Pohl A: The UCSC Genome Browser Database: update 2009. Nucleic Acids Res. 2009, 37: D755-761. 10.1093/nar/gkn875.
Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39: 1278-1284. 10.1038/ng2135.
Yu G, Xiao CL, Bo X, Lu CH, Qin Y, Zhan S, He QY: A new method for measuring functional similarity of microRNAs. JIOMICS. 2011, 1: 49-54.
Wang JZ, Du Z, Payattakool R, Yu PS, Chen CF: A new method to measure the semantic similarity of GO terms. Bioinformatics. 2007, 23: 1274-1281. 10.1093/bioinformatics/btm087.
Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S: GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics. 2010, 26: 976-978. 10.1093/bioinformatics/btq064.
Suzuki R, Shimodaira H: Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006, 22: 1540-1542. 10.1093/bioinformatics/btl117.
Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association. Bioinformatics. 2007, 23: 257-258. 10.1093/bioinformatics/btl567.
Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: D355-360. 10.1093/nar/gkp896.
Li C, Li X, Miao Y, Wang Q, Jiang W, Xu C, Li J, Han J, Zhang F, Gong B, Xu L: SubpathwayMiner: a software package for flexible identification of pathways. Nucleic Acids Research. 2009, 37: e131-10.1093/nar/gkp667.
Wickham H: ggplot2: Elegant Graphics for Data Analysis. 2009, Springer, 1
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Schliep KP: phangorn: phylogenetic analysis in R. Bioinformatics. 2011, 27: 592-593. 10.1093/bioinformatics/btq706.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.
Robinson DF, Foulds LR: Comparison of phylogenetic trees. Mathematical Biosciences. 1981, 53: 131-147. 10.1016/0025-5564(81)90043-2.
Park HJ, Williams TL: A Fitness Distance Correlation Measure for Evolutionary Trees. Proceedings of the 1st International Conference on Bioinformatics and Computational Biology. 2009, 331-342.
Tian E, Hagen KGT: O-linked glycan expression during Drosophila development. Glycobiology. 2007, 17: 820-827. 10.1093/glycob/cwm056.
Ludford-Menting MJ, Crimeen-Irwin B, Oliaro J, Pasam A, Williamson D, Pedersen N, Guillaumot P, Christansen D, Manie S, Gaus K, Russell SM: The Reorientation of T-Cell Polarity and Inhibition of Immunological Synapse Formation by CD46 Involves Its Recruitment to Lipid Rafts. Journal of Lipids. 2011, 2011: 1-10.
Amoroso R, Fitzsimmons L, Thomas WA, Kelly GL, Rowe M, Bell AI: Quantitative Studies of Epstein-Barr Virus-Encoded MicroRNAs Provide Novel Insights into Their Regulation. J Virol. 2011, 85: 996-1010. 10.1128/JVI.01528-10.
Umbach JL, Cullen BR: The role of RNAi and microRNAs in animal virus replication and antiviral immunity. Genes & Development. 2009, 23: 1151-1164. 10.1101/gad.1793309.
Umbach JL, Kramer MF, Jurak I, Karnowski HW, Coen DM, Cullen BR: MicroRNAs expressed by herpes simplex virus 1 during latent infection regulate viral mRNAs. Nature. 2008, 454: 780-783.
Portis T, Longnecker R: Epstein-Barr virus LMP2A interferes with global transcription factor regulation when expressed during B-lymphocyte development. J Virol. 2003, 77: 105-114. 10.1128/JVI.77.1.105-114.2003.
Miller DM, Cebulla CM, Rahill BM, Sedmak DD: Cytomegalovirus and transcriptional down-regulation of major histocompatibility complex class II expression. Semin Immunol. 2001, 13: 11-18. 10.1006/smim.2001.0291.
Samols MA, Skalsky RL, Maldonado AM, Riva A, Lopez MC, Baker HV, Renne R: Identification of Cellular Genes Targeted by KSHV-Encoded MicroRNAs. PLoS Pathog. 2007, 3: e65-10.1371/journal.ppat.0030065.
Narizhneva NV, Razorenova OV, Podrez EA, Chen J, Chandrasekharan UM, DiCorleto PE, Plow EF, Topol EJ, Byzova TV: Thrombospondin-1 up-regulates expression of cell adhesion molecules and promotes monocyte binding to endothelium. FASEB J. 2005, 19: 1158-1160.
Scaria V, Jadhav V: microRNAs in viral oncogenesis. Retrovirology. 2007, 4: 82-10.1186/1742-4690-4-82.
Feederle R, Linnstaedt SD, Bannert H, Lips H, Bencun M, Cullen BR, Delecluse HJ: A Viral microRNA Cluster Strongly Potentiates the Transforming Properties of a Human Herpesvirus. PLoS Pathog. 2011, 7: e1001294-10.1371/journal.ppat.1001294.
Philippe H, Douady CJ: Horizontal gene transfer and phylogenetics. Current Opinion in Microbiology. 2003, 6: 498-505. 10.1016/j.mib.2003.09.008.
Scaria V, Hariharan M, Maiti S, Pillai B, Brahmachari S: Host-virus interaction: a new role for microRNAs. Retrovirology. 2006, 3: 68-10.1186/1742-4690-3-68.
Pan X, Zhang B, San Francisco M, Cobb GP: Characterizing viral microRNAs and its application on identifying new microRNAs in viruses. J Cell Physiol. 2007, 211: 10-18. 10.1002/jcp.20920.
Weiler J, Hunziker J, Hall J: Anti-miRNA oligonucleotides (AMOs): ammunition to target miRNAs implicated in human disease?. Gene Ther. 2005, 13: 496-502.
Acknowledgements and Funding
This work was partially supported by the 2007 Chang-Jiang Scholars Program, "211" Projects, National "973" Projects of China (2011CB910700), National Natural Science Foundation of China (20871057), Guangdong Natural Science Research Grant (32209003), and the Fundamental Research Funds for the Central Universities (21611303 to G Yu and 11610101 to QY He).
The authors declare that they have no competing interests.
G Yu conceived and designed the prototype of the study, conducted the data analyses and drafted the manuscript. QY He supervised the study. All authors approved the final manuscript.
About this article
Cite this article
Yu, G., He, QY. Functional similarity analysis of human virus-encoded miRNAs. J Clin Bioinform 1, 15 (2011). https://doi.org/10.1186/2043-9113-1-15