key: cord-295750-0gpyi4ii
authors: Raev, Sergei; Yuzhakov, Anton; Bulgakov, Alexandr; Kostina, Ludmila; Gerasianinov, Alexei; Verkhovsky, Oleg; Zaberezhny, Alexei; Aliper, Taras
title: An Outbreak of a Respiratory Disorder at a Russian Swine Farm Associated with the Co-Circulation of PRRSV1 and PRRSV2
date: 2020-10-15
journal: Viruses
DOI: 10.3390/v12101169
sha: 
doc_id: 295750
cord_uid: 0gpyi4ii

We conducted a cross-sectional study to identify the major respiratory pathogen responsible for an outbreak of respiratory disease at a swine farm in West Siberia in 2019. We discovered that the peak of morbidity and mortality coincided with a high level of porcine reproductive and respiratory syndrome virus (PRRSV) 1 and 2-related viremia. Based on longer PRRSV2 viremia, the dominant role of PRRSV2 over PRRSV1 in the outbreak was assumed. Phylogenetic analysis revealed that the PRRSV1 strain belonged to sub-genotype 2—one of the predominant groups of genotype 1 PRRSVs in Russia. A partial open reading frame 7 sequence of the PRRSV2 isolate demonstrated a high identity with modified live vaccine-related strains from Denmark (93%) and wild-type VR2332 (92%). We identified the first instance of PRRSV1/PRRSV2 mixed infection in Russia. This finding indicates that further field investigations are needed to access PRRSV2 epidemiology in eastern Europe.

Porcine reproductive and respiratory syndrome (PRRS) has been simultaneously described in the USA and Europe and remains a major problem for the pig industry. Economic losses associated with this disease are estimated at US$664 million per year globally. In sows, reproductive failure is characterized by abortions, fetus mummification, stillbirths, and birth of weak offspring. Growing piglets mainly demonstrate respiratory symptoms such as dyspnea, coughing, and fever [1] .

Further investigations have revealed that the causative agent of this disease is porcine reproductive and respiratory syndrome virus (PRRSV), for which two genotypes have been described: the European (PRRSV1) and North American (PRRSV2) genotypes. Significant genetic (50-60% identity) and antigenic differences eventually led to their divergence into two separate species: Betaarterivirus suid 1 (the European type, or PRRSV1) and Betaarterivirus suid 2 (the North American type, or PRRSV2) [2] . Both viruses belong to genus Betaarterivirus in the Variartevirinae subfamily and Arteriviridae family. PRRSV's genome comprises about 15,000 nucleotides and contains ten open reading frames (ORFs) including ORF7, which encodes the nucleocapsid protein [3] .

Shortly after PRRSV's discovery, both PRRSV1 and PRRSV2 were isolated beyond their original places of detection, namely Europe and North America, respectively. PRRSV1 strains are currently present in Europe, North America, and Asia, and PRRSV2 strains are predominant in North America

The blood samples were taken from piglets ages 0, 3, 7, 12, 16, 20 , and 26 weeks (8 to 10 animals per age group). Serum samples were collected on the same day and stored frozen at −70 • C before analysis.

All serum samples were tested for the presence of antibodies against PRRSV capsid protein using two commercial ELISA test kits. The first kit, INgezim PRRS Universal (Ingenasa, Madrid Spain), is prepared on the basis of both PRRSV1 and PRRSV2 recombinant nucleoproteins. This kit detects antibodies against both PRRSV1 and PRRSV2 designated as "pan-PRRSV antibodies". The second assay, RRSS-SEROTEST (Vetbiochem, Moscow, Russia), is capable of detecting antibodies against PRRSV1 only. Both ELISA test kits were used in accordance with the manufacturers' instructions.

PRRSV, PCV2, swine influenza virus (SIV), porcine parvovirus virus (PPV), and porcine respiratory coronavirus (PRCV) in serum samples were detected using commercial polymerase chain reaction (PCR) kits (Vetbiochem, Russia) according to the manufacturer's instructions. In all serum samples, detection was performed using commercial real-time polymerase chain reaction (PCR) kits (Vetbiochem, Russia) in accordance with the manufacturer's instructions. According to previous data, both ORF5 and ORF7 sequences might be chosen for PRRSV1 classification [8] . We used the following primers for the amplification and sequencing of PRRSV1 ORF7, kindly provided by Ivan Trus (Ghe University, Ghent, Belgium): 5 -TGGCCCCTGCCCAICACGT-3 (PRRSV1-ORF7-F) and

Viruses 2020, 12, 1169 3 of 10 5 -TGATCGCCCTAATTGAATAGGTGACT-3 (PRRSV1-ORF7-R). For the amplification and sequencing of PRRSV2 ORF7, we used two pairs of nested primers. For the first round of PCR, we used the following primers: 5 -TTCTGGCCCCTGCCCAYC-3 (F7) and 5 -CGCCCTAATTGAATAGGTGAC-3 (R7). For the second round of PCR as well as for sequencing, we used 5'-CCAAATAACAACGGCAAG-3' (AmerF7_12) and 5'-TCATGCTGAGGGTGATG-3' (AmerR7_12). The expected size of PCR products of PRRSV1 ORF7 and PRRSV2 ORF7 was around 644 and 368 nucleotides, respectively. A one-tube, real-time RT-PCR kit (Alpha Ferment, Moscow, Russia) was used to perform reverse-transcription polymerase chain reaction (RT-PCR). The PCR amplification program consisted of 10 min at 50 • C, 5 min at 95 • C, 35 cycles at 95 • C for 15 s, 55 • C for 15 s, and 72 • C for 30 s, as well as 1 cycle of 5 min at 72 • C. The amplification products were visualized in 1% agarose gel in a buffer containing a mixture of Tris base, acetic acid and EDTA (x1 TAE) after which purification was performed using Monarch DNA Gel Extraction Kit (New England Biolabs, MA, USA) in accordance with the manufacturer's instructions. Sanger sequencing was performed using the Big Dye®Terminator v.3.1 Cycle Sequencing Kit (Applied Biosystems, CA, USA) in accordance with the manufacturer's instructions, and the primers specified above. The nucleotide sequences of the genome fragments were determined using an AB3130 genomic automated analyzer (Applied Biosystems, USA). The nucleotide sequences were analyzed using Lasergene 11.1.0. (DNASTAR, WI, USA). Multiple alignment was performed using MUSCLE (MEGA 7.0.18). Phylogenetic dendrograms were plotted using the maximum likelihood method, GTR model (MEGA 7.0.18). The topology of the trees was confirmed following 1000 bootstrapping replication steps.

Viral RNA was not detected in sera from either newborn or 3-week-old piglets, whereas in 7-week-old piglets, PRRSV1 and PRRSV2, and viral RNA were present in 50% and 40% of serum samples, respectively ( Figure 1 ). The mixed RNA belonging to both viruses was detected in only one 7-week-old piglet. In sharp contrast to PRRSV1, where the RNA was not present in samples from 12-week-old animals, PRRSV2 viremia was detected in 60% of 12-week-old piglets. Simultaneously, PCV2 DNA was detected in 20% of 12-week-old piglets, and viremia rates increased with age until the age of 20 weeks, with 40% of the piglets tested being viremic. All serum samples were negative for SIV, PPV, and PRCV. Whereas the total concentration of PRRSV1-specific antibodies showed a clear tendency to rise from 3 weeks (when all samples were seronegative) to 20 weeks of age, a minor numerical decrease in the total level of antibodies against both strains from 3 to 7 weeks of age was detected ( Figure 2 ).

Comprehensive analysis of total PRRSV antibodies (to both genotypes) and PRRSV1-specific antibodies (Table 1 ) did not reveal any certain tendency. PRRSV2 RNA was found in an only numerically higher (three out of nine) proportion of PRRSV1 IgG negative/total PRRSV Ig-positive samples compared to PRRSV1 (one out of nine). 

