Corrigendum: Genetic diversity of imported PRRSV-2 strains, 2005–2020, Hungary
- 1Veterinary Medical Research Institute, Budapest, Hungary
- 2National Laboratory for Infectious Animal Diseases, Antimicrobial Resistance, Veterinary Public Health and Food Chain Safety, Budapest, Hungary
- 3Department of Pharmacology and Toxicology, University of Veterinary Medicine, Budapest, Hungary
- 4Veterinary Diagnostic Directorate, National Food Chain Safety Office, Budapest, Hungary
- 5Hungarian Association for Porcine Health Management, Budapest, Hungary
- 6National PRRS Eradication Committee, Budapest, Hungary
Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) remains sporadic in Europe. In this study, we investigated the molecular epidemiology of PRRSV-2 infections encompassing 15 years in Hungary. Partial (423 bp long) ORF5 sequences (n = 44) from 20 Hungarian pig herds were analyzed. The study strains fell into two genetic lineages, L1 and L5, being L5 strains more prevalent (88.6 vs. 11.4%). Pairwise sequence identities within Hungarian representative PRRSV-2 strains ranged between 84.7 to 100% (nucleotide, nt) and 85 to 100% (amino acid, aa). When compared with reference strains, identity values fell between 87 and 100% (L1, nt 87–91%, aa 87–93%, reference strain IAF-exp91; L5, nt 87–100%, aa 88–100%, reference strain Ingelvac MLV). Epidemiologic examination implied that the majority of L5 strains were imported repeatedly from other European countries where Ingelvac MLV was approved for routine use. The emergence of L1 strains was thought to be associated with a single introduction and subsequent dissemination between pig farms of a large integrator. Results presented here contribute to a better understanding of the epizootiology of PRRSV-2 infections and shed light on the genetic diversity of viral strains in non-endemic countries.
Introduction
Porcine reproductive and respiratory syndrome (PRRS) has become a major cause of economic losses for the global swine industry since its first description in Europe and North America (1, 2). PRRS virus (PRRSV), the causative agent, is divided into two species within the Arteriviridae family, Betaarterivirus suid 1 and Betaarterivirus suid 2 (PRRSV-1 and -2, respectively) (3). PRRSVs are enveloped, positive-sense, single-stranded RNA viruses with approximately 15 kb long genome which encodes at least 10 open reading frames (4).
Reproductive failure and respiratory problems in affected sows and piglets are commonly seen irrespectively of the virus species (5), and comparable virulence has been observed in the reproductive disorder among sows (6). However, the respiratory disease is typically more severe in infections caused by PRRSV-2 (7–9).
PRRSV-2 strains have been classified into nine lineages (L1-9) and numerous sub-lineages based on the ORF5 gene (10). This gene is a highly variable region of the genome making it the main target for evaluating phylogenetic relationships among strains. The majority of genetic lineages have been reported to occur in the Americas and Asia (11). On the contrary, in Europe where the enzootic species is PRRSV-1, the diversity of PRRSV-2 is lower and the circulating European PRRSV-2 strains predominantly belong to the sub-lineage 5.1 within L5 (12–18). This sub-lineage is represented by the PRRSV-2 prototype strain VR2332 and its attenuated form, the Ingelvac MLV (also known as RespPRRS) vaccine virus along with their descendants. Beside the prevalent L5 genotype, a few genotype L1 isolates were identified in Slovakia and Hungary, which originated presumably from Canada or the United States (19–22).
In Hungary, the control of PRRS is limited to disease caused by PRRSV-1 and is supported by routine use of different live, attenuated and/or inactivated vaccines, such as Porcilis PRRS, Unistrain PRRS, Reprocyc PRRS EU, Ingelvac PRRS Flex EU, Suvaxyn PRRS and Progressis (23). PRRSV-2 MLV vaccines (including Ingelvac PRRS MLV) have not been registered to use in Hungary.
PRRSV-2 is not widespread in Hungary, and as a consequence, published data about the genetic diversity of Hungarian PRRSV-2 strains is limited (18, 19, 22). This study reports the ORF5 gene based characterization of all Hungarian origin PRRSV-2 strains in a period encompassing 2005–2020.
Materials and methods
Samples
From 2005 to 2020, serum and organ samples were collected for routine PRRSV monitoring, or from herds suspected to sustain PRRS disease, as described in detail previously (18, 24, 25). Samples were tested within 48 h of arrival to laboratory by PRRSV-specific PCR or ELISA methods. In case of positive test results, the import licenses and TRACES certificates (provided by the Hungarian veterinary authorities) were inspected to determine the origin of imported pigs or the pig movement between Hungarian farms.
The recorded clinical symptoms associated with PRRSV-2 infection ranged from asymptomatic to mild, but in four outbreaks, severe problems such as abortions, stillborn piglets and respiratory disorders among pre-fatteners with high mortality were documented.
Following PRRSV-specific diagnostic testing the remainders of the serum and organ samples were placed at −20°C and then thawed when molecular epidemiologic investigation was initiated.
RNA extraction and ORF5 sequencing
RNA was isolated from samples using QIAmp Viral RNA Mini Kit. The amplification of partial ORF5 was performed by nested RT-PCR assay by means of previously published protocols, which are routinely utilized for diagnostic purposes (22, 26). In brief, the first round included reverse transcription followed by PCR using the Qiagen One-Step RT-PCR Kit and the second round was carried out with Top Taq DNA polymerase to amplify a 432 bp long region of the ORF5. The second-round PCR products were sequenced by Sanger method on ABI PRISM 3100 automatic sequencer, or by semiconductor sequencing technology on Ion Torrent PGM (27).
Sequence analysis
Sequence data were collected from various sources. For the determination of lineages a large ORF5 gene dataset (n = 847), which was kindly provided by Dr. Mang Shi, was utilized as reference (10). Furthermore, 86 PRRSV-2 sequences were collected from GenBank deposited by European countries other than Hungary. A total of 37 sequences were newly determined in this study (see Supplementary material, fasta file) and an additional seven GenBank records from Hungary (accession numbers, EF406336, DQ366650, KM514315, MN150539, MN150535, MN150538, MN150540) were added to the analyses. GenBank accession numbers were retrieved for a subset of Hungarian PRRSV-2 strains, which represented unique sequences (accession numbers, ON939096-ON939121).
Multiple alignments of nucleotide and deduced amino acid sequences (corresponding to nt and aa positions 91-522 and 31-174 relative to the ORF5 gene of strain VR2332) were generated by the MUSCLE algorithm implemented in AliView software (28). Pairwise sequence identities were calculated for a subset of sequences in Mega X (29) (such as identities between all Hungarian sequences and the Quebec reference L1 strain, IAF-exp91, the prototype L5 strain, VR2332, and the lineage L5 vaccine virus strain, Ingelvac PRRS MLV (accession numbers, L40898, EF536003, and AF066183). Genetic recombination events were analyzed with program RDP v5.5 (30). Phylogenetic tree of partial ORF5 sequences was built on the IQ-TREE web-server applying the best fit maximum likelihood model (TVM+F+R8) with 1000 ultrafast bootstrap replications (31). Alignment of deduced GP5 amino acid sequences was visually inspected and sequences were compared with relevant reference strains, focusing primarily on the antigenic regions (32–36) and aa sites contributing to virus rebound or escaping from the host immune system (37, 38). Potential N-glycosylation sites (N-X-S/T) in the GP5 were predicted using the NetGlyc 1.0 Server (https://services.healthtech.dtu.dk/service.php?NetNGlyc-1.0).
Results and discussion
Classification of Hungarian PRRSV-2 strains
In general, the PRRSV positivity rate between 2005 and 2020 was 10-11% in each year (in total 40,785 positives out of 363,721 tested samples). With the progression of the eradication program this rate reduced gradually (18). Between 2005 and 2020, the presence of PRRSV-2 was confirmed in 11 out of 19 counties involving a total of 20 pig farms (Table 1). PRRSV-2 was detected in nearly all years, although the spatiotemporal distribution of the cases was uneven (Supplementary Figure 1). Following the launch of the National PRRS Eradication Program in 2014, the more intense surveillance activity resulted in an increased detection rate of PRRSV-2.
Altogether, 44 partial ORF5 gene sequences were obtained from Hungarian pig holdings over a period of 15 years, including the previously characterized isolates (18, 19, 22). The Hungarian sequences shared 84.7–100% nucleotide and 85–100% amino acid identities with each other. By comparing our samples to the IAF-exp91 strain and to the Ingelvac PRRS MLV vaccine virus, 87–91% and 87–100% nucleotide, 87–93% and 88–100% amino acid identities, respectively, were seen. Sequence analyses revealed neither indel mutations, nor recombination events.
Phylogenetic tree reconstruction was performed in order classify the newly determined PRRSV-2 sequences and to verify the relationships of the previously reported records. The nucleotide sequences of this study and the reference dataset of lineages were aligned along with ORF5 sequences from other European countries (Supplementary Table 2). In the dataset 71, 6, 3, 3, 1, 1, and 1 sequences originated from Denmark, Germany, Poland, Austria, Spain, Slovakia and Lithuania, respectively. The strain classification was performed according to Shi et al. (10). Our analyses indicated that all Hungarian PRRSV-2 sequences could be unequivocally classified; i.e. the majority (39/44, 88.6%) typed as L5, while the remainder (5/44, 11.4%) typed as L1 (Figure 1). The ORF5 based phylogenetic tree showed that most of the Hungarian lineage L5 strains are closely related to Danish MLV-related sequences across Europe. The Hungarian L1 variants formed a distinct monophyletic cluster including a Slovakian strain, 36M (39), with no other sequences from additional neighboring countries, which emphasize their unique position among European PRRSV-2 field strains.
Figure 1. (A) Phylogenetic analysis based on global ORF5 sequence dataset including European sequences (n = 977) of PRRSV-2. Maximum likelihood method with TVM+F+R8 model and 1,000 ultrafast bootstrap were utilized. (B) Lineage 5 and a (C) part of lineage 1 sub-trees were expanded for further presentation. The different geographic origin of the sequences were marked with distinct colors. Red circles () denote the partial ORF5 sequences determined in this study; each dot represents an individual sequence. Prototype strain VR2332 () and the vaccine strain Ingelvac MLV () are indicated separately.
Molecular epidemiology
L5 is a predominant PRRSV-2 lineage worldwide along with L1, L9 and L8. The main geographic locations of lineage L5 strains include Canada, China, USA, Mexico, Korea and Denmark (11). Within lineage L5, it is of great importance to discriminate the Ingelvac PRRS MLV or its descendants from the parental, wild-type VR2332 strain. To date, a particular strain is considered vaccine-related if its ORF5 sequence shares lower than 5% sequence distance from the ORF5 sequence of Ingelvac MLV vaccine virus (10, 11, 40).
Concerning the Hungarian L5 PRRSV-2 strains (39 samples from 17 pig farms), their partial ORF5 sequences shared 95-100% nucleotide identity with the vaccine virus suggesting that they were vaccine derived viruses. Additionally, as confirmed by observations of field veterinarians, several PRRSV-2 strains were detected after the settling of fattening pigs from abroad—e.g., Belgium, Czech Republic, Germany, Slovakia or Denmark—or could be connected in some way to herd(s) abroad (Figure 2; Table 1). In this regard, it is important to emphasize that Ingelvac PRRS MLV is registered to use in those respective countries, but not in Hungary. PRRSV-2 sequence data are currently scarce in Europe (except Denmark). Yet, sequence and epidemiologic data (together with available import licenses and TRACES certificates) indicated that the imported L5 PRRSV-2 isolates could have originated mainly from Denmark and Slovakia. Although, based on the Hungarian legislation, import solely has been allowed from PRRSV-free herds since 2017, in various cases the animals were reported to arrive from collection stations. In collection stations, pigs from PRRS-infected and -free herds may have been kept in the same air space, a prerequisite that could have favored viral spread between infected and naïve animals. Due to varying regulation of PRRSV-2 vaccination in Europe the Ingelvac MLV vaccine strain could be periodically imported to PRRSV-2 free countries by trading of live animals. Such incidents impose significant risk to Hungarian herds currently declared PRRSV-free, since the PRRS eradication program will come to a completion in 2022. So far, several regulations were made to reduce the risk of reinfection at PRRS free farms of Hungary; the most significant achievement might be that only officially determined PRRS-free stocks are admitted to the country from 2017 onward.
Figure 2. Spatial distribution of all pig farms that were infected with PRRSV-2 from 2005 to 2020. The location of the farms were designated by different colored markers according to the detected lineage of PRRSV-2. Origin of imported fattening stocks were also displayed (for further details the reader is referred to Table 1).
Lineage L1 PRRSV-2 is the most abundant group worldwide; it is mainly dispersed in USA, Canada and China (11). In Hungary, five lineage L1 sequences (from four pig farms) were identified; all five sequences were linked to a single recognized event of introduction dating back to 2005 to a farm of a large integrator and spread between other locations of this integrator between 2006 and 2013 (index farm, A; other farms with subsequent transmission events, B, C, D, and G; Figure 2; Table 1) (19, 22). Epidemiological data suggested that the most likely country of origin for the first recorded L1 strain in Hungary was Slovakia (Figure 2; Table 1). By the end of 2014, the year when the national eradication program in Hungary was launched, PRRSV-2 was eliminated in all farms where L1 strains were previously in circulation (Supplementary Figure 1). Since then no additional L1 strains were identified. In Europe, this particular lineage is represented only by our few characterized sequences and a single Slovakian sample. The phylogenetic analysis further confirmed the hypothesis that these PRRSV-2 strains are wild-type and originated from the USA or Canada (19). The Hungarian and Slovakian lineage L1 PRRSV-2 strains display a unique cluster within Europe and could constitute a country-specific group of L1, a finding also recognized in Korea (41, 42). A more recent extensive study that focused on lineage diversity in USA expanded the representative dataset of L1 sub-lineages (43). According to results of BLAST search, Hungarian L1 sequences shared the greatest sequence identity with a newly defined sub-lineage, L1D-alpha (93%). L1D-alpha contains progenitor isolates and it was predicted to have emerged before 1990 (43, 44). In addition, this sub-lineage presumably predominated the viral population from 1990 to 2004 in the USA (43).
Analysis of GP5 amino acid sequences
In order to evaluate the diversity of GP5, the non-redundant deduced amino acid sequences were selected from our dataset, resulting in 32 unique variants for subsequent genetic analyses (Figure 3). The Hungarian PRRSV-2 samples showed a relatively high sequence identity within L5 and L1 (proportion of conservative aa sites: 122/144 and 139/144) according to the amino acid alignment of partial GP5 (a protein fragment encompassing aa31-174). Certain aa residues such as aa56, aa57, aa102, and aa104 have already been identified in association with escape mechanisms from homologous neutralization (36, 37). In our case, this finding is particularly relevant for the vaccine related L5 samples. Some of the characterized sequences (such as 25004, 47206, 46978, and 40533) that originated from pig farms with high mortality rate and increased number of stillborn piglets carried mutations compared to Ingelvac MLV, V102L/F and G104E, which may represent immune escape variants of the vaccine virus. However, these mutations appear in samples of symptom-free pig farms and, for example, not all Danish L5 isolates manifesting clinical signs showed substitutions at these positions (13). Thus, our understanding is limited whether these residues play a role in the evolution of immune escape variants and in the mechanism of virus rebound of the Ingelvac MLV vaccine. Another substantial residue, aa151, is considered to be involved in the attenuation process of VR2332 strain (R → G) (45). We identified numerous substitutions (G151K/V/R) in our L5 sequences and half of these reverted to the amino acid of parental strain VR2332 (Figure 3). In fact, high mutation rate and positive selective pressure were found at this position (46, 47).
Figure 3. Multiple alignment of the partial GP5 amino acid sequences of all Hungarian isolates of PRRSV-2 from 2005 to 2020 together with reference strain IAF-exp91 and VR2332 (The reader is referred to Supplementary Figure 1 for the colored signs). Epitope sites and highly variables regions are represented by blank boxes and gray boxes, respectively. The site of attenuation of VR2332 is indicated in the red box.
GP5 is a dominant protein participating in the neutralization process, also it has been acknowledged as a crucial factor contributing to virus rebound (32, 36, 38, 48, 49). Among the recognized antigenic regions, epitope B and C (aa37-45 and aa52-61) were proved to be neutralizing (32–35, 50). Epitope B was completely conserved (37SHLQLIYNL45) among all Hungarian PRRSV-2 while epitope C, which spans the HVR-2, was variable at position aa57-59. Analyzing the T cell epitopes (aa117-131 and aa149-163) they carried a few substitutions at different amino acid positions.
N-glycosylation of GP5 is a significant mechanism contributing to the infectivity and the antigenicity of PRRSVs and affects the immune response of the host (51–53). Potential N-glycosylation sites were observed at several positions in the GP5 of Hungarian sequences; these include N32, N33, N34, N35, N44 and N51. Two glycosylation sites corresponding to amino acid N44 and N51 seem to be essential for PRRSV infection, therefore these modifications are usually highly conserved in reference (51, 52) as well as in Hungarian PRRSV-2 strains. We found six different N-glycosylation patterns (A-F) composed of three or four asparagines (Supplementary Table 1). The most abundant pattern A (33-44-51) and B (34-44-51) seemed to be widely distributed in the USA (46). N30 is another common glycosylation site (13, 46, 50) which, however, could not be examined from the partial ORF5 sequences of Hungarian PRRSV-2 strains. The N-glycosylation pattern of the full length GP5 regarding Ingelvac MLV vaccine virus strain is 30-33-44-51. Apart from the N30 the Hungarian L5 sequences carry several different changes compared to the vaccine virus, which could imply a trend to evade the host immune response.
Conclusions
In this study, classification, genetic diversity and distribution of circulating PRRSV-2 strains have been investigated in Hungary, a European country where PRRSV-2 is not enzootic. A dataset containing partial PRRSV-2 ORF5 sequences as well as information about PRRSV-2 positive pig farms were collated and analyzed. Despite the low sample size, the data presented here are thought to give a true overall picture of the epizootiology of PRRSV-2 infections recorded in Hungary between 2005 and 2020, with over 360,000 swine origin samples being tested for PRRSV in this period.
In our study, PRRSV strain surveillance was based on sequencing of the partial ORF5 gene. Although the analysis of partial ORF5 alone may hide some important features of viral evolution, this genomic region is still preferably used in molecular epidemiologic investigations (40, 41, 43). Furthermore, this single-gene based strain characterization approach appeared to be an economically reasonable choice for a pioneering national PRRSV eradication program. Sequencing of the structural region or the whole genome of PRRSV strains had not been a primary aim of Hungarian authorities when the PRRS eradication program was launched, however, sequencing the whole genome in some cases nicely complemented the current methodology (19, 54). Thus, as massively parallel sequencing technologies become part of routine diagnostic laboratory workflows and the amount of data these technologies generate is easier to analyze, it is expected that such new approaches will replace old ones over time. Ideally, these technologies will permit whole genome sequencing and bioinformatics pipelines at reasonable costs to be performed that will be an added value when compared to current technologies and protocols used in routine viral strain characterization.
PRRSV-2 is not a major cause of porcine reproductive and respiratory syndrome in Europe and with few exceptions PRRSV-2 lineages circulating in this continent have been known only from isolated outbreaks. This epidemiologic landscape markedly influences the number of viral sequences reported from different countries; thus, one need to keep in mind that conclusions of this study concerning the country of origin of isolates may have some limitations, even if the route of animal transport and available sequence information mutually increase the plausibility of our findings. Our study indicate that key features about the epizootiology and evolution of PRRSV-2 could be missing in Europe as the virus falls out the focus of intensive surveillance (12, 13). Hence, it seems a realistic scenario that new virus variants that may cause local epizootics in the future may remain overlooked.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
SJ, EK, SM, KB, ÁB, IN, and IS: collection and analysis of data. SJ, KB, ÁB, IN, and IS: interpretation of results. SJ, KB, ÁB, and IS: writing the draft manuscript. All authors have agreed to the published version of the manuscript.
Funding
This work was supported by the National Laboratory for Infectious Animal Diseases, Antimicrobial Resistance, Veterinary Public Health and Food Chain Safety, RRF-2.3.1-21-2022-00001. Prepared with the professional support of the Doctoral Student Scholarship Program of the Co-operative Doctoral Program of the Ministry of Innovation and Technology financed from the National Research, Development and Innovation Fund.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.986850/full#supplementary-material
References
1. Wensvoort G, Terpstra C, Pol JM, ter Laak EA, Bloemraad M, de Kluyver EP, et al. Mystery swine disease in the Netherlands: the isolation of Lelystad virus. Vet Q. (1991) 13:121–30. doi: 10.1080/01652176.1991.9694296
3. Brinton MA, Gulyaeya A, Balasuriya UBR, Dunowska M, Faaberg KS, Goldberg T, Expansion of the Rank Structure of the Family Arteriviridae Renaming Its Taxa. (2017). Available online at: https://talk.ictvonline.org/taxonomy/p/taxonomy-history?taxnode_id=20186087 (accessed May 2022).
5. Maclachlan NJ, Dubovi EJ, editors. Fenner's Veterinary Virology. 5th ed. London: Academic Press (2016). p. 602.
6. Han K, Seo HW, Park C, Kang I, Youn S-K, Lee SY, et al. Comparative virulence of reproductive diseases caused by type 1 (European-like) and type 2 (North American-like) porcine reproductive and respiratory syndrome virus in experimentally infected pregnant gilts. J Comp Pathol. (2014) 150:297–305. doi: 10.1016/j.jcpa.2013.11.205
7. Han K, Seo HW, Oh Y, Kang I, Park C, Chae C. Comparison of the virulence of European and North American genotypes of porcine reproductive and respiratory syndrome virus in experimentally infected pigs. Vet J. (2013) 195:313–8. doi: 10.1016/j.tvjl.2012.06.035
8. Martínez-Lobo FJ, Díez-Fuertes F, Segalés J, García-Artiga C, Simarro I, Castro JM, et al. Comparative pathogenicity of type 1 and type 2 isolates of porcine reproductive and respiratory syndrome virus (PRRSV) in a young pig infection model. Vet Microbiol. (2011) 154:58–68. doi: 10.1016/j.vetmic.2011.06.025
9. Choi K, Lee J, Park C, Jeong J, Chae C. Comparison of the pathogenesis of single or dual infections with type 1 and type 2 porcine reproductive and respiratory syndrome virus. J Comp Pathol. (2015) 152:317–24. doi: 10.1016/j.jcpa.2015.03.002
10. Shi M, Lam TT-Y, Hon C-C, Murtaugh MP, Davies PR, Hui RK-H, et al. Phylogeny-based evolutionary, demographical, and geographical dissection of North American type 2 porcine reproductive and respiratory syndrome viruses. J Virol. (2010) 84:8700–11. doi: 10.1128/JVI.02551-09
11. Brar MS, Shi M, Murtaugh MP, Leung FC-C. Evolutionary diversification of type 2 porcine reproductive and respiratory syndrome virus. J Gen Virol. (2015) 96:1570–80. doi: 10.1099/vir.0.000104
12. Stadejek T, Stankevicius A, Murtaugh MP, Oleksiewicz MB. Molecular evolution of PRRSV in Europe: current state of play. Vet Microbiol. (2013) 165:21–8. doi: 10.1016/j.vetmic.2013.02.029
13. Kvisgaard LK, Hjulsager CK, Brar MS, Leung FCC, Larsen LE. Genetic dissection of complete genomes of Type 2 PRRS viruses isolated in Denmark over a period of 15 years. Vet Microbiol. (2013) 167:334–44. doi: 10.1016/j.vetmic.2013.09.023
14. Bøtner A, Strandbygaard B, Sørensen KJ, Have P, Madsen KG, Madsen ES, et al. Appearance of acute PRRS-like symptoms in sow herds after vaccination with a modified live PRRS vaccine. Vet Rec. (1997) 141:497–9. doi: 10.1136/vr.141.19.497
15. Indik S, Schmoll F, Sipos W, Klein D. Genetic variability of PRRS virus in Austria: consequences for molecular diagnostics and viral quantification. Vet Microbiol. (2005) 107:171–9. doi: 10.1016/j.vetmic.2005.01.024
16. Greiser-Wilke I, Fiebig K, Drexler C, grosse Beilage E. Genetic diversity of Porcine reproductive and respiratory syndrome virus (PRRSV) in selected herds in a pig-dense region of North-Western Germany. Vet Microbiol. (2010) 143:213–23. doi: 10.1016/j.vetmic.2009.12.006
17. Papatsiros V, Stylianaki I, Papakonstantinou G, Tsekouras N, Bitchava D, Christodoulopoulos G, et al. Histopathological lesions accompanied with first-time isolation of a PRRSV-2 strain in Greece. Viral Immunol. (2020) 33:1–6. doi: 10.1089/vim.2020.0087
18. Szabó I, Molnár T, Nemes I, Abonyi T, Terjék Z, Bálint Á, et al. eradication on large-scale fattening pig farms in Hungary between 2014 and 2019. Acta Vet Hung. (2019) 67:529–42. doi: 10.1556/004.2019.052
19. Balka G, Wang X, Olasz F, Bálint Á, Kiss I, Bányai K, et al. Full genome sequence analysis of a wild, non-MLV-related type 2 Hungarian PRRSV variant isolated in Europe. Virus Res. (2015) 200:1–8. doi: 10.1016/j.virusres.2015.01.014
20. Dirbáková Z, MojŽišová G, Jadud P, MojŽiš M. Genetic diversity and evolution of PRRS virus strains circulated in pig herds during last 15 years in Slovakia. The 2011 International Porcine Reproductive and Respiratory Syndrome Symposium. Chicago, IL, USA. (2011). Available online at: https://www.researchgate.net/publication/271837642 (accessed October 1, 2021).
21. Vilcek S, Jackova A, Leskova V, Vlasakova M. Genetic Analysis of PRRSV Isolates From Pig farms in Slovakia. EuroPRRS2011. Understanding and combating PRRS in Europe. Novi Sad: Scientific Veterinary Institute Novi Sad (2011). p. 44–6.
22. Balka G, Hornyák Á, Bálint Á, Kiss I, Kecskeméti S, Bakonyi T, et al. Genetic diversity of porcine reproductive and respiratory syndrome virus strains circulating in Hungarian swine herds. Vet Microbiol. (2008) 127:128–35. doi: 10.1016/j.vetmic.2007.08.001
23. Bálint Á, Molnár T, Kecskeméti S, Kulcsár G, Soós T, Szabó PM, et al. Genetic variability of PRRSV vaccine strains used in the National Eradication Programme, Hungary. Vaccines. (2021) 9:849. doi: 10.3390/vaccines9080849
24. Nemes I, Molnár T, Abonyi T, Terjék Z, Bálint Á, Szabó I. Eradication of PRRS from backyard swine herds in Hungary between 2012 and 2018. Acta Vet Hung. (2019) 67:543–52. doi: 10.1556/004.2019.053
25. Szabó I, Bognár L, Molnár T, Nemes I, Bálint Á, PRRS. eradication from swine farms in five regions of Hungary. Acta Vet Hung. (2020) 68:257–62. doi: 10.1556/004.2020.00043
26. Stadejek T, Stankevicius A, Storgaard T, Oleksiewicz MB, Belák S, Drew TW, et al. Identification of radically different variants of porcine reproductive and respiratory syndrome virus in Eastern Europe: towards a common ancestor for European and American viruses. J Gen Virol. (2002) 83:1861–73. doi: 10.1099/0022-1317-83-8-1861
27. Jakab S, Marton S, Szabó I, Kecskeméti S, Bálint Á, Bányai K. Preliminary observations concerning the integration of amplicon deep sequencing in the PRRS Eradication Program. Hungarian Vet J. (2022) 144:115–28.
28. Larsson A. AliView: a fast and lightweight alignment viewer and editor for large data sets. Bioinformatics. (2014) 30:3276–8. doi: 10.1093/bioinformatics/btu531
29. Kumar S, Stecher G, Li M, Knyaz C, Tamura K, MEGA X. Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:1547–9. doi: 10.1093/molbev/msy096
30. Martin DP, Varsani A, Roumagnac P, Botha G, Maslamoney S, Schwab T, et al. RDP5: a computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. (2021) 7:87. doi: 10.1093/ve/veaa087
31. Trifinopoulos J, Nguyen LT, von Haeseler A, Minh BQ, W-IQ-TREE. a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. (2016) 44:W232–5. doi: 10.1093/nar/gkw256
32. Ostrowski M, Galeota JA, Jar AM, Platt KB, Osorio FA, Lopez OJ. Identification of neutralizing and nonneutralizing epitopes in the porcine reproductive and respiratory syndrome virus GP5 ectodomain. J Virol. (2002) 76:4241–50. doi: 10.1128/JVI.76.9.4241-4250.2002
33. de Lima M, Pattnaik AK, Flores EF, Osorio FA. Serologic marker candidates identified among B-cell linear epitopes of Nsp2 and structural proteins of a North American strain of porcine reproductive and respiratory syndrome virus. Virology. (2006) 353:410–21. doi: 10.1016/j.virol.2006.05.036
34. Plagemann PGW. The primary GP5 neutralization epitope of North American isolates of porcine reproductive and respiratory syndrome virus. Vet Immunol Immunopathol. (2004) 102:263–75. doi: 10.1016/j.vetimm.2004.09.011
35. Vashisht K, Goldberg TL, Husmann RJ, Schnitzlein W, Zuckermann FA. Identification of immunodominant T-cell epitopes present in glycoprotein 5 of the North American genotype of porcine reproductive and respiratory syndrome virus. Vaccine. (2008) 26:4747–53. doi: 10.1016/j.vaccine.2008.06.047
36. Popescu LN, Trible BR, Chen N, Rowland RRR. GP5 of porcine reproductive and respiratory syndrome virus (PRRSV) as a target for homologous and broadly neutralizing antibodies. Vet Microbiol. (2017) 209:90–6. doi: 10.1016/j.vetmic.2017.04.016
37. Fan B, Liu X, Bai J, Zhang T, Zhang Q, Jiang P. The amino acid residues at 102 and 104 in GP5 of porcine reproductive and respiratory syndrome virus regulate viral neutralization susceptibility to the porcine serum neutralizing antibody. Virus Res. (2015) 204:21–30. doi: 10.1016/j.virusres.2015.04.015
38. Chen N, Trible BR, Kerrigan MA, Tian K, Rowland RRR. ORF5 of porcine reproductive and respiratory syndrome virus (PRRSV) is a target of diversifying selection as infection progresses from acute infection to virus rebound. Infect Genet Evol. (2016) 40:167–75. doi: 10.1016/j.meegid.2016.03.002
39. Jacková A, Vlasáková M, Mandelík R, Vilček S. Genetic typing of porcine reproductive and respiratory syndrome virus isolates from central European countries. Acta Virol. (2013) 57:363–8.
40. Kikuti M, Sanhueza J, Vilalta C, Paploski IAD, VanderWaal K, Corzo CA. Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) genetic diversity and occurrence of wild type and vaccine-like strains in the United States swine industry. PLoS ONE. (2021) 16:1–10. doi: 10.1371/journal.pone.0259531
41. Kim S-C, Jeong C-G, Park G-S, Park J-Y, Jeoung H-Y, Shin G-E, et al. Temporal lineage dynamics of the ORF5 gene of porcine reproductive and respiratory syndrome virus in Korea in 2014–2019. Arch Virol. (2021) 166:2803–15. doi: 10.1007/s00705-021-05169-w
42. Kwon T, Yoo SJ, Sunwoo SY, Lee D-U, Je SH, Park JW, et al. Independent evolution of porcine reproductive and respiratory syndrome virus 2 with genetic heterogeneity in antigenic regions of structural proteins in Korea. Arch Virol. (2019) 164:213–24. doi: 10.1007/s00705-018-4048-7
43. Paploski IAD, Pamornchainavakul N, Makau DN, Rovira A, Corzo CA, Schroeder DC, et al. Phylogenetic structure and sequential dominance of sub-lineages of PRRSV type-2 lineage 1 in the United States. Vaccines. (2021) 9:1–17. doi: 10.3390/vaccines9060608
44. Shi M, Lemey P, Singh Brar M, Suchard MA, Murtaugh MP, Carman S, et al. The spread of type 2 porcine reproductive and respiratory syndrome virus (PRRSV) in North America: A phylogeographic approach. Virology. (2013) 447:146–54. doi: 10.1016/j.virol.2013.08.028
45. Yuan S, Mickelson D, Murtaugh MP, Faaberg KS. Complete genome comparison of porcine reproductive and respiratory syndrome virus parental and attenuated strains. Virus Res. (2001) 74:99–110. doi: 10.1016/S0168-1702(00)00250-1
46. Delisle B, Gagnon CA, Lambert M-È, D'Allaire S. Porcine reproductive and respiratory syndrome virus diversity of Eastern Canada swine herds in a large sequence dataset reveals two hypervariable regions under positive selection. Infect Genet Evol. (2012) 12:1111–9. doi: 10.1016/j.meegid.2012.03.015
47. Storgaard T, Oleksiewicz M, Bøtner A. Examination of the selective pressures on a live PRRS vaccine virus. Arch Virol. (1999) 144:2389–401. doi: 10.1007/s007050050652
48. Robinson SR Li J, Nelson EA, Murtaugh MP. Broadly neutralizing antibodies against the rapidly evolving porcine reproductive and respiratory syndrome virus. Virus Res. (2015) 203:56–65. doi: 10.1016/j.virusres.2015.03.016
49. Young JE, Dvorak CMT, Graham SP, Murtaugh MP. Isolation of porcine reproductive and respiratory syndrome virus GP5-specific, neutralizing monoclonal antibodies from hyperimmune sows. Front Immunol. (2021) 12:1–10. doi: 10.3389/fimmu.2021.638493
50. Zhou YJ Yu H, Tian ZJ Li GX, Hao XF, Yan LP, et al. Genetic diversity of the ORF5 gene of porcine reproductive and respiratory syndrome virus isolates in China from 2006 to 2008. Virus Res. (2009) 144:136–44. doi: 10.1016/j.virusres.2009.04.013
51. Wei Z, Lin T, Sun L, Li Y, Wang X, Gao F, et al. N-Linked glycosylation of GP5 of porcine reproductive and respiratory syndrome virus is critically important for virus replication in vivo. J Virol. (2012) 86:9941–51. doi: 10.1128/JVI.07067-11
52. Ansari IH, Kwon B, Osorio FA, Pattnaik AK. Influence of N-Linked glycosylation of porcine reproductive and respiratory syndrome virus GP5 on virus infectivity, antigenicity, and ability to induce neutralizing antibodies. J Virol. (2006) 80:3994–4004. doi: 10.1128/JVI.80.8.3994-4004.2006
53. Li J, Tao S, Orlando R, Murtaugh MP. N-glycosylation profiling of porcine reproductive and respiratory syndrome virus envelope glycoprotein 5. Virology. (2015) 478:86–98. doi: 10.1016/j.virol.2015.02.013
Keywords: Betaarterivirus suid 2, molecular epidemiology, phylogenetic analysis, ORF5, vaccine virus
Citation: Jakab S, Kaszab E, Marton S, Bányai K, Bálint Á, Nemes I and Szabó I (2022) Genetic diversity of imported PRRSV-2 strains, 2005–2020, Hungary. Front. Vet. Sci. 9:986850. doi: 10.3389/fvets.2022.986850
Received: 05 July 2022; Accepted: 23 September 2022;
Published: 11 October 2022.
Edited by:
Igor Paploski, University of Minnesota Twin Cities, United StatesReviewed by:
Michael Zeller, Iowa State University, United StatesJuan Manuel Sanhueza, Catholic University of Temuco, Chile
Ruwini Rupasinghe, University of California, Davis, United States
Copyright © 2022 Jakab, Kaszab, Marton, Bányai, Bálint, Nemes and Szabó. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Krisztián Bányai, bkrota@hotmail.com