Evaluation of the capability of the PCV2 genome to encode miRNAs: lack of viral miRNA expression in an experimental infection
Veterinary Research volume 46, Article number: 48 (2015)
Porcine circovirus type 2 (PCV2) is a ssDNA virus causing PCV2-systemic disease (PCV2-SD), one of the most important diseases in swine. MicroRNAs (miRNAs) are a new class of small non-coding RNAs that regulate gene expression post-transcriptionally. Viral miRNAs have recently been described and the number of viral miRNAs has been increasing in the past few years. In this study, small RNA libraries were constructed from two tissues of subclinically PCV2 infected pigs to explore if PCV2 can encode viral miRNAs. The deep sequencing data revealed that PCV2 does not express miRNAs in an in vivo subclinical infection.
Introduction, methods, and results
Porcine circovirus type 2-systemic disease (PCV2-SD) is a devastating disease that causes important economic losses . The disease is essentially caused by PCV2, a single stranded DNA, non enveloped virus belonging to the Circoviridae family . The PCV2 genome encodes four ORFs. ORF1 encodes for two proteins (Rep and Rep’) which are involved in replication. ORF2 encodes for the Cap protein, which constitutes the unique structural protein. The two proteins coded by ORF3 and ORF4 have been related with cellular apoptosis while the function of both still needs further research .
MicroRNAs (miRNAs) comprise a class of small noncoding RNAs that post-transcriptionally regulate the expression of many genes by mRNA degradation or translation inhibition . miRNAs are involved in the modulation of gene expression and replication of many viruses and play a pivotal role in host-virus interactions. In addition, many viruses encode miRNAs that can play a role in the infection process. To date, 308 hairpin precursors and 502 mature miRNAs have been discovered in several viruses as shown in the miRBase , most of them encoded by herpesviruses, which is the virus family with the highest miRNAs encoding capacity [6,7]. Other viruses belonging to the families Polyomaviridae, Adenoviridae, Papillomaviridae, Baculoviridae and Ascoviridae encode miRNAs in low numbers [8-12]. Recently, a Human Torque Teno virus, a small, ssDNA virus from the Anelloviridae family, that encodes a miRNA involved in interferon modulation has been described . All of them share in common to be DNA viruses, with an essential nuclear phase in their replication cycle, necessary for initial miRNA biogenesis. miRNAs encoded by RNA viruses is a matter of controversy [14-16].
In the present study, the expression of miRNAs in subclinically PCV2 infected pigs was analysed using high throughput sequencing. Firstly, in silico prediction was carried out in order to check if the PCV2 genome encodes possible miRNA precursors. The Vmir prediction algorithm  was used to predict the possible presence of hairpin structures in the PCV2 genome compatible with the existence of miRNAs. Computational prediction of viral miRNAs indicated that 41 miRNA candidates could be identified, 16 of them with a score among 100-150 and two with a score >150 (Figure 1). In order to explore whether these candidates were present in the virus, next-generation sequencing (NGS) of small RNAs was carried out from tonsil and mediastinal lymph node of animals subclinically infected with PCV2.
We have reported previously modification of the expression pattern of host miRNAs due to the PCV2 infection . For both purposes, four animals were inoculated with 7 × 104.8TCID50 of PCV2 isolate Sp-10-7-54-13  and two animals were inoculated with PBS. Animal experiments were carried out at the CReSA facilities, all procedures were performed under the supervision of the Ethical and Animal Welfare Committee of the UAB in compliance with national guidelines and EU regulations. At 21 days post-inoculation samples were taken and total RNA extractions were carried out in order to construct small RNA libraries as described in  with some modifications. A total of 12 small RNA libraries were constructed in a two-step ligation procedure with the 3′ and 5′ adaptors from IDT technologies. Amplification by RT-PCR was carried out with fusion primers containing sequences complementary to the 3′ and 5′ adaptors used for miRNA library construction and sequences complementary to the A and B adaptors used for high-throughput (HT) sequencing with the GS FLX 454 device (Roche) at CRAG (Centre de Recerca Agrogenòmica, Universitat Autònoma de Barcelona, Spain). From the total reads obtained (1 106 437), primer sequences were trimmed and only those insert sequences between 15 and 29 nucleotides and with total number of sequences ≥3 were kept for further analysis. This procedure resulted in a total of 796 710 reads.
The analysis of the reads that aligned to the Sus scrofa genome are reported elsewhere because these sequences constitute porcine miRNAs . For viral miRNA discovery, sequences were blasted to the PCV2 isolate Sp-10-7-54-13 genome (NCBI Reference Sequence: GU049342) considering only sequences with 100% of alignment and identity (perfect match). In order to search potential viral miRNAs, a blast against the viral genome was done with an increased number of mismatches in the extremes due to miRNA variability (isomiRs) [21,22]. The presence of isomiRs entails differences in length and point mutations in both extremes, with the 3′-terminus having a higher proportion of such mutations. In some cases, a blast was done allowing internal variations with a <100% of alignment in order to consider the variability of the viral genome . Also, sequences were blasted to the output of Vmir hairpin PCV2 structures.
No hits were found in the viral PCV2 genome sequence (GU049342) with a 100% homology. One candidate of 18 nt with 58 copies was identified when allowing <100% of alignment and identity in the ORF2, positions 1189 to 1206. In accordance with the Vmir prediction, these positions were included in the pre-miRNA candidate MD18, (the hairpin with the highest score), presenting one internal mismatch, and 94.4% sequence identity. The secondary structure and the minimum free energy of the pre-miRNA candidate was analysed with MFold software  (Figure 2). As this candidate presented a point mutation, and due to the described variability of the virus, the corresponding viral DNA fragment from the mediastinal lymph node of an infected animal was amplified and Sanger sequenced. Because the exact candidate sequence was not found in the viral sequence (94.4% sequence identity), it was not considered a miRNA encoded by the virus. The miRNA candidate sequence was compared to the porcine genome Sscrofa10.2 (GCA_000003025.4) using blast and showed 100% sequence identity with the ssc-miR-29a hairpin precursor located in chromosome 18. Also, the 18 nt miRNA candidate sequence was compared to the miRBase (v.21) showing 100% sequence identity with miR-29a-5p. The alignment of this region was carried out including all PCV2 sequences available in the databases (Figure 3). miR-29a-5p has been described in many species like human, bovine, mouse, but it has not been described in pigs, where the precursor miR-29 has been included in miRBase along with the mature ssc-miR-29a-3p. From the total reads obtained in this study, 1276 sequences of 22 nt, comprising the 18 nt of the candidate, blasted to the miR-29a-5p. From all of the above data, miR-29a-5p can be considered as a miRNA encoded by the porcine genome and the 18 nt candidate as one of the miR-29a-5p isomiRs.
This is the first study that tries to identify if PCV2 can encode miRNAs. PCV2 is a ssDNA virus with a nuclear phase in its replicating cycle . This step is considered essential for the production of viral miRNAs, considering that the first steps in miRNA maturation take place in the cellular nucleus . Nevertheless, high throughput sequencing has failed to identify any miRNA encoded by the viral genome in the natural host after an experimental infection. Thus, PCV2 may increase the list of DNA viruses not encoding miRNAs, as Cowpox virus . Notwithstanding, the capacity of PCV2 to encode miRNAs has to be evaluated in cell culture, in a different clinical form or at a different time points.
Only one miRNA candidate was initially identified, but the posterior analysis indicated that it was a host miRNA, miR-29a-5p. The homology of the viral sequence with miR-29a-5p, with only one internal mismatch, led us to analyze with more detail if miR-29a-5p can regulate the expression of the Cap protein. A miR-29a-5p target prediction was evaluated using Miranda software  in order to explore if the Cap gen constitutes a target. As expected, a predicted target was identified in ORF2, with a low free energy of -21.65 Kcal/mol and a high complementarity score of 187. The alignment of the region with all PCV2 isolates available in the database indicated that A, C and T, but not G, has been found at position 1202 respecting Sp-10-7-54-13. All substitutions lead to a non synonymous change. The lack of the presence of G at position 1202 could be due to restriction in the protein conformation or due to the pressure exerted by miR-29a-5p in order to avoid regulation by this miRNA. If this is a consequence of viral evolution as  proposed, it needs further investigation. How viruses can evolve to avoid the inhibition by host miRNAs is a critical question. Some authors indicate that this evolution allows viruses to replicate without being targeted by host miRNAs by encoded viral miRNAs or modifying the expression of host miRNAs . On the other hand, the secondary structure of MD18 could avoid the regulation by miR-29a-5p due to the inaccessibility of the target sequence as has been proposed for HIV .
In a previous study , in subclinically infected pigs, we demostrated that PCV2 can alter the miRNA expression pattern of the host. If PCV2 can express miRNAs in vitro has to be determined, but in its natural host, in a subclinical infection, NGS failed to indentify viral miRNAs. The exploration of the possible capacity of PCV2 to encode miRNAs could contribute to the understanding of the pathogenesis of PCV2, especially for the candidate with highest score identified by the in silico prediction. In addition, further studies on the similarity of this candidate with miR-29a-5p, and its significance, could shed light on how miRNAs affect viral evolution .
Segalés J (2012) Porcine circovirus type 2 (PCV2) infections: clinical signs, pathology and laboratory diagnosis. Virus Res 164:10–19
Segalés J, Allan GM, Domingo M (2005) Porcine circovirus diseases. Anim Health Res Rev 6:119–142
Gao Z, Dong Q, Jiang Y, Opriessnig T, Wang J, Quan Y, Yang Z (2014) ORF4-protein deficient PCV2 mutants enhance virus-induced apoptosis and show differential expression of mRNAs in vitro. Virus Res 183:56–62
Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116:281–297
miRBase: the microRNA database [http://www.mirbase.org]
Grundhoff A, Sullivan CS (2011) Virus-encoded microRNAs. Virology 411:325–343
Skalsky RL, Cullen BR (2010) Viruses, microRNAs, and host interactions. Annu Rev Microbiol 64:123–141
Sullivan CS, Grundhoff AT, Tevethia S, Pipas JM, Ganem D (2005) SV40-encoded microRNAs regulate viral gene expression and reduce susceptibility to cytotoxic T cells. Nature 435:682–686
Sano M, Kato Y, Taira K (2006) Sequence-specific interference by small RNAs derived from adenovirus VAI RNA. FEBS Lett 580:1553–1564
Qian K, Pietila T, Ronty M, Michon F, Frilander MJ, Ritari J, Tarkkanen J, Paulin L, Auvinen P, Auvinen E (2013) Identification and validation of human papillomavirus encoded microRNAs. PLoS One 8:e70202
Hussain M, Taft RJ, Asgari S (2008) An insect virus-encoded microRNA regulates viral replication. J Virol 82:9164–9170
Singh J, Singh CP, Bhavani A, Nagaraju J (2010) Discovering microRNAs from Bombyx mori nucleopolyhedrosis virus. Virology 407:120–128
Kincaid RP, Burke JM, Cox JC, de Villiers EM, Sullivan CS (2013) A human torque teno virus encodes a microRNA that inhibits interferon signaling. PLoS Pathog 9:e1003818
Ouellet DL, Plante I, Landry P, Barat C, Janelle ME, Flamand L, Tremblay MJ, Provost P (2008) Identification of functional microRNAs released through asymmetrical processing of HIV-1 TAR element. Nucleic Acids Res 36:2353–2365
Schopman NC, Willemsen M, Liu YP, Bradley T, van Kampen A, Baas F, Berkhout B, Haasnoot J (2012) Deep sequencing of virus-infected cells reveals HIV-encoded small RNAs. Nucleic Acids Res 40:414–427
Whisnant AW, Bogerd HP, Flores O, Ho P, Powers JG, Sharova N, Stevenson M, Chen CH, Cullen BR (2013) In-depth analysis of the interaction of HIV-1 with cellular microRNA biogenesis and effector mechanisms. MBio 4:e000193
Grundhoff A (2011) Computational prediction of viral miRNAs. Methods Mol Biol 721:143–152
Nunez-Hernandez F, Perez LJ, Munoz M, Vera G, Tomas A, Egea R, Cordoba S, Segales J, Sanchez A, Nuñez JI (2015) Identification of microRNAs in PCV2 subclinically infected pigs by high throughput sequencing. Vet Res 46:18
Fort M, Sibila M, Nofrarias M, Perez-Martin E, Olvera A, Mateu E, Segalés J (2010) Porcine circovirus type 2 (PCV2) Cap and Rep proteins are involved in the development of cell-mediated immunity upon PCV2 infection. Vet Immunol Immunopathol 137:226–234
Timoneda O, Nunez-Hernandez F, Balcells I, Munoz M, Castello A, Vera G, Perez LJ, Egea R, Mir G, Cordoba S, Rosell R, Segalés J, Tomas A, Sanchez A, Nunez JI (2014) The role of viral and host microRNAs in the Aujeszky’s disease virus during the infection process. PLoS One 9, e86965
Ebhardt HA, Tsang HH, Dai DC, Liu Y, Bostan B, Fahlman RP (2009) Meta-analysis of small RNA-sequencing errors reveals ubiquitous post-transcriptional RNA modifications. Nucleic Acids Res 37:2461–2470
Morin RD, O’Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, Eaves CJ, Marra MA (2008) Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res 18:610–621
Firth C, Charleston MA, Duffy S, Shapiro B, Holmes EC (2009) Insights into the evolutionary history of an emerging livestock pathogen: porcine circovirus 2. J Virol 83:12813–12821
Zuker M (2003) Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 31:3406–3415
Tischer I, Peters D, Pociuli S (1995) Occurrence and role of an early antigen and evidence for transforming ability of porcine circovirus. Arch Virol 140:1799–1816
Cullen BR (2010) Five questions about viruses and microRNAs. PLoS Pathog 6:e1000787
Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS (2003) MicroRNA targets in Drosophila. Genome Biol 5:R1
Bogerd HP, Skalsky RL, Kennedy EM, Furuse Y, Whisnant AW, Flores O, Schultz KL, Putnam N, Barrows NJ, Sherry B, Scholle F, Garcia-Blanco MA, Griffin DE, Cullen BR (2014) Replication of many human viruses is refractory to inhibition by endogenous cellular microRNAs. J Virol 88:8065–8076
Cullen BR (2013) How do viruses avoid inhibition by endogenous cellular microRNAs? PLoS Pathog 9:e1003694
Guo YE, Steitz JA (2014) Virus meets host microRNA: the destroyer, the booster, the hijacker. Mol Cell Biol 34:3780–3787
The authors thank all the personnel of the animal facilities at CReSA. This work was supported by the projects AGL2007-66371-C02 and AGL2010-22358-C02 from MINECO (Spanish government). FN-H is the recipient of a PhD fellowship FI-DGR (AGAUR) from Generalitat de Catalunya. The authors thank Dr. Kevin Dalton for editing the manuscript.
The authors declare that they have no competing interests.
JIN and AS conceived, designed and coordinated this study. FNH, and LJP performed the experiments. JIN and FNH analyzed data. JS, and SC, contributed reagents/material/analysis tools. JIN and FNH wrote the paper. GV performed the bioinformatics analysis. JS critically read the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Núñez-Hernández, F., Pérez, L.J., Vera, G. et al. Evaluation of the capability of the PCV2 genome to encode miRNAs: lack of viral miRNA expression in an experimental infection. Vet Res 46, 48 (2015). https://doi.org/10.1186/s13567-015-0181-4
- Torque Teno Virus
- Porcine Circovirus Type
- PCV2 Infection
- Viral miRNAs
- Porcine Genome