Genetic variation and recombination analysis of PRRSV-2 GP3 gene in China from 1996 to 2023

Porcine reproductive and respiratory syndrome virus (PRRSV) has become widespread in China particularly the highly pathogenic porcine reproductive and respiratory syndromes (HP-PRRSV), NADC30, and NADC34 strains, and has posed a threat to the swine industry for over 20 years. To monitor genetic variation in PRRSV-2 GP3 strains in China, we analyzed 618 strains isolated between 1996 to 2023 and constructed phylogenetic trees. Additionally, 60 selected strains were used to analyze nucleotide and amino acid homology. PRRSV GP3 gene exhibited nucleotide identity ranging from 78.2% to 100.0% and amino acid similarity ranging from 74.9% to 99.6%. The GP3 gene in the 60 selected strains consisted of 254 amino acids, and amino acid mutations in the strains primarily occurred in B-cell epitopes, T-cell epitopes, and highly variable regions. The glycosylation sites of the strains used for amino acid sequence comparisons remained unaltered, except for the N²⁹ site in the GD20220303-2022 strain. PRRSV-2 strains in China belong to lineages 1, 3, 5, and 8. Recombination analysis detected two recombination events, involving lineages 1 and 8. In conclusion, this study investigated multiple strains of the PRRSV-2 GP3 gene to explore the prevalence and genetic diversity of the GP3 gene in China from a gene family perspective. The results of the analyses provide a basis for clinical prevention strategies and vaccine development.

1 Introduction

Porcine reproductive and respiratory syndrome (PRRS) is an acute and virulent infectious disease caused by the porcine reproductive and respiratory syndrome virus (PRRSV) (Yu et al., 2022). The typical symptom of PRRS is the cyanosis of the skin at the end of the ears, commonly known as the “blue-ear disease.” Pregnant sows and piglets are the most susceptible to this virus. The symptoms in pregnant sows after infection are complicated and manifest primarily as reproductive disorders, abortion, weak fetuses, still-births, and mummified fetuses, accompanied by immunosuppressive diseases and certain respiratory diseases. In piglets, the syndrome is manifested as congenital dysplasia, hyperthermia, and interstitial pneumonitis. PRRS was initially discovered in the United States in 1986 and was called a “mystery disease”. Subsequently, it spread rapidly worldwide. This disease was initially reported in China in 1996 and subsequently spread widely and mutated into highly pathogenic strains (HP-PRRSV) in 2006 and NADC-like strains in 2012 (Snijder and Meulenberg, 1998; Gao et al., 2017; Zhang et al., 2018a; Xie et al., 2020; Liu et al., 2021; Xu et al., 2022). PRRS is highly contagious with a high mortality rate, causing huge economic losses to large-scale piggeries worldwide, and is considered one of the most important infectious diseases threatening the sustained and healthy development of the global animal breeding industry, together with various virulent infectious diseases of animals (Allende et al., 1999; Tian et al., 2007; Zhou et al., 2018; Cai et al., 2022; Liang et al., 2022; Goatley et al., 2023).

PRRSV is classified within the family Arterivirus, under the order Nidoviridae (Cavanagh, 1997). Other viruses in this family include equine arteritis virus (EAV), lactate dehydrogenase-elevating virus (LDV), and simian hemorrhagic fever virus (SHFV) (Shi et al., 2020). PRRSV is a single-stranded, positive-stranded RNA virus with spherical particles, a capsid, a smooth, ciliated surface, and a diameter of approximately 50–60 nm. The nucleocapsid is cubic in shape and has a diameter of 25–35 nm. The virus is temperature-sensitive and can be inactivated by heating to 56°C for 20 min. Organic chemicals can also destroy the virus, and chloroform, ether, and trichloromethane can inhibit its reproduction by destroying its capsid (Tong et al., 2007). PRRSV can be classified into two genotypes, namely, type 1 (European) and type 2 (American). The two genome sequences are similar by approximately about 54–60% (Meng et al., 1995), and the prevalent strains in Europe are mostly PRRSV-1 (symbolized by LV strains extracted from the Netherlands). In contrast, in China, PRRSV strains are mainly PRRSV-2 (symbolized by the VR2332 strain from the United States). Per the recent classification scheme of the International Committee for the Classification of Viruses (ICVC), PRRSV can be categorized into Betaarterivirus suid 1 (former PRRSV-1 type) and Betaarterivirus suid 2 (former PRRSV-2 type) under Betaarterivirus (Adams et al., 2017).

The genome of PRRSV spans approximately 15 kb, including at least 11 open reading frames (ORFs): ORF1a, ORF1b, ORF2a, ORF2b, ORF3-7, ORF5a, and TF. TF overlaps with the coding region of non-structural protein 2 (NSP2) of ORF1a (Fang et al., 2012). PRRSV encodes a structural protein (SP) and a non-structural protein (NSP), which can be classified into two types: major structural proteins (GP5, M, and N) and minor structural proteins (GP2, GP2a, GP3, GP4 and E) (Ke and Yoo, 2017). The remaining are non-structural proteins that constitute approximately about 75% of the genome (Ma et al., 2021). ORF1a and ORF1b are responsible for encoding two important viral replicase proteins (pp1a and pp1ab), which can be translated into 16 non-structural proteins (NSPs): NSP1α, NSP1β, NSP2, NSP2TF, NSP2N, NSP3-6, NSP7α, NSP7β and NSP8-12 (Sun et al., 2022). ORF2-7 are responsible for encoding GP2, E, GP3, GP4, GP5a, GP5, M and N (Snijder et al., 1999). These proteins play an important role in activating replicate, shear, and protease activities during replication, transcription, and translation, and some NSPs play important roles in PRRSV infection, anti-body production, and immune escape regulation (Oleksiewicz et al., 2001).

GP3, a protein encoded by ORF3, has a molecular mass of approximately 27-29 kDa and is composed of a segmented signal peptide, a domain with extensive glycosylation, a brief hydrophobic segment, and a C-terminal structural domain that is not glycosylated and varies greatly. GP3 is the most glycosylated protein in PRRSV, with seven highly conserved glycosylation sites, and is the second most conserved structural protein after GP5. Differences in antigenic epitopes and the secondary structures of GP3 sequences between different typologies and strains may be responsible for this phenomenon. GP3 also has high immunogenicity and antigenicity, which can induce neutralizing antibody production and interact with CD163, a key receptor for PRRSV entry into the cell, to increase infectivity (Vu et al., 2011; Zhang et al., 2020). GP3 is found in relatively small amounts of viral particles. It is mainly produced in cellular secretions and may act as a decoy to disperse antibodies to keep susceptible cells away from viral particles (Zhang et al., 2018b). Antibodies produced against GP3 have been experimentally demonstrated to prevent viral infections, but the ability of GP3 to generate antibodies is not as strong as that of GP5 (Guarino et al., 1999). Furthermore, the fusion of heat shock protein 70 with GP3 of PRRSV induced the production of interferon-gamma (IFN-γ) and interleukin-4 (IL-4) in porcine sera (Du et al., 2012). IL-4 production in the serum enhances the body’s immune response, protects swine, and reduces the infection rate of PRRSV (Du et al., 2012). The high degree of GP3 glycosylation may influence viral immune evasion by the virus, GP3 can also form a heterotrimer with GP2a and GP4 to dope the virion, and viral particles deficient in these three proteins lose their infectivity (Murtaugh et al., 1995; Wieringa et al., 2002; Vu et al., 2011).

2 Materials and methods

2.1 Data collection

A comprehensive genetic variation analysis of the GP3 gene in Chinese PRRSV-2 was conducted by selectively screening 618 PRRSV-2 strains (610 isolates from China and 8 classical strains from the United States) from the GenBank database available on the NCBI website¹ (Supplementary Table 1). The study included a selection of strains from various years and lineages (lineages 1, 3, 5, and 8). A total of 60 representative strains with GP3 sequences were meticulously chosen from a larger pool of 618 strains, including representatives from various sources such as representative, vaccine, and prevalent strains. This thorough selection process was carried out to ensure the inclusiveness and comprehensiveness of the genetic variation analysis of the GP3 gene in Chinese PRRSV-2 (Table 1).

TABLE 1

Table 1. The GP3 reference sequences of 60 PRRS strains.

2.2 PRRSV-2 GP3 gene sequence analysis

GP3 nucleotide homology was examined using the Clustal W tool within the MegaAlign feature of DNAStar software (version 7.0, Madison, WI, USA), and the nucleotide sequences were then translated into the corresponding amino acid sequences using the EditSeq program to analyze GP3 amino acid homology. Amino acid sequences were analyzed using the DNAStar and MEGA software (version 5.1, MegaLimited, Auckland, New Zealand).

2.3 Phylogenetic analysis

Phylogenetic analysis of the 618 PRRSV strains was performed using the GP3 sequences. Phylogenetic trees were constructed using the maximum likelihood (ML) and neighbor-joining (NJ) methods in the MEGA software with 1000 times self-expansion analyses and default settings for other parameters. The obtained phylogenetic tree results were imported into the online visualization and editing software The Interactive Tree of Life² for annotation and editing.

2.4 Recombination analysis

The RDP4 software employs seven distinct algorithms (RDP, GENECONV, BootScan, MaxChi, Chimera, SiScan, and 3Seq) to forecast possible recombination occurrences. Strains whose analyses were positive (+) for more than four algorithms with p < 0.05 were considered recombinant, assuming that recombination events had occurred. Strains identified as having undergone recombination were validated using SimPlot software (version 3.5.1).

3 Results

3.1 Nucleotide and amino acid homology analysis of GP3 gene

To investigate genetic variation in the PRRSV GP3 gene, 60 strains from lineages 1, 3, 5, and 8 were selected. Nucleotide and amino acid homology were analyzed using DNAStar software to fully comprehend the genetic evolutionary relationships among these different lineages. The findings indicated that the nucleotide similarity among the GP3 genes of the selected strains varied between 78.2 to 100% whereas the amino acid similarity ranged from 74.9 to 99.6%.

The lowest nucleotide homology (78.2%) was found in GD20220303-2022 and R98-2006 strains, while strains GD20220303-2022 and R98-2006 had the lowest amino acid homology (74.9%). The highest nucleotide homology (100%) was observed between the following set of strains: ZJqz21-2021 and NADC30-2008, YN-2011-2011 and BJ-4-2000, GS2002-2009, GS2003-2009 and GS2004-2009, YN-2011-2011, BJ-4-2000 and MLV RespPRRS Repro-1999, S1-2006 and PRRSV01-2008, ATCC VR-2332-1993 and PRRSV03-2008, CH2002-2009, CH2003-2009 and CH2004-2009, rJXA1-R-2020 and JXA1 P80-2008, YN-DZ-2023 and JXA1-2006, QY2010-2010 and QYYZ-2011 strains. Additionally, the amino acid similarity between strains NADC30-2008 and ZJqz21-2021, BJ-4-2000 and MLV RespPRRS Repro-1999, S1-2006 and PRRSV01-2008, ATCC VR-2332-1993 and PRRSV03-2008, GS2002-2009, GS2003-2009 and GS2004-2009, BJ-4-2000, MLV RespPRRS Repro-1999 and YN-2011-2011, TJ-2006 and BJ-2007, CH2002-2009, CH2003-2009 and CH2004-2009, QY2010-2010 and QYYZ-2011 were 99.6% (Table 2).

TABLE 2

Table 2. Nucleotide and amino acid homology analysis of GP3 gene (%).

3.2 Comparison of GP3 amino acid sequence

MegAlign feature in DNAStar software was utilized to analyze the 60 amino acid mutation sites in PRRSV GP3. Across all 60 strains selected for sequence comparison, the GP3 protein consistently had 254 amino acids. Notably, mutations were primarily concentrated around residues 2–32 and 200–230, encompassing the B-cell epitope region, T-cell epitope region, and high variant region (HVR). It is worth mentioning that the GD20220303 strain exhibited an N²⁹ → D²⁹ mutation at the glycosylation site of amino acid N²⁹, while the remaining 59 strains showed no mutations in the seven glycosylation sites (N²⁹, N⁴², N⁵⁰, N¹³¹, N¹⁵², N¹⁶⁰, and N¹⁹⁵). The 7 strains of lineage 1 showed the highest concentration of mutations at positions 2-32 compared with those of lineages 3, 5, and 8, and the 19 strains of lineage 5 showed the lowest concentration of mutations compared with those of lineages 1, 3, and 8. All 60 strains showed the F³² → Y³² mutation at position 32, with the lowest concentration of mutations at positions 110–130 and 170–200; the highest concentration of mutations was found in the epitope region of the B-cells. The B-cell epitope region exhibited the highest variation. Some strains had characteristic mutations according to the lineage, such as T⁶ → A⁶, V²⁵ → A²⁵, and S⁸³ → E⁸³ in lineage 1, R⁷¹ → G⁷¹ in lineages 5, and R⁷¹ → K⁷¹, L⁷³ → F⁷³, and Q²¹⁴ → R²¹⁴ and 8 (Figure 1).

FIGURE 1

Figure 1. SimPlot validation of the recombination events in the GP3 gene. (A) Virulent strain GDDX-2018-2018 in which a recombination event occurred. (B) Virulent strain FJ1402-2014 in which a recombination event occurred. (C) Virulent strain CH-SCMY-2-2019-2019 in which a recombination event occurred. The blue line represents the primary parental strain, the yellow line represents the secondary parental strain, and the red line represents the control strain CH-1a.

3.3 Phylogenetic analysis

According to the worldwide PRRSV categorization framework, 618 PRRSV-2 strains were identified from the NCBI website. Analysis of GP3 sequences sourced from the GenBank database revealed that the common PRRSV-2 strains in China belong to lineages 1, 3, 5, and 8, with lineage 8 being closer to lineages 3 and 5 and farthest from lineage 1 by the ML method, but lineage 8 is more closely related to lineage 1 and more distantly related to lineages 3 and 5 by the NJ method. The recent dominant strains and mostly newly emerged strains in the last three years belong to lineages 1, 5, and 8 (Figures 2, 3).

FIGURE 2

Figure 2. Phylogenetic analysis of the PRRSV GP3 gene of 618 strains. The phylogenetic tree was constructed by the ML method using MEGA software (version 5.1, MegaLimited, Auckland, New Zealand), and self-expansion analysis was performed 1000 times.

FIGURE 3

Figure 3. Phylogenetic analysis of the PRRSV GP3 gene of 618 strains. The phylogenetic tree was constructed by the NJ method using MEGA software (version 5.1, MegaLimited, Auckland, New Zealand), and self-expansion analysis was performed 1000 times.

3.4 Recombination analysis

The RDP4 software was employed to predict potential recombination occurrences, which were then confirmed using SimPlot. Three possible recombination events were identified in GP3, two of which were verified by SimPlot with p < 0.05, proving that the results were statistically significant (Figure 4). Event A occurred in lineages 3 and 8, indicating recombination between lineage 3 and 8, events C occurred in lineages 1 and 3, indicating a recombination event between these two lineages 1 and 3; and event B occurred exclusively in lineage 8, suggesting recombination within lineage 8 (Table 3). RDP4 detected recombination in the CH-SCMY-2-2019-2019 and GDDX-2018-2018 strains belonging to lineages 1 and 8 respectively. The recombinant GDDX-2018-2018 strain had HEB-108-2017 strain as the primary parental strain, and PRRSV2-CN-F2-2019-2019 as the secondary parental strain. For the CH-SCMY-2-2019-2019 strain, the primary and secondary parental strains were SCLS-2-2020 and PRRSV2-CN-N0-2021-2021, respectively; SimPlot confirmed the occurrence of these two recombination events (Figure 5).

FIGURE 4

Figure 4. Prediction of genetic recombination events in GP3 using RDP4 software. (A) Strain GDDX-2018-2018 in which recombination events occurred. (B) Strain FJ1402-2014 in which recombination events occurred. (C) Strain CH-SCMY-2-2019-2019 in which recombination events occurred (R, RDP, G, GENECONV, B, BootScan, M, MaxChi, C, Chimaera, S, SiScan, T, 3Seq).

TABLE 3

Table 3. Recombination analysis of the GP3 gene.

FIGURE 5

Figure 5. Amino acid sequences of the 60 PRRSV GP3 strains were compared with each other in lineages 1, 3, 5, and 8, with the red region representing B-cell epitopes, the blue region representing T-cell epitopes, the purple region representing glycosylation sites, and the green region representing amino acid sequences inserted.

4 Discussion

Following the initial identification of PRRSV in China in 1996, continuous mutation and genetic recombination events have resulted in its extensive proliferation, presenting a significant hurdle for the swine farming sector. GP3, with its high immunogenicity and antigenicity, plays an important role in PRRSV infection and the production of neutralizing antibodies (Wang et al., 2014). To study more deeply the variation and genetic evolution of the PRRSV-2 GP3 gene in China, we analyzed the nucleotide and amino acid homology of 60 representative strains from different lineages, as well as their amino acid sequence comparisons. Our findings showed that the nucleotide homology of GP3 ranged from 78.2% to 100.0%, while the amino acid homology ranged from 74.9 to 99.6%. The Chinese classical CH-1a strain had the highest nucleotide and amino acid homology with the MLV RespPRRS Repro strain compared to the eight classical strains from the United States, with 91.6% and 88.6%, respectively, and also had high homology with the newly emerged YN-DZ strain from China in 2023, with 95.8% of the nucleotides and 92.5% of the amino acids, respectively 92.5%. Through the homology analysis of these 60 strains (Table 2), it was observed that lineage 5 strains displayed the most significant differences in nucleotide and amino acid homology. This suggests notable recombination and variation within lineage 5, indicating the presence of genetic mutations in the GP3 gene. Lineage 3 strains emerged later in China, with a lower clinical detection rate. In contrast, lineage 5 strains, represented by strains such as BJ-4, had a comparatively longer period of circulation in China than lineage 3. However, lineage 5 strains exhibited weaker infectivity in natural environments, resulting in a lower clinical detection rate (Guo et al., 2018). Lineage 8 strains, represented by strains such as CH-1a, were the earliest to emerge and have become prevalent in China. Furthermore, lineage 8 strains were further classified into sublineages 8.1, 8.7, and 8.9. In contrast, strains from lineages 1, 3, and 5 became more prevalent in China later than those from lineage 8. Most strains detected in China belong to lineage 8, making it the dominant circulating lineage in the country.

Amino acid deletions, insertions, and substitutions are the key factors affecting mutations in PRRSV. In this study, a sequence comparison of amino acids was performed, revealing several mutations in the PRRSV-2 type GP3 protein. These mutations mainly occurred in the B-cell epitope region, the T-cell epitope region, and the HVR. Although all strains across the spectrum exhibited varying degrees of mutations in these three regions, none of them showed any additions or deletions. The GP3 protein is known to be the most heavily glycosylated protein in PRRSV. Among the 60 strains analyzed, only the first glycosylation site of the GD20220303 strain had a mutation from N²⁹ to D²⁹, while the remaining six glycosylation sites and all glycosylation sites of the other 59 strains remained unmutated. This suggests a high degree of conservation in the glycosylation sites of GP3 across different PRRSV subtypes. Studies have demonstrated functional similarities between the glycosylation sites on GP3 and GP5 proteins, highlighting their crucial role in PRRSV immune evasion (Liu et al., 2011), therefore, further research on these glycosylation sites is required.

Phylogenetic analysis revealed that Chinese native PRRSV-2 strains belong to lineages 1, 3, 5, and 8, with lineages 1 and 8 exhibiting a more pronounced pattern of occurrence. In the constructed phylogenetic tree, strain CH-1a was most closely homologous to HB-1-sh strain from China and shared lineage 8 with the YN-DZ strain; it was most closely homologous to ATCC VR-2332 strain in lineage 5 from the United States. The genetic variations: between the GP3 and the GP5 (GP5a, which is prone to variation among structural proteins) genes, as well as among the NSP2 gene, which undergoes frequent high-level variations and is susceptible to deletions, were compared (Luo et al., 2023a,b; Zhang et al., 2023). In the evolutionary tree generated by the ML method for GP3, lineages 1 and 5 showed a closer relationship, as did lineages 3 and 8. Lineage 5 could be divided into two parts, one part of the strain was located in the same branch as that of lineage 1 and is further away from lineage 8, and the other part of the strain is closer to lineage 3, but both parts are closer to each other than to the other lineages and to both lineages 1. Strains in lineage 3 were closer to lineage 8, suggesting that recombination might have occurred in lineages 3 and 8, a possibility that was also successfully predicted by RDP4 and identified by SimPlot. In contrast, the NJ evolutionary tree indicated that lineages 1 and 5 were closely related, lineages 1 and 3 were more distantly related, and lineages 3 and 8 were genetically distant from each other. Taken together, it can be inferred that the evolutionary tree constructed using the ML method had the highest degree of conformity with the actual situation. In addition, lineages 1 and 8 had a greater predominance in the prevalence of PRRS and have been the most prevalent lineages in China in recent years, which is consistent with the results of the two studies mentioned above (Zhang et al., 2022, 2023; Luo et al., 2023a). Variations in the phylogenetic tree could be linked to the quantity, timing, and diversity of the selected aggressive strains. In this study, we posited that variances in genetic distances could stem from factors such as the number of aggressive strain samples, duration, specific genes, and targeted algorithms, along with the genes and algorithms used for target analysis. In the last three years, most PRRSV-2 strains belonged to lineages 1, 5, and 8. 5 of the 11 new strains that emerged in 2023 belonged to lineage 1 along with the currently prevalent NADC30-like and NADC34-like strains in China, and the numbers of strains belonging to lineages 5 and 8 were 4 and 2, respectively, with no strains belonging to lineage 3, which is the same as the current predominance of PRRSV in China. This is the same as the current predominant PRRSV lineage 1 and 8 results in China. Based on our findings, we hypothesize that the prevalence trend of lineage 1 may continue to expand in the future and become the most prevalent lineage in China (Tian et al., 2007).

Since its identification in China, PRRSV has experienced a swift evolutionary process and has swiftly proliferated across the nation. Initial identification of the PRRSV NADC30 strain occurred in 2008 in the United States. At the end of 2012, the discovery of HENAN-XINX, HNyc15, JL580, and HENAN-HEB strains, which are homology to the NADC30 strain, was reported in China; these were named NADC30-like (NADC30-like) strains, which rapidly spread to several regions of the country and formed independent branches with an increased percentage of clinical detection (Tian, 2017). NADC30-like strains are less pathogenic than HP-PRRSV strains and are prone to recombination with classical strains such as CH-1a and HP-PRRSV strains (Zhou et al., 2018). In 2015, the combined presence of NADC30-like and HP-PRRSV strains increased becoming leading strains in China. In 2017, NADC34-like strains were initially identified in Liaoning Province, and by 2021, NADC34-like strains were reported in several other Chinese provinces (Xu et al., 2022). Reports suggest that NADC34-like strains may emerge as the leading PRRSV epidemic strains in China (Zhao et al., 2022). PRRSV undergoes mutations and recombination, resulting in more complex PRRSV strain types of higher virulence, making PRRS prevention, control, and decontamination more difficult. Natural genetic recombination between strains used in vaccines and those in the wild is possible, and recombination events of the GM2 strain can occur between the recombinant vaccine strain MLV RespPRRS Repro and the wild strain QYYZ, suggesting that recombination events are key to the evolution of PRRSV (Liu et al., 2021).

In this study, three sets of recombination events were identified by performing recombination analysis, and the presence of two of these sets was verified using the SimPlot software. Lineages 1, 3, and 8 had a higher likelihood of recombination, and the recombination strains CH-SCMY-2-2019-2019 and GDDX-2018-2018, both of which had major and minor parental strains, belonged to lineages 1 and 8, respectively. The analysis of the ML evolutionary tree, which aligns well with real-world data, suggests a close relationship between lineages 1, 3, and 8. Specifically, lineages 3 and 8 exhibit a tendency towards recombination mutations, confirming the predictions made by RDP4. Recombination is a reason why PRRSV continues to evolve and adapt to the environment, and its high prevalence is one of its characteristics. In recent years, strains born as a result of continuous recombination have greatly reduced the effectiveness of existing vaccines, which makes the prevention and control of blue-ear disease more challenging. Presently, in China, the predominant strains are NADC30-like and NADC34-like, followed by the HP-PRRSV strains. There are three distinct stages in the epidemiological timeline of PRRS in China, each of which includes representative strains and clinical symptoms. The first phase was from 1996 to 2006, during which the dominant epidemic strain was the C-PRRSV strain represented by the classical strain CH-1a, which was the earliest strain to be identified in China and may have led to the termination of pregnancy in sows and respiratory illnesses in swine bars (Zhao et al., 2022). In the second phase, from 2006 to 2012, HP-PRRSV strains replaced the C-PRRSV strains as the new dominant strains in China (Tian et al., 2007), which caused higher morbidity and mortality rates, more severe clinical manifestations, and affected fattening swine and gilts of all ages, adversely impacting the Chinese economic swine industry. In the final phase, from 2012 to the present, the prevalence of lineage 1 strains increased, rapidly spreading to all parts of the country, forming an independent branch and gradually becoming the new dominant prevalent spectrum in China, represented by the NADC30-like strain (Liu et al., 2021). Furthermore, the management and prevention of PRRS have become more complex owing to the ineffectiveness of existing market vaccines in hindering recombinant strains. In 2017, the intermediate virulence strain NADC34-like was widely prevalent, with clinical manifestations of reproductive disorders in sows, and the prevalence of the NADC34-like strain has increased significantly since 2020, and it may become the dominant prevalent strain in the future (Zhao et al., 2022), further complicating the clinical prevention and treatment of PRRS more difficult in China.

Our research focused on examining the genetic diversity and genetic mixing of 618 PRRSV-2 GP3 genes across China, spanning the years 1996 to 2023, and investigated the similarities and differences in nucleotides and amino acids among various spectrum strains, PRRSV generates new dominant prevalent strains to adapt to the changes in the transmission environment through continuous mutations. Existing vaccines are unable to effectively prevent PRRS infection and transmission, which poses significant challenges to prevention and control efforts. Although the GP3 gene exhibits relatively low conservation, it is presently considered moderately conserved. However, the potential for mutations and recombination events during future transmissions should not be overlooked. Through the analysis of GP3 genes from diverse strains and over different years, our objective is to provide a theoretical framework and roadmap for developing novel vaccines and treatments. This research aims to mitigate the impact of porcine PRRS on the pig farming industry.

5 Conclusion

Between 1996 and 2023, the prevalent PRRSV-2 strains in China have been categorized into lineages 1, 3, 5, and 8. Notably, lineages 1 and 8 have shown extensive distribution, prolonged prevalence, high recombination rates, and increased mutation probabilities. These lineages are also considered the predominant strains among the commonly found variants in China. GP3 is not well conserved and is vulnerable to mutations. It plays a role in viral infection and antibody production, making it a promising research focus for developing new vaccines. By studying the gene sequence of GP3, it is possible to better understand its prevalence and mutation trends in China and screen for key mutation sites, as well as to guide the clinical selection of effective vaccine strains close to the clinic for the prevention and control of PRRS, to reduce all kinds of losses caused by PRRS.

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

CL: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing. YZ: Data curation, Formal analysis, Methodology, Writing – original draft. KL: Data curation, Writing – original draft. GL: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft. QL: Writing – review & editing. HZ: Writing – review & editing. HS: Writing – review & editing. RW: Writing – review & editing. WK: Writing – review & editing. MZ: Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of the article. This research was funded by the Guangdong Provincial Department of Agriculture and Rural Affairs-Guangdong Agricultural Technical Service “Light Cavalry” Project 2024 (NITG20240253), Guangdong Provincial Department of Education’s distinctive innovation initiative (2023KTSCX128) and the National Natural Science Foundation of China (31902279) back this study.

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

Footnotes

1. ^ www.ncbi.nlm.nih.gov
2. ^ https://itol.embl.de

References

Adams, M. J., Lefkowitz, E. J., King, A. M. Q., Harrach, B., Harrison, R., Knowles, N., et al. (2017). Changes to taxonomy and the international code of virus classification and nomenclature ratified by the international committee on taxonomy of viruses (2017). Arch. Virol. 162, 2505–2538. doi: 10.1007/s00705-017-3358-5

PubMed Abstract | Crossref Full Text | Google Scholar

Allende, R., Lewis, T. L., Lu, Z., Rock, D., Kutish, G., Ali, A., et al. (1999). North American and European porcine reproductive and respiratory syndrome viruses differ in non-structural protein coding regions. J. Gen. Virol. 80, 307–315. doi: 10.1099/0022-1317-80-2-307

PubMed Abstract | Crossref Full Text | Google Scholar

Cai, H., Qi, N., Li, J., Lv, M., Lin, X., Hu, J., et al. (2022). Research progress of the avian coccidiosis vaccine. Vet. Vaccine 1:100002. doi: 10.1016/j.vetvac.2022.100002

Crossref Full Text | Google Scholar

Cavanagh, D. (1997). Nidovirales: A new order comprising coronaviridae and arteriviridae. Arch. Virol. 142, 629–633.

Google Scholar

Du, Y. J., Qi, J., Lu, Y., Wu, J., Yoo, D., Liu, X., et al. (2012). Evaluation of a DNA vaccine candidate co-expressing GP3 and GP5 of porcine reproductive and respiratory syndrome virus (PRRSV) with interferon α/γ in immediate and long-lasting protection against HP-PRRSV challenge. Virus Genes 45, 474–487. doi: 10.1007/s11262-012-0790-1

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, Y., Treffers, E. E., Li, Y. H., Tas, A., Sun, Z., van der Meer, Y., et al. (2012). Efficient-2 frameshifting by mammalian ribosomes to synthesize an additional arterivirus protein. Proc. Natl. Acad. Sci. U.S.A. 109, E2920–E2928. doi: 10.1073/pnas.1211145109

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, J. C., Xiong, J. Y., Ye, C., Chang, X., Guo, J., Jiang, C., 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

Goatley, L. C., Tng, P., and Al-Adwani, L. (2023). ASFV antigens selected from genotype I immunised pigs are immunogenic, but do not protect against genotype II challenge. Vet. Vaccine 2:4. doi: 10.1016/j.vetvac.2023.100042

Crossref Full Text | Google Scholar

Guarino, H., Goyal, S. M., Murtaugh, M. P., Morrison, R., and Kapur, V. (1999). Detection of porcine reproductive and respiratory syndrome virus by reverse transcription-polymerase chain reaction using different regions of the viral genome. J. Vet. Diagn. Invest. 11, 27–33. doi: 10.1177/104063879901100104

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Z. H., 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. J. 15:2. doi: 10.1186/s12985-017-0910-6

PubMed Abstract | Crossref Full Text | Google Scholar

Ke, H. Z., and Yoo, D. (2017). The viral innate immune antagonism and an alternative vaccine design for PRRS virus. Vet. Microbiol. 209, 75–89. doi: 10.1016/j.vetmic.2017.03.014

PubMed Abstract | Crossref Full Text | Google Scholar

Liang, W., Meng, X. G., Zhen, Y. R., Zhang, Y., Hu, X., Zhang, Q., et al. (2022). Integration of transcriptome and proteome in lymph nodes reveal the different immune responses to PRRSV between PRRSV-resistant Tongcheng pigs and PRRSV-Susceptible large white pigs. Front. Genet. 13:800178. doi: 10.3389/fgene.2022.800178

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Zhou, Y., and Yin, J. (2011). Functional analysis of glycosylation of GP3 protein of porcine reproductive and respiratory syndrome virus. Chin. J. Anim. Infect. 19:8. doi: 10.3969/j.issn.1674-6422.2011.03.001

Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, Q., Zheng, Y. J., He, Y. X., Li, G., Zhang, H., Sha, H., et al. (2023a). Genetic variation and recombination analysis of the GP5 (GP5a) gene of PRRSV-2 strains in China from 1996 to 2022. Front. Microbiol. 14:1238766. doi: 10.3389/fmicb.2023.1238766

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, Q., Zheng, Y. J., Zhang, H., Yang, Z., Sha, H., Kong, W., et al. (2023b). Research progress on glycoprotein 5 of porcine reproductive and respiratory syndrome virus. Animals 13:813. doi: 10.3390/ani13050813

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, J., Ma, L. L., Yang, M. T., Wu, W., Feng, W., and Chen, Z. (2021). The function of the PRRSV-host interactions and their effects on viral replication and propagation in antiviral strategies. Vaccines 9:364. doi: 10.3390/vaccines9040364

PubMed Abstract | Crossref Full Text | Google Scholar

Meng, X. J., Paul, P. S., Halbur, P. G., and Morozov, I. (1995). Sequence comparison of open reading frames 2 to 5 of low and high virulence United States isolates of porcine reproductive and respiratory syndrome virus. J. Gen. Virol. 76, 3181–3188. doi: 10.1099/0022-1317-76-12-3181

PubMed Abstract | Crossref Full Text | Google Scholar

Murtaugh, M. P., Elam, M. R., and Kakach, L. T. (1995). Comparison of the structural protein-coding sequences of the VR-2332 and Lelystad virus-strains of the PRRS virus. Arch. Virol. 140, 1451–1460. doi: 10.1007/bf01322671

PubMed Abstract | Crossref Full Text | Google Scholar

Oleksiewicz, M. B., Botner, A., Toft, P., Normann, P., and Storgaard, T. (2001). Epitope mapping porcine reproductive and respiratory syndrome virus by phage display: The nsp2 fragment of the replicase polyprotein contains a cluster of B-cell epitopes. J. Virol. 75, 3277–3290. doi: 10.1128/jvi.75.7.3277-3290.2001

PubMed Abstract | Crossref Full Text | Google Scholar

Shi, Y. J., Tong, X. H., Ye, G., Xiu, R., Li, L., Sun, L., et al. (2020). Structural characterization of the Helicase nsp10 encoded by porcine reproductive and respiratory syndrome virus. J. Virol. 94:15. doi: 10.1128/jvi.02158-19

PubMed Abstract | Crossref Full Text | Google Scholar

Snijder, E. J., and Meulenberg, J. J. M. (1998). The molecular biology of arteriviruses. J. Gen. Virol. 79, 961–979. doi: 10.1099/0022-1317-79-5-961

PubMed Abstract | Crossref Full Text | Google Scholar

Snijder, E. J., van Tol, H., Pedersen, K. W., Raamsman, M., and de Vries, A. (1999). Identification of a novel structural protein of arteriviruses. J. Virol. 73, 6335–6345. doi: 10.1128/jvi.73.8.6335-6345.1999

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, R. Q., Guo, Y. Y., Li, X. Y., Li, R., Shi, J., Tan, Z., et al. (2022). PRRSV non-structural proteins orchestrate porcine E3 ubiquitin ligase RNF122 to promote PRRSV proliferation. Viruses 14:424. doi: 10.3390/v14020424

PubMed Abstract | Crossref Full Text | Google Scholar

Tian, K. (2017). NADC30-like porcine reproductive and respiratory syndrome in China. Open Virol. J. 11, 59–65. doi: 10.2174/1874357901711010059

PubMed Abstract | Crossref Full Text | Google Scholar

Tian, K. G., Yu, X. L., Zhao, T. Z., 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., Zhou, Y., and Hao, X. (2007). Identification and molecular epidemiology of thevery virulent porcine reproductive and respiratory syndromeviruses emerged in China. Chin. J. Prev. Vet. Med. 29:5.

Google Scholar

Vu, H. L. X., Kwon, B., Yoon, K. J., Laegreid, W., Pattnaik, A., Osorio, F., et al. (2011). Immune evasion of porcine reproductive and respiratory syndrome virus through glycan shielding involves both glycoprotein 5 as well as glycoprotein 3. J. Virol. 85, 5555–5564. doi: 10.1128/jvi.00189-11

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X. L., Dang, R. Y., Liu, W. K., Yang, Z., Du, E., and Zhang, S. (2014). Antigenic characteristics of glycosylated protein 3 of highly pathogenic porcine reproductive and respiratory syndrome virus. Virus Res. 189, 24–28. doi: 10.1016/j.virusres.2014.04.016

PubMed Abstract | Crossref Full Text | Google Scholar

Wieringa, R., de Vries, A. A. F., Raamsman, M. J. B., and Rottier, P. (2002). Characterization of two new structural glycoproteins, GP3 and GP4, of equine arteritis virus. J. Virol. 76, 10829–10840. doi: 10.1128/jvi.76.21.10829-10840.2002

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, S. S., Liang, W., Wang, X. Y., Chen, H., Fan, J., Song, W., et al. (2020). Epidemiological and genetic characteristics of porcine reproduction and respiratory syndrome virus 2 in mainland China, 2017-2018. Arch. Virol. 165, 1621–1632. doi: 10.1007/s00705-020-04661-z

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, H., Li, C., Li, W. S., 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

Yu, Y., Zhang, Q. Y., Cao, Z., Tang, Y., Xia, D., Wang, G., et al. (2022). Recent advances in porcine reproductive and respiratory syndrome virus NADC30-Like research in China: Molecular characterization, pathogenicity, and control. Front. Microbiol. 12:791313. doi: 10.3389/fmicb.2021.791313

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, H. L., Zhang, W. L., Xiang, L. R., Leng, C., Tian, Z., Tang, Y., et al. (2018a). 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

Zhang, M. Z., Krabben, L., Wang, F. K., and Veit, M. (2018b). Glycoprotein 3 of porcine reproductive and respiratory syndrome virus exhibits an unusual Hairpin-like membrane topology. J. Virol. 92, e660–18. doi: 10.1128/jvi.00660-18

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, H., Luo, Q., Zheng, Y. J., Sha, H., Li, G., Kong, W., et al. (2023). Genetic variability and recombination of the NSP2 gene of PRRSV-2 strains in China from 1996 to 2021. Vet. Sci. 10:325. doi: 10.3390/vetsci10050325

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, H., Sha, H. Y., Qin, L. M., Wang, N., Kong, W., and Huang, L. (2022). Research progress in porcine reproductive and respiratory syndrome virus-host protein interactions. Animals 12:1381. doi: 10.3390/ani12111381

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y. J., Zhang, K., Zheng, H., Liu, C., Jiang, Y., Du, N., et al. (2020). Development of a monoclonal antibody against porcine CD163 SRCR5 domain which partially blocks infection of PRRSV. Front. Vet. Sci. 7:597843. doi: 10.3389/fvets.2020.597843

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, J., Xu, L., Xu, Z. W., Deng, H., Li, F., Sun, X., et al. (2022). Emergence and spread of NADC34-like PRRSV in southwest China. Transbound. Emerg. Dis. 69, E3416–E3424. doi: 10.1111/tbed.14463

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, L., Kang, R., Xie, B., Tian, Y., Yang, X., Yu, J., et al. (2018). Complete genome sequence of a recombinant NADC30-Like Strain, SCnj16, of porcine reproductive and respiratory syndrome virus in southwestern China. Genome Announc. 6:e00004–18. doi: 10.1128/genomeA.00004-18

PubMed Abstract | Crossref Full Text | Google Scholar