Genomic variation in macrophage-cultured European porcine reproductive and respiratory syndrome virus Olot/91 revealed using ultra-deep next generation sequencing
© Lu et al.; licensee BioMed Central Ltd. 2014
Received: 3 December 2013
Accepted: 24 February 2014
Published: 4 March 2014
Porcine Reproductive and Respiratory Syndrome (PRRS) is a disease of major economic impact worldwide. The etiologic agent of this disease is the PRRS virus (PRRSV). Increasing evidence suggest that microevolution within a coexisting quasispecies population can give rise to high sequence heterogeneity in PRRSV.
We developed a pipeline based on the ultra-deep next generation sequencing approach to first construct the complete genome of a European PRRSV, strain Olot/9, cultured on macrophages and then capture the rare variants representative of the mixed quasispecies population. Olot/91 differs from the reference Lelystad strain by about 5% and a total of 88 variants, with frequencies as low as 1%, were detected in the mixed population. These variants included 16 non-synonymous variants concentrated in the genes encoding structural and nonstructural proteins; including Glycoprotein 2a and 5.
Using an ultra-deep sequencing methodology, the complete genome of Olot/91 was constructed without any prior knowledge of the sequence. Rare variants that constitute minor fractions of the heterogeneous PRRSV population could successfully be detected to allow further exploration of microevolutionary events.
KeywordsPRRSV Microevolution Variant spectra Ultra-deep next generation sequencing
Porcine Reproductive and Respiratory Syndrome virus (PRRSV) is the causative agent of a significant disease of the domestic pig (Sus scrofa) with global consequences. The severity of PRRSV infection ranges from subclinical to lethal and it affects pigs in both growing and reproductive stages. The virus has a positive-sense 15 kb RNA genome and its genetic diversity has been well characterised within and between European and North American strains . Extensive viral genetic heterogeneity may have contributed towards the observed variations between PRRSV isolates and clones in term of virulence, interactions with the immune system, and antigenic properties of viral proteins. Such a broad diversity indeed poses serious challenges to diagnostics and control measures.
Most previous studies of PRRSV genetic diversity have been restricted to the ORF5 and ORF7 sequences of type 2 “North American-like” viruses that also include the Asian variants. Only 14 of the 303 completed PRRSV genomes in Genbank belong to genotype 1. Furthermore, studies have shown that PRRSV mutates rapidly and multiple intra-strain variants can coexist in individually infected pigs . The extensive genetic diversity displayed by PRRSV and other RNA viruses such as HIV and influenza reflects the error prone nature of RNA polymerases, which lack a proofreading function [3, 4].
To identify PRRSV quasispecies, previous studies have employed conventional methodologies including reverse-transcription, PCR, cloning and Sanger sequencing of a subset of PRRSV structural and non-structural proteins [2, 5]. More recently next-generation sequencing (NGS) of fragments generated by long range RT-PCR has been used to characterise multiple PRRSV genomes . However, this approach relies upon prior knowledge of the target sequence and the assumption that the PCR primer binding sites are non-variable. Here we describe an approach which requires no prior knowledge of the target sequences and which should enable the detection of low-frequency nucleotide variants and hence provides a snapshot of the microevolution in the entire viral population.
We analysed the intra-strain sequence diversity of low passage PRRSV Olot/91 strain, passaged exclusively on primary porcine alveolar macrophages (PAM). First reported in Spain in 1991, Olot/91 is the parent strain of the commercial Suvaxyn PRRSV inactivated vaccine which is used in Spain and Portugal. Only a partial sequence of 3,383 nt [GenBank:X92942], that covers ORFs 2-7 and the 3'-UTR of this strain has previously been published .
While variants in the major Olot/91 strain do congregate on known highly hypervariable B- and T-cell epitopes of type I PRRSV (Figure 2 and Additional file 3: Figure S2) [19–23], most of the variants representative of the probable quasispecies population map to the nsp12 and gp2a genes of ORF1b and ORF2 respectively (Figure 2); suggesting a higher microevolutionary rate at this region under the investigated culture conditions. In addition, the predicted short signal peptide of GP2a was found to harbour a high concentration of non-synonymous variants and a potential N-glycosylation site was created from a residue change (S37N) on GP5. Further analysis, like molecular modelling, may be necessary to decipher the potential impacts of these evolving variants have on such functions as interactions between these proteins and other viral or host cellular proteins.
Using ultra-deep NGS we have constructed the complete sequence of the PRRSV Olot/91 genome cultured in their natural host - porcine alveolar macrophages - and identified single-nucleotide variants present in the associated viral population. Further investigation using this methodology will help to establish if a link exists between the microevolutionary dynamics and pathogenesis of PRRSV viral strains. However, it is also important to note that the exact nature of the pathogenesis cannot be truly identified until the viral haplotypes within the quasispecies population can be reconstructed with confidence .
We are indebted to Lise K. Kvisgaarda and Lars E. Larsen for kindly providing the sequence of the Marc 145-adapted Olot/91 strain. TAA, ALA and ZHL were supported by UK Biotechnology and Biological Sciences Research Council (BBRSC) Institute Strategic Programme Grants (BB/J004235/1). ADW was supported by EC FP7 PoRRScon Grant (Agreement number 245141). Edinburgh Genomics is funded by BBSRC National Capacity Grants (BB/J004243/1). We are also grateful for support from the COST Action FA0902.
- Murtaugh MP, Stadejek T, Abrahante JE, Lam TT, Leung FC: The ever-expanding diversity of porcine reproductive and respiratory syndrome virus. Virus Res 2010, 154: 18-30. 10.1016/j.virusres.2010.08.015PubMedView Article
- Goldberg TL, Lowe JF, Milburn SM, Firkins LD: Quasispecies variation of porcine reproductive and respiratory syndrome virus during natural infection. Virology 2003, 317: 197-207. 10.1016/j.virol.2003.07.009PubMedView Article
- Lauring AS, Andino R: Quasispecies theory and the behavior of RNA viruses. PLoS Pathogens 2010, 6: e1001005. 10.1371/journal.ppat.1001005PubMedPubMed CentralView Article
- Domingo E, Sheldon J, Perales C: Viral quasispecies evolution. Microbiol Mol Biol Rev 2012, 76: 159-216. 10.1128/MMBR.05023-11PubMedPubMed CentralView Article
- Schommer SK, Kleiboeker SB: Use of a PRRSV infectious clone to evaluate in vitro quasispecies evolution. Adv in Exp Med and Biol 2006, 581: 435-438. 10.1007/978-0-387-33012-9_78View Article
- Kvisgaard LK, Hjulsager CK, Fahnøe U, Breum SØ, Ait-Ali T, Larsen LE: A fast and robust method for full genome sequencing of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Type 1 and Type 2. J Virol Methods 2013, 193: 697-705. 10.1016/j.jviromet.2013.07.019PubMedView Article
- Urniza A, Climent I, Duran J, Corte´s E, Sarraseca J, Casal JI, Vela C: Baculovirus expression of proteins of porcine reproductive and respiratory syndrome virus strain olot/91. Involvement of ORF3 and ORF5 proteins in protection. Virus Genes 1997, 14: 19-29. 10.1023/A:1007931322271PubMedView Article
- Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, Rogel-Gaillard C, Park C, Milan D, Megens HJ, Li S, Larkin DM, Kim H, Frantz LA, Caccamo M, Ahn H, Aken BL, Anselmo A, Anthon C, Auvil L, Badaoui B, Beattie CW, Bendixen C, Berman D, Blecha F, Blomberg J, Bolund L, Bosse M, Botti S, Bujie Z, et al.: Analyses of pig genomes provide insight into porcine demography and evolution. Nature 2012, 491: 393-398. 10.1038/nature11622PubMedPubMed CentralView Article
- BMTagger: Best match tagger for removing human reads from metagenomics datasets. ftp://ftp.ncbi.nlm.nih.gov/pub/agarwala/bmtagger/
- Sickle - a windowed adaptive trimming tool for FASTQ files using quality http://github.com/najoshi/sickle
- Scythe - a bayesian adapter trimmer http://github.com/vsbuffalo/scythe
- Meulenberg JJM, Hulst MM, de Meijer EJ, Moonen PLJM, den Besten A, de Kluyver EP, Wensvoort G, Moormann RJM: Lelystad Virus, the Causative Agent of Porcine Epidemic Abortion and Respiratory Syndrome (PEARS), Is Related to LDV and EAV. Virology 1993, 192: 62-72. 10.1006/viro.1993.1008PubMedView Article
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324PubMedPubMed CentralView Article
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing Subgroup: The sequence alignment/Map format and SAMtools. Bioinformatics 2078–2079, 2009: 25.
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 2008, 18: 821-829. 10.1101/gr.074492.107PubMedPubMed CentralView Article
- Wilm A, Aw PP, Bertrand D, Yeo GH, Ong SH, Wong CH, Khor CC, Petric R, Hibberd ML, Nagarajan N: LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res 2012, 40: 11189-11201. 10.1093/nar/gks918PubMedPubMed CentralView Article
- Delrue I, Van Gorp H, Van Doorsselaere J, Delputte PL, Nauwynck HJ: Susceptible cell lines for the production of porcine reproductive and respiratory syndrome virus by stable transfection of sialoadhesin and CD163. BMC Biotechnol 2010, 10: 48. 10.1186/1472-6750-10-48PubMedPubMed CentralView Article
- Ait-Ali T, Wilson AD, Carre W, Westcott DG, Frossard JP, Mellencamp MA, Mouzaki D, Matika O, Waddington D, Drew TW, Bishop SC, Archibald AL: Host inhibits replication of European porcine reproductive and respiratory syndrome virus in macrophages by altering differential regulation of type-I interferon transcriptional response. Immunogenetics 2011, 63: 437-448. 10.1007/s00251-011-0518-8PubMedView Article
- Díaz I, Pujols J, Ganges L, Gimeno M, Darwich L, Domingo M, Mateu E: In silico prediction and ex vivo evaluation of potential T-cell epitopes in glycoproteins 4 and 5 and nucleocapsid protein of genotype-I (European) of porcine reproductive and respiratory syndrome virus. Vaccine 2009, 27: 5603-5611. 10.1016/j.vaccine.2009.07.029PubMedView Article
- Oleksiewicz MB, Bøtner A, Normann P: Porcine B-cells recognize epitopes that are conserved between the structural proteins of American- and European-type porcine reproductive and respiratory syndrome virus. J Gen Virol 2002, 83: 1407-1418.PubMedView Article
- Costers S, Vanhee M, Van Breedam W, Van Doorsselaere J, Geldhof M, Nauwynck HJ: GP4-specific neutralizing antibodies might be a driving force in PRRSV evolution. Virus Res 2010, 154: 104-113. 10.1016/j.virusres.2010.08.026PubMedView Article
- Parida R, Choi I-S, Peterson DA, Pattnaik AK, Laegreid W, Zuckermann FA, Osorio FA: Location of T-cell epitopes in nonstructural proteins 9 and 10 of type-II porcine reproductive and respiratory syndrome virus. Virus Res 2012, 169: 13-21. 10.1016/j.virusres.2012.06.024PubMedView Article
- Liao YC, Lin HH, Lin CH, Chung WB: Identification of cytotoxic T lymphocyte epitopes on swine viruses: multi-epitope design for universal T cell vaccine. PLoS One 2013, 8: e84443. 10.1371/journal.pone.0084443PubMedPubMed CentralView Article
- Schirmer M, Sloan WT, Quince C: Benchmarking of viral haplotype reconstruction programmes: an overview of the capacities and limitations of currently available programmes. Brief Bioinform 2012. doi:10.1093/bib/bbs081
- Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies. Bioinformatics 2005, 21: 676-679. 10.1093/bioinformatics/bti079PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.