BRIEF RESEARCH REPORT article

Front. Vet. Sci., 11 October 2022

Sec. Veterinary Epidemiology and Economics

Volume 9 - 2022 | https://doi.org/10.3389/fvets.2022.986850

Genetic diversity of imported PRRSV-2 strains, 2005–2020, Hungary

  • 1. Veterinary Medical Research Institute, Budapest, Hungary

  • 2. National Laboratory for Infectious Animal Diseases, Antimicrobial Resistance, Veterinary Public Health and Food Chain Safety, Budapest, Hungary

  • 3. Department of Pharmacology and Toxicology, University of Veterinary Medicine, Budapest, Hungary

  • 4. Veterinary Diagnostic Directorate, National Food Chain Safety Office, Budapest, Hungary

  • 5. Hungarian Association for Porcine Health Management, Budapest, Hungary

  • 6. National PRRS Eradication Committee, Budapest, Hungary

Abstract

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 (, ). 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) (). PRRSVs are enveloped, positive-sense, single-stranded RNA viruses with approximately 15 kb long genome which encodes at least 10 open reading frames ().

Reproductive failure and respiratory problems in affected sows and piglets are commonly seen irrespectively of the virus species (), and comparable virulence has been observed in the reproductive disorder among sows (). However, the respiratory disease is typically more severe in infections caused by PRRSV-2 ().