Serum samples from 7-week-old piglets were used for PRRSV1 (sample #4) and for PRRSV2 (sample #6) sequencing. The ORF7 sequences were deposited into the GenBank sequence database under the accession numbers ORF7_Kem19EU-MT344126 and ORF7_Kem19NA-MT344127. The ORF7_Kem19EU-MT344126 (PRRSV1) has 378 nucleotides, and further phylogenetic analysis revealed that this isolate belongs to subtype 2 of PRRSV1 ( Figure 3 ). This strain is most closely related to strains Eig иAus from Lithuania and strain Gk from the European part of Russia [8] .

The partial nucleotide sequence of ORF7_Kem19NA-MT344127 (PRRSV2) comprises 337 nucleotides ( Figure 4 ). The comparative analysis of partial sequence ORF7_Kem19NA-MT344127 against those of reference strains demonstrated a high nucleotide identity with wild-type VR2332 (92%) and MLV-related strains from Denmark: KC577602.1_DK-2004-3-1 and AF095479.1_strain_17704A (93%) [14] . The level of identity between the Hungarian non-MLV-related strain KM514315.1:14861-15232_Hungary_102_2012 and ORF7_Kem19NA-MT344127 was only 81%. Nucleotide identity with GU256774.1, the only one PRRSV2 strain detected in Russia, was 88%. 

JQ308798.1_strain_QYYZ KC618159_strain_GXLSF02-2011  ....-....................................-.......................................................................  EF536003_strain_VR2332  ..N.-....................................-.......................................................................  AF184212_strain_SP  ...K-....................................-.............Y.........................................................  EF536000_strain_MN30100  ...K-...N................................-.....................................................................A.  EF442777_strain_MN184  ...K-................R................S..-.S...................................................A.................  KC577603_isolate_DK-2008-16-4-....-....................................-...................................................?...................  KC577602_DK-2004-3-1  ....-....................................-.......................................................................  KF183964_DK-2010-10-13-5  ....-....................................-.......................................................................  AF121131_MD-001_tw91  ..NK-......................S...V....R....-...........................................I...........................  AY881994_strain_FJ-1  ..R.-.....................H........E..I..-...........................................I...........................  AY032626_strain_CH-1a  .R.K-...N................................-.S.....................................................................  AF095479_strain_17704A  ....-....................................-.......................................................................  KM514315_isolate_Hungary_102_2 ...K-....................................-..L........................................I...........................  EF112445_strain_JXA1  ...K-...N..............................R.-...........................................A.................Q.......A.  GU256774_isolate_Irkutsk/RUS/2 ...K-...N..............................R.-...........................................A.................Q.......A.  JQ308798_strain_QYYZ  ...K-..........................A.......N.-......................V.........R......................................  JN654459_strain_NADC30  ..NK-......................S............N-R...............................R......................................  DQ459471_isolate_S1  ....-....................................-.......................................................................  AF396833_isolate_NADC-8_(E) . ...-....................................-.......................................................................  MN073140_isolate_lung23199-S4-....-....................................-.......................................................................  AF396842_isolate_JA-142_(E) . ..K-.R.N.............................I.N-.......................................................................  AY656994_isolate_17198-6  .L..-..............P.........P...........-.............................................................?.........  MK860181_isolate_394-1  ..NK-...G..................S............N-G...............................R............................R.......A.  AF066183_strain_RespPRRS_MLV  ....-....................................-.......................................................................  DQ988080_isolate_Ingelvac_ATP  ...K-.R.N.............................I.N-.......................................................................  AF066384_strain_PrimePac PRRS  ...K-....................................-.............Y.........................................................  KP262345_isolate_BACU_UNAM  ....I...N................................-.....................................................................A.  AY209224_isolate_29D1  ...KK.Q.N................................-.....................................................A.................  EF441808_isolate_05K175 . KC862536_isolate_DK-2011-05-11 ..............................................-...............................................................  DQ324710_strain_Porcilis_PRRS  ..............................................-...............................................................  DQ324712_strain_Pyrsvac-187  ............N.................................-..R............................................................  MT008024_strain_Tyu16  ......R....EN..................K.M......SKR..T-.R.R.......................................A....P..............  JF802085_strain_Lena  ......R.FR..N..................R...T..........-.................V.....G........................P..............  KP889243_strain_SU1-Bel  ......R.VR..N..................K...T..........-...R.............V.............................................  strain_Bor59*  ...R..R..RN.N......D...........K.M............N.R...............V.....A........................P.......... N.N......D...........K.M..........RVN.R...............V.....S.......................................  KC668221_isolate_WestSib13  ...N..R...N.NA.....D...........K.......A....R.N...R.........T...V.....S........................A...I..........  AF438363_strain_Sid  ...R..R...NNN......D........M..R.M.....R......N.R.R.............V.....LA..................A...................  AF438362__strain_Aus  ...R..R..RN.N......D...........K.............VN.................V.....S........................P..............  KC713969_strain_Eig  ...R.PR...N.NI.....D...........K..............S.R...............V.....S.......................................  EU071308_strain_GK . 

The results of this study indicate that the respiratory disease outbreak correlated with increased PRRSV-related viremia rates, while it did not correlate with the prevalence of another major respiratory pathogen, PCV2. Remarkably, PRRSV2-related viremia peaked in a more pronounced manner and lasted longer compared to PRRSV1-related viremia. This might be an indication of the predominant role of PRRSV2 in respiratory outbreak development.

We assume that whereas total PRRSV antibodies from piglets aged 12 weeks might be considered as post-infectious, at 3 weeks of age, these antibodies are of maternal origin. Due to differences in sensitivity and specificity between the test kits used in this study and simultaneous PRRSV1 and PRRSV2 viremia detection, our effort to distinguish between PRRSV1 and PRRSV2 antibodies by using universal (PRRSV1/PRRSV2) and PRRSV1-specific ELISA test kits cannot describe the real situation. After the initial introduction of PRRSV2 into Europe, differentiation between antibodies against PRRSV of different genotypes became crucial to PRRSV diagnosis. An ELISA assay for simultaneous detection and differentiation between antibodies to PRRSV1 and PRRSV2 was developed [15] , but no such ELISA test kit is commercially available in Russia. Hence, the use of universal PCR test kits for simultaneous detection of both viruses is the only option to prevent PRRSV2 from remaining undetected.

Data on the comparative pathogenicity of PRRSV1, PRRSV2, and mixed infection are scarce. In one study, the predominant role of PRRSV2 over PRRSV1 in the course of PRRSV1-1/PRRSV2 (lineage 1) mixed infection was demonstrated on the basis of replication and pathogenicity levels [16] . Significant differences in the pathogenicity of different genotypes/lineages of PRRSV1 and PRRSV2 should be considered in the case of mixed infection. Further experimental studies will be required to establish the impact of at least PRRSV1-2 and PRRSV2 (lineage 5) presence. Despite recombination never being demonstrated between PRRSV1 and PRRSV2 strains, the simultaneous presence of PRRSV1 and PRRSV2 might also be considered a possible risk for the emergence of recombinant strains. The use of MLV vaccines under such circumstances may only increase this risk [17, 18] .

The idea that PRRSV1 is a homogenous group of viruses was revised after the discovery of genetically distinct isolates from eastern Europe, especially from Belarus and Russia [9] . Thus, the detection of a new member of PRRSV1-2 in the present study offers additional proof of the wide distribution of this genotype in Russia. Isolates of PRRSV1-1 (including the so-called Russian group of viruses), PRRSV1-2, and PRRSV1-3 differ significantly in pathogenicity [11, 12, 19, 20] . Marked differences in the nucleotide sequence between PRRSV1-2 isolates and PRRSV1 vaccine strains, which all belong to PRRSV1-1, including differences in the ORF5 sequence (which encodes a major epitope: the envelope glycoprotein), could potentially hamper vaccine efficacy [12] . The benefits of using currently available vaccines against PRRSV1-2, be they live, inactivated, or subunit, remain unknown. The circulation of PRRSV2 further complicates the search for working control strategies, as no vaccines against PRRSV2 are currently available in Russia.

Phylogenetic analysis of a partial ORF7 sequence of the PRRSV2 strain isolated in the course of this study revealed that this strain was very dissimilar (88% identity) to the only known PRRSV2 strain (JXA1-related, sublineage 8.7) detected in East Siberia, Russia, in 2007. This suggests that a different instance of PRRSV2 introduction into Russia occurred at a certain point in addition to the one resulting in the outbreak of 2007. Further genetic analysis (including an ORF5 sequence) should be performed to classify this isolate more precisely.

High homology within lineage 5 of PRRSV2 makes the distinction between the Ingelvac PRRS MLV vaccine strain, its parental strain VR-2332, VR-2332-related wild-type isolates, and vaccine-derived isolates very complicated and requires at least ORF5 sequence analysis [18] . The Boehringer MLV vaccine has not been licensed in Russia. However, there are known cases of the detection of an MLV-like PRRSV1 strain in Vietnam, regardless of the corresponding type of vaccine not being licensed in that country [21] . Notably, one locally-produced MLV vaccine against PRRSV2 has been licensed but is currently commercially unavailable in Russia, and the sequence of the vaccine virus has not yet been published. It is also conceivable that the virus originated from breeder animals infected with an MLV-related PRRSV2 strain. It is well known that MLV-related strains (including strains that have a high identity level with ORF7_Kem19NA-MT344127) are widely spread in Denmark, one of the major exporters of breeding pigs [14] . For example, the average number of pigs imported from Denmark to Russia has been more than 3000 animals a year over the past three years. Finally, an independent introduction of a wild-type VR2332-related strain could result in the emergence of the ORF7_Kem19NA-MT344127 virus.

In the course of infection, NLS interacts with nuclear transport proteins, leading to the penetration of the N protein into the nucleus and nucleolus of infected cells [22] . A critical role of this element of the N protein in pathogenesis has been shown [23] . A similar insertion is known to be found Viruses 2020, 12, 1169 9 of 10 in a Mexican field isolate DMZC144/2016 (GenBank ID MF352184.1) [24] . This might indicate that both Kem19NA and DMZC144/2016 have a nucleoprotein consisting of 124 amino acids in length. Complete ORF7 sequencing should be performed to confirm this assumption. Interestingly, Kem19NA and DMZC144/2016 are genetically distant from each other, their nucleotide identity being only 87%, suggesting that this insertion may not be lineage-specific.

To the best of our knowledge, this is the first report of PRRSV1/PRRSV2 co-circulation in Russia. Considered alongside data from field veterinary laboratories, our observations indicate that PRRSV2 circulates in both the European and Asian parts of the country. The prevalence and distribution of PRRSV2 in Russia therefore warrants further investigation. 

Porcine reproductive and respiratory syndrome: Clinical disease, pathology and immunosuppression

The molecular biology of arteriviruses

Genetic diversity of porcine reproductive and respiratory syndrome virus strains circulating in Hungarian swine herds

Lelystad-like strain of porcine reproductive and respiratory syndrome virus (PRRSV) identified in Canadian swine

Characterization of Emerging European-Like Porcine Reproductive and Respiratory Syndrome Virus Isolates in the United States

Phylogeny-Based Evolutionary, Demographical, and Geographical Dissection of North American Type 2 Porcine Reproductive and Respiratory Syndrome Viruses

Definition of subtypes in the European genotype of porcine reproductive and respiratory syndrome virus: Nucleocapsid characteristics and geographical distribution in Europe

Molecular evolution of PRRSV in Europe: Current state of play

Geographic distribution and molecular analysis of porcine reproductive and respiratory syndrome viruses circulating in swine farms in the Republic of Korea between

Full-genome analysis and pathogenicity of a genetically distinct Russian PRRSV-1 Tyu16 strain

Aliper, T.I. Genetic and pathogenic characterization of a Russian subtype 2 PRRSV-1 isolate

First Outbreak of Atypical Porcine Reproductive and Respiratory Syndrome in Russia Caused by Highly Pathogenic Chinese-Like Prrs Virus

Examination of the selective pressures on a live PRRS vaccine virus

Blocking ELISA's for the distinction between antibodies against European and American strains of porcine reproductive and respiratory syndrome virus

Comparison of the Pathogenesis of Single or Dual Infections with Type 1 and Type 2 Porcine Reproductive and Respiratory Syndrome Virus

Appearance of acute PRRS-like symptoms in sow herds after vaccination with a modified live PRRS vaccine

Comparison of Molecular and Biological Characteristics of a Modified Live Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Vaccine (Ingelvac PRRS MLV), the Parent Strain of the Vaccine (ATCC VR2332), ATCC VR2385, and Two Recent Field Isolates of PRRSV

Pathogenesis and antigenic characterization of a new East European subtype 3 porcine reproductive and respiratory syndrome virus isolate

Pathogenicity of three genetically diverse strains of PRRSV Type 1 in specific pathogen free pigs

Genetic Analysis of ORF5 Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Isolated in Vietnam

The localization of porcine reproductive and respiratory syndrome virus nucleocapsid protein to the nucleolus of infected cells and identification of a potential nucleolar localization signal sequence

Mutations within the nuclear localization signal of the porcine reproductive and respiratory syndrome virus nucleocapsid protein attenuate virus replication

Phylogenetic analysis of ORF5 and ORF7 of porcine reproductive and respiratory syndrome (PRRS) virus and the frequency of wild-type PRRS virus in México

The authors would like to thank Alexander Mishin and Alexander Skrylev for their administrative support and Yana Streltsova for her excellent technical assistance. We thank Tomasz Stadejek and Sara Botti for providing the sequences of the strains Bor-59 and ILI-6.

The authors declare no conflict of interest.