Next Article in Journal
Autonomic Nervous System Adaptation and Circadian Rhythm Disturbances of the Cardiovascular System in a Ground-Based Murine Model of Spaceflight
Next Article in Special Issue
Histopathological Pulmonary Lesions in 1st-Day Newborn Piglets Derived from PRRSV-1 MLV Vaccinated Sows at the Last Stage of Gestation
Previous Article in Journal
Bladder Neck Obstruction: Experience and Management in a Sperm Bank
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 13
Issue 3
10.3390/life13030843
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Characteristics and Pathogenicity of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) 1 in Taiwan during 2019–2020

by
Fu-Chun Hsueh
1,
Kun-Lin Kuo
1,2,
Feng-Yang Hsu
1,2,
Sheng-Yuan Wang
1,
Hsien-Jen Chiu
1,2,
Meng-Tien Wu
1,
Chuen-Fu Lin
1,2,
Yu-Han Huang
1,
Ming-Tang Chiou
1,2,* and
Chao-Nan Lin
1,2,*
1
Animal Disease Diagnostic Center, College of Veterinary Medicine, National Pingtung University of Science and Technology, Pingtung 91201, Taiwan
2
Department of Veterinary Medicine, College of Veterinary Medicine, National Pingtung University of Science and Technology, Pingtung 91201, Taiwan
*
Authors to whom correspondence should be addressed.
Life 2023, 13(3), 843; https://doi.org/10.3390/life13030843
Submission received: 23 February 2023 / Revised: 16 March 2023 / Accepted: 19 March 2023 / Published: 21 March 2023
(This article belongs to the Special Issue Porcine Reproductive and Respiratory Syndrome Virus (PRRSV))
Browse Figures
Versions Notes

Abstract

:
Two variants of porcine reproductive and respiratory syndrome virus (PRRSV), PRRSV 1 and PRRSV 2, have caused abortion in pregnant sows and respiratory distress in nursery pigs worldwide. PRRSV 2 has been thoroughly researched in Taiwan since 1993; however, the first case of PRRSV 1 was not reported until late 2018. To decipher the genetic characteristics of PRRSV 1 in Taiwan, open reading frame 5 (ORF5) genes of PRRSV 1 strains collected from 11 individual pig farms in 2019–2020 were successfully sequenced. All Taiwanese ORF5 sequences were closely related to Spanish-like PRRSV strains, which are considered to share a common evolutionary origin with the strain used for the PRRSV 1 vaccine. Analyses of amino acid (aa) and non-synonymous substitutions showed that genetic variations resulted in numerously specific codon mutations scattered across the neutralizing epitopes within the ORF5 gene. The PRRSV 1 challenge experiment disclosed the pathogenetic capability of the NPUST2789 isolate in nursery pigs. These findings provide comprehensive knowledge of the molecular diversity of the PRRSV 1 variant in local Taiwanese fields and facilitate the development of suitable immunization programs against this disease.

1. Introduction

Since porcine reproductive and respiratory syndrome virus (PRRSV) was first discovered in Europe and the USA [1,2], it has resulted in middle- or late-term reproductive failure in pregnant sows and respiratory disorders in nursery pigs. This porcine disease, called porcine reproductive and respiratory syndrome (PRRS), also devastates the swine industry, causing great economic losses worldwide [3]. As an enveloped, positive-stranded RNA virus, PRRSV belongs to the order Nidovirales, family Arteriviridae, genus Betaarterivirus and can be classified into two species, Betaarterivirus suid 1 (PRRSV 1, formerly called the European type) and Betaarterivirus suid 2 (PRRSV 2, formerly called the North American type), with approximately 44% genetic differences between each other [4,5]. Both PRRSV strains manifest the same clinical signs in naïve sows and gilts [6]. In Taiwan, the first outbreak of PRRSV 2 was reported in 1993 [7]. Since then, PRRSV 2 has uninterruptedly circulated in swine herds in the absence of effective management or strategies against it. Later, PRRSV 1 emerged in Taiwanese fields in late 2018, hindering the birth of piglets and undermining the confidence in vaccination [8].
Porcine reproductive and respiratory syndrome virus is approximately 15.1 kb in length and encodes 10 open reading frames (ORFs), ORF1a, 1b, 2a, 2b, 3, 4, 5, 5a, 6 and 7 [9]. Among them, the highly mutated ORF5 gene, which encodes glycoprotein 5 (GP5), plays a pivotal role in inducing viral neutralizing antibodies and imparting cross-protection to pigs against heterologous PRRSV strains [10,11]. GP5 consists of three regions, the ectodomain, the hydrophobic region and the putative endodomain [12]. Specifically, two neutralizing epitopes have been identified within the ectodomain, the N-terminal ectodomain (amino acid (aa) 29–35) [12] and the ectodomain epitope (aa 38–54) [13], which are recognized by neutralizing monoclonal antibodies. Additionally, an immunodominant T-epitope, also named the major histocompatibility complex class I (MHC-I) epitope, has been located in the region of aa 117–128. This epitope has the capability of activating cell-mediated immunity against PRRSV by stimulating IFN-γ secreting cells [14]. On the other hand, a deceptive “decoy epitope” situated in the hypervariable area of GP5 plays an important role in evading neutralizing responses by the N-glycan-shielding mechanism [15,16]. The variation in N-glycosylation sites may interfere with the recognition by neutralizing antibodies and further compromise the immunogenicity of neutralizing epitopes since specific glycosylation sites could differ during viral replication in vivo [17,18].
As a highly mutated virus, PRRSV 1 was initially separated into four different genotypes based on the ORF5 sequence [19]. A few years later, the novel phylogenetic classification system further divided PRRSV 1 into 12 sub-clades to explicitly accentuate the genomic diversity [20]. Following the incorporation of nucleotide sequences from central and eastern European countries, a more solid classification system proposed that PRRSV 1 should be categorized into three lineages [21]. Among them, Lineage 1 is composed of the prototype PRRSV 1, Lelystad strain as well as PRRSV DV (Merck, Sharp and Dohme Animal Health, Rahway, NJ, USA), Pyrsvac-183 (Laboratorios Syva, S.A., León, Spain), Amervac PRRSV (Laboratorios Hipra, S.A., Girona, Spain) and Unistrain PRRSV (Laboratorios Hipra, S.A., Girona, Spain) vaccine strains [21]. Currently, live-attenuated or modified live vaccines have similar concerns over safety, meager heterologous cross-protection and the reversion of virulence [22,23]. This is the common warning of all commercial PRRSV modified live vaccines, which emphasizes the importance of adjustable immune programs. Thus, the elucidation of the molecular diversity of PRRSV 1 has become an essential task to illuminate the actual epidemic in the field and further facilitate more appropriate vaccination strategies for disease control. Here, this study records the molecular characteristics of the ORF5 gene and pathogenicity of PRRSV 1 in Taiwan following the first outbreak in 2018. These results provide the in situ status of currently circulating PRRSV 1 and will promote the advancement of effective procedures to deal with this virus.

2. Materials and Methods

2.1. Sample Collection

In 2019–2020, clinical specimens (serum or tissues) were collected from 15,034 pigs for PRRSV detection from northern, central, southern and eastern Taiwan. Severe reproductive and/or respiratory epidemics took place on 11 different pig farms in central and southern Taiwan. On those farms, pregnant sows manifested abortion, prematurity and stillbirth, and nursery pigs had high fever, severe dyspnea and inappetence accompanied by edematous eyelids and discolored ears. Detailed information on vaccination and background history is summarized in Table 1. Blood samples and lung tissues were submitted to the Animal Disease Diagnostic Center at National Pingtung University of Science and Technology for molecular diagnosis. All specimens were confirmed as PRRSV 1-positive by quantitative reverse transcription polymerase chain reaction (RT-qPCR) performed as previously mentioned [8].

2.2. Amplification of the ORF5 Gene and Phylogenetic Analysis

Viral nucleic acids were first extracted using the MagNA Pure LC total nucleic acid isolation kit (Roche Diagnostics, Mannheim, Germany) and then used to synthesize complementary DNA (cDNA) using the PrimeScriptTM RT reagent kit (Takara Bio Inc., Kusatsu, Shiga, Japan). Two microliters of template cDNA was amplified utilizing KAPA HiFi HotStart ReadyMix (Roche), and two primer pairs were designed outside of the ORF5 gene (Table 2) under the following thermal conditions: initial denaturation at 95 °C for 3 min; 40 cycles of 98 °C for 20 s, 56 °C for 30 s and 72 °C for 50 s; and a final extension at 72 °C for 50 s. PCR products were subjected to 2% agarose gel electrophoresis, visualized by ultraviolet illumination after ethidium bromide staining, and the TA cloning of ORF5 gene as described in a previous study [8] was performed and confirmed by nucleotide sequencing (Mission Biotech Inc., Taipei, Taiwan). Multiple sequence alignment of historical global PRRSV 1 reference strains was established using Clustal W in the Molecular Evolutionary Genetics Analysis (MEGA), version 11 software program to construct the phylogenetic trees by the maximum likelihood estimation (MLE) method in the Kimura 2-parameter (K2P) mode using a discrete gamma distribution and bootstrap methods estimated for 1000 replications.

2.3. Identification of Selected Sites of the ORF5 Gene

To precisely deduce the positive/diversifying and negative/purifying selection sites, four algorithms, Fixed Effects Likelihood (FEL), Single-Likelihood Ancestor Counting (SLAC), Fast Unconstrained Bayesian AppRoximation (FUBAR) and Mixed Effects Model of Evolution (MEME), were performed. The ratio of non-synonymous substitutions (beta) and synonymous substitutions (alpha) per site based on the maximum likelihood tree was calculated using the FEL, SLAC and MEME methods at a significance level of 0.05 [24], while the posterior probability of positive or negative selection at a single site was evaluated using the FUBAR method with a threshold set at 0.95 [25]. All of the abovementioned methods were implemented on the Datamonkey web server of the HyPhy software package (http://www.datamonkey.org (accessed on 14 February 2022)) [26].

2.4. Animal Experiment

To evaluate the pathogenicity of the PRRSV 1 strain in pigs, we performed a previous isolate, NPUST2789 (GeneBank accession number: MN242825) [8] as the PRRSV 1 challenged variant in this experiment. A total of 21 4-week-old Landrace × Yorkshire × Duroc crossed-breed pigs were automatically separated into five groups: the high-inoculated (HIN) group (n = 3), the high-inoculated contacted (HC) group (n = 6), the low-inoculated (LIN) group (n = 3), the low-inoculated contacted (LC) group (n = 6) and the vehicle group (n = 3). No pigs received the PRRSV vaccines in this study. All pigs were first confirmed to be PRRSV 1-seronegative negative and then intranasally challenged with 2 mL of 105 and 103 TCID50/mL NPUST2789 in the HIN and LIN group, respectively. As for the HC and LC groups, those pigs were herded in the same batch as pigs from the HIN and LIN groups without being challenged. Body weights and body temperatures of pigs were monitored weekly and daily, respectively. Serum samples from pigs were collected in Vacutainer® Plus Plastic SSTTM tubes with polymer gel (BD Medical, East Rutherford, NJ, USA) at 0-, 3-, 7-, 10- and 14-days post-inoculation (DPI) for evaluating both PRRSV 1 antibodies and viral loads by using the IDEXX PRRS X3 Ab Test kit (IDEXX, Westbrook, ME, USA) and RT-qPCR, respectively. After euthanasia for all PRRSV 1-challenged pigs, all lobes of the lung were excised, fixed in 10% neutral-buffered formalin, paraffin-embedded and stained with hematoxylin for examining the PRRSV 1 invasion in lung tissue sections.

2.5. Statistics

Body weights, body temperatures and viral loads were depicted and statistically measured using GraphPad Prism 9.5.0 (GraphPad Software, San Diego, CA, USA) with two-way analysis of variance (ANOVA). Statistically significant differences were represented by p-values < 0.001.

3. Results

3.1. PRRSV 1 Nucleotide Sequencing and Phylogenetic Analysis of the ORF5 Gene

A total of 5464 (36.3%) samples from 15,034 pigs were positive for PRRSV. Among them, 158 (2.9%) pigs were PRRSV 1-positive, 5275 (96.5%) were PRRSV 2-positive, and 31 (0.6%) were positive for both PRRSV 1 and 2 (Table 3). Severe reproductive and/or respiratory epidemics that took place on 11 different pig farms were further analyzed in the present study. A total of 11 full-length ORF5 genes were successfully amplified and sequenced. The location of the herds encompassed five separate counties, Pingtung County (four herds), Kaohsiung City (two herds), Miaoli County (one herd), Yunlin County (three herds) and Hsinchu County (one herd), with details on case history and vaccination information listed in Table 1. The nucleotide lengths of the ORF5 gene of the PRRSV 1 strains in this study were all 606 bp. The ORF5 nucleotide sequences of these 11 strains (GenBank accession no. OM677752-OM6777762) showed 93.7–99.7% identity to each other. Compared with the rest of the PRRSV 1 strains in the blue cluster and the remaining unmarked viral variants (Figure 1), our Taiwan PRRSV 1 strains manifested 87.9–99.8% and 79.9–88.2% nucleotide identities, respectively. Phylogenetic analysis of all Taiwanese PRRSV 1s (Figure 1) revealed that they belonged to Lineage 1 based on the previous classification system [21], which was established from prototype PRRSV isolates in Western Europe in 1991. These local variants were closely related to the PRRSV 1 strains Vp-046_bis and All-183 (GenBank accession nos. DQ345725 and DQ345726), with nucleotide identities of 99.84% and 99.67%, respectively.

3.2. Amino Acid Analysis of ORF5 of PRRSV 1

Multiple amino acid alignment of ORF5 showed that all 11 Taiwanese PRRSV isolates in 2019–2020 displayed prolific posterity aa sites analogous to the Spanish-like PRRSV 1 strains [27]. However, the PRRSV 1 variants in Taiwan still exhibited numerous aa mutations at a percentage of 4.6–20.3% in comparison with all the other historical global viral strains. These aa mutations scattered across neutralizing epitopes of the ORF5 gene comprised A32V in the N-terminal ectodomain, D37N/S in the ectodomain and C117F, F119L and F122L in the T-cell epitope (Figure 2). Furthermore, some unique aa mutations were also identified in the majority of Taiwanese PRRSV strains, including G8E, A32V, G36D/V, D56E/K, G63D, A101T/I, C111S, V154I, D172G and N174D (Figure 2). Of note, 3 aa mutations were discovered at the N-glycosylation site, 2 of which, S60N and G63N, were discovered in the PRRSV 1/3599 12w-2/TW-2019 strain, and 1 of them, D56N, was noticed in the PRRSV 1/4260 8w-1/TW-2020 variant (Figure 2).

3.3. Estimation of Non-Synonymous and Synonymous Substitutions of the ORF5 Gene

Non-synonymous (beta) and synonymous (alpha) estimates of each codon among 11 Taiwanese variants and 6 historical variants, including Vp-046_bis, All-183, DV, Lelystad virus, AF0611 and AS1210 (Figure 1), were analyzed. The average value of beta/alpha was 0.287. The results of FEL analysis indicated that 1 specific site at codon 90 was identified under diversifying positive selection (Figure 3) with a likelihood ratio test (LRT) at 5.063 and p value = 0.024. Using the FUBAR method, a similar result was obtained at codon 90 by positive selection, but 5 codon sites were also highlighted by negative selection (Table 4). The MEME method indicated that a distinct site at codon 36 was computed by positive selection with a non-synonymous substitution rate of 220.51, an LRT of 5.78, and a p value of 0.03. No positive or negative selection was found by the SLAC method.

3.4. Pathogenicity of Virulent NPUST2789 Challenge in Nursery Pigs

The variations of both body weights and body temperatures in each group are recorded in Figure 4A,B. No statistically significant difference in weekly body weights and daily body temperatures was observed among different groups. During this study, no sign of fever was discovered in a single pig. Regarding the PRRSV antibodies, all pigs in the HIN group and 1 pig in the LIN group showed positive after 3 DPI. The levels of PRRSV 1 viral loads in each group were depicted in Figure 4C. In the HIN group, the average viral titers in serum were 1.94 ± 0.45 log10 copies/µL at 3 DPI, and gradually climbing at 7, 10 and 14 DPI with viral titers of 2.51 ± 0.72, 3.80 ± 1.01 and 4.10 ± 0.58 log10 copies/µL, respectively. As to the LIN group, only 1 pig showed positive viral titers starting from 0.24 ± 0.34 log10 copies/µL at 3 DPI, 0.67 ± 0.94 log10 copies/µL at 7DPI, to a viral peak of 0.81 ± 1.14 log10 copies/µL at 10 DPI and finally 0.81 ± 1.14 log10 copies/µL at 14 DPI. In the HC group, the viral titer of 0.80 ± 1.13 log10 copies/µL was discovered at 14 DPI. Significantly statistical differences (p < 0.001) were noticed at 14 DPI. Moreover, the histopathological examination manifested that the PRRSV-associated lesions of severe interstitial pneumonia were only noted in pigs at the HIN group (Figure 4D).

4. Discussion

The earliest case of field-circulating PRRSV 1 in Taiwan was discovered in 2018, which caused numerous abortions and stillbirths in sows [8]. Thereafter, no study focused on the phylogenetic variation of PRRSV 1 in Taiwan. To fill this gap, 11 full-length ORF5 genes of PRRSV 1 were sequenced from 11 separate pig farms across 5 different geographic areas in Taiwan. All Taiwanese strains were closely related to historic Spanish PRRSV 1s, Vp-046_bis and All-183, which had been designated Amervac PRRSV (Laboratorios Hipra) and Pyrsvac-183 (Laboratorios Syva), respectively. Multitudinous amino acid mutations scattered across the neutralizing epitopes of Taiwanese PRRSV 1 glycoprotein 5 were identified and compared with each other or with historical global viral strains. In terms of sample sizes, this limitation of only 11 sequences is due the low detection rate of PRRSV 1 (2.9%) during 2019–2020 in Taiwan. To our knowledge, this is the first study to demonstrate the phylogenetic characteristics of PRRSV 1 in Taiwan. These findings can serve as essential information for developing a useful strategy for disease control.
Porcine reproductive and respiratory syndrome virus 1 has been recognized as a pathogen with high evolutionary rates, especially in the ORF5 gene [28], which resulted in recombinant or mutation events [29]. Among the ectodomain and the hydrophobic region, Taiwanese PRRSV 1 variants exhibited numerous aa mutations scattered across the N-terminal domain (1 site), the ectodomain epitope (1 site) and the rest of the region (14 sites). Of note, unique aa mutations at specific N-glycosylation sites, such as D37N, D56N, S60N and G63N, are closely associated with immune evasion, inhibition of viral recognition by the host and loss of neutralizing capabilities [27,30]. Another aa mutation (R9C) at a conserved cysteine residue, which might alter the folding and function of glycoprotein [31], was noticed in three of our sequences. Furthermore, amid the T-cell epitope, four aa mutations, including C117F at the conserved cysteine site and F119L, F122L and V126A, have also been discovered in vaccine-derived wild-type PRRSV 1 [31,32]. Collectively, these results might be a feasible explanation for the incompletion of PRRSV vaccines; in addition, the highest titer of PRRS 1 RNA was found in the diseased pigs (Table 1). Therefore, challenge trials of local wild-type PRRSV 1 isolate are crucial to demonstrate its actual pathogenicity. In this study, we evaluated the previous PRRSV 1 isolate, NPUST2789, challenged in nursery pigs and verified the PRRSV 1 viral loads in serum and the PRRSV 1-related lesion in the lung tissue sections, which point outs that Taiwanese wild-type PRRSV 1s have potentially high infectious capabilities.
The primary concern of live-attenuated vaccines is the recovery of virulence [23]. In Asia, cases of vaccine-derived wild-type PRRSV 1 infection are not infrequent [29,33]. One Thai study documented that most of their isolates shared over 99% identities in the ORF5 gene with the European vaccine, even though modified live virus (MLV) vaccines for PRRSV 1 have not been used in local fields [33]. Interestingly, there were comparable events on Farms E and K, which did not receive the PRRSV 1 vaccine. This might be due to the international trade of breeding pigs or semen from foreign countries and the close interactions between farms that have administered the live-attenuated PRRSV 1 vaccine. In Korea, two studies demonstrated that the Korean viral strains were derived from a common ancestor of the Spanish type 1 PRRSV [32] and the VP-046-like viral strain [34]. A similar study in China also indicated that their local isolates might have originated from Amervac PRRSV vaccine-related isolates and undergone further genetic evolution [35]. In European countries, two studies in Denmark [31] and the Netherlands [36] mentioned that both their wild-type isolates strongly resembled the PRRSV vaccine strain DV. In the current study, our PRRSV 1 variants were comparable to the Spanish-like PRRSV 1 vaccine strain. Furthermore, potential recombinant PRRSV strains derived from vaccines and wild-type viral strains had also discovered what was mentioned above. PRRSV MLV vaccines do confer late but effective protection against homologous wild-type virus, while only limited protection or no protection against heterologous strains [37]. Therefore, both efficacy and safety should be concerns for the use of MLV vaccines [38]. Taken together, authorized vaccines are still the optimal solution against PRRSV, but they must be administered only under the proper evaluation by veterinarians, and the vaccine schedule must be adjusted for adequate immunization. Moreover vaccination and strict external biosecurity protocols are still the only choice for the control and prevention of PRRSV in Taiwan in the future.

5. Conclusions

In the present study, we provide up-to-date knowledge with regard to the phylogenetic diversity of circulating PRRSV 1 in Taiwan after the first outbreak in late 2018. Miscellaneous amino acid mutations in the glycoprotein 5 might be ascribed to the significant evolutionary capabilities of PRRSV 1. Moreover, the challenge experiment proved that the PRRSV 1 isolate, NPUST2789, has the ability to induce viremia and to damage lungs in nursery pigs, and accordingly, this highlights the urgency of the appropriate immunization program using live-attenuated vaccines. Our results can serve as the basis for developing suitable strategies for PRRS control.

Author Contributions

Conceptualization, M.-T.C. and C.-N.L.; methodology, F.-C.H., K.-L.K., F.-Y.H., S.-Y.W., H.-J.C. and Y.-H.H.; software, F.-C.H., K.-L.K., M.-T.W. and C.-F.L.; validation, F.-C.H., K.-L.K., C.-F.L., M.-T.C. and C.-N.L.; formal analysis, F.-C.H., H.-J.C. and K.-L.K.; investigation, F.-C.H., K.-L.K. and H.-J.C.; resources, M.-T.C. and C.-N.L.; data curation, F.-C.H. and K.-L.K.; writing—original draft preparation, F.-C.H. and K.-L.K.; writing—review and editing, M.-T.C. and C.-N.L.; visualization, M.-T.C. and C.-N.L.; supervision, M.-T.C. and C.-N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministry of Science and Technology, Taiwan (MOST 110-2313-B-020-009-MY2).

Institutional Review Board Statement

All clinical samples were submitted to the ADDC of NPUST for diagnostic purposes. The animal experimental procedure was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of the National Pingtung University of Science and Technology with approval No. NPUST-109-123.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank all members in Animal Disease Diagnostic Center of National Pingtung University of Science and Technology for assisting the diagnostic works. We would like to thank Yung-Cheng Lin for comments that greatly improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wensvoort, G.; Terpstra, C.; Pol, J.M.; ter Laak, E.A.; Bloemraad, M.; de Kluyver, E.P.; Kragten, C.; van Buiten, L.; den Besten, A.; Wagenaar, F. Mystery Swine Disease in The Netherlands: The Isolation of Lelystad Virus. Vet. Q. 1991, 13, 121–130. [Google Scholar] [CrossRef] [PubMed]
  2. Collins, J.E.; Benfield, D.A.; Christianson, W.T.; Harris, L.; Hennings, J.C.; Shaw, D.P.; Goyal, S.M.; Mccullough, S.; Morrison, R.B.; Joo, H.S.; et al. Isolation of Swine Infertility and Respiratory Syndrome Virus (Isolate ATCC VR-2332) in North America and Experimental Reproduction of the Disease in Gnotobiotic Pigs. J. Vet. Diagn. Invest. 1992, 4, 117–126. [Google Scholar] [CrossRef] [PubMed]
  3. Rossow, K.D. Porcine Reproductive and Respiratory Syndrome. Vet. Pathol. 1998, 35, 1–20. [Google Scholar] [CrossRef] [PubMed]
  4. Darwich, L.; Gimeno, M.; Sibila, M.; Diaz, I.; de la Torre, E.; Dotti, S.; Kuzemtseva, L.; Martin, M.; Pujols, J.; Mateu, E. Genetic and Immunobiological Diversities of Porcine Reproductive and Respiratory Syndrome Genotype I Strains. Vet. Microbiol. 2011, 150, 49–62. [Google Scholar] [CrossRef]
  5. Adams, M.J.; Lefkowitz, E.J.; King, A.M.Q.; Harrach, B.; Harrison, R.L.; Knowles, N.J.; Kropinski, A.M.; Krupovic, M.; Kuhn, J.H.; Mushegian, A.R.; et al. Changes to Taxonomy and the International Code of Virus Classification and Nomenclature Ratified by the International Committee on Taxonomy of Viruses (2017). Arch. Virol. 2017, 162, 2505–2538. [Google Scholar] [CrossRef] [Green Version]
  6. Yang, S.; Kang, I.; Jeong, J.; Oh, T.; Park, K.H.; Park, S.-J.; Ham, H.J.; Jin, G.-R.; Lee, B.-H.; Chae, C. A Comparison of Two Commercially Available Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Modified-Live Virus Vaccines Analyzing the Growth Performance in 1-Day-Old Vaccinated Swine Located on Endemic Farms Co-Circulating PRRSV-1 and PRRSV-2. J. Vet. Med. Sci. 2020, 82, 224–228. [Google Scholar] [CrossRef] [Green Version]
  7. Chang, C.C.; Chung, W.B.; Lin, W.M.; Yang, P.C.; Weng, C.N.; Chiu, Y.T.; Chang, W.F.; Chu, R.M. Porcine Reproductive and Respiratory Syndrome (PRRS) in Taiwan. I: Viral Isolation. J. Chin. Soc. Vet. Sci. 1993, 19, 268–276. [Google Scholar]
  8. Lin, W.H.; Kaewprom, K.; Wang, S.Y.; Lin, C.F.; Yang, C.Y.; Chiou, M.T.; Lin, C.N. Outbreak of Porcine Reproductive and Respiratory Syndrome Virus 1 in Taiwan. Viruses 2020, 12, 316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Kuhn, J.H.; Lauck, M.; Bailey, A.L.; Shchetinin, A.M.; Vishnevskaya, T.V.; Bào, Y.; Ng, T.F.F.; LeBreton, M.; Schneider, B.S.; Gillis, A.; et al. Reorganization and Expansion of the Nidoviral Family Arteriviridae. Arch. Virol. 2016, 161, 755–768. [Google Scholar] [CrossRef] [Green Version]
  10. Mardassi, H.; Mounir, S.; Dea, S. Molecular Analysis of the ORFs 3 to 7 of Porcine Reproductive and Respiratory Syndrome Virus, Quebec Reference Strain*. Arch. Virol. 1995, 140, 1405–1418. [Google Scholar] [CrossRef]
  11. Kim, W.-I.; Kim, J.-J.; Cha, S.-H.; Wu, W.-H.; Cooper, V.; Evans, R.; Choi, E.-J.; Yoon, K.-J. Significance of Genetic Variation of PRRSV ORF5 in Virus Neutralization and Molecular Determinants Corresponding to Cross Neutralization among PRRS Viruses. Vet. Microbiol. 2013, 162, 10–22. [Google Scholar] [CrossRef] [PubMed]
  12. Wissink, E.H.J.; van Wijk, H.A.R.; Kroese, M.V.; Weiland, E.; Meulenberg, J.J.M.; Rottier, P.J.M.; van Rijn, P.A. The Major Envelope Protein, GP5, of a European Porcine Reproductive and Respiratory Syndrome Virus Contains a Neutralization Epitope in Its N-Terminal Ectodomain. J. Gen. Virol. 2003, 84, 1535–1543. [Google Scholar] [CrossRef]
  13. Plagemann, P.G.W. GP5 Ectodomain Epitope of Porcine Reproductive and Respiratory Syndrome Virus, Strain Lelystad Virus. Virus Res. 2004, 102, 225–230. [Google Scholar] [CrossRef] [PubMed]
  14. Díaz, I.; Pujols, J.; Ganges, L.; Gimeno, M.; Darwich, L.; Domingo, M.; Mateu, E. In Silico Prediction and Ex Vivo Evaluation of Potential T-Cell Epitopes in Glycoproteins 4 and 5 and Nucleocapsid Protein of Genotype-I (European) of Porcine Reproductive and Respiratory Syndrome Virus. Vaccine 2009, 27, 5603–5611. [Google Scholar] [CrossRef] [PubMed]
  15. Vu, H.L.X.; Kwon, B.; Yoon, K.-J.; Laegreid, W.W.; Pattnaik, A.K.; Osorio, F.A. Immune Evasion of Porcine Reproductive and Respiratory Syndrome Virus through Glycan Shielding Involves Both Glycoprotein 5 as Well as Glycoprotein 3. J. Virol. 2011, 85, 5555–5564. [Google Scholar] [CrossRef] [Green Version]
  16. Ansari, I.H.; Kwon, B.; Osorio, F.A.; Pattnaik, A.K. Influence of N-Linked Glycosylation of Porcine Reproductive and Respiratory Syndrome Virus GP5 on Virus Infectivity, Antigenicity, and Ability To Induce Neutralizing Antibodies. J. Virol. 2006, 80, 3994–4004. [Google Scholar] [CrossRef] [Green Version]
  17. Wissink, E.H.J.; Kroese, M.V.; Maneschijn-Bonsing, J.G.; Meulenberg, J.J.M.; van Rijn, P.A.; Rijsewijk, F.A.M.; Rottier, P.J.M. Significance of the Oligosaccharides of the Porcine Reproductive and Respiratory Syndrome Virus Glycoproteins GP2a and GP5 for Infectious Virus Production. J. Gen. Virol. 2004, 85, 3715–3723. [Google Scholar] [CrossRef]
  18. Wei, Z.; Lin, T.; Sun, L.; Li, Y.; Wang, X.; Gao, F.; Liu, R.; Chen, C.; Tong, G.; Yuan, S. N-Linked Glycosylation of GP5 of Porcine Reproductive and Respiratory Syndrome Virus Is Critically Important for Virus Replication In Vivo. J. Virol. 2012, 86, 9941–9951. [Google Scholar] [CrossRef] [Green Version]
  19. Stadejek, T.; Oleksiewicz, M.B.; Potapchuk, D.; Podgórska, K. Porcine Reproductive and Respiratory Syndrome Virus Strains of Exceptional Diversity in Eastern Europe Support the Definition of New Genetic Subtypes. J. Gen. Virol. 2006, 87, 1835–1841. [Google Scholar] [CrossRef]
  20. Shi, M.; Lam, T.T.Y.; Hon, C.C.; Hui, R.K.H.; Faaberg, K.S.; Wennblom, T.; Murtaugh, M.P.; Stadejek, T.; Leung, F.C.C. Molecular Epidemiology of PRRSV: A Phylogenetic Perspective. Virus Res. 2010, 154, 7–17. [Google Scholar] [CrossRef]
  21. Balka, G.; Podgórska, K.; Brar, M.S.; Bálint, Á.; Cadar, D.; Celer, V.; Dénes, L.; Dirbakova, Z.; Jedryczko, A.; Márton, L.; et al. Genetic Diversity of PRRSV 1 in Central Eastern Europe in 1994-2014: Origin and Evolution of the Virus in the Region. Sci. Rep. 2018, 8, 7811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kimman, T.G.; Cornelissen, L.A.; Moormann, R.J.; Rebel, J.M.J.; Stockhofe-Zurwieden, N. Challenges for Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Vaccinology. Vaccine 2009, 27, 3704–3718. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, L.; Ge, X.; Yang, H. Porcine Reproductive and Respiratory Syndrome Modified Live Virus Vaccine: A “Leaky” Vaccine with Debatable Efficacy and Safety. Vaccines 2021, 9, 362. [Google Scholar] [CrossRef] [PubMed]
  24. Kosakovsky Pond, S.L.; Frost, S.D.W. Not So Different After All: A Comparison of Methods for Detecting Amino Acid Sites Under Selection. Mol. Biol. Evol. 2005, 22, 1208–1222. [Google Scholar] [CrossRef] [Green Version]
  25. Murrell, B.; Moola, S.; Mabona, A.; Weighill, T.; Sheward, D.; Kosakovsky Pond, S.L.; Scheffler, K. FUBAR: A Fast, Unconstrained Bayesian AppRoximation for Inferring Selection. Mol. Biol. Evol. 2013, 30, 1196–1205. [Google Scholar] [CrossRef] [Green Version]
  26. Weaver, S.; Shank, S.D.; Spielman, S.J.; Li, M.; Muse, S.V.; Kosakovsky Pond, S.L. Datamonkey 2.0: A Modern Web Application for Characterizing Selective and Other Evolutionary Processes. Mol. Biol. Evol. 2018, 35, 773–777. [Google Scholar] [CrossRef] [Green Version]
  27. Prieto, C.; Vázquez, A.; Núñez, J.I.; Álvarez, E.; Simarro, I.; Castro, J.M. Influence of Time on the Genetic Heterogeneity of Spanish Porcine Reproductive and Respiratory Syndrome Virus Isolates. Vet. J. 2009, 180, 363–370. [Google Scholar] [CrossRef]
  28. Nguyen, V.G.; Kim, H.K.; Moon, H.J.; Park, S.J.; Chung, H.C.; Choi, M.K.; Park, B.K. A Bayesian Phylogeographical Analysis of Type 1 Porcine Reproductive and Respiratory Syndrome Virus (PRRSV). Transbound Emerg. Dis. 2014, 61, 537–545. [Google Scholar] [CrossRef]
  29. Sun, Q.; Xu, H.; Li, C.; Gong, B.; Li, Z.; Tian, Z.-J.; Zhang, H. Emergence of a Novel PRRSV-1 Strain in Mainland China: A Recombinant Strain Derived from the Two Commercial Modified Live Viruses Amervac and DV. Front. Vet. Sci. 2022, 9, 974743. [Google Scholar] [CrossRef]
  30. Smith, N.; Power, U.F.; McKillen, J. Phylogenetic Analysis of Porcine Reproductive and Respiratory Syndrome Virus Isolates from Northern Ireland. Arch. Virol. 2018, 163, 2799–2804. [Google Scholar] [CrossRef] [Green Version]
  31. Kvisgaard, L.K.; Hjulsager, C.K.; Kristensen, C.S.; Lauritsen, K.T.; Larsen, L.E. Genetic and Antigenic Characterization of Complete Genomes of Type 1 Porcine Reproductive and Respiratory Syndrome Viruses (PRRSV) Isolated in Denmark over a Period of 10 Years. Virus Res. 2013, 178, 197–205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lee, C.; Kim, H.; Kang, B.; Yeom, M.; Han, S.; Moon, H.; Park, S.; Kim, H.; Song, D.; Park, B. Prevalence and Phylogenetic Analysis of the Isolated Type I Porcine Reproductive and Respiratory Syndrome Virus from 2007 to 2008 in Korea. Virus Genes 2010, 40, 225–230. [Google Scholar] [CrossRef] [PubMed]
  33. Thanawongnuwech, R.; Amonsin, A.; Tatsanakit, A.; Damrongwatanapokin, S. Genetics and Geographical Variation of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) in Thailand. Vet. Microbiol. 2004, 101, 9–21. [Google Scholar] [CrossRef] [PubMed]
  34. Kang, H.; Yu, J.E.; Shin, J.-E.; Kang, A.; Kim, W.-I.; Lee, C.; Lee, J.; Cho, I.-S.; Choe, S.-E.; Cha, S.-H. Geographic Distribution and Molecular Analysis of Porcine Reproductive and Respiratory Syndrome Viruses Circulating in Swine Farms in the Republic of Korea between 2013 and 2016. BMC Vet. Res. 2018, 14, 160. [Google Scholar] [CrossRef] [PubMed]
  35. Zhai, S.-L.; Lin, T.; Zhou, X.; Pei, Z.-F.; Wei, Z.-Z.; Zhang, H.; Wen, X.-H.; Chen, Q.-L.; Lv, D.-H.; Wei, W.-K. Phylogeographic Analysis of Porcine Reproductive and Respiratory Syndrome Virus 1 in Guangdong Province, Southern China. Arch. Virol. 2018, 163, 2443–2449. [Google Scholar] [CrossRef]
  36. Dortmans, J.C.F.M.; Buter, G.J.; Dijkman, R.; Houben, M.; Duinhof, T.F. Molecular Characterization of Type 1 Porcine Reproductive and Respiratory Syndrome Viruses (PRRSV) Isolated in the Netherlands from 2014 to 2016. PLoS ONE 2019, 14, e0218481. [Google Scholar] [CrossRef]
  37. Martín-Valls, G.E.; Cortey, M.; Allepuz, A.; Illas, F.; Tello, M.; Mateu, E. Introduction of a PRRSV-1 Strain of Increased Virulence in a Pig Production Structure in Spain: Virus Evolution and Impact on Production. Porc. Health Manag. 2023, 9, 1. [Google Scholar] [CrossRef]
  38. Nan, Y.; Wu, C.; Gu, G.; Sun, W.; Zhang, Y.-J.; Zhou, E.-M. Improved Vaccine against PRRSV: Current Progress and Future Perspective. Front. Microbiol. 2017, 8, 1635. [Google Scholar] [CrossRef] [Green Version]
Life 13 00843 g001 550
Figure 1. A phylogenetic tree based on the ORF5 gene of PRRSV 1 sequences was constructed by maximum likelihood estimation (MLE) using the Kimura 2-parameter (K2P) mode under the gamma distribution with 1000 bootstrapping replicates. Horizontal branch lengths indicate genetic distances among different viral strains. The red triangle represents the Taiwan PRRSV 1 strains isolated in 2019–2020 within the blue cluster. The scale bar indicates nucleotide substitutions per site.
Figure 1. A phylogenetic tree based on the ORF5 gene of PRRSV 1 sequences was constructed by maximum likelihood estimation (MLE) using the Kimura 2-parameter (K2P) mode under the gamma distribution with 1000 bootstrapping replicates. Horizontal branch lengths indicate genetic distances among different viral strains. The red triangle represents the Taiwan PRRSV 1 strains isolated in 2019–2020 within the blue cluster. The scale bar indicates nucleotide substitutions per site.
Life 13 00843 g001
Life 13 00843 g002 550
Figure 2. Alignment of amino acid sequences of glycoprotein 5 of 11 Taiwanese PRRSV 1 variants along with historical global viral strains. Potential neutralizing epitopes are marked above. Ellipses represent the amino acids identical to the reference sequence.
Figure 2. Alignment of amino acid sequences of glycoprotein 5 of 11 Taiwanese PRRSV 1 variants along with historical global viral strains. Potential neutralizing epitopes are marked above. Ellipses represent the amino acids identical to the reference sequence.
Life 13 00843 g002
Life 13 00843 g003 550
Figure 3. Maximum likelihood estimates of non-synonymous rates (beta; β) and synonymous rates (alpha; α) at each site are shown as bars. The line shows the estimates under the null model (α = β). Estimates above 10 at codon 90 are censored at the site under diversifying positive selection at p < 0.05.
Figure 3. Maximum likelihood estimates of non-synonymous rates (beta; β) and synonymous rates (alpha; α) at each site are shown as bars. The line shows the estimates under the null model (α = β). Estimates above 10 at codon 90 are censored at the site under diversifying positive selection at p < 0.05.
Life 13 00843 g003
Life 13 00843 g004 550
Figure 4. Body weights (A), body temperature (B), PRRSV loads in serum (C) and PRRSV-challenged lesions suggestive of diffuse and severe interstitial pneumonia (D) in pigs following the NPUST2789 challenge. The PRRS viral load was detected by a probe−based quantitative reverse transcription PCR (RT−qPCR) and transformed into the log10 copies/µL. The error bars indicate the standard deviation. Statistically significant differences are marked as **** (p < 0.001).
Figure 4. Body weights (A), body temperature (B), PRRSV loads in serum (C) and PRRSV-challenged lesions suggestive of diffuse and severe interstitial pneumonia (D) in pigs following the NPUST2789 challenge. The PRRS viral load was detected by a probe−based quantitative reverse transcription PCR (RT−qPCR) and transformed into the log10 copies/µL. The error bars indicate the standard deviation. Statistically significant differences are marked as **** (p < 0.001).
Life 13 00843 g004
Table
Table 1. Complete information of PRRSV 1 in each pig farm.
Table 1. Complete information of PRRSV 1 in each pig farm.
FarmCollection TimeLocationScale of HerdType of Vaccine 1Age in WeeksAnimal Condition 2Quantitation of PRRSV 1/PRRSV 2 Loads 3Sample
Farm AJanuary 2019Pingtung1400 sowsLive attenuated PRRSV 13Sick1.10 × 106/− 4NPUST3064
Farm BMarch 2019Pingtung170 sowsLive attenuated PRRSV 18Sick5.85 × 106/−NPUST108-355
Farm CApril 2019Kaohsiung200 sowsLive attenuated PRRSV 16–7Sick1.87 × 106/−NPUST108-603
Farm DAugust 2019Kaohsiung140 sowsLive attenuated PRRSV 1 and Modified live PRRSV 24Normal1.90 × 104/−NPUST3554
Farm EAugust 2019Pingtung480 sowsModified live PRRSV 212Normal4.40 × 104/−NPUST3599
Farm FNovember 2019Miaoli400 sowsLive attenuated PRRSV 15Sick1.77 × 106/−NPUST108-2275
Farm GFebruary 2020Yunlin900 sowsLive attenuated PRRSV 15–6Sick2.26 × 105/−NPUST4028
Farm HFebruary 2020Pingtung250 sowsLive attenuated PRRSV 17Normal7.1 × 104/5.2 × 103NPUST4035
Farm IMarch 2020Yunlin200 sowsLive attenuated PRRSV 17Normal1.7 × 104/−NPUST4178
Farm JMarch 2020Hsinchu200 sowsNon provided8Normal4.7 × 104/−NPUST4260
Farm KApril 2020Yunlin1100 hogsModified live PRRSV 25Sick6.64 × 107/−NPUST109-920
1 Vaccination program of the PRRSV vaccine: every 3 months for sows and at 7–10 days of age for piglets; 2 definition of animal condition: “Sick” means that PRRS-related signs were obviously noticed in pigs of the farm; 3 units of PRRSV loads: copies/µL. 4—means PRRSV qPCR negative.
Table
Table 2. The list of primers used for sequencing in this study.
Table 2. The list of primers used for sequencing in this study.
Name of PrimersSequences
EU-ORF5-F15′-TGCATTTCYTGACACCATC-3′
EU-ORF5-R15′-CCCBARRAGTCGGCCRCGWGA-3′
EU-ORF5-F25′-ATCYRCAATGAGGTGGGCTAC-3′
EU-ORF5-R25′-ACYTTNAGGGCRTADATCAT-3′
Table
Table 3. Positive rates of PRRSV 1 and PRRSV2 during 2019–2020.
Table 3. Positive rates of PRRSV 1 and PRRSV2 during 2019–2020.
YearNumberPositive (%)PRRSV 1 Only (%)PRRSV 2 Only (%)PRRSV 1 + PRRSV 2 (%)
201988522864 (32.4)108 (3.8)2747 (95.9)9 (0.3)
202061822600 (42.1)50 (1.9)2528 (97.2)22 (0.9)
Total15,0345464 (36.3)158 (2.9)5275 (96.5)31 (0.6)
Table
Table 4. Results of positive and negative selection by the FUBAR method.
Table 4. Results of positive and negative selection by the FUBAR method.
Codon SiteMean Posterior Beta-AlphaNegative SelectionPositive Selection
Prob [Alpha > Beta] aProb [Alpha < Beta] bBayes Factor [Alpha < Beta] c
77−12.9550.9740.0170.023
909.5780.0180.96131.945
108−17.8700.9840.0090.012
123−14.8160.9790.0140.018
164−20.9860.9880.0070.010
181−13.2730.9720.0180.024
a Posterior probability ≥ 0.95 of negative selection; b posterior probability ≥ 0.95 of positive selection; c empirical Bayes factor for positive selection. Bold font represents the probabilities are above 0.95 or meaningful.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hsueh, F.-C.; Kuo, K.-L.; Hsu, F.-Y.; Wang, S.-Y.; Chiu, H.-J.; Wu, M.-T.; Lin, C.-F.; Huang, Y.-H.; Chiou, M.-T.; Lin, C.-N. Molecular Characteristics and Pathogenicity of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) 1 in Taiwan during 2019–2020. Life 2023, 13, 843. https://doi.org/10.3390/life13030843

AMA Style

Hsueh F-C, Kuo K-L, Hsu F-Y, Wang S-Y, Chiu H-J, Wu M-T, Lin C-F, Huang Y-H, Chiou M-T, Lin C-N. Molecular Characteristics and Pathogenicity of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) 1 in Taiwan during 2019–2020. Life. 2023; 13(3):843. https://doi.org/10.3390/life13030843

Chicago/Turabian Style

Hsueh, Fu-Chun, Kun-Lin Kuo, Feng-Yang Hsu, Sheng-Yuan Wang, Hsien-Jen Chiu, Meng-Tien Wu, Chuen-Fu Lin, Yu-Han Huang, Ming-Tang Chiou, and Chao-Nan Lin. 2023. "Molecular Characteristics and Pathogenicity of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) 1 in Taiwan during 2019–2020" Life 13, no. 3: 843. https://doi.org/10.3390/life13030843

APA Style

Hsueh, F. -C., Kuo, K. -L., Hsu, F. -Y., Wang, S. -Y., Chiu, H. -J., Wu, M. -T., Lin, C. -F., Huang, Y. -H., Chiou, M. -T., & Lin, C. -N. (2023). Molecular Characteristics and Pathogenicity of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) 1 in Taiwan during 2019–2020. Life, 13(3), 843. https://doi.org/10.3390/life13030843

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop