key: cord-0949604-cqzfvfoo authors: Nguyen Thi, Thu Hien; Chen, Chi-Chih; Chung, Wen-Bin; Chaung, Hso-Chi; Huang, Yen-Li; Cheng, Li-Ting; Ke, Guan-Ming title: Antibody Evaluation and Mutations of Antigenic Epitopes in the Spike Protein of the Porcine Epidemic Diarrhea Virus from Pig Farms with Repeated Intentional Exposure (Feedback) date: 2022-03-07 journal: Viruses DOI: 10.3390/v14030551 sha: a6310e56dd9f7fd2fbde7a063f282055ba7a879c doc_id: 949604 cord_uid: cqzfvfoo The feedback strategy, or controlled exposure of pig herd to the porcine epidemic diarrhea virus (PEDV), significantly decreased losses during a severe outbreak in late 2013 in Taiwan. However, some pig farms still suffered from recurrent outbreaks. To evaluate the association between antibody titers and clinical manifestations, sera and colostra were analyzed from one pig farm that employed the feedback strategy. Furthermore, spike (S) gene full sequences from six positive samples of two farms with and without using feedback were compared to investigate the evolution of PEDV variants circulating in pig herds. The results in this study showed that high PEDV antibody titers do not correlate with the high rate of protection from PEDV infection. In addition, repeated feedback generated the emergence of PEDV variants with unique substitutions of N537S and Y561H in the COE domain and S769F in the SS6 epitopes. These mutations indicated the pathogenetic evolution of PEDV strains existing in the cycle of the feedback method. A very strict biosecurity practice to block the routes of pathogen transfer should be followed to achieve successful control of PEDV infections in pig herds. Porcine epidemic diarrhea (PED) is a highly transmissible disease that causes substantial economic losses to pork producers all over the world [1] [2] [3] . PED clinical signs include watery diarrhea, vomiting, dehydration, high morbidity, and mortality in the suckling piglets. The causative agent, porcine epidemic diarrhea virus (PEDV), is an enveloped, single-stranded, positive-sense RNA virus, belonging to the family Coronaviridae, subfamily Coronavirinae, genus Alphacoronavirus. The genome of PEDV is approximately 28 kb in length and comprises 2 overlapping open reading frames (ORFs) 1a/1b encoding for 16 non-structural proteins (nsps), then a series of downstream ORFs encoding the 4 structural proteins, spike (S), envelope (E), membrane (M), and nucleocapsid (N), plus an accessory protein gene ORF3 [4] . The PEDV S protein is a homotrimeric membrane glycoprotein that is composed of two subunits, S1 and S2, mediating receptor binding and cell membrane fusion, respectively [5, 6] . The essential roles of S protein in host-virus The concentrations of IgG in serum or IgA in colostra were detected by using a commercial PEDV ELISA kit (GENTAUR, Kampenhout, Belgium) and AniGen PED IgA ELISA (BioNote Inc., Hwaseong-si, South Korea), respectively, and performed following the manufacturer's instructions. The neutralizing antibody titers of PEDV in sera and colostra were evaluated following a previous study [18] , with some modifications. The sera and milk of the sows were diluted in two-fold serial dilution in serum-free DMEM medium before mixing with an equal volume of 200 TCID 50 PEDV and incubated for 1 h at 37 • C, then transferred to 96-well plates containing confluent Vero cell monolayers. After 2 days, CPE was observed using an inverted microscope. Samples were identified as PEDV-positive by the conventional RT-PCR method, which targets the ORF3 gene, and confirmed again by quantitative reverse transcriptionpolymerase chain reaction (RT-qPCR), performed as previously described [9, 19] . From the positive samples, RNA of the virus was extracted by the AxyPrep Body Fluid Viral DNA/RNA Miniprep Kit (Bioscience, Unioncity, CA, USA). The cDNA was immediately synthesized with random primers by the High-Capacity cDNA Reverse Transcription Kit (ThermoFisher, Carlsbad, CA, USA). Sequences encoding the S gene (4161 bp) were amplified by PCR using GDP-HiFi DNA polymerase (GeneDireX, Las Vegas, NV, USA), with specific primers designed based on the sequences of the PEDV reference (Table 2) , following a previous study [20] under the following conditions of adjusted thermal cycles: initial denaturation at 95 • C for 5 min, 35 cycles of denaturation at 95 • C for 30 s, annealing at 60 • C for 30 s, extension at 72 • C for 4 min 30 s, and a final extension at 72 • C for 5 min. The amplified PCR products were purified by 1% agarose gel electrophoresis. The PCR product was excised from the gel then purified using FavorPrep GEL/PCR Extraction (Favorgen Biotech Corporation, Pingtung, Taiwan). The purified PCR product was then ligated into the pGEM T-Easy vector system (Promega, Fichburg, WI, USA). The identity of clones was confirmed by automated DNA sequencing (Genomics, New Taipei City, Taiwan). All data of the SN titers in serum were log-transformed as the base of 2 to show potential linear associations. Data are expressed as mean ± SEM. A Student's t-test was performed to compare treatment means between groups, and the serum antibody titers of each sow before and after the feedback were analyzed by the paired t-test. The serum antibody levels in piglets were categorized at three levels: low (L: SN ≤ 8 or ELISA IgG ≤ 25, IgA ≤ 1.5), medium (M: 8 < SN ≤ 32, 25 < ELISA IgG ≤ 40, 1.5 < ELISA IgA ≤ 2.5), and high (H: SN ≥ 32, or ELISA IgG ≥ 40, IgA ≥ 2.5). The associations of the categorical variables with the subsequent incidence of diarrhea were analyzed using the χ 2 test. Differences in which p < 0.05 were considered statistically significant. To investigate the immune responses of sows after intentional exposure to PEDV, sera were collected from 15 sows before and after feedback in farm SL. Results showed that a significant (p < 0.05) increase in antibody titers was observed either examined by the SN test (3.7 ± 0.6 and 21.5 ± 5.8 for before and after feedback) or ELISA (20.1 ± 2.7 and 27.5 ± 5.0 for before and after feedback) ( Table 3 ). Sera were also collected from two offspring of each sow seven days after farrowing. The SN titers against PEDV of piglet serum (77.2 ± 16.5) were significantly (p < 0.05) higher than those of sow serum (21.5 ± 5.8), and with a moderate positive correlation (r = 0.42) between each other. To determine the associations between antibody titers and clinical manifestations, the antibody titer distributions in pigs with and without diarrhea were compared. Sera were collected from 2 piglets 2 days after birth from each of 16 and 15 sows in farms SL and RS, respectively. Colostra were also collected from 85 sows on the day of farrowing in the SL farm. Piglets with clinical signs of diarrhea caused by PEDV were confirmed by RT-PCR. Results showed that there was no association of the serum antibody levels with diarrhea in piglets either determined by SN (χ 2 p-values of 0.15 and 0.78 for farms SL and RS, respectively) or ELISA IgG (χ 2 p-values of 0.30 and 0.39 for farms SL and RS, respectively) ( Figure 1 ). There were also no associations of the IgA antibody titers in sow colostra with the protection of piglets from PEDV infection observed with χ 2 p-values of 0.40 for farm SL (Figure 2 ). ELISA 15 20.1 ± 2.7 27.5 ± 5.0 * Data are presented as mean ± SEM. a,b Means with different superscript letters in a row are significantly different (p < 0.05). To determine the associations between antibody titers and clinical manifestations, the antibody titer distributions in pigs with and without diarrhea were compared. Sera were collected from 2 piglets 2 days after birth from each of 16 and 15 sows in farms SL and RS, respectively. Colostra were also collected from 85 sows on the day of farrowing in the SL farm. Piglets with clinical signs of diarrhea caused by PEDV were confirmed by RT-PCR. Results showed that there was no association of the serum antibody levels with diarrhea in piglets either determined by SN (χ 2 p-values of 0.15 and 0.78 for farms SL and RS, respectively) or ELISA IgG (χ 2 p-values of 0.30 and 0.39 for farms SL and RS, respectively) ( Figure 1 ). There were also no associations of the IgA antibody titers in sow colostra with the protection of piglets from PEDV infection observed with χ 2 p-values of 0.40 for farm SL (Figure 2 ). Figure 1 . Comparison of the distribution of serum neutralization and ELISA antibody titers in piglets with or without diarrhea in farms SL (a,b) and RS (c,d). Each point represents one litter and sera were collected from 2 piglets 2 days after birth from each sow. For the SL farm, n = 9 and 7 for litters with or without diarrhea, respectively, and for the RS farm, n = 10 and 5 for litters with or without diarrhea, respectively. Comparison of the distribution of serum neutralization and ELISA antibody titers in piglets with or without diarrhea in farms SL (a,b) and RS (c,d). Each point represents one litter and sera were collected from 2 piglets 2 days after birth from each sow. For the SL farm, n = 9 and 7 for litters with or without diarrhea, respectively, and for the RS farm, n = 10 and 5 for litters with or without diarrhea, respectively. The full-length S gene of PEDV strains collected from SL and CL farms was analyzed to determine their phylogenetic and genetic relationship with 34 sequences of PEDV strains from GenBank (Table A1 ). The phylogenetic tree was divided into two major groups (G1 and G2) ( Figure 3 ). Phylogenetic analysis shows that the six strains in this study cluster within the G2b group, closely related to the Japanese strain (IBR-7/JPN/2014), but less related to other Asian strains, such as Chinese PEDV strains (CH/GSTS/2016), Vietnamese PEDV strains (HUA-PED45), and Korean PEDV strains (CNUP1F-2019). However, they differ genetically from European strains (including CV777). The six strains analyzed in this study form a different sub-group from the previous virus strains isolated in Taiwan The full-length S gene of PEDV strains collected from SL and CL farms was analyzed to determine their phylogenetic and genetic relationship with 34 sequences of PEDV strains from GenBank (Table A1 ). The phylogenetic tree was divided into two major groups (G1 and G2) ( Figure 3 ). Phylogenetic analysis shows that the six strains in this study cluster within the G2b group, closely related to the Japanese strain (IBR-7/JPN/2014), but less related to other Asian strains, such as Chinese PEDV strains (CH/GSTS/2016), Vietnamese PEDV strains (HUA-PED45), and Korean PEDV strains (CNUP1F-2019). However, they differ genetically from European strains (including CV777). The six strains analyzed in this study form a different sub-group from the previous virus strains isolated in Taiwan between 2013-2014 and 2016-2018. The full length of the S gene of all five strains from farm SL is 4161 bp (1386 aa). The nucleotide sequences are separated into two closely related lineages: Group A (PT01, 02, 05) and Group B (PT 03, 04). When aligning to the CV777 vaccine strain, the identification of five strains had several point mutations, showing a low nucleotide sequence identity at 93.39% to 93.63%. In comparison with the previous Taiwan PED strains circulating between 2013-2014 and 2016-2018, the five strains collected in our study shared between 99.15% to 99.61% nucleotide sequence identities. Bootstrapping with 1000 replicates was performed to determine the percentage reliability for each internal node. Horizontal branch lengths were correlated to genetic distances among PEDV strains, with different genogroups labeled on the right. Black solid circles represent PEDV strains isolated from the SL farm that used feedback. The hollow square represents CY4 strains isolated from the CL farm with no feedback. The scale bar represents nucleotide substitutions per site. The S protein comprises a total of 1386 amino acids for all six PEDV strains obtained from SL and CL farms. In this study, we focused on analyzing the antigenic regions containing epitopes capable of inducing neutralizing antibodies against PEDV. The amino acid sequence corresponding to the COE domain and epitopes SS2, SS6, and 2C10, were located in the residues 502-641, 751-758, 767-774, and 1371-1377, respectively [21] [22] [23] . Corresponding to the phylogenetic analysis shown in Figure 3 , the five PEDV strains collected from the SL farm can be categorized into two groups when analyzed using protein sequence alignment software. Complete identities in the S protein were observed amongst the PT01, PT02, and PT05 strains (Group A) and between the PT03 and PT04 strains (Group B). A total of 13 substitutions (P27S, A169D, A174T, S182H, M214T, L299F, S769F, P885S, A1054V, P1094L, R1104K, G1225S, and E1382K) in the S protein were observed when comparing the amino acid sequences of Group B PEDV strains with those of Group A. The two groups share 99.06% sequence identity. The Group A and Group B strains exhibit 98.2% and 98.7% identity with CY4, respectively. Overall, the sequence alignment of the S proteins demonstrated that the Taiwanese topotype strain TW-Pingtung-41 (GenBank no. KP276251.1) and the six strains identified in this study share 92.22% to 93.01% of amino acid identity with the PEDV prototype CV777. Both CY4 and Group B strains exhibit 98.85% of amino acid identity with Taiwanese topotype strain TW-Pingtung-41, but Group A strain shares a slightly lower identity of 98.34%. Multiple alignments of amino acid sequences of each antigenic region in S proteins among the identified field PEDV strains in this study and those of 35 Taiwanese strains published in GenBank were conducted ( Table 4 ). The amino acid sequence of the COE domain of CY4 shares 100% identity with that of strain TW-Pingtung-41, while both Group A and Group B strains share only 98.57%. As compared with CY4 and 35 Taiwanese PEDV strains, the 5 strains from SL farms were uniquely characterized by substitution mutations of N537S and Y561H in the COE domain. Amino acid variations at residues D509A, H524P/S, K566N, L603I, Y606H, E608D/A, G612V/A/D, and E636V in the COE region were occasionally identified among those of published Taiwanese PEDV strains. Another unique substitution (S769F) in epitope motif SS6 of Group B strains but not of Group A was also observed. No mutations were identified in epitopes SS2 and 2C10. Intentional exposure of sows and gilts to PEDV via feedback of homogenized intestines could offer protective maternal immunity to the most susceptible suckling piglets and thus shorten the outbreak period and minimize economic losses on pig farms. However, determined by the method of feedback, the farm types and sizes, the herd management systems, and biosecurity measures implemented, the intentional virus-exposure practices might result in a negative or positive impact on pig farms [24, 25] . Therefore, our research at periodically used feedback farms aimed to: (i) evaluate the relationship between passive humoral immunity and the time showing clinical signs and (ii) investigate the evolution of the PEDV variants that circulated in this farm. The lactogenic immunity to protect neonatal suckling piglets via colostrum and milk is critical in the prevention and control of PEDV infection [26] . Several studies support that exposure to infected PEDV tissues can induce a significant increase in specific IgA and IgG levels [17] [18] [19] . Our data from the SL farm are consistent with these reports. Results from the SN and ELISA IgA analyses in this study showed significantly higher antibody levels in sow sera after feedback. In the SL farm, sows were given oral exposure to PEDVinfected materials. This practice might provide mucosal immunization. The previous studies showed that boosting inactivated virus in previously exposed sows may increase protection in nursing piglets against PEDV and maternally derived PEDV-neutralizing activity. The specific memory B cells could develop subsequent intentional exposure [27, 28] . Memory B cells might be critically related to the ability to induce and secrete mucosalspecific antibodies [29, 30] . In our study, PEDV-specific antibodies (IgA and IgG) were detected in colostrum and milk, implying that memory B cells responsible for IgA and IgG synthesis existed at the intestinal mucosa and can respond to this PEDV stimulation. Colostrum/milk antibody levels associated with a protective immune response in piglets may be a more reliable diagnostic tool than systemic antibody responses for detecting vaccination efficacy [27, 31] . In piglets, humoral immunity is important for protection against PED. In one recent study at a Midwestern United States sow farm, the sows were given oral immunization by feedback after an initial PEDV outbreak. Two months after this feedback approach, all weaned pigs in that farm tested negative for PEDV [32] . It has been shown that the immune response of pregnant sows vaccinated with the PEDV S1 protein provided lactogenic immunity to the piglets through colostrum and milk [33] . However, the correlation between the antibody levels and the time of showing clinical signs in piglets is not clearly described. Few researchers have considered IgA antibody responses as a correlate of protection [28, 34, 35] . The previous report has also demonstrated that anti-PEDV antibodies through intraperitoneal administration could contribute to the protection of piglets against PEDV infection and the passively transferred circulating antibodies partially improved the course of infection [36] . On the other hand, in our study, there was no association of the serum antibody levels with the diarrhea symptom either determined by the SN test or ELISA IgG. Our results also indicated no associations of the IgA antibody titers in sow colostrum with the protection of piglets from PEDV infection. Moderate or even severe PED occurred in piglets with high titers of maternal antibody in the SL pig farm, indicating that the viral loads in the barn environment might exceed the protection threshold provided by the passively transferred lactogenic immunity. Current PEDV intentional exposure is frequently ineffective, and it should be noted that: (i) there is no standard, optimized, feedback protocol, (ii) the potential of contamination of the feedback material with other infectious pathogens, (iii) mutations altering the antigenic nature of feedback viruses, and (iv) the possibility of continuous re-infection. Indeed, several studies in recent years have indicated that hazard to the herd arises from the re-infecting virus (using feedback or live vaccines) from reversion or recombination to virulence [17, 37, 38] . Therefore, the feedback protocols should only be applied in urgent explosive outbreaks to initiate and maintain a high lactogenic immunity. Disinfection and strict biosecurity measures should be thoroughly followed for a successful control and eradication of PEDV in pig herds. PEDV phylogenetic investigation in Taiwan recently confirmed that the primary PEDV strain is related to Genogroup 2b, closely related to the strains USA/Texas39/2013 and Thailand P1915-NPF-01511A [14, 20] . Although the strains isolated in this study are closely related with the US Genogroup 2b, they have greater similarity with the Japanese strain IBR7/JPN/2014 than with Chinese or earlier Taiwanese strains. Currently, very limited research has been conducted on the evolution of the PEDV circulating on the farm using feedback. While there have been multiple reports on the incidence and genetic diversity of the PED virus in Taiwan, no studies of the virus circulation in closed farms that used feedback regularly have been conducted. This research displayed several point mutations in the full spike gene of five virus strains-notably potential significant mutations of antigenic epitopes in the S gene. From our study, PEDV strains sequenced in that farm in 2020 have several mutations within the S gene compared to earlier Taiwanese strains (2013-2014 and 2016-2018) . Phylogenetic analysis categorized the five PEDV strains identified from the SL farm into two groups, showing that two variants coexist in the pig population. The coexistence of different PEDV variants within pig herds or even individual pigs with high frequency has been reported [18, 39] . The repeated feedback practices implemented in the SL farm might promote the emergence of new PEDV variants. The CY4 strain identified from the CL farm without frequent feedback exhibited 100% identity on the COE domain of S protein as compared with the topotype strain TW-Pingtung-41, but the five strains identified from the SL farm shared only 98.57%. These results revealed the fact that a complicated, constantly mutated, mixed genotype infection of PEDV might occur in pig herds utilizing a repeated feedback program. The antigenic epitopes inducing neutralizing antibodies against PEDV infection in pigs are mainly located in the S glycoprotein. The mutations therefore more frequently occur in the S protein of PEDV [40] [41] [42] . Mutations in the S protein can lead to a variety of structural and functional consequences and contribute to the lack of efficacy of PEDV vaccines [42] [43] [44] [45] . In this study, the two particular substitutions, N537S and Y561H, within the COE domain of the S protein could be significant in the pathogenetic nature of the PEDV strain isolated from the SL farm. It is interesting to note that lethal strains identified in China, South Korea, Japan, Taiwan, Canada, Mexico, and the United States all have common substitutions (-LQDGQVKI-to -SQSGQVKI-) within the SS6 epitope region of the S1 gene, which may play a role in PEDV virulence [19] . However, in our study, the unique substitution S→F (-SQSGQVKI-to -SQFGQVKI-) on the SS6 region of Group B strains could be an important mutation point to analyze the B cell epitopes on the S gene. Interestingly, no such above-mentioned mutations were found in PEDV strain CY4 identified from the CL farm without frequent feedback, and only 10 km distant from the SL farm during the same period. Furthermore, the SL farm was a tightly closed herd, and the replacement breeding stock came from only one source farm free from PEDV infection. The PEDV strains with substitution mutations would be unlikely to originate from other pig herds. The above information suggests that the amino acids at the positions of 537, 561, and 769 in the antigenic regions of S protein might be the most susceptible residues subject to substitution mutation in pig herds under immense immune pressure of the host and the high viral loads in the environments. Further examination of the molecular structure of these mutations should be carried out to understand the significance of such changes within the antigenic epitopes of PEDV in more detail. Similar to another report [46] , the sequences of SS2 (aa 751-758) and 2C10 (aa 1371-1377) motifs were quite conserved, and no substitution in these motifs was found in this study. In the case of PEDV infection, passively transferred maternal antibodies cannot provide full protection against the virus in piglets under field conditions with high viral loads. A very strict biosecurity program should be followed to achieve successful control of PEDV infections in pig herds with feedback practice. PEDV variants with unique substitution mutations of N537S and Y561H in the COE domain and S769F in SS6 epitopes emerged in pig herds under repeated intentional exposure of PEDV to sows and gilts. Informed Consent Statement: Not applicable. The data that support the findings of this study are available in the main manuscript of this article. We are grateful to Martin Ryan, School of Biology, University of St. Andrews, UK, for proofreading this paper, and our colleagues at RCAB for their excellent assistance. The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Porcine epidemic diarrhea virus: An overview of current virological and serological diagnostic methods Porcine epidemic diarrhea virus: Current insights An economic evaluation of intervention strategies for Porcine Epidemic Diarrhea (PED) Structure and Immune Recognition of the Porcine Epidemic Diarrhea Virus Spike Protein Receptor Usage and Cell Entry of Porcine Epidemic Diarrhea Coronavirus Porcine epidemic diarrhea virus S1 protein is the critical inducer of apoptosis Genomic and epidemiological characteristics provide new insights into the phylogeographical and spatiotemporal spread of porcine epidemic diarrhea virus in Asia Cellular entry of the porcine epidemic diarrhea virus Analysis of the spike, ORF3, and nucleocapsid genes of porcine epidemic diarrhea virus circulating on Thai swine farms A new coronavirus-like particle associated with diarrhea in swine Porcine epidemic diarrhea virus: An emerging and re-emerging epizootic swine virus Evolution, antigenicity and pathogenicity of global porcine epidemic diarrhea virus strains Swine enteric coronavirus disease: A review of 4 years with porcine epidemic diarrhoea virus and porcine deltacoronavirus in the United States and Canada US-like strain of porcine epidemic diarrhea virus outbreaks in Taiwan Host factors affecting generation of immunity against porcine epidemic diarrhea virus in pregnant and lactating swine and passive protection of neonates Passive immunisation, an old idea revisited: Basic principles and application to modern animal production systems Infective Material, Concepts and Procedures for Intentional Sow Herd Exposure to Porcine Epidemic Diarrhea Virus Deletion in the S1 Region of Porcine Epidemic Diarrhea Virus Reduces the Virulence and Influences the Virus-Neutralizing Activity of the Antibody Induced Molecular characterization of the spike and ORF3 genes of porcine epidemic diarrhea virus in the Philippines Phylogenetic Analysis of the Spike (S) Gene of the New Variants of Porcine Epidemic Diarrhoea Virus in Taiwan Identification of the epitope region capable of inducing neutralizing antibodies against the porcine epidemic diarrhea virus The S2 glycoprotein subunit of porcine epidemic diarrhea virus contains immunodominant neutralizing epitopes Identification of two novel B cell epitopes on porcine epidemic diarrhea virus spike protein Analysis of the effect of feedback feeding on the farm-level occurrence of porcine epidemic diarrhea in kagoshima and miyazaki prefectures, japan Successful Eradication of Porcine Epidemic Diarrhea in an Enzootically Infected Farm: A Two-Year Follow-Up Study Lactogenic immunity and vaccines for porcine epidemic diarrhea virus (PEDV): Historical and current concepts A heterologous 'prime-boost' anti-PEDV immunization for pregnant sows protects neonatal piglets through lactogenic immunity against PEDV Evaluation of humoral immune status in porcine epidemic diarrhea virus (PEDV) infected sows under field conditions Effect of Booster Vaccination with Inactivated Porcine Epidemic Diarrhea Virus on Neutralizing Antibody Response in Mammary Secretions Field trials evaluating the efficacy of porcine epidemic diarrhea vaccine, RNA (Harrisvaccine) in the Philippines Vaccination of sows with a dendritic cell-targeted porcine epidemic diarrhea virus S1 protein-based candidate vaccine reduced viral shedding but exacerbated gross pathological lesions in suckling neonatal piglets Characterization of anti-porcine epidemic diarrhea virus neutralizing activity in mammary secretions Immunogenicity and protective efficacy of recombinant S1 domain of the porcine epidemic diarrhea virus spike protein Cross protective immune responses in nursing piglets infected with a US spike-insertion deletion porcine epidemic diarrhea virus strain and challenged with an original US PEDV strain Mucosal and systemic isotype-specific antibody responses and protection in conventional pigs exposed to virulent or attenuated porcine epidemic diarrhoea virus Does circulating antibody play a role in the protection of piglets against porcine epidemic diarrhea virus? Porcine epidemic diarrhea virus (PEDV): An update on etiology, transmission, pathogenesis, and prevention and control Evaluation of responses to both oral and parenteral immunization modalities for porcine epidemic diarrhea virus in production units Novel porcine epidemic diarrhea virus (PEDV) variants with large deletions in the spike (S) ene Coexist with PEDV strains possessing an intact S gene in domestic pigs in Japan: A new disease situation Genetic characterization of the spike gene of porcine epidemic diarrhea viruses (PEDVs) circulating in Vietnam from The spike glycoprotein genes of porcine epidemic diarrhea viruses isolated in China A porcine epidemic diarrhea virus strain with distinct characteristics of four amino acid insertion in the COE region of spike protein A single V672F substitution in the spike protein of field-isolated PEDV Promotes Cell-Cell fusion and replication in Vero E6 cells Identification and pathogenicity of a variant porcine epidemic diarrhea virus field strain with reduced virulence The four of structural genes sequences of a porcine epidemic diarrhea virus from domestic piglet in Fujian Bioinformatics insight into the spike glycoprotein gene of field porcine epidemic diarrhea strains during 2011-2013 in Guangdong, China