Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Vet. Sci., 08 October 2025

Sec. Veterinary Epidemiology and Economics

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1667063

Molecular characteristics and potential antigenic epitope analysis of porcine epidemic diarrhea virus in China from 2022 to 2025

Zunbao Wang,&#x;Zunbao Wang1,2Kai Yang&#x;Kai Yang1Mingfang BiMingfang Bi1Kaijie LiKaijie Li1Wen WangWen Wang2Yan SongYan Song2Xiaomei PanXiaomei Pan2Tianzeng Li
Tianzeng Li2*Xiaobing Mo
Xiaobing Mo1*
  • 1State Key Laboratory for Diagnosis and Treatment of Severe Zoonotic Infectious Diseases, Key Laboratory for Zoonosis Research of the Ministry of Education, College of Veterinary Medicine, Institute of Zoonosis, Jilin University, Changchun, China
  • 2Tecon Bio-Pharmaceuticals Co., Ltd., Urumqi, China

Porcine epidemic diarrhea virus (PEDV) causes varying degrees of diarrhea, vomiting, dehydration, and emaciation in pigs, and severe cases can result in the death of suckling piglets. PEDV is one of the significant viruses affecting the global swine industry. In this study, we investigated the prevalence of PEDV infection in 2,346 pig samples collected from 20 provinces in China, with an overall positive rate of 43.3% (1,015/2,346). Phylogenetic analysis of 15 newly sequenced PEDV strains revealed that 1 strain belonged to the G1c (S-INDEL) genotype, 4 strains to G2b, and the remaining 10 strains to G2c, indicating that G2c is currently the dominant subtype. Compared with the classical CV777 strain, the nucleotide and amino acid sequence identities of the S gene of these 15 PEDV strains ranged from 92.71 to 94.83% and 92.89 to 94.99%, respectively. Amino acid alignment identified mutations of varying degrees within key neutralizing epitopes, including COE, SS2, SS6, and 2C10. Using the server, a new potential N-glycosylation site at position 302 (NKTI) in the S protein of the PEDV/XinJiang/2 strain was identified. Furthermore, linear and discontinuous antigenic epitopes of resolved PEDV proteins were predicted using ElliPro, revealing conserved potential antigenic regions at amino acid 31–54, 90–103, 1,170–1,177, 1,179–1,197, 1,201–1,211, 1,221–1,236, and 1,243–1,254. These findings provide new epidemiological data on PEDV circulating in China and offer valuable insights for the development of future vaccines and molecular diagnostics based on these antigenic epitopes, thereby contributing to improved PED control strategies.

1 Introduction

Porcine epidemic diarrhea virus (PEDV) is a single-stranded positive-sense RNA virus belonging to the genus Alphacoronavirus within the family Coronaviridae, which causes porcine epidemic diarrhea (1). Once inside the host, the virus infects the small intestinal mucosa, causing varying degrees of diarrhea, vomiting, dehydration, and emaciation in pigs of all ages, with the highest mortality observed in suckling piglets (2, 3). PEDV was first identified in the 1970s in the United Kingdom and Belgium (4). Subsequently, PEDV was first isolated in 1984 in China (5). The impact of PEDV on pigs in China was mitigated due to the use of inactivated vaccines. However, this relative control was disrupted in 2010, a highly virulent PEDV variant with high mortality emerged (6).

At the molecular level, the full-length PEDV genome is approximately 28 kb in size and consists of a 5′ end, a 5′ cap, a 5′ untranslated region (UTR), seven open reading frames (ORFs), a 3′ untranslated region (UTR), and a polyadenylated (polyA) tail (1). These ORFs encode several structural proteins, including the spike protein (S protein), the membrane protein (M protein), the envelope protein (E protein), and the nucleocapsid protein (N protein) (7). Among them, the S protein is encoded by the S gene and plays a crucial role in PEDV invasion of host cells, mediating attachment to and fusion with the host cell membrane (8). Consequently, the S gene is considered a key indicator of PEDV genetic evolution. Based on this gene, PEDV can be classified into genotypes G1a, G1b, G1c (S-INDEL), G2a, G2b, and G2c (9). Structurally, the PEDV S protein is a homotrimeric glycoprotein composed of an N-terminal S1 subunit, which contains the receptor-binding domain, and a C-terminal S2 subunit, which mediates membrane fusion (10). Importantly, certain epitopes of the S protein are involved in inducing the host to produce neutralizing antibodies (11). Among them, the antigenic epitopes such as The Collagenase Equivalence, SS2, SS6, and 2C10 have been identified to elicit the production of neutralizing antibodies in the host (1214). Given their functional significance, monitoring these important antigenic epitopes will aid in the prevention and control of PEDV infection. Nevertheless, studies have shown that these epitopes have undergone varying degrees of mutation (15). Therefore, it remains necessary to identify potential antigenic epitopes on the S protein that can induce the host cells to produce neutralizing antibodies.

In this study, an epidemiological investigation of PEDV was conducted across 20 provinces in China from 2022 to 2025. The S gene of the virus was sequenced, and a phylogenetic tree was constructed. Meanwhile, the changes in antigenic epitopes between the sequenced strains and the vaccine strains were comparatively analyzed. Additionally, ElliPro was used to predict antigenic epitopes located on the surface of the naturally occurring trimeric PEDV S protein. This research will provide valuable guidance for PEDV prevention and control in pig farms, as well as vaccine development.

2 Materials and methods

2.1 Collection of samples

From 2022 to 2025, a total of 2,346 samples from pigs suspected of dying from diarrhea were collected from 20 provinces in China. The samples included pig small intestinal tissues, feces, and anal swabs. The collected samples were placed in sterile 10 mL centrifuge tubes and transported to the laboratory for testing.

2.2 Nucleic acid extraction and detection of PEDV

Viral nucleic acids were extracted using the YALEPIC® Viral RNA Isolation Kit (YALI, Jiangsu) following the manufacturer’s instructions. Subsequently, PEDV detection was performed by RT-qPCR according to the previously established method (16). The RT-qPCR assay was carried out using the HiScript III U + One Step qRT-PCR Probe 5 × Master Mix (Vazyme, Nanjing). Negative and positive plasmid controls were included in each run. Samples with Ct value below 35 were considered positive. Primer and probe sequences are listed in Table 1.

Table 1
www.frontiersin.org

Table 1. The nucleotide sequences used in this study.

2.3 Geographical analysis of PEDV infection

A geographical distribution analysis of PEDV infection was conducted. The visualization of results was performed using GraphPad Prism 10.0, Adobe Illustrator 2022, and the ggplot2 package. Geographic data were obtained from DataV (Accessed on June 11, 2025).1

2.4 Sequencing of the PEDV S gene

Nucleic acid samples with low Ct values from qPCR detection were selected for amplification by RT-PCR. The amplification products were analyzed using 1.5% agarose gel electrophoresis. Bands with high brightness were excised and purified, then sent to Sangon Biotech (Shanghai, China) for sequencing. The sequencing primers are listed in Table 1.

2.5 Sequence similarity analysis and phylogenetic tree construction of PEDV

Using the Distance method in MEGA 11.0, nucleotide and amino acid sequence similarity analyses were performed on 15 sequenced PEDV S genes (Supplementary Table 1). Meanwhile, together with 32 reference sequences selected from the GenBank database (Supplementary Table 3), a phylogenetic tree was constructed using the Neighbor-Joining (NJ) method in MEGA 11.0 (nucleotide sequence similarity >90%). The robustness of the phylogenetic tree was assessed by 1,000 bootstrap replicates.

2.6 Comparison of major antigenic epitopes of PEDV S protein

The amino acid sequences of the major antigenic epitopes COE, SS2, SS6, and 2C10 within the PEDV S protein were aligned using MegAlign software to analyze whether mutations had occurred in these epitopes.

2.7 Prediction of potential N-glycosylation sites specific to the S protein of PEDV virulent strains versus vaccine strains

Use the NetNGlyc-1.0 online server2 to predict potential N-glycosylation sites in the S proteins of the sequenced PEDV field strains and vaccine strains. Only accept results with a score > 0.5 and a Jury agreement of 9/9.

2.8 Prediction of antigenic epitopes in the PEDV S protein

Potential linear and discontinuous antigenic epitopes of the resolved PEDV S protein (PDB: 6VV5) were predicted using ElliPro in IEDB. Peptides with scores greater than 0.7 were selected, and PyMOL was used for the structural visualization of the PEDV S protein.

3 Results

3.1 The nucleotide sequences used in this study

3.2 Amplification status of the PEDV S gene

In this study, 15 virus strains with low Ct values and representative geographic origins were selected for sequencing. The RT-PCR results are shown in Figure 1, where three clear fragments corresponding to S1 (1972 bp), S2 (1,554 bp), and S3 (1,437 bp) were amplified, consistent with the expected sizes.

Figure 1
Gel electrophoresis image showing DNA bands in five lanes. The lanes are labeled DNA Ladder, NTC, S1, S2, S3, and DNA Ladder. The DNA ladder bands are visible on both ends, indicating fragment sizes ranging from 100 to 2000 base pairs. The S1, S2, and S3 lanes display distinct bright bands, while the NTC lane lacks noticeable bands, serving as a control.

Figure 1. The PEDV S gene fragment on a 1.5% agarose gel. The nine lanes represent the corresponding amplified fragments of S1, S2, and S3, and three biological replicates were conducted.

3.3 PEDV infection and geographic distribution

We conducted a survey of PEDV infection in pigs that died of diarrhea across 20 provinces in China. Among the collected samples, the positive rate of PEDV was 43.3% (95% CI: 41.3–45.3%; Figure 2A; Supplementary Table 2). Notably, the positive rates in Liaoning, Anhui, Guangxi, Inner Mongolia, Xinjiang, and Yunnan exceeded 60%. The positive rate in Shaanxi was below 10%, and no PEDV infections were detected in Chongqing. The positive rates were relatively high in Xinjiang, southern regions, and most parts of Northeast China; the coastal areas exhibited moderate positive rates, while the positive rates in inland areas were relatively low (Figure 2B).

