Molecular Characteristics and Pathogenicity of Porcine Epidemic Diarrhea Virus Isolated in Some Areas of China in 2015–2018

Abstract

Since 2010, Porcine epidemic diarrhea virus (PEDV) has caused severe diarrhea disease in piglets in China, resulting in large economic losses. To understand the genetic characteristics of the PEDV strains that circulated in some provinces of China between 2015 and 2018, 375 samples of feces and small intestine were collected from pigs and tested. One hundred seventy-seven samples tested positive and the PEDV-positive rate was 47.20%. A phylogenetic tree analysis based on the entire S gene showed that these strains clustered into four subgroups, GI-a, GI-b, GII-a, and GII-b, and that the GII-b strains have become dominant in recent years. Compared with previous strains, these strains have multiple variations in the SP and S1-NTD domains and in the neutralizing epitopes of the S protein. We also successfully isolated and identified a new virulent GII-b strain, GDgh16, which is well-adapted to Vero cells and caused a high mortality rate in piglets in challenge experiments. Our study clarifies the genetic characteristics of the prevalent PEDV strains in parts of China, and suggests that the development of effective novel vaccines is both necessary and urgent.

Introduction

Porcine epidemic diarrhea virus (PEDV) is the etiological agent of porcine epidemic diarrhea (PED), a severe diarrhea disease in piglets that is characterized by severe watery diarrhea, vomiting, dehydration, weight loss, and nearly 100% mortality (1). PED occurred sporadically around the world in 1990–2009, but in 2010, an acute and severe outbreak of PED in piglets occurred in China and spread to other Asian countries, causing large economic losses (2–8). In April 2013, PED suddenly erupted in the United States, causing many piglets to die, and the mortality rate in suckling piglets reached 100% (9, 10). The disease was shown to be caused by a highly pathogenic PEDV variant. The S genes of the classical CV777 strain and the new strain OH851 have the same insertions and deletions (S-INDEL strains), unlike those of the variant PEDV strains (11, 12).

The genome of PEDV is ~28 kb in length and contains seven open reading frames (ORFs), which encode four structural proteins and three non-structural proteins (13). S is the largest structural protein, and contains neutralizing antibody epitopes and a specific receptor-binding site for viral entry (14). At present, four antigenic epitopes have been characterized in the S protein, including the CO equivalent (COE) domain (amino acids 499–638), the epitopes SS2 (amino acids 748–755) and SS6 (amino acids 764–771), and epitope 2C10 (1368-GPRLQPY-1374) (15, 16). Because the S protein plays a vital role and the S gene is extensively mutated, it is often used as the target gene in the analysis of viral genetic variation. Based on whether the S gene contains the INDEL sequence or not, PEDV strains can be classified into genogroup II (GII) or genogroup I (GI), respectively. GI is further divided into two subgroups (GI-a, GI-b) according to INDEL sequence differences. At present, most isolates recovered in China belong to GII (17). A new mutation in the S gene of PEDV has recently been reported (18). Different GI-b strains have also been reported in different areas of China (19, 20). Studies have shown that PEDV strains of different genotypes can coexist, in one province in particular. These findings indicate that PEDV has continued to spread widely to most areas of China and has caused serious economic losses in the pig industry, reflecting the complex evolution of the virus. Therefore, extensive research into the evolutionary pathogenic mechanism of these strains in China is essential.

To control the spread of PEDV, a classical-CV777-derived vaccine has been widely used in many areas of China. However, it does not provide adequate protection against PEDV invasion (6, 21). In contrast, the wide-scale use of vaccines has increased the environmental stress upon the virus, causing PEDV to mutate to escape its host's immune defenses. To further and fully understand the prevalence and evolution of PEDV in southern China, diarrhea samples were collected from piglets in this study, and the variation of the S genes of the PEDV-positive samples were analyzed with sequence alignment and a phylogenetic tree.

Materials and Methods

Sample Collection

A total of 375 diarrheic samples from the small intestine tissues or feces were collected from suckling piglets on pig farms in eight provinces of China (Fujian, Guangdong, Guangxi, Guizhou, Jiangxi, Shandong, Hubei, Hu'nan, and Hainan) between June 2015 and October 2018. The piglets suffered severe watery diarrhea and dehydration. The diarrheic feces were resuspended in 1 mL of phosphate-buffered saline (PBS) in 1.5 mL Eppendorf tubes. After centrifugation at 10,000 × g for 5 min, 200 μL of each supernatant was transferred to a new tube for RNA extraction and virus isolation.

RNA Extraction and Sequencing

The total RNA from the collected supernatants was extracted with TRIzol Reagent (TaKaRa), according to the manufacturer's instructions. The extracted RNA was subjected to reverse transcription (RT–PCR) with three pairs of newly designed primers to amplify and detect the PEDV S gene (Table 1). The three overlapping PCR products were identified with 1.5% agarose gel electrophoresis. The positive PCR products were sequenced by Sangon Biological Engineering Co. Ltd, and the entire sequence of the S gene was determined with the DNAStar software. The complete S gene sequences were submitted to GenBank, under the accession numbers shown in Table 2.

S Gene Sequence Analysis

The complete genome sequences of reference strains available in GenBank were downloaded and used in a phylogenetic analysis (Table 3). A phylogenetic tree was constructed from all the S genes of the representative strains and isolates, using the neighbor joining method with 1,000 bootstrap replicates, with the Molecular Evolutionary Genetics Analysis (MEGA, version 6.0) software (22).

Virus Isolation

Vero cells grown in a 24-well cell culture plate were infected with the previously prepared supernatants and maintained in Dulbecco's modified Eagle's medium (Thermo Scientific) containing 7 μg/mL trypsin without EDTA (Thermo Scientific). The cells were monitored daily for a cytopathic effect (CPE). When the CPE appeared in 70% of the cells, the cells were fixed with anhydrous ethanol. An immunofluorescence assay (IFA) was then performed with an anti-N protein monoclonal antibody (mAb; cat. # PEDV12-F, Alpha Diagnostic International Inc., USA) diluted 1:1,000 and an Alexa-Fluor®-488-conjugated Affinipure goat anti-mouse IgG(H+L) secondary antibody (SA00013-1; Proteintech, USA) diluted 1:400.

Titer Determination for the Viral Proliferation Curve

Vero cells cultured in a 24-well cell culture plate were infected with PEDV at a multiplicity of infection (MOI) of 0.01. The cells and supernatants were collected at 12, 24, 36, 48, 60, 72, and 96 h post-infection (hpi). The cells were then frozen and thawed three times. After centrifugation at 10,000 × g for 5 min at 4°C, the supernatants were collected and the median tissue culture infective dose (TCID₅₀) was determined with a microtitration infection assay.

Piglet Challenge Experiment

To determine the virulence of the third-generation isolated strain GDgh16, six healthy 4-day-old colostrum-deprived suckling piglets were artificially fed bovine milk from birth. The colostrum-deprived piglets were randomly divided into two groups, with three piglets in each group. One group was challenged orally with 0.5 mL of PEDV at 10^5.0 TCID₅₀/mL. The other group received cell-culture medium. Duplicate samples of small intestine were collected in from all piglets, which had been euthanized at 48 h postchallenge. One of the duplicate samples was crushed in a grinder with 2 mL of PBS. The crushed intestine was then centrifuged at 10,000 × g for 10 min at 4°C. The supernatant was collected and their RNA extracted. The PEDV N gene copies in the small intestine were detected with real-time quantitative PCR (qPCR). RT–qPCR was performed with the PowerUp™ SYBR® Green Master Mix (A25742; Thermo Fisher) in a 20 μL reaction containing 1 ng of cDNA as the template, in a CFX96 thermal cycler, under the following cycling conditions: 50°C for 2 min; 95°C for 2 min; and 40 cycles of 95°C for 15 s, 60°C for 15 s, and 72°C for 60 s. The other sample was stained with anti-N protein mAb (diluted 1:1,000) for immunohistochemical (IHC) examination.

Statistical Analysis

The numerical data are expressed as means ± standard deviations (SD), and all data were analyzed with the GraphPad Prism software (version 5.02 for Windows; GraphPad Software Inc.).

Results

PEDV Detection and Phylogenetic Analysis Based on the S Gene

As shown in Table 4, of the 375 feces and small -intestine samples tested between 2015 and 2018, 177 were positive for PEDV (47.20%). The positivity rates in 2015, 2016, 2017, and 2018 were 48.57% (34 positive samples and 70 test samples), 67.14%% (47 positive samples and 70 test samples), 53.33% (32 positive samples and 60 test samples), and 36.57% (64 positive samples and 175 test samples), respectively. The positive rate was highest in 2016 and lowest in 2018. A sequence alignment showed that these strains shared 92.9–100% nucleotide homology and 91–100% amino acid identity. They also shared 93.1–96.8% nucleotide homology and 91.5–96.8% amino acid identity with reference strain CV777, and 93.8–99% nucleotide homology and 92.5–98.9% amino acid identity with the reference strains isolated from China.

Sixty-two S genes from the test strains and representative strains downloaded from GenBank were analyzed with a phylogenetic tree. As shown in Figure 1, the phylogenetic analysis divided these strains into two groups, GI and GII, based on whether the S gene contained the S-INDEL (23). GI included the classical strains (CV777 and SM98) and some isolates from China, the USA, and Japan collected after 2010. Therefore, GI was further divided into three subgroups: GI-a, GI-b. GI-a contained classical S-INDEL strains. GI-b contained a new S-INDEL strain. GII contained non-S-INDEL strains and was also divided into two subgroups, GII-a and GII-b, which consisted of a number of extremely virulent strains from all over the world, isolated since 2010. The strains isolated in the present study belonged to GI-a, GI-b, GII-a, and GII-b. GDjm18-2, was categorized as subtype GI-a, which also included the classical vaccine strains CV777-attenuated and JS2008. GDjm17-1 was categorized in GI-b cluster. The other strains identified in the present study formed eight clusters. Of these strains, 25 isolates from Guangdong, three isolates from Fujian, and one isolate from Jiangxi formed three clusters and belonged to GII-b, with strong similarity to GD-A and CH-GXNN-2012. The other 34 isolates formed five clusters and belonged to GII-a. Among these 34 strains, JXyc15 was closely related to the C4 cluster (North American strains), whereas the other strains showed closer identity to CH-ZMDZ-11, CH-HNAY-2015, and CH-HNCDE-2016L. As shown in Table 5, all the strains isolated in 2015 belonged to GII-a (100%). In 2016 and 2017, 46.15% and 43.75% of the isolated strains belonged to GII-a, respectively. Compared with GII-a, the rate slightly increased, and in 2016 and 2017, 53.84 and 50% of the isolated strains belonged to GII-b, respectively. However, in 2018, 72.22% of the isolated strains belonged to GII-b, which was much higher than the proportion that belonged to GII-a in 2017 (22.22%). The comparison result show: Variation of PEDV S gene is continuously occurring and GII-b strains may be the dominant strains in China in the future.

Figure 1

Amino Acid Sequence Analysis of Neutralizing Epitopes in the S Protein

Neutralizing antibodies play an important role in the prevention and control of viral infections. Therefore, it is important to identify and analyze the amino acid sequences of the neutralizing epitopes in viral proteins. To analyze the genetic characteristics of the South China PEDV strains, the deduced amino acid sequences of the S proteins detected in our study were aligned and compared with those of representative PEDV strains, including strains from GI-a (CV777 and DR13 virulent), GI-b (OH851 and CH-ZWZBa-01-2015), GII-a (CH-HNQX-3-14, CH-HNAY-2015, CH-ZMDZY-11), and GII-b (CH-GXNN-2012, CD-A). As shown in Figure 2, compared with strain CV777, the GI-a strain GDjm18-2 had three amino acid substitutions in the COE domain and one amino acid substitution in epitope SS6. The GII-a strains had amino acid substitutions at 35 positions in the COE domain, at two positions in epitope SS2, and at five positions in epitope SS6. Many new amino acid substitutions were detected in the COE regions of the GII-a strains, at positions 502 (S→ P), 507 (P→ M), 510 (N→ S), 516 (N→ D), 522 (S→ A), 527 (S→ G), 533 (A→ V), 535 (D→ E), 547 (D→ E), 559 (V→ I or A), 562 (S→ D), 567 (S→ A), 568 (K→ T or N), 570 (Q→ H), 571 (D→ N or Y), 575 (P→ L), 580 (S→ A), 588 (S→ G), 594 (T→ R or C), 608 (Y→ H), 613 (S→ I or G), 614 (G→ V), 626 (K→ E or S), and 637 (L→ F or S). Epitope 2C10 was conserved in all GII-a strains. Among the GII-a strains, GDhz16 had four continuous amino acid mutations in epitope SS6, which differed from the epitope sequence in the other strains and the reference strains. Compared with strain CV777, the GII-b strains had amino acid substitutions at 17 positions in the COE domain, and amino acid substitutions at one position in three epitopes (SS2, SS6, and 2C10). As well as the common amino acid mutations that were similar to those in the GII-b reference strains, there were novel amino acid substitutions at eight positions in the COE region: 504 (V→ L), 510 (N→ D), 535 (D→ H), 542 (S→ H), 567 (S→ Y), 614 (G→ V), 626 (K→ T), and 637 (L→ V). These amino acid sequences demonstrate that the neutralizing epitopes of the PEDV S protein are constantly mutating. This phenomenon increases the difficulty of preventing and controlling PEDV infections because existing vaccines cannot effectively protect against PEDV. However, these finding may facilitate the development of effective novel vaccines in the future.

Figure 2

Numbers of Mutated Amino Acid in Different Domains of the S Protein

To further analyze the amino acid mutations in the different domains of the S protein in these isolates, the different domains of the S protein were aligned with those of CV777, and the average number of amino acid mutations present in each year was calculated. The S protein can be divided into the S1 protein and the S2 protein. The S1 protein contains four domains: SP (amino acids 1–18), S1-NTD (amino acids 19–233), COE and RBD (amino acids 501–629), whereas the S2 protein contains five domains: SS6 (amino acids 764–771), HR1 (amino acids 978–1117), HR2 (amino acids 1274–1313), TM (amino acids 1324–1346), and 2C10 (amino acids 1368–1374). Previous data have indicated that 2C10 is conserved, so we did not analyze the 2C10 domain. As shown in Figure 3, in these strains, the S1 sequence had more amino acid mutations than the S2 sequence. From 2015 to 2018, the number of mutated amino acids in S1 remained at a high level, whereas that in S2 decreased. Furthermore, the numbers of mutated amino acids in SP (amino acids 1–18) and S1-NTD (amino acids 19–233) increased slightly, whereas the numbers of mutated amino acids in the COE and RBD domains decreased. SS6, HR1, HR2, and TM in the S2 protein did not change obviously from 2015 to 2018.

Figure 3

Pathogenicity of GDgh16

Because the samples of GDgh16 came from a scale pig farm that had experienced high mortality, and it displayed a high viral titer in Vero cells, we investigated its pathogenicity in vivo. As shown in Figure 4A, PEDV-infected cells showed a characteristic green color, indicating that PEDV was isolated successfully. The viral proliferation curve indicated that the titer of strain GDgh16 increased to 10^6.33 TCID₅₀/mL at 36 hpi but decreased to 10^4.99 TCID₅₀/mL at 96 hpi (Figure 4B). Six piglets were divided into two groups; one group was challenged orally with GDgh16 and the other group was inoculated with cell culture medium. All the challenged piglets showed classical clinical signs, including vomiting, watery diarrhea, and dehydration, at 16 hpi. The challenged piglets began to die at 24 hpi, and all had died by 48 hpi (Figure 5A). The control pigs remained healthy, with no detectable PEDV shedding. The control piglets were euthanized and necropsy was performed on all the piglets. The viral copy numbers in different parts of the intestine were determined with RT–qPCR. The duodenum, jejunum, ileum, cecum, and colon had higher viral copy numbers than the rectum (Figure 5B). The duodenums, jejunums, and ileums of the piglets were subjected to an IHC assay. As shown in Figure 5C, the tissues from the piglets in the challenged group showed remarkable levels of viral antigens compared to those in the control group. The results of GDgh16 challenge test indicate that the variant strains are a large threat to the pig industry and that the control of PEDV spread has become a critical issue.

Figure 4

Figure 5

Discussion

PEDV has become an important diarrhea virus, causing extensive damage to pig farms worldwide. Because there is no effective vaccine against the emerging prevalent strains in China, the variant PEDV strains occur frequently on many farms in different areas (24). Because there are extensive viral variants and the protection afforded by commercial vaccines is limited, it is necessary to fully understand the genetic variations and epidemiology of PEDV to facilitate the development of next-generation vaccines.

In the present study, the genetic variations in PEDV in parts of China in 2015–2018 were analyzed.

The S gene encodes the largest structural protein of PEDV and stimulates the host body to produce neutralizing antibodies against the virus. Because its variants are extensive, the S gene is commonly used as the target gene in studies of the genomic characteristics of PEDV (25). A phylogenetic analysis showed that strains from four subgroups of PEDV were present from 2015 to 2018, and that GII-a and GII-b were the two most prevalent subgroups in China at that time. From 2015 to 2018, eight strains belonging to four subgroups (GI-a, GI-b, GII-a, and GII-b) were epidemic in Jiangmen (Guangdong), which suggests that PEDV had mutated widely and the PEDV epidemic was becoming more complex. These results are consistent with those of Wen et al. (26). In 2015, all the isolated strains belonged to GII-a, whereas in 2018, 72.22% of strains belonged to GII-b, and only 22.22% of strains belonged to GII-a. Interestingly, unlike GII-a, which includes strains from other countries, such as America, South Korea, and Japan, the GII-b subgroup only contains Chinese-isolated strains. Combined with previous studies, these results suggest that GII-b strains may be the dominant strains in China in the future (27, 28).

The S protein is highly variable, and many studies have shown that amino acid changes in the S protein can affect the virulence and pathogenicity of PEDV. Our study has shown that the numbers of amino acid mutations in the SP1 and S1-NTD domains of PEDV increased in 2017 and 2018. It had been suggested that S1-NTD is a vital domain related to viral virulence (29, 30) and that conformational changes in S1-NTD are related to the high pathogenicity of PEDV strain FJzz1 (18, 27). Increasing numbers of more-virulent PEDV strains have recently emerged (18, 27, 31). Whether the mutations identified in this study alter the major conformation and thus the pathogenicity of these strains will be investigated further in the future. Our data show that the PEDV positivity rate in the provinces tested increased from 2015 to 2016, but decreased from 2016 to 2018, which might be attributable to improvements in disease prevention and control strategies. Many pig farms use the “feed-back” mode to ensure sow immunity to PEDV and to protect piglets against PEDV infection. This is an effective measure to prevent PED, but there is also a risk of virus dispersal, which is responsible for the many GI-b strains reported to date (19, 20, 32–34).

Four neutralizing epitopes of the PEDV S protein have been determined: the COE domain (499–638), epitope SS2 (748–755), epitope SS6 (764–771), and epitope 2C10 (1368–1374) (15, 16). In the present study, we detected amino acid changes at 35 positions in the COE domain. Moreover, one strain, GDhz16, had four continuous amino acid mutations in epitope SS6. Epitopes SS2 and 2C10 also contained amino acid substitutions. The antigenicity, pathogenicity, and neutralization properties of isolated strains are altered by such mutations, especially some insertions and deletions in the S protein (35, 36). Therefore, the vaccine derived from prototype strain CV777 protects against the disease induced by classical strains but not the disease caused by variant strains (24, 37). Whether these amino acid changes affect the antigenicity and neutralization properties of the four neutralizing epitopes warrants investigation in future studies.

Based on previous epidemiological and clinical observations of field strains since 2010, the emerging GII strains are highly pathogenic (38). To investigate the pathogenicity of the isolated variant strains, three piglets were infected orally with GDgh16. The piglets in the infected group began to show clinical signs of diarrhea at 12 h, and developed the typical symptoms of PED at 16 h. Morbidity reached 100%. The piglets began to die at 24 hpi, and all had died by 48 hpi. Moreover, their small intestines contained high viral copies and many viral antigens, indicating that GDgh16 was a highly pathogenic strain. Other researchers have demonstrated that different types of pigs infected with variant PEDV strains shared consistent outcomes (39–42). These results indicate that the variant strains are a large threat to the pig industry, and that the control of PEDV spread has become a critical issue.

In conclusion, the PEDV strains circulating in parts of China between 2015 and 2018 clustered into four subgroups: GI-a, GI-b, GII-a, and GII-b. The GII-b strains became dominant in 2018. Compared with previous strains, these strains displayed multiple variations in the SP and S1-NTD domains and the neutralizing epitopes of the S protein. We successfully isolated and identified a new virulent GII-b strain, GDgh16, which is well-adapted to Vero cells and causes a high mortality rate in piglets. Our study provides insight into the genetic characteristics of the prevalent PEDV strains in parts of China, and suggests that the development of effective novel vaccines is both necessary and urgent.

References

• 1.

  PensaertMBde BouckP. A new coronavirus-like particle associated with diarrhea in swine. Arch Virol. (1978) 58:243–7. 10.1007/bf01317606

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2.

  BoniottiMBPapettiALavazzaAAlboraliGSozziEChiapponiCet al. Porcine epidemic diarrhea virus and discovery of a recombinant swine enteric coronavirus, Italy. Emerg Infect Dis. (2016) 22:83–7. 10.3201/eid2201.150544

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  KimYKChoYYAnBHLimSILimJAChoIS . Molecular characterization of the spike and ORF3 genes of porcine epidemic diarrhea virus in the Philippines. Arch Virol. (2016) 161:1323–8. 10.1007/s00705-016-2758-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  PasickJBerhaneYOjkicDMaxieGEmbury-HyattCSweklaKet al. Investigation into the role of potentially contaminated feed as a source of the first-detected outbreaks of porcine epidemic diarrhea in Canada. Transbound Emerg Dis. (2014) 61:397–410. 10.1111/tbed.12269

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  SteinriglAFernandezSRStoiberFPikaloJSattlerTSchmollF. First detection, clinical presentation and phylogenetic characterization of Porcine epidemic diarrhea virus in Austria. Bmc Vet Res. (2015) 11:310. 10.1186/s12917-015-0624-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  SunRQCaiRJChenYQLiangPSChenDKSongCX. Outbreak of porcine epidemic diarrhea in suckling piglets, China. Emerg Infect Dis. (2012) 18:161–3. 10.3201/eid1801.111259

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  TianYYuZChengKLiuYHuangJXinYet al. Molecular characterization and phylogenetic analysis of new variants of the porcine epidemic diarrhea virus in Gansu, China in 2012. Viruses. (2013) 5:1991–2004. 10.3390/v5081991

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  SungMHDengMCChungYHHuangYLChangCYLanYCet al. Evolutionary characterization of the emerging porcine epidemic diarrhea virus worldwide and 2014 epidemic in Taiwan. Infect Genet Evol. (2015) 36:108–15. 10.1016/j.meegid.2015.09.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  StevensonGWHoangHSchwartzKJBurroughERSunDMadsonDet al. Emergence of Porcine epidemic diarrhea virus in the United States: clinical signs, lesions, and viral genomic sequences. J Vet Diagn Invest. (2013) 25:649–54. 10.1177/1040638713501675

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  JungKSaifLJ. Porcine epidemic diarrhea virus infection: etiology, epidemiology, pathogenesis and immunoprophylaxis. Vet J. (2015) 204:134–43. 10.1016/j.tvjl.2015.02.017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  MesquitaJRHakze-vanDHRAlmeidaALourencoMvan der PoelWHNascimentoMS. Outbreak of porcine epidemic diarrhea virus in Portugal, 2015. Transbound Emerg Dis. (2015) 62:586–8. 10.1111/tbed.12409

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  WangLByrumBZhangY. New variant of porcine epidemic diarrhea virus, United States, 2014. Emerg Infect Dis. (2014) 20:917–9. 10.3201/eid2005.140195

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  JarvisMCLamHCZhangYWangLHesseRAHauseBMet al. Genomic and evolutionary inferences between American and global strains of porcine epidemic diarrhea virus. Prev Vet Med. (2016) 123:175–84. 10.1016/j.prevetmed.2015.10.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  WangDGeXChenDLiJCaiYDengJet al. The S gene is necessary but not sufficient for the virulence of porcine epidemic diarrhea virus novel variant strain BJ2011C. J Virol. (2018) 92:e00603-18. 10.1128/JVI.00603-18

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  ChangSHBaeJLKangTJKimJChungGHLimCWet al. Identification of the epitope region capable of inducing neutralizing antibodies against the porcine epidemic diarrhea virus. Mol Cells. (2002) 14:295–9.

  • Pubmed Abstract
  • Google Scholar

• 16.

  SunDFengLShiHChenJCuiXChenHet al. Identification of two novel B cell epitopes on porcine epidemic diarrhea virus spike protein. Vet Microbiol. (2008) 131:73–81. 10.1016/j.vetmic.2008.02.022

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  ZhangQLiuXFangYZhouPWangYZhangY. Detection and phylogenetic analyses of spike genes in porcine epidemic diarrhea virus strains circulating in China in 2016-2017. Virol J. (2017) 14:194. 10.1186/s12985-017-0860-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  LiuXZhangQZhangLZhouPYangJFangYet al. A newly isolated Chinese virulent genotype GIIb porcine epidemic diarrhea virus strain: biological characteristics, pathogenicity and immune protective effects as an inactivated vaccine candidate. Virus Res. (2019) 259:18–27. 10.1016/j.virusres.2018.10.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  WangPZhuJLiuXGuoJGuXRuanW. Isolation and recombinant analysis of variants of porcine epidemic diarrhea virus strains from Beijing, China. Virus Dis. (2019) 30:294–301. 10.1007/s13337-019-00513-w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  ChenNLiSZhouRZhuMHeSYeMet al. Two novel porcine epidemic diarrhea virus (PEDV) recombinants from a natural recombinant and distinct subtypes of PEDV variants. Virus Res. (2017) 242:90–5. 10.1016/j.virusres.2017.09.013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  WangXNiuBYanHGaoDYangXChenLet al. Genetic properties of endemic Chinese porcine epidemic diarrhea virus strains isolated since 2010. Arch Virol. (2013) 158:2487–94. 10.1007/s00705-013-1767-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  KumarSStecherGLiMKnyazCTamuraK. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:1547–9. 10.1093/molbev/msy096

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  WangDFangLXiaoS. Porcine epidemic diarrhea in China. Virus Res. (2016) 226:7–13. 10.1016/j.virusres.2016.05.026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24.

  LiWLiHLiuYPanYDengFSongYet al. New variants of porcine epidemic diarrhea virus, China, 2011. Emerg Infect Dis. (2012) 18:1350–3. 10.3201/eid1808.120002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  ChenJLiuXShiDShiHZhangXLiCet al. Detection and molecular diversity of spike gene of porcine epidemic diarrhea virus in China. Viruses. (2013) 5:2601–13. 10.3390/v5102601

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26.

  WenZLiJZhangYZhouQGongLXueCet al. Genetic epidemiology of porcine epidemic diarrhoea virus circulating in China in 2012-2017 based on spike gene. Transbound Emerg Dis. (2018) 65:883–9. 10.1111/tbed.12825

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27.

  ChenPWangKHouYLiHLiXYuLet al. Genetic evolution analysis and pathogenicity assessment of porcine epidemic diarrhea virus strains circulating in part of China during 2011-2017. Infect Genet Evol. (2019) 69:153–65. 10.1016/j.meegid.2019.01.022

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  YuJChaiXChengYXingGLiaoADuLet al. Molecular characteristics of the spike gene of porcine epidemic diarrhoea virus strains in Eastern China in 2016. Virus Res. (2018) 247:47–54. 10.1016/j.virusres.2018.01.013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  HouYLinCMYokoyamaMYountBLMarthalerDDouglasALet al. Deletion of a 197-amino-acid region in the N-terminal domain of spike protein attenuates porcine epidemic diarrhea virus in piglets. J Virol. (2017) 91:e00227-17. 10.1128/JVI.00227-17

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  SuYHouYPraratMZhangYWangQ. New variants of porcine epidemic diarrhea virus with large deletions in the spike protein, identified in the United States, 2016-2017. Arch Virol. (2018) 163:2485–9. 10.1007/s00705-018-3874-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31.

  ZhangLLiuXZhangQZhouPFangYDongZet al. Biological characterization and pathogenicity of a newly isolated Chinese highly virulent genotype GIIa porcine epidemic diarrhea virus strain. Arch Virol. (2019) 164:1287–95. 10.1007/s00705-019-04167-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  LiBLiuHHeKGuoRNiYDuLet al. Complete genome sequence of a recombinant porcine epidemic diarrhea virus strain from eastern china. Genome Announc. (2013) 1:e10513. 10.1128/genomeA.00105-13

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33.

  LiKSongDZhangFGongWGuoNLiAet al. Complete genome sequence of a recombinant porcine epidemic diarrhea virus strain, CH/JXJA/2017, isolated in jiangxi, China, in 2017. Genome Announc. (2018) 6:e01590-17. 10.1128/genomeA.01590-17

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  LiRQiaoSYangYGuoJXieSZhouEet al. Genome sequencing and analysis of a novel recombinant porcine epidemic diarrhea virus strain from Henan, China. Virus Genes. (2016) 52:91–8. 10.1007/s11262-015-1254-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  ParkSKimSSongDParkB. Novel porcine epidemic diarrhea virus variant with large genomic deletion, South Korea. Emerg Infect Dis. (2014) 20:2089–92. 10.3201/eid2012.131642

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  ZhangXPanYWangDTianXSongYCaoY. Identification and pathogenicity of a variant porcine epidemic diarrhea virus field strain with reduced virulence. Virol J. (2015) 12:88. 10.1186/s12985-015-0314-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  PuranavejaSPoolpermPLertwatcharasarakulPKesdaengsakonwutSBoonsoongnernAUrairongKet al. Chinese-like strain of porcine epidemic diarrhea virus, Thailand. Emerg Infect Dis. (2009) 15:1112–5. 10.3201/eid1507.081256

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  LinCMSaifLJMarthalerDWangQ. Evolution, antigenicity and pathogenicity of global porcine epidemic diarrhea virus strains. Virus Res. (2016) 226:20–39. 10.1016/j.virusres.2016.05.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  MadsonDMArrudaPHMagstadtDRBurroughERHoangHSunDet al. Characterization of porcine epidemic diarrhea virus isolate US/Iowa/18984/2013 infection in 1-day-old cesarean-derived colostrum-deprived piglets. Vet Pathol. (2016) 53:44–52. 10.1177/0300985815591080

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  ChenQGaugerPCStafneMRThomasJTMadsonDMHuangHet al. Pathogenesis comparison between the United States porcine epidemic diarrhoea virus prototype and S-INDEL-variant strains in conventional neonatal piglets. J Gen Virol. (2016) 97:1107–21. 10.1099/jgv.0.000419

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  LinCMAnnamalaiTLiuXGaoXLuZEl-TholothMet al. Experimental infection of a US spike-insertion deletion porcine epidemic diarrhea virus in conventional nursing piglets and cross-protection to the original US PEDV infection. Vet Res. (2015) 46:134. 10.1186/s13567-015-0278-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42.

  ThomasJTChenQGaugerPCGimenez-LirolaLGSinhaAHarmonKMet al. Effect of porcine epidemic diarrhea virus infectious doses on infection outcomes in naive conventional neonatal and weaned pigs. PLoS ONE. (2015) 10:e139266. 10.1371/journal.pone.0139266

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  YuLLiuYWangSZhangLLiangPWangLet al. Molecular characteristics and pathogenicity of porcine epidemic diarrhea virus in some areas of China from 2015 to 2018. Res Square. (2020). 10.21203/rs.3.rs-41418/v1

  • CrossRef
  • Google Scholar