PRRSV-2 strains have been classified into nine lineages (L1-9) and numerous sub-lineages based on the ORF5 gene (). 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 (). 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 (). 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 ().

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 (). 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 (, , ). 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 (, , ). 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 (, ). 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 ().

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 (). 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 (). Pairwise sequence identities were calculated for a subset of sequences in Mega X () (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 (). 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 (). Alignment of deduced GP5 amino acid sequences was visually inspected and sequences were compared with relevant reference strains, focusing primarily on the antigenic regions () and aa sites contributing to virus rebound or escaping from the host immune system (, ). 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 (). 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.

Table 1

Pig farms#Farm typeSample IDsYear (and lineage)Clinical symptomsMortalityOrigin (and year) of animals positive for PRRSV-2§
SeverityRespiratoryReproduction
ALarge-scale farrowing finishing farmHU12
7703_NEBIH_2013_HU*
2005 (L1)
2013 (L1)
Mild
Mild to moderate
YesNot visible3% during prefattening periodPig import from Slovakia
(2005, 2013)
BLarge-scale fattening farmHU21
PRRSV-2_Hungary_102_2012
2006 (L1)
2012 (L1)
MildYesNot applicable13% during fattening periodCirculation among herds within Hungary
CLarge-scale fattening farm4020_NEBIH_2013_HU*2013 (L1)MildYesNot applicable8% during fattening periodCirculation among herds within Hungary
DLarge-scale farrowing prefattening farm19678_NEBIH_2010_HU*
28663_NEBIH_2011_HU*
27125_NEBIH_2011_HU*
2010 (L1)
2011 (L5)
Mild
None
YesNot visible2% during prefattening periodPig import from Slovakia (2011)
ELarge-scale fattening farm273_NEBIH_2009_HU*2009 (L5)MildYesNot applicable10% during fattening periodNot available
FLarge-scale farrowing finishing farm23090_NEBIH_2010_HU*2010 (L5)MildYesNot visible10% from weaning to slaughterCirculation among herds within Hungary
GLarge-scale fattening farm32708_NEBIH_2011_HU*2011 (L5)NoneYes12.7 weaned piglets/farrowing3.5% during prefattening periodPig import from Slovakia (2011)
HLarge-scale fattening farm50708_NEBIH_2014_HU2014 (L5)NoneYesNot applicable6% during fattening periodPig import from Slovakia (2014)
IBreeding farm11033_NEBIH_2015_HU*2015 (L5)Stillborn pigletsOnly PRRSV infection90% died newborn piglets rate, at normal farrowing timeNo visible effectCirculation among herds within Hungary
JLarge scale fattening farm6576_NEBIH_2015_HU*
42487_NEBIH_2016_HU*
60499_56-60_NEBIH_2017_HU* 35565_96-100_NEBIH_2017_HU*
3158_11-15_NEBIH_2017_HU*
2015 (L5)
2016 (L5)
2017 (L5)
Mild
Mild
Mild
YesNot applicableMortality rate doubles, from 30/month to 75/monthPig import from Czech Republic (2015),
circulation among herds within Hungary (2016, 2017)
KLarge scale fattening farm59074_NEBIH_2015_HU*2015 (L5)NoneYesNot applicable5% during fattening periodPig import from Denmark (2015)
LLarge-scale breeding farm61748_NEBIH_2015_HU*2015 (L5)NoneYesWeaned piglets: 11/farrow farrowing rate: 84.78%11.2% during lactation 2.1% during prefattening and 2% during fattening periodPig import from Belgium (2015)
MLarge-scale fattening farm8727_NEBIH_2016_HU2016 (L5)NoneYesNot applicable6.5% during fattening periodPig import from Slovakia (2016)
NLarge-scale fattening farm25004_NEBIH_2016_HU*2016 (L5)MortalityYesNot applicableMortality rate doublesPig import from Czech Republic (2016)
OLarge-scale breeding farm23920_28_NEBIH_2017_HU*2017 (L5)NoneYesNot visibleMortality rate increased among weaned piglets significantlyCirculation among herds within Hungary
PFattening farm32607_NEBIH_2017_HU2017 (L5)NoneYesNot applicable5.5% during fattening periodPig import from Denmark (2017)
QLarge-scale fattening farm55352_NEBIH_2018_HU
55352_5_NEBIH_2018_HU*
55352_10_NEBIH_2018_HU*
55352_15_NEBIH_2018_HU*
55352_20_NEBIH_2018_HU*
55352_16-20_NEBIH_2018_HU*
2018 (L5)NoneOnly PRRSV infectionNot applicableNo visible effectPig import from Denmark (2018)
RLarge-scale breeding farm64196_2_NEBIH_2020_HU*2020 (L5)NoneYesFarrowing rate 89%, pregnancy rate 91%Liveborn piglets 13.8/farrowCirculation among herds within Hungary
SFattening farm47206_125_NEBIH_2020_HU*
47206_150_NEBIH_2020_HU*
46978_4_NEBIH_2020_HU*
46978_1_NEBIH_2020_HU*
46978_11_NEBIH_2020_HU*
40533_NEBIH_2020_HU*
32023_6-10_NEBIH_2020_HU*
107f_NEBIH_2020_HU*
122f_NEBIH_2020_HU*
156f_NEBIH_2020_HU*
2020 (L5)Stillborn pigletsOnly PRRSV infectionLiveborn/farrow decreased from 14.28 to 12.53Mortality during prefattening period increased from 21 to 30 piglets/monthSeed import from Ireland (2020)
TFattening farm56612_81-85_NEBIH_2020_HU*
56612_40-44_NEBIH_2020_HU*
56612_45-49_NEBIH_2020_HU*
2020 (L5)NoneOnly PRRSV infectionNot applicable5 died piglets out of 4,000Pig import from Germany (2020)

PRRSV-2 associated epidemiological and clinical data, Hungary, 2005–2020.

#

Location of pig farms can be seen in Figure 2.

§

Based on inspection of import licenses and TRACES certificates (provided by the Hungarian veterinary authorities).

*

Determined in this study.

Altogether, 44 partial ORF5 gene sequences were obtained from Hungarian pig holdings over a period of 15 years, including the previously characterized isolates (, , ). 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. (). 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 (), with no other sequences from additional neighboring countries, which emphasize their unique position among European PRRSV-2 field strains.

Figure 1

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 (). 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 (, , ).

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

Lineage L1 PRRSV-2 is the most abundant group worldwide; it is mainly dispersed in USA, Canada and China (). 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) (, ). 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 (). 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 (, ). A more recent extensive study that focused on lineage diversity in USA expanded the representative dataset of L1 sub-lineages (). 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 (, ). In addition, this sub-lineage presumably predominated the viral population from 1990 to 2004 in the USA ().

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 (, ). 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 (). 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) (). 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 (, ).

Figure 3

GP5 is a dominant protein participating in the neutralization process, also it has been acknowledged as a crucial factor contributing to virus rebound (, , , , ). Among the recognized antigenic regions, epitope B and C (aa37-45 and aa52-61) were proved to be neutralizing (, ). 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 (). 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 (, ) 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 (). N30 is another common glycosylation site (, , ) 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 (, , ). 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 (, ). 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 (, ). Hence, it seems a realistic scenario that new virus variants that may cause local epizootics in the future may remain overlooked.

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.

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.

Statements

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.

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.

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.

    WensvoortGTerpstraCPolJMter LaakEABloemraadMde KluyverEPet al. Mystery swine disease in the Netherlands: the isolation of Lelystad virus. Vet Q. (1991) 13:12130. 10.1080/01652176.1991.9694296

  • 2.

    KeffaberK. Reproductive failure of unknown etiology. Am Assoc Swine Pract Newsl. (1989) 1:110.

  • 3.

    BrintonMAGulyaeyaABalasuriyaUBRDunowskaMFaabergKSGoldbergTExpansion 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).

  • 4.

    MeulenbergJJMM. PRRSV the virus. Vet Res. (2000) 31:1121. 10.1051/vetres:2000103

  • 5.

    MaclachlanNJDuboviEJeditors. Fenner's Veterinary Virology. 5th ed. London: Academic Press (2016). p. 602.

  • 6.

    HanKSeoHWParkCKangIYounS-KLeeSYet 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:297305. 10.1016/j.jcpa.2013.11.205

  • 7.

    HanKSeoHWOhYKangIParkCChaeC. 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:3138. 10.1016/j.tvjl.2012.06.035

  • 8.

    Martínez-LoboFJDíez-FuertesFSegalésJGarcía-ArtigaCSimarroICastroJMet 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:5868. 10.1016/j.vetmic.2011.06.025

  • 9.

    ChoiKLeeJParkCJeongJChaeC. 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:31724. 10.1016/j.jcpa.2015.03.002

  • 10.

    ShiMLamTT-YHonC-CMurtaughMPDaviesPRHuiRK-Het al. Phylogeny-based evolutionary, demographical, and geographical dissection of North American type 2 porcine reproductive and respiratory syndrome viruses. J Virol. (2010) 84:870011. 10.1128/JVI.02551-09

  • 11.

    BrarMSShiMMurtaughMPLeungFC-C. Evolutionary diversification of type 2 porcine reproductive and respiratory syndrome virus. J Gen Virol. (2015) 96:157080. 10.1099/vir.0.000104

  • 12.

    StadejekTStankeviciusAMurtaughMPOleksiewiczMB. Molecular evolution of PRRSV in Europe: current state of play. Vet Microbiol. (2013) 165:218. 10.1016/j.vetmic.2013.02.029

  • 13.

    KvisgaardLKHjulsagerCKBrarMSLeungFCCLarsenLE. Genetic dissection of complete genomes of Type 2 PRRS viruses isolated in Denmark over a period of 15 years. Vet Microbiol. (2013) 167:33444. 10.1016/j.vetmic.2013.09.023

  • 14.

    BøtnerAStrandbygaardBSørensenKJHavePMadsenKGMadsenESet al. Appearance of acute PRRS-like symptoms in sow herds after vaccination with a modified live PRRS vaccine. Vet Rec. (1997) 141:4979. 10.1136/vr.141.19.497

  • 15.

    IndikSSchmollFSiposWKleinD. Genetic variability of PRRS virus in Austria: consequences for molecular diagnostics and viral quantification. Vet Microbiol. (2005) 107:1719. 10.1016/j.vetmic.2005.01.024

  • 16.

    Greiser-WilkeIFiebigKDrexlerCgrosse BeilageE. 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:21323. 10.1016/j.vetmic.2009.12.006

  • 17.

    PapatsirosVStylianakiIPapakonstantinouGTsekourasNBitchavaDChristodoulopoulosGet al. Histopathological lesions accompanied with first-time isolation of a PRRSV-2 strain in Greece. Viral Immunol. (2020) 33:16. 10.1089/vim.2020.0087

  • 18.

    SzabóIMolnárTNemesIAbonyiTTerjékZBálintÁet al. eradication on large-scale fattening pig farms in Hungary between 2014 and 2019. Acta Vet Hung. (2019) 67:52942. 10.1556/004.2019.052

  • 19.

    BalkaGWangXOlaszFBálintÁKissIBányaiKet al. Full genome sequence analysis of a wild, non-MLV-related type 2 Hungarian PRRSV variant isolated in Europe. Virus Res. (2015) 200:18. 10.1016/j.virusres.2015.01.014

  • 20.

    DirbákováZMojŽišováGJadudPMojŽ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.

    VilcekSJackovaALeskovaVVlasakovaM. 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. 446.

  • 22.

    BalkaGHornyákÁBálintÁKissIKecskemétiSBakonyiTet al. Genetic diversity of porcine reproductive and respiratory syndrome virus strains circulating in Hungarian swine herds. Vet Microbiol. (2008) 127:12835. 10.1016/j.vetmic.2007.08.001

  • 23.

    BálintÁMolnárTKecskemétiSKulcsárGSoósTSzabóPMet al. Genetic variability of PRRSV vaccine strains used in the National Eradication Programme, Hungary. Vaccines. (2021) 9:849. 10.3390/vaccines9080849

  • 24.

    NemesIMolnárTAbonyiTTerjékZBálintÁSzabóI. Eradication of PRRS from backyard swine herds in Hungary between 2012 and 2018. Acta Vet Hung. (2019) 67:54352. 10.1556/004.2019.053

  • 25.

    SzabóIBognárLMolnárTNemesIBálintÁPRRS. eradication from swine farms in five regions of Hungary. Acta Vet Hung. (2020) 68:25762. 10.1556/004.2020.00043

  • 26.

    StadejekTStankeviciusAStorgaardTOleksiewiczMBBelákSDrewTWet 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:186173. 10.1099/0022-1317-83-8-1861

  • 27.

    JakabSMartonSSzabóIKecskemétiSBálintÁBányaiK. Preliminary observations concerning the integration of amplicon deep sequencing in the PRRS Eradication Program. Hungarian Vet J. (2022) 144:11528.

  • 28.

    LarssonA. AliView: a fast and lightweight alignment viewer and editor for large data sets. Bioinformatics. (2014) 30:32768. 10.1093/bioinformatics/btu531

  • 29.

    KumarSStecherGLiMKnyazCTamuraKMEGAX. Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:15479. 10.1093/molbev/msy096

  • 30.

    MartinDPVarsaniARoumagnacPBothaGMaslamoneySSchwabTet al. RDP5: a computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. (2021) 7:87. 10.1093/ve/veaa087

  • 31.

    TrifinopoulosJNguyenLTvon HaeselerAMinhBQW-IQ-TREE. a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. (2016) 44:W2325. 10.1093/nar/gkw256

  • 32.

    OstrowskiMGaleotaJAJarAMPlattKBOsorioFALopezOJ. Identification of neutralizing and nonneutralizing epitopes in the porcine reproductive and respiratory syndrome virus GP5 ectodomain. J Virol. (2002) 76:424150. 10.1128/JVI.76.9.4241-4250.2002

  • 33.

    de LimaMPattnaikAKFloresEFOsorioFA. 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:41021. 10.1016/j.virol.2006.05.036

  • 34.

    PlagemannPGW. The primary GP5 neutralization epitope of North American isolates of porcine reproductive and respiratory syndrome virus. Vet Immunol Immunopathol. (2004) 102:26375. 10.1016/j.vetimm.2004.09.011

  • 35.

    VashishtKGoldbergTLHusmannRJSchnitzleinWZuckermannFA. 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:474753. 10.1016/j.vaccine.2008.06.047

  • 36.

    PopescuLNTribleBRChenNRowlandRRR. GP5 of porcine reproductive and respiratory syndrome virus (PRRSV) as a target for homologous and broadly neutralizing antibodies. Vet Microbiol. (2017) 209:906. 10.1016/j.vetmic.2017.04.016

  • 37.

    FanBLiuXBaiJZhangTZhangQJiangP. 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:2130. 10.1016/j.virusres.2015.04.015

  • 38.

    ChenNTribleBRKerriganMATianKRowlandRRR. 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:16775. 10.1016/j.meegid.2016.03.002

  • 39.

    JackováAVlasákováMMandelíkRVilčekS. Genetic typing of porcine reproductive and respiratory syndrome virus isolates from central European countries. Acta Virol. (2013) 57:3638.

  • 40.

    KikutiMSanhuezaJVilaltaCPaploskiIADVanderWaalKCorzoCA. 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:110. 10.1371/journal.pone.0259531

  • 41.

    KimS-CJeongC-GParkG-SParkJ-YJeoungH-YShinG-Eet al. Temporal lineage dynamics of the ORF5 gene of porcine reproductive and respiratory syndrome virus in Korea in 2014–2019. Arch Virol. (2021) 166:280315. 10.1007/s00705-021-05169-w

  • 42.

    KwonTYooSJSunwooSYLeeD-UJeSHParkJWet 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:21324. 10.1007/s00705-018-4048-7

  • 43.

    PaploskiIADPamornchainavakulNMakauDNRoviraACorzoCASchroederDCet al. Phylogenetic structure and sequential dominance of sub-lineages of PRRSV type-2 lineage 1 in the United States. Vaccines. (2021) 9:117. 10.3390/vaccines9060608

  • 44.

    ShiMLemeyPSingh BrarMSuchardMAMurtaughMPCarmanSet al. The spread of type 2 porcine reproductive and respiratory syndrome virus (PRRSV) in North America: A phylogeographic approach. Virology. (2013) 447:14654. 10.1016/j.virol.2013.08.028

  • 45.

    YuanSMickelsonDMurtaughMPFaabergKS. Complete genome comparison of porcine reproductive and respiratory syndrome virus parental and attenuated strains. Virus Res. (2001) 74:99110. 10.1016/S0168-1702(00)00250-1

  • 46.

    DelisleBGagnonCALambertM-ÈD'AllaireS. 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:11119. 10.1016/j.meegid.2012.03.015

  • 47.

    StorgaardTOleksiewiczMBøtnerA. Examination of the selective pressures on a live PRRS vaccine virus. Arch Virol. (1999) 144:2389401. 10.1007/s007050050652

  • 48.

    Robinson SR LiJNelsonEAMurtaughMP. Broadly neutralizing antibodies against the rapidly evolving porcine reproductive and respiratory syndrome virus. Virus Res. (2015) 203:5665. 10.1016/j.virusres.2015.03.016

  • 49.

    YoungJEDvorakCMTGrahamSPMurtaughMP. Isolation of porcine reproductive and respiratory syndrome virus GP5-specific, neutralizing monoclonal antibodies from hyperimmune sows. Front Immunol. (2021) 12:110. 10.3389/fimmu.2021.638493

  • 50.

    Zhou YJ YuHTian ZJ LiGXHaoXFYanLPet 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:13644. 10.1016/j.virusres.2009.04.013

  • 51.

    WeiZLinTSunLLiYWangXGaoFet 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:994151. 10.1128/JVI.07067-11

  • 52.

    AnsariIHKwonBOsorioFAPattnaikAK. 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:39944004. 10.1128/JVI.80.8.3994-4004.2006

  • 53.

    LiJTaoSOrlandoRMurtaughMP. N-glycosylation profiling of porcine reproductive and respiratory syndrome virus envelope glycoprotein 5. Virology. (2015) 478:8698. 10.1016/j.virol.2015.02.013

  • 54.

    BálintÁBalkaGHorváthPKecskemétiSDánÁFarsangAet al. Full-length genome sequence analysis of a Hungarian porcine reproductive and respiratory syndrome virus isolated from a pig with severe respiratory disease. Arch Virol. (2015) 160:41722. 10.1007/s00705-014-2265-2

Summary

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

Volume

9 - 2022

Edited by

Igor Paploski, University of Minnesota Twin Cities, United States

Reviewed by

Michael Zeller, Iowa State University, United States; Juan Manuel Sanhueza, Catholic University of Temuco, Chile; Ruwini Rupasinghe, University of California, Davis, United States

Updates

Copyright

*Correspondence: Krisztián Bányai

This article was submitted to Veterinary Epidemiology and Economics, a section of the journal Frontiers in Veterinary Science

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics