ORIGINAL RESEARCH article

Front. Vet. Sci., 10 September 2020

Sec. Veterinary Epidemiology and Economics

Volume 7 - 2020 | https://doi.org/10.3389/fvets.2020.00603

Full Genomic Analysis of New Variants of Porcine Reproductive and Respiratory Syndrome Virus Revealed Multiple Recombination Events Between Different Lineages and Sublineages

  • JW

    Jinglong Wang

  • SL

    Siyuan Lin

  • DQ

    Dongqun Quan

  • HW

    Hao Wang

  • JH

    Jiabin Huang

  • YW

    Yuxu Wang

  • TR

    Tongwei Ren

  • KO

    Kang Ouyang

  • YC

    Ying Chen

  • WH

    Weijian Huang

  • TL

    Tingrong Luo *

  • ZW

    Zuzhang Wei *

  • Laboratory of Animal infectious Diseases and Molecular Immunology, College of Animal Science and Technology, Guangxi University, Nanning, China

Abstract

Porcine reproductive and respiratory syndrome virus (PRRSV) has had a devastating impact on the pig industry in China, and monitoring its genetic diversity is important for epidemiological surveillance and understanding its evolution. Here, we determine the complete genome sequences of two PRRSV strains, GXYL1403 and GXNN1839. Comparative, phylogenetic, and recombination detection program analyses show that the two isolates are recombinant strains with large-fragment amino acid deletions in nsp2. GXYL1403 possesses a unique deletion region of 124 amino acids in nsp2, and GXNN1839 contains a deletion of 131 amino acids in nsp2 as compared with VR2332. Further analysis of the full-length sequence suggests that GXYL1403 is a natural recombinant between sublineages 8.1 (CH-1a like) and 8.3 (JXA1-like). The recombination site of GXYL1403 is located in nsp9–nsp12 (8961nt−11181nt). GXNN1839 is a natural recombinant between the lineage 5 (VR-2332-like) and lineage 1 (NADC30-like) strains. The recombination events occurred in nsp9 (7872nt-8162nt) and in ORF2 (12587nt−13282nt) in the genome of GXNN1839. These results provide new evidence that PRRSV strains circulating in the environment have undergone recombination among the different lineages or sublineages of field strains, and these add to our understanding of RNA combination events that occur in PRRSV.

Introduction

Porcine reproductive and respiratory syndrome (PRRS) is an important infectious disease caused by porcine reproductive and respiratory syndrome virus (PRRSV), and it poses a serious threat to the global pig industry (1). PRRSV is a small, enveloped, positive-stranded RNA virus whose genome consists of at least 10 open reading frames (ORFs), namely ORF1a, ORF1b, ORF2a, ORF2b, ORF3, ORF4, ORF5, ORF5a, ORF6, and ORF7 (2). ORF1a and ORF1b are situated at the 5′ end of the genome and encode polyproteins that are self-cleaved to at least 16 non-structural proteins (nsps) by viral proteases (2, 3). ORF2a, ORF2b, ORF3, ORF4, ORF5, ORF5a, ORF6, and ORF7 are situated at the 3′ end of the genome and encode eight structural proteins, namely glycoprotein (GP) 2a, 2b, GP3, GP4, GP5, ORF5a protein, matrix protein (M), and nucleocapsid (N), respectively (3–5). PRRSV is divided into type 1 and 2 PRRSV genotypes, which only have about 60% nucleotide sequence identity (6, 7). Phylogenetic analyses of the ORF5 of type 2 PRRSV strains show that type 2 PRRSV could be divided into nine distinct lineages with each lineage containing several sublineages (8, 9). Type 2 PRRSV strains are prevalent in China and are distributed into four lineages and sublineages: lineages 1, 3, 5, and 8 are representatives of the NADC-30, QYYZ-like, VR-2332, and JXA1 strains, respectively (9). In 2006, a devastating outbreak of PRRS, caused by a highly pathogenic PRRSV (HP-PRRSV) representative with the JXA1 strain (sublineage 8.3) emerged in China, and this caused substantial economic losses for swine farming (10, 11). The lineage 3 PRRSV represented by the QYYZ strain emerged in 2010, and it is mainly circulating in southern China (12, 13). The NADC30-like PRRSV as represented by the WUH6 strain was first isolated in 2011 and is mainly found in northeast and central China (14–17).

The increasing diversity of the type 2 PRRSV is related to point mutations and genetic recombination among the field strains or in between the field and modified live vaccine (MLV) strains, which can lead to some new emerging antigenic variant strains (6, 18–24). In the present study, we identify two novel recombinant PRRSV strains circulating in swine farms in Guangxi, a province of southern China. We perform genome-wide sequence analysis of the two isolates and find new combination patterns in PRRSV evolution.

Materials and Methods

Cell Culture and Virus Isolation

