- Short report
- Open Access
Genetic diversity and multiple introductions of porcine reproductive and respiratory syndrome viruses in Thailand
Virology Journal volume 8, Article number: 164 (2011)
Porcine reproductive and respiratory syndrome virus (PRRSV) is prevalent in Thailand, causing a huge impact on the country's swine industry. Yet the diversity and origin of these Thai PRRSVs remained vague. In this context, we collected all the Thai PRRSV sequences described earlier and incorporated them into the global diversity. The results indicated that PRRSVs in Thailand were originated from multiple introductions involving both Type 1 and Type 2 PRRSVs. Many of the introductions were followed by extensive geographic expansion, causing regional co-circulation of diverse PRRSV variants in three major pig-producing provinces. Based on these results, we suggest (1) to avoid blind vaccination and to apply vaccines tailor-made for target diversity, (2) to monitor pig importation and transportation, and (3) to implement a better biosecurity to reduce horizontal transmissions as three potentially effective strategies of controlling PRRS in Thailand.
PRRSV is a major swine disease virus causing economic losses in swine industry worldwide including Thailand. This disease was first reported during almost concurrent epidemics in the North American countries (late 1980s) [1, 2] and in the European countries (early 1990) . The causative agent, porcine reproductive and respiratory syndrome virus (PRRSV) belongs to the family Arteriviridae in the order Nidovirales . It is an enveloped virus containing a positive-sense RNA genome of approximately 15 kb in length, encoding at least nine open reading frames (ORFs) including ORF1a, 1b, 2a, 2b, and 3-7 . Among them, ORF5 encoding the major envelope protein is often used for phylogenetic analysis because of its high variability. PRRSV can be divided into two major genotypes: Type 1 (European strains) and Type 2 (North American strains). Both genotypes are found to be genetically and antigenically heterogeneous [6, 7]. According to the recent reports for PRRSV classification, Type 1 PRRSV is divided into 3 subtypes. Among them, the cosmopolitan Subtype I was further divided into 12 different clades [8, 9]. For Type 2 PRRSV, 9 well-defined lineages have been described [8, 10].
PRRSV has been circulating in Thailand for a long time. Seropositive animals could be traced back to as early as 1989, and seropositive rate increased annually from 1991 to 2002 . In 1996, the first Thai PRRSV was isolated and was identified as Type 2 PRRSV . A few years later, Type 1 PRRSV was also reported in Thailand, and some was found to be co-circulating with Type 2 PRRSV within the same herd . Despite of the effort on characterizing PRRSV diversity in Thailand, the origin and epidemiological history of the viruses remain unknown. This study re-analyzes previously characterized Thai sequences in the context of global PRRSV diversity, which helps shed lights on the potential origin and prevalence of different variants of PRRSVs in Thailand.
2. Materials and Methods
2.1 Viral isolation
All Thai PRRS viruses were isolated from nursery pigs (3-8 week of age) showing sign of PRDC in different regions of Thailand during 2000 to 2008. All pig farms are multi-site farms and practicing partial all in-all out and mostly continuous-flow system. The infected pigs or tissue samples or sera were submitted to the Chulalongkorn University-Veterinary Diagnosis Laboratory (CU-VDL). Virus isolation was performed by using MARC-145 cells and PRRSV positive result was confirmed by PCR and Immunoperoxidase monolayer assay (IPMA) using SDOW17 [11, 13].
2.2 Data sets
Four hundred and twenty two complete Type 1 and five hundred and eighty one complete Type 2 PRRSV ORF5 gene sequences were downloaded from Gene Bank and PRRSV database http://prrsvdb.org/ as reference backbone. The sequences are comprised of sequences representative of diversity of Type 1  and Type 2  and all available Thai PRRSV ORF5 sequences (Table 1). All Thai PRRSVs were sequenced at CU-VDL during 2000 to 2008 [11, 13].
2.3 Phylogenetic reconstruction
The ORF5 nucleotide sequences were aligned in MUSCLE v3.6 , using default settings with minor manual adjustments. The alignments were then screened for recombination using Recombination Detection Program v.2 . No recombination was found in the alignment. Then, a phylogeny was constructed using a Bayesian Markov Chain Monte Carlo (BMCMC) method implemented in MrBayes v.3.2 . A general time-reversible nucleotide substitution model with 4 categories of gamma-distributed rate heterogeneity and a proportion of invariant sites (GTR+ Γ4+I) were used. The posterior distribution of trees and model parameters were summarized from Markov Chain Monte Carlo sampling over 10 million generations, during which trees were sampled every 1000 generations. We adopted the classification system from the previous report [8, 10] for lineage diversity of Thai sequences. Average pairwise genetic distances between two lineages were calculated to evaluate the interlineage genetic diversity.
3.1 Type 1 PRRSV in Thailand
Phylogenetic analysis indicated that there were at least four independent introductions of Type 1 PRRSV into Thailand during 2000 to 2008. Among these introductions, three of them were from Clade A (Lelystad-like) while the other from Clade H. The definite sources of these introductions are difficult to identify due to complex geographic constituents for both clades (Figure 1A). Nevertheless, one Thai sequence (02CB12, 2002) was found to be closely related to the prototype (LV) with 98.5% of similarity as well as Porcilis® PRRS with 99.1% of similarity (Intervet, The Netherlands) in Clade A (Figure 1A) indicating virus sample 02CB12 is probably a vaccine re-isolate or a vaccine descendant. Despite multiple introductions, the prevalence and diversity of Type 1 PRRSV in Thailand was quite limited except for those within Clade H (n = 6), which had an average genetic distance of 7% and a maximum pair of genetic distance of 10%.
3.2 Type 2 PRRSV in Thailand
Type 2 PRRSV in Thailand was originated from at least three independent introductions, which is evidenced by three monophyletic clusters: one (n = 9) within Lineage 1 while another one (n = 1) within Lineage 5 and the other (n = 3) as an orphan cluster that is beyond the description of current classification system by Shi et al (2010b). The phylogenetic backbone suggested this orphan cluster had closer relation to Lineage 6-9 forming a monophyl supported by high posterior probability (0.93). Nevertheless, further classification for this cluster remains unattainable. Since this orphan cluster had larger than 10% average genetic distances with those established lineages, we suggest this cluster as an independent lineage. Moreover, Thai prototype of Type 2 PRRSV (01NP1.2) has 100% nucleotide identity with MLV vaccine in Lineage 5 indicating an apparent vaccine re-isolate or vaccine descendant.
The major proportion of Thai Type 2 PRRSVs formed a single cluster within Lineage 1 (Figure 1B). This cluster had the largest diversity comparing to other Thai clusters, with an average genetic distance of 9% and a maximum pair of genetic distance of 14%. In addition, our estimation of the tMRCA of this cluster suggested an introduction time of around early 1990s if not earlier (Tun, personal communication).
3.3 Geographic distribution of PRRSV in Thailand
In order to examine the distribution spread of PRRSVs after their initial introduction into Thailand, geographical locations with different shapes were marked as previously described (in section 3.1 and 3.2) as Type 1 and 2 Thai clusters. Each represented an independent introduction event (Figure 2). For some provinces (Chonburi, Ratchaburi and Nakhon Pathom) in central Thailand, co-circulation of multiple clusters from the two types of PRRSV was observed whereas others had only one cluster of PRRSV. In addition, it is worth mentioning that the Thai cluster within Lineage 1 had the widest distribution, spreading to at least 8 provinces in Thailand (Figure 2).
The appearance of the vaccine-related strains in Thailand is of interest since the Porcilis® PRRS is not commercially available in Thailand. One possibility is the smuggling or unauthorized use of this modified live vaccine. Alternatively, it might originate from previous introduction of vaccinated breeder animals. In Thailand, most of the swine breeders and frozen semen have been continuously imported from both European and North American countries. Evidences of seropositivity have been reported in the imported pigs and semen from previous retrospective study . Despite of the presence of vaccine-related PRRSVs in Thailand, their impact to the swine industry was most likely insignificant. It is because for each vaccine-related variant, only one sequence was identified in the phylogeny. And no transmission chain has been identified to be seeded by these vaccine strains..
Interestingly, three Type 2 PRRSV isolates collected from 2007 and 2008 belong to an orphan lineage which has limited diversity but is well isolated from any of the established lineages. This is also supported by NSP2 phylogeny, where the Thai orphan lineage (08NP147) is distantly related to other Type 2 PRRSV . The missing diversity between the orphan lineage and other lineages are most likely due to incomplete sampling in Thailand. Indeed, the sampling in this study is limited by not only quantity, but also time span, compared to that of North America. Therefore, we suggest more intensive sampling and characterization of PRRSV in Thailand to fill in the diversity gaps.
Our recent biogeographic analyses on Type 2 PRRSV suggest that Lineage 6-9 were mostly endemic in the United States whereas Lineage 1 and 2 were originally circulating in Canada and later emerged in the United States due to intensive pig flow from Canada to the United States (Shi, unpublished data). Since the Thai orphan lineage is clustered within Lineage 6-9, it was most likely introduced from the United States. On the other hand, the Thai cluster within Lineage 1 was potentially introduced from Canada or United States. The Canadian origin of this cluster is more favored, because its common ancestor was estimated as early 1990s, which is earlier than the emergences of Lineage 1 and 2 PRRSV in the United States. Yet this perspective remains to be examined by incorporating earlier Canadian samples into the phylogenetic analyses.
Chonburi, Ratchaburi and Nakhon Pathom are three major Thai swine producing provinces located in the central region. Evidence based on our phylogenetic analyses show co-circulation of multiple clusters of both genotypes of PRRSV in these three provinces. This co-circulation of multiple PRRS viruses renders huge burden on the swine industry in these areas, because it might cause succession of PRRSV infection (due to lack of cross-protection for heterologous infection) and recombination. Since our analyses are based on ORF5 sequence alignment, no recombination events were detected but potential recombination breakpoint may exist in other part of the genome. In addition, we are able to identify a widely distributed Thai cluster (Lineage 1 of Type 2) originated from an introduction of approximately 20 years ago. Additionally, it is likely the movement of the PRRSV infected pigs between the herds that maintain the transmission chain of this particular cluster.
Our study demonstrates multiple introductions of PRRSV strains into Thailand, which result in co-circulation of diverse PRRSV strains in some areas. This posts difficulty on the vaccine selection since current monovalent vaccine has been shown to provide no or only partial protection against heterologous PRRSV strains [6, 18, 19]. Therefore, correct choice of vaccine is essential for prevention and control of PRRSV in Thailand. In addition to vaccination, the virus control measures should also aim for pig importation and transportation. First, an upgraded biosecurity system is required to prevent not only the introduction of new types of PRRSV into Thailand, but also the subsequent dissemination of the viruses among major pig-producing provinces. Second, it is necessary to establish a surveillance system which monitors the pig flows. This helps provide timely remedy in cases of outbreak.
porcine reproductive and respiratory syndrome virus
open reading frame
porcine respiratory disease complex
Chulalongkorn University-Veterinary Diagnosis Laboratory
Immunoperoxides monolayer assay
: Bayesian Markov Chain Monte Carlo
Hill H: Overview and history of mystery swine disease (swine infertility respiratory syndrome). Proceedings of the Mystery Swine Disease Committee Meeting; October 6; Denver Colorado 1990, 29-30.
Keffaber K: Reproductive failure of unknown etiology. Am Assoc Swine Pract Newsletter 1989, 1: 1-9.
Wensvoort G, Terpstra C, Pol JM, ter Laak EA, Bloemraad M, de Kluyver EP, Kragten C, van Buiten L, den Besten A, Wagenaar F, et al.: Mystery swine disease in The Netherlands: the isolation of Lelystad virus. Vet Q 1991, 13: 121-130.
Cavanagh D: Nidovirales: a new order comprising Coronaviridae and Arteriviridae. Arch Virol 1997, 142: 629-633.
Meulenberg JJ, Hulst MM, de Meijer EJ, Moonen PL, den Besten A, de Kluyver EP, Wensvoort G, Moormann RJ: 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.1008
Meng XJ: Heterogeneity of porcine reproductive and respiratory syndrome virus: implications for current vaccine efficacy and future vaccine development. Vet Microbiol 2000, 74: 309-329. 10.1016/S0378-1135(00)00196-6
Murtaugh MP, Elam MR, Kakach LT: Comparison of the structural protein coding sequences of the VR-2332 and Lelystad virus strains of the PRRS virus. Arch Virol 1995, 140: 1451-1460. 10.1007/BF01322671
Shi M, Lam TT, Hon CC, Hui RK, Faaberg KS, Wennblom T, Murtaugh MP, Stadejek T, Leung FC: Molecular epidemiology of PRRSV: a phylogenetic perspective. Virus Res 2010, 154: 7-17. 10.1016/j.virusres.2010.08.014
Stadejek T, Oleksiewicz MB, Scherbakov AV, Timina AM, Krabbe JS, Chabros K, Potapchuk D: Definition of subtypes in the European genotype of porcine reproductive and respiratory syndrome virus: nucleocapsid characteristics and geographical distribution in Europe. Arch Virol 2008, 153: 1479-1488. 10.1007/s00705-008-0146-2
Shi M, Lam TT, Hon CC, Murtaugh MP, Davies PR, Hui RK, Li J, Wong LT, Yip CW, Jiang JW, Leung FC: Phylogeny-based evolutionary, demographical, and geographical dissection of North American type 2 porcine reproductive and respiratory syndrome viruses. J Virol 2010, 84: 8700-8711. 10.1128/JVI.02551-09
Thanawongnuwech R, Amonsin A, Tatsanakit A, Damrongwatanapokin S: Genetics and geographical variation of porcine reproductive and respiratory syndrome virus (PRRSV) in Thailand. Vet Microbiol 2004, 101: 9-21. 10.1016/j.vetmic.2004.03.005
Damrongwatanapokin S, Arsayuth K, Kongkrong C, Parchariyanon S, Pinyochon W, Tantaswasdi U: Serological studies and isolation of porcine reproductive and respiratory syndrome (PRRS) virus in Thailand. J Thai Vet Med Assoc 1996, 47: 19-30.
Amonsin A, Kedkovid R, Puranaveja S, Wongyanin P, Suradhat S, Thanawongnuwech R: Comparative analysis of complete nucleotide sequence of porcine reproductive and respiratory syndrome virus (PRRSV) isolates in Thailand (US and EU genotypes). Virol J 2009, 6: 143. 10.1186/1743-422X-6-143
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32: 1792-1797. 10.1093/nar/gkh340
Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics 2005, 21: 260-262. 10.1093/bioinformatics/bth490
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180
Kedkovid R, Nuntawan Na Ayudhya S, Amonsin A, Thanawongnuwech R: NSP2 gene variation of the North American genotype of the Thai PRRSV in central Thailand. Virol J 2010, 7: 340. 10.1186/1743-422X-7-340
Labarque G, Reeth KV, Nauwynck H, Drexler C, Van Gucht S, Pensaert M: Impact of genetic diversity of European-type porcine reproductive and respiratory syndrome virus strains on vaccine efficacy. Vaccine 2004, 22: 4183-4190. 10.1016/j.vaccine.2004.05.008
van Woensel PA, Liefkens K, Demaret S: Effect on viraemia of an American and a European serotype PRRSV vaccine after challenge with European wild-type strains of the virus. Vet Rec 1998, 142: 510-512. 10.1136/vr.142.19.510
The authors would like to thank Drs. Suphasawatt Puranaveja and Roongtham Kedkovid from Chulalongkorn University-Veterinary Diagnosis Laboratory (CU-VDL) for samples and providing information. This work was partially supported by the Strategic Research Theme of Infection and Immunology, The University of Hong Kong.
The authors declare that they have no competing interests.
HMT carried out the advanced phylogenetic studies and drafted manuscript. HMT, MS and CLYW participated in data analysis. SNNA did sampling and sequencing. AA and RT provided sequence data and helped to draft the manuscript. FCCL and HMT conceived the study, participated in its design and coordination. All authors read and approved the final manuscript.
About this article
Cite this article
Tun, H.M., Shi, M., Wong, C.L. et al. Genetic diversity and multiple introductions of porcine reproductive and respiratory syndrome viruses in Thailand. Virol J 8, 164 (2011). https://doi.org/10.1186/1743-422X-8-164
- Multiple Introduction
- Average Genetic Distance
- Swine Industry
- Bayesian Markov Chain Monte Carlo
- ORF5 Gene Sequence