Genomic characteristics of a novel emerging PRRSV branch in sublineage 8.7 in China

Porcine reproductive and respiratory syndrome virus (PRRSV) has caused serious economic losses to the pig industry worldwide. During the continuous monitoring of PRRSV, a new PRRSV strain type with novel characteristics was first identified in three different regions of Shandong Province. These strains presented a novel deletion pattern (1 + 8 + 1) in the NSP2 region and belonged to a new branch in sublineage 8.7 based on the ORF5 gene phylogenetic tree. To further study the genomic characteristics of the new-branch PRRSV, we selected a sample from each of the three farms for whole-genome sequencing and sequence analysis. Based on the phylogenetic analysis of the whole genome, these strains formed a new independent branch in sublineage 8.7, which showed a close relationship with HP-PRRSV and intermediate PRRSV according to nucleotide and amino acid homology but displayed a completely different deletion pattern in NSP2. Recombinant analysis showed that these strains presented similar recombination patterns, all of which involved recombination with QYYZ in the ORF3 region. Furthermore, we found that the new-branch PRRSV retained highly consistent nucleotides at positions 117–120 (AGTA) of a quite conserved motif in the 3’-UTR; showed similar deletion patterns in the 5’-UTR, 3’-UTR and NSP2; retained characteristics consistent with intermediate PRRSV and exhibited a gradual evolution trend. The above results showed that the new-branch PRRSV strains may have the same origin and be similar to HP-PPRSV also evolved from intermediate PRRSV, but are distinct strains that evolved simultaneously with HP-PRRSV. They persist in some parts of China through rapid evolution, recombine with other strains and have the potential to become epidemic strains. The monitoring and biological characteristics of these strains should be further studied.

Introduction

Porcine reproductive and respiratory syndrome (PRRS) is one of the most economically important infectious diseases of swine globally (Gao et al., 2017). The etiologic agent of PRRS has been identified as porcine reproductive and respiratory syndrome virus (PRRSV), which is currently considered to consist of two different species in the Arteriviridae family: Betaarterivirus suid 1 (PRRSV-1) and Betaarterivirus suid 2 (PRRSV-2) (Brinton et al., 2021). The disease causes reproductive failure in sows, mainly premature parturition, late-term abortions and farrowing of stillborn and nonviable piglets. Respiratory dysfunction can develop predominantly in younger pigs (Albina, 1997). Currently, the most widely prevalent strain is PRRSV-2 in North America and Asia, which is classified into nine lineages (1–9) based on phylogenetic analyses of ORF5 gene sequences collected worldwide (Shi et al., 2010), among which five main sublineages (1.5 and 1.8, 3.5, 5.1, 8.7) occur in China (Guo et al., 2018; Zhang et al., 2018; Xu et al., 2020; Yu et al., 2020).

Sublineage 8.7 was the earliest strain found in mainland China (Liu et al., 2017; Guo et al., 2018) and mainly includes three subgroups: Classical PRRSV [representing strain CH-1a (Guo et al., 1996)], intermediate PRRSV [representing strains HB-1(sh)/2002 and HB-2(sh)/2002 (Gao et al., 2004)] and highly pathogenic PRRSV [HP-PRRSV; representing strains JXA1 (Tian et al., 2007) and HuN4 (Tong et al., 2007)]. Research by An showed that HP-PRRSV appeared in China in 2006 originated from the Classical CH-1a-like PRRSV. Its evolution was based on the gradual accumulation of mutations in intermediate PRRSV strains (Zhou et al., 2009a; An et al., 2010). The pathogenicity of the viruses shows a gradually increasing trend from Classical PRRSV to intermediate PRRSV to HP-PRRSV (Zhou and Yang, 2010). Classical PRRSV and intermediate PRRSV show moderate pathogenicity (Yang et al., 2008; Li et al., 2009; Ma et al., 2009). HP-PRRSV has the strongest pathogenicity and has been one of the main circulating strains in Chinese pig herds since 2006, causing serious economic losses to the pig industry in China and Southeast Asian countries (Tian et al., 2007; Tong et al., 2007). Sublineage 1.5 and 1.8 PRRSVs have gradually replaced sublineage 8.7 PRRSV as major epidemic strains in China (Xu et al., 2022), but sublineage 8.7 strains have not disappeared, and their genetic evolution has become increasingly complex (Jiang et al., 2020). We have been paying attention to whether this branch will evolve into another new branch similar to HP-PRRSV.

In this study, a new PRRSV strain type with novel characteristics was first found in three different regions of Shandong Province during the continuous monitoring of PRRSV. These strains belonged to a new independent branch in sublineage 8.7 and showed a novel deletion pattern in the NSP2 region. Further analysis of their whole-genome characteristics alongside all strains on this branch revealed that the new branch strains presented similar characteristics to intermediate PRRSV. Furthermore, multiple unique deletions and point mutations were first found in the PRRSV genome for the first time herein, which providing an important data reference for the epidemiology of PRRSV in China.

Materials and methods

Sample collection

In 2020–2022, 31 clinical samples of pigs suspected of PRRSV infection, including lung, lymph node and serum samples, were collected from farms in three different regions (Qingdao, Weihai and Yantai) of Shandong Province in northern China. The morbidity caused by suspected PRRSV infection on these pig farms was approximately 10–20%, and most of the affected pigs had mixed bacterial infections. These samples were transported at low temperature and stored at −20°C. Sample names, collection times, types, isolation regions and accession numbers are summarized in Table 1.

TABLE 1

Table 1. Information of new-branch PRRSV.

RNA extraction and RT–PCR

Tissue samples were ground with 1 ml PBS containing 2% penicillin streptomycin. Total RNA extraction from lapping fluid, cDNA preparation, RT-PCR and viral genome sequencing were performed as described previously (Xu et al., 2022). The primers used to detect PRRSV have been previously reported (Zhang et al., 2019). Seven pairs of overlapping primers were designed based on sublineage 8.7 and used to amplify the complete genome (Table 2).

TABLE 2

Table 2. Primers for complete genome amplification.

Sequence alignment and phylogenetic analysis

A total of 2,127 ORF5 sequences, 26 NSP2 sequences and 26 complete genomic sequences were obtained from GenBank as reference strain sequences. The nucleotide (nt) and amino acid (aa) similarity comparison of each gene fragment and the deduced amino acid sequences of the full-length genomic sequence alignment with reference strain sequences were assessed using the ClustalW method in Lasergene software (DNASTAR Inc., Madison, WI, United States) (Xu et al., 2022). Multiple sequence alignments of ORF5, NSP2, ORF3 and the whole genome were conducted using ClustalW in MEGA version 7.0 (Song et al., 2020), and phylogenetic trees were also constructed using the maximum likelihood method with 1,000 replicates and the Kimura 2-parameter substitution model (Kumar et al., 2016).

Recombination analysis

To test for recombination, SimPlot (version 3.5.1) was applied with a 200-bp window, sliding along the genome alignments with a 20-bp step size. We also further checked for potential recombination events by using seven different algorithms (RDP, GENECONV, Boot Scan, Max Chi, chimera, Si Scan, and 4Seq) in RDP4 software (Li et al., 2022).

Results and discussion

Clinical sample detection

Thirty-one clinical samples from animals with suspected PRRSV infection were collected from farms in three different regions of Shandong Province in China between 2020 and 2022. Among these samples, 18 (58.1%) samples were positive for PRRSV by RT–PCR. Eighteen ORF5 sequences, 6 NSP2 gene sequences and 3 complete genomic sequences of these samples were amplified and sequenced for further analysis with reference strain sequences. Strain names, areas and accession numbers are summarized in Table 1.

Phylogenetic analysis of ORF5 gene

In the construction phylogenetic trees of PRRSV, the ORF5 sequence is often used as the ideal target (Shi et al., 2010). Sequencing of GP5 revealed that these strains belonged to sublineage 8.7 PRRSV (data not shown). To understand the evolutionary relationships of these newly obtained PRRSV strains with representative strains in China, phylogenetic trees based on the sublineage 8.7 PRRSV ORF5 gene sequence set (2127) collected from GenBank was constructed using the maximum likelihood method in MEGA version 7.0. The results demonstrated that all 18 new PRRSV strains in this study formed a new independent branch within sublineage 8.7 and belonged to transitional PRRSV strains between intermediate PRRSV and HP-PRRSV (Figure 1A). To visualize the results, all strains of the new branch are displayed in the ORF5 phylogenetic tree in Figure 1B. It is worth noting that some sequences from GenBank were also assigned to this new branch, and only the ZJXS1412, HB2104 and JS18-3 strains had complete genome information. To explore the origin and evolution of the new-branch strains, we also included these strains from GenBank in our subsequent analysis. Information of these strains is shown in Table 1.

FIGURE 1

Figure 1. Phylogenetic analysis of GP5 and NSP2 genes of new-branch PRRSV. (A) Phylogenetic tree constructed based on the GenBank dataset (n = 2,127) of the sublineage 8.7 PRRSV ORF5 gene from China; (B) Phylogenetic tree constructed based on the ORF5 gene of new-branch PRRSV isolates and reference strains from different sublineages; (C) Phylogenetic tree constructed based on the NSP2 gene of new-branch PRRSV isolates and reference strains from different sublineages. The new-branch PRRSV strains obtained in this study are labeled with red triangles (▲). The strains obtained from GenBank belonging to the new branch are marked with blue squares (■).

Phylogenetic and sequence alignment analyses of NSP2 gene

In recent years, an increasing number of studies have reported that the recombination events continue to increase in PRRSV in China (Yu et al., 2020; Chen et al., 2021). Therefore, the PRRSV classification criteria only based on the ORF5 gene alone may not be rigorous. The Nsp2 gene is one of the most variable regions in the entire PRRSV genome (Nelsen et al., 1999), which can tolerate some amino acid deletion, insertion or mutation and it is commonly usually used as a basis for classifying different strains. According to our phylogenetic analysis based on the NSP2 gene, ZJXS1412 from GenBank and 6 new PRRSV strains (SDWH86, SDQD95, SDYT91, SDQD94, SDHW3, and SDHW8) were classified into sublineage 8.7 and formed a relatively independent branch parallel to the HP-PRRSV branch, while HB2104 and JS18-3 from GenBank belonged to sublineage 1.8 (Figure 1C). However, amino acid alignment based on the Nsp2 gene results showed that the new-branch strains displayed a completely different deletion pattern in the Nsp2 region from other sublineage 8.7 strains, retained characteristics consistent with intermediate PRRSV, and presented a gradual evolutionary trend. Most of the new-branch PRRSV, intermediate PRRSV and HP-PRRSV strains showed a deletion at the 481th aa of NSP2 relative to the reference strain VR-2332 (Figure 2). JS18-3 and HB2104 showed the same deletion pattern with the NADC30 strain (111 aa+1 aa+19 aa) due to recombination with the NADC30 strain (Han et al., 2020). ZJXS1412 maintained the same deletions at positions 481 and 533–561 (1 aa+29 aa) as HP-PRRSV. SDWH86, SDQD95 and SDYT91 showed further deletions at positions 15, 482–488 and 583 on the basis of the 481st aa deletion. To the best of our knowledge, this is a novel deletion pattern (1 aa+8 aa+1 aa) that has not been previously reported in sublineage 8.7 PRRSV in China. Although studies have shown that the 30 aa deletion in Nsp2 of the HP-PRRSV virus observed in China is not related to its virulence (Zhou et al., 2009b), some new discoveries about the function of NSP2 have been reported recently (Kappes et al., 2013; Morgan et al., 2013; Song et al., 2019; Chen et al., 2022). Therefore, the effect of the new deletion pattern of NSP2 found in this study needs to be further explored based on an infectious clone platform. At present, many major epidemic PRRSV strains in China have formed relatively stable and unique characteristics in the NSP2 region, such as a discontinuous 30 aa (1 aa+29 aa) deletion in Nsp2 that is used as the gene marker distinguishing HP-PRRSV from other strains (Tian et al., 2007). Similarly, NADC30-like PRRSV and NADC34-like PRRSV have discontinuous 131 aa deletion (111 aa+1 aa+19 aa) and 100 aa deletion in Nsp2, respectively (Xu et al., 2020). Therefore, excluding recombinant strains such as JS18-3 and HB2104, this new-branch PRRSV with a 1 + 8 + 1 aa stable deletion in NSP2 may represent a necessary condition for becoming a mainstream strain. This further illustrates that the Nsp2 region of sublineage 8.7 PRRSV displays more complex deletion patterns in China.

FIGURE 2

Figure 2. The deduced NSP2 amino acid sequences of the new-branch strains compared with the reference PRRSV strains. The position of the deduced Nsp2 amino acid sequence relative to the VR-2332 strain. Purple indicates the common aa deletion of intermediate PRRSV, new-branch PRRSV and HP-PRRSV. Red indicates the deletion pattern of HP-PRRSV; green indicates the novel characteristic discontinuous deletion pattern of new-branch PRRSV; yellow indicates the deletion pattern of NADC30 PRRSV (not fully displayed); cyan indicates other deletion patterns of intermediate PRRSV strains. The new-branch PRRSV strains obtained in this study are labeled with red triangles (▲). The strains obtained from GenBank belonging to the new branch are marked with blue squares (■).

Homology analysis of full-genome

In view of the unique characteristics of GP5 and NSP2, we further analyzed the full-length genomic characteristics of these PRRSV strains belonging to the new branch. The full-length genome sequences of 3 new PRRSV strains (SDWH86, SDQD95, and SDYT91) consisted of 15,378 nt, a 188-nt 5′-untranslated region (5’-UTR) and a 148-nt 3′-UTR. The full-length genomic nucleotide sequence homology percentages of these sequences with representative PRRSV strains HuN4, CH-1a, HB-1(sh)/2002, VR-2332, NADC30 and QYYZ were 94.0–94.2%, 91.6–91.7%, 92.8–92.9%, 87.5–87.7%, 83.6–83.8% and 85.1–86.1%, respectively. The 5′-UTR, ORF1a, ORF1b, ORF5, and ORF7 shared the highest homology with the HP-PRRSV HuN4 strain; ORF2a, ORF4, ORF5a and ORF6 shared the highest homology with the intermediate PRRSV HB-1(sh)/2002 strain; ORF2b shared the highest homology with the Classical PRRSV CH-1a strain; ORF3 shared the highest homology with the QYYZ strain; and the 3′-UTR shared the highest homology with the NADC30 strain (Table 3). Considering the high homology resulting from the evolutionary relationships among Classical PRRSV, intermediate PRRSV and HP-PRRSV strains (An et al., 2010), we speculated that these strains may show recombination in the ORF3 region.

TABLE 3

Table 3. Nucleotide and amino acid sequence similarity between new PRRSV and reference strains.

Recombination and phylogenetic analysis of full-genome

Indeed, recombination analysis based on RDP4 and SimPlot software verified our conjecture. The results of RDP4 analysis showed that the new-branch strains (SDWH86, SDQD95, and SDYT91) were recombinant viruses with two different sources, among which HP-PRRSV HuN4 was the major parental strain, while QYYZ was the minor parental strain (Table 4). It is worth noting that the RDP4 software analysis showed that HP-PRRSV was the major parental strain, but the homology comparison results showed that these new PRRSV strains presented some degree of homology with the intermediate strains, and their NSP2 region did not have the same characteristics as that of HP-PRRSV. Therefore, their major parents may have evolved from intermediate strains and been homologous to HP-PRRSV, but we have not yet found them. And SimPlot results were similar to those of RDP4. Taking the SDWH86 strain as an example, two recombination breakpoints identified based on RDP4 were used to divide the whole genome into three regions: 1–12,885 (a), 12,885–13,056 (b), and 13,056–15,378 (c) (Figure 3A). Phylogenetic analysis based on breakpoints provided consistent results. The a + c segment belonged to sublineage 8.7 (Figure 3B), and the b segment belonged to sublineage 3.5 (Figure 3C). Therefore, SDWH86, SDQD95 and SDYT91 are recombinant viruses positioned between sublineage 8.7 and 3.5. However, the ZJXS1412, HB2104 and JS18-3 strains on the new branch have more complex recombination patterns (Supplementary Table S1). ZJXS1412 is a sublineage 8.7, 3.5 and 5.1 recombinant; JS18-3 is a sublineage 8.7, 3.5 and 1.8 recombinant; HB2104 is a sublineage 8.7, 3.5, 1.8 and 5.1 recombinant (Supplementary Figure S1). Further phylogenetic analysis based on full-length genomes showed that the SDWH86, SDQD95, SDYT91, ZJXS1412, HB2104, and JS18-3 isolates formed independent branches in sublineage 8.7 (Figure 4A). Additionally, phylogenetic analysis based on the GP3 gene showed that they all belonged to sublineage 3.5, which was a common characteristic of the new-branch strains (Figure 4B). Our previous study described the possible mechanism of complex recombination pattern formation (Xiang et al., 2022). Therefore, we speculated that the complex recombination pattern of the new-branch PRRSV ZJXS1412, HB2104, and JS18-3 strains evolved from the parental strain with sublineage 8.7 and 3.5 recombination, which provided skeletons or fragments for recombination and continuously produced more complex recombination patterns.

TABLE 4

Table 4. Information of recombination events of three new PRRSV detected by RPD4 software.

FIGURE 3

Figure 3. Genome recombination analysis of the SDWH86 isolate. (A) Similarity plots generated with the representative strains HuN4, VR-2332 and QYYZ by Simplot v3.5.1; (B) Phylogenetic trees based on recombinant fragment a + c within SDWH86 and reference PRRSV strains; (C) Phylogenetic trees based on recombinant fragment b within SDWH86 and reference PRRSV strains. The blue region (a + c fragment) indicates that the major parent was HP-PRRSV, and the red region (b fragment) indicates that the minor parent was QYYZ. The new-branch PRRSV strains obtained in this study are labeled with red triangles (▲). The strains obtained from GenBank belonging to the same branch are marked with blue squares (■).

FIGURE 4

Figure 4. Phylogenetic analysis of the whole genome and ORF3 gene of new-branch PRRSV. (A) Phylogenetic tree constructed based on the whole genome of new-branch PRRSV isolates and reference strains from different sublineages; (B) Phylogenetic tree constructed based on the ORF3 gene of new-branch PRRSV isolates and reference strains from different sublineages. The new branch PRRSV strains obtained in this study are labeled with red triangles (▲). The strains obtained from GenBank belonging to the new branch are marked with blue squares (■).

Sequence alignment analysis of 5’-UTR and 3’-UTR

Sequence alignment analysis showed that the new-branch PRRSV strains had similar characteristics in their 5′-UTR and 3′-UTR. They not only retained the same deletion pattern characteristics as intermediate PRRSV strains but also showed a trend of gradual evolution (Figure 5). In the 5′-UTR, new-branch PRRSV, intermediate PRRSV and HP-PRRSV all presented a deletion at the 120th nucleotide (Figure 5A). JS18-3 and HB2104 maintained the same deletion at the 120th nucleotide position as most intermediate PRRSV and HP-PRRSV strains, while SDWH86, SDQD95, SDYT91 and ZJXS1412 showed a continuous 2-nucleotide deletion at positions 119 and 120. The continuous deletion of two nucleotides from the 119th to 120th positions in the 5′-UTR has been previously reported (Zhou et al., 2011; Shi et al., 2013). Since the PRRSV 5′-UTR plays a crucial role in replication, mRNA transcription and protein translation (Gao et al., 2012), the effect of additional deletions in the 5′-UTR needs further study. In the 3′-UTR, new-branch PRRSV, intermediate PRRSV and HP-PRRSV all had a deletion at the 19th nucleotide (Figure 5B). Even more interestingly, SDWH86, SDQD95 and SDYT91 showed further deletion of the 20th and 40th nucleotides. This is the first report of a novel deletion pattern in the 3′-UTR of sublineage 8.7 PRRSV in China. JS18-3, HB2104 and ZJXS1412 showed further deletion of the 32nd nucleotide on the basis of the 19-20th and 40th nucleotide deletions. In addition, the new-branch PRRSV strains presented a completely different motif (ATGA) from the 117-120th nucleotides of the 3′-UTR compared to the sublineage 8.7 strains (Figure 5B). The most recent research indicated that the motif of the 117-120th nucleotides in the 3′-UTR is quite conserved within each sublineage of PRRSV; for example, in most sublineage 8.7 strains, this motif is AAAG (in HP-PRRSVs) or GAGA (in Classical PRRSVs) (Xiong et al., 2022). Some mutations in this motif have been confirmed to enhance the replication ability of HP-PRRSV or L1 PRRSV (Xiong et al., 2022). Therefore, we speculated that the new-branch strains with the same motif in the 3′-UTR may have the same origin; although this motif is inconsistent with those of the sublineage 8.7 strains, it may be a result of evolution. This finding contradicts our previous understanding of sublineage 8.7 PRRSV once again.

FIGURE 5

Figure 5. The deduced 5’-UTR and 3’-UTR nucleotide sequences of the new-branch strains compared with the reference PRRSV strains. (A) The positions of the deduced 5’-UTR nucleotide sequences relative to the VR-2332 strain. (B) The positions of the deduced 3’-UTR nucleotide sequences relative to the VR-2332 strain. Purple indicates the common nt deletions of intermediate PRRSV, new-branch PRRSV and HP-PRRSV, while green indicates new deletions on this basis; a red box indicates a highly conserved motif at positions 117–120 in the 3’-UTR, yellow represents the motif of the new-branch PRRSV, and cyan represents the motif of the Classical strain, intermediate strain and HP-PRRRSV. The new-branch PRRSV strains obtained in this study are labeled with red triangles (▲). The strains obtained from GenBank belonging to the new branch are marked with blue squares (■).

The origin and evolution analysis of the new-branch PRRSV

To deepen the understanding of the origin and genetic relationships of the new-branch PRRSV strains and identify common points in their genome evolution shared with sublineage 8.7 PRRSV, we further collated and analyzed all available the information of the new-branch strains (Supplementary Table S1) and compared the characteristics of the new-branch PRRSV and sublineage 8.7 PRRSV strains (Classical PRRSV, intermediate PRRSV, and HP-PRRSV) from China (Table 5). Previous studies have shown that HP-PRRSV originated from the Classical CH-1a-like PRRSV. This process is formed on the gradual accumulation of intermediate PRRSV strains (An et al., 2010). The new-branch PRRSV strains shared many common characteristics with intermediate PRRSV in the 5′-UTR, 3′-UTR and NSP2 and displayed a gradual evolution trend. The regular deletions in these regions indicated that the new-branch PRRSV strains may have evolved from intermediate PRRSV, and these deletions might have evolved in a step-by-step manner. Furthermore, the results of nucleotide and amino acid homology analyses showed that the new-branch PRRSV strains presented higher homology with the HP-PPRSV strains. However, the NSP2 region of most new-branch strains did not have the same characteristics as that of HP-PPRSV strains except for ZJXS1412, while ORF2a, ORF4, ORF5a and ORF6 of the new-branch strains shared the highest homology with the intermediate PRRSV HB-1(sh)/2002 strain. Therefore, we speculated that the new-branch PRRSV strains may have the same origin and be similar to HP-PRRSV also evolved from intermediate PRRSV, but they constituted separate strains that evolved simultaneously with HP-PRRSV. The new-branch strains showed three different evolutionary directions: some strains, such as ZJXS1412, evolved the same 1 aa+29 aa deletion pattern as HP-PRRSV; some strains, such as SDWH86, SDQD95, and SDYT91, evolved a completely different deletion pattern (1 aa+8 aa+1 aa) from HP-PRRSV but also maintained a close relationship with HP-PRRSV; and some strains, such as HB2104 and JS18-3, displayed a complex recombination pattern. Based on the limited background information obtained from GenBank, these strains of this new-branch could be traced back to 2011. They persist in some areas of China, continue to evolve based on new mutation and recombination events and have the potential to become epidemic strains. At present, new-branch strains have appeared in six provinces, including Shandong, Zhejiang, Jiangsu, Shanghai, Guangdong and Hubei. Unfortunately, the new-branch strains found in this study could not be isolated from primary PAM cells or passaged Marc-145 cells, and their pathogenicity needs further verification. To date, only the JS18-3 strain has been systematically evaluated. It exhibits a relatively high replication ability and severe cytopathic effects in vitro and presents higher virulence in pigs than low-pathogenicity CH-1a-like field strains (Han et al., 2020). This result shows similarity to the characteristics of reported intermediate strains in terms of pathogenicity and cell tropism (Yang et al., 2008; Li et al., 2009; Ma et al., 2009). Therefore, the possibility of the new-branch strains further producing new strains with stronger pathogenicity based on mutation and recombination cannot be ruled out. This highlights the importance of the continuous monitoring of sublineage 8.7 PRRSV in China, especially these strains with novel characteristics.

TABLE 5

Table 5. Comparison of characteristics between new-branch PRRSV and sublineage 8.7 PRRSV in China.

Conclusion

In summary, we identified a new type of PRRSV strain with novel deletion patterns in the 5′-UTR, 3′-UTR and NSP2 regions, which showed the same origin and formed a new independent branch in sublineage 8.7. They may be similar to HP-PRRSV also evolved from intermediate PRRSV, but represent separate strains that evolved simultaneously with HP-PRRSV. In addition, they have undergone rapid evolution and recombined with other strains and have the potential to become epidemic strains, which should receive more attention.

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.

Ethics statement

This study was approved by the Animal Ethics Committee of the School of Harbin Veterinary Research Institute of the Chinese Academy of Agricultural Sciences and was performed in accordance with animal ethics guidelines and approved protocols. The Animal Ethics Committee approval number was SYXK (Hei) 2011022.

Author contributions

HZ, QW, and ZT: conceived and designed the experiments. WL, CLi, and ZG: performed the experiments. HX, BG, QS, JZ, LX, CLe, JP, GZ, YT, and HL: sample collection and data curation. WL, HZ, and QW: contributed to the writing of the manuscript. TA, XC, and ZT: reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by grants from the National Key Research and Development Program (grant no. 2022YFF0711004), the Natural Science Foundation of Heilongjiang Province (grant no. YQ2022C042), the National Natural Science Foundation of China (grant nos. 32002315 and 32172890), the State Key Laboratory of Veterinary Biotechnology Foundation (grant no. SKLVBF202208) and the National Center of Technology Innovation for Pigs.

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/fmicb.2023.1186322/full#supplementary-material

References

Albina, E. (1997). Epidemiology of porcine reproductive and respiratory syndrome (PRRS): an overview. Vet. Microbiol. 55, 309–316. doi: 10.1016/s0378-1135(96)01322-3

PubMed Abstract | CrossRef Full Text | Google Scholar

An, T. Q., Tian, Z. J., Xiao, Y., Li, R., Peng, J. M., Wei, T. C., et al. (2010). Origin of highly pathogenic porcine reproductive and respiratory syndrome virus, China. Emerg Infect Dis. 16, 365–367. doi: 10.3201/eid1602.090005

PubMed Abstract | CrossRef Full Text | Google Scholar

Brinton, M. A., Gulyaeva, A. A., Balasuriya, U. B. R., Dunowska, M., Faaberg, K. S., Goldberg, T., et al. (2021). ICTV virus taxonomy profile: arteriviridae 2021. J. Gen. Virol. 102:001632. doi: 10.1099/jgv.0.001632

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, P., Tan, X., Lao, M., Wu, X., Zhao, X., Zhou, S., et al. (2021). The novel PRRSV strain HBap4-2018 with a unique recombinant pattern is highly pathogenic to piglets. Virol. Sin. 36, 1611–1625. doi: 10.1007/s12250-021-00453-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Yu, L., Zhou, Y., Yang, S., Bai, Y., Wang, Q., et al. (2022). Nonstructural protein 2 is critical to infection efficiency of highly pathogenic porcine reproductive and respiratory syndrome virus on PAMs and influence virulence in vivo. Viruses 14:2613. doi: 10.3390/v14122613

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Z. Q., Guo, X., and Yang, H. C. (2004). Genomic characterization of two Chinese isolates of porcine respiratory and reproductive syndrome virus. Arch. Virol. 149, 1341–1351. doi: 10.1007/s00705-004-0292-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, F., Lu, J., Yao, H., Wei, Z., Yang, Q., and Yuan, S. (2012). Cis-acting structural element in 5' UTR is essential for infectivity of porcine reproductive and respiratory syndrome virus. Virus Res. 163, 108–119. doi: 10.1016/j.virusres.2011.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, J. C., Xiong, J. Y., Ye, C., Chang, X. B., Guo, J. C., Jiang, C. G., et al. (2017). Genotypic and geographical distribution of porcine reproductive and respiratory syndrome viruses in mainland China in 1996-2016. Vet. Microbiol. 208, 164–172. doi: 10.1016/j.vetmic.2017.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, B. Q., Chen, Z. S., Liu, X. W., and Cui, Y. Z. (1996). Isolation and identification of porcine reproductive and respiratory syndrome (PRRS) virus. Chin. J. Animal Poult. Infect. Dis. 25, 1–5.

Google Scholar

Han, G., Lei, K., Xu, H., and He, F. (2020). Genetic characterization of a novel recombinant PRRSV2 from lineage 8, 1 and 3 in China with significant variation in replication efficiency and cytopathic effects. Transbound. Emerg. Dis. 67, 1574–1584. doi: 10.1111/tbed.13491

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Li, G., Yu, L., Li, L., Zhang, Y., Zhou, Y., et al. (2020). Genetic diversity of porcine reproductive and respiratory syndrome virus (PRRSV) from 1996 to 2017 in China. Front. Microbiol. 11:618. doi: 10.3389/fmicb.2020.00618

PubMed Abstract | CrossRef Full Text | Google Scholar

Kappes, M. A., Miller, C. L., and Faaberg, K. S. (2013). Highly divergent strains of porcine reproductive and respiratory syndrome virus incorporate multiple isoforms of nonstructural protein 2 into virions. J. Virol. 87, 13456–13465. doi: 10.1128/JVI.02435-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Fang, L., Xu, Z., Liu, S., Gao, J., Jiang, Y., et al. (2009). Recombination in vaccine and circulating strains of porcine reproductive and respiratory syndrome viruses. Emerg. Infect. Dis. 15, 2032–2035. doi: 10.3201/eid1512.090390

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Xu, H., Zhao, J., Gong, B., Sun, Q., Xiang, L., et al. (2022). Epidemiological investigation and genetic evolutionary analysis of PRRSV-1 on a pig farm in China. Front. Microbiol. 13:1067173. doi: 10.3389/fmicb.2022.1067173

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. K., Zhou, X., Zhai, J. Q., Li, B., Wei, C. H., Dai, A. L., et al. (2017). Emergence of a novel highly pathogenic porcine reproductive and respiratory syndrome virus in China. Transbound. Emerg. Dis. 64, 2059–2074. doi: 10.1111/tbed.12617

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, P., Cai, X., Liu, Y., Shi, W., Wang, S., Wang, H., et al. (2009). Isolation and characterization of PRRSV GD3 strain isolated in China. Chin. J. Prev. Vet. Med. 31, 28–33.

Google Scholar

Morgan, S. B., Graham, S. P., Salguero, F. J., Sanchez Cordon, P. J., Mokhtar, H., Rebel, J. M., et al. (2013). Increased pathogenicity of European porcine reproductive and respiratory syndrome virus is associated with enhanced adaptive responses and viral clearance. Vet. Microbiol. 163, 13–22. doi: 10.1016/j.vetmic.2012.11.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelsen, C. J., Murtaugh, M. P., and Faaberg, K. S. (1999). Porcine reproductive and respiratory syndrome virus comparison: divergent evolution on two continents. J. Virol. 73, 270–280. doi: 10.1128/JVI.73.1.270-280.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, M., Lam, T. T., Hon, C. C., Hui, R. K., Faaberg, K. S., Wennblom, T., et al. (2010). Molecular epidemiology of PRRSV: a phylogenetic perspective. Virus Res. 154, 7–17. doi: 10.1016/j.virusres.2010.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, K., Tao, C., Lin, C., Li, G., and Li, J. (2013). Complete genome sequence of a porcine reproductive and respiratory syndrome virus variant with a new deletion in the 5′ untranslated region. Genome Announc. 1:e00090-12. doi: 10.1128/genomeA.00090-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, J., Gao, P., Kong, C., Zhou, L., Ge, X., Guo, X., et al. (2019). The nsp2 Hypervariable region of porcine reproductive and respiratory syndrome virus strain JXwn06 is associated with viral cellular tropism to primary porcine alveolar macrophages. J. Virol. 93:e01436-19. doi: 10.1128/JVI.01436-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, S., Xu, H., Zhao, J., Leng, C., Xiang, L., Li, C., et al. (2020). Pathogenicity of NADC34-like PRRSV HLJDZD32-1901 isolated in China. Vet. Microbiol. 246:108727. doi: 10.1016/j.vetmic.2020.108727

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, K., Yu, X., Zhao, T., Feng, Y., Cao, Z., Wang, C., et al. (2007). Emergence of fatal PRRSV variants: unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark. PLoS One 2:e526. doi: 10.1371/journal.pone.0000526

PubMed Abstract | CrossRef Full Text | Google Scholar

Tong, G. Z., Zhou, Y. J., Hao, X. F., Tian, Z. J., An, T. Q., and Qiu, H. J. (2007). Highly pathogenic porcine reproductive and respiratory syndrome, China. Emerg Infect Dis. 13, 1434–1436. doi: 10.3201/eid1309.070399

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, L., Xu, H., Li, C., Tang, Y. D., An, T. Q., Li, Z., et al. (2022). Long-term genome monitoring retraces the evolution of novel emerging porcine reproductive and respiratory syndrome viruses. Front. Microbiol. 13:885015. doi: 10.3389/fmicb.2022.885015

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, J., Cui, X., Zhao, K., Wang, Q., Huang, X., Li, D., et al. (2022). A novel motif in the 3'-UTR of PRRSV-2 is critical for viral multiplication and contributes to enhanced replication ability of highly pathogenic or L1 PRRSV. Viruses 14:166. doi: 10.3390/v14020166

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H., Li, C., Li, W., Zhao, J., Gong, B., Sun, Q., et al. (2022). Novel characteristics of Chinese NADC34-like PRRSV during 2020-2021. Transbound. Emerg. Dis. 69, e3215–e3224. doi: 10.1111/tbed.14485

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H., Song, S., Zhao, J., Leng, C., Fu, J., Li, C., et al. (2020). A potential endemic strain in China: NADC34-like porcine reproductive and respiratory syndrome virus. Transbound. Emerg. Dis. 67, 1730–1738. doi: 10.1111/tbed.13508

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, M., Xing, G., Sun, W., Huang, J., Jiang, X., and Zhou, J. (2008). Isolation and pathogenicity analysis of porcine reproductive and respiratory syndrome virus. Chin. J Prev. Vet. Med. 30, 846–851.

Google Scholar

Yu, F., Yan, Y., Shi, M., Liu, H. Z., Zhang, H. L., Yang, Y. B., et al. (2020). 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. 94:e01813-19. doi: 10.1128/JVI.01813-19

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H. L., Zhang, W. L., Xiang, L. R., Leng, C. L., Tian, Z. J., Tang, Y. D., et al. (2018). Emergence of novel porcine reproductive and respiratory syndrome viruses (ORF5 RFLP 1-7-4 viruses) in China. Vet. Microbiol. 222, 105–108. doi: 10.1016/j.vetmic.2018.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Chen, S., Zhang, J., Zeng, J., Guo, X., Ge, X., et al. (2009a). Molecular variation analysis of porcine reproductive and respiratory syndrome virus in China. Virus Res. 145, 97–105. doi: 10.1016/j.virusres.2009.06.014

CrossRef Full Text | Google Scholar

Zhou, Z., Ni, J., Cao, Z., Han, X., Xia, Y., Zi, Z., et al. (2011). The epidemic status and genetic diversity of 14 highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) isolates from China in 2009. Vet. Microbiol. 150, 257–269. doi: 10.1016/j.vetmic.2011.02.013

CrossRef Full Text | Google Scholar

Zhou, L., and Yang, H. (2010). Porcine reproductive and respiratory syndrome in China. Virus Res. 154, 31–37. doi: 10.1016/j.virusres.2010.07.016

CrossRef Full Text | Google Scholar

Zhou, L., Zhang, J., Zeng, J., Yin, S., Li, Y., Zheng, L., et al. (2009b). 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. 83, 5156–5167. doi: 10.1128/JVI.02678-08

PubMed Abstract | CrossRef Full Text | Google Scholar