- 1Key Laboratory of Animal Diseases and Human Health of Sichuan Province, College of Veterinary Medicine, Sichuan Agricultural University, Chengdu, China
- 2Livestock and Poultry Biological Products Key Laboratory of Sichuan Province, Sichuan Animal Science Academy, Chengdu, China
- 3Animal Breeding and Genetics Key Laboratory of Sichuan Province, Sichuan Animal Science Academy, Chengdu, China
Pseudorabies virus is a major pathogen in the pig industry, causing substantial economic losses. The emergence of pseudorabies virus variant strains in China has led to extensive spread, raising concerns about their potential impact. However, the differences in pathogenicity between the classical strains and the variant strains of genotype II are not well understood. In this study, we isolated three pseudorabies virus strains to evaluate their replication characteristics and to examine the differences in virulence genes among various subgenotypes strains. Additionally, a piglet infection model was utilized to investigate the clinical features of infection, tissue tropism, and the inflammatory responses induced by these strains. Our results showed that the genotype II variant strains (MS, XJ, LS, and CZ) had significantly larger plaque sizes and higher replication capacities than the genotype II classical strain Fa. The animal experiments revealed significant differences in pathogenicity among the pseudorabies virus subgenotype strains, with the variant strains showing higher mortality rates, more severe clinical symptoms, increased nasal virus shedding, and a more robust inflammatory response compared to the genotype II classical strain. There were also notable differences in tissue tropism among the strains. In terms of tissue viral loads, the genotype II variant strains did not exhibit a significant advantage over the genotype I classical strain. Furthermore, our findings indicate that antibodies against the genotype II classical strains have a reduced neutralizing capacity against the genotype II variant strains. On the other hand, antibodies against the genotype II variant strains displayed similar neutralizing abilities against both classical and variant strains. Overall, these findings offer important insights into the distinctions among pseudorabies virus subgenotypes and their implications for the clinical control of pseudorabies virus infections in pig farming.
1 Introduction
Pseudorabies (PR), also known as Aujeszky’s disease (AD), is an acute infectious disease caused by the pseudorabies virus (PRV), which belongs to the order Herpesvirales, family Herperviridae, subfamily alphaherpersvirinae, and genus Varicellovirus (1).1 PRV is transmitted primarily through direct contact with infected suids or by consuming raw offal from infected pigs, whether wild or domestic, or pig carcasses (2). PRV strains are currently classified into two genotypes: genotype I and genotype II. Genotype I is commonly found in Argentina, Spain, Italy, the United States, and China (3–6), whereas genotype II PRV isolates are predominantly present in China. Genotype II is further categorized into two subgenotypes: the classic strains and variant strains (7). Infection with PRV leads to acute and fatal infections in piglets, abortion in sows and growth retardation or weight loss in growing pigs (8, 9). While PRV can infect various domestic and wild animals, pigs are its only natural host (10, 11). In recent years, several cases have confirmed that PRV can also cause encephalitis in humans, leading to symptoms such as fever, sweating, weakness, seizures, and even death (12). The recovery prospects for patients are extremely poor, posing a severe threat to public health (13).
Before 2011, China had effectively controlled PRV through vaccination using the Bartha-Hungary-1961 strain-based vaccine (Bartha-K61). Some regions, like parts of the United States, Canada, and Europe, even achieved complete eradication of PRV (14). However, from October 2011, a rapid spread of the PRV epidemic occurred in China, causing substantial economic losses to the pig industry. The outbreak was caused by the genotype II PRV variant strain, which differs from classical strains due to deletions in amino acid positions 75–77 of the glycoprotein B (gB) (15). It is worth noting that the Bartha-K61 vaccine has been found ineffective against the genotype II variant strain (16). Since its emergence, the PRV variant strain has caused significant economic losses to the pig industry in China (17).
Variant strains of PRV, in comparison to the classical PRV strains, display enhanced virulence, increased infectivity, faster transmission, and present with more severe clinical symptoms. Recent studies have even indicated increased virulence in recombinant strains (18). However, the relationship between PRV genetic features and virulence is still not fully understood. In this study, we isolated variant strains of LS, MS, and CZ. The early classical strain Fa and the variant strain XJ were used as reference strains to analyze the genetic evolution and virulence changes of the current popular variant strains, focusing on the main glycoprotein genes gB, glycoprotein C (gC), glycoprotein D (gD), and glycoprotein E (gE) and PRV pathogenicity to piglets.
2 Materials and methods
2.1 Virus and cells
Baby hamster kidney (BHK-21) cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, Waltham, MA, United States), supplemented with 10% fetal bovine serum (FBS, Gibco, Waltham, MA, United States), and maintained in an incubator with 5% CO2. The PRV XJ strain was isolated and preserved in our laboratory, as previously described (19). The PRV Fa strain was also preserved in our laboratory (20). Pig anti-Fa hyperimmune serum and pig anti-XJ strain hyperimmune serum were prepared and stored in our laboratory.
2.2 Isolation and identification of PRV
The polymerase chain reaction (PCR) method was utilized for the detection of PRV in brain, liver, and lung tissues collected from pigs from farms located in different regions of Sichuan Province. The PCR was conducted using the following primers: gE forward primer (5′-ATCTGGACGTTCCTGCCC-3′) and gE reverse primer (5′-GTAGATGCAGGGCTCGTACA-3′) (21). The PRV-positive samples were minced and homogenized in 1 mL of DMEM. Subsequently, the PRV-positive samples underwent three freeze–thaw cycles and were centrifuged at 5000 rpm for 10 min at 4°C to obtain the supernatant. Utilize the PCR method to confirm the presence of PRV in the supernatant (21). The supernatant from PRV-positive samples was then filtered through a 0.22 μm membrane. The filtered supernatant was used to inoculate confluent BHK-21 cells monolayers. Following a 1-h incubation at 37°C, the supernatant was aspirated, and the cells were further cultured in DMEM with 2% FBS for 48 h. Cytopathic effects (CPE) were monitored daily. Infected cells were fixed using 4% paraformaldehyde and then characterized by immunofluorescence assay (IFA) with a monoclonal antibody against the PRV gB protein sourced from Beijing TianTech Biotechnology (19).
2.3 Transmission electron microscopic observation of virions
At 48 h post-infection (hpi), the cells were harvested and fixed using a 3% glutaraldehyde fixative. Following fixation, the cells were processed for transmission electron microscopy by Chengdu Lilai Biotechnology Co., Ltd., enabling the capture of virions images (21).
2.4 One-step growth curve
As described previously (22), one-step growth kinetics assays were conducted. BHK-21 cells were plated in 96-well plates at a density sufficient to form confluent monolayer. The cells were then infected with the PRV strains MS, CZ, XJ, LS, or Fa at a multiplicity of infection (MOI) of 0.1. Post-infection, cells and supernatants were harvested at intervals of 6, 12, 24, 36, 48, and 72 hpi. Virus titers were calculated using the 50% tissue culture infectious dose (TCID50) assay, and these data were utilized to generated the one-step growth curves for each PRV strain.
To evaluate plaque sizes of PRV strains from different subgenotype, BHK-21 cells were infected with PRV at a titer of 100 TCID50, employing plaque assays according to the previous study protocol (21).
2.5 Phylogenetic analysis
Total DNA from the XJ, MS, CZ, LS and Fa strain was extracted using a universal genomic DNA kit (Cofitt, Jiangsu, China). PCR amplification was conducted using PrimeSTAR HS DNA Polymerase and GC Buffer (TaKaRa Bio, Shiga, Japan), with primers for each gene listed in Table 1. Purified positive PCR products were cloned into the pUCm-T vector and submitted to Shanghai Bioengineering Co., Ltd. for sequencing. Full-length sequences of the gB, gC, gD, and gE genes from the MS, CZ, XJ, LS, and Fa strains were aligned using DNAStar version 15.7 software. Phylogenetic trees were constructed based on the major virulence genes (gB, gC, gD, and gE) using Molecular Evolutionary Genetics Analysis (MEGA; version 5.1) software, employing the maximum likelihood method and 1,000 bootstrap replicates for validation.
2.6 Animal experiments
Thirty two-week-old female piglets were screened and confirmed to be negative for PRV gB-specific antibodies via an enzyme-linked immunosorbent assay (ELISA) kit (ID. vet, France) and negative for PRV DNA using PCR with gE-specific primers. They were then randomly assigned into 6 groups, each group consisting of 5 piglets. Piglets in groups 1 to 5 (n = 5 each) were intranasally infected with 2 mL of MS, CZ, Fa, XJ, or LS strains, at a titer of 107 TCID50/mL, respectively. Group 6 served as the uninfected control and was intranasally inoculated with 2 mL of DMEM. Rectal temperatures and clinical signs were recorded daily following inoculation. Daily nasal swabs collections were conducted, and viral shedding was assessed using qRT-PCR. Clinical symptoms were scored according to previous descriptions as follows (23): (1) temperature elevation above 40°C; (2) temperature elevation above 40°C with respiratory distress; (3) ataxia; (4) convulsions; and (5) moribund or dead.
Following the challenge, piglets exhibiting distress symptoms, including respiratory difficulties, rapid or irregular breathing, inability to stand, or anorexia, were humanely euthanized using a lethal dose of pentobarbital sodium (40 mg/kg). On day 14, all remaining piglets, including the survivors and control group piglets, were euthanized via an intrajugular injection of pentobarbital sodium (40 mg/kg). After euthanasia, tissues including the cerebellum, brainstem, tonsils, inguinal lymph nodes, mesenteric lymph nodes, lungs, liver, brain, intestines, spleen, and kidneys were collected for further analysis.
2.7 Real-time fluorescence quantitative PCR (qRT-PCR)
Fresh tissues samples and nasal swabs were collected for viral loads analysis. Genomic DNA was extracted from the samples using a universal genomic DNA kit (Cofitt, Jiangsu, China). We performed qRT-PCR using TB Green® Premix Ex Taq™ (Tli RNaseH Plus) (RR420Q, Takara, Japan) to quantify the PRV DNA copies, specifically for the gE gene with the following primers: Forward (5’-CTTCCACTCGCAGCTCTTCT-3′), Reverse (5’-TAGATGCAGGGCTCGTACAC-3′). The gene copies per gram of tissue sample were expressed as log10 copies.
2.8 Cross-neutralization assay
Neutralizing titers of antisera obtained from Fa vaccine and rPRV-gE/gI/TK− strain-inoculated animals were determined using methods previously described (24). The sera were first heat-inactivated at 56°C for 30 min prior to serial dilution. Subsequently, the diluted sera were incubated with an equal volume of XJ, LS, Fa, CZ, or MS, each at a titer of 200 TCID50/0.1 mL, at 37°C for 1 h. Next, 100 μL of BHK-21 cells, prepared at a concentration of 106 cells/mL in DMEM with 2% FBS, were added to each well of the 96-well plates. The plates were incubated at 37°C with 5% CO2, and CPE were monitored daily. Neutralizing antibody titers were calculated using the Reed-Muench method, defining the serum dilution preventing 50% of cells from CPE as the neutralization titer.
2.9 Quantification of inflammatory markers
Concentration of IL-1β, IL-6, IL-8, and TNF-α were measured in lung, brain, spleen, and intestine tissues using Elisa kits from Thermo Fisher, according to the manufacturer’s protocol.
2.10 Statistical analysis
All data were analyzed using GraphPad 7.04 software with one-way analysis of variance (ANOVA). All results are presented as the mean ± standard deviation. A significance level of p < 0.05 was considered statistically significant.
3 Results
3.1 In vitro proliferation characteristics of PRV
Following PRV infection, a distinct green fluorescence signal, indicative of the PRV gB protein, was detected within the CPE regions (Figure 1). Transmission electron microscopy showed virus particles, devoid of a capsule structure, within the nucleus, measuring approximately 100 nm in diameter (Figure 2A). Furthermore, mature virus particles, enclosed by a capsule, were identified in the cytoplasm, measuring about 160 nm in diameter. These observations are consistent with the known surface characteristics of PRV. The isolated strains, MS, CZ, and LS, were named according to their geographical origins. The one-step growth curve demonstrated that the MS, CZ, LS, XJ, and Fa strains exhibited similar growth characteristics in BHK-21 cells, peaking at 36 hpi with virus titers of 108.18/mL for MS, 108.33/mL for CZ, 108.08/mL for XJ, 107.81/mL for LS, and 107.83/mL for Fa, as illustrated in Figure 2B. Plaque assay results showed that the LS strain had a significantly higher plaque generation ability compared to the MS, XJ, CZ, and Fa strains (Figure 2C). The MS, XJ, and CZ strains exhibited similar plaque generation abilities, each significantly exceeding that of the Fa strain (Figure 2C). This indicates that the infectivity of the LS, MS, XJ, and CZ strains was significantly stronger than that of the Fa strain.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438354/fvets-11-1438354-HTML-r1/image_m/fvets-11-1438354-g001.jpg)
Figure 1. IFA to identify PRV strains. 4′,6-diamidino-2-phenylindole (DAPI): DAPI was utilized for staining cell nuclei. Fluorescein Isothiocyanate (FITC): following binding of the anti-gB protein mouse monoclonal antibody to the PRV gB protein, FITC-labeled goat anti-mouse immunoglobulin G (IgG) is utilized to bind to the complex, resulting in the emission of green fluorescence. Merge: combine DAPI staining pictures and FITC pictures. Bright: capture images of cells under bright field.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438354/fvets-11-1438354-HTML-r1/image_m/fvets-11-1438354-g002.jpg)
Figure 2. Biological characteristics of different strains of PRV in vitro. (A) Transmission electron microscopy images of PRV. The red arrow points to the mature virion, and the black arrow points to the noncapsular nucleocapsid. (B) One-step growth curves of the LS, XJ, MS, CZ, and Fa strains in BHK-21 cells. In the growth curves, the x-axis represents the infection time, and the y-axis represents the virus titer. (C) BHK-21 cells were infected with the LS, XJ, MS, CZ and Fa strains. The plaque size of each virus was assessed by performing the plaque assay. ** p < 0.01; *** p < 0.001; **** p < 0.0001; ns, not significant.
3.2 Phylogenetic analysis
The gB, gC, gD, and gE gene sequences of the MS, CZ, and LS strains were uploaded to NCBI (accession numbers: OR338844-OR338846 and OR334583-OR334586).
Analysis of the gB and gC gene sequences indicated that the MS, CZ, and LS strains are significantly distinct from the genotype II classical strains in China and genotype I strains from other countries, clustering more closely with the variant strain XJ. Within this cluster, the MS and CZ strains showed a close relationship with HN1201, whereas the LS strain showed greater similarity to JS2012, all of which are distinct from the classical strains Fa and Ea, as illustrated in Figures 3A,B. Analysis of the gD sequences revealed that the MS strain shared similarity with HN1201, the CZ strain had a close relationship with HeNLH/2017, and the LS strain was closely associated with SD2018 (Figure 3C). Subsequentanalysis of gE gene sequences indicated that the MS, CZ, and LS formed a distinct cluster within the genotype II variant strains (Figure 3D). Altogether, the phylogenetic analysis of the gB, gC, gD, and gE genes suggests that the MS, CZ, and LS strains are classified as genotype II variant strains.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438354/fvets-11-1438354-HTML-r1/image_m/fvets-11-1438354-g003.jpg)
Figure 3. Phylogenetic analysis of the gB, gC, gD and gE genes of PRV. Phylogenetic trees based on the gB (A), gC (B), gD (C) and gE (D) genes were constructed using the maximum likelihood method. The PRV strains used in this study are indicated by a red dot. The PRV genotype II classical strains are indicated by green, and the PRV genotype II variant strains are indicated by blue.
3.3 Clinical symptoms and viral excretion
Piglets infected with the variant viruses LS, XJ, CZ, and MS developed severe clinical symptoms such as muscle spasms, depression, decreased appetite, and reduced water intake, resulting in rapid death. In contrast, piglets in the Fa strain-infected group showed mild shivering, depression, and decreased appetite. From the 9th day post-infection, the symptoms of the surviving piglets stabilized and did not progress (Figure 4A). Within 2 days post-infection (dpi), piglets in the variant virus-infected groups developed a fever exceeding 40°C that persisted until their demise. By the 10th day post-infection, surviving piglets in the Fa-infected group had returned to a normal body temperatures range of 38.5°C to 39.5°C (Figure 4B). Piglets infected with the variant viruses LS, XJ, CZ, and MS began to die 72 h post-infection, and all had succumbed within 8 days. In the Fa-infected group, 2 out of the 5 piglets survived (Figure 4C). Daily nasal swab collections were conducted to assess viral shedding, followed by analysis with qRT-PCR. PRV was detected in the nasal secretions of all infected piglets, with peak viral shedding observed on the 6th day post-infection. Compared to other PRV-infected groups, the MS-infected group exhibited higher viral copy numbers from the initial day infection (Supplementary Figure S2). In contrast, the Fa strain showed a declining trend in viral copy numbers following the 6th day post-infection (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438354/fvets-11-1438354-HTML-r1/image_m/fvets-11-1438354-g004.jpg)
Figure 4. Clinical symptoms and shedding virus of newborn piglets challenged with the PRV strains. Newborn piglets were infected with different PRV strains or infected with DMEM as a mock group. The survival rate, clinical symptoms and shedding virus of each group were observed and recorded daily. (A) Clinical symptom scores of piglets after challenge. Clinical symptoms, including body temperature, respiratory distress, ataxia, convulsion, and moribund or dead, were scored as previously described (25). (B) The rectal temperature of all piglets after challenge with PRV or DMEM. The newborn piglets were challenged with different PRV strains or DMEM, and then the rectal temperature was recorded daily. Clinical fever was defined as a rectal temperature ≥ 40°C. (C) Survival curves of newborn piglets infected with different PRV strains were generated. (D) The shedding virus of newborn piglets after challenge. Nasal swabs were collected daily, and the PRV copy number was determined by qRT-PCR. Copies/μg represents the number of DNA copies of the PRV genome per 1 μg of total DNA.
3.4 Pathogenicity analysis in piglets
Quantitative PCR results showed no significant differences in viral loads in the cerebellum and ileum of piglets infected with any of the five different PRV strains. However, strains XJ and LS showed significantly higher viral loads within the brain. Additionally, strains XJ and Fa strains exhibited significantly higher viral loads in the brainstem compared to the other strains (Figure 5A). The CZ strain was observed to have a significant increase in viral loads in lymphoid organs, including the spleen, inguinal lymph nodes, and mesenteric lymph nodes, surpassing those of other strains. In digestive organs, such as the liver, ileum, and duodenum, the MS and Fa strains showed significantly higher viral loads compared to the other strains. In lung tissue, the MS and CZ strains showed significantly higher viral loads compared to the other strains, showing variations in tissue tropism among different PRV strains (Figure 5B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438354/fvets-11-1438354-HTML-r1/image_m/fvets-11-1438354-g005.jpg)
Figure 5. Viral load and histopathological analysis in PRV-challenged piglets. (A,B) Viral load of different PRV strains in various tissues of PRV-challenged piglets. Brain, cerebellum, brainstem, ileum, jejunum, duodenum, liver, spleen, kidney, lung, inguinal lymph nodes, mesenteric lymph nodes and tonsil from PRV-challenged piglets were tested to determine the copy number of PRV by qRT-PCR. copies/g represents the number of DNA copies per microgram of tissue. The data in (A,B) represent the mean ± SD. Experiments were performed independently with at least three biological replicates. One-way ANOVA was used for analysis. * p < 0.5; ** p < 0.01; *** p < 0.001; **** p < 0.0001; ns, not significant.
3.5 Cross-neutralization responses between variant strains and classical strains
Subsequently, we performed cross-neutralization assays to evaluate the serological cross-reactivity between different PRV strains. Antisera derived from vaccination with the classical Fa strain showed a significantly reduced neutralizing capacity against the variant strains, including XJ, as compared to its neutralizing capacity against the Fa strain itself (Figure 6A). In contrast, antisera from the XJ strain displayed neutralizing titers comparable to those against the Fa strain for all variant strains assessed (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1438354/fvets-11-1438354-HTML-r1/image_m/fvets-11-1438354-g006.jpg)
Figure 6. Cross-neutralization tests and inflammatory factor levels in the lungs, liver, brain and intestine. (A) Neutralizing titers of Fa strain antisera against the 5 PRV strains. The significance of differences between the Fa group and other groups was analyzed with one-way ANOVA. **** p < 0.0001; ns, not significant. (B) Neutralizing titers of XJ strain antisera against the 5 PRV strains. The inflammatory factor levels of the tissues in all groups, including IL-1β (C), IL-6 (D), IL-8 (E) and TNF-α (F). The significance of differences between the mock group and other groups was analyzed with one-way ANOVA. ** p < 0.01; *** p < 0.001; ****, p < 0.0001. The data in this figure represent the mean ± SD. Experiments were performed independently with at least three biological replicates.
3.6 Induction of high levels of inflammatory cytokines in tissues after PRV infection
Inflammatory cytokine assays demonstrate that compared to the mock group, IL-1β, IL-6, IL-8, and TNF-α were significantly upregulated in the lungs, spleen, intestines, and brain post-infection. Specifically, the LS strain induced significantly higher levels of IL-1β in the lungs, and the XJ strain induced significantly higher levels of IL-1β in the brain compared to other PRV strains (Figure 6C). The MS strain induced significantly higher levels of IL-6 in the lungs, and the CZ strain induced significantly higher levels of IL-6 in the spleen compared to other strains. In both the brain and intestines, elevated IL-6 levels were observed with the XJ strain, surpassing those induced by the other strains (Figure 6D). The CZ strain induced significantly higher levels of IL-8 in the lungs, and both the Fa and CZ strains induced significantly higher levels of IL-8 in the intestines compared to other strains (Figure 6E). Furthermore, the LS strain induced significantly higher levels of TNF-α in the lungs compared to other strains (Figure 6F).
4 Discussion
PRV, the causative agent of Aujeszky’s disease, has caused significant economic losses in the swine industry since the emergence of variant PRV strains in China in 2011. The three isolated PRV variant strains exhibited indistinguishable in vitro proliferation characteristics in BHK-21 cells, with all strains reaching peak viral titers within 36 h, reflecting their comparable replication abilities in vitro. But plaque assay results showed that the LS, MS, XJ, and CZ strains exhibited significantly higher plaque generation ability compared to the Fa strain, which may be due to the weaker cell-to-cell spread ability of the Fa strain. Previous studies have found that the deletion of amino acids 124–495 in the HSV-2 gE protein severely inhibits the cell-to-cell spread efficiency of HSV-2 (26). The gE protein, as a membrane protein, is closely related to the cell-to-cell spread ability of PRV, and is also involved in the transport and localization of viral particles, affecting viral evasion and dissemination (27). Distinct from the Fa strain, the MS, CZ, XJ, and LS strains have mutations at aa54 (G to D), aa403 (P to A), aa448 (V to L), aa492 (G to D), and aa518 (S to P) (Supplementary Figure S1). These specific amino acid substitutions could potentially influence the efficiency of cell-to-cell spread in the variant strains. To validate these findings, further studies involving the construction of mutant and revertant strains are required. This will be a focus of our laboratory’s future research efforts.
In assessing the neutralizing abilities of sera, we observed that animals immunized with the XJ strain exhibited similar neutralizing capacities against all variant strains. This suggests that vaccines based on variant strains can effectively prevent infection by prevalent strains. This neutralizing ability may be attributed to the consistency of neutralizing epitopes on the gB protein between variant strains and the classical strain (28). In contrast, neutralizing antibodies derived from the classical strain showed limited effectiveness against variant strains, possibly because the variant strains have evolved immune escape mechanisms, including the evasion of neutralizing antibodies induced by the classical strain. Ren et al. (29) found that gC and gD proteins contribute to evading neutralizing antibodies induced by the Bartha-K61 vaccine. Deng et al. (30) discovered differences in the usage of immunoglobulin heavy chain variable region genes and the combination of immunoglobulin heavy chain variable region and IGHJ genes between the Bartha-K61 and XJ strain infection groups. Additionally, PRV has developed mechanisms to evade host innate immunity, such as UL13 suppressing the expression of RIG-I and MDA5 by inhibiting the activation of the transcription factor NF-κB (31), and gE counteracting cGAS/STING-mediated IFN production (32). These mechanisms may differ between classical and variant strains. In conclusion, we hypothesize that PRV variant strains may have evolved multiple mechanisms to escape neutralization by antibodies against classical strains, but further research is needed to confirm this.
The different strains of PRV induced varying levels of cytokines in specific tissues, reflecting their unique pathogenic profiles. For instance, the LS strain’s induction of high IL-1β and TNF-α levels in the lungs suggests it may cause more severe pulmonary pathology compared to other strains. In addition, IL-1β is crucial for the activation of inflammatory responses, and TNF-α, a key mediator of inflammation, should exhibit trends consistent with IL-1β. This is consistent with our findings. The XJ strain’s significant induction of IL-1β and IL-6 in the brain implies a propensity for causing severe neurological damage, aligning with observed clinical symptoms like muscle spasms and depression. Several studies have indicated that encephalitis is a major cause of mortality in infected animals (33–36). However, in the surviving piglets of the Fa group, elevated levels of inflammatory cytokines were also detected in the lungs, brain, intestines, and spleen compared to the control group. Therefore, the mechanisms underlying acute mortality in piglets caused by PRV remain to be further elucidated.
Our study has several limitations. First, the limited number of PRV strains, all from specific geographical locations, may not fully represent the genetic diversity of PRV. For instance, we cannot determine whether there are significant differences in the onset, course, and pathogenicity of diseases caused by early variant strains compared to currently prevalent variant strains. Second, our in vivo model was restricted to piglets, which may not comprehensively capture the disease dynamics in other pig populations. Future research should address these limitations by including a more diverse set of PRV strains from different regions. Additionally, exploring other animal models and conducting longitudinal studies will provide a more comprehensive understanding of PRV pathogenesis.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The animal studies were approved by the Ethics Committee of Sichuan Agricultural University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
LX: Conceptualization, Data curation, Formal analysis, Software, Writing – original draft, Writing – review & editing. QT: Formal analysis, Software, Writing – review & editing. TX: Formal analysis, Software, Writing – review & editing. YY: Supervision, Writing – review & editing. YaZ: Supervision, Writing – review & editing. ZL: Formal analysis, Writing – review & editing. YuZ: Software, Writing – review & editing. LZ: Conceptualization, Funding acquisition, Project administration, Writing – review & editing. ZX: Conceptualization, Resources, Validation, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Chongqing Municipal Technology Innovation and Application Development Project (No. cstc2021jscx-dxwt BX0007), the Key K&D Program of Sichuan Science and Technology Plan (No. 2022YFN0007), the Porcine Major Science and Technology Project of Sichuan Science and Technology Plan (No. 2021ZDZX0010-3), the Sichuan Science and Technology Program Projects (Key R&D Projects) (No. 2023YFN0021) and the Agricultural Industry Technology System of Sichuan Provincial Department of Agriculture (No. CARS-SVDIP).
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.2024.1438354/full#supplementary-material
Footnotes
References
1. Xu, J, Cheng, X, Liu, Y, Fu, X, Tong, W, Zheng, H, et al. Pseudorabies virus UL16 protein influences the inhibition of LRPPRC for the viral proliferation. Vet Microbiol. (2022) 265:109327. doi: 10.1016/j.vetmic.2021.109327
2. Laval, K, and Enquist, LW. The neuropathic itch caused by pseudorabies virus. Pathogens. (2020) 9:254. doi: 10.3390/pathogens9040254
3. Moreno, A, Musto, C, Gobbi, M, Maioli, G, Menchetti, M, Trogu, T, et al. Detection and molecular analysis of pseudorabies virus from free-ranging Italian wolves (Canis lupus italicus) in Italy - a case report. BMC Vet Res. (2024) 20:9. doi: 10.1186/s12917-023-03857-0
4. Moreno, A, Sozzi, E, Grilli, G, Gibelli, LR, Gelmetti, D, Lelli, D, et al. Detection and molecular analysis of pseudorabies virus strains isolated from dogs and a wild boar in Italy. Vet Microbiol. (2015) 177:359–65. doi: 10.1016/j.vetmic.2015.04.001
5. Serena, MS, Metz, GE, Lozada, MI, Aspitia, CG, Nicolino, EH, Pidone, CL, et al. First isolation and molecular characterization of Suid herpesvirus type 1 from a domestic dog in Argentina. Open Vet J. (2018) 8:131–9. doi: 10.4314/ovj.v8i2.3
6. Verin, R, Varuzza, P, Mazzei, M, and Poli, A. Serologic, molecular, and pathologic survey of pseudorabies virus infection in hunted wild boars (Sus scrofa) in Italy. J Wildl Dis. (2014) 50:559–65. doi: 10.7589/2013-01-004
7. Lin, Y, Tan, L, Wang, C, He, S, Fang, L, Wang, Z, et al. Serological investigation and genetic characteristics of pseudorabies virus in Hunan Province of China from 2016 to 2020. Front Vet Sci. (2021) 8:762326. doi: 10.3389/fvets.2021.762326
8. An, T-Q, Peng, J-M, Tian, Z-J, Zhao, H-Y, Li, N, Liu, Y-M, et al. Pseudorabies virus variant in Bartha-K61–vaccinated pigs, China, 2012. Emerg Infect Dis J. (2013) 19:1749–55. doi: 10.3201/eid1911.130177
9. Gu, Z, Dong, J, Wang, J, Hou, C, Sun, H, Yang, W, et al. A novel inactivated gE/gI deleted pseudorabies virus (PRV) vaccine completely protects pigs from an emerged variant PRV challenge. Virus Res. (2015) 195:57–63. doi: 10.1016/j.virusres.2014.09.003
10. Liu, A, Xue, T, Zhao, X, Zou, J, Pu, H, Hu, X, et al. Pseudorabies virus associations in wild animals: review of potential reservoirs for cross-host transmission. Viruses. (2022) 14:2254. doi: 10.3390/v14102254
11. Verpoest, S, Cay, AB, and De Regge, N. Molecular characterization of Belgian pseudorabies virus isolates from domestic swine and wild boar. Vet Microbiol. (2014) 172:72–7. doi: 10.1016/j.vetmic.2014.05.001
12. Ou, J, Cai, S, Zheng, F, Lu, G, and Zhang, G. Human pseudorabies virus infection: a new threat in China. J Infect. (2020) 80:578–606. doi: 10.1016/j.jinf.2019.12.018
13. Fan, S, Yuan, H, Liu, L, Li, H, Wang, S, Zhao, W, et al. Pseudorabies virus encephalitis in humans: a case series study. J Neurovirol. (2020) 26:556–64. doi: 10.1007/s13365-020-00855-y
14. Sun, Y, Luo, Y, Wang, CH, Yuan, J, Li, N, Song, K, et al. Control of swine pseudorabies in China: opportunities and limitations. Vet Microbiol. (2016) 183:119–24. doi: 10.1016/j.vetmic.2015.12.008
15. Sun, Y, Liang, W, Liu, Q, Zhao, T, Zhu, H, Hua, L, et al. Epidemiological and genetic characteristics of swine pseudorabies virus in mainland China between 2012 and 2017. PeerJ. (2018) 6:e5785. doi: 10.7717/peerj.5785
16. Hu, D, Zhang, Z, Lv, L, Xiao, Y, Qu, Y, Ma, H, et al. Outbreak of variant pseudorabies virus in Bartha-K61-vaccinated piglets in Central Shandong Province, China. J Vet Diagn Invest. (2015) 27:600–5. doi: 10.1177/1040638715593599
17. Tan, L, Yao, J, Yang, Y, Luo, W, Yuan, X, Yang, L, et al. Current status and challenge of pseudorabies virus infection in China. Virol Sin. (2021) 36:588–607. doi: 10.1007/s12250-020-00340-0
18. Qin, Y, Qin, S, Huang, X, Xu, L, Ouyang, K, Chen, Y, et al. Isolation and identification of two novel pseudorabies viruses with natural recombination or TK gene deletion in China. Vet Microbiol. (2023) 280:109703. doi: 10.1016/j.vetmic.2023.109703
19. Yin, Y, Xu, Z, Liu, X, Li, P, Yang, F, Zhao, J, et al. A live gI/gE-deleted pseudorabies virus (PRV) protects weaned piglets against lethal variant PRV challenge. Virus Genes. (2017) 53:565–72. doi: 10.1007/s11262-017-1454-y
20. Zhu, L, Yi, Y, Xu, Z, Cheng, L, Tang, S, and Guo, W. Growth, physicochemical properties, and morphogenesis of Chinese wild-type PRV fa and its gene-deleted mutant strain PRV SA215. Virol J. (2011) 8:272. doi: 10.1186/1743-422X-8-272
21. Tao, Q, Zhu, L, Xu, L, Yang, Y, Zhang, Y, Liu, Z, et al. The construction and immunogenicity analyses of a recombinant pseudorabies virus with Senecavirus a VP2 protein Coexpression. Microbiol Spectr. (2023) 11:e0522922. doi: 10.1128/spectrum.05229-22
22. Xu, L, Wei, JF, Zhao, J, Xu, SY, Lee, FQ, Nie, MC, et al. The immunity protection of central nervous system induced by pseudorabies virus DelgI/gE/TK in mice. Front Microbiol. (2022) 13:862907. doi: 10.3389/fmicb.2022.862907
23. Hu, RM, Zhou, Q, Song, WB, Sun, EC, Zhang, MM, He, QG, et al. Novel pseudorabies virus variant with defects in TK, gE and gI protects growing pigs against lethal challenge. Vaccine. (2015) 33:5733–40. doi: 10.1016/j.vaccine.2015.09.066
24. Zhao, J, Zhu, L, Xu, L, Li, F, Deng, H, Huang, Y, et al. The construction and immunogenicity analyses of recombinant pseudorabies virus with NADC30-like porcine reproductive and respiratory syndrome virus-like particles co-expression. Front Microbiol. (2022) 13:846079. doi: 10.3389/fmicb.2022.846079
25. Ren, X, Cao, N, Tian, L, Liu, W, Zhu, H, Rong, Z, et al. A self-assembled nanoparticle vaccine based on pseudorabies virus glycoprotein D induces potent protective immunity against pseudorabies virus infection. Vet Microbiol. (2023) 284:109799. doi: 10.1016/j.vetmic.2023.109799
26. Wang, F, Zumbrun, EE, Huang, J, Si, H, Makaroun, L, and Friedman, HM. Herpes simplex virus type 2 glycoprotein E is required for efficient virus spread from epithelial cells to neurons and for targeting viral proteins from the neuron cell body into axons. Virology. (2010) 405:269–79. doi: 10.1016/j.virol.2010.06.006
27. Ning, Y, Huang, Y, Wang, M, Cheng, A, Yang, Q, Wu, Y, et al. Alphaherpesvirus glycoprotein E: a review of its interactions with other proteins of the virus and its application in vaccinology. Front Microbiol. (2022) 13:970545. doi: 10.3389/fmicb.2022.970545
28. Wang, Y, Qiao, S, Li, X, Xie, W, Guo, J, Li, Q, et al. Molecular epidemiology of outbreak-associated pseudorabies virus (PRV) strains in Central China. Virus Genes. (2015) 50:401–9. doi: 10.1007/s11262-015-1190-0
29. Ren, J, Wang, H, Zhou, L, Ge, X, Guo, X, Han, J, et al. Glycoproteins C and D of PRV strain HB1201 contribute individually to the escape from Bartha-K61 vaccine-induced immunity. Front Microbiol. (2020) 11:323. doi: 10.3389/fmicb.2020.00323
30. Deng, L, Yang, F, Xu, Z, Li, F, Zhao, J, Deng, H, et al. Characterization of B cell receptor H-CDR3 repertoire of spleen in PRV-infected mice. BMC Vet Res. (2022) 18:228. doi: 10.1186/s12917-022-03340-2
31. Zhao, N, Wang, F, Kong, Z, and Shang, Y. Pseudorabies virus tegument protein UL13 suppresses RLR-mediated antiviral innate immunity through regulating receptor transcription. Viruses. (2022) 14:1465. doi: 10.3390/v14071465
32. Lu, M, Qiu, S, Zhang, L, Sun, Y, Bao, E, and Lv, Y. Pseudorabies virus glycoprotein gE suppresses interferon-β production via CREB-binding protein degradation. Virus Res. (2021) 291:198220. doi: 10.1016/j.virusres.2020.198220
33. Engelhardt, S, Schneider, S, Buder, A, Aupperle-Lellbach, H, and Pfeil, I. MRI in a dog with confirmed pseudorabies infection. Tierarztl Prax Ausg K Klientiere Heimtiere. (2019) 47:272–81. doi: 10.1055/a-0948-8760
34. Yang, H, Han, H, Wang, H, Cui, Y, Liu, H, and Ding, S. A case of human viral encephalitis caused by pseudorabies virus infection in China. Front Neurol. (2019) 10:534. doi: 10.3389/fneur.2019.00534
35. Liu, Q, Wang, X, Xie, C, Ding, S, Yang, H, Guo, S, et al. A novel human acute encephalitis caused by pseudorabies virus variant strain. Clin Infect Dis. (2020) 73:e3690–700. doi: 10.1093/cid/ciaa987
Keywords: pseudorabies, pathogenicity, variant strains, isolation, inflammatory factor
Citation: Xu L, Tao Q, Xu T, Yang Y, Zhang Y, Liu Z, Zhou Y, Zhu L and Xu Z (2024) Pathogenicity characteristics of different subgenotype pseudorabies virus in newborn piglets. Front. Vet. Sci. 11:1438354. doi: 10.3389/fvets.2024.1438354
Edited by:
Jianzhu Liu, Shandong Agricultural University, ChinaReviewed by:
Giuseppe Marruchella, University of Teramo, ItalyYan-Dong Tang, Chinese Academy of Agricultural Sciences, China
Copyright © 2024 Xu, Tao, Xu, Yang, Zhang, Liu, Zhou, Zhu and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ling Zhu, abtczl72@126.com; Zhiwen Xu, abtcxzw@126.com
†These authors have contributed equally to this work