Porcine epidemic diarrhea virus (PEDV) is an enveloped, single-stranded, positive-sense RNA virus in the genus Alphacoronavirinae of the family Coronaviridae. PED, which is caused by PEDV, is characterized by severe diarrhea, dehydration, and high mortality rates in the affected swine. PEDV variants belonging to genogroup 2 (G2) have been emerging in China in a large-scale outbreak characterized by approximately 80–100 % morbidity rates and high mortality rates among suckling piglets since late 2010. Highly virulent PEDV strains, phylogenetically related to the G2 Chinese PEDV variants (shared ≥99.5 % nt identity), suddenly emerged in the United States in May 2013 and rapidly spread throughout the country, causing severe economic losses [3, 4]. The second PEDV variant in the US, designated S INDEL PEDV, with insertions and deletions in the N terminal region (S1 subunit) of the spike (S) protein same as the G1 PEDV, was identified subsequently. The S INDEL strains caused reportedly milder disease in the field, indicating that the S gene contains major virulent determinants. The other large S-deletion PEDV strains were also identified, including one Korean field strain with a 204-aa deletion, and a US cell-adapted strain with a 197-aa deletion. Whether these S-deletion PEDV strains have distinct pathogenic characteristics, have not been described.
In July of 2013, fecal samples from diarrheal sows, collected in a pig farm in Jiangsu province of east China, were submitted to our lab for routine laboratory diagnostics. PEDV positivity was confirmed by reverse transcription PCR (RT-PCR) according to the method reported previously. One of the PEDV-positive samples, from a sow with very mild clinical sign, was used to inoculate into Vero cells to isolate the virus, as described previously. This PEDV strain, designated FL2013 strain, was successfully passaged and propagated, as characterized by typical PEDV-induced cytopathic effects (CPE), such as cell fusion and syncytia formation (data not shown), which was similar to what we observed for the CHGD-01 strains. The other field PEDV samples collected from pigs with severe clinical sign could not adapt to Vero cells. The complete genomic cDNA of the FL2013 strain was further determined using the 3rd cell culture passaging virus by amplification of twelve regions covering the PEDV genome as described previously. The sequences were assembled and analyzed using DNASTAR program and MEGA6.0 program.
The FL2013 PEDV genomic sequence had the size of 28,044 nt excluding the polyadenosine tail (GenBank accession no. KP765609). Twenty available PEDV sequences with complete genomes were used for multiple alignment and phylogenetic analysis (Table 1). The phylogenetic trees based upon either the complete genome (Fig. 1a) or S gene (Fig. 1b) showed that the FL2013 strain belongs to the genogroup G2b as proposed by Huang et al.. The S gene harbors two significant insertions at aa 58 to 61 (QGVN) and 142 (N), and a deletion of two aa (DI) between aa positions 167 and 168 at the N-terminus in comparison with the prototype CV777 strain (Fig. 2), which is the unique sequence signature of the virulent PEDV strains in the genogroup G2 [2–4]. However, the extreme C-terminus of the FL2013 S gene has a unique 21-nt deletion, leading to a 7-aa deletion (FEKVHVQ) in comparison with the other G2 PEDV sequences (Fig. 2). Interestingly, this unique deletion was also found in a G1 Korean strain SM98 (GenBank accession no. GU937797; Fig. 2). The complete S gene was sequenced in the original fecal sample, the 5th, 10th and 20th cell culture passages, respectively, and no sequence alterations were found, indicating that the 7-aa deletion was naturally present in the PEDV field strain and was stable in cell culture passaging. In addition, we did not find significant mutations on the other structural genes, non-structural genes and untranslated regions (UTR) potentially associated with attenuation.
Since the FL2013 strain was isolated from a sow with very mild clinical sign, we were interested to determine its virulence by experimental infection in newborn piglets in comparison with the virulent G2 strain. Thirty-six 3-day-old piglets, negative for PEDV RNA were assigned into three groups with 12 in each. Piglets in each group were housed with their mothers (PEDV RNA and antibody negative) with no artificial supply of colostrum and milk. Group A serves as negative control, whereas piglets in groups B and C were challenged orally with the FL2013 strain and the virulent CHGD-01 strain at 1.0 × 105 50 % tissue culture infectious doses (TCID50)/3 mL, respectively. Two piglets from each group were euthanized and necropsied at 3 days post-inoculation (dpi) for histopathologic examinations and immunohistochemistry staining according to Stevenson et al.. The remaining 10 piglets in each group were examined monitored for clinical signs and recorded to evaluate survival rate for 25 dpi. The result showed that Group-B piglets, together with the negative control group-A piglets did not show clinical sign, and all 10 remaining piglets in each group survived during the observation period (Fig. 3). Group-C piglets had PED signs characterized by acute vomiting and watery diarrhea, and immunohistochemistry staining showed abundant viral antigen in the severely atrophic villi (Fig. 3). In contrast, less PEDV positive staining was found in the intestinal villi of group-B samples, and the villi was only slightly damaged (Fig. 3). 10 of 10 remaining piglets in group C died within 7 dpi. The comparative pathology study indicated that the FL2013 strain had reduced virulence to newborn piglets.
In summary, we identified an apparently attenuated G2 PEDV field strain with 7-aa deletion at the 3′-end of the S gene. Genetic characterizations of various PEDV field strains are of importance for demonstration of PEDV outbreaks and development of efficacy vaccines. Findings from present study reveals that variant PEDV strains circulating in swine-producing areas in China, and provides new molecular epidemiological data for a better understanding of this disease. It is highly possible that the sequence deletions (especially the 7-aa deletion), and other unknown mutations found in the variant strain FL2013 might have contributed to the reduced severity of the clinical disease in the piglets. Whether the changes in the genomic sequence, insertions and deletion, could collectively alter the efficiency of viral replication and RNA synthesis in the FL2013 strain leading to reduced pathogenicity will be investigated in further. More studies will be conducted to test this hypothesis. In addition, the FL2013 isolate may serve as a potential vaccine candidate that could protect the Chinese piglets from the infection caused by the virulent strain of PEDV.