key: cord-0701605-6zz0o89w authors: Zhao, Pengwei; Wang, Song; Chen, Zhi; Yu, Jiang; Tang, Rongzhi; Qiu, Wenbin; Zhao, Lu; Liu, Yueyue; Guo, Xiaozhen; He, Hongbin; Xu, Guanlong; Li, Jinxiang; Wu, Jiaqiang title: Successive Passage In Vitro Led to Lower Virulence and Higher Titer of A Variant Porcine Epidemic Diarrhea Virus date: 2020-04-01 journal: Viruses DOI: 10.3390/v12040391 sha: dc61db3f0571d3a11a3451ebd58997cada9d988a doc_id: 701605 cord_uid: 6zz0o89w A highly virulent porcine epidemic diarrhea virus (PEDV) appeared in China and spread rapidly to neighbor countries, which have led to great economic losses to the pig industry. In the present study, we isolated a PEDV using Vero cells and serially propagated 100 passages. PEDV SDSX16 was characterized in vitro and in vivo. The viral titers increased to 10(7.6) TCID(50)/mL (100th) by serial passages. The spike (S) gene and the whole gene of the SDSX16 virus was fully sequenced to assess the genetic stability and relatedness to previously identified PEDV. Along with successive passage in vitro, there were 18 nucleotides (nt) deletion occurred in the spike (S) gene resulting in a deletion of six amino acids when the SDSX16 strain was passaged to the 64th generation, and this deletion was stable until the P100. However, the ORF1a/b, M, N, E, and ORF3 genes had only a few point mutations in amino acids and no deletions. According to growth kinetics experiments, the SDSX16 deletion strain significantly enhanced its replication in Vero cells since it was passaged to the 64th generation. The animal studies showed that PEDV SDSX16-P10 caused more severe diarrhea and vomiting, fecal shedding, and acute atrophic enteritis than SDSX16-P75, indicating that SDSX16-P10 is enteropathogenic in the natural host, and the pathogenicity of SDSX16 decreased with successive passage in vitro. However, SDSX16-P10 was found to cause lower levels of cytokine expression than SDSX16-P75 using real-time PCR and flow cytometry, such as IL1β, IL6, IFN-β, TNF-α, indicating that SDSX16-P10 might inhibit the expression of cytokines. Our data indicated that successive passage in vitro resulted in virulent attenuation in vivo of the PEDV variant strain SDSX16. Porcine epidemic diarrhea (PED) is a highly contagious swine disease. The main clinical signs include severe vomiting, diarrhea, and dehydration in suckling piglets [1] [2] [3] . The disease was first Intestine samples and faeces were collected from piglets with diarrhea in a pig farm during a large-scale outbreak in Shandong province, China in 2016. The samples were PEDV-positive, but negative to transmissible gastroenteritis virus (TGEV), porcine rotavirus (PoRV), and porcine delta virus (PDCoV) was by using multi-RT-PCR. Small-intestine tissues were used to generate a 10% (w/v) homogenate in Earle's balanced salt solution (Sigma-Aldrich, St. Louis, MO, USA). The suspension was centrifuged at 4000× g for 10 min at 4 • C. The supernatant was filtered through a 0.22 µm-pore-size filter and used for virus isolation. A one-tenth gram of feces was suspended in 1 mL phosphate-buffered saline (PBS), vortexed for 5 min, and then centrifuged at 4000× g for 5 min. The supernatant was filtered through a 0.22 µm-pore size filter and used as an inoculum for virus isolation. The virus was isolated in Vero cells (ATCC CCL-81), as described previously with some modifications [27, 28] . Vero cells were seeded at 70% confluency in 6-well plates (Corning, NY, USA) and incubated overnight. After washing three times with PBS (pH = 7.4), cells were inoculated with 400 µL inoculum diluted in 200 µL of maintenance medium with trypsin (MMT) containing DMEM, 0.3% tryptose phosphate broth (TPB, Sigma-Aldrich, St. Louis, MO, USA), 3 µg/mL of trypsin and 1% penicillin/streptomycin, were added to each well [7, 28, 29] . After incubation 2 h at 37 • C, 5% CO 2 , the inoculated cells were maintained at 37 • C under 5% CO 2 and monitored daily for cytopathic effects (CPE). When 80% of CPE appeared, cells were subjected to three rounds of freeze and thaw. The culture supernatants were then collected and centrifuged at 500× g for 10 min. The supernatants were aliquoted and stored at −80 • C until use. This was designated as 'passage 1 (P1)' viral stock, then confluent Vero E6 cells in T25 flasks were infected by PEDV strain at an MOI of 0.01. When extensive CPE was observed, the cell culture supernatant was harvested and used for the next passage in Vero E6 cells. To conduct growth kinetics experiments, Vero cells cultured in 6-well-plates were infected with the SDSX16 isolate from different passages at an MOI of 0.01 and incubated at 37 • C under 5% CO 2 . The culture supernatants were harvested from at different time points (6, 12, 18, 24 and 48 h) and stored at −80 • C. Virus titers were measured in 96-well plates by 10-fold serial dilutions of samples, in triplicate per dilution. The Vero E6 cells were seeded in a 96-well plate at a density of 1 × 10 5 cells/well in 200 µL and incubated for 24 h at 37 • C in 5% CO 2 . Ten-fold serial dilutions of the virus were added to the prepared wells with 100 µL of diluted virus per well. The tissue culture with 50% infective dose (TCID 50 ) was expressed as the reciprocal of the highest dilution showing CPE by the Reed-Muench method [30] . The Vero cells infected by PEDV (MOI = 0.01) in 24-well plates for 48 h were fixed with 4% paraformaldehyde for 15 min at room temperature (RT). Cells were permeabilized with 0.5% Triton X-100 in PBS at 4 • C for 15 min and were blocked with 1% bovine serum albumin (BSA) in PBS for 30 min at RT. The monoclonal antibody specific against PEDV spike protein was added to the cells as primary antibodies and incubated for 1 h at 37 • C. After washing three times in PBS, cells were incubated for 1 h at RT with dilution 1:1000 of the FITC-conjugated donkey anti-mouse IgG antibody (ThermoFisher, Waltham, MA, USA), followed by counter-staining with 4, 6-diamidino-2-phenylindole (DAPI, Sigma-Aldrich, St. Louis, MO, USA). The coverslips were mounted on microscope glass slides in mounting buffer and staining was visualized using a fluorescence Leica DM IL LED microscope (Leica, Wetzlar, Germany). ImageJ software version 1.51p was used to calculate the percentage of CPE cells per 200 cells in the visual field. Vero cells were seeded in 6-well plates and incubated overnight until the cell density reached 80% confluency. After washing with PBS, cells were inoculated with SDSX16-P5 at a multiplicity of infection (MOI = 1) diluted in MMT. After 2 h of incubation at 37 • C, cells were washed with PBS and then cultured with MMT. Cells were observed daily to check for any cytopathic effect (CPE), particularly syncytium formation. When CPE was obvious, the supernatant was discarded and the cells were collected by gently blowing with a cell scraper. The cells were collected in sterile tubes by centrifugation at 3000× g for 10 min and then supernatants were discarded. Cell pellets were washed with PBS and fresh 3% glutaraldehyde fixative was added and incubated for 3 h at 4 • C. The conventional TEM sample preparation method was followed for rinsing, 1% osmium acid (OsO4) fixed, drift-washed, dehydrated, soaked, Epon 812 embedded, semi-thin sliced with LKB-V-type ultra-thin Slicer line ultra-thin slices, lead citrate and uranium acetate dyes, and examined by JEOL-1200EX electron microscopy. Images were recorded using MORADA-G2. Virus RNA was extracted from the supernatants of homogenized samples using Trizol reagent according to the manufacturer protocols (ThermoFisher, Waltham, MA, USA). Complementary DNA (cDNA) was produced through reverse transcription, using PrimeScript™ One-Step RT-PCR Kit Ver.2 according to the manufacturing protocols (TaKaRa Bio, Otsu, Shiga, Japan). Targeted sequence primers of the PEDV S gene and the whole gene were amplified by PCR, with reference to the PEDV CV777 strain (GenBank: AF353511.1). The S1 fragment was amplified using primers PEDV-S1F: 5 -GTGGCTTTTCTAATCATTTGGTC-3 and PEDV-S1R: 5 -CTGGGTGAGTAATTGTTTACA-3 , under the following conditions: reverse transcription at 50 • C for 30 min, denaturation at 95 • C for 2 min, 35 cycles of 94 • C for 30 s, 50 • C for 30 s, and 72 • C for 3 min, and 72 • C for 10 min. The product fragment was 2508 bp [31] . The S2 fragment was amplified using primers PEDV-S2F: 5 -GGCCAAGTCAAGATTGCACC-3 and PEDV-S2R: 5 -AGCTCCAACTCTTGGACAGC-3 , under the following conditions: reverse transcription at 50 • C for 30 min, denaturation at 95 • C for 2 min, 35 cycles of 94 • C for 30 s, 52 • C for 30 s, and 72 • C for 2 min, and 72 • C for 10 min. The product fragment was 1850 bp. The whole genome of PEDV was amplified using primers in Table 1 , the reaction conditions are as follows: Reverse transcription at 50 • C for 30 min, denaturation at 95 • C for 2 min, 35 cycles of 94 • C for 30 s, 52 • C for 30 s, and 72 • C for 1 min 30 s, and 72 • C for 10 min. The RT-PCR products were identified by electrophoresis on a 1% agarose gel, and cloned into the pMD19/T vector (TaKaRa Bio, Otsu, Shiga, Japan). The positive clones were sequenced at the Shanghai Sangon Biological Engineering Technology & Services Co., Ltd. (China). Sequence assembly was carried out using Lasergene Software version 7.1. Moreover, sequence alignments and phylogenetic trees of the S genes and the whole SDSX16 genomes were also performed using distance-based neighbor-joining method of MEGA4.0 software. The numbers at each branch represent bootstrap values greater than 50% of 1000 replicates. The three-day-old Duroc crossbred piglets were obtained from the Shandong Academy of Agricultural Science (Jinan, Shandong, China). All piglets were confirmed to be free of PEDV, TGEV, PoRV, and PDCoV infections using a viral serum neutralizing test for the detection of antibodies against PEDV, TGEV, PoRV, and PDCoV, as well as RT-PCR to detect viral nucleic acids. To determine the virulence and pathogenesis of the PEDV isolate and successive passage strains, SDSX16 variants at five passage levels (SDSX16-P10, -P30, -P63, -P64, -P75) were selected for pig study. Thirty three-day-old Duroc crossbred piglets were randomly divided into six groups (P10, P30, P63, P64, P75, and Mock), each group contained five piglets. The infected group (P10, P30, P63, P64, and P75) were orally challenged with P10, P30, P63, P64, and P75, respectively (10 6 TCID 50 /piglet), while the uninfected group (Mock) was orally given maintenance medium in the same route. These piglets were observed and recorded daily for signs of vomiting, diarrhea, and body weight until the end of the experiment and humanely euthanized following the appearance of severe diarrhea or after becoming moribund. At the end of the study, all surviving piglets were euthanized humanely. At necropsy, intestinal tissue specimens (<0.3 cm thick) were collected from each piglet and were fixed with 10% formalin for 24 h at RT and embedded in paraffin according to the standard laboratory procedures. The formalin-fixed paraffin-embedded tissues were cut in 5 µm thick on a microtome (Leica), floated on a 40 • C water bath containing distilled water, and transferred onto glass slides. The tissues were then deparaffinized in xylene for 5 min and washed in decreasing concentrations of ethanol (100%, 95%, 90%, 80%, and 70%) for 3 min each. The fixed intestine sections were evaluated for PEDV antigen by immunohistochemistry (IHC) using 200× diluted PEDV-specific monoclonal antibodies (1:200) with the antigen retrieval method described previously [32] . This study was approved by the local animal ethics committee. Porcine intestinal epithelial cells were cultured in 6-well plates at a concentration of 10 6 cells per well, were infected with SDSX16-P10, -P63, -P64, -P75 (MOI = 1) respectively. The porcine intestinal epithelial cells in 6-well plates were transfected with poly(I:C) (5 µg/mL) (Sigma-Aldrich, St. Louis, MO, USA) were served as a positive control group by using Lipofectamine 3000 transfection reagent (ThermoFisher, Waltham, MA, USA). Hours post-infection, the cells were collected for real-time qPCR and Flow cytometry to detect cytokines. To detect mRNA levels of cytokines, total RNA was extracted with Trizol reagent according to the manufacturer protocols (ThermoFisher, Waltham, MA, USA) and subjected to reverse transcription with PrimeScript RT Master Mix (TaKaRa Bio, Otsu, Shiga, Japan). Measurements were performed using the SYBR Premix Ex Taq (TaKaRa Bio, Otsu, Shiga, Japan ) on a LightCycler 480 II with a 96-well plate block (Roche, Switzerland). Amplification was performed for 30 s at 95 • C, followed by 40 cycles of 95 • C for 5 s, 60 • C for 30 s. The relative expression of mRNA was quantified within each sample using the 2 −∆∆Ct method, and the cut-off Ct value was 35. β-actin mRNA expression levels were used as a housekeeping gene for internal control. Primers for each cytokine and gene were listed in Table 2 . Statistically significant differences were determined by Student's t-test using GraphPad Prism software (version 8.0; San Diego, CA, USA). Two-way ANOVA was used to calculate statistical significance, and multiple samples were compared by one-way ANOVA. Values are expressing as mean ± SEM of at least independently repeated three times unless otherwise noted, and a p-value below 0.05 was considered statistically significant (*p < 0.05, **p < 0.01, ***p < 0.001). We tried to isolate PEDV from RT-PCR-positive clinical samples in Vero cells. One PEDV isolate designated as SDSX16 was successfully isolated from the feces of a naturally infected piglet. This sample was originated from a commercial swine farm located in Shandong Province and was obtained on September 16, 2016. This SDSX16 virus and it could produce distinct CPEs from passages 5 (P5), typical of PEDV infection, such as cell fusion, syncytia formation, and detachment in infected Vero cells. We then investigated whether the isolate could be efficiently cultivated and maintained in cultured cells. To do this, the isolated PEDV strain SDSX16 was serially passaged in Vero cells for a total of 100 passages (P1 to P100). The time for the onset of CPE was typically 24 h and prominent CPE was observed within 48 h post-inoculation. The development of CPE was faster in later passages. The visible CPE appeared at 18 h and became predominant by 36 h. Virus propagation was confirmed by detecting PEDV antigens by IFA using an anti-PEDV spike-specific monoclonal antibody. The distinct staining pattern for typical syncytia formation was distributed in the cytoplasm. In contrast, neither CPE nor spike-specific monoclonal antibody staining was detected in mock-inoculated cells. Examples of CPE, and corresponding IFA images in cells inoculated with passages SDSX16-P10 and -P64, are shown in Figure 1 . The titer of the SDSX16 was measured every ten generations. The virus titers ranged from 10 4.15 to 10 4.50 TCID 50 /mL up to SDSX16-P5, and the peak viral titer reached 10 7.6 TCID 50 /mL after passage 100. The titer of SDSX16-P70 (7.44 log 10 TCID 50 /mL) was approximately 10-fold higher than that of SDSX16-P60 (6.22 log 10 TCID 50 /mL). By retrospective testing, we found that the titer of SDSX16-P64 (7.22 log 10 TCID 50 /mL) was significantly higher than that of SDSX16-P63 (Figure 2A ). The growth kinetics study further showed that SDSX16 replicated rapidly and efficiently in Vero cells, reaching a maximum titer >10 7.0 TCID 50 /mL by 48 h of inoculation, and the titer of SDSX16-P64 began to be significantly higher than the titer of SDSX16-P63 by 18 h of inoculation ( Figure 2B ). The morphology of PEDV SDSX16-P5 was characterized in Vero cells by electron microscopy. The virus-induced intracellular membrane structures were captured with a bright domain in the cytoplasm. The presence of clustered spherical particles in loose membrane vesicles was found in the cytoplasm of Vero cells infected by PEDV SDSX16-P5, showing a typical shape of coronaviruses (Figure 3) . Sequencing of the spike gene of the SDSX16 isolate (GenBank number: KY679906.1) yielded a sequence of 4161 bp, which was analyzed through a comparison with 23 typical strains available in NCBI. The phylogenetic tree ( Figure 4A ) revealed that these PEDV strains could be divided into Subgroup I and Subgroup II. Amino acid analysis indicated that the SDSX16 isolate shared 93.5%~94.9%, 97.5%~98.4%, and 98.5%~99% identity to Subgroup I, Subgroup IIa, and Subgroup IIb, respectively. The major variable region of the S gene was the S1 domain. Compared with Subgroup I (CV777 and DR13) and Subgroup IIb (SDSX16, SD2014, and PC22A-P10 et al.) [14, 18, 19, 33] isolates exhibited an insertion of five amino acids at 59-62 position, 140 position, as well as the deletion of two amino acid at 163-164 position ( Figure 4B ). The same insertions and deletions have been found in other variant strains isolated, indicating that epidemics are caused by variant strains in many countries. Moreover, in comparison with the CV777 strain, some amino acid substitutions in the neutralizing region of the core of epitope (COE) (499-638aa) were observed [34, 35] and the regions of the aligned sequences, corresponding to these regions were COE (514-642aa, Figure 4C ). The sequence analysis revealed that the other variant strains had the same COE mutations. The S gene of SDSX16 at different passage levels (SDSX16 P10-P100 GenBank number: MN241452-MN241464) were analyzed. The results showed that the spike gene had an 18 nt deletion resulting in the deletion of 6 amino acids when the SDSX16 was passaged to the P64, and this deletion was stable until the P100 ( Figure 4D Supplementary Table S1 . Interestingly, the amino acid sequence analysis results of SDSX16-P63, -P64, -P75 showed no mutations in ORF1a/b, M, N, E, and ORF3 genes. Also, the S gene of SDSX16-P64 mainly lost six amino acids than the -P63 S gene and the S gene of -P75 had no amino acid deletion or addition compared to the S gene of -P64. These results suggest that the biggest difference in the genomes of SDSX16-P64 and -P63 is the 6 amino acid deletion in -P64. To determine the virulence and pathogenesis of the PEDV isolate and different passage strains, we have conducted oral inoculation of nursery pigs with PEDV SDSX16-P10, -P30, -P63, -P64, -P75, respectively. We found that piglets (group A) orally challenged with of SDSX16-P10 (10 6 TCID 50 /piglet), developed diarrhea at 18 h, furthermore, the infected piglets presented with anorexia and depression, and the piglets began to exhibit watery diarrhea at 2 days post-challenge ( Figure 5A ). Two piglets died at 2 days post-challenge and all piglets died at 5 days post-challenge in group A ( Figure 5B) . Whereas, the piglets (group B) started to vomit 24 h after orally inoculated with the -P30 and showed diarrhea at 30 h and some animal's clinical signs of vomiting after suckling. However, there was no significant difference in the clinical examination between groups C and D, infected by -P63, and -P64, respectively, and piglets were found to have slight diarrhea at 48 h. The piglets' group E challenged with P75 virus was not showing any significant diarrhea symptom. No, any clinical signs were observed in group F. In this animal experiment study, there was a significant decrease in weight post-challenge in the P10, P30 infected group, whereas the uninfected piglets gained weight ( Figure 5C ). No obvious temperature changes were observed between these groups. Immunohistochemistry of PEDV confirmed the presence of the virus in the cytoplasm of epithelial cells on atrophic villi in some segments of the small intestine ( Figure 6A-E) . PEDV antigen was detected in some cells of ileum, jejunum, duodenum, and cecum of piglets infected by SDSX16-P10, P30, P64, P65 or P75. SDSX16-P10 and -P30 caused severe villous atrophy in the cytoplasm of villous enterocytes in most jejunum and ileum. P64 and P65 caused moderate villous atrophy in some jejunum and ileum. P75 caused mild villous atrophy in a few jejunum or ileum. PEDV IHC antigen was not detected in the intestines of mock piglets ( Figure 6F ). The piglets infected with SDSX16-10 developed severe acute enteritis typical of PEDV throughout the study, demonstrating that the isolate was virulent and enteropathogenic in neonatal piglets. The piglets were infected by the higher passage of SDSX16, which had slighter diarrhea, vomiting and fecal shedding. The results indicated that the mutations in the S gene, caused by the serial passages of SDSX16, might be related to the attenuation of virulence in vivo. Pro-inflammatory cytokines and type I interferons are produced at the porcine intestinal epithelial cells as part of the innate immune response during the infection process. We evaluated the changes in gene expression of pro-inflammatory cytokines (IL1β, IL6, and TNF-α) and type I interferon (IFN-β) in porcine intestinal epithelial cells infected with different SDSX16 strains. The mRNA levels of IL1β, IL6, TNF-α and IFN-β in porcine intestinal epithelial cells infected with SDSX16-P10 strain were much lower than SDSX16-P75 at 18 h post-infection (hpi). However, there was no significant difference in the examination of cytokine expression between porcine intestinal epithelial cells infected by P63 and P64 respectively (Figure 7 ). Next, we used flow cytometry to detect IL6 and TNF-α protein levels in cells infected with SDSX16-P10 or -P75. The results showed that SDSX16-P10 strain-stimulated IL6 and TNF-α expression of the protein was markedly lower than SDSX16-P75 at 18 hpi ( Figure 8 ). Collectively, these data suggest that there was no significant difference in cytokines expression between cells infected with SDSX16-P63 and -P64, and SDSX16-P10 had a stronger ability to inhibit cytokines expression than SDSX16-P75. The PEDV variant strains have been detected in many provinces in China and have become one of the most economically important viral causes for diarrhea in piglets [36, 37] . At present, piglet diarrhea, caused by PEDV infection, leads to huge losses in many pig farms [38] . We collected clinical samples from pig breeding farms in Shandong Province where severe acute diarrhea was associated with high mortality in suckling piglets. Isolation of variant strains is important for effective control of piglet PEDV infection and vaccine development. In this study, we successfully isolated a strain of PEDV from the intestinal tissues of a PED-infected piglet. The PEDV SDSX16 isolate showed obvious CPE (syncytium formation) in Vero cells from after passage 5, the viral titers started from roughly 10 4.15 TCID 50 /mL and increased after serial passages to 10 7.6 TCID 50 /mL. Phylogenetic analysis based on the S gene sequences showed the SDSX16 isolate belonged to subgroup II. More recently, it has been reported that PEDV G2 strain-based inactivated vaccine candidates are candidates for developing an effective vaccine against the virulent pandemic PEDV strains [39] . These results suggest the SDSX16 strain could be a candidate strain for vaccine development against PEDV variant strains. To develop the PEDV vaccine, we attenuate the virulence of the SDSX16 strain by the successive passage in vitro. Considering that S protein plays a crucial role in mediating virus-cell fusion and inducing neutralizing antibodies, we focus on analyzing to vary the S gene during the successive passage. The result found the S gene had greater diversity, while other genes (ORF1a/b, M, N, E, and ORF3) occurred acid point mutants. We further determined the serial passage of SDX16 in Vero cells resulted in the viral growth adaptation and promoting its replication in vitro by growth kinetics and virus titration assay. The serial passage of PEDV and other coronaviruses often results in the attenuation of virus virulence [40] [41] [42] [43] [44] . To study whether the serial passages of SDSX16 have influenced the viral pathogenicity, PEDV-seronegative, three-day-old pigs were inoculated orally with SDSX16 different viruses. The results suggested that successive passage resulted in attenuation of the SDSX16 strain, and deleting six amino acids in the S gene might affect the pathogenicity of the virus. Although data indicating pivotal roles of S protein in determining phenotypes of coronaviruses have been accumulated, it remains to be determined whether the mutations in the S gene of SDSX16 contribute to the viral attenuation and virus virulence. Therefore, we will need to apply PEDV reverse genetics systems to further studies to define the genetic mechanisms that underlie the attenuation of PEDV. It has been widely reported that type I interferons (IFN-α/β) are two essential cytokines that can control viral infection [45] . TNF-α has an antiviral function, and IL6 is associated with improving humoral and mucosal immune response [21, 46] . Variant PEDV (genogroup G2) exhibits higher resistance to the interferon than the existing PEDV and avoids the neutralizing antibody against PEDV [47] [48] [49] . In order to understand whether there were differences in the cytokine levels between the SDSX16-P10 and -P75 strain. We compared SDSX16-P10 with -P75 and found that IFN-β mRNA levels could be up-regulated in porcine intestinal epithelial cells. However, the mRNA and protein levels of pro-inflammatory cytokines, in cells induced by SDSX16-P75, were much higher than SDSX16-P10. Additionally, it has been reported that highly pathogenic PEDV strain suppressed the induction of pro-inflammatory cytokines and type 1 interferon production through the down-regulation of Toll-like receptors (TLRs) [25, 26, [50] [51] [52] . SDSX16-P75 induces higher levels of cytokine production indicating that the pathogenicity of the SDSX16 strain is weakened by the successive passage in vitro, which is consistent with the results of animal experiments. Our study demonstrated that the ability of PEDV to induce type I IFNs is strain-specific. Thus, Type I IFNs may play an important role in PEDV replication and pathogenesis. In summary, this study explored the virulent PEDV SDSX16 strain that was isolated in Shandong Province, China. The mutants of the S gene maybe improve the virus replication capability in Vero cells it was passaged to the 64th generation and was not related to the virus virulence. SDSX16-P75 induced higher cytokines levels in porcine intestinal epithelial cells and highly attenuated compared to -P10. SDSX16-P75 may be a suitable candidate virus to develop a variant strain-based PEDV vaccine to control current PEDV outbreaks. Phylogenetic analysis of porcine epidemic diarrhea virus field strains prevailing recently in China Porcine epidemic diarrhea virus variants with high pathogenicity Epidemic strain YC2014 of porcine epidemic diarrhea virus could provide piglets against homologous challenge Chinese-like strain of porcine epidemic diarrhea virus New variants of porcine epidemic diarrhea virus, China Porcine Epidemic Diarrhea Virus among Farmed Pigs Isolation and molecular characterization of porcine epidemic diarrhea viruses collected in Japan in 2014 Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States The first case of porcine epidemic diarrhea in Canada Distinct characteristics and complex evolution of PEDV strains Emergence of Porcine epidemic diarrhea virus in the United States: Clinical signs, lesions, and viral genomic sequences Complete Genome Sequence of Porcine Epidemic Diarrhea Virus Strain USA/Colorado/2013 from the United States Porcine epidemic diarrhoea virus: A comprehensive review of molecular epidemiology, diagnosis, and vaccines Detection and molecular diversity of spike gene of porcine epidemic diarrhea virus in China Immunogenicity and protective efficacy of recombinant S1 domain of the porcine epidemic diarrhea virus spike protein Immunogenicity of a recombinant parapoxvirus expressing the spike protein of Porcine epidemic diarrhea virus Detection of antibodies against porcine epidemic diarrhea virus in serum and colostrum by indirect ELISA Molecular characterization and phylogenetic analysis of membrane protein genes of porcine epidemic diarrhea virus isolates in China Molecular characterization and phylogenetic analysis of new variants of the porcine epidemic diarrhea virus in Gansu Molecular characterization of pig epidemic diarrhoea viruses isolated in Japan from Tumour necrosis factors α and β inhibit virus replication and synergize with interferons Tumor necrosis factor-alpha and interleukin-1 beta play a critical role in the resistance against lethal herpes simplex virus encephalitis Pattern recognition receptors and inflammation Interferons and viruses: An interplay between induction, signalling, antiviral responses and virus countermeasures Suppression of type I interferon production by porcine epidemic diarrhea virus and degradation of CREB-binding protein by nsp1 Inhibition of NF-kappaB activity by the porcine epidemic diarrhea virus nonstructural protein 1 for innate immune evasion Quantitation, biological and physicochemical properties of cell culture-adapted porcine epidemic diarrhea coronavirus (PEDV) Isolation and characterization of a Korean porcine epidemic diarrhea virus strain KNU-141112 Pathogenesis of porcine epidemic diarrhea virus isolate (US/Iowa/18984/2013) in 3-week-old weaned pigs A simple method of estimating fifty per cent endpoints12 Isolation and characterization of porcine epidemic diarrhea viruses associated with the 2013 disease outbreak among swine in the United States Pathology of US porcine epidemic diarrhea virus strain PC21A in gnotobiotic pigs Attenuation of an original US porcine epidemic diarrhea virus strain PC22A via serial cell culture passage Identification of the epitope region capable of inducing neutralizing antibodies against the porcine epidemic diarrhea virus The GPRLQPY motif located at the carboxy-terminal of the spike protein induces antibodies that neutralize Porcine epidemic diarrhea virus A molecular epidemiological investigation of PEDV in China: Characterization of co-infection and genetic diversity of S1-based genes Genetic evolution analysis and pathogenicity assessment of porcine epidemic diarrhea virus strains circulating in part of China during Outbreak of porcine epidemic diarrhea in suckling piglets Evaluation and comparison of immunogenicity and cross-protective efficacy of two inactivated cell culture-derived GIIaand GIIb-genotype porcine epidemic diarrhea virus vaccines in suckling piglets The group-specific murine coronavirus genes are not essential, but their deletion, by reverse genetics, is attenuating in the natural host attenuated coronavirus vaccines through the directed deletion of group-specific genes provide protection against feline infectious peritonitis Mutations in the spike gene of porcine epidemic diarrhea virus associated with growth adaptation in vitro and attenuation of virulence in vivo Complete genome sequences of a Korean virulent porcine epidemic diarrhea virus and its attenuated counterpart Attenuation of replication by a 29 nucleotide deletion in SARS-coronavirus acquired during the early stages of human-to-human transmission Identification and characterization of MAVS, a mitochondrial antiviral signaling protein that activates NF-kappaB and IRF 3 Evaluation on the efficacy and immunogenicity of recombinant DNA plasmids expressing spike genes from porcine transmissible gastroenteritis virus and porcine epidemic diarrhea virus Resistance of field-isolated porcine epidemic diarrhea virus to interferon and neutralizing antibody Vaccines for porcine epidemic diarrhea virus and other swine coronaviruses Generation and protective efficacy of a cold-adapted attenuated genotype 2b porcine epidemic diarrhea virus The Interplay Between Host Innate Immune Responses and Porcine Enteric Coronaviruses Porcine epidemic diarrhea virus inhibits dsRNA-induced interferon-beta production in porcine intestinal epithelial cells by blockade of the RIG-I-mediated pathway Porcine epidemic diarrhea virus nucleocapsid protein antagonizes beta interferon production by sequestering the interaction between IRF3 and TBK1 We particularly thank the animal house staff for regularly monitoring and preparing the piglets for the study. The authors declare no conflict of interest.