Characterization of two Austrian porcine reproductive and respiratory syndrome virus (PRRSV) field isolates reveals relationship to East Asian strains
Veterinary Research volume 47, Article number: 17 (2016)
Porcine reproductive and respiratory syndrome virus (PRRSV) causes major problems for the swine industry worldwide. Due to Austria’s central location in Europe, a large number of animals are transported through the country. However, little is known about current PRRSV strains and epidemiology. We determined full-length genome sequences of two Austrian field isolates (AUT13-883 and AUT14-440) from recent PRRSV outbreaks and of a related German isolate (GER09-613). Phylogenetic analysis revealed that the strains belong to European genotype 1 subtype 1 and form a cluster together with a South Korean strain. Remarkably, AUT14-440 infected the simian cell line MARC-145 without prior adaptation. In addition, this isolate showed exceptional deletions in nonstructural protein 2, in the overlapping region of glycoprotein 3 and 4 and in the 3′ untranslated region. Both Austrian isolates caused similar lung lesions but only pigs infected with AUT14-440 developed clear clinical signs of infection. Taken together, the genetic and biological characterization of two novel Austrian PRRSV field isolates revealed similarities to East Asian strains. This stresses the necessity for a more detailed analysis of current PRRSV strains in Europe beyond the determination of short ORF5 and ORF7 sequences.
Porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most economically important pathogens in the swine industry worldwide [1, 2]. It is the etiological agent of porcine reproductive and respiratory syndrome (PRRS), which is characterized by respiratory disorders as well as by growth retardation in growing pigs and reproductive failure in late gestation sows [3, 4]. The disease emerged in the late eighties in North America, with the first outbreaks in Europe recorded in 1990 .
The causative agent, PRRSV, was first isolated in 1991 in the Netherlands . This strain, Lelystad virus (LV), is regarded as the prototype strain of European PRRSV type 1 (PRRSV-1), whereas VR2332 represents the North American PRRSV type 2 (PRRSV-2). The two genotypes share only about 60% identity at the nucleotide level . Due to high mutation and recombination rates , variability is also high within the genotypes, especially in type 1. Therefore, three subtypes have been proposed for PRRSV-1 based the size of the nucleocapsid protein (N): pan-European subtype 1 and Eastern European subtypes 2 and 3 .
PRRSV is a small, enveloped virus with a single-stranded positive-sense RNA genome that belongs to the family Arteriviridae, order Nidovirales . The 5′-capped and 3′-polyadenylated genome of PRRSV is about 15 kb in length and contains ten open reading frames (ORF) [9, 10]. ORF1a and 1b constitute over 75% of the viral genome and encode two polyproteins, which are cleaved into at least 14 nonstructural proteins (nsp) that are responsible for genome replication and transcription . ORF2–4 encode the minor structural proteins, including glycoprotein (GP) 2, 3 and 4, which form a hetero-trimer that is believed to be involved in virus entry . The major structural proteins GP5, membrane protein (M) and N protein are encoded by ORF5–7. Nucleotide sequences of ORF5 (603–606 bp) and ORF7 (372–387 bp) are widely used for phylogenetic studies. However, the analysis is limited as these short genomic sequences might be subject to recombination and immunological selection pressure . Phylogenetic analysis based on ORF5 and ORF7 may also lead to different subtyping of PRRSV-1 strains . For these reasons, the use of further methodologies for genetic subtyping has been suggested, such as whole genome sequencing [6, 14]. The lack of published full-length sequences is exemplified by the situation in Austria, where the limited data on current strains consists solely of ORF5 or ORF7 sequences [15–17]. This gap in knowledge complicates the identification of the source of outbreaks, which is crucial as the most important PRRSV problems in Austria occur as a consequence of (re-) introduction into formerly free herds. An awareness of recent isolates and genomic changes that possibly lead to lower protection is also important in terms of assessing the effectiveness of vaccines.
The aim of this study was to determine the etiological agents of two PRRSV-1 outbreaks in Austria and to characterize them genetically and biologically. The results reveal that the isolated PRRSV-1 subtype 1 field strains (AUT13-883 and AUT14-440) have genetic and phylogenetic similarities to East Asian strains.
Materials and methods
Viruses and cells
Wild-type PRRSV strains AUT13-883 and AUT14-440 were isolated from sera of two farms in Upper Austria in 2013 and 2014. A cell culture adapted PRRSV-1 subtype 1 strain (GER09-613) that was isolated from the field in Germany in 2009 was also included in the analysis (kindly provided by L. Haas, TiHo Hannover). MARC-145 cells (CCLV RIE 277 ) were obtained from the Friedrich-Loeffler-Institute in Germany. Preparation of porcine alveolar macrophages (PAM) was essentially performed as described previously . Briefly, lungs of 6- to 12 week-old PRRSV-free, healthy pigs were washed 2–3 times with PBS (pH 7.2). Aliquots of lavage fluid were pooled and centrifuged for 5 min at 1000 g. The cells were washed in PBS and repelleted twice. For storage, PAM were resuspended in FCS supplemented with 10% dimethyl sulfoxide and frozen at −150 °C. MARC-145 cells were cultured in Dulbeccos’s modified Eagle’s medium (Life Technologies, Waltham, USA) supplemented with 10% fetal calf serum (Bio&Sell, Nürnberg, Germany) and antibiotics (100 U/mL penicillin and 100 µg/mL streptomycin). Two additional antibiotics (enrofloxacin and kanamycin) were added for the culture of PAM. All cells were maintained at 37 °C and 5% CO2 and observed daily for cytopathic effects.
Serum samples were collected from two medium-sized piglet producers in Upper Austria in November 2013 (farm 883) and March 2014 (farm 440). Both farms were experiencing PRRSV-related problems such as pneumonia, wasting, variance in growth, conjunctivitis and coughing in rearing piglets, as well as reproductive disorders such as mummification and abortions in sows. Whereas one farm (883) had experienced minor problems related to PRRSV over a period of 2 years, the other (440) was formerly free of PRRSV and had losses of up to 50% in one farrowing badge. Serum samples from both farms had been tested positive for PRRSV by qRT-PCR by the Animal Health Service Upper Austria. Neither of the farms had used vaccines against PRRSV.
Virus isolation and titration
For virus isolation, 50 µL serum of diseased piglets from farms 440 and 883 were used to inoculate 5 × 106 PAM on a six-well format. After 2 days, the cell culture supernatant was clarified by centrifugation (5 min, 3000 g) and passaged on PAM and MARC-145 cells. Virus infection was generally analysed 2 days post-infection by immunofluorescence staining. Virus titres of supernatants were confirmed on PAM from a single animal to ensure comparability and expressed as tissue culture infectious dose of 50% (TCID50). In the animal trial a viral stock from the 6th to 8th passage was used for inoculation.
Detection of PRRSV-positive cells
Cells were fixed with 1:2 methanol-acetone for 2 min at room temperature and air-dried. Anti-PRRSV-N-protein monoclonal antibodies (clone P10/b1)  (kindly provided by A. Saalmüller, Vienna) were employed for the detection of infected cells. Goat anti-mouse conjugated with Cy3 (Dianova/Jackson, Hamburg, Germany) was used as a secondary antibody.
Determination and analysis of full-length genome sequences
Total RNA was extracted from PAM 5 days after inoculation and from serum samples using the RNeasy Mini Kit and the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. RNA was eluted in 30 µL distilled water and either used immediately for RT-PCR or stored at −80 °C for subsequent analysis.
RT-PCR was carried out using the LongAmpKit (NEB, Ipswich, USA) or the One Step RT-PCR Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Seven primer pairs that were highly conserved in PRRSV-1 strains were designed (available upon request) based on published sequences available in NCBI GenBank. The resulting PCR amplicons were overlapping and covered the whole genome.
The DNA fragments were subjected to gel electrophoresis, purified by the peqGOLD Gel Extraction Kit (Peqlab, Erlangen, Germany) and sequenced by commercial laboratories (Microsynth Austria, Vienna, Austria and Eurofins Genomics, Ebersberg, Germany). Primer selection for sequencing was based on available PRRSV-1 sequences. The determined full-length genome sequences of GER09-613, AUT14-440 and AUT13-883 were submitted to GenBank (KT344816, KT334375 and KT326148).
Initial phylogenetic analysis was carried out with NCBI’s Basic Local Alignment Search Tool for nucleotides (BLASTn). The two closest neighbours of the isolates presented in this study were determined for ORF5 and ORF7. Pairwise comparison and identity calculations were carried out with CLC Main Workbench 7.6 (CLCBIO, Aarhus, Denmark). Alignments and phylogenetic trees were generated with the software CLC Sequence Viewer 7.6 (CLCBIO, Aarhus, Denmark) with bootstrap values based on 1000 replicates. All PRRSV-1 strains with full genome sequences deposited in GenBank were used to construct the phylogenetic trees. For the trees based on ORF5 and ORF7, the two closest neighbours—as determined by NCBI BLASTn—were added. The PRRSV-2 prototype VR-2332 was used as an out-group. Recombination analysis was performed with the recombination analysis tool (RAT)  using the full genome alignment of PRRSV-1 strains.
Twenty-seven-week-old conventional pigs (#1–20) were obtained from a PRRSV-negative herd and housed in a biosafety level 2 facility. They had been vaccinated against PCV-2 and M. hyopneumoniae in their third week of life. To confirm their PRRSV status nasal swabs taken at the time of arrival were analysed by qRT-PCR and serum samples were tested for antibodies against PRRSV with the commercial IDEXX X3 ELISA (IDEXX laboratories, Westbrook, USA). The animal experiments were approved by the ethics committee of the University of Veterinary Medicine, Vienna and the Federal Ministry of Science, Research and Economy (BMWF-68.205/0196-WF/V/3b/2014).
Tenfold dilution series of sera were prepared and transferred to naïve PAM. For each time point and each animal, two replicates were performed. Three days later cells were fixed and stained and the TCID50 was calculated.
Quantitative reverse transcription-PCR (qRT-PCR)
To detect viremia, viral RNA was extracted from 140 µL serum with QIAamp Viral RNA Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. RNA was eluted in 60 µL distilled water and 1 µL was directly used for amplification with the KAPA™ SYBR® FAST One-Step qRT-PCR Kit (Peqlab, Erlangen, Germany) on an ABI 7500 cycler (Applied Biosystems, Foster City, USA). Published primers Pesch PLS (5′-ATGGCCAGCCAGTCAATC-3′) and Pesch PLR (5′-TCGCCCTAATTGAATAGGTG-3′) were used for amplification . A recombinant cDNA clone of a PRRSV-1 strain was purified using the QIAGEN Plasmid Midi Kit (Qiagen, Hilden, Germany) and spectrophotometrically quantified. The copy number of recombinant plasmids was calculated with the formula: N (molecules per µL) = [C (DNA concentration in µg/µL)/K (fragment size in bp)] × 185.5 × 1013 (factors derived from DNA weight, volume and the Avogadro constant). To obtain a standard curve, a tenfold dilution series of cDNA of the PRRSV-1 clone was included in the qRT-PCR setup. Cycling conditions were 42 °C 5:00, 95 °C 5:00 and 42 cycles of 95 °C 0:03, 60 °C 1:00 and 80 °C 0:33 (fluorescence detection step) followed by a dissociation step (95 °C 0:15, 60 °C 1:00, 95 °C 0:15). The numbers of genome copies were calculated with the 7500 System SDS Software (Applied Biosystems, Foster City, USA) based on the standard curve and projected to 1 mL serum by multiplication by 428.4.
To detect viral shedding, nasal swabs were submerged in 2 mL of sterile saline (0.7% NaCl) and thoroughly shaken for 30 s. 200 µL of supernatant were used for viral RNA extraction with the High Pure Viral RNA Kit (Roche Diagnostics, Vienna, Austria) according to the manufacturer’s instructions. RNA was eluted in 50 µL of the kit’s elution buffer and 13.6 µL were directly used for amplification with the LightCycler® RNA Amplification Kit SYBR Green I with a LightCycler® 480-II System (Roche Diagnostics, Vienna, Austria). The nested RT-PCR from Pesch  was adapted for detection with SYBR Green I®. Reverse transcription and activation was performed according to the manufacturer’s instructions, followed by 32 cycles of 60 °C 0:10, 72 °C 0:22 and 95 °C 0:10. The nested amplification step was carried out with LightCycler® FastStart DNA Master(plus) SYBR Green I (Roche Diagnostics, Vienna, Austria) using 2 µL of the first PCR reaction. Cycling conditions were 25 cycles of 60 °C 0:10, 72 °C 0:20 and 95 °C 0:20. The fluorescence detection step was carried out at 84 °C. The PCR was followed by a dissociation step. Results were considered positive when a clear amplification curve was observed and the melting temperature of the amplicons (calculated by LightCycler® 480 System Software) was within the expected range.
Serum samples were tested for antibodies against PRRSV with the commercial IDEXX X3 ELISA (IDEXX laboratories, Westbrook, USA) according to the manufacturer’s instructions. Ratios of optical densities of sample and positive control (P/PK ratios) were calculated and considered positive when >0.4.
Data were analysed with the SPSS statistics software (SPSS Ltd, Quarry Bay, Hong Kong). Oneway ANOVA with Posthoc LSD test was applied for the analysis of all results except for those from the ELISA, which were analysed with the t test with the test value 0. A p value < 0.05 was considered statistically significant.
Two PRRSV-1 outbreaks in 2013 (farm 883) and 2014 (farm 440) were brought to our attention by the responsible veterinarians. Case 883 was characterized by mild respiratory disorders and a low virus load that persisted over 2 years. In case 440, suckling piglet mortality reached 50%, which is severe for PRRSV outbreaks in Austria, and abortions were observed. Both farms had been found to be PRRSV-1 positive by qRT-PCR before we obtained any samples. To investigate the PRRSV strains responsible for the outbreaks, we isolated two viruses from the serum of diseased piglets after serial inoculation of PAM (termed AUT13-883 and AUT14-440). MARC-145 cells were also susceptible for AUT14-440 without adaption but were not infected by AUT13-883.
The sequences of the ORF5 and 7 genes revealed that both isolates belong to PRRSV-1 subtype 1  but the results for ORF5 and ORF7 were not consistent. In the phylogenetic tree based on ORF5, AUT14-440 clustered together with German and South Korean PRRSV-1 isolates (among them KNU-07) (Additional file 1). AUT13-883 represented a different branch of the tree and was most closely related to German strain EU-2a from 1992  and to Austrian strain 2888 . In contrast, when ORF7 sequences were used for phylogenetic analysis, the two Austrian isolates formed a group together with Korean strain KNU-07 , a Croatian strain  and German strains  (Additional file 2).
Since the results of the phylogenetic trees based on ORF5 and ORF7 were inconclusive, complete genome sequences for AUT13-883 and AUT14-440 were determined (GenBank: KT326148 and KT334375). A phylogenetic tree based on the full genome sequences of the isolates presented in this study, 38 published PRRSV-1 strains and PRRSV-2 prototype VR-2332 (Figure 1) confirmed the relationship between the two Austrian isolates, GER09-613 and KNU-07. All four strains cluster together, with no other strain in the same branch. AUT13-883 and GER09-613 are the most closely related strains.
Distinct deletions in the genome of AUT14-440
Assembly of the overlapping sequences resulted in complete genomes consisting of 15022 nucleotides (nt) for AUT14-440 and 15095 nt for GER09-613 and AUT13-883, excluding the 3′ poly(A) tails. Detailed comparison to the European prototype LV revealed marked differences in the region of nsp2. Both Austrian virus isolates and the German strain GER09-613 possess a shorter ORF1a due to deletions in nsp2 (Figure 2). The difference of three nucleotides in ORF1a between AUT13-883 (7188 nt) and LV (7191 nt) results in the deletion of proline 182 in nsp2 (Figure 3A, dotted line box). The deletion is also present in the South Korean strain KNU-07 and has been described for strains from Hong Kong. GER09-613 and AUT14-440 (ORF1a lengths of 7188 and 7152 nt) share a deletion of three nt, resulting in the loss of glutamine 183 in nsp2, one amino acid (aa) downstream of proline 182 (Figure 3A, dotted line box). In comparison to LV, AUT14-440 has further deletions in nsp2 at aa positions 320–323, aa 359–364 and aa 699–700, all in hypervariable regions of nsp2 where deletions are frequently found (Figure 3A, solid line boxes). ORF3 and ORF4 of AUT14-440 were also truncated by 36 nt when compared to LV, resulting in a deletion of 12 aa in the overlapping region of GP3 and GP4 (Figures 3B and C solid line boxes). Shorter deletions (up to 8 aa) in this area have been described for several strains, including HK5 and the Chinese isolate BJEU06-1. AUT14-440 also shows a single nucleotide deletion at position 87 in the 3′ untranslated region (UTR) (Figure 3D, solid line box). All deletions in virus isolate AUT14-440 were also found when RNA from the original serum sample was used for RT-PCR and PCR products were directly sequenced.
Pairwise nucleotide comparison of the complete genome resulted in identities of 87.7, 87.9 and 88.6% with LV for AUT14-440, AUT13-883 and GER09-613. Identities between the isolates were 86.4% (GER09-613 to AUT14-440), 85.4% (AUT14-440 to AUT13-883) and 89.9% (AUT13-883 to GER09-613) (Figure 2).
Determination of virulence
Genetic and phylogenetic analysis showed a striking relationship between the Austrian isolates, the German strain GER09-613 and the South Korean strain KNU-07. The reported severity of the clinical symptoms associated with AUT13-883 and AUT14-440 differed markedly, while the clinical symptoms of GER09-613 have not been described. However, the particular clinical signs in an affected herd provide no more than an indication of the pathogenicity of a PRRSV isolate, as other pathogens may be involved. To assess the virulence and to test Koch’s postulates, an animal trial was performed with the Austrian isolates AUT13-883 and AUT14-440 and the German isolate GER09-613, which clusters in the same branch of the phylogenetic tree.
PRRSV-seronegative pigs were distributed equally into four groups according to their weight and sex. Each group of five animals was housed in a separate compartment to avoid cross-contamination. After 1 week of adaptation (day 7 to day 0), pigs from groups GER09-613, AUT13-883 and AUT14-440 were inoculated intranasally with 3 mL medium (1.5 mL in each nostril) containing 1 × 105 TCID50 of the respective strain. Pigs from the control group received the same amount of medium without virus. Clinical signs and rectal temperature were monitored daily. A scoring system was applied that rates liveliness, dyspnoea, coughing, nasal and ocular discharge, conjunctivitis and cyanosis with scores from 0 (physiological) to 3 (severe clinical signs). Daily scores for each animal were added to obtain a value for the overall health status and the mean clinical scores for the four groups were calculated (Figure 4A). All pigs infected with AUT14-440 showed mild to moderate dyspnoea on three or more study days. Other clinical signs that occurred frequently in all animals in this group were apathy and conjunctivitis. Two waves of clinical signs could be distinguished: the first wave occurred around day 4 post infection (pi) and the second, with more prominent signs of illness, took place around day 12 pi. In the other infection groups, only one (AUT13-883) or two (GER09-613) animals showed mild dyspnoea for no more than 2 days in a row. Pigs in the control group showed no relevant clinical signs in the trial, except for animal #4 that had to be euthanized on study day 10 due to a severe lameness in the hind legs. There was no significant difference in rectal temperature between the four groups.
The pigs were weighed at −7, 0, 3, 7, 10 and 13 days post infection (dpi). For each week of the trial (study days −7 to 0, 0 to 7, 7 to 13) daily weight gain was calculated and is shown for each group as a box plot with maximum and minimum values represented as whiskers (Figure 4B). From the first week to the second week of the trial, all groups showed an average increase of daily weight gain between 0.08 and 0.16 kg. From the second to the third week, the control group and group GER09-613 showed a mean increase of daily weight gain of 0.05 and 0.11 kg, whereas in the groups infected with the Austrian field isolates daily weight gain did not increase (mean values of −0.07 kg for AUT14-440 and −0.02 kg for AUT13-883). Despite the apparent differences between the groups, the figures are not statistically significant due to the small group size and the high variation within the groups (p = 0.055 between GER09-613 and AUT14-440 for the third week).
Gross pathology and histopathology
At 13 and 14 dpi, pigs were euthanized and necropsied with special emphasis on the lungs. No other organs than lungs showed pathological lesions. Samples from each lung lobe were fixed in formalin, embedded in paraffin, stained with hematoxylin and eosin and examined histologically. Lesions were hypertrophy and hyperplasia of pneumocytes, septal infiltration by mononuclear cells, perivascular and intraalveolar accumulation of inflammatory cells and necrotic debris. Depending on their occurrence in different lung lobes and on the severity of the lesions, a total score per lung was calculated, applying a PRRSV specific scoring system that has been described elsewhere . The mean values for the Austrian isolates were 52.4 ± 20.9 for AUT14-440 and 53.6 ± 9.4 for AUT13-883 (Figure 5A). In these two groups all types of lesions mentioned above were present in mild to severe manifestations with no visible differences between the groups. The results differed significantly (p < 0.01) from both the control group with a mean score of 5.6 ± 3.0 and group GER09-613, which mainly showed minor accumulations of inflammatory cells (mean histological score of 15 ± 6.0). The scores for each group are displayed as a boxplot in Figure 5B.
Viremia and virus shedding
To exclude the possibility of cross-contamination between the groups, 50 µL of serum taken 3 dpi from one animal per group (#10, 14, 19) were inoculated onto naïve PAM and virus was re-isolated. RT-PCR and subsequent sequencing of the 700 bp PCR amplicon covering ORF4 confirmed the identity of the viral strains. Compared to the parental strain, two point mutations were detected in each strain. One of them was the same for all three strains and resulted in an amino-acid change from aspartic acid to asparagine at position 33 in GP4. In addition melting temperatures of amplicons from positive nasal swabs confirmed that no cross-contamination had occurred (data not shown).
Viral nucleic acids in serum could be detected with qRT-PCR as early as day 3 pi in all animals from group AUT14-440 and in two animals from group AUT13-883 (Figure 6). In the GER09-613 group, the first virus-positive animals were detected 10 dpi. On day 13 pi all animals from this group except one (#8) tested positive for viral RNA; animal #8 remained negative for qRT-PCR as well as for titration on PAM (data not shown) on all sampling days. The number of genomic copies in the serum of animals infected with AUT14-440 declined towards the end of the trial. In the AUT13-883 group, the onset of viremia differed between the pigs. The maximum number of genome copies per mL serum was up to 100 times higher in the AUT14-440 group than in the other groups. On days 3 and 7 pi, the difference between group AUT14-440 and all other groups was significant (p < 0.01), whereas on day 10 pi only the difference to the control group and the AUT13-883 group was significant (p < 0.05). Serum titres on PAM were also highest in the AUT14-440 group (up to 2.3 × 103 TCID50/mL, data not shown). Viral RNA or infectious virus could be detected in the serum of all infected animals, except for #8, with qRT-PCR and titration on PAM on one or more sampling days. The animals from the control group gave negative results with both assays throughout the trial.
To examine virus shedding from the upper respiratory tract, nasal swabs were taken at −7, 3, 7, 10 and 13 dpi and analysed with qRT-PCR for PRRSV-specific nucleic acids (data not shown). All animals from infected groups were virus-positive 3 dpi. Group GER09-613 tested completely negative on day 7 pi and one pig (#8) remained negative until the end of the trial, while the other animals tested positive on day 10 and 13 pi. Virus shedding could be determined for four pigs each from groups AUT14-440 and AUT13-883 on study days 7 and 10, whereas only two animals from group AUT14-440 and all animals from group AUT13-883 were positive on the last sampling day, day 13. The nasal swabs from the control group tested negative for PRRSV-specific nucleic acids at every time point.
Antibody titres against PRRSV were determined at the beginning of the trial to ensure the PRRSV-free status of the animals. At this time point all pigs tested negative for N-specific antibodies with the commercial IDEXX X3 ELISA (Figure 7). Pigs were further screened for seroconversion 7, 10 and 13 dpi. First detection of PRRSV-specific antibodies was possible 7 dpi in the sera of two pigs from groups AUT14-440 and GER09-613. On days 10 and 13 pi all infected animals but one (#19) tested positive in ELISA. The clearest antibody response with the highest P/PK ratios was seen in the AUT14-440 group. The animals from the control group were negative for PRRSV-specific antibodies at every time point.
Austria is surrounded by seven neighbouring countries and lies on the border between Western and Eastern Europe. Due to its central location it has become a major European corridor for travel and the movement of goods, including live pigs. In fact, one of the isolates (AUT440-14) originated from a farm directly adjacent to the major west-east transit road. The reintroduction of PRRSV into virus-free herds is the most urgent yet least understood problem for the swine industry in Austria. Epidemiological evidence is often insufficient as only few ORF5 and ORF7 sequences from Austria are accessible at GenBank, most of them originating from only three publications [15–17]. In none of these cases were sequence data correlated to the biological properties of the isolates.
In this study, the initial question about the phylogenetic background of two field strains not only led to complete genome sequences of two recent Austrian field isolates but also to their characterization in vivo. (Phylo-) genetic information could be interpreted in the light of the phenotypic properties of the isolates. This is especially important because genetic relationships do not give evidence for the virulence of PRRSV strains . To go beyond the descriptive determination of selected sequences, it is necessary to address the pathogen’s properties by carefully characterizing the isolate.
The initial virus isolation in this study was performed on primary PAM, which are believed to be the main target cells of PRRSV in vivo and are widely used for PRRSV-1 isolation and propagation . The much more convenient simian cell line MARC-145 can usually not be used for unadapted PRRSV-1 virus strains, although it is readily infected by PRRSV-2 without adaptation [6, 28]. Interestingly, it also supported growth of the first passage of AUT14-440, a PRRSV-1 isolate. This observation encouraged us to investigate the phylogenetic classification of the strain.
Virus isolation is rarely performed in routine PRRSV diagnostics and sequences are normally obtained directly from clinical samples. Generally, only ORF5 or ORF7 sequences are determined. However, there are doubts about the validity of phylogenetic trees based solely on ORF5 or ORF7 sequences because of the recombination potential and high rate of nucleotide substitution shown by these genes [6, 14]. As a matter of fact, in the phylogenetic trees we determined based on ORF5 and 7 the relationship between the two Austrian isolates and closely related strains remained unclear as AUT13-883 grouped differently in the two trees. The difference might relate to recombination events, which are believed to have a major impact on PRRSV diversity [6, 29]. No evidence for this phenomenon was found by recombination analysis, probably due to the small number of full PRRSV-1 genomes available in GenBank (data not shown). This finding underlines the need for more full-length genome sequencing. The phylogenetic tree based on full genomes clearly showed relatedness between the two Austrian isolates, the German strain GER09-613 and the South Korean isolate KNU-07. However, it has to be considered that the four strains cluster together in a phylogenetic tree despite not being closely related (identities between 85.4 and 90.2%) and that only 38 PRRSV-1 strains were available for comparison, most of them from earlier than 2010. The picture might be less clear if the full genome sequences of more (recent) strains were available for inclusion in the analysis. At present we can only speculate whether the South Korean strain was imported to Europe or whether the cluster is the result of convergent evolution. It is also possible that all four strains originate from older strains from Germany (e.g., BH_95_10-08_EU) that show close relationships in ORF5 and ORF7 but for which complete sequences are not available. Detailed epidemiological statements will not be possible without a broader range of full-length PRRSV genomes.
Nevertheless, comparing not only a small part of the 15 kb PRRSV genome but full sequences revealed remarkable characteristics of the strains and provided indications that the four strains have a common origin. AUT13-883 and KNU-07 have a proline deletion in nsp2 that has only been noted in several strains from Hong Kong . AUT14-440 and GER09-613 also share a single aa deletion in nsp2, a single position after the proline. Additionally, on either side of the 12 aa deletion in GP3 in AUT14-440 the four strains share areas of aa sequences that are very similar to one another but clearly differ from those of other PRRSV-1 strains, providing further indications of a possible relationship.
A closer look at the full genome sequences revealed additional deletions in strain AUT14-440. In nsp2, three more deletions of 2–6 aa support the conclusion from the phylogenetic analysis that the isolate is not as closely related as the three other strains. The 12 aa deletion in the overlapping region of GP3 and GP4 is unique to this strain. Other published isolates that carry a deletion in this highly variable area only miss 1–8 aa compared to LV and the concomitant occurrence of deletions in nsp2 and ORF3/4 has only been described for a few Chinese strains [31–33]. The first strains carrying a deletion in the overlapping region of ORF3 and ORF4 were isolated in Denmark as long ago as 1992, shortly after the discovery of PRRSV. In subsequent evaluations, deletion mutants were reported to evolve more quickly than non-deleted viruses . This might be due to immunological pressure, since neutralizing antibodies against this region in GP4 have been described in vitro for LV  and two Belgian strains . This region and the corresponding region in GP3 were also shown to be immunogenic for other isolates, including a Danish deletion mutant with an 8 aa deletion [34, 37]. Therefore, it is very likely that the highly variable region in the overlap of GP3 and GP4 is under negative selection pressure, causing deletions or mutations to occur . This is in line with the existence of very similar aa sequences next to the deletion in the Austrian isolates and their nearest relatives GER09-613 and KNU-07, indicating that the four strains evolved in a similar direction, probably due to immunological pressure. The 12 aa deletion in AUT14-440 might prevent the binding of GP4-specific antibodies and thereby result in an evolutionary advantage. This idea will be tested in further experiments.
Another interesting aspect of the deletion is the fact that the heterotrimeric minor glycoproteins GP2, GP3 and GP4 are abundant on the viral surface and probably determine cell tropism . A deletion of 12 aa in two of these proteins might lead to an altered cell tropism, as seen in this study for AUT14-440, which is, despite its type 1 genotype, capable of infecting MARC-145 cells. The single nt deletion in the 3′UTR of AUT14-440 has not been described in other PRRSV-1 strains, although there is a report of a deletion at a different position within the 3′UTR in highly pathogenic Chinese strains . As the 3′UTR is essential for PRRSV replication , deletions in this area might have an implication for the replication efficacy of the virus.
AUT14-440 differed from the other Austrian isolate AUT13-883 and the German isolate GER09-613 also in other aspects. It was the only strain to cause clear respiratory distress in pigs in the animal trial. Animals infected with AUT14-440 ceased to increase daily weight gain in the second week after infection, highlighting the economic relevance of the strain. Although pigs from the other groups remained relatively healthy throughout the trial, moderate to severe lung lesions were found in pigs infected with AUT13-883, which also showed varying developments of daily weight gain. It is conceivable that clinical signs would have been more prominent in all infection groups if a less natural route of infection, e.g., intramuscular, or higher viral titres had been chosen.
Viremia was already detected 3 dpi in pigs infected with the Austrian strains, which is in accordance with other experimental infections [42, 43], whereas the GER09-613 group showed the first viral RNA and titres in serum 10 dpi. A likely explanation for the late onset and the lack of clinical signs is that strain GER09-613 was cell-culture adapted on MARC-145 cells and hence potentially attenuated and thus had to re-adapt to the host. Future studies will be necessary to address the sequence differences between inoculated GER09-613 and the re-isolated virus.
In the present study, high amounts of viral RNA in blood were associated with severity of clinical signs, with titres of AUT14-440 exceeding those of the other strains by up to 100-fold. In contrast, the amounts of viral RNA in nasal swabs were lowest for pigs infected with AUT14-440, indicating that the amount of virus shedding via the respiratory route is not necessarily correlated with viremia. The detection of virus shedding 3 dpi is probably an artefact of intranasal infection, as pigs from group GER09-613 do not show any viral RNA in nasal swabs 7 dpi, which corresponds to the late onset of viremia in this group.
The early detection of non-neutralizing N-specific antibodies is in accordance with previous studies [42, 43]. The IDEXX X3 ELISA is not quantitative but the high P/PK ratios of group AUT14-440 are in agreement with the high number of viral genome copies in the blood of animals infected with this strain.
In summary, we have determined the complete genome sequence of two recent Austrian field isolates of PRRSV and interpret the nucleic acid sequences and the phylogeny in the light of the strains’ in vivo characteristics. Both isolates cluster with a German field isolate from 2009 and the South Korean KNU-07 strain in whole genome phylogenetic analysis. One of the field isolates, AUT14-440, caused clinical signs in an animal experiment, grew on MARC-145 cells and showed exceptional deletions. These well characterized isolates represent an excellent basis for further studies on the implications of the described molecular properties on virus entry, replication and pathogenicity. It remains unclear whether the similarities to East Asian strains are unique to European strains because of the lack of full genome sequences of strains currently present in the field. As PRRSV strains are subject to constant evolutionary pressure, both from vaccination and from instruments of modern herd management, it is of the utmost importance to keep up with the increasing diversity of PRRSV by studying current isolates.
Neumann EJ, Kliebenstein JB, Johnson CD, Mabry JW, Bush EJ, Seitzinger AH, Green AL, Zimmerman JJ (2005) Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States. J Am Vet Med Assoc 227:385–392
Nieuwenhuis N, Duinhof TF, van Nes A (2012) Economic analysis of outbreaks of porcine reproductive and respiratory syndrome virus in nine sow herds. Vet Rec 170:225
Goyal SM (1993) Porcine reproductive and respiratory syndrome. J Vet Diagn Invest 5:656–664
Wensvoort G, Terpstra C, Pol JM, ter Laak EA, Bloemraad M, de Kluyver EP, Kragten C, van Buiten L, den Besten A, Wagenaar F, Broekhuijsen JM, Moonen PL, Zetstra T, de Boer EA, Tibben HJ, de Jong MF, van’t Veld P, Groenland GJ, van Gennep JA, Voets MT, Verheijden JH, Braamskamp J (1991) Mystery swine disease in The Netherlands: the isolation of Lelystad virus. Vet Q 13:121–130
Nelsen CJ, Murtaugh MP, Faaberg KS (1999) Porcine reproductive and respiratory syndrome virus comparison: divergent evolution on two continents. J Virol 73:270–280
Murtaugh MP, Stadejek T, Abrahante JE, Lam TT, Leung FC (2010) The ever-expanding diversity of porcine reproductive and respiratory syndrome virus. Virus Res 154:18–30
Stadejek T, Oleksiewicz MB, Scherbakov AV, Timina AM, Krabbe JS, Chabros K, Potapchuk D (2008) Definition of subtypes in the European genotype of porcine reproductive and respiratory syndrome virus: nucleocapsid characteristics and geographical distribution in Europe. Arch Virol 153:1479–1488
Fields BN, Knipe DM, Howley PM (2013) Fields virology. Wolters Kluwer Health/Lippincott Williams & Wilkins, Philadelphia
Firth AE, Zevenhoven-Dobbe JC, Wills NM, Go YY, Balasuriya UB, Atkins JF, Snijder EJ, Posthuma CC (2011) Discovery of a small arterivirus gene that overlaps the GP5 coding sequence and is important for virus production. J Gen Virol 92:1097–1106
Johnson CR, Griggs TF, Gnanandarajah J, Murtaugh MP (2011) Novel structural protein in porcine reproductive and respiratory syndrome virus encoded by an alternative ORF5 present in all arteriviruses. J Gen Virol 92:1107–1116
Fang Y, Snijder EJ (2010) The PRRSV replicase: exploring the multifunctionality of an intriguing set of nonstructural proteins. Virus Res 154:61–76
Wissink EH, Kroese MV, van Wijk HA, Rijsewijk FA, Meulenberg JJ, Rottier PJ (2005) Envelope protein requirements for the assembly of infectious virions of porcine reproductive and respiratory syndrome virus. J Virol 79:12495–12506
Shi M, Lam TT, Hon CC, Hui RK, Faaberg KS, Wennblom T, Murtaugh MP, Stadejek T, Leung FC (2010) Molecular epidemiology of PRRSV: a phylogenetic perspective. Virus Res 154:7–17
Stadejek T, Stankevicius A, Murtaugh MP, Oleksiewicz MB (2013) Molecular evolution of PRRSV in Europe: current state of play. Vet Microbiol 165:21–28
Indik S, Schmoll F, Sipos W, Klein D (2005) Genetic variability of PRRS virus in Austria: consequences for molecular diagnostics and viral quantification. Vet Microbiol 107:171–178
Balka G, Hornyak A, Dan A, Ladinig A, Biksi I, Rusvai M (2010) PriProET based melting point analyses on PRRSV positive field samples. Mol Cell Probes 24:411–414
Jackova A, Vlasakova M, Mandelik R, Vilcek S (2013) Genetic typing of porcine reproductive and respiratory syndrome virus isolates from central European countries. Acta Virol 57:363–368
Kim HS, Kwang J, Yoon IJ, Joo HS, Frey ML (1993) Enhanced replication of porcine reproductive and respiratory syndrome (PRRS) virus in a homogeneous subpopulation of MA-104 cell line. Arch Virol 133:477–483
Zeman D, Neiger R, Yaeger M, Nelson E, Benfield D, Leslie-Steen P, Thomson J, Miskimins D, Daly R, Minehart M (1993) Laboratory investigation of PRRS virus infection in three swine herds. J Vet Diagn Invest 5:522–528
Weiland E, Wieczorek-Krohmer M, Kohl D, Conzelmann KK, Weiland F (1999) Monoclonal antibodies to the GP5 of porcine reproductive and respiratory syndrome virus are more effective in virus neutralization than monoclonal antibodies to the GP4. Vet Microbiol 66:171–186
Etherington GJ, Dicks J, Roberts IN (2005) Recombination analysis tool (RAT): a program for the high-throughput detection of recombination. Bioinformatics 21:278–281
Fetzer C, Pesch S, Ohlinger VF (2006) High risk of false positive results in a widely used diagnostic test for detection of the porcine reproductive and respiratory syndrome virus (PRRSV). Vet Microbiol 115:21–31
Martinez-Lobo FJ, Diez-Fuertes F, Simarro I, Castro JM, Prieto C (2011) Porcine reproductive and respiratory syndrome virus isolates differ in their susceptibility to neutralization. Vaccine 29:6928–6940
Nam E, Park CK, Kim SH, Joo YS, Yeo SG, Lee C (2009) Complete genomic characterization of a European type 1 porcine reproductive and respiratory syndrome virus isolate in Korea. Arch Virol 154:629–638
Novosel D, Stadejek T (2011) Genetic diversity of European genotype PRRSV in Central and Eastern Europe—an update. GenBank. http://ncbi.nlm.nih.gov/nuccore/JQ043208.1
Wernike K, Bonilauri P, Dauber M, Errington J, LeBlanc N, Revilla-Fernandez S, Hjulsager C, Isaksson M, Stadejek T, Beer M, Hoffmann B (2012) Porcine reproductive and respiratory syndrome virus: interlaboratory ring trial to evaluate real-time reverse transcription polymerase chain reaction detection methods. J Vet Diagn Invest 24:855–866
Balka G, Ladinig A, Ritzmann M, Saalmuller A, Gerner W, Kaser T, Jakab C, Rusvai M, Weissenbock H (2013) Immunohistochemical characterization of type II pneumocyte proliferation after challenge with type I porcine reproductive and respiratory syndrome virus. J Comp Pathol 149:322–330
OIE (2012) Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. Office international des épizooties, Paris
van Vugt JJ, Storgaard T, Oleksiewicz MB, Botner A (2001) High frequency RNA recombination in porcine reproductive and respiratory syndrome virus occurs preferentially between parental sequences with high similarity. J Gen Virol 82:2615–2620
Brar MS, Shi M, Hui RK, Leung FC (2014) Genomic evolution of porcine reproductive and respiratory syndrome virus (PRRSV) isolates revealed by deep sequencing. PLoS One 9:e88807
Chen N, Cao Z, Yu X, Deng X, Zhao T, Wang L, Liu Q, Li X, Tian K (2011) Emergence of novel European genotype porcine reproductive and respiratory syndrome virus in mainland China. J Gen Virol 92:880–892
Li B, Gao S, Zhou T, Liu X, Lu H, Feng F (2014) Complete genome sequence of European genotype porcine reproductive and respiratory syndrome virus strain LNEU12 in Northern China. Genome Announc 2:e00957–01014
Zhou Z, Liu Q, Hu D, Yue X, Yu X, Zhang Q, Gu X, Ni J, Li X, Zhai X, Tian K (2013) Complete genome sequence of a European genotype porcine reproductive and respiratory syndrome virus in china. Genome Announc 1:e00175–00213
Oleksiewicz MB, Botner A, Toft P, Grubbe T, Nielsen J, Kamstrup S, Storgaard T (2000) Emergence of porcine reproductive and respiratory syndrome virus deletion mutants: correlation with the porcine antibody response to a hypervariable site in the ORF 3 structural glycoprotein. Virology 267:135–140
Meulenberg JJ, van Nieuwstadt AP, van Essen-Zandbergen A, Langeveld JP (1997) Posttranslational processing and identification of a neutralization domain of the GP4 protein encoded by ORF4 of Lelystad virus. J Virol 71:6061–6067
Vanhee M, Costers S, Van Breedam W, Geldhof MF, Van Doorsselaere J, Nauwynck HJ (2010) A variable region in GP4 of European-type porcine reproductive and respiratory syndrome virus induces neutralizing antibodies against homologous but not heterologous virus strains. Viral Immunol 23:403–413
Oleksiewicz MB, Botner A, Toft P, Normann P, Storgaard T (2001) Epitope mapping porcine reproductive and respiratory syndrome virus by phage display: the nsp2 fragment of the replicase polyprotein contains a cluster of B-cell epitopes. J Virol 75:3277–3290
Costers S, Vanhee M, Van Breedam W, Van Doorsselaere J, Geldhof M, Nauwynck HJ (2010) GP4-specific neutralizing antibodies might be a driving force in PRRSV evolution. Virus Res 154:104–113
Tian D, Wei Z, Zevenhoven-Dobbe JC, Liu R, Tong G, Snijder EJ, Yuan S (2012) Arterivirus minor envelope proteins are a major determinant of viral tropism in cell culture. J Virol 86:3701–3712
Li B, Fang L, Liu S, Zhao F, Jiang Y, He K, Chen H, Xiao S (2010) The genomic diversity of Chinese porcine reproductive and respiratory syndrome virus isolates from 1996 to 2009. Vet Microbiol 146:226–237
Verheije MH, Olsthoorn RC, Kroese MV, Rottier PJ, Meulenberg JJ (2002) Kissing interaction between 3′ noncoding and coding sequences is essential for porcine arterivirus RNA replication. J Virol 76:1521–1526
Frydas IS, Trus I, Kvisgaard LK, Bonckaert C, Reddy VR, Li Y, Larsen LE, Nauwynck HJ (2015) Different clinical, virological, serological and tissue tropism outcomes of two new and one old Belgian type 1 subtype 1 porcine reproductive and respiratory virus (PRRSV) isolates. Vet Res 46:37
Han K, Seo HW, Park C, Oh Y, Kang I, Han HJ, Kim SH, Chae C (2013) Comparative pathogenicity of three Korean and one Lelystad type 1 porcine reproductive and respiratory syndrome virus (pan-European subtype 1) isolates in experimentally infected pigs. J Comp Pathol 149:331–340
The authors declare that they have no competing interests.
LJS, HK, BL, CR, IH-P and TR planned and conceived the experiments. LJS, LZ, HK and IH-P performed the animal trial. GM and BJ provided field samples and executed qRT-PCR on nasal swabs. LJS, LZ and BL performed all other laboratory work and analysed the data. LJS, BL and TR interpreted the results and designed the figures. LJS and TR wrote the manuscript. All authors read and approved the manuscript.
The authors thank Herbert Weißenböck for careful examination of histological lung samples, Graham Tebb for detailed correction of the manuscript text, Alexander Tichy for valuable help with the statistical analysis and Ludwig Haas for providing the isolate GER09-613.
Additional file 1. Phylogenetic analysis of ORF5. Phylogenetic tree based on ORF5 nucleotide sequences of 41 PRRSV-1 strains and PRRSV-2 prototype VR2332 as an out-group as well as the two closest related sequences in NCBI BLASTn phylogenetic tree for each isolate. The PRRSV strains presented in this study are marked in red. The tree was constructed using the neighbour joining method with the numbers at the nodes representing bootstrap values in % of 1000 replicates. Scale bar: number of substitutions per site.
Additional file 2. Phylogenetic analysis of ORF7. Phylogenetic tree based on ORF7 nucleotide sequences of 41 PRRSV-1 strains and PRRSV-2 prototype VR2332 as an out-group as well as the two closest related sequences in NCBI BLASTn phylogenetic tree for each isolate. The PRRSV strains presented in this study are marked in red. The tree was constructed using the neighbour joining method with the numbers at the nodes representing bootstrap values in % of 1000 replicates. Scale bar: number of substitutions per site
About this article
Cite this article
Sinn, L.J., Zieglowski, L., Koinig, H. et al. Characterization of two Austrian porcine reproductive and respiratory syndrome virus (PRRSV) field isolates reveals relationship to East Asian strains. Vet Res 47, 17 (2016). https://doi.org/10.1186/s13567-015-0293-x