Porcine epidemic diarrhea virus (PEDV) is a member of genus Alphacoronavirus
within the family Coronaviridae
and is a single-stranded positive-sense RNA virus. PEDV was firstly identified in Belgium in 1978 [1
]. The disease caused by PEDV was firstly observed in feeder pigs and fattening swine in England [2
] and was proposed the name “porcine epidemic diarrhea” (PED) [3
]. PED is an acute and highly contagious enteric disease characterized by severe enteritis, vomiting and watery diarrhea in swine.
In China, PED was firstly occurred in Shanghai in 1973, however the causative agent was PEDV confirmed in 1984 [4
]. Since first identified there have been occasional outbreaks of PEDV in China. However, outbreaks of PEDV have been observed on most swine breeding farms in most provinces since late 2010 in China. The piglets within one week of age and sometimes within only a few hours showed vomiting, severe watery diarrhea, severe dehydration and high mortality (up to 100%), but few sows or fattening swine showed any clinical signs during the outbreaks [5
]. The economic loss caused by PEDV infection has been continuous and serious in China.
The Spike (S) gene is considered the most useful in revealing the genetic diversity of coronavirus isolates [7
]. Like other coronaviruses, PEDV S gene also plays an important role in the molecular epidemiology and in the genetic variation of PEDV field strains [10
]. In order to control and prevent PEDV infection, it is necessary for us to further investigate the prevalence of PEDV and the molecular characteristics of the S genes of Chinese PEDV field strains during 2011–2012. In this study, we firstly detected PEDV in the samples collected from 177 farms, then analyzed the sequence characteristics of S genes of field strains and performed the phylogenetic analysis of the field strains.
Since December 2010, outbreaks of viral diarrhea have been observed on most swine-breeding farms in most of the provinces of China. The clinical appearance was similar to that of PEDV infection. In this study, 79.66% (141 of 177) pig farms in 29 provinces (excluding Tibet and Hainan, China) were confirmed as infected with PEDV, and 72.27% (417 of 577) clinical samples were positive for PEDV. The results indicated that PEDV has a high prevalence rate in pig herds. TGEV was also detected by RT-PCR [16
] in our lab, and the overall detection rate was 29.12% in clinical samples. Two other research groups in China reported that 40.8% (20 of 49) and 22.73% (5 of 22) sow milk samples were positive for PEDV in 2011 [5
], and hypothesized that sow milk could represent a possible route for the vertical transmission of PEDV from sow to suckling piglets [5
]. We also detected PEDV in 76.71% (56 of 73) sow milk samples in 2012. These results indicated that newborn piglets could be infected by PEDV in sow milk.
Our previously study indicated that the S gene of CH/FJND-3/2011) has the lowest sequence homology (93.75%) with that of CH/S and is 9 nt longer than those of CV777 and CH/S [17
]. A more precise study with the full-length S genes of recently Chinese field strains will help to best understand the genetic relationships among prevailing PEDV in China. In this study, the entire S genes of 33 representative Chinese field strains during 2011–2012 were determined and compared with those of CV777-attenuated, CH/S and CHGD-1.
The S genes of attenuated strains (CV777-attenuated, DR13-attenuated, 83p-5-100) are 4149 nt in length [18
], however those of classical virulent strains are 4152 nt in length, such as CV777 [20
], Br1/87 [21
] and Chinese classical strain CH/S [22
]. The S genes of 33 recently field strains show the diversity in length. The S genes of 13 recently field strains in G1-1 are shorter than those of classical strains in length because of the nucleotides deletion, however, the S genes of 20 recently field strains in G2 are longer than those of classical strains in length because of the nucleotides insertion and deletion.
The unique characteristic of 20 field strains in G2 in the present study was the presence of the nucleotides insertion and deletion in the S genes. The alignment of S genes reveals that one domain which locates within the 5' most 1100 nt of the S gene exhibited an increased divergence compared to the remaining part of the sequence. The insertion and deletion regions are mainly clustered in this region. Furthermore, the largest number of nucleotide differences is also clustered in this domain. Most of the diversity between the recent field strains in G1-1 and in G2 is located within the 1100 nt-domain. The phenomenon is also found in the S gene of CH/FJND-3/2011 [17
Our findings showed the S genes of 13 field strains in G1-1 exhibited higher homologies with those of CH/S and CV777-attenuated than with that of CHGD-1. However, the S genes of 20 field strains in G2 exhibited higher sequence identities with that of CHGD-1 than with those of CH/S and CV777-attenuated. These results further indicated the heterogeneity in S genes of PEDV field strains prevailing in China.
According to the phylogenetic analyses, the recently Chinese field strains were divided into two distinct clusters. 13 field strains and reference strains including Korean field strains, Japanese strains, PEDV-attenuated strains and Chinese classical strains were determined to belong to G1, whereas 20 field strains and 2 reference strains (CH/FJND-3/2011, CHGD-1) were clustered together within G2. These data demonstrate that prevalent PEDV field strains in China have two types (classical strains and variant strains), and most field strains are variant strains. These variants are more closely related to each other and to other previously identified Chinese variants CH/FJND-3/2011 and CHGD-1 rather than Chinese classical strains (CH/S, DX, JS-2004-2, LJB/03, LZC) and attenuated strains.
The antigen epitope motif (1368
] and the COE domain (aa499-638) [14
] of S protein can induce neutralizing antibodies against PEDV. According to the alignment of amino acid sequences of S proteins, the antigen epitope motif (1368
) is conserved in all the recent field strains. The COE domain consists of 140 aa. Compared to CV777-attenuated, 27 field strains have 1-8 different aa, 7 field strains have no different aa. A previously study reported that three of these mutations in the recent isolates influenced the hydrophobicity of the S protein as compared with that for CV777 [5
4.1. Clinical Samples Detection
A total of 577 samples (504 small intestines of dead neonatal piglets and 73 sow milk) were collected from 177 farms in 29 provinces from February 2011 to November 2012. The collection procedures of milk were as follows. The sows’ teats were washed three times with sterile physiological saline. The workers milked the teats with sterile gloves, and took the middle section of the milk in 1.5mL sterile microcentrifuge tubes. The milk samples were sent to our lab with ice bags or dry ice. All the samples were initially detected by reverse transcription-polymerase chain reaction (RT-PCR) using previously described methods [23
]. The animal experiments were approved by Harbin Veterinary Research Institute and performed in accordance with animal ethics guidelines and approved protocols. The approval number of Animal Ethics Committee is Heilongjiang-SYXK-2006-032.
4.2. Sequencing of S Gene
To obtain the full-length S gene, four primers (S1-1U, 5'- ATCGTCAGAGGCATTTTTAA-3'; S1-1L, 5'-ATCCATCACCATTAAACGAA-3'; and S1-2U, 5'- ATGTTGTGTTAGGCTTGTTG-3'; S1-2L, 5'- CACTAACAGGCGTGTTGTAA-3') and three primers (S2-U, 5'- CTGATTCTGGACAGTTGTTA-3'; S2-1L, 5'- TTGGACAGCATCCAAAGACA-3'; S2-2L, 5'- CTTCGAGACATCTTTGACAA-3') were designed and synthesized according to the corresponding region of CV777 (AF353511) to amplify S1 domain and S2 domain by PCR or nested-PCR, respectively. The complementary DNA (cDNA) of 33 representative field strains was synthesized using S2-2L as the reverse transcription primer according to the previously described methods [23
PCR was carried out in a two-step reaction, first with a pair of primers flanking the region to be amplified and then using a pair of primers within the amplified sequence in a total volume of 25 µL. Exactly, 1µL (10 ng) cDNA was mixed with a reaction mixture containing 12.5 µL 2× EmeraldAmpTM PCR Master Mix (TaKaRa, Dalian, China), 0.5 µL each specific primer (10 µM), and 10.5 µL sterile deionized water. The first round of amplification was performed under the reaction conditions (pre-denaturation at 94 °C for 5 min followed by 30 cycles of denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, extension at 72 °C for 3 min, and a final extension at 72 °C for 7 min). The second round of amplification was performed: a 100-fold dilution in distilled water of the first PCR products were used as the template under the reaction conditions of the first round of amplification (pre-denaturation at 94 °C for 5 min followed by 25 cycles of denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, extension at 72 °C for 2.5 min, and a final extension at 72 °C for 7 min).
PCR products were excised from 1.0% agarose gels, extracted using QIAquick® Gel Extraction Kit (QIAGEN GmbH, Hilden, German), cloned into a T-tailed vector, pMD18-T (TaKaRa, Dalian, China) and transformed using JM109 competent cells (TaKaRa, Dalian, China). Three to five recombinant DNA clones of each PEDV strains were sequenced by Beijing Genomics Institute (Beijing, China).
4.3. Multiple Alignments and Phylogenetic Analysis
The sequences of full-length S genes of 33 field strains (Table 1
) and 30 reference strains (Table 4
) were used in sequence alignments and phylogenetic analysis. The multiple sequence alignments, nucleotide and amino acid sequences divergences were analyzed by Jotun Hein method in the software program MegAlign [24
]. Phylogenetic tree was constructed from the aligned nucleotide sequences by the maximum likelihood method in MEGA 5.05 program [25
Reference strains used in this study
Reference strains used in this study
|Reference strains||Countries||S genes (nt)||S proteins (aa)||Accession no.|
|DR13 attenuated ||South Korea||4149||1382||DQ462404|