Meat Animal Research Center-145 (MARC-145) cells and porcine alveolar macrophages (PAMs) were maintained in cell cultures as described in our previous study (25). Field samples (sera) were collected from clinically diseased pigs in Guangxi province, China. Samples were detected by RT-PCR using specific primers for the ORF5 gene of genotype 2 PRRSV (26). One positive sample was from 56 dying piglets exhibiting severe respiratory syndrome on a pig farm located in YuLing city, Guangxi Province, which has a total pig population of 823. One positive sample was from 18 pregnant sows suffering reproductive failure in a farm located in Nanning city, Guangxi Province in South China, which has a total sow population of 312. The positive samples were used for virus isolation using MARC-145 cells or PAMs. Cells were observed daily for the appearance of cytopathic effects (CPE). When 75% CPE was observed, the culture media were harvested and used for blinding passage in MARC-145 cells or PAMs. The PRRSV strain (GXYL1403) isolated from MARC-145 cells was then plaque-purified three times on MARC-145 cells for subsequent complete genome sequencing. The third passage of each virus strain was identified by RT-PCR and an indirect immunofluorescence assay (IFA).

RT-PCR Amplification and Complete Genome Determination

The viral genomic RNA was extracted from PRRSV-infected PAMs or MARC-145 cells and reverse transcribed into cDNA using M-MLV reverse transcriptase (Takara) according to the manufacturer's instructions. PCR was performed to amplify the complete genomic sequences of PRRSV isolates GXYL1403 and GXNN1839 by using PRRSV-specific primers as described previously (27, 28). The PCR products were purified using an E.Z.N.A.TM Gel Extraction Kit (OMEGA, USA) and cloned into a pMD18-T vector (TaKaRa, Dalian, China) for nucleotide sequencing. The complete genomic sequences of PRRSV isolates were assembled using the SeqMan program of DNAstar software, version 7.0.

Immunofluorescence Assay (IFA)

IFA was used to detect PRRSV N protein expression in PAMs and MARC-145 cells as described previously (25). At 72 h infection, the isolate-infected PAMs or MARC-145 cells were washed with phosphate-buffered saline (PBS) twice and fixed using cold methanol for 15 min at room temperature and then blocked with 3% BSA (fraction V bovine serum albumin) (Roche, Mannheim, Germany) for 30 min. The cells were treated with a monoclonal antibody (mAb) against PRRSV N protein (SR30-A) (Rural Technologies, Inc., Brookings, SD, USA) as the primary antibody and incubated for 2 h at 37°C. Then, the cells were washed five times with PBS and incubated with goat antirabbit IgG H&L (Alexa Fluor® 488, Abcam) for 1 h at 37°C. Following five further washes with PBS, the cells were observed under a fluorescence microscope.

Sequence Comparison and Phylogenetic Analysis

Differences of sequences from the two isolates and reference strains from all over the world were analyzed and aligned using the MegAlign program in DNAstar 7.0 software (DNASTAR Inc., Madison, WI, USA). Phylogenetic trees, including the complete genomes and ORF5 from this study and the representative strains from China and other countries, were constructed by neighbor-joining in MEGA 6.0, using the maximum composite likelihood and bootstrap confidence values from 1,000 replicates. Detailed information on the PRRSV reference sequences is shown in Table 1.

Table 1

No.Virus strainAccession no.CountryLineagePhylogenetic tree
1NADC30JN654459American1GP5, Genome
2MN184EF442777American1GP5, Genome
3CHsx1401KP861625China1GP5, Genome
4HNhxKX766379China1GP5, Genome
5PRRSV0000008659EU758687American2GP5
6PRRSV0000008973EU758940American2GP5
7PRRSV0000031DQ474791American2GP5
8GD-KPKU978619China3GP5, Genome
9QYYZJQ308798China3GP5, Genome
10GM2JN662424China3GP5, Genome
11Ibaraki08–5AB546113Japan4GP5
12Miyagi08–2AB546105Japan4GP5
13Miyagi08–3AB546106Japan4GP5
14EDRD-1AB288356Japan4GP5, Genome
15VR-23332AY150564American5GP5, Genome
16NADC-8AF396833American5GP5
17PA8AH006184Canada5GP5
18NVSL-14AF396839American6GP5
19Aichi N20AB175715Japan7GP5
20Neb-1EU755263American7GP5, Genome
21PrimePacAF066384American7GP5, Genome
22CH-1aAY032626China8.1GP5, Genome
23CH-1REU807840China8.1GP5, Genome
24HBJM2EU399826China8.2GP5
25HBSZEU399825China8.2GP5
26JXA1EF112445China8.3GP5, Genome
27JXA1 P80FJ548853China8.3GP5, Genome
28JXwn06EF641008China8.3GP5, Genome
29TJEU860248China8.3GP5, Genome
30JXZX2HQ832215China8.4GP5
31HeN-2FJ237419China8.5GP5
32AHW01EU399828China8.5GP5
33GXLSN06–2012KC618172China8.6GP5
34Yunnan-08EU819086China8.6GP5
35HK1KF287132China8.7GP5, Genome
36HK4KF287134China8.7GP5, Genome
37JA-142AF396842American9GP5, Genome
38SDSU73AY656993American9GP5, Genome
39LVM96262NetherlandsGP5, Genome

Information regarding the reference strains used in this study.

Recombination Analyses

The complete genomic sequences of GXYL1403 and GXNN1839 were compared with PRRSV reference strains JXA1, CH1a, VR-2332, and NADC30 using the Recombination Detection Program (RDP) v4.66 with six different algorithms (RDP, GENECONV, MaxChi, BootScan, SiScan, and 3Seq) (29). The results were presented using the RDP method. Recombinant events were confirmed by Bootscan analysis in Simplot software (v3.5.1, JHK University, Baltimore, MD, USA) with the default parameters. A series of phylogenetic trees based on each of the sequence regions were constructed to further identify these putative recombination events.

Results

Isolation and Identification of PRRSV Strains

Two PRRSV-positive sera were used to inoculate the MARC-145 cells and PAMs for PRRSV isolation. Typical CPE, characterized by cell fusion and shedding, was induced in both MARC-145 cells and PAMs 48 h post-inoculation (hpi) (data not shown). To confirm the isolation of the PRRSV strains, IFA was conducted using a mAb against PRRSV N (SR-30A). PAMs and MARC-145 cells infected with PRRSV isolates reacted with the specific antibody against the PRRSV N protein (data not shown). These results indicate successful isolation of the infectious PRRSV strains. These PRRSV isolates were designated as GXYL1403 and GXNN1839, respectively.

Genomic Characteristics and Phylogenetic Analysis of the Two PRRSV Isolates

The full-length genomes of GXYL1403 and GXNN1839 were determined and submitted to the GenBank database with accession numbers MN660069 and MN660070. The entire genomes of GXYL1403 and GXNN1839 are 15,018 and 15,019 nucleotides in length, respectively, excluding the poly (A) tails as shown in Tables 2, 3. The nucleotide similarities of GXYL1403 with PRRSV representative strains NADC30, VR-2332, QYYZ, CH-1a, and JXA1 were 82.5, 88.8, 87.2, 94.2, and 98.0%, respectively, with only 60.6% nucleotide similarity with the type 1 PRRSV LV strain (Table 2). The nucleotide similarities of GXNN1839 with NADC30, VR-2332, QYYZ, CH-1a, and JXA1 were 92.2, 84.5, 81.9, 83.6, and 82.9%, respectively, with only 59.8% nucleotide similarity with the LV strain (Table 3). This indicates that the GXYL1403 and GXNN1839 strains belong to type 2 PRRSV. Each region of the GXYL1403 and GXNN1839 genomes were further compared with those of the reference strains as shown in Tables 2, 3. The results show that GXYL1403 exhibits the highest nucleotide similarity (92.8–100%) and amino acid similarity (93.4–100%) with JXA1. The 5′UTR and nsp8-nsp12 of GXYL1403 share the highest nucleotide similarity with that of CH-1a. The other regions of GXYL1403 exhibit the highest nucleotide and amino acid similarities with those of JXA1. GXNN1839 shares the highest nucleotide (86.4–97.9%) and amino acid (83.9–98.9%) similarities with the NADC30 strain (Table 3). ORF2–ORF4 of GXNN1839 share higher (92.6–96.5%) nucleotide and amino acid similarities (93–96.9%) with VR-2332 but exhibit lower nucleotide identity with the other reference strains. The other regions of the GXNN1839 genome exhibit the highest nucleotide and amino acid similarities with NADC30. These results indicate that GXYL1403 and GXNN1839 may be made up of mosaic isolates.

Table 2

