Evolutionary Analysis of Four Recombinant Viruses of the Porcine Reproductive and Respiratory Syndrome Virus From a Pig Farm in China

The porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most important pathogens causing substantial economic losses to the Chinese swine industry. In this study, we analyzed the complete genome sequences of four PRRSV isolates (PRRSV2/CN/SS0/2020, PRRSV2/CN/SS1/2021, PRRSV2/CN/L3/2021, and PRRSV2/CN/L4/2020) isolated from a single pig farm from 2020 to 2021. The genomes of the four isolates were 14,962–15,023 nt long, excluding the poly (A) tails. Comparative analysis of the genome sequences showed that the four isolates shared 93.2–98.1% homology and they had no close PRRSV relatives registered in the GenBank (<92%). Furthermore, PRRSV2/CN/SS0/2020 and PRRSV2/CN/SS1/2021 had characteristic 150-aa deletions (aa481+aa537-566 +aa628–747) that were identical to the live attenuated virus vaccine strain TJM-F92 (derived from the HP-PRRSV TJ). Further analysis of the full-length sequences suggests that the four isolates were natural recombinant strains between lineages 1 (NADC30-like), 3 (QYYZ-like), and 8.7 (JXA1-like). Animal experiments revealed discrepancies in virulence between PRRSV2/CN/SS0/2020 and PRRSV2/CN/L3/2021. The strain with high homology to HP-PRRSV demonstrates higher pathogenicity for pigs than the other isolate with low homology to HP-PRRSV. Taken together, our findings suggest that PRRSVs have undergone genome evolution by recombination among field strains/MLV-like strains of different lineages.

Introduction

The porcine reproductive and respiratory syndrome (PRRS) is caused by the porcine reproductive and respiratory syndrome virus (PRRSV). It was first reported in 1987 and has spread rapidly since to become a worldwide menace (1, 2). The PRRSV genome is an ~15 kb-long positive single-stranded RNA virus that has been classified under the Nidovirales order and Arteriviridae family (https://talk.ictvonline.org/taxonomy). The viral RNA contains at least ten open overlapping reading frames (ORFs): ORF1a, ORF1b, ORF2a, ORF2b, ORF3, ORF4, ORF5a, and ORF5–ORF7 (3–6). Among these ORFs, ORF1a and ORF1b comprise almost three-quarters of the viral genome and encode at least 16 nonstructural proteins (such as nsp1α, nsp1β, nsp2-6, nsp2TF, nsp2N, nsp7α, nsp7β, and nsp8–12), whereas other ORFs located at the 3' terminus code for eight structural proteins: GP2, E, GP3, GP4, GP5a, GP5, M, and N, respectively (3, 5, 7, 8).

All PRRSV strains isolated to date have been classified into two major genotypes, namely, the European (type 1) and North American (type 2) genotypes, with type 2 PRRSV being predominant in China since its initial report there in 1996 (9, 10). PRRS is considered one of the most important infectious diseases in Chinese swineherds as it causes severe economic losses in the swine industry every year. Currently, the Chinese type 2 PRRSV strains can be classified into different lineages, including lineage 1 (NADC30-like), 3 (QYYZ-like), 5 (VR2332-like), and 8 (JXA1-like and CH-1a-like) based on phylogenetic analysis of ORF5 sequences (11, 12). The diversity of the Chinese type 2 PRRSV has been increasing due to recombination events among the different PRRSV lineages since the emergence of NADC30-like PRRSVs in China in 2012 (12–23). In the present study, we identify four recombinant PRRSV strains circulating in a swine farm in Fujian province, China. Phylogenetic and molecular evolutionary analyses indicated that these strains evolved from natural recombination among the NADC30-, QYYZ-, and JXA1-like variants. To further understand the four strains, we genetically characterized the complete genomes of the viruses.

Materials and Methods

Sample Collection and Viral Isolation

From 2020 to 2021, a severe reproductive and respiratory disease was observed in pigs from four independent pens in a swine farm in Fujian province, China. This farm created a PRRS-positive stable herd by inoculating 0.05 dose/pig of MLV TJM-F92 vaccine (derived from the HP-PRRSV TJ) from 2018. The affected pigs from the two pens exhibited high fever (40.3–41.8°C), severe respiratory syndrome, and high morbidity (30%) and mortality (20%), whereas those from the other two pens exhibited fever, severe respiratory syndrome and reproductive problems, and especially high abortion rate (10%) was observed in the sows. PRRSV was detected using an RT-PCR kit (Beijing Anheal Laboratories Co., Ltd., Beijing, China) according to the manufacturer's instructions. PRRSV was positive when the CT value was ≤30. MARC-145^CD163 cell lines (a stable cell line highly expressing porcine CD163) were used for PRRSV isolation from the positive samples (12). The tissues (lung, lymph nodes, and serum samples) were homogenized with DMEM containing antibiotics and antimycotics to obtain a 50% (w/v) suspension. After freeze-thawing thrice, the samples were centrifuged at 10,000 × g for 10 min and the supernatants were sieved through a 0.22-μm filter. MARC-145^CD163 cells were mixed with 100 μL of supernatant and incubated at 37 °C with 5% CO₂ for 3–5 days. PRRSVs were passaged three times in MARC-145^CD163 cells for subsequent analysis after being plaque-purified three times.

RNA Isolation and RT-PCR

RNA was extracted from the positive tissues (lungs, serum, and lymph nodes) and the virus using a Viral RNA Mini Kit (TIANGEN, Beijing, China) according to the manufacturer's instructions. cDNA was generated using the HiScript^® III 1st Strand cDNA Synthesis Kit (Vazyme Biotech Co., Ltd, Nanjing, China) following the supplier's instructions. The full-length viral genomes were amplified using six viral-specific primers (Supplementary Table 1). Subsequently, we confirmed the complete genome using other PCR primers as described previously (24, 25). The purified PCR products were cloned into the pEASY^®-Blunt Simple Cloning Vector (TransGen Biotech, Beijing, China). To determine a consensus sequence of each fragment, at least three recombinant clones were sent to the Ruibo Life Technologies Corporation (Beijing, China) for sequencing by using the Sanger approach in both directions for each fragment.

Complete Genomic Sequence and Recombination Analysis

Forty-four representative PRRSV strains available in GenBank, including the type 1 and type 2 strains, were used for the comparative sequence analysis in this study (Table 1). Multiple sequence alignments and genome analyses were performed using the MEGA 7.0 software and the DNAStar 7.0 package. Phylogenetic trees were constructed using the neighbor-joining method in MEGA 7.0 and bootstrap confidence values from 1,000 replicates.

TABLE 1

Table 1. Representative PRRSV strains used in this study.

Potential recombination within the whole genome sequences was determined using seven methods (RDP, BootScan, GENECONV, Chimera, Maxchi, SiScan, and 3Seq) via the recombination detection program 4.10 (RDP 4.10) (26). A recombination event was identified when at least five of the seven methods reported recombination signals in RDP 4.10 with the highest acceptable p of 0.05 (27). The possible recombination event was further confirmed using the SIMPLOT software (version 3.5.1) with a 200-bp window width and a 20-bp step size (28).

Challenge Experiment

To evaluate the pathogenicity of the recombinant strains, two strains PRRSV2/CN/SS0/2020 (98.3% homology with JXA1) and PRRSV2/CN/L3/2021 (86.3% homology with JXA1) were selected for animal experiments. Fifteen 4-week-old pigs confirmed to be free of PRRSV and PCV2 were randomly divided into 3 groups (5 pigs/group). The pigs in group 1 were intranasally administered with 2 mL of PRRSV2/CN/SS0/2020 containing 2 × 10⁵ TCID50 and those in group 2 were intranasally administered with 2 mL of PRRSV2/CN/SS0/2020 containing 2 × 10⁵ TCID50. The pigs in the negative control group were inoculated with 2 mL of Dulbecco's Modified Eagle Medium (DMEM). Rectal temperature was recorded daily from 0 to 14 days post-challenge (dpc). Serum was collected at 0, 4, 7, 11, and 14 dpc to detect PRRSV-specific antibodies using PRRSV using ELISA kit (IDEXX Laboratories Inc., Westbrook, ME, USA). Viral load in the sera of each group was detected by an IFA-microtitration infectivity assay as previously described (29). All animals were euthanized at 14 dpc for necropsy. In addition, lung tissues of pigs were collected at necropsy and fixed in 10% neutral-buffered formalin and routinely processed for histological examination. The procedures of animal handling and experimentation performed in this study were approved by the Longyan University Animal Ethics Committee (Permit number: permit no.LY20210010X).

Data Analysis

Statistical analysis in this study was performed using a one- or two-way ANOVA analysis of variance in GraphPad Prism software (version 6.0), and the results were considered statistically significant when p < 0.05.

Results

Complete Genome Sequence Analysis of PRRSV

Four strains of PRRSV were isolated using MARC-145^CD163 cells and designated as PRRSV2/CN/SS0/2020 (Accession number: ON365556), PRRSV2/CN/SS1/2021 (Accession number: ON093974), PRRSV2/CN/L3/2021 (Accession number: OL416130), and PRRSV2/CN/L4/2020 (Accession number: OL422822). Typical PRRSV CPE characterized by cell fusion and shedding was observed in MARC-145CD¹⁶³ cells (Supplementary Figure 1). The genomes of the four isolates were 14,962–15,023 nt in length, excluding the poly (A) tails at the 3' end. Comparative analyses of genome sequences showed that the four isolates shared 93.2–98.1% homology and had a 83.7–86.8% identity with VR2332-like PRRSVs (VR-2332, RespPRRSV MLV, and BJ-4), 86.1–91.5% with JXA1-like PRRSVs (JXA1, HuN4, and TJd14-1), 81.9–86.6% with NADC30-like PRRSVs (NADC30, FJZ03, and Chsx1401), 82.9–86.5% with QYYZ-like PRRSVs (QYYZ, GM2, and FJFS), and only 59.2–59.8% with LV (Supplementary Table 2).

To evaluate the genomic characteristics of the four PRRSV isolates, each region of the genomes of the four strains was further compared with four viruses from different lineages, including NADC30-like strains, QYYZ-like strains, VR2332-like strains (VR-2332 and BJ-4), JXA1-like strains, and the LV (the prototype of type 1 PRRSV) strain (Table 2). The 5'-UTR of the four PRRSV isolates was 188–191 nt long and they had an identity of 87.3–91.5% with NADC30-like strains, 93.7–95.7% with QYYZ-like strains, 89.4–91.0% with VR2332-like strains, 97.4–98.4% with JXA1-like strains, and 61.5–62.2% with LV. The 3'-UTR of the four PRRSVs was found to be ~145–154 nt long, excluding the poly (A) tail, and their sequence was 86.1–89.9% homologous with that of the NADC30-like strains, 89.4–92.1% with that of the QYYZ-like strains, 90.7–91.4% with that of the VR2332-like strains, 85.3–88.7% with that of the JXA1-like strains, and 74.3–75.2% with that of LV (Table 2).

TABLE 2

Table 2. Sequence distance of four isolates in this study with reference strains.

ORF1a and ORF1b encode the non-structural proteins (Nsp) of PRRSV. A nucleotide sequence comparison of ORF1 showed that in the four strains, ORF1b was relatively conserved as compared to ORF1a. Of all the Nsp sequences, the most variable ones were from Nsp1β and Nsp2 within ORF1a, which had 72.4–94.9% and 63.4–86.9% amino acid identity with the reference strains, respectively (Table 2).

Open reading frames 2a to 7 encode the PRRSV structural proteins. Nucleotide sequence comparison of this region showed that the four strains shared 95.2–99.1% nucleotide identity, 84.9–86.4% identity with NADC30-like strains, 89.1–90.7% identity with QYYZ-like strains, 87.6–88.8% identity with VR2332-like strains, and 86.6–87.5% identity with JXA1-like strains. In contrast, these ORF regions only shared 65.2–65.3% identity with LV.

Phylogenetic Analysis of PRRSV

Four respective phylogenetic trees were constructed based on the full-length genome sequences, Nsp2, ORF2-7, and ORF5 nucleotide sequences of the four PRRSV isolates and 27 representative strains of type 2 PRRSV.

Interestingly, the tree based on the full-length genome sequences showed that the four PRRSVs formed a separate branch that was in the middle of lineages 3 and 8.7, which are represented by QYYZ and JXA1, respectively (Figure 1A). The tree based on the ORF2-7 or ORF5 sequences indicates that the four strains were grouped in lineage 3, together with the QYYZ-like strains (Figures 1C,D). The tree constructed based on Nsp2 indicates that PRRSV2/CN/L3/2021 and PRRSV2/CN/L4/2020 were clustered into lineage 1, together with the NADC30-like strains, and PRRSV2/CN/SS0/2020 was clustered into lineage 8.7, together with the JXA1-like strains; however, PRRSV2/CN/SS1/2021 was in a separate branch (Figure 1B).

FIGURE 1

Figure 1. Phylogenetic tree based on the full length (A), Nsp2 (B), ORF2-7 (C), and ORF5 (D) of the four PRRSV strains from this study and reference strains. The four PRRSV strains isolated in this study are labeled with a red circle. The reliability of the tree was evaluated by bootstrapping using Mega 7.0 with 1000 replicates.

Amino Acid Analysis of NSP2 and ORF5

The Nsp2-coding region is recognized as one of the most variable proteins with different deletions and insertions. Therefore, Nsp2 is often used to analyze the genetic variations and molecular epidemiology of PRRSV. As shown in Figure 2, two isolates, namely PRRSV2/CN/L3/2021 and PRRSV2/CN/L4/2020, had a unique discontinuous deletion of 131-aa in the Nsp2-coding region, which was identical to that of the NADC30 and NADC30-like strains. Meanwhile, PRRSV2/CN/L4/2020 had an additional 1-aa deletion at the 15th aa. Interestingly, both PRRSV2/CN/L3/2021 and PRRSV2/CN/L4/2020 had two additional 2-aa insertions at amino acid positions 224–225 (Figure 2A). Notably, PRRSV2/CN/SS0/2020 and PRRSV2/CN/SS1/2021 had characteristic 150-aa deletions (aa481+aa537-566+aa628–747) that were identical to those of the live attenuated virus vaccine strain TJM-F92 (derived from the HP-PRRSV TJ) (30) and TJbd14–1 (13), which is an MLV-like strain that evolved from the HP-PRRSV vaccine virus TJM-F92 (Figure 2B). Notably, comparative sequence analysis showed that PRRSV2/CN/SS0/2020 and PRRSV2/CN/SS1/2021 exhibited the highest nucleotide similarity (89.5–95.3%) and amino acid similarity (87.0–94.0%) with TJbd14-1 and TJM-F92.

FIGURE 2

Figure 2. Amino acids alignment for Nsp2 of the four PRRSV isolates with representative strains. (A) The 131-aa discontinuous deletions (aa324–434, aa485, and aa504–522) highlighted in gray boxes show the deletion signature of PRRSV2/CN/L3/2021, PRRSV2/CN/L4/2021, and NADC30-like strains. Additional 2-aa insertion in PRRSV2/CN/L3/2021 and PRRSV2/CN/L4/2021 is marked in the green box. (B) The 150-aa discontinuous deletions (aa481+aa537-566+aa628–747) highlighted in gray boxes show the deletion signature of PRRSV2/CN/SS0/2020, PRRSV2/CN/SS1/2021, and MLV TJM-F92 strain.

The GP5 protein encoded by ORF5 is the most variable structural protein of the virus. At least six antigenic regions (ARs) were reported within GP5 (AR1-15, AR27-35, AR37-51, AR149-156, AR166-181, and AR192-200) (31, 32). Three ARs (AR27-35, AR37-51, and AR192-200) at the N-terminus of the four strains were the most similar to the QYYZ-like strains but differed from the JXA1-like, VR-2332-like, and NADC30-like strains. Additionally, seven unique amino acids distributed in GP5 were only identified in QYYZ-like strains and the four strains (I²⁶, Y³⁸, C⁶⁶, S⁹², F¹¹⁷, I¹⁵², and H¹⁹⁹) (Figure 3). Studies have shown that the 13th and 151st amino acids of GP5 are related to the virulence of the virus (33). In the present study, four strains had both R¹³Q and R¹⁵¹G mutation as compared with the VR2332 strain.

FIGURE 3

Figure 3. Analysis and comparison of amino acid mutations in GP5 compared to reference viruses. The signal peptide and transmembrane (TM) domains in GP5 are highlighted by rectangles. Important amino acid changes between four strains (red cicrlce) and QYYZ-like PRRSVs in GP5 are indicated in red.

Recombination Analysis

The recombination analysis using the RDP version 4.10 software revealed that the four isolates are the result of recombination between JXA1-like, QYYZ-like, and NADC30-like strains circulating in China (Table 3, Figure 4, Supplementary Figure 2). Additionally, the putative recombinant events and statistically incongruent phylogenetic trees were further confirmed using SimPlot v3.5.1 and statistically incongruent phylogenetic trees (Figure 4, Supplementary Figure 2). From the similarity plot, the breakpoints of recombination events of four strains are mainly located in nsp1, nsp2, nsp3, nsp11, and nsp12 (Figure 4, Supplementary Figure 2).

TABLE 3

Table 3. Information of recombination events of four isolates.

FIGURE 4

Figure 4. Similarity plot analysis (A) and phylogenetic trees analysis based on different regions (B) of PRRSV2/CN/SS0/2020. For similarity plot analysis, the y-axis indicates the percentage identity between the parental sequences and the query sequence. The complete genome RRSV2/CN/SS0/2020 was chosen as the query sequence, respectively.

PRRSV2/CN/SS0/2020 and PRRSV2/CN/L3/2021 Exhibited Pathogenicity for Pigs

Pigs infected with PRRSV2/CN/SS0/2020 exhibited obvious clinical signs, including respiratory distress, anorexia, and coughing from 3 dpc. The body temperatures of PRRSV2/CN/SS0/2020-inoculated group began to reach above 40°C at 3–14 dpc, with a peak (>41 °C) at 8–9 dpc. In contrast, the pigs infected with PRRSV2/CN/L3/2021 had less fever and less severe clinical signs (Figure 5A). Control pigs showed normal rectal temperature and behavior throughout the experiment. One of five PRRSV2/CN/SS0/2020-infected pigs died at 10 dpc, whereas no mortality was observed in PRRSV2/CN/L3/2021-infected group and negative control group. In addition, pigs in PRRSV2/CN/SS0/2020-infected group had significantly more severe interstitial pneumonia than did pigs inoculated with PRRSV2/CN/L3/2021 (Figure 5B), which indicates that PRRSV2/CN/SS0/2020 is a highly virulent PRRSV isolate for pigs.

FIGURE 5

Figure 5. The rectal temperature, microscopic lung lesions, viremia, and antibody levels of pigs during the challenge experiment. (A) Rectal temperatures of pigs inoculated with PRRSV2/CN/SS0/2020, PRRSV2/CN/L3/2021, and DMEM. The clinical fever cut-off value was set at 40 °C. (B) Microscopic lung lesions of the inoculated pigs. (C) Dynamics of viral loads of the inoculated pigs during the experiment. (D) PRRSV-specific antibodies response in serum of challenged pigs. The cut-off value for seroconversion was set at a sample-to-positive (s/p) ratio of 0.4. Data are expressed as mean ± standard deviation (S.D). Asterisk indicates significant differences between the PRRSV2/CN/SS0/2020- and PRRSV2/CN/L3/2021-infected groups (*p < 0.05).

The PRRSV was isolated from the serum of pigs challenged with the two strains at 4 dpc. The mean virus titers of PRRSV2/CN/SS0/2020 or PRRSV2/CN/L3/2021 reached peak levels at 7 dpc (10^6.2 TCID50/mL and 10^5.0 TCID50/mL, respectively) and high viremia persisted until 14 dpc in every PRRSV-infected group (Figure 5C). Furthermore, the mean virus loads for PRRSV2/CN/SS0/2020-infected group were significantly higher than those of the PRRSV2/CN/L3/2021-infected group at 7 and 11 dpc (p < 0.05). Meanwhile, no PRRSV was detected in the control group throughout the experiment. The specific antibodies of PRRSV N protein in the sera of pigs were detected by a commercial ELISA kit. All pigs in the PRRSV2/CN/SS0/2020-infected group became seroconverted at 7 dpc, whereas PRRSV2/CN/L3/2021-infected were seroconverted at 11 dpc (Figure 5D). In contrast, PRRSV-specific antibodies were not detected in the control group throughout the experiment.

Discussion

Porcine reproductive and respiratory syndrome virus is one of the most important pathogens causing substantial economic losses to the Chinese swine industry. According to the global PRRSV classification system, four different lineages of PRRSV (lineages 1, 3, 5, and 8) co-existed in Chinese swine herds, and two of them (lineages 1 and 8) are currently the predominant strains circulating in China. With the emergence of NADC30-like PRRSVs in China since 2012, recombinant strains have been frequently observed in the field. Notably, multiple studies have shown that PRRSV variants evolved from recombination events between lineage 1 and one or two other lineages (12–23). In this study, four PRRSV strains (PRRSV2/CN/SS0/2020, PRRSV2/CN/SS1/2021, PRRSV2/CN/L3/2021, and PRRSV2/CN/L4/2020) were also identified as recombinant viruses from three lineages of type 2 PRRSV (NADC30-like, QYYZ-like, and JXA1-like viruses). Importantly, the four strains in this study had a low whole-genome nucleotide similarity with other PRRSV sequences in GenBank (<92%), which raised the concern of the introduction of new viruses into other large-scale pig farms. Therefore, to further our knowledge about these recombinant PRRSV strains, we characterized them genetically and phylogenetically to obtain insights that link genotypic with phenotypic data.

The natural genomic recombination between different lineages of PRRSV has been shown to play an important role in the generation of novel strains (34–36). Several studies have shown that the patterns of recombination of PRRSV in the field are becoming increasingly complex; for instance, 14LY01-FJ, 14LY02-FJ, 15LY01-FJ, 15LY02-FJ, JL580, GXNN1839, GXYL1403, HeN1401, HeN1601, FJ1402, FJXS15, SC-d, and TJnh1501 have evolved from the recombination between two lineages (17, 18, 20, 22, 36–39), whereas SD17–38, PRRSV2/CN/N9185/2018, PRRSV2/CN/X4833/2018, SDhz1512, and SCcd16 were natural recombinant viruses among three lineages (12, 14, 16, 23). More surprisingly, FJLIUY-2017 was derived from recombinant strains of four lineages (15). Notably, the diversity of recombinant PRRSV strains poses a major obstacle to the effective control of viral transmission (38, 40). In the present study, a comparison of whole-genome sequences and recombination analysis revealed that the four recombinant isolates (PRRSV2/CN/L3/2021, PRRSV2/CN/L4/2020, PRRSV2/CN/SS02020, and PRRSV2/CN/SS1/2021) were also natural recombinant viruses among the three lineages (lineage 1, 3, and 8.7). However, full-length sequencing of the virus showed that the entire genome of the four strains had low identity with the PRRSV strains available in GenBank, suggesting that PRRSVs have undergone evolution via natural recombination in recent years. Previous studies showed that Nsp1, Nsp2, Nsp3, Nsp9, Nsp11, and ORF2 are the hot spots for PRRSV RNA recombination (12, 13, 15, 17, 22, 23, 25, 36). In the present study, breakpoints of the four recombinant strains were also mainly located in Nsp1, Nsp2, Nsp3, and Nsp11, indicating that PRRSV gains genetic diversity by increasing recombination events at specific regions.

Nsp2 and ORF5 are the most variable proteins in the viral genome and are usually used as target genes for the molecular epidemiological surveillance of PRRSV (27). Compared with VR2332, the Nsp2 gene of PRRSV contains different patterns of amino acid insertions and deletions. For example, an outbreak in China with highly pathogenic PRRSVs (HP-PRRSV, JXA1-like) in 2006 had a unique 30-aa deletion in the Nsp2 region, whereas a unique 131-aa (111aa+1aa+19aa) deletion pattern in Nsp2 is the footprint of NADC30 and NADC30-like strains (11, 41). In the present study, the Nsp2 gene of the two strains (PRRSV2/CN/L3/2021 and PRRSV2/CN/L4/2020) had the 131-aa deletion, which is identical to NADC30-like strains, suggesting that these two strains may belong to lineage 1. Additionally, 2-aa insertions were also found in PRRSV2/CN/L3/2021 and PRRSV2/CN/L4/2020. Previous studies showed that MLV TJM-F92 has a unique 150-aa deletion (aa481+aa537-566+aa628–747) signature in its Nsp2-coding region (30). Interestingly, two strains identified in our study (PRRSV2/CN/SS0/2020 and PRRSV2/CN/SS1/2021) had characteristic 150-aa deletions that were identical to TJbd14-1 (an MLV-like strain that evolved from TJM-F92) and TJM-F92, suggesting that the PRRSV genome may have evolved more compactly by eliminating dispensable genomic regions (42). Furthermore, the phylogenetic analysis based on the ORF5 gene demonstrated that the four strains belonged to the QYYZ-like virus. In addition, the GP5 antigenic regions in the four strains were similar to those of the related lineage 3 representative strains but different from other lineage representative strains. The discordance between Nsp2 patterns and ORF5 lineages resulted from the recombination of PRRSVs. In other words, there is limited understanding of the recombinant PRRSV genomic data, especially if quantifying PRRSV diversity has only focused on ORF5 and Nsp2 analysis. In addition, the 137th aa of the GP5 protein is assumed to differentiate the attenuated vaccine strain (A¹³⁷) and the wild strain (S¹³⁷). Residue S¹³⁷ was present in all four strains, suggesting that the four strains may be wild viruses (43). Additionally, a retrospective survey found that a massive vaccination with HP-PRRSV live attenuated vaccine (TJM-F92) was effective in preventing and controlling PRRS in this farm since 2018. Thus, we propose that PRRSV2/CN/SS0/2020 and PRRSV2/CN/SS1/2021 may have evolved from recombination among MLV TJM-F92 and NADC30-like or other PRRSV strains. Whether this strain directly evolved from MLV TJM-F92 needs to be further investigated.

Viremia and severity of fever are closely related to PRRSV virulence (44). In the present study, PRRSV2/CN/SS0/2020-infected pigs had higher viremia than PRRSV2/CN/L3/2021 from 4 dpc to 14 dpc, suggesting that PRRSV2/CN/SS0/2020 has higher virulence for pigs. In addition, pigs inoculated with the PRRSV2/CN/SS0/2020 strain had persistently higher fever (>40°C for 12 days) and interstitial pneumonia compared to the PRRSV2/CN/L3/2021-infected group. Notably, one pig in the PRRSV2/CN/SS0/2020 group died within 2 weeks, whereas the pigs in the PRRSV2/CN/L3/2021 survived throughout the experiment. These data showed that PRRSV2/CN/SS0/2020 demonstrates higher pathogenicity than PRRSV2/CN/L3/2021. This finding is consistent with previous results that recombinant PRRSV strains derived from field PRRSV strains and vaccine-like strains could result in higher pathogenicity for pigs, such as FJXS15 and TJnh150 (13, 18).

In conclusion, we determined the complete genome sequences of four recombinant PRRSV isolates (PRRSV2/CN/SS0/2020, PRRSV2/CN/SS1/2021, PRRSV2/CN/L3/2021, and PRRSV2/ CN/L4/2020). Four isolates were natural recombinant among three lineages (lineage 1, 3, and 8) and have the absence of close relatives to PRRSV sequence in GenBank (<92%). PRRSV2/CN/SS0/2020 exhibits higher pathogenicity than PRRSV2/CN/L3/2021. Our findings suggest that PRRSVs have undergone genome evolution by recombination among field strains/MLV-like strains of different lineages.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

JLiu and CW performed the experiments. YX, CL, and YY conceived and designed the experiments. LL, CH, and JLi performed the data analyses. JLiu completed the writing of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by Leading Project Foundation of Science Department of Fujian Province (2021N0032), Research Foundation of Engineering Research Center for the Prevention and Control of Animal Original Zoonosis, Fujian Province University, (Grant No. 2021ZN002), Qimai Foundation of Shanghang District of Longyan City, Fujian Province (2020SHQM07) and Longyan University Young Talents Program.

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.

Publisher's Note

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.933896/full#supplementary-material

References

1. Benfield DA, Nelson E, Collins JE, Harris L, Goyal SM, Robison D, et al. Characterization of Swine Infertility and Respiratory Syndrome (SIRS) virus (isolate ATCC VR-2332). J Vet Diagn Invest. (1992) 4:127–33. doi: 10.1177/104063879200400202

CrossRef Full Text | Google Scholar

2. Wensvoort G, Terpstra D, Pol JM, ter Laak EA, Bloemraad M, de Kluyver EP, et al. Mystery swine disease in the Netherlands: the isolation of Lelystad virus. Vet Q. (1991) 13:121–30. doi: 10.1080/01652176.1991.9694296

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Firth AE, Zevenhoven-Dobbe JC, Wills NM, Go YY, Balasuriya UB, Atkins JF, et al. Discovery of a small arterivirus gene that overlaps the GP5 coding sequence and is important for virus production. J Gen Virol. (2011) 92:1097–106. doi: 10.1099/vir.0.029264-0

PubMed Abstract | CrossRef Full Text | Google Scholar

4. 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:1107–16. doi: 10.1099/vir.0.030213-0

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lunney JK, Fang Y, Ladinig A, Chen N, Li Y, Rowland B, et al. Porcine reproductive and respiratory syndrome virus (PRRSV): pathogenesis and interaction with the immune system. Annu Rev Anim Biosci. (2016) 4:129–54. doi: 10.1146/annurev-animal-022114-111025

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Muhlenberg JJ. PRRSV, the virus. Vet Res. (2000) 31:11–21. doi: 10.1051/vetres:2000103

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Dokland T. The structural biology of PRRSV. Virus Res. (2010) 154:86–97. doi: 10.1016/j.virusres.2010.07.029

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wu WH, Fang Y, Rowland RR, Lawson SR, Christopher-Hennings J, Yoon KJ, et al. The 2b protein as a minor structural component of PRRSV. Virus Res. (2005) 114:177–81. doi: 10.1016/j.virusres.2005.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Nelsen CJ, Murtaugh MP, Faaberg KS. Porcine reproductive and respiratory syndrome virus comparison: divergent evolution on two continents. J Virol. (1999) 73:270–80. doi: 10.1128/JVI.73.1.270-280.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zhou YJ, Hao XF, Tian ZJ, Tong GZ, Yoo D, An TQ, et al. Highly virulent porcine reproductive and respiratory syndrome virus emerged in China. Transbound Emerg Dis. (2008) 55:152–64. doi: 10.1111/j.1865-1682.2008.01020.x

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Guo Z, Chen XX, Li R, Qiao S, Zhang G. The prevalent status and genetic diversity of porcine reproductive and respiratory syndrome virus in China: a molecular epidemiological perspective. Virol J. (2018) 15:2. doi: 10.1186/s12985-017-0910-6

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Liu J, Xu Y, Lin Z, Fan J, Dai A, Deng X, et al. Epidemiology investigation of PRRSV discharged by faecal and genetic variation of ORF5. Transbound Emerg Dis. (2021) 68:2334–44. doi: 10.1111/tbed.13894

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Bian T, Sun Y, Hao M, Zhou L, Ge X, Guo X, et 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. doi: 10.1016/j.meegid.2017.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chen N, Ye M, Li S, Huang Y, Zhou R, Yu X, et 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. doi: 10.1111/tbed.12952

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Liu J, Wei C, Lin Z, Fan J, Xia W, Dai A, et 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. doi: 10.1016/j.meegid.2018.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Liu Y, Li J, Yang J, Zeng H, Guo L, Ren S, et al. Emergence of different recombinant porcine reproductive and respiratory syndrome viruses, China. Sci Rep. (2018) 8:4118. doi: 10.1038/s41598-018-22494-4

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Liu JK, Zhou X, Zhai JQ, Li B, Wei CH, Dai AL, et al. Emergence of a novel highly pathogenic porcine reproductive and respiratory syndrome virus in China. Transbound Emerg Dis. (2017) 64:2059–74. doi: 10.1111/tbed.12617

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Liu J, Zhou X, Zhai J, Wei C, Dai A, Yang X, et al. Recombination in JXA1-R vaccine and NADC30-like strain of porcine reproductive and respiratory syndrome viruses. Vet Microbiol. (2017) 204:110–20. doi: 10.1016/j.vetmic.2017.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sun YK, Li Q, Yu ZQ, Han XL, Wei YF, Ji CH, et al. Emergence of novel recombination lineage 3 of porcine reproductive and respiratory syndrome viruses in Southern China. Transbound Emerg Dis. (2019) 66:578–87. doi: 10.1111/tbed.13067

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wang HM, Liu YG, Tang YD, Liu TX, Zheng LL, Wang TY, et al. A natural recombinant PRRSV between HP-PRRSV JXA1-like and NADC30-like strains. Transbound Emerg Dis. (2018) 65:1078–86. doi: 10.1111/tbed.12852

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zhang Q, Xu X, You S, Li Y, Wang H, Bai J, et al. Emerging of two new subgenotypes of porcine reproductive and respiratory syndrome viruses in Southeast China. Microb Pathog. (2016) 97:27–33. doi: 10.1016/j.micpath.2016.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zhao K, Ye C, Chang XB, Jiang CG, Wang SJ, Cai XH, et 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. doi: 10.1128/JVI.01446-15

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhou L, Kang R, Zhang Y, Ding M, Xie B, Tian Y, et 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. doi: 10.3390/v10060328

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Zhang H, Leng C, Ding Y, Zhai H, Li Z, Xiang L, et al. Characterization of newly emerged NADC30-like strains of porcine reproductive and respiratory syndrome virus in China. Arch Virol. (2019) 164:401–11. doi: 10.1007/s00705-018-4080-7

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhou L, Kang R, Ji G, Tian Y, Ge M, Xie B, et 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. doi: 10.1007/s11262-017-1519-y

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. (2010) 26:2462–3. doi: 10.1093/bioinformatics/btq467

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Yin B, Qi S, Sha W, Qin H, Liu L, Yun J, et al. Molecular characterization of the Nsp2 and ORF5 (ORF5a) Genes of PRRSV strains in nine provinces of China during 2016–2018. Front Vet Sci. (2021) 8:605832. doi: 10.3389/fvets.2021.605832

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lole KS, Bollinger RC, Paranjape RS, Gafkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. (1999) 73:152–60. doi: 10.1128/JVI.73.1.152-160.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zhou L, Yang B, Xu L, Jin H, Ge X, Guo X, et al. Efficacy evaluation of three modified-live virus vaccines against a strain of porcine reproductive and respiratory syndrome virus NADC30-like. Vet Microbiol. (2017) 207:108–16. doi: 10.1016/j.vetmic.2017.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Leng X, Li Z, Xia M, Li X, Wang F, Wang W, et 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. doi: 10.1016/j.vetmic.2011.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

31. de Lima M, Pattnaik AK, Flores EF, Osorio FA. Serologic marker candidates identified among B-cell linear epitopes of nsp2 and structural proteins of a North American strain of porcine reproductive and respiratory syndrome virus. Virology. (2006) 353:410–21. doi: 10.1016/j.virol.2006.05.036

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhou YJ, Yua H, Tian ZJ, Liu JX, An TQ, Peng JM, et al. Monoclonal antibodies and conserved antigenic epitopes in the C terminus of GP5 protein of the North American type porcine reproductive and respiratory syndrome virus. Vet Microbiol. (2009) 138:1–10. doi: 10.1016/j.vetmic.2009.01.041

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Allende R, Kutish GF, Laegreid W, Lu Z, Lewis TL, Rock DL, et al. Mutations in the genome of porcine reproductive and respiratory syndrome virus responsible for the attenuation phenotype. Arch Virol. (2000) 145:1149–61. doi: 10.1007/s007050070115

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kappes MA, Faaberg KS. PRRSV structure, replication and recombination: origin of phenotype and genotype diversity. Virology. (2015) 479:475–86. doi: 10.1016/j.virol.2015.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Shi M, Holmes EC, Brar MS, Leung FC. Recombination is associated with an outbreak of novel highly pathogenic porcine reproductive and respiratory syndrome viruses in China. J Virol. (2013) 87:10904–7. doi: 10.1128/JVI.01270-13

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zhang Z, Zhou L, Ge X, Guo X, Han J, Yang H. 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. doi: 10.1016/j.meegid.2018.09.024

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhang Q, Jiang P, Song Z, Lv L, Li L, Bai J. 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. doi: 10.1016/j.vetmic.2016.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wang J, Lin S, Quan D, Wang H, Huang J, Wang Y, et al. 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. (2020) 7:603. doi: 10.3389/fvets.2020.00603

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Sun YF, Zhou L, Bian T, Tian XX, Ren WK, Lu C, et al. Efficacy evaluation of two commercial modified-live virus vaccines against a novel recombinant type 2 porcine reproductive and respiratory syndrome virus. Vet Microbiol. (2018) 216:176–82. doi: 10.1016/j.vetmic.2018.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Meng XJ. Heterogeneity of porcine reproductive and respiratory syndrome virus: implications for current vaccine efficacy and future vaccine development. Vet Microbiol. (2000) 74:309–29. doi: 10.1016/S0378-1135(00)00196-6

PubMed Abstract | CrossRef Full Text | Google Scholar

41. 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:2256–7. doi: 10.3201/eid2112.150360

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Han J, Wang Y, Faaberg KS. Complete genome analysis of RFLP 184 isolates of porcine reproductive and respiratory syndrome virus. Virus Res. (2006) 122:175–82. doi: 10.1016/j.virusres.2006.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Wesley RD, Mengeling WL, Lager KM, Clouser DF, Landgraf JG, Frey ML. Differentiation of a porcine reproductive and respiratory syndrome virus vaccine strain from North American field strains by restriction fragment length polymorphism analysis of ORF5. J Vet Diagn Invest. (1998) 10:140–4. doi: 10.1177/104063879801000204

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Stadejek T, Larsen LE, Podgórska K, Bøtner A, Botti S, Dolka I, et al. Pathogenicity of three genetically diverse strains of PRRSV Type 1 in specific pathogen free pigs. Vet Microbiol. (2017) 209:13–9. doi: 10.1016/j.vetmic.2017.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar