key: cord-0697333-dm7rx4ls authors: Wang, Xiao-meng; Niu, Bei-bei; Yan, He; Gao, Dong-sheng; Yang, Xia; Chen, Lu; Chang, Hong-tao; Zhao, Jun; Wang, Chuan-qing title: Genetic properties of endemic Chinese porcine epidemic diarrhea virus strains isolated since 2010 date: 2013-06-25 journal: Arch Virol DOI: 10.1007/s00705-013-1767-7 sha: 9e7a8a7cddeaff40a45a8a7d3965d34d52b8f617 doc_id: 697333 cord_uid: dm7rx4ls Acute diarrhea outbreaks caused by porcine epidemic diarrhea virus (PEDV) have been observed in various pig-breeding provinces of China since December 2010. Endemic strains of PEDV were isolated from different areas, and the complete genome sequences of 10 isolates were determined. Our objective in this study was to genetically characterize current Chinese field isolates of PEDV to better understand their epidemiology and genetic diversity. Sequence analysis showed that 10 post-2010 isolates shared high homology with each other and were always clustered together with the virulent DR13 strains (South Korea) and/or one earlier Chinese strain, CH-S, in phylogenetic analysis. All post-2010 isolates possessed common sequence changes in each gene. Our results suggest that current Chinese PEDV isolates originated from either South Korean and/or Chinese ancestors that underwent some genetic variation, thereby forming a new PEDV genotype in China. Porcine epidemic diarrhea virus (PEDV) is a positive-sense single-stranded RNA virus belonging to group 1 of the genus Alphacoronavirus, within the family Coronaviridae. The virus is the etiological agent of porcine epidemic diarrhea (PED), an acute and highly contagious enteric disease characterized by watery diarrhea, vomiting, dehydration, and high mortality in nursery piglets [1] . PEDV was first recognized as the causative agent of PED in Belgium [2] and the United Kingdom in 1978 [3] . The presence of PED in China was confirmed, by immunofluorescence assays and serum neutralization tests, in 1984. Bivalent killed or attenuated vaccines have been used to prevent transmissible gastroenteritis (TGE) and PED in China since then. Clinical PED outbreaks have reemerged in the major pigproducing areas of China since December 2010, resulting in major losses of suckling piglets. Morbidity rates between 90 and100 % have been observed, with mortalities of 70-100 % on affected pig farms [4] . Heterogeneity of genomic sequences is believed to be responsible for the diversity in PEDV pathogenicity. A number of molecular investigations have been conducted, revealing different levels of variation in nucleotide sequences among PEDV isolates [5] [6] [7] . Because of incomplete genomic information, most molecular epidemiological analyses have targeted a single gene of PEDV. To date, complete genome sequences of 18 PEDV strains, including 10 endemic Chinese isolates, have been submitted to the GenBank database [4, [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] . To better understand the diversity, epidemiology and molecular characteristics of PEDV isolates in China, we sequenced the complete genome of a virulent PEDV strain from Henan Province, China. This isolate was from a suckling piglet that exhibited severe diarrhea. We compared its genomic sequence with those of other Chinese PEDV isolates and reference strains. Isolation and identification of PEDV strain CH/ZMDZY/11 PEDV strain CH/ZMDZY/11 was isolated from 2-to 3-day-old suckling piglets suffering from watery diarrhea and dehydration. These piglets were from a single farm in Henan Province, China, which contained 1500 breeding sows. Cases were characterized by watery diarrhea, dehydration with milk curd vomitus, thin-walled intestines, 100 % morbidity and 90 % mortality. Pigs of other ages were affected, exhibiting diarrhea and loss of appetite to different degrees, with severity of symptoms dependent on age. Affected sows were characterized by anorexia, depression and transient diarrhea, with most affected sows recovering within one week. Production records showed that breeding herds had been immunized with a killed CV777-based vaccine that had been produced in China. Virus isolation was conducted using Vero E6 cells (American Type Culture Collection, VA, USA). Samples from the small intestines of pigs that were positive for PEDV but negative for TGE virus (TGEV), porcine rotavirus (PRV), classical swine fever virus (CSFV), and porcine circovirus 2 (PCV2) were homogenized by adding phosphate-buffered saline (PBS) to achieve a 10 % (w/v) homogenate. These were then centrifuged (10,0009g, 15 min, 4°C) in a Sorvall ST 16R microcentrifuge (Fisher Scientific, MA, USA) and filtered through a 0.22-lm filter (Millipore, MA, USA) for use as inocula. Vero E6 cells were seeded onto 21-cm 2 tissue culture plates 24 h prior to infection. The growth medium was Dulbecco's modified Eagle's medium (DMEM) (Life Technologies, NY, USA) supplemented with 10 % heat-inactivated fetal calf serum, 100 U/ml penicillin G, and 100 lg/ml streptomycin. When confluent, growth medium was removed, and cell monolayers were washed once with PBS and then with DMEM containing 0.3 % tryptose phosphate broth and 0.02 % yeast extract. The washed cells were infected with the filtered intestinal preparation that had been resuspended in 'infection medium' (DMEM, 0.3 % tryptose phosphate broth, 0.02 % yeast extract, and 10 lg/ml trypsin). After adsorption at 37°C for 1 h, an equal volume of infection medium was added. Every 24 h, the culture supernatants were supplied with trypsin at 80 % of the initial concentration. PEDV strain CH/ZMDZY/11 was identified by amplifying structural and non-structural genes of the virus using PEDV-specific primers in a reverse transcription polymerase chain reaction (RT-PCR). Other Chinese PEDV isolates and reference strains Another 11 Chinese PEDV isolates from 3-to 10-day-old suckling piglets suffering from watery diarrhea, dehydration and vomiting were examined. For these isolates, from different farms across seven provinces within China, morbidity of affected piglets was 70-100 %, and mortality, 60-90 %. Some of the breeding herds from these farms had been immunized with CV777-based vaccines. Information pertaining to the reference strains and 11 PEDV isolates we analyzed is summarized in Table 1 . Nucleotide sequences were edited, aligned and analyzed with the MegAlign program, a part of the Lasergene sequence analysis software package. Phylogenetic trees were constructed using Molecular Evolutionary Genetics Analysis (MEGA) software version 5.02, applying the neighbor-joining method, which was based on the nucleotide sequence of each gene. Cytopathic effects, such as cell fusion and syncytia formation, typical of coronavirus infection, were observed after 11 passages in Vero-E6 cells. Genes encoding the spike protein (S), envelope protein (E), membrane protein (M), and nucleoprotein (N), along with open reading frame 3 (ORF3), were successfully amplified by RT-PCR ( Fig. 1 ). Our sequencing results showed that the amplified fragments shared high homology to other reported PEDV strains. The genome of the CH/ZMDZY/11 strain comprises 28,038 nucleotides (nt), excluding the poly (A) tail. The organization of genes is similar to that for other reported PEDVs (5 0 -ORF1a/1b-S-ORF3-E-M-N-3 0 ), with seven coding regions identified [1, 20] . The 5 0 and 3 0 ends of the genome contain untranslated regions (UTRs) of 292 and 334 nt, respectively. ORF1a/1b encodes two large nonstructural precursor polyproteins, replicases 1a and 1b. The S gene encodes the S protein (150-220 kDa), E encodes the E protein (7 kDa), M encodes the M protein (20-30 kDa), and the N gene encodes N (58 kDa). ORF3, between the S and E genes, encodes a non-structural protein. The nucleotide and deduced amino acid (aa) sequences of certain reference strains and 11 Chinese PEDV isolates were examined and compared with that for the CH/ ZMDZY/11 strain (Table 1) . Of these isolates, 10 were identified after 2010 from PEDV-affected swine farms of major swine-breeding provinces in China. The remaining two isolates, CH-S and LZC, were isolated in 1986 and 2006, respectively [4] . The length of ORF1a/1b for the 10 Phylogenetic analysis of the nucleotide and amino acid sequences of ORF1a/1b confirmed that all field strains belonged to three groups (Fig. 3a) . Group 1 comprised CV777 (Europe), LZC and SM98 (South Korea). Group 2 only contained attenuated DR13 (South Korea). The 10 Chinese PEDV strains isolated after 2010 were grouped together with virulent DR13 and CH-S (China) in group 3. Pairwise alignment and comparison of the S gene from the 10 post-2010 Chinese PEDV isolates showed high levels of nucleotide sequence identity (97.4-100 %). When the S gene sequence for these isolates was compared with the European prototype CV777, we found 92.9-94.6 % nucleotide sequence identity; when they were compared against SM98, attenuated and virulent DR13 strains, nucleotide sequence identities of 92.9-94.6 % were seen. The homology of the S gene of CH-S to CV777, SM98, attenuated and virulent DR13 was 95.2-97 %, while that of LZC to CV777, SM98, attenuated and virulent DR13 was 95.2-99.1 %. Amino acid sequence analysis of the S glycoprotein gene revealed that the 10 post-2010 PEDV isolates displayed sequence homology ranging from 97 to 100 % when compared with each other. These isolates shared relative low levels of identity with earlier Chinese isolates and those from Europe and South Korea, at 92. 6 Phylogenetic analysis of the S gene and protein confirmed that all field strains could be classified into three and D?N at 168. There was a continuous 17-aa deletion in ORF3 of attenuated DR13 (82YCPLLYYCGALLD-ATII98); however, such a deletion did not exist in CH-ZMDZY-11 and other Chinese PEDV strains. Phylogenetic analysis based on ORF3 nucleotide and amino acid sequences showed that PEDV Chinese field strains clustered into three groups (Fig. 3c) Table 1 Porcine epidemic diarrhea virus strains in China 2491 identity and 94.3-100 % amino acid similarity) when the E sequences from post-2010 isolates were compared with earlier Chinese strains and reference strains. The following nucleotide changes in post-2010 PEDV strains and virulent DR13 were observed at the following nucleotide positions: C?T at 90, 165, and 166, and G?A at 194. The change at 194 resulted in an altered amino acid at residue 65 (R?Q). Attenuated DR13 has a 21-nt deletion from 67-87, which corresponds to a 7-aa deletion. Phylogenetic analysis showed that Chinese field isolates belonged to three groups (Fig. 3d) . Group Phylogenetic analysis of the N gene confirmed that all field strains could be classified into three groups (Fig. 3f) . Group 1 contained CV777 (Europe), LZC (China) and SM98 (South Korea). Group 2 comprised CH-S (China) and attenuated DR13 (South Korea), and group3 contained the 10 post-2010 Chinese PEDV strains and virulent DR13 (South Korea). In the present study, a virulent PEDV strain was isolated from a suckling piglet, and its genome was sequenced and compared with those of other PEDV strains. To the best of our knowledge, we are the first to report the genetic characteristics of viruses causing PED outbreaks in China. Comparative analysis of the whole genome sequence indicated that CH/ZMDZY11 shared high nucleotide sequence identity and amino acid sequence similarity with PEDV strains isolated between 2010 and 2012. The PEDV strain prevailing in China has developed several mutations in its genome. The majority of the 10 post-2010 Chinese PEDV strains were isolated from farms where breeding herds were immunized with inactivated or attenuated PEDV vaccines. These were manufactured with the CV777 vaccine strain; however, the morbidity and mortality of suckling piglets suffering from PED remains high, indicating low efficacy of the current vaccine in use. The polymerase gene comprises two large ORFs, 1a and 1b, which encode two replicases that determine the replication efficacy of the virus. Sequence analysis of these ORFs showed that the post-2010 Chinese PEDV isolates and virulent DR13 shared multiple mutations at the same location and were clustered in one group. The 10 post-2010 PEDV strains appear to be related with the earlier prevailing Chinese strain, CH-S. The S protein is an antigen on the surface of the virion that plays a pivotal role in regulating interactions with specific host-cell receptor glycoproteins to mediate viral entry. This stimulates induction of neutralizing antibodies in the natural host [21, 22] . It is also associated with growth adaptation in vitro and attenuation of virulence in vivo [23, 24] . The S genes of the 10 post-2010 Chinese PEDV isolates were highly homologous. The S genes in all 10 isolates contained some altered amino acids, with a discontinuous 5-aa insertion and a 2-aa deletion at the same sites of the N-terminal region. The S gene has been considered a target gene for analyzing genetic variation of PEDV strains in the field [4, 6] . Our findings indicate that the current PEDV isolate prevailing in China is a novel strain with a variant S gene that is distinct from previous Chinese and reference strains. This could explain the lesseffective immunological protection against PEDV strains induced by the current vaccine in China, and the effectiveness of feeding sows feces or intestinal contents of infected piglets. Further investigation of the antigenic differences among current isolates and the vaccine strain used in China is necessary for the development of an effective vaccine against PED in China. ORF3 is an accessory gene, located between the S and E genes, and its function has rarely been examined. Sequence analysis revealed a continuous 17-aa deletion in ORF3 of attenuated DR13 and the CV777 vaccine strain. However, this deletion does not exist in current Chinese PEDV strains isolated after 2010. Differences in ORF3 sequences between cell-adapted and field viruses could be a marker of their adaptation to environmental conditions. A recently study suggested that PEDV ORF3 is related to viral infectivity and pathogenicity and encodes an ion channel protein that regulates virus production [25] . The post-2010 isolates were clustered together with virulent DR13, implying that current Chinese PEDV strains retain high infectivity and pathogenicity. The effects of this deletion on virulence and other properties requires further investigation. The E protein plays an important role during coronavirus budding. It transiently resides in a pre-Golgi compartment before progressing to the Golgi apparatus [26] . Recently, many studies have focused on the genetic and phylogenetic analysis of PEDV S, M and ORF3 genes, whereas a comparison of E genes has rarely been conducted. Our results show a 21-nt deletion in the E gene of attenuated DR13, which does not exist in post-2010 Chinese PEDV isolates and virulent DR13. All post-2010 Chinese PEDV strains exhibited a close relationship to the South Korean virulent DR13 and its derivative. The M protein plays an important role in the viral assembly process and also induces antibodies that neutralize the virus [27, 28] . Heterogeneity of the M gene sequence is known to confer diverse PEDV pathogenicities [5] . Sequence comparisons indicate that the 10 post-2010 Chinese PEDV strains contain an amino acid change (V46A) compared with earlier Chinese isolates and other reference strains. The post-2010 isolates were clustered together with virulent DR13. The N protein, which binds to virion RNA and provides a structural basis for the helical nucleocapsid, is a basic phosphoprotein associated with the genome [29] . Alignment with reference strains indicated that the N gene is highly conserved, with the exception of some point mutations. Phylogenetic analysis confirmed that all 10 post-2010 Chinese PEDV isolates were grouped with virulent DR13. Taken together, our results show that a novel PEDV strain is the etiological agent contributing to the re-emergence of PED in China. Prevailing PEDV field isolates share common molecular characteristics in their S, E, and M genes. Phylogenetic analyses of ORF1a/b, ORF3, E, M, and N genes confirmed that the 10 post-2010 Chinese PEDV isolates have a close relationship to a virulent DR13 strain from South Korea. The differences between the current isolates and the vaccine strain must be defined to develop an effective vaccine that better prevents and controls PED epidemics. Porcine epidemic diarrhea A new coronavirus-like particle associated with diarrhea in swine Virus-like particles associated with porcine epidemic diarrhoea Complete genome sequence of a Chinese virulent porcine epidemic diarrhea virus strain Chinese-like strain of porcine epidemic diarrhea virus Heterogeneity in spike protein genes of porcine epidemic diarrhea viruses isolated in Korea Sequence analysis of the partial spike glycoprotein gene of porcine epidemic diarrhea viruses isolated in Korea Complete genome sequence of porcine epidemic diarrhea virus strain AJ1102 isolated from a suckling piglet with acute diarrhea in China Complete genome sequence of a variant porcine epidemic diarrhea virus strain isolated in China Complete genome sequence of a porcine epidemic diarrhea virus variant Complete genome sequence of a novel porcine epidemic diarrhea virus in south China Complete genome sequence of a highly prevalent isolate of porcine epidemic diarrhea virus in south China Complete genome sequence of novel porcine epidemic diarrhea virus strain GD-1 in China Complete genome sequence of a vero cell-adapted isolate of porcine epidemic diarrhea virus in eastern China Complete genome sequence of a virulent porcine epidemic diarrhea virus strain Porcine epidemic diarrhea virus strains in China 2493 Complete genome sequences of a Korean virulent porcine epidemic diarrhea virus and its attenuated counterpart Phylogenetic analysis of porcine epidemic diarrhea virus field strains prevailing recently in China Complete genome sequence of a variant porcine epidemic diarrhea virus strain isolated in central China Isolation and characterization of a variant porcine epidemic diarrhea virus in China Porcine epidemic diarrhoea virus: a comprehensive review of molecular epidemiology, diagnosis, and vaccines 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 Cloning and further sequence analysis of the spike gene of attenuated porcine epidemic diarrhea virus DR13 Mutations in the spike gene of porcine epidemic diarrhea virus associated with growth adaptation in vitro and attenuation of virulence in vivo PEDV ORF3 encodes an ion channel protein and regulates virus production Porcine epidemic diarrhea virus E protein causes endoplasmic reticulum stress and up-regulates interleukin-8 expression Protein interactions during coronavirus assembly Organization of two transmissible gastroenteritis coronavirus membrane protein topologies within the virion and core Heterologous gene expression from transmissible gastroenteritis virus replicon particles Acknowledgments This work was supported by the Program for Science & Technology Innovation Talents in the Universities of Henan Province (Program no. 2011HASTIT009).