Genomic regionPairwise % identity (nt/aa)
NADC30VR-2332QYYZCH-1aJXA1Lelystad virus
Complete genome82.588.887.394.298.060.6
5′UTR90.692.695.894.990.862.2
Nsp183.8/85.187.9/89.087.7/89.692.6/91.496.8/96.157.6/53.0
Nsp268.6/58.584.0/79.083.4/78.690.8/86.497.0/96.134.8/29.5
Nsp382.5/90.490.7/94.682.1/88.695.1/98.298.9/99.159.6/58.9
Nsp483.5/92.689.7/93.184.2/91.795.4/95.699.0/98.562.7/63.1
Nsp586.5/91.289.4/92.482.4/90.094.3/94.798.8/10064.0/71.8
Nsp693.8/10095.8/93.897.9/10097.9/100100/10068.8/75.0
Nsp781.6/83.488.5/89.292.8/94.695.0/96.198.6/99.257.7/46.5
Nsp888.1/91.295.6/97.894.1/97.899.3/97.898.5/97.863.7/66.7
Nsp986.2/93.191.5/97.790.8/97.095.4/98.197.9/99.167.2/75.0
Nsp1083.7/87.289.3/96.589.4/97.093.1/97.796.7/98.960.4/64.5
Nsp1189.3/96.489.4/95.587.5/94.696.2/98.794.5/93.466.1/77.2
Nsp1289.4/71.988.9/70.686.8/62.795.4/88.892.8/93.644.1/45.8
ORF284.0/83.191.3/92.089.0/88.794.9/97.997.6/96.465.4/62.4
ORF382.6/80.889.0/87.190.7/87.895.9/93.399.1/97.665.2/57.1
ORF485.3/87.189.9/91.193.7/95.596.5/98.396.5/96.666.3/72.1
ORF584.1/84.588.4/87.183.1/80.694.5/91.599.2/98.564.4/57.4
ORF689.0/93.795.0/97.790.5/96.697.1/97.799.6/10069.2/80.5
ORF790.3/90.192.7/93.589.0/91.196.5/96.099.2/99.269.2/63.6
3′UTR8992.690.092.795.265.4

Nucleotide and amino acid identities of different regions of GXYL1403 compared with other PRRSV reference strains.

Table 3

Genomic regionPairwise % identity (nt/aa)
NADC30VR-2332QYYZCH-1aJXA1Lelystad virus
Complete genome92.284.581.983.982.959.8
5′UTR93.692.687.994.990.862.4
Nsp192.0/90.183.9/85.480.5/80.785.6/85.18771.7/81.556.9/53.8
Nsp289.6/83.976.6/68.672.7/63.377.7/69.973.2/64.633.3/29.9
Nsp390.9/92.884.3/94.180.0/86.383.3/90.881.2/93.660.2/57.5
Nsp493.6/97.585.5/93.183.7/91.784.2/92.284.6/93.661.2/61.6
Nsp592.2/92.987.5/89.481.6/89.486.2/92.985.5/91.863.5/70.0
Nsp689.6/93.887.5/10085.4/93.889.6/90.887.5/93.872.9/81.2
Nsp792.3/90.384.0/85.779.5/83.482.8/85.780.4/83.455.2/48.0
Nsp891.1/98.988.1/88.983.7/88.984.4/88.985.9/86.760.0/68.9
Nsp992.9/96.987.8/96.685.1/95.286.6/86.886.4/88.967.8/74.2
Nsp1093.8/97.785.6/94.883.9/93.484.7/85.584.2/86.960.8/64.5
Nsp1194.6/96.087.2/95.186.0/94.287.7/86.489.4/92.965.8/76.3
Nsp1293.7/83.787.0/90.284.2/90.888.0/90.888.0/92.244.3/42.7
ORF291.2/92.496.5/96.983.7/84.885.5/87.284.4/85.667.5/60.4
ORF392.6/92.794.8/93.085.1/82.787.5/79.686.0/79.266.0/57.5
ORF486.40/87.892.6/93.283.4/84.985.8/89.484.0/89.467.8/68.7
ORF593.0/95.087.0/85.684.8/83.086.7/87.686.7/87.661.9/53.3
ORF697.7/96.189.9/91.588.8/92.088.4/88.089.3/88.170.782.8
ORF794.9/97.490.9/92.486.6/88.790.1/93.589./91.965.8/60.3
3′UTR9892.685.787.889.164.9

Nucleotide and amino acid identities of different regions of GXNN1839 compared with other type 2 PRRSV reference strains.

Nsp2 is one of the most heterogeneous proteins in PRRSV and contains different patterns of amino acid insertions and deletions. Compared with VR-2332, GXYL1403 contains a discontinuous 124-amino acid deletion in the hypervariable region (HVR) of nsp2. GXNN1839 has a discontinuous 131-amino acid deletion in the HVR of nsp2, which has the same amino acid deletion pattern as the NADC30-like strain (Figure 1). This data suggest that novel PRRSV strains with amino acid deletions in nsp2 are emerging.

Figure 1

The phylogenetic tree based on GP5 shows that the type 2 PRRSV strains can be classified into nine different lineages. GXYL1403 belongs to sublineage 8.3, which is represented by the reference strain JXA1. GXNN1839 belongs to lineage 3, which is represented by the reference strain NADC30 (Figure 2A). The phylogenetic tree based on the complete genome also shows that GXYL1403 belongs to sublineage 8.3, and GXNN1839 belongs to lineage 3 (Figure 2B).

Figure 2

Recombination Analysis

To determine the recombinant events in the genomes of the GXYL1403 and GXNN1839 isolates, we performed recombinant detection using RDP4 and SimPlot by comparing GXYL1403 and GXNN1839 with the representative PRRSV strains from different lineages in the phylogenetic tree, including strains from lineage 1 (NADC30), lineage 3 (QYYZ), lineage 5 (VR2332), and lineage 8 (CH-1a and JX1A). The results show that putatively recombinant breakpoints located at positions 8961nt−11181nt in the GXYL1403 genome were detected with a high degree of reliability (Table 4). The breakpoints in the GXYL1403 genome separated its genome into three regions (Figure 3A). The phylogenetic trees also paralleled and confirmed the recombinant events as shown in Figure 3B. The region 1nt−8961nt is closely genetically related to the JXA1-like strain. The region 8961nt−11181nt is closely genetically related to the CH-1a-like strain, and the region 11182nt−5018nt is closely related to the JXA1-like strain (Figure 3B). It was suggested that GXYL1403 is a natural recombinant between sublineage 8.1 (CH-1a-like) and sublineage 8.3 (JXA1-like). Two putatively major recombinant breakpoints in the GXNN1839 genome were detected with high reliability, and they are located at the nsp9 coding region (7872nt−8162nt) and the ORF2 coding region (12587nt−13282nt), respectively (Table 4). The breakpoints in the GXNN1839 genome separated its genome into five regions (Figure 4A). The phylogenetic trees also paralleled and confirmed the recombinant events as shown in Figure 4B. This suggests that GXNN1839 is a natural recombinant between linage 5 (VR-2332-like) and lineage 1 (NADC30-like).

Table 4

IsolateBreakpointsMajor parentMinor parentp-Value of the detection methods
BeginningEndingRDPGENECONVBootScanMaxChiChimeraSiScan3Seq
GXNN18397,8728,162NADC30VR-23323.373 × 10−113.971 × 10−132.189 × 10−171.190 × 10−142.497 × 10−153.534 × 10−161.545 × 10−14
12,58713,282NADC30VR-23325.253 × 10−266.773 × 10−211.725 × 10−211.316 × 10−111.655 × 10−122.546 × 10−182.22 × 10−15
GXYL14038,96111,181JXA1CH-1a3.408 × 10−186.05 × 10−232.648 × 10−161.146 × 10−242.228 × 10−173.782 × 10−141.515 × 10−17

Recombination breakpoints identified in GXYL1403, GXNN1839, and associated parental strains of PRRSV.

Figure 3

Figure 4

Discussion

PRRSV has had a large impact on the Chinese swine industry since its initial outbreak in China in 1996. In 2006, a large outbreak of PRRS caused by HP-PRRSV (10, 11, 26) emerged in most areas of China and caused major economic losses for swine farming. With the emergence of lineages 3 and 1 PRRSV in some areas in China in 2010 and 2015 (12, 14–17), respectively, the genetic diversity of RPRSV has significantly expanded and has further complicated vaccine development. In the present study, we isolated two PRRSV strains (GXYL1403 and GXNN1839) and determined their complete genome sequences. Consistent with previous studies, phylogenetic trees based on GP5 show type 2 PRRSV strains isolated in China are divided into four lineages (1, 3, 5, and 8) (9). The GXYL1403 and GXNN1839 strains reported in this study are grouped into lineages 1 and 3, respectively. The nsp2 is the most heterogeneous protein in PRRSV (26, 30, 31). Compared to strain VR-2332, the nsp2 of the two isolates in this study contain different amino acid deletions, which suggest that novel PRRSV strains with amino acid deletions in nsp2 are emerging.

Recombination plays a critical role in the genetic evolution of PRRSV and presents great challenges for the prevention and control of the disease it causes. Recombination can occur during the virus-replication process. It has been shown that there are different patterns of recombination among the different lineages and sublineages of PRRSV strains (14, 32–39). In the present study, whole genome sequence comparison and recombination analysis reveals that the two isolates (GXYL1403 and GXNN1839) show novel recombination patterns. GXNN1839 is a recombinant virus that originated from recombination events that occurred in nsp9 and GP2 between strains of lineages 1 (NADC30-like) and 5 (VR2332-like). Many recombinant PRRSV strains that originated from recombination events between NADC30-like and VR2332-like strains have also been identified in previous studies (23, 32, 38, 40). Indeed, the NADC30-like strains seem to be prone to recombine with other PRRSV strains. Numerous recombinant PRRSV strains originated from recombination events between NADC30-like strains and other lineages of PRRSV strains that were found in the field (14, 23, 24, 34, 40–42). GXYL1403 was detected as a recombinant virus from strains of sublineages 8.1 (Ch-1a-like) and 8.3 (JX1A-like) that have circulated in China in recent years. Recombinant PRRSV strains that originated from recombination events that occurred in nsp9 between sublineages 8.1 (Ch-1a-like) and 8.3 (JX1A-like) are also reported in a previous study (42). The emergence of new PRRSV recombination variants poses a major obstacle to the effective control of the spread of virus infection and complicates vaccine development for adequate preventative measures (43). The exact mechanism of genetic recombination in PRRSV remains unknown. The breakpoints of recombination sites in the genome among different PRRS strains can be found at each region of the PRRSV genome and seem to appear randomly. A recently published study shows that nsp9 and GP2–GP3 are the hot spots for PRRSV RNA recombination (44). In this study, we find the two strains (GXYL1403 and GXNN1839) both display a breakpoint in nsp9, and GXNN1839 has an additional breakpoint in GP2, indicating that PRRSV gains genetic diversity through the frequency of recombination events at specific regions with PRRSV strains located in China.