Figure 2
Graph A shows a bar chart of positive rates across Chinese regions, ranging from zero to over seventy percent, with upper and lower limits. Graph B displays a map of China colored by positive rate percentages, with a gradient from blue to red. Regions like Xinjiang and Liaoning are highlighted.

Figure 2. PEDV infection and geographic distribution in pigs that died of diarrhea. (A) The bar chart represents the positive rate of PEDV in each province, while the line chart shows the upper and lower limits of the 95% confidence interval. (B) Geographic distribution of PEDV. Colors ranging from green to red indicate increasing positive rates.

3.4 Nucleotide similarity of the PEDV S gene and amino acid similarity of the S protein

To evaluate the protective efficacy of existing vaccines, we performed nucleotide and amino acid similarity analyses between the vaccine strain CV777, AJ1102 and the sequences obtained in this study (Figure 3). The results showed that the nucleotide and amino acid similarities of the S gene from the 15 sequenced PEDV strains ranged from 94.26 to 100% and 93.67 to 100%, respectively. Compared with CV777, the nucleotide and amino acid similarities ranged from 92.71 to 94.83% and 92.89 to 94.99%, respectively. Compared with AJ1102, the nucleotide and amino acid similarities ranged from 95.04 to 99.05% and 95.29 to 98.83%, respectively.

Figure 3
Table comparing two groups labeled A and B, showing percentage similarities between various regions such as AnHui, GanSu, and ShanXi. Colored circles indicate different percentage ranges across both tables.

Figure 3. Nucleotide and amino acid similarity. (A) Nucleotide similarity of the PEDV S gene. (B) Amino acid similarity of the PEDV S protein. Green, yellow, and red circles represent similarity levels from high to low.

3.5 Phylogenetic tree of PEDV

A phylogenetic tree was constructed based on the S gene sequences of 15 sequenced PEDV strains and 32 reference strains retrieved from GenBank (Figure 4). The tree included genotypes G1a, G1b, G1c, G2a, G2b, and G2c, classified according to S gene typing. The results showed that among the 15 sequenced strains, 1 strain (PEDV/HeNan/2) belonged to the G1c (S-INDEL) genotype, 4 strains clustered within G2b, and the remaining 10 strains were classified as G2c, indicating that G2c is currently the predominant subtype in most regions.

Figure 4
Phylogenetic tree illustrating relationships among various PEDV strains. The tree is circular, grouped into five segments labeled G1b, G1a, G1c, G2b, G2a, and G2c. Each segment contains multiple strain labels with G2a and G2c highlighted in shades of blue. Red stars mark specific strains, emphasizing them within the tree.

Figure 4. Phylogenetic tree of the PEDV S gene. Red pentagrams indicate the sequences obtained in this study. Different colors represent different genotypes.

3.6 Comparison of key antigenic epitopes of PEDV S protein

PEDV continues to evolve under the selective pressure of vaccines, with significant variations in its S protein. In this study, 15 newly sequenced strains were compared with the reference strains CV777 and AJ1102. Additionally, the previously identified antigenic epitopes COE, SS2, SS6, and 2C10 were analyzed (Figure 5).

Figure 5
Protein sequence alignment chart displaying amino acid variations across different samples. The chart highlights differences using colored columns for specific positions. Dots represent conserved amino acids, and letters indicate substitutions or mutations at those positions. Sequence positions and sample identifiers are provided along the top and left sides.

Figure 5. Comparison of key antigenic epitopes of the PEDV S protein. The green, orange, blue, and yellow regions correspond to the antigenic epitopes COE, SS2, SS6, and 2C10, respectively.

The results showed that, relative to CV777, the newly sequenced strains exhibited 3 to 5 amino acid mutations in the COE region. In the SS2 epitope, a mutation (Asn755 to Asp755) was identified exclusively in the JiangXi/2 strain. Within the SS6 epitope, all newly sequenced strains shared a common amino acid substitution (Ser771 to Tyr771); furthermore, the GanSu and JiangXi/2 strains each harbored an additional mutation (Val774 to Ala774 and Gln773 to Leu773, respectively). In the 2C10 epitope, the GanSu and JiangXi/2 strains carried unique amino acid changes (Pro1375 to Ser1375 and Ile1379 to Cys1379, respectively).

When compared to AJ1102, the new strains showed 0 to 4 amino acid mutations in the COE region. The mutations observed in the SS2 epitope were consistent with those found relative to CV777. For the SS6 epitope, only the GanSu and JiangXi/2 strains exhibited amino acid substitutions (Val774 to Ala774 and Gln773 to Leu773, respectively). In the 2C10 epitope, the mutation was the same as that observed in comparison to CV777.

3.7 Prediction results of potential specific N-glycosylation sites in the S protein of 15 PEDV strains compared with the vaccine strain

Most variants at position 62 (NSTW) and position 118 (NATA) are similar to AJ1102 (Figure 6; Supplementary Table 5). However, CV777 has NSSW (Thr to Ser) at position 62 and a deletion at position 118, which may result in partial loss of protective efficacy. Notably, the PEDV/GuangXi strain, classified as G2c subtype, shows deletions of N-glycosylation sites at positions 353 (NSSD), 1,263 (NRTG), and 1275 (NATY) compared to AJ1102, which might lead to immune evasion. Compared to CV777, the glycosylation patterns in regions 62 and 118 differ from those of other strains, which is associated with the ability to evade traditional vaccines. The PEDV/XinJiang/2 strain exhibits a novel potential N-glycosylation site at position 302 (NKTI), which could alter the antigenic activity and structure of the S protein.

Figure 6
Grid showing amino acid changes across different strains. Strains are listed vertically, amino acid positions horizontally. Red cells indicate changes such as NSSW, NKTL, and NSSN, surrounded by teal indicating no change. Positions include 62, 118, 131, and others up to 1275.

Figure 6. Prediction results of potential specific N-glycosylation sites in the S protein of 15 PEDV strains compared with the vaccine strain. The vertical axis represents the strain names, and the horizontal axis represents the amino acid sites. The blue areas indicate the same N-glycosylation at the same amino acid sites, while the red areas indicate differences.

3.8 Prediction of antigenic epitopes of the PEDV S protein

In this study, the antigenic epitopes of the PEDV S protein were predicted using ElliPro. Seven peptide segments with antigenicity scores greater than 0.7 were selected (Supplementary Table 4) and visualized (Supplementary Figures 1–7). Under natural conditions, the PEDV S protein exists in a trimeric state, and the structure revealed potential antigenic epitopes on the S protein in its trimeric conformation. Conserved regions among the sequenced strains and the reference strains CV777 and AJ1102 were visualized (Figures 7, 8). The potential conserved antigenic epitopes were identified in the following regions: amino acids 31–54, 90–103, 1,170–1,177, 1,179–1,197, 1,201–1,211, 1,221–1,236, and 1,243–1,254. Additionally, discontinuous (conformational) epitope prediction based on the PEDV S protein structure was performed. Two groups of discontinuous epitopes with antigenicity scores above 0.7 were selected and visualized (Supplementary Table 4; Supplementary Figures 8, 9).

Figure 7
Diagram showing A) amino acid sequence alignment of different strains, highlighting variations; B-E ) Visualization of conserved antigenic epitopes. The blue, green, and white regions represent the A, B, and C chains of the S protein trimer, respectively. Yellow spheres indicate potential antigenic epitope regions with scores greater than 0.7. (F) An enlarged view of the key amino acid sites of the S protein monomer, where the red areas indicate the key amino acid sites of the S protein monomer.

Figure 7. Visualization of conserved regions within continuous epitopes of the S protein in its native trimeric state. (A) Amino acid alignment of the sequenced strains and vaccine strains in the predicted epitope regions. (B–E) Visualization of conserved antigenic epitopes. The blue, green, and white regions represent the A, B, and C chains of the S protein trimer, respectively. Yellow spheres indicate potential antigenic epitope regions with scores greater than 0.7. (F) An enlarged view of the key amino acid sites of the S protein monomer, where the red areas indicate the key amino acid sites of the S protein monomer.

Figure 8
A panel of images labeled A to F. A: An alignment chart of amino acid sequences for different strains, displaying sequence conservation with some variations marked. B-E:Visualization of conserved antigenic epitopes. The blue, green, and white regions represent the A, B, and C chains of the S protein trimer, respectively. Yellow spheres indicate potential antigenic epitope regions with scores greater than 0.7. (F) An enlarged view of the key amino acid sites of the S protein monomer, where the red areas indicate the key amino acid sites of the S protein monomer.

Figure 8. Visualization of conserved regions within continuous epitopes of the S protein in its native trimeric state. (A) Amino acid alignment of the sequenced strains and vaccine strains in the predicted epitope regions. (B–E) Visualization of conserved antigenic epitopes. The blue, green, and white regions represent the A, B, and C chains of the S protein trimer, respectively. Yellow spheres indicate potential antigenic epitope regions with scores greater than 0.7. (F) An enlarged view of the key amino acid sites of the S protein monomer, where the red areas indicate the key amino acid sites of the S protein monomer.

4 Discussion

PEDV can cause diarrhea, vomiting, and dehydration in piglets, and severe cases may result in death (17). It is necessary to continuously update disease prevention and control strategies as well as vaccine development plans to effectively control and prevent PEDV infection. The primary transmission route of PEDV is fecal-oral and nasal transmission (18). From 2011 to 2014, the PEDV positive rate in 29 provinces of China ranged from approximately 61.1 to 78.5% (19). Epidemiological investigations from 2015 to 2018 showed that the positive rate of PEDV in China was approximately 66.9% (363/543) (20). From 2018 to 2021, the positive rate of PEDV in the southern region of China was 55.8% (3,986/7107) (21). Epidemiological data from seven regions in China between 2020 and 2022 showed that the positive rate of PEDV was approximately 57.4% (199/347) (22). This study investigated the PEDV positive rate in 20 provinces of China from 2022 to 2025, which was 43.3% (1,015/2346). Regions with high positive rates were mainly concentrated in southern and northeastern China, while regions with low positive rates were primarily located in the western and northern parts of the country. These areas may have maintained lower positive rates due to factors such as lower pig farm density, geographical isolation, limited mobility, or effective prevention and control measures.

The similarity results of the PEDV S gene and amino acids compared with CV777 and AJ1102 indicate a significant difference from CV777, suggesting that the vaccine’s protective efficacy may be weakened. In contrast, there is a relatively high similarity with AJ1102, but a certain degree of mutation has also occurred, which may affect the efficacy of the vaccine in future evolution.

The S protein of PEDV contains numerous antigenic epitopes that can induce the host to produce neutralizing antibodies, and most of these epitopes are located in the S1 region (23, 24). Under the selective pressure exerted by vaccines, mutations occur in the PEDV S protein (25). When mutations occur in the antigenic epitopes on the S protein, PEDV’s ability to bind antibodies is significantly reduced (26). Therefore, monitoring whether antigenic epitopes on the PEDV S protein undergo changes is of great significance for the prevention and control of PEDV and the development of vaccines. Previous studies have identified several important antigenic epitopes. The Core Neutralizing Epitope (COE), combined with an adjuvant and administered into mice, can induce the production of a large amount of PEDV neutralizing antibodies (27, 28). A PEDV vaccine based on the COE antigenic epitope has already been developed (27). In addition, the SS2, SS6, and 2C10 epitopes have also been identified (29, 30).

This study analyzed the previously identified COE, SS2, SS6, and 2C10 epitopes, finding that the newly sequenced strains exhibit varying degrees of mutations in these epitope regions. These mutations may potentially affect the protective efficacy of current mainstream vaccines, indicating the need for further monitoring.

The extensive N-glycosylation of the S protein is a prerequisite for the virus to enter host cells, complete membrane fusion, and carry out subsequent replication. An intact N-glycosylation pathway is essential for PEDV proliferation (31). The addition, loss, or positional changes of N-glycosylation sites can create a glycan shield on the surface of the S protein, masking neutralizing epitopes and thereby helping the virus evade host antibody recognition. The GII subtype (such as GIIc) has acquired new glycosylation sites (N62, N118) in the Domain 0 region of the S protein, which is believed to be associated with its high pathogenicity and ability to escape traditional vaccines (32). Therefore, monitoring the N-glycosylation sites of the PEDV S protein is also of certain significance for future vaccine development.

Linear and discontinuous epitopes hold significant importance in immunology. Linear antigen epitopes are composed of consecutively arranged amino acids and generally do not depend on the three-dimensional structure of the antigen. B cells can recognize linear epitopes exposed on the surface of antigens and, upon activation, produce specific antibodies targeting these epitopes. Discontinuous epitopes, which are also one of the main types of epitopes recognized by B cells, consist of residues that are brought into close proximity through the spatial folding of the protein, forming a specific three-dimensional structure that is recognized by antibodies or B cell receptors (33).

ElliPro is an online tool that can be used to identify antigenic epitopes in structural proteins. ElliPro is considered a potentially useful research tool for identifying antigens epitopes (34). Previously, researchers used ElliPro to predict linear B-cell epitopes on the envelope protein of the Zika virus (35). In 2023, 24 linear epitopes and 7 discontinuous epitopes of avian-origin Escherichia coli were predicted using ElliPro, and two of these peptide segments were selected as candidate peptides (36). In addition, other researchers used the IEDB ElliPro server to predict discontinuous and linear epitope-specific B-cell epitopes of the SARS-CoV-2 spike protein (37).

This study used ElliPro to predict 7 linear epitopes and 2 discontinuous epitopes of the PEDV S protein in its trimeric state, which has already been resolved. Conserved epitopes were identified among those with high scores. However, these antigenic epitopes based on bioinformatics analysis still require further validation through in vitro and in vivo experiments.

5 Conclusion

In summary, this study revealed that among 2,346 pig samples from pigs that died of diarrhea collected across 20 provinces in China during 2022–2025, the PEDV positive rate was 43.3% (1,015/2,346). Three subgroups were identified, among which 1 strains belonged to G1c, 4 strains to G2b, and the remaining 10 strains to G2c, indicating that G2c is currently the dominant subtype. The nucleotide and amino acid similarities between the sequenced strains and the CV777 strain ranged from 92.71 to 94.83% and 92.89 to 94.99%, respectively. Mutations of varying degrees were observed in the neutralizing epitopes COE, SS2, SS6, and 2C10. Prediction of antigenic epitopes identified aa31–54, 90–103, 1,170–1,177, 1,179–1,197, 1,201–1,211, 1,221–1,236, and 1,243–1,254 as conserved potential antigenic regions. These findings provide important references and insights for developing effec-tive PED control strategies, as well as for the research and development of vaccines and diagnostic methods.

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 author/s.

Ethics statement

Ethical approval was not required for the studies involving animals in accordance with the local legislation and institutional requirements because all samples in this study (feces from diseased pigs and intestinal tissues from deceased pigs) were proactively submitted to our laboratory by veterinarians from pig farms experiencing outbreaks of piglet diarrhea for the purpose of identifying the causative agents. The samples were collected from naturally diseased or deceased individuals at clinical outbreak farms. Fecal collection was a non-invasive procedure (using sterile swabs or fecal scoops to collect only naturally excreted feces) and did not cause additional stress or harm to the animals. Intestinal tissues from deceased pigs were collected after natural death. This study is classified as non-interventional and therefore does not require additional committee approval or exemption were required. Written informed consent was not obtained from the owners for the participation of their animals in this study because since the pig farms voluntarily submitted diarrheic deceased pigs to the laboratory for clinical diagnosis, no additional certification is required.

Author contributions

ZW: Investigation, Resources, Conceptualization, Validation, Project administration, Methodology, Writing – original draft, Formal analysis, Data curation. KY: Writing – review & editing, Formal analysis, Writing – original draft, Methodology, Project administration, Visualization, Investigation, Software, Conceptualization, Validation. MB: Investigation, Software, Writing – original draft, Project administration, Methodology. KL: Methodology, Investigation, Writing – original draft, Software, Visualization. WW: Writing – original draft, Investigation, Resources, Methodology. YS: Methodology, Investigation, Resources, Writing – original draft. XP: Writing – original draft, Resources, Software, Investigation, Methodology. TL: Supervision, Data curation, Conceptualization, Software, Investigation, Writing – review & editing, Writing – original draft, Methodology, Resources, Funding acquisition, Project administration, Visualization, Formal analysis, Validation. XM: Methodology, Conceptualization, Investigation, Visualization, Validation, Writing – review & editing, Supervision, Formal analysis, Writing – original draft, Resources, Software, Data curation, Project administration, Funding acquisition.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

ZW, WW, YS, XP, and TL are employed by the company Tecon Bio-Pharmaceuticals Co., Ltd.

The remaining 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.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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

Footnotes

References

1. Kocherhans, R, Bridgen, A, Ackermann, M, and Tobler, K. Completion of the porcine epidemic diarrhoea coronavirus (Pedv) genome sequence. Virus Genes. (2001) 23:137–44. doi: 10.1023/A:1011831902219

PubMed Abstract | Crossref Full Text | Google Scholar

2. Lee, C. Porcine epidemic diarrhea virus: an emerging and re-emerging epizootic swine virus. Virol J. (2015) 12:193. doi: 10.1186/s12985-015-0421-2

PubMed Abstract | Crossref Full Text | Google Scholar

3. Song, D, Moon, H, and Kang, B. Porcine epidemic diarrhea: a review of current epidemiology and available vaccines. Clin Exp Vaccine Res. (2015) 4:166–76. doi: 10.7774/cevr.2015.4.2.166

PubMed Abstract | Crossref Full Text | Google Scholar

4. Antas, M, and Wozniakowski, G. Current status of porcine epidemic diarrhoea (ped) in European pigs. J Vet Res. (2019) 63:465–70. doi: 10.2478/jvetres-2019-0064

PubMed Abstract | Crossref Full Text | Google Scholar

5. Zhang, H, Zou, C, Peng, O, Ashraf, U, Xu, Q, Gong, L, et al. Global dynamics of porcine enteric coronavirus Pedv epidemiology, evolution, and transmission. Mol Biol Evol. (2023) 40:msad052. doi: 10.1093/molbev/msad052

PubMed Abstract | Crossref Full Text | Google Scholar

6. Chen, J, Wang, C, Shi, H, Qiu, HJ, Liu, S, Shi, D, et al. Complete genome sequence of a Chinese virulent porcine epidemic diarrhea virus strain. J Virol. (2011) 85:11538–9. doi: 10.1128/JVI.06024-11

PubMed Abstract | Crossref Full Text | Google Scholar

7. Duarte, M, Tobler, K, Bridgen, A, Rasschaert, D, Ackermann, M, and Laude, H. Sequence analysis of the porcine epidemic diarrhea virus genome between the nucleocapsid and spike protein genes reveals a polymorphic Orf. Virology. (1994) 198:466–76. doi: 10.1006/viro.1994.1058

PubMed Abstract | Crossref Full Text | Google Scholar

8. Chambers, P, Pringle, CR, and Easton, AJ. Heptad repeat sequences are located adjacent to hydrophobic regions in several types of virus fusion glycoproteins. J Gen Virol. (1990) 71:3075–80.

Google Scholar

9. Han, X, Liu, Y, Wang, Y, Wang, T, Li, N, Hao, F, et al. Isolation and characterization of porcine epidemic diarrhea virus with a novel continuous mutation in the S1(0) domain. Front Microbiol. (2023) 14:1203893. doi: 10.3389/fmicb.2023.1203893

PubMed Abstract | Crossref Full Text | Google Scholar

10. Kirchdoerfer, RN, Bhandari, M, Martini, O, Sewall, LM, Bangaru, S, Yoon, KJ, et al. Structure and immune recognition of the porcine epidemic diarrhea virus spike protein. Structure. (2021) 29:e5:385–92. doi: 10.1016/j.str.2020.12.003

Crossref Full Text | Google Scholar

11. Chang, SH, Bae, JL, Kang, TJ, Kim, J, Chung, GH, Lim, CW, et al. Identification of the epitope region capable of inducing neutralizing antibodies against the porcine epidemic diarrhea virus. Mol Cell. (2002) 14:295–9. doi: 10.1016/S1016-8478(23)15106-5

PubMed Abstract | Crossref Full Text | Google Scholar

12. Godet, M, Grosclaude, J, Delmas, B, and Laude, H. Major receptor-binding and neutralization determinants are located within the same domain of the transmissible gastroenteritis virus (coronavirus) spike protein. J Virol. (1994) 68:8008–16.

Google Scholar

13. Sun, D, Feng, L, Shi, H, Chen, J, Cui, X, Chen, H, et al. Identification of two novel B cell epitopes on porcine epidemic diarrhea virus spike protein. Vet Microbiol. (2008) 131:73–81. doi: 10.1016/j.vetmic.2008.02.022

PubMed Abstract | Crossref Full Text | Google Scholar

14. Tian, Y, Yang, X, Li, H, Ma, B, Guan, R, Yang, J, et al. Molecular characterization of porcine epidemic diarrhea virus associated with outbreaks in Southwest China during 2014-2018. Transbound Emerg Dis. (2021) 68:3482–97. doi: 10.1111/tbed.13953

PubMed Abstract | Crossref Full Text | Google Scholar

15. Zhao, X, Li, Z, Zeng, X, Zhang, G, Niu, J, Sun, B, et al. Sequence analysis of the spike gene of porcine epidemic diarrhea virus isolated from South China during 2011-2015. J Vet Sci. (2017) 18:237–43. doi: 10.4142/jvs.2017.18.2.237

PubMed Abstract | Crossref Full Text | Google Scholar

16. Huang, X, Chen, J, Yao, G, Guo, Q, Wang, J, and Liu, G. A TaqMan-probe-based multiplex real-time rt-qpcr for simultaneous detection of porcine enteric coronaviruses. Appl Microbiol Biotechnol. (2019) 103:4943–52. doi: 10.1007/s00253-019-09835-7

PubMed Abstract | Crossref Full Text | Google Scholar

17. Le Dafniet, M, Blumberg-Tick, J, Gozlan, H, Barret, A, Joubert Bression, D, and Peillon, F. Altered balance between thyrotropin-releasing hormone and dopamine in prolactinomas and other pituitary tumors compared to normal pituitaries. J Clin Endocrinol Metab. (1989) 69:267–71.

Google Scholar

18. Jung, K, Saif, LJ, and Wang, Q. Porcine epidemic diarrhea virus (Pedv): an update on etiology, transmission, pathogenesis, and prevention and control. Virus Res. (2020) 286:198045. doi: 10.1016/j.virusres.2020.198045

PubMed Abstract | Crossref Full Text | Google Scholar

19. Wang, D, Fang, L, and Xiao, S. Porcine epidemic diarrhea in China. Virus Res. (2016) 226:7–13. doi: 10.1016/j.virusres.2016.05.026

PubMed Abstract | Crossref Full Text | Google Scholar

20. Su, M, Li, C, Qi, S, Yang, D, Jiang, N, Yin, B, et al. A molecular epidemiological investigation of Pedv in China: characterization of co-infection and genetic diversity of S1-based genes. Transbound Emerg Dis. (2020) 67:1129–40. doi: 10.1111/tbed.13439

PubMed Abstract | Crossref Full Text | Google Scholar

21. Li, C, Lu, H, Geng, C, Yang, K, Liu, W, Liu, Z, et al. Epidemic and evolutionary characteristics of swine enteric viruses in south-Central China from 2018 to 2021. Viruses. (2022) 14:1420. doi: 10.3390/v14071420

PubMed Abstract | Crossref Full Text | Google Scholar

22. Peng, Q, Fu, P, Zhou, Y, Lang, Y, Zhao, S, Wen, Y, et al. Phylogenetic analysis of porcine epidemic Diarrhea virus (Pedv) during 2020-2022 and isolation of a variant recombinant Pedv strain. Int J Mol Sci. (2024) 25:10878. doi: 10.3390/ijms252010878

PubMed Abstract | Crossref Full Text | Google Scholar

23. Fernandez, A. Sars-CoV-2 glycosylation suggests that vaccines should have adopted the S1 subunit as antigen. Acs Pharmacol Transl Sci. (2021) 4:1016–7. doi: 10.1021/acsptsci.1c00036

PubMed Abstract | Crossref Full Text | Google Scholar

24. Follis, KE, York, J, and Nunberg, JH. Furin cleavage of the Sars coronavirus spike glycoprotein enhances cell-cell fusion but does not affect virion entry. Virology. (2006) 350:358–69. doi: 10.1016/j.virol.2006.02.003

PubMed Abstract | Crossref Full Text | Google Scholar

25. Sun, J, Li, Q, Shao, C, Ma, Y, He, H, Jiang, S, et al. Isolation and characterization of Chinese porcine epidemic diarrhea virus with novel mutations and deletions in the S gene. Vet Microbiol. (2018) 221:81–9. doi: 10.1016/j.vetmic.2018.05.021

PubMed Abstract | Crossref Full Text | Google Scholar

26. Li, C, Li, W, Lucio De Esesarte, E, Guo, H, Van Den Elzen, P, Aarts, E, et al. Cell attachment domains of the porcine epidemic Diarrhea virus spike protein are key targets of neutralizing antibodies. J Virol. (2017) 91:10–1128. doi: 10.1128/JVI.00273-17

PubMed Abstract | Crossref Full Text | Google Scholar

27. Wang, X, Wang, L, Huang, X, Ma, S, Yu, M, Shi, W, et al. Oral delivery of probiotics expressing dendritic cell-targeting peptide fused with porcine epidemic Diarrhea virus Coe antigen: a promising vaccine strategy against Pedv. Viruses. (2017) 9:312. doi: 10.3390/v9110312

PubMed Abstract | Crossref Full Text | Google Scholar

28. Yao, X, Zhu, Y, Qiao, WT, Lu, WH, Zhang, YQ, and Li, JL. Based on the results of PEDV phylogenetic analysis of the most recent isolates in China, the occurrence of further mutations in the antigenic site S1 degrees and coe of the S protein which is the target protein of the vaccine. Transbound Emerg Dis. (2023) 2023:1227110. doi: 10.1155/2023/1227110

Crossref Full Text | Google Scholar

29. Li, Z, Ma, Z, Li, Y, Gao, S, and Xiao, S. Porcine epidemic diarrhea virus: molecular mechanisms of attenuation and vaccines. Microb Pathog. (2020) 149:104553. doi: 10.1016/j.micpath.2020.104553

PubMed Abstract | Crossref Full Text | Google Scholar

30. Okda, FA, Lawson, S, Singrey, A, Nelson, J, Hain, KS, Joshi, LR, et al. The S2 glycoprotein subunit of porcine epidemic diarrhea virus contains immunodominant neutralizing epitopes. Virology. (2017) 509:185–94. doi: 10.1016/j.virol.2017.06.013

PubMed Abstract | Crossref Full Text | Google Scholar

31. Zhu, H, Lou, J, Yang, Z, Bai, J, Jiang, P, Wang, X, et al. Stt3B promotes porcine epidemic diarrhea virus replication by regulating N-glycosylation of Pedv S protein. J Virol. (2025) 99:e0001825. doi: 10.1128/jvi.00018-25

PubMed Abstract | Crossref Full Text | Google Scholar

32. Fu, Y, Wang, Y, Dai, L, Cheng, B, Xiao, S, and Yin, Y. Evolutionary dynamics and antigenic diversity of porcine epidemic diarrhea virus (Pedv) in China: phylogenetic and recombination analyses based on large-scale S gene sequences. BMC Vet Res. (2025) 21:426. doi: 10.1186/s12917-025-04873-y

PubMed Abstract | Crossref Full Text | Google Scholar

33. Kumar, N, Bajiya, N, Patiyal, S, and Raghava, GPS. Multi-perspectives and challenges in identifying B-cell epitopes. Protein Sci. (2023) 32:e4785. doi: 10.1002/pro.4785

PubMed Abstract | Crossref Full Text | Google Scholar

34. Ponomarenko, J, Bui, HH, Li, W, Fusseder, N, Bourne, PE, Sette, A, et al. ElliPro: a new structure-based tool for the prediction of antibody epitopes. Bmc Bioinformatics. (2008) 9:514. doi: 10.1186/1471-2105-9-514

PubMed Abstract | Crossref Full Text | Google Scholar

35. Srivastava, K, and Srivastava, V. Prediction of conformational and linear B-cell epitopes on envelop protein of zika virus using Immunoinformatics approach. Int J Pept Res Ther. (2023) 29:17. doi: 10.1007/s10989-022-10486-y

PubMed Abstract | Crossref Full Text | Google Scholar

36. Grakh, K, Mittal, D, Prakash, A, Kumar, R, and Jindal, N. uspA gene-based phylogenetic analysis and antigenic epitope prediction for Escherichia coli strains of avian origin. Front Vet Sci. (2023) 10:1183048. doi: 10.3389/fvets.2023.1183048

PubMed Abstract | Crossref Full Text | Google Scholar

37. Moustafa, RI, Faraag, AHI, El-Shenawy, R, Agwa, MM, and Elsayed, H. Harnessing immunoinformatics for developing a multiple-epitope peptide-based vaccination approach against Sars-CoV-2 spike protein. Saudi J Biol Sci. (2023) 30:103661. doi: 10.1016/j.sjbs.2023.103661

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: PEDV, S protein, molecular characteristics, mutation analysis, N-glycosylation, antigen prediction

Citation: Wang Z, Yang K, Bi M, Li K, Wang W, Song Y, Pan X, Li T and Mo X (2025) Molecular characteristics and potential antigenic epitope analysis of porcine epidemic diarrhea virus in China from 2022 to 2025. Front. Vet. Sci. 12:1667063. doi: 10.3389/fvets.2025.1667063

Received: 16 July 2025; Accepted: 22 September 2025;
Published: 08 October 2025.

Edited by:

Kimberly VanderWaal, University of Minnesota Twin Cities, United States

Reviewed by:

Xiaoliang Hu, Yibin University, China
Xin Yao, Northeast Agricultural University, China

Copyright © 2025 Wang, Yang, Bi, Li, Wang, Song, Pan, Li and Mo. 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: Tianzeng Li, Mjg5NTkzMjM2QHFxLmNvbQ==; Xiaobing Mo, bW94QGpsdS5lZHUuY24=

These authors have contributed equally to this work

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