1. Introduction
Porcine epidemic diarrhea (PED) is caused by the PED virus (PEDV), which belongs to the
Alphacoronavirus genus of the family
Coronaviridae. PEDV contains a positive-sense, single-stranded RNA genome of approximately 28 kb in size, which encodes four structural proteins, namely, spike (S), envelope (E), membrane (M), and nucleocapsid (N) proteins as well as three non-structural proteins (NSP), namely, replicases 1a and 1b and open reading frame 3 (ORF3). The first outbreak of PED was recognized in England in 1971. During the 1980s and 1990s, PED caused severe economic losses in the swine industry in Europe and Asia, but it later turned into a sporadic disease [
1]. However, since late 2010, devastating outbreaks of PED that affect all pigs, regardless of their age, with a mortality rate of up to 95% among suckling pigs, regardless of their vaccination status, have been reported in China [
2]. The disease has rapidly spread to Asia and, subsequently, to North America, and has resulted in the death of millions of pigs, and has severely affected the swine industry [
3,
4].
Compared to traditional and CV777-based vaccine strains, the high genetic diversity of the genogroup 2b (G2b) PEDV suggests that PEDV CV777- or DR13-based attenuated vaccines may not completely protect pigs against the infection or control the disease progression [
2,
5]. Therefore, a new generation of PEDV vaccine against G2b PEDV is urgently required. In the past, successful attenuation of traditional PEDV strains such as SM98 [
6], 83P-5 [
3], KPED-9 [
7], and DR13 [
8] by in vitro serial passages in Vero cells has been achieved [
9]. It has been demonstrated that conventional 11-day-old pigs orally inoculated with 2.55 × 10
5 fluorescent focus forming units of the highly-passaged cell culture strain CV777 presents transient or no fecal viral shedding without clinical symptoms [
10]. In addition, sows intramuscularly immunized with 10
7 50% tissue culture infective dose (TCID
50)/mL of the highly cell culture-adapted CV777 strain twice every two weeks induces serum and colostrum antibody titers and reduces the mortality rate from 100 to 20% in neonatal piglets born to these vaccinated sows challenged with 10 50% lethal dose (LD
50) of wild PEDV compared with those born to unvaccinated sows [
7]. By using the cell culture passage strategy, evident attenuations of a high-passage PEDV YN Chinese strain (passage 144; YN-144) in 10-day-old conventional pigs [
11] or highly-passaged PC21A-derived strains (passage 120 and passage 160) in 4-day-old, caesarian-derived, colostrum-deprived (CD/CD) or conventional suckling piglets [
12] have been established. However, the detailed immune responses or the efficacy of cross-protection induced by these attenuated G2b PEDV strains against the parental virus or heterologous viruses are either unknown or only partially effective [
12].
In Taiwan, new strains of PEDV have been detected since late 2013. Phylogenetic analysis of these field strains indicate that their
S gene sequences are highly identical to each other and are closely related to those of G2b PEDV strains [
5,
13]. In the present study, a Taiwan PEDV Pingtung 52 (PEDVPT) strain was successfully isolated from a suckling pig during a farm outbreak in 2014 and serially passaged in Vero cells. To evaluate and compare the pathogenicity of the G2b Taiwan PEDVPT field isolate with its cell culture-adapted strain, we orally inoculated 5-week-old, PEDV-seronegative conventional pigs with a low-passage (passage 5; PEDVPT-P5) or high-passage (passage 96; PEDVPT-P96) PEDV strain. After confirming complete viral clearance in feces, these pigs were then challenged or re-challenged with PEDVPT-P5 to compare the efficacy of immune protection between the groups previously inoculated with PEDVPT-P5 or PEDVPT-P96.
2. Materials and Methods
2.1. Virus Isolation and Serial Passage of PEDVPT Strain
The parental PEDVPT strain (GenBank Accession No. KP276252) was isolated in early 2014 from the intestinal homogenate of a 7-day-old suckling pig in Taiwan and adapted to Vero cells as previously described [
7,
14]. Viral infection and propagation were confirmed by daily observation of cytopathic effects (CPE), real-time reverse transcription quantitative PCR (RT-qPCR), immunofluorescence assay (IFA), and immunocytochemistry (ICC) staining. In brief, Vero C1008 cells (American Type Culture Collection (ATCC) No. CRL-1586) were cultured overnight (O/N) to 80% confluence and were washed twice with modified post-inoculation medium (PI medium) containing Dulbecco’s modified Eagle’s medium (DMEM, Gibco, Grand Island, NY, USA) supplemented with tryptose phosphate broth (0.3%), yeast extract (0.02%), and 10 μg/mL of trypsin. After two blind passages, PEDV-inoculated Vero cells showing CPE, characterized by cell fusion, syncytial cell formation, and cell detachment, were subjected to three rounds of plaque purification. PEDV-infected Vero cells showing more than 90% CPE were subjected to one freeze-and-thaw cycle, and the supernatants were harvested. Viral stocks of the low-passage PEDV (passage 5; PEDVPT-P5) strain and the high-passage PEDV (passage 96; PEDVPT-P96) strain were prepared by serial passaging of the culture supernatants in Vero cells. The virus stocks were titrated by performing a 10-fold serial dilution of the supernatants in 96-well plates in triplicates for each dilution. The viral titers of the PEDVPT-P5 and PEDVPT-P96 stocks were 10
5 and 10
6 TCID
50/mL, respectively.
2.2. Detection of PEDV Antigens by IFA and ICC
Supernatants of the virus-inoculated Vero cells showing typical CPE were eluted, 200 µL of 80% ice-cold acetone (Sigma-Aldrich, St. Louis, MO, USA) was added to each well, and the cells were fixed at −20 °C for 10 min. After fixation, acetone was removed, and the plate was air dried at room temperature for 30 min. The air-dried cells were then incubated with an in-house anti-PEDV N protein monoclonal antibody diluted 1000-fold (for ICC stain) or with PEDV hyperimmune pig serum diluted 200-fold (for IFA) and incubated at room temperature for 1 h. After incubation, the wells were washed three times with phosphate-buffered saline (PBS). For IFA, 100 μL of fluorescein isothiocyanate (FITC)-labeled goat anti-mouse immunoglobulin G (IgG) (Bethyl Laboratories, Montgomery, TX, USA) diluted to 200-fold in PBS was added to each well and incubated in dark at room temperature for 1 h. For ICC, a polyclonal anti-rabbit/mouse immunoglobulin EnVision-DAB+ system (Dako, Carpinteria, CA, USA) was used according to the manufacturer’s protocol. After incubation with the secondary antibody, each well was washed three times with PBS. For ICC staining, the labeled cells were exposed to 3,3′-diaminobenzidine (DAB) chromogen using a peroxidase DAB substrate kit (Dako, Carpinteria, CA, USA), according to the manufacturer’s instructions. For IFA, the labeled cells in the 96-well plate were mounted with 1% glycerol (Sigma-Aldrich) in PBS. The stained cells were detected using an inverted light/fluorescence microscope.
2.3. Comparative Pathogenicity and Assessment of Immunogenicity in PEDVPT-P5 and PEDVPT-P96-Inoculated 5-Week-Old Conventional Pigs
Eighteen 5-week-old, PEDV-seronegative, PEDV-fecal RNA-negative, Large White × Duroc crossbred pigs were selected from a conventional pig farm. Selected pigs were negative for porcine circovirus type II, porcine respiratory and reproductive virus, porcine rotavirus, transmissible gastroenteritis virus, porcine respiratory coronavirus, and porcine deltacoronavirus as tested by routine conventional or real-time PCR detection. These pigs were randomly assigned to four groups, viz., the mock-inoculated group with PEDVPT-P5 challenge (pigs A1 to A3; n = 3), mock-inoculated group without PEDVPT-P5 challenge (pigs A4 to A6; n = 3), the PEDVPT-P5-inoculated group with PEDVPT-P5 challenge (pigs B1 to B6; n = 6), and the PEDVPT-P96-inoculated group with PEDVPT-P5 challenge (pigs C1 to C6; n = 6). All pigs were labeled with ear tags. Three pigs of the same group were housed in a separate room. Each pig in the PEDVPT-P5- and PEDVPT-P96-inoculated groups was orally inoculated with 5 mL of 105 TCID50/mL of the respective virus diluted in PI medium. Each pig in the mock-inoculated group received 5 mL of PI medium orally. Rectal swabs were collected every day to monitor the duration of viral shedding, and clinical signs were recorded and scored daily. To monitor seroconversion and mucosal immunity, 10 mL of ethylenediaminetetraacetic acid (EDTA)-anticoagulated blood and fecal swabs were collected from inoculated pigs every two weeks for detecting PEDV-specific plasma IgG and mucosal IgA, respectively. All experiments procedures performed on the animals were reviewed and approved by the Institutional Animal Care and Use Committee of National Taiwan University (Taipei, Taiwan, NTU105-EL-00087).
2.4. Re-Challenging with PEDVPT-P5 to Assess the Protection in PEDVPT-P5- and PEDVPT-P96-Inoculated Pigs
After confirming that all PEDVPT-P5- and PEDVPT-P96-inoculated pigs were seroconverted and negative for fecal virus shedding by RT-qPCR, all pigs were orally inoculated with 5 mL of 105 TCID50/mL of PEDVPT-P5 to evaluate immune protection. At this time point, all pigs were 9-weeks old.
2.5. Scoring for Clinical Signs of Infection
After each inoculation, animals were monitored daily for clinical signs of infection until the end of the experiment. PEDV-associated diarrhea was scored based on the fecal consistency as follows: 0, normal; 1, loose; 2, semi-fluid; and 3, watery [
15]. The weight of individual pigs in each group was monitored and recorded weekly.
2.6. RNA Extraction, cDNA Synthesis, and SYBR Green-Based Quantitative Real-Time PCR
Fecal swabs were diluted with 700 μL of PBS, and RNA was extracted using a QIAamp Viral RNA Mini Kit (Qiagen, Chatsworth, CA, USA), according to the manufacturer’s instructions. The total amount of feces contained in each fecal swab ranged from 0.25 to 0.35 g [
15,
16]. Complementary DNA (cDNA) synthesis was performed by reverse transcription using the QuantiTect Reverse Transcription Kit (Qiagen), according to the manufacturer’s instructions. The synthetized cDNA was amplified by real-time qPCR using SYBR Advantage qPCR Premix (Clontech, Palo Alto, CA, USA) under the following conditions: 45 cycles of 95 °C for 5 s, followed by 60 °C for 30 s. The primers used for detection of PEDV in the fecal samples have been previously published [
17]. All real-time qPCR reactions were performed in duplicates, and the results were expressed as genomic equivalents (GE), as previously stated [
17]. Ten-fold serial dilutions of a known amount of plasmid (pCR-XL-TOPO DNA; Thermo Fisher Scientific, Waltham, MA, USA) containing the PEDV
N gene was used as the positive control to generate a standard curve. According to this standard curve, the detection limit of RT-qPCR was approximately 60 GE of DNA (data not shown). The amplification product for each reaction was also confirmed by performing the melting curve assay. Samples with a melting temperature of 84.7–86.1 °C were considered PEDV positive. The efficiency of the qPCR ranged from 90.42 to 96.92%.
2.7. Antibody Responses
For concentration and purification of PEDV virions, the cell component of the supernatants collected from PEDV-infected Vero cells was first freeze–thawed three times, and the cell debris was eluted by centrifugation at 6000× g for 30 min. To concentrate the virions, viral supernatants were pelleted by centrifugation at 75,000× g for 2.5 h using an Avanti J-25 centrifuge (Beckman, Fullerton, CA, USA) and were resuspended in PBS. Viral protein concentrations were determined using a Pierce BCA (bicinchoninic acid) Protein Assay Kit (Thermo Fisher Scientific Inc.). Ninety-six-well, flat-bottom plates (Corning Life Sciences, Corning, NY, USA) were coated O/N with 2 μg/well of concentrated PEDV virions. The PEDV-coated plates were washed six times with 100 mL of PBST (PBS containing 0.05% Tween 20), and wells were incubated for 1 h at room temperature with plasma samples diluted to 20-fold in PBS or with syringe-filtered (0.22 μm pore size) fecal samples, which was diluted two-fold in PBS; each fecal swab was diluted in 700 μL PBS. Each sample was added to the plate in duplicates. After incubation with samples, the wells were washed six times with PBST, and antibodies were detected using horseradish peroxidase (HRP) conjugated goat anti-pig IgG (Kirkegaard & Perry Laboratories, Gaithersburg, MD, USA) diluted to 1:1000 or goat-anti-pig IgA (Abcam, Cambridge, UK) diluted to 1:5000. After washing six times with PBST, 100 µL of tetramethylbenzidine (TMB) substrate solution (Kirkegaard & Perry Laboratories) was added to each well and incubated at room temperature for 20 min. The reaction was stopped using TMB stop solution (Kirkegaard & Perry Laboratories), and the optical density (OD) at 405 nm was recorded on an enzyme-linked immunosorbent assay (ELISA) reader. Antibody titers were expressed as sample-to-positive control ratio (S/P ratio) values. Using a panel of positive and negative controls, the cutoff values of plasma IgG and fecal IgA detection were 0.1 (95% confidence level) and 0.3 (95% confidence level), respectively.
2.8. Neutralizing Antibody Assay
For neutralization assay, 100 µL of Vero cell suspension at a density of 3 × 105/mL were seeded onto 96-well culture plates (Corning Life Sciences) and incubated at 37 °C, 5% CO2 overnight to reach 70% confluence. Prior to performing neutralizing antibody assay, plasma samples were heated in water bath at 56 °C for 30 min to inactivate the complement and diluted at 1:10, 1:20, 1:40, 1:80, and 1:160 in PI medium. Fifty microliters of diluted plasma samples were admixed and incubated with 200 TCID50/mL PEDVPT-P5 at 37 °C and 5% CO2 for 1 h. After removing the medium from the 70% confluent Vero cells and washing the cells twice with PI medium, mixtures containing PEDVPT-P5 and different diluted plasma were added to Vero cells in duplicates and were incubated at 37 °C, 5% CO2 for 1 h. After 1 h of incubation, the supernatants were replaced by fresh PI medium. The cells were maintained at 37 °C and 5% CO2, and the CPE were monitored for the next three days. The neutralizing titer was determined as the last dilution without CPE.
2.9. Complete Viral Genome Sequencing and Full-Length Sequence Analysis
Viral RNA extraction and subsequent reverse transcription of both PEDVPT-P5 and PEDVPT-P96 were performed as previously described [
5] using sequence-specific reverse primers (see
Supplementary Table S1). The complete genomes of both viruses were obtained by using seven pairs of oligonucleotide primers that flank different regions of the novel strains of PEDV genome with at least 200-base pair overhangs (see
Supplementary Table S1). The PCR amplicons of each fragment were gel-purified, cloned into plasmid vector (pJET1.2; Thermo Fisher Scientific), and subjected for sequencing after plasmid isolation. To acquire sequences of the extreme 3′ and 5′ ends, SMART RACE cDNA Amplification Kit was used according to the manufacture’s instruction (Clontech, Tokyo, Japan). Sequence assembly and comparison of complete genomes of both PEDVPT-P5 and PEDVPT-P96 were performed by using Lasergene 7.1 (DNASTAR, Madison, WI, USA) and MEGA7 [
18] , respectively.
2.10. Statistical Analysis
The results of body weights, antibody titers, and fecal viral shedding were analyzed statistically using SAS 9.4 (Statistical Analysis System, SAS Institute Inc., Cary, NC, USA). The variables among groups were compared using one way analysis of variance (ANOVA) combined with Scheffe’s method. A p-value of <0.05 was interpreted as statistically significant.
4. Discussion
In the present study, we successfully isolated and characterized the pathogenicity of a new Taiwan PEDVPT strain, PEDVPT-P5. We have demonstrated that this G2b PEDVPT-P5 strain can induce typical watery diarrhea in 5-week-old pigs but milder symptoms and a shorter viral shedding period in 9-week-old pigs. By serial passaging of PEDVPT in Vero cells, a PEDVPT-P96 viral stock was also prepared. Compared to PEDVPT-P5-inoculated pigs, PEDVPT-P96-inoculated pigs showed a delayed, mild, and transient diarrhea, or in some cases undetectable, providing evidence of viral attenuation. Of importance, inoculation of both PEDVPT-P5 and PEDVPT-P96 induced PEDV-specific plasma IgG, mucosal IgA, and neutralizing antibodies and elicited protection against a subsequent challenge with the low passage PEDVPT-P5 strain. This study provides valuable information regarding the viral shedding period, immunogenicity, and pathogenicity of different passaged PEDV in conventional post-weaning pigs. Moreover, the evidence of viral attenuation and the efficacy of protection induced by PEDVPT-P96 in conventional pigs highlights PEDVPT-P96 as a potential, live-attenuated PEDV vaccine candidate.
It has been shown that PEDV can infect seronegative pigs of all ages. However, the susceptibility to PEDV infection is age-dependent and is significantly higher in neonatal pigs than in weaning pigs [
20]. In the present study, we demonstrated that compared to 5-week-old pigs, which exhibited PEDV-associated clinical signs for 7–9 days and long viral shedding period of 26 days, older 9-week-old pigs exhibited milder clinical symptoms for a shorter period lasting of 4–7 days and viral shedding lasting for only nine days. Furthermore, we demonstrated that even without displaying clinical signs, the viral shedding periods could be longer than three weeks in PEDV-infected post-weaning pigs. It has been reported that oral fluid samples are suitable for screening PEDV in commercial growing pig herds [
21]. These findings suggest that routine screening for oral and fecal PEDV shedding in clinically healthy pigs when purchasing feeder pigs from swine producers is necessary to reduce the risk of disease transmission from farm to farm.
In the present study, we demonstrated that 5-week-old conventional pigs inoculated with the PEDVPT-P5 strain exhibited fecal viral shedding, peaked at 2–4 DPI, and persisted for 26 days. An increase in PEDV-specific IgG titer in these pigs started at 14 DPI and reached a plateau at 21 DPI. The result is similar to a previous finding in 4-week-old conventional pigs challenged with 1.42 × 10
5 TCID
50 of a USA/KS/2013 PEDV isolate. In this study, USA/KS/2013 isolate-inoculated pigs displayed fecal viral shedding primarily during the first two weeks with a peak at 5–6 DPI and fecal swabs that were positive for PEDV nucleic acid at 21 and 28 DPI in some pigs. Seroconversion was also noted at 14 DPI and reached a peak level at 21 DPI [
22]. When the spike sequence of the USA/KS/2013 isolate was compared to the Taiwan PEDVPT isolate, we found that these two viruses share high sequence identity (99%).
In Asia, the attenuated CV777-based and DR13-based PEDV vaccines have been used for decades and, as a consequence, PED in Asia has been generally under control with low prevalence and mortality rates in the past [
9]. These live-attenuated vaccines are derived from cell culture-adapted PEDV strains and have proven to be effective in controlling PED by reducing the severity and duration of diarrhea, shortening the period of viral shedding, stimulating the immune response of sows, and protecting the delivered piglets through elevated maternal antibodies [
7,
23,
24]. Unfortunately, the exposure to traditional PEDV strains or the above-mentioned live-attenuated vaccine strains provide limited cross-protection against the epidemics of the novel virulent PEDV strains [
15]. Therefore, a novel vaccine against these new strains is urgently needed. In the present study, we have developed a high passage PEDVPT-P96 viral stock by serial passaging in Vero cells and have provided evidence of its viral attenuation. In addition, PEDVPT-P96 inoculation of pigs induced plasma IgG, mucosal IgA, and neutralizing antibodies and provided complete immunization against a subsequent low passage PEDVPT-P5 challenge, suggesting that PEDVPT-P96 may serve as a potential live-attenuated vaccine candidate.
By comparing the sequences of the parental virulent PEDVPT-P5 and attenuated PEDVPT-P96 strains to identify potential genetic determinants for viral attenuation, we demonstrated that most non-silent mutations are located in the
S gene. On comparing the result of sequence analysis with other previous studies, two mutations (S887R and C1354F) identified in the attenuated late passages of YN1 strain were also identified in our PEDVPT-P96 strain. However, the deletion at the position of aa 144 of the
S gene, suggested by the same study as a potential attenuation marker for PEDV [
10] was not found in PEDVPT-P96; alternatively, an aa substitution (T144I) occurred in this particular region. The clinical significance of these mutations requires further elucidation. Of note, one unique aa substitution (F554S) was observed in the COE neutralizing domain in PEDVPT-P96, but the mechanism of how it interferes with antibody recognition remains to be elucidated. However, the results of in vivo animal challenge performed in this study suggest that this mutation does not result in the loss of the cross-protection against the parental PEDVPT-P5 strain. In the present study, aa substitutions of F669S in the
NSP3 gene and M252I in the
NSP15 gene were identified in PEDVPT-P96. The NSP of coronaviruses are well known to be responsible for viral replication, and many of them are reported to be interferon (IFN) antagonists. In PEDV, NSP1, NSP3, NSP7, NSP14, NSP15, and NSP16 have been demonstrated to suppress the IFN-β and IRF3 promoter activities [
23]. The effects of mutations in
NSP3 and
NSP15 genes of PEDVPT-P96 on immune modulation need further investigation. Unlike other reports wherein the
NSP3 and
ORF3 genes tend to comprise more mutations through PEDV attenuation, only one non-silent mutation was observed in each gene in the present study. To date, many studies have been conducted aiming to map attenuating mutation(s) [
10,
11,
24]; however, there is little consistency among different studies, indicating that different strains may harbor different mutation patterns, thereby suggesting that the attenuation may be multifactorial and introduction of attenuating mutation(s) may have implications in generating a safe live attenuated vaccine. This information is invaluable for further quick and rational development of effective and genetically stable vaccines via reverse genetics.
It has been reported that the gut–mammary secretory IgA axis is the most promising and effective pathway against enteric diseases, including PED. For PEDV vaccine development, it is crucial for neonatal piglets to receive sufficient maternal antibodies against PEDV via colostrum or milk from previously immunized sows [
8,
25]. Furthermore, the amount and the duration of maternal antibodies is closely related to the antibody titers of the immunized sows [
9]. Immune protection against PEDV infection, including pathogenicity, immunogenicity, and fecal viral shedding, in nursing pigs born to PEDVPT-P96-inoculated sows as well as the potential for reversal of virulence of the PEDVPT-P96 strain by serial animal passages will be further evaluated.