In summary, we determined the complete genome sequences of two novel PRRSV isolates (GXYL1403 and GXNN1839). Both of these are natural recombinant strains that contain amino acid deletions in nsp2. Our findings suggest that PRRSV strains circulating in southern China have undergone recombination among the different lineages or sublineages of field strains, and RNA recombination contributes to their genetic diversity.

Statements

Data availability statement

The datasets generated for this study can be found in the nucleotide sequences generated in present study were submitted to the NCBI GenBank with accession numbers, MN660069 and MN660070.

Author contributions

ZW, TL, and WH conceived and designed the experiments. JW, SL, DQ, WH, and JH performed the experiments. TR, YC, and KO analyzed the data. JW wrote the manuscript. ZW and TL revised the manuscript. All the authors read and approved the final version of the manuscript.

Funding

This study was supported by the National Key Research and Development Program of China (No. 2018YFD0500100), Natural Sciences Foundation of China (No. 31660716), Natural Science Foundation of Guangxi province (No. 2018GXNSFDA281021), and Guangxi Science and Technology Bureau (No. AA17204057).

Acknowledgments

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

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    NeumannEJKliebensteinJBJohnsonCDMabryJWBushEJSeitzingerAHet al. 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:385–92. 10.2460/javma.2005.227.385

  • 2.

    FangYSnijderEJ. The PRRSV replicase: exploring the multifunctionality of an intriguing set of nonstructural proteins. Virus Res. (2010) 154:61–76. 10.1016/j.virusres.2010.07.030

  • 3.

    SnijderEJMeulenbergJJ. The molecular biology of arteriviruses. J Gen Virol. (1998) 79:961–79. 10.1099/0022-1317-79-5-961

  • 4.

    JohnsonCRGriggsTFGnanandarajahJMurtaughMP. Novel structural protein in porcine reproductive and respiratory syndrome virus encoded by an alternative ORF5 present in all arteriviruses. J Gen Virol. (2011) 92:1107–16. 10.1099/vir.0.030213-0

  • 5.

    WuWHFangYFarwellRSteffen-BienMRowlandRRChristopher-HenningsJet al. A 10-kDa structural protein of porcine reproductive and respiratory syndrome virus encoded by ORF2b. Virology. (2001) 287:183–91. 10.1006/viro.2001.1034

  • 6.

    MurtaughMPStadejekTAbrahanteJELamTTLeungFC. The ever-expanding diversity of porcine reproductive and respiratory syndrome virus. Virus Res. (2010) 154:18–30. 10.1016/j.virusres.2010.08.015

  • 7.

    NelsenCJMurtaughMPFaabergKS. Porcine reproductive and respiratory syndrome virus comparison: divergent evolution on two continents. J Virol. (1999) 73:270–80. 10.1128/JVI.73.1.270-280.1999

  • 8.

    ShiMLamTTHonCCMurtaughMPDaviesPRHuiRKet al. Phylogeny-based evolutionary, demographical, and geographical dissection of North American type 2 porcine reproductive and respiratory syndrome viruses. J Virol. (2010) 84:8700–11. 10.1128/JVI.02551-09

  • 9.

    GaoJCXiongJYYeCChangXBGuoJCJiangCGet al. Genotypic and geographical distribution of porcine reproductive and respiratory syndrome viruses in mainland China in 1996–2016. Vet Microbiol. (2017) 208:164–72. 10.1016/j.vetmic.2017.08.003

  • 10.

    TianKYuXZhaoTFengYCaoZWangCet al. Emergence of fatal PRRSV variants: unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark. PLoS ONE. (2007) 2:e526. 10.1371/journal.pone.0000526

  • 11.

    TongGZZhouYJHaoXFTianZJAnTQQiuHJ. Highly pathogenic porcine reproductive and respiratory syndrome, China. Emerg Infect Dis. (2007) 13:1434–36. 10.3201/eid1309.070399

  • 12.

    LuWHTunHMSunBLMoJZhouQFDengYXet al. Re-emerging of porcine respiratory and reproductive syndrome virus (lineage 3) and increased pathogenicity after genomic recombination with vaccine variant. Vet Microbiol. (2015) 175:332–40. 10.1016/j.vetmic.2014.11.016

  • 13.

    WenhuiLZhongyanWGuanqunZZhiliLJingYunMQingmeiXet al. Complete genome sequence of a novel variant porcine reproductive and respiratory syndrome virus (PRRSV) strain: evidence for recombination between vaccine and wild-type PRRSV strains. J Virol. (2012) 86:9543. 10.1128/JVI.01341-12

  • 14.

    ZhaoKYeCChangXBJiangCGWangSJCaiXHet al. Importation and recombination are responsible for the latest emergence of highly pathogenic porcine reproductive and respiratory syndrome virus in China. J Virol. (2015) 89:10712–6. 10.1128/JVI.01446-15

  • 15.

    ZhouLWangZDingYGeXGuoXYangH. NADC30-like strain of porcine reproductive and respiratory syndrome virus, China. Emerg Infect Dis. (2015) 21:2256–7. 10.3201/eid2112.150360

  • 16.

    ZhangQJiangPSongZLvLLiLBaiJ. Pathogenicity and antigenicity of a novel NADC30-like strain of porcine reproductive and respiratory syndrome virus emerged in China. Vet Microbiol. (2016) 197:93–101. 10.1016/j.vetmic.2016.11.010

  • 17.

    LiCZhuangJWangJHanLSunZXiaoYet al. Outbreak Investigation of NADC30-Like PRRSV in South-East China. Transbound Emerg Dis. (2016) 63:474–9. 10.1111/tbed.12530

  • 18.

    MartonSSzalayDKecskemetiSForroBOlaszFZadoriZet al. Coding-complete sequence of a vaccine-derived recombinant porcine reproductive and respiratory syndrome virus strain isolated in Hungary. Arch Virol. (2019) 164:2605–8. 10.1007/s00705-019-04338-2

  • 19.

    LiuJWeiCLinZFanJXiaWDaiAet al. Recombination in lineage 1, 3, 5 and 8 of porcine reproductive and respiratory syndrome viruses in China. Infect Genet Evol. (2019) 68:119–26. 10.1016/j.meegid.2018.12.006

  • 20.

    ZhangQXuXYouSLiYWangHBaiJet al. Emerging of two new subgenotypes of porcine reproductive and respiratory syndrome viruses in Southeast China. Microb Pathog. (2016) 97:27–33. 10.1016/j.micpath.2016.05.011

  • 21.

    SunYKLiQYuZQHanXLWeiYFJiCHet al. Emergence of novel recombination lineage 3 of porcine reproductive and respiratory syndrome viruses in Southern China. Transbound Emerg Dis. (2019) 66:578–87. 10.1111/tbed.13067

  • 22.

    ShiMHolmesECBrarMSLeungFC. Recombination is associated with an outbreak of novel highly pathogenic porcine reproductive and respiratory syndrome viruses in China. J Virol. (2013) 87:10904–7. 10.1128/JVI.01270-13

  • 23.

    LiuYLiJYangJZengHGuoLRenSet al. Emergence of different recombinant porcine reproductive and respiratory syndrome viruses, China. Sci Rep. (2018) 8:4118. 10.1038/s41598-018-22494-4

  • 24.

    ZhouLKangRZhangYDingMXieBTianYet al. Whole genome analysis of two novel type 2 porcine reproductive and respiratory syndrome viruses with complex genome recombination between lineage 8, 3, and 1 strains identified in Southwestern China. Viruses. (2018) 10:328. 10.3390/v10060328

  • 25.

    WeiZLinTSunLLiYWangXGaoFet al. N-linked glycosylation of GP5 of porcine reproductive and respiratory syndrome virus is critically important for virus replication in vivo. J Virol. (2012) 86:9941–51. 10.1128/JVI.07067-11

  • 26.

    HongSWeiYLinSHuangJHeWYaoJet 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. (2019) 15:171. 10.1186/s12917-019-1906-9

  • 27.

    LengXLiZXiaMLiXWangFWangWet al. Mutations in the genome of the highly pathogenic porcine reproductive and respiratory syndrome virus potentially related to attenuation. Vet Microbiol. (2012) 157:50–60. 10.1016/j.vetmic.2011.12.012

  • 28.

    ZhangHLengCDingYZhaiHLiZXiangLet al. Characterization of newly emerged NADC30-like strains of porcine reproductive and respiratory syndrome virus in China. Arch Virol. (2019) 164:401–11. 10.1007/s00705-018-4080-7

  • 29.

    MartinDPLemeyPLottMMoultonVPosadaDLefeuvreP. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. (2010) 26:2462–3. 10.1093/bioinformatics/btq467

  • 30.

    KedkovidRNuntawan Na AyudhyaSAmonsinAThanawongnuwechR. NSP2 gene variation of the North American genotype of the Thai PRRSV in central Thailand. Virol J. (2010) 7:340. 10.1186/1743-422X-7-340

  • 31.

    YoshiiMOkinagaTMiyazakiAKatoKIkedaHTsunemitsuH. Genetic polymorphism of the nsp2 gene in North American type–porcine reproductive and respiratory syndrome virus. Arch Virol. (2008) 153:1323–34. 10.1007/s00705-008-0098-6

  • 32.

    LiYJiGWangJTanFZhuangJLiXet al. Complete genome sequence of an NADC30-like porcine reproductive and respiratory syndrome virus characterized by recombination with other strains. Genome Announc. (2016) 4:e00330–16. 10.1128/genomeA.00330-16

  • 33.

    ZhangZZhouLGeXGuoXHanJYangH. Evolutionary analysis of six isolates of porcine reproductive and respiratory syndrome virus from a single pig farm: MLV-evolved and recombinant viruses. Infect Genet Evol. (2018) 66:111–9. 10.1016/j.meegid.2018.09.024

  • 34.

    ZhouLKangRJiGTianYGeMXieBet al. Molecular characterization and recombination analysis of porcine reproductive and respiratory syndrome virus emerged in southwestern China during 2012–2016. Virus Genes. (2018) 54:98–110. 10.1007/s11262-017-1519-y

  • 35.

    ZhaoHHanQZhangLZhangZWuYShenHet al. Emergence of mosaic recombinant strains potentially associated with vaccine JXA1-R and predominant circulating strains of porcine reproductive and respiratory syndrome virus in different provinces of China. Virol J. (2017) 14:67. 10.1186/s12985-017-0735-3

  • 36.

    DongJWangYYuLZhangPLiuXZhangLet al. Pathogenicity of a newly emerged recombined porcine reproductive and respiratory syndrome virus strain (subgenotype III) in China. Vet Microbiol. (2017) 210:162–6. 10.1016/j.vetmic.2017.09.018

  • 37.

    HanGXuHWangKHeF. Emergence of Two different recombinant PRRSV strains with low neutralizing antibody susceptibility in China. Sci Rep. (2019) 9:2490. 10.1038/s41598-019-39059-8

  • 38.

    ZhouLKangRXieBTianYWuXLvXet al. Identification of a novel recombinant type 2 porcine reproductive and respiratory syndrome virus in China. Viruses. (2018) 10:151. 10.3390/v10040151

  • 39.

    ChenNYeMLiSHuangYZhouRYuXet al. Emergence of a novel highly pathogenic recombinant virus from three lineages of porcine reproductive and respiratory syndrome virus 2 in China 2017. Transbound Emerg Dis. (2018) 65:1775–85. 10.1111/tbed.12952

  • 40.

    ZhouLKangRYuJXieBChenCLiXet al. Genetic characterization and pathogenicity of a novel recombined porcine reproductive and respiratory syndrome virus 2 among Nadc30-like, Jxa1-like, and Mlv-like strains. Viruses. (2018) 10:551. 10.3390/v10100551

  • 41.

    BianTSunYHaoMZhouLGeXGuoXet al. A recombinant type 2 porcine reproductive and respiratory syndrome virus between NADC30-like and a MLV-like: genetic characterization and pathogenicity for piglets. Infect Genet Evol. (2017) 54:279–86. 10.1016/j.meegid.2017.07.016

  • 42.

    ZhouLKangRZhangYYuJXieBChenCet al. Emergence of two novel recombinant porcine reproductive and respiratory syndrome viruses 2 (lineage 3) in Southwestern China. Vet Microbiol. (2019) 232:30–41. 10.1016/j.vetmic.2019.01.026

  • 43.

    HuangYWMengXJ. Novel strategies and approaches to develop the next generation of vaccines against porcine reproductive and respiratory syndrome virus (PRRSV). Virus Res. (2010) 154:141–9. 10.1016/j.virusres.2010.07.020

  • 44.

    YuFYanYShiMLiuHZZhangHLYangYBet al. Phylogenetics, genomic recombination, and NSP2 polymorphic patterns of porcine reproductive and respiratory syndrome virus in China and the United States in 2014–2018. J Virol. (2020) 94:e01813–19. 10.1128/JVI.01813-19

Summary

Keywords

PRRSV, phylogenetic analysis, nsp2, deletion, recombination

Citation

Wang J, Lin S, Quan D, Wang H, Huang J, Wang Y, Ren T, Ouyang K, Chen Y, Huang W, Luo T and Wei Z (2020) Full Genomic Analysis of New Variants of Porcine Reproductive and Respiratory Syndrome Virus Revealed Multiple Recombination Events Between Different Lineages and Sublineages. Front. Vet. Sci. 7:603. doi: 10.3389/fvets.2020.00603

Received

12 February 2020

Accepted

27 July 2020

Published

10 September 2020

Volume

7 - 2020

Edited by

Francisco Ruiz-Fons, Consejo Superior de Investigaciones Científicas (CSIC), Spain

Reviewed by

Ping Jiang, Nanjing Agricultural University, China; Yifeng Jiang, Chinese Academy of Agricultural Sciences (CAAS), China

Updates

Copyright

*Correspondence: Zuzhang Wei Tingrong Luo

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

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics