Skip to main content

Genetic analysis of porcine productive and respiratory syndrome virus between 2013 and 2014 in Southern parts of China: identification of several novel strains with amino acid deletions or insertions in nsp2

Abstract

Background

Porcine respiratory and reproductive syndrome virus (PRRSV) is one of the most economically significant pathogens in the Chinese swine industry. ORF5 and nsp2 are highly variable regions of the PRRSV genome. Therefore, nsp2 and GP5 are often selected for investigation of variations and phylogenetic analyses for their genetic diversities. Knowledge of the molecular evolution of PRRSV field strains may contribute to the control of PRRS in China.

Results

The results of multiple sequence alignments of GP5 showed that there is 84.5–100% aa identity among the 56 strains in this study. These strains shared 84.5–99.0% aa identity with the prototypical type 2 PRRSV VR-2332 and 56.6–59.2% with strain LV, prototypical type 1 PRRSV. Phylogenetic analysis showed there is considerable diversity among PRRSV ORF5 and the existence of two lineages (5 and 8). Most of the strains were classified into lineage 8 with multiple sub-lineages (3, 4 and 6). Moreover, PRRSV strains with 5 novel patterns of deletions or insertions in the nsp2 region were found.

Conclusions

Phylogenetic analysis based on ORF5 sequences indicated the diversity of PRRSV in southern parts of China and the strains with 30 aa deletion in nsp2 are dominant in the porcine population. Also, new PRRSV strains with different patterns of deletions or insertions in nsp2 are emerging. The data presented here constitute a useful basis for further epidemiological studies regarding the heterogeneity of PRRSV strains in China and provide a basis for the prevention of PRRS in southern parts of China.

Background

Porcine reproductive and respiratory syndrome virus (PRRSV) is acknowledged as one of the most economically important diseases for the swine industry worldwide [1]. PRRSV, the etiological agent of porcine reproductive and respiratory syndrome (PRRS), is a single-stranded, enveloped, RNA virus. The PRRSV genome consists of approximately15.4 kb and contains a 5′-untranslated region (UTR), open reading frames (ORFs), a 3′-UTR and a 3′-poly(A) tail. The 5′ two-thirds of the genome encodes polyproteins that are processed by viral protease to 14 nonstructural proteins (nsps) [2]. The 3′ one-third region of genome encodes the structural proteins that are translated from a 3′–5′ co-terminal, nested set of subgenomic mRNAs. In addition to the three major structural proteins, GP5, M and N, the genome of PRRSV encodes minor structural proteins, GP2, 2b, GP3, GP4 and ORF5a [3,4,5].

It has been shown that PRRSV is continuously evolving through point mutations and genome recombination, which can lead to some new emerging antigenic variant strains [6]. According to the genetic diversity, PRRSV has been classified as two separate species: type 1 (European) PRRSV and type 2 (North American) PRRSV. The two genotypes share about 60% identity at the nucleotide level [7]. GP5 is highly variable and contains important immunological domains associated with viral neutralization [8]. Nsp2 is the most variable region of PRRSV genome and substitutions, deletions and insertions have been observed in the nsp2 coding region [9, 10]. Therefore, nsp2 and GP5 are often selected for investigation of variations and phylogenetic analyses for their genetic diversities. The genetically extensive variation of PRRSV is likely to pose a major obstacle for the effective control of the most economically significant disease that affects the swine industry [11].

PRRS outbreaks were documented in an intensive pig farm in China at the end of 1995, and it has become one of the most important swine diseases in the Chinese swine industry. In 2006, a large outbreak of porcine high fever syndrome (PHFD), caused by a highly pathogenic PRRSV (HP-PRRSV), emerged in China and affected over 20 million pigs with about 400,000 fatal cases [12]. The disease was subsequently confirmed in southeastern Asian countries including Vietnam, Laos and the Philippines, and this has caused disastrous damage to the swine industry [13, 14]. Whole-genome analysis of the isolated viruses revealed that these PRRSV strains could be grouped into genotype 2 and were highly homologous to a Chinese isolate HB-1 (96.5% nucleotide identity) [15]. It was furthermore observed that these emerging strains contained a single amino acid deletion at position 481 and a 29-amino acid deletion from position 532 to 560 in nsp2 [12, 16]. The study with regards to the origin of the HP-PRRSV did not find recombination or large fragment replacement, which suggests that all HP-PRRSVs originated from the same Chinese ancestor by gradual evolution [15].

Guangxi is one of the biggest breeding regions in the southern parts of China. The genotype prevalence of PRRSV in South China is not currently known. The aim of this study is to investigate the genetic variation of PRRSV using strains isolated in 2013–2014 from pigs that exhibited symptoms of the disease.

Results

Prevalence of PRRSV in Guangxi Province, China from 2013 to 2014

Of the 475 filed samples collected from clinically diseased pigs found between 2013 and 2014 in Guangxi Province of China, 133 samples (28%) were positive for PRRSV, as determined by specific PCR. These results indicated that PRRSV was widely distributed among swine populations in the southern parts of China.

Sequence analyses of the ORF5 gene and nsp2 hypervariable regions of PRRSV

One hundred and thirty-three PRRSV positive samples were used for ORF5 gene and nsp2 HVR amplification. Fifty-six ORF5 sequences and 35 nsp2 HVR sequences were selected for sequencing and analysis (Table 1). To investigate the amino acid difference among PRRSV strains, the GP5 amino acid sequences of 56 PRRSV strains were aligned, together with some North American genotype strains and those from China and other countries. The results showed that all 56 strains encoded a GP5 protein of 200 amino acid residues, but substitutions were extensive (Fig. 1). Sequence alignments revealed that there was 84.5–100% amino acid identity between the 56 Guangxi PRRSV strains and shared 84.5–99.0% amino acid identity with the prototypical type 2 PRRSV VR-2332, 87.0–99.0% with JX-A1, 89.5–92.5% with CH-1a and 56.6–59.2% with prototypical type 1 PRRSV strain LV (Additional file 1: Table S1).

Table 1 Geographic origin and amplified sequence size of target genes from clinical samples in this study
Fig. 1
figure 1

The alignment of GP5 of PRRSV. A multiple alignment of PRRSV GP5 was performed by Clustal W. The sequence of VR-2332 is shown at the top; the residues conserved with it are hidden and substitutions are indicated by the amino acid letter codes. The functional domains are shown within boxes. HRV: highly variable regions, DCE: decoy epitope, PNE: primary neutralizing epitope

To investigate variation in the deduced amino acid sequences of ORF5 gene products, the amino acid sequences of 56 PRRSV strains including some representative strains were aligned. As shown in Fig. 1, critical amino acid variations in some motifs and regions such as the peptide signal, HRV, the decoy epitope (DCE), the primary neutralizing epitope (PNE) and T epitopes were found in GP5 of these strains. Decoy epitopes of VR-2332 is 27VLAN30 and of JXA1 is 27VLVN30. Strains GXLB1403, GXGL1305a/b and GXNN1310b/c contained the same aa in their decoy epitopes (27VLAN30) as those of VR2332. Specific substitution at position 27 (27V → 27A) was found in strain GXWZ1401b, while strains GXGG1305a, GXGG1306 and GXBH1311b contained variations at positions 28 (28L → 28P) and strains GXBS1401a, GXWZ1301a and GXYL1308b have a specific substitution at position 30 (30N → 30D) in the decoy epitope compared to those of VR2332 and JXA1. The other strains have the same aa at the decoy epitope as those of JXA1. Great diversities in HVR1 and HVR2 were found at amino acid positions 32–37 and 57–62, respectively. Specifically, substitutions (S32 → N32, S32 → G32, N34 → S34, N35 → G35 and S36 → G36) in HVR1 of some strains, resulted in loss or gain of the N-linked glycosylation site at specific position in GP5 (Table 2). Amino acid substitutions in the primary neutralizing epitope were also found. Strains GXLB1403, GXGL1305a/b and GXNN1310b/c contained the same aa in the primary neutralizing epitope (37SHLQLIYNL45) as those of VR2332. All PRRSV strains in sub-genotypes 3, representative of JXA1, contained an amino acid mutation (39L → 39I) when compared with VR-2332, while strains GXLB1403 and GXYL1403a contain variations at positions 39 (39L → 39F). Substitutions (N58 → Q/G/K/R58, K59 → H/N59 and D60 → N60) in HVR2 were observed in some strains compared with VR2332. In T cell epitopes, strains GXGG1305a, GXGG1305b and GXGG1306 carried substitutions at positions 126(I126 → T126). Some strains carried substitutions at position at 151(R151 → K151) (Fig. 1).

Table 2 The potential N-glycosylation sites in GP5 of different strains in this study

Identification of several novel strains with amino acids deletions or insertions in nsp2

The nsp2 gene has the highest genetic diversity in the genomes of PRRSV field strains and also was used as an epidemiological genetic marker. To investigate the amino acid differences among PRRSV strains, a predicted 1862-bp DNA fragment containing nsp2 HVR from 35 positive samples was amplified, cloned and sequenced. As shown in Table 1, the amplified nsp2 HVR exhibited various lengths. Compared to strains VR-2332, 1 of 35 nsp2 sequences was 1722 nucleotides in length which is the same as that of VR-2332. 29 out of 35 nsp2 HV region sequences had the same length of 1682 nucleotides, containing the same 30-aa deletion as JXA1 and other HP-PRRSV strains, suggesting that strains with a 30 aa deletion in nsp2 is the dominating strain circulating in the southern parts of China.

Strains GXYL1310, GDHZ1401 and GXBH1404 had the same length of 1322 nucleotides and were found to contain the same 30-aa deletion as JXA1. They also have an extra continuous 120 aa deletion in nsp2. Strains GXBS1401a, GXNN1396 and GXYL1403e contained a discontinuous 31, 49 and 123 aa deletion in their HVR, respectively, compared with strainVR-2332. We also found that one isolated strain (GXNN1407a) contained a 30 aa deletion and 1 aa insertion compared with VR-2332 and JXA1 (Fig. 2). Pairwise comparisons revealed that 86.4–100% nucleotide identity and 84.5–100% amino acid identity between the 35 Guangxi PRRSV strains and shared 65.9%~ 99.2% amino acid identity with the prototypical type 2 PRRSV VR-2332, 87.0–99.0% with JX-A1, 89.5–92.5% with CH-1a and 56.6–59.2% with strain LV of the European type. The data suggested that nsp2 is highly variable and novel HP-PRRSV strains with aa deletions and insertions in the nsp2 are emerging.

Fig. 2
figure 2

Identification of PRRSV strains with amino acid deletions or insertions in nsp2. A multiple alignment of PRRSV NSP2 HV was performed by Clustal W. The sequence of VR-2332 is shown at the top; the residues conserved with it are hidden. The deleted or inserted amino acids are labeled with boxes

Phylogenetic analyses of the ORF5 gene and nsp2 hypervariable region of PRRSV

To gain a better understanding of the genetic relationship, the phylogenetic analysis based on deduced amino acid sequences of ORF5 gene products was conducted by using the 56 ORF5 sequences obtained in this study together with 39 downloaded referenced PRRSV sequences (Table 3). As shown in Fig. 3, the results showed that the PRRSV strains in this study could be divided into 2 different subgroups. Among the 56 GP5 sequences, two strains (GXBL1403 and GXYL1403a) belonged to lineage 5, as represented by VR-2332. Fifty-four strains belonged to lineage 8, with six strains being classified as sub-lineage 8.4 and three strains being classified as sub-lineage 8.6. The other 45 strains formed a large cluster being classified as sub-lineage 8.3 with the representative strains being JXA1, JXwn06 and HUN4.

Table 3 The information of reference strains used in this study
Fig. 3
figure 3

Phylogenetic tree based on a comparison of 95 complete PRRSV ORF5 sequences, including the 56 sequences from this study and 39 ORF5 sequences from PRRSV strains originating from China and other countries. The tree was constructed using the neighbor-joining algorithm based on the p-distance model. The numbers indicate percentage bootstrap values calculated from 1000 bootstrap replicates. The PRRSV ORF5 sequences collected in this study are marked with a black square. The ORF5 sequences from prototype viruses of PRRSV are marked with a black square

Discussion

Since its emergence in China at the end of 1995, PRRSV has spread widely in all areas of China and is continuously evolving. This has led to the emergence of some new antigenic variant strains [17]. In 2006, a large outbreak of porcine high fever syndrome (PHFD), caused by a highly pathogenic form of PRRSV, emerged in China and Southeast Asian countries and caused major economic losses for swine farming [12, 15, 16] . In our previous study, 133 of 475 samples (28%) were positive for PRRSV, indicating that PRRSV is widely distributed among swine populations in southern parts of China. Fifty-six ORF5 sequences and 35 nsp2 HVR sequences were selected for investigation of variations and phylogenetic analyses for their genetic diversities. Sequences alignments of GP5 and nsp2 showed that there was extensive genetic variability between them (84.5–100% amino acid identity) or with the representative strain, VR-2332 (84.5–99.0% amino acid identity). GP5 based phylogenetic trees showed all these strains belonged to the type 2 PRRSV which are scattered in 2 lineages (lineages 5 and 8). But most of the strains belonged to a large cluster in sub-lineage 8.3 with the representative HP-PRRSV strains being JXA1, JXwn06 and HUN4. This is consistent with other studies showing that the dominant PRRSV seen in Guangdong Province, which is also located in Southern China, was HP-PRRSV between 2007 and 2014 [18, 19]. The sub-lineage 8.3 PRRSV was also the predominant virus at the country-wide scale in the subsequent years since 2007 [17, 20, 21].

As a transmembrane protein, GP5 possesses two to four potential N-linked glycosylation sites that are located in a small ectodomain [22]. The N-linked glycosylation of GP5 have been shown to be involved in diverse functions such as the proper folding protein, receptor binding, virus infectivity and induction of immune response [8, 23,24,25]. The amino acids in the proximal region of the ectodomain of GP5 are highly variable. In this study, we showed that substitutions at each consensus sequence of N-linked glycosylation site, N-X-T/S, in GP5 of some strains, resulted in loss or gain of N-linked glycosylation sites at specific positions in GP5. The N44 and N51-linked glycosylation sites were well conserved. The DCE upstream of the PNE was speculated to elicit a great abundance of the non-neutralized antibodies against GP5 and delay the production of neutralizing antibodies stimulated by PNE [8]. The alignment of GP5 showed that variations in DCE were observed and key aa substitutions in PNE were found among the strains in this study. As a result, a different number of N-glycosylation sites among the strains and the key aa variations in DCE and PNE might allow these field strains to escape neutralization by the antibodies induced by current vaccines.

The results of alignment and phylogenetic tree studies based on GP5 and HV of nsp2 also showed that a great number of emergences of PRRSV might be related to the extensive use of the attenuated modified live virus (MLV) PRRS vaccine in China. Three strains, GXYL1310, GDHZ1401 and GXBH1404, which have a specific 120 aa deletion in nsp2 were identified. It was suggested that there is the possibility that these three strains were derived directly from the widely used commercial vaccine strain, TJM, which is characterized by a 120 aa deletion in nsp2 and derived originally from the TJ strain by serial passage in MARC-145 cells of up to 92 times [26]. VR2332-derived MLV strains were also found. Two strains were clustered into this minor branch shared a high identity with MLV vaccine and its parent virus VR-2332, with amino acid similarities of 99.2 and 99.7%, respectively. Several studies showed that some prevalent PRRSV strains may be related to the reversion of commercial MLVs and the recombination between the vaccine virus and field viruses [27,28,29,30,31]. It has been suggested that more attention should be paid to MLV-like strains which have undergone evolutionary changes and have since circulated widely in the field.

The nsp2 of PRRSV is a highly heterogeneous protein. Remarkably, natural deletions and insertions have continued to occur in the HV of nsp2, and these have led to genome size differences among PRRSV strains [9, 32, 33]. Type 2 PRRSV with 1–150 aa deletions and 1–36 aa insertions in the nsp2 coding region has been identified in USA, China, Japan, Denmark and Thailand [9, 10, 33, 34]. In this study, most of the strains have a discontinuous 30 aa deletion, suggesting that the strain with a 30 aa deletion in nsp2 is the dominant virus prevalent in the southern parts of China. It is noted that several novel strains (GXNN1396, GXYL1403e, GXBS1401a and GXNN1407a) with additional aa deletions or insertions are also found, suggesting that strains with other types of aa deletions in nsp2 may have also been prevalent in this region. GXNN1396 has a discontinuous 30 aa deletion and a continuous 19 aa deletion at position 499–517 in nsp2 compared to VR-2332. Similar 19 aa deletions in nsp2 were also observed in a Japanese PRRSV strain, Jyc, and several USA PRRSV strains (MN184, NADC30 and NADC31), which have 19 aa deletions at position 495–513. GXYL1403e had a novel discontinuous 124 aa deletion at positions 481 and 496–619 in the nsp2-coding region in contrast to the VR2332 reference strain. Similar aa deletions at this region were also observed in a Chinese PRRSV strain, Em2007, which has a 68 aa deletion at position 499–566. Compared to JXA1, strains GXBS1401a and GXNN1407a have one aa deletion at position 816 and one aa insertion at position 830, respectively. One aa insertion in nsp2 was also identified in a Japanese strain, Jtg, which have one aa insertion at position 554 (Additional file 2: Table S2). The emergence of HP-PRRSV strains with 30 aa deletions in the nsp2 was once speculated to be related to its virulence. However, the following study showed that the discontinuous deletion of 30 amino acids in nsp2 was not related to the virulence of the emerging HP- PRRSV [35]. But recent studies showed that PRRSV strains with deletions in the nsp2 were more likely to be pathogenic [33, 35, 36]. Overall, the mechanism underlying spontaneous deletions in nsp2 during viral passages in vivo and their effect on viral replication and pathogenicity remains unclear.

Nsp2 is also a highly immunogenic protein. It has been shown that nsp2 contains several putative B-cell and T-cell epitopes. Antibodies against nsp2 were generated as early as 1 week after PRRSV infection [37, 38]. Most of these epitopes were found to be mapped to the HV of nsp2 which usually occur by substitutions, natural deletions and insertions. It has been shown that a natural deletion or an engineered deletion in nsp2 of PRRSV plays an important role in modulating the induction of inflammatory cytokines in vitro [36, 39, 40]. The biological and immunological characteristics of the strains with specific aa deletions in nsp2 remain topics for further studies.

Conclusions

In this study, we showed PRRSV is widely distributed among swine populations in the southern parts of China. GP5 based phylogenetic trees and sequence alignments showed that extensive genetic variability exists compared with the representative stains and the PRRSV strains with 30 amino acid deletions in nsp2 and these are dominant in the porcine population. In addition, more and more PRRSV strains with different patterns of deletions or insertions in nsp2 are emerging. This study expands the existing knowledge of the genetic diversity and evolution of PRRSV in southern parts of China and can potentially help to better control the spread of PRRSV.

Methods

Sample collection, viral RNA extraction and PRRSV detection

Field samples (n = 475) (sera, lungs, lymph nodes and spleens) from clinically diseased pigs between 2013 and 2014 in different regions of Guangxi, China were submitted to the Laboratory of Animal infectious Diseases and Molecular Immunology, Guangxi University, Nanning for PRRSV testing. A summary of the samples studied is presented in Table 1. Total RNA was extracted using TRizol reagent (Invitrogen, Grand Island, NY, USA) according to the manufacturer’s instructions and then used for synthesis of cDNA with random hexamers (Fermentas, Glen Burnie, MD, USA). All the samples were screened for PRRSV by PCR using the forward and reverse primers, (5′-AAGCTGTTAAACAGGGAGTGG-3′) and (5′-CCAAAGAATACCAGCCCATCA-3′), respectively. Thermal cycling conditions were 95 °C for 3 min, followed by 35 cycles of 95 °C for 1 min, 59 °C for 40s, 72 °C for 1 min, and a final elongation step at 72 °C for 10 min. Finally, the PCR products were analyzed on 1.0% agarose gel electrophoresis ultraviolet imaging. Positive samples were determined by the presence of 443 bp amplified products.

Cells and virus

MARC-145 cells were grown at 37°Cin minimum essential medium (MEM) supplemented with 10% fetal bovine serum (FBS). Sera or supernatants of tissue homogenates from PRRSV-positive samples were used to inoculate the MARC-145 cells for PRRSV isolation.

Amplification of ORF5 and nsp2 hypervariable region and sequence determination

PRRSV positive samples were used for amplification of complete ORF5 and nsp2 hypervariable regions (HVR). The forward (5′-AGGTGGGCAACCGTTTTA-3′) and reverse primers (5′-CATCACTGGCGTGTAGGTAAT-3′) were used for amplification of the complete ORF5. PCR reaction conditions were 95 °C for 3 min, followed by 30 cycles of 95 °C for 1 min, 59 °C for 1 min, 72 °C for 1.5 min, and a final elongation step at 72 °C for 10 min. The forward (5′-AATGTTGTTCTTCCTGGGGTTGAG-3′) and reverse primers (5′-AAGCTGCAAAACCCCAATCACCCG-3′) were used for amplification of the nsp2 HVR. PCR reaction condition was 95 °C for 3 min, followed by 30 cycles of 95 °C for 40 s, 57 °C for 40 s, 72 °C for 2 min, and a final elongation step at 72 °C for 10 min. The PCR products were purified with an E.Z.N.A.TM Gel Extraction Kit (OMEGA, USA) and cloned into pBST-II vector (TIANGEN Inc., Beijing, China). Positive clones were sequenced in both directions using universal primers T7 and SP6 promoter-specific primers.

Amino acid mutation analysis of Nsp2 HV and ORF5

To further characterize the amino acid mutation in Nsp2 HV and GP5, differences of the amino acid sequences derived from ORF5 gene and nsp2 HVR of these strains and other representative strains from China and other countries were analyzed and aligned using the MegAlign program (version 5.01) of the DNASTAR package. (DNASTAR Inc., Madison, WI, USA).

Phylogenetic tree analysis

The multiple sequence alignment of the nucleotide sequences of ORF5 or nsp2 HVR were performed by using the Clustal W method in MEGA5.2. MEGA version 5.2 with the p-distance model was used to evaluate phylogenetic relationships by the neighbor-joining method with 1000 bootstrap replicates. The sequences obtained in this study were submitted to the GenBank database under the accession numbers (MG604994 - MG605049 for ORF5 and MG604959 - MG604993 for nsp2) and the reference strains from China and other countries (lineages 1 to 9) used in this study are listed in Table 2. The classification of lineages and sub-lineages was according to their description in recent studies [41, 42].

Abbreviations

DCE:

decoy epitope

HP:

Highly pathogenic

HVR:

hypervariable region

nsp:

nonstructural proteins

ORF:

Open reading frame

PRRS:

Porcine reproductive and respiratory syndrome

PRRSV:

Porcine reproductive and respiratory syndrome virus

References

  1. Neumann EJ, Kliebenstein JB, Johnson CD, Mabry JW, Bush EJ, Seitzinger AH, Green AL, Zimmerman JJ. Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States. J Am Vet Med Assoc. 2005;227(3):385–92.

    Article  Google Scholar 

  2. Fang Y, Snijder EJ. The PRRSV replicase: exploring the multifunctionality of an intriguing set of nonstructural proteins. Virus Res. 2010;154(1–2):61–76.

    Article  CAS  Google Scholar 

  3. Johnson CR, Griggs TF, Gnanandarajah J, Murtaugh MP. Novel structural protein in porcine reproductive and respiratory syndrome virus encoded by an alternative ORF5 present in all arteriviruses. J Gen Virol. 2011;92(Pt 5:1107–16.

    Article  CAS  Google Scholar 

  4. Wu WH, Fang Y, Farwell R, Steffen-Bien M, Rowland RR, Christopher-Hennings J, Nelson EA. A 10-kDa structural protein of porcine reproductive and respiratory syndrome virus encoded by ORF2b. Virology. 2001;287(1):183–91.

    Article  CAS  Google Scholar 

  5. Snijder EJ, Meulenberg JJ. The molecular biology of arteriviruses. The Journal of general virology. 1998;79(Pt 5):961–79.

    Article  CAS  Google Scholar 

  6. Murtaugh MP, Stadejek T, Abrahante JE, Lam TT, Leung FC. The ever-expanding diversity of porcine reproductive and respiratory syndrome virus. Virus Res. 2010;154(1–2):18–30.

    Article  CAS  Google Scholar 

  7. Nelsen CJ, Murtaugh MP, Faaberg KS. Porcine reproductive and respiratory syndrome virus comparison: divergent evolution on two continents. J Virol. 1999;73(1):270–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Ostrowski M, Galeota JA, Jar AM, Platt KB, Osorio FA, Lopez OJ. Identification of neutralizing and nonneutralizing epitopes in the porcine reproductive and respiratory syndrome virus GP5 ectodomain. J Virol. 2002;76(9):4241–50.

    Article  CAS  Google Scholar 

  9. Yoshii M, Okinaga T, Miyazaki A, Kato K, Ikeda H, Tsunemitsu H. Genetic polymorphism of the nsp2 gene in north American type--porcine reproductive and respiratory syndrome virus. Arch Virol. 2008;153(7):1323–34.

    Article  CAS  Google Scholar 

  10. Kedkovid R, Nuntawan Na Ayudhya S, Amonsin A, Thanawongnuwech R. NSP2 gene variation of the north American genotype of the Thai PRRSV in Central Thailand. Virol J. 2010;7:340.

    Article  Google Scholar 

  11. Huang YW, Meng XJ. Novel strategies and approaches to develop the next generation of vaccines against porcine reproductive and respiratory syndrome virus (PRRSV). Virus Res. 2010;154(1–2):141–9.

    Article  CAS  Google Scholar 

  12. Tian K, Yu X, Zhao T, Feng Y, Cao Z, Wang C, Hu Y, Chen X, Hu D, Tian X, et al. Emergence of fatal PRRSV variants: unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark. PLoS One. 2007;2(6):e526.

    Article  Google Scholar 

  13. Normile D. Virology. China, Vietnam grapple with 'rapidly evolving' pig virus. Science. 2007;317(5841):1017.

    Article  CAS  Google Scholar 

  14. Ni J, Yang S, Bounlom D, Yu X, Zhou Z, Song J, Khamphouth V, Vatthana T, Tian K. Emergence and pathogenicity of highly pathogenic porcine reproductive and respiratory syndrome virus in Vientiane, Lao People's Democratic Republic. J Vet Diagn Invest. 2012;24(2):349–54.

    Article  Google Scholar 

  15. An TQ, Tian ZJ, Xiao Y, Li R, Peng JM, Wei TC, Zhang Y, Zhou YJ, Tong GZ. Origin of highly pathogenic porcine reproductive and respiratory syndrome virus, China. Emerg Infect Dis. 2010;16(2):365–7.

    Article  Google Scholar 

  16. Tong GZ, Zhou YJ, Hao XF, Tian ZJ, An TQ, Qiu HJ. Highly pathogenic porcine reproductive and respiratory syndrome, China. Emerg Infect Dis. 2007;13(9):1434–6.

    Article  Google Scholar 

  17. Zhou L, Chen S, Zhang J, Zeng J, Guo X, Ge X, Zhang D, Yang H. Molecular variation analysis of porcine reproductive and respiratory syndrome virus in China. Virus Res. 2009;145(1):97–105.

    Article  CAS  Google Scholar 

  18. Xie J, Cui T, Cui J, Chen Y, Zhang M, Zhou P, Deng S, Su S, Zhang G. Epidemiological and evolutionary characteristics of the PRRSV in southern China from 2010 to 2013. Microb Pathog. 2014;75:7–15.

    Article  Google Scholar 

  19. Wu Q, Li Z, Zhang G, Niu J, Zeng X, Sun B, Ma J. Genetic diversity and phylogenetic analysis of porcine reproductive and respiratory syndrome virus in southern China from 2007 to 2014. J Vet Sci. 2017;18(3):317–26.

    Article  Google Scholar 

  20. Li B, Fang L, Liu S, Zhao F, Jiang Y, He K, Chen H, Xiao S. The genomic diversity of Chinese porcine reproductive and respiratory syndrome virus isolates from 1996 to 2009. Vet Microbiol. 2010;146(3–4):226–37.

    Article  CAS  Google Scholar 

  21. Zhou L, Yang X, Tian Y, Yin S, Geng G, Ge X, Guo X, Yang H. Genetic diversity analysis of genotype 2 porcine reproductive and respiratory syndrome viruses emerging in recent years in China. Biomed Res Int. 2014;2014:748068.

    PubMed  PubMed Central  Google Scholar 

  22. Dea S, Gagnon CA, Mardassi H, Pirzadeh B, Rogan D. Current knowledge on the structural proteins of porcine reproductive and respiratory syndrome (PRRS) virus: comparison of the north American and European isolates. Arch Virol. 2000;145(4):659–88.

    Article  CAS  Google Scholar 

  23. Ansari IH, Kwon B, Osorio FA, Pattnaik AK. Influence of N-linked glycosylation of porcine reproductive and respiratory syndrome virus GP5 on virus infectivity, antigenicity, and ability to induce neutralizing antibodies. J Virol. 2006;80(8):3994–4004.

    Article  CAS  Google Scholar 

  24. Helenius A, Aebi M. Intracellular functions of N-linked glycans. Science. 2001;291(5512):2364–9.

    Article  CAS  Google Scholar 

  25. 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(18):9941–51.

    Article  CAS  Google Scholar 

  26. Leng X, Li Z, Xia M, Li X, Wang F, Wang W, Zhang X, Wu H. Mutations in the genome of the highly pathogenic porcine reproductive and respiratory syndrome virus potentially related to attenuation. Vet Microbiol. 2012;157(1–2):50–60.

    Article  CAS  Google Scholar 

  27. Lu WH, Tun HM, Sun BL, Mo J, Zhou QF, Deng YX, Xie QM, Bi YZ, Leung FC, Ma JY. Re-emerging of porcine respiratory and reproductive syndrome virus (lineage 3) and increased pathogenicity after genomic recombination with vaccine variant. Vet Microbiol. 2015;175(2–4):332–40.

    Article  CAS  Google Scholar 

  28. Zhao K, Ye C, Chang XB, Jiang CG, Wang SJ, Cai XH, Tong GZ, Tian ZJ, Shi M, An TQ. Importation and recombination are responsible for the latest emergence of highly pathogenic porcine reproductive and respiratory syndrome virus in China. J Virol. 2015;89(20):10712–6.

    Article  CAS  Google Scholar 

  29. Jiang YF, Xia TQ, Zhou YJ, Yu LX, Yang S, Huang QF, Li LW, Gao F, Qu ZH, Tong W, et al. Characterization of three porcine reproductive and respiratory syndrome virus isolates from a single swine farm bearing strong homology to a vaccine strain. Vet Microbiol. 2015;179(3–4):242–9.

    Article  CAS  Google Scholar 

  30. Zhou L, Wang Z, Ding Y, Ge X, Guo X, Yang H. NADC30-like strain of porcine reproductive and respiratory syndrome virus, China. Emerg Infect Dis. 2015;21(12):2256–7.

    Article  CAS  Google Scholar 

  31. Nielsen HS, Oleksiewicz MB, Forsberg R, Stadejek T, Botner A, Storgaard T. Reversion of a live porcine reproductive and respiratory syndrome virus vaccine investigated by parallel mutations. J Gen Virol. 2001;82(Pt 6):1263–72.

    Article  CAS  Google Scholar 

  32. Shen S, Kwang J, Liu W, Liu DX. Determination of the complete nucleotide sequence of a vaccine strain of porcine reproductive and respiratory syndrome virus and identification of the Nsp2 gene with a unique insertion. Arch Virol. 2000;145(5):871–83.

    Article  CAS  Google Scholar 

  33. Brockmeier SL, Loving CL, Vorwald AC, Kehrli ME Jr, Baker RB, Nicholson TL, Lager KM, Miller LC, Faaberg KS. Genomic sequence and virulence comparison of four type 2 porcine reproductive and respiratory syndrome virus strains. Virus Res. 2012;169(1):212–21.

    Article  CAS  Google Scholar 

  34. Kvisgaard LK, Hjulsager CK, Brar MS, Leung FC, Larsen LE. Genetic dissection of complete genomes of type 2 PRRS viruses isolated in Denmark over a period of 15 years. Vet Microbiol. 2013;167(3–4):334–44.

    Article  CAS  Google Scholar 

  35. Zhou L, Zhang J, Zeng J, Yin S, Li Y, Zheng L, Guo X, Ge X, Yang H. The 30-amino-acid deletion in the Nsp2 of highly pathogenic porcine reproductive and respiratory syndrome virus emerging in China is not related to its virulence. J Virol. 2009;83(10):5156–67.

    Article  CAS  Google Scholar 

  36. Choi HW, Nam E, Lee YJ, Noh YH, Lee SC, Yoon IJ, Kim HS, Kang SY, Choi YK, Lee C. Genomic analysis and pathogenic characteristics of type 2 porcine reproductive and respiratory syndrome virus nsp2 deletion strains isolated in Korea. Vet Microbiol. 2014;170(3–4):232–45.

    Article  CAS  Google Scholar 

  37. Yan Y, Guo X, Ge X, Chen Y, Cha Z, Yang H. Monoclonal antibody and porcine antisera recognized B-cell epitopes of Nsp2 protein of a Chinese strain of porcine reproductive and respiratory syndrome virus. Virus Res. 2007;126(1–2):207–15.

    Article  CAS  Google Scholar 

  38. Oleksiewicz MB, Botner A, Toft P, Normann P, Storgaard T. 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. 2001;75(7):3277–90.

    Article  CAS  Google Scholar 

  39. Liu X, Bai J, Wang H, Fan B, Li Y, Jiang P. Effect of amino acids residues 323-433 and 628-747 in Nsp2 of representative porcine reproductive and respiratory syndrome virus strains on inflammatory response in vitro. Virus Res. 2015;208:13–21.

    Article  CAS  Google Scholar 

  40. Faaberg KS, Kehrli ME Jr, Lager KM, Guo B, Han J. In vivo growth of porcine reproductive and respiratory syndrome virus engineered nsp2 deletion mutants. Virus Res. 2010;154(1–2):77–85.

    Article  CAS  Google Scholar 

  41. Shi M, Lam TT, Hon CC, Murtaugh MP, Davies PR, Hui RK, Li J, Wong LT, Yip CW, Jiang JW, et al. Phylogeny-based evolutionary, demographical, and geographical dissection of north American type 2 porcine reproductive and respiratory syndrome viruses. J Virol. 2010;84(17):8700–11.

    Article  CAS  Google Scholar 

  42. Gao JC, Xiong JY, Ye C, Chang XB, Guo JC, Jiang CG, Zhang GH, Tian ZJ, Cai XH, Tong GZ, et al. Genotypic and geographical distribution of porcine reproductive and respiratory syndrome viruses in mainland China in 1996-2016. Vet Microbiol. 2017;208:164–72.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Dev Sooranna, Imperial College London, for editing the manuscript.

Funding

This study was supported by the Natural Sciences Foundation of China (Nos. 31372444 and 31660716), Guangxi Natural Science Foundation (No. 2016GXNSFAA380159 and 2018GXNSFDA281021), Guangxi Science and Technology Bureau (No. AA17204057), and Guangxi University Foundation (No. XGZ130959).

Availability of data and materials

All data were included in the manuscript as tables and figures. The sequences obtained in this study were submitted to the GenBank database under the accession numbers (MG604994 - MG605049 for ORF5 and MG604959 - MG604993 for nsp2) and the reference strains used in this study are listed in Table 2.

Author information

Authors and Affiliations

Authors

Contributions

HS, WY, LS and HJ conducted the experiments and analyzed the data. HW1 and YJ assisted with sample preparation and experiments. CY and KO shared ideas and discussed the research data. HW2 and WZ contributed to supervision, had the idea for the project and directed the research. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Weijian Huang or Zuzhang Wei.

Ethics declarations

Competing interest

The authors declare that they have no competing interests.

Ethics approval and consent to participate

For all porcine clinical samples used in this study, written consents were obtained from farm owners and all procedures were carried out in strict accordance with the Animal Ethics Procedures and Guidelines of the People’s Republic of China. All the animal protocols in this study were approved by the Ethics Committee of Guangxi University.

Consent for publication

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Comparison of the GP5 sequences of the different PRRSV strains examined in this study. (DOCX 14 kb)

Additional file 2:

Table S2. The positions and sizes of aa insertions and deletions in nsp2 of PRRSV strains compared to VR2332. (DOCX 26 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hong, S., Wei, Y., Lin, S. et al. Genetic analysis of porcine productive and respiratory syndrome virus between 2013 and 2014 in Southern parts of China: identification of several novel strains with amino acid deletions or insertions in nsp2. BMC Vet Res 15, 171 (2019). https://doi.org/10.1186/s12917-019-1906-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12917-019-1906-9

Keywords