key: cord-0872605-rqlb8cjg authors: Fan, Jing-Hui; Zuo, Yu-Zhu; Li, Jian-Hui; Pei, Li-Hua title: Heterogeneity in membrane protein genes of porcine epidemic diarrhea viruses isolated in China date: 2012-05-15 journal: Virus Genes DOI: 10.1007/s11262-012-0755-4 sha: 41cbc4bd8b333f0225d4050b8f012a650681314a doc_id: 872605 cord_uid: rqlb8cjg Since late 2010, porcine epidemic diarrhea virus (PEDV) has been re-emerging in immunized swine herds with devastating impact in the Hebei province of China. Seven prevailing strains of PEDV were isolated from fecal samples out of piglets suffering from severe diarrhea. The M gene of the seven PEDV isolates encompasses an open reading frame of 681 nucleotides, encoding a protein of 226 amino acids. The seven PEDV isolates showed 99.4–99.9 % nucleotide sequence identity and 98.2–99.1 % deduced amino acid identity. When compared with other Chinese isolates and foreign isolates, the seven isolates showed high nucleotide identity with the Thailand isolate M-NIAH1005 (99.6–99.9 %) and Korea isolate PFF188 (99.7–100 %), but low identity with other Chinese isolates (96.6–99.1 %) and with the vaccine strain CV777 used in China (97.8–98.2 %). Phylogenetic analyses showed that all seven Chinese field isolates were grouped together in the same cluster. Although CV777 was also separated into the same cluster with the seven isolates, they were belonged to different sub-cluster. These results showed that the seven prevailing isolates in China are closely related phylogenetically to each other and have close relationships with the Korean strain PFF188 and Thailand strain M_NIAH1005. However, they differ genetically from other Chinese isolates and the vaccine strain CV777. Therefore, a more efficient vaccine strain should be chosen to prevent outbreaks of PEDV in China. Porcine epidemic diarrhea virus (PEDV), a member of the family Coronaviridae, is an enveloped single-stranded RNA virus [7, 8, 17] . Its genome contains genes for the spike (S), membrane (M), small membrane (sM), open reading frame (ORF) 3, and nucleocapsid (N) proteins [2, 4] . The M protein is a structural membrane glycoprotein and the most abundant of all the envelope proteins. It has a short amino-terminal domain on the outside of the virion and a long carboxy-terminal domain on the inside. The M protein plays a central role in PEDV assembly. Together with the sM protein, it is responsible for assembly of the PEDV envelope [1, 9, 21] . In addition, the M protein directs the incorporation of the S protein [13, 15] and the nucleocapsid [12] into the budding particle. PEDV causes severe entero-pathogenic diarrhea in piglets, especially in neonates, and the disease has a high mortality rate which can reach 80 % in certain situations [19] . PEDV was first reported in Belgium and the United Kingdom in 1978. Since then, outbreaks of the disease have been reported in several swine-raising countries, most notably in Europe and Asia, including Japan, China, and Korea [11] . It causes heavy economic losses in the European and Asian swine industry and is an important gastrointestinal disease [7, 16] . Several vaccines have been applied to the regions with a high prevalence of PEDV, but these vaccines have not been completely effective in preventing the disease [6, 18] . Since early 2006, PEDV has been re-emerging in immunized swine herds in China [5, 20] . In late 2010, outbreaks of PED occurred in HeBei province. Pigs of all ages were affected and losses from this disease outbreak were extensive. Most of the affected farms lost 100 % of their newborn piglets. In order to identify the genetic heterogeneity of the prevailing strains in China, PEDV strains were isolated from fecal samples from piglets infected with PEDV in 2010. The M genes were cloned and sequenced for analysis of genetic diversity, and the phylogenetic relationships of the field isolates with other PEDV reference strains were also analyzed. Eighty-two fecal samples were collected from piglets of 10 farms where the piglets presented clinical signs related to the PED definition and were from affected herds with severe diarrhea. Seven PEDVs were isolated from the feces that presented positive for PEDV. These isolates were used for nucleotide sequence analysis, amino acid sequence analysis, and phylogenetic analysis in this study, and are described in Table 1 . Restriction endonucleases, polymerases, and ligases were purchased from TaKaRa Biotechology (Dalian, China) Co., Ltd. Extraction of genomic PEDV RNA Each fecal sample was centrifuged using a Beckman F3602 rotor at 2,000 g at 4°C for 5 min. The supernatant was used for the extraction of the viral genomic RNA using Trizol reagent (Invitrogen, USA) based on the manufacturer's instructions. The extracted RNA was resuspended in RNasefree water and was stored at -80°C. Primers were designed on the basis of the published sequence of the M gene of PEDV. The forward and reverse primers designed to amplify the whole sequence of M gene of PEDV isolates were MP1 5 0 -GGATCCATGTCTAACGGT-3 0 and MP2 5 0 -AAGCTTTCTGTTTAGACTAAAT-3 0 . The primers were synthesized commercially. The conditions for amplification included a reverse transcription process. Reverse transcription was carried out using the primer Oligo (dT) 18. The mixture of 1 ll virus RNA, 1 ll of 50 pM Oligo (dT)18, and 5 ll of RNase-free dH 2 O was incubated at 70°C for 5 min and placed on ice for at least 2 min. Subsequently, 2 ll of 59 RT buffer, 0.5 ll of 10 mM dNTP, 0.25 ll of RNase inhibitor (40 U/ll), and 0.25 ll of M-MLV reverse transcriptase (200 U/ll) were added to initiate the reaction at 42°C for 60 min. The reaction was stopped by heating at 70°C for 15 min. The cDNA was either stored at -20°C or amplified immediately. The M genes of the seven field isolates were amplified by PCR using primers MP1 and MP2. The PCR conditions Cloning and sequencing The M gene fragments obtained using RT-PCR were separated by gel electrophoresis and purified using a Gel Extraction Mini Kit (TIANGEN BIOTECH (Beijing) Co., Ltd.) according to the manufacturer's recommendation. The purified product was cloned into a pMD18-T plasmid DNA vector, and the resulting plasmid was named pMD-M. The recombinant plasmid was used to transform competent Escherichia coli JM109. The positive clones were selected by their lacZ-negative phenotypes and verified by restriction enzyme digestion, PCR, and DNA sequencing. Nucleotide sequence accession number The complete nucleotide sequences of the M genes of the seven PEDV isolates have been deposited in the GenBank Database. Their accession numbers are given in Table 1 . Comparison and phylogenetic analysis of the nucleotide sequences of the M genes of the seven PEDV isolates with those of other PEDV strains were performed with PHYLIP v3.69 software. The reference strains used for the sequence analysis are described in Table 1 . Phylogenetic relationships were determined by neighbor-joining method with PHYLIP v3.69 and reliability of topologies was estimated by performing bootstrap analysis with 1,000 replicates. TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/services/ TMHMM-2.0/) was used to predict the transmembrane helices in M proteins. Antigenic index was evaluated by Jameson-Wolf method. The 82 fecal samples were submitted to pathogen detection and isolation. Detection results showed that 51 out of 82 (62. The nucleotide and deduced amino acid sequences of our isolates were compared with those of the reference strains. The sequence homology among these seven isolates was 99.4-99.9 % at the nucleotide level. When compared with the Chinese reference strains, the seven isolates had 96.6-99.1 % sequence homology with them: JS/2004 had 98.7-99.1 % sequence homology with the seven isolates, while LZC had 96.6-97.1 % homology with them. All seven isolates had low nucleotide identity (97.8-98.2 %) with CV777. HB/GY and HB/LF had only 97.8 % sequence homology with CV777. When the field isolates were compared with reference strains from other countries, the nucleotides identity values were 96.2-100 %. The highest nucleotides identity was observed between the seven field isolates and PFF188 and M-NIAH1005. The seven isolates had 99.7-100 % sequence homology with PFF188 and 99.6-99.9 % sequence homology with M-NIAH1005. HB-HS had 100 % sequence homology with PFF188 and 99.9 % with M-NIAH1005, while HB-GY and M-NIAH380 only showed a low nucleotide identity of 96.2 %. The deduced amino acid sequences of the field isolates showed 98.2-99.1 % homology, and they showed homology of 94.7-100 % with the deduced amino acid sequences of the reference strains. To analyze the phylogenetic relationships of the seven PEDV isolates with the reference strains isolated in various parts of the world, we constructed a neighbor-joining phylogenetic tree using the M gene sequences. A representative tree for the M gene is shown in Fig. 1 . The seven PEDV isolates and reference strains were separated into three potential clusters. All seven field isolates were segregated into the same cluster. In terms of the topology of the phylogenetic tree, the seven PEDV isolates might have a common origin. Porcine epidemic diarrhea virus (PEDV) is a member of the coronavirus group causing watery diarrhea, dehydration, and high mortality in suckling pigs. Although the bi-combined attenuated vaccine against TGEV and PEDV infection is authorized for use on swine farms, PED still occurs in the immunized swine herds in China. Since late 2010, outbreaks of PED have occurred in Hebei province. Most of the affected farms lost 100 % of their newborn piglets. PEDV is known to have only one serotype, but the genome shows wide genetic diversity [3, 10, 22] . The structural M protein of PEDV plays an important role in the assembly process and induces virus neutralization antibodies [14] . Therefore it is necessary to investigate the prevalence of PEDV on the basis of phylogenetic analysis of the M gene of the field strains of PEDV in China. This study was conducted to characterize PEDVs isolated from the fecal samples of piglets suffering from severe diarrhea in Hebei province of the People's Republic of China. The initial aim of this study was to isolate and investigate the genetic diversity among field strains of PEDV in China. A large amount of sequence data on PEDVs isolated throughout the world has been published over the years and is now available for sequence comparison and phylogenetic analysis. In this study, seven Chinese field strains of PEDV were isolated. The complete nucleotide and deduced amino acid sequences of the M gene were determined and compared with those of the reference strains. None of the seven isolates were found to have sequence insertions or deletions in the M genes, although all the PEDV isolates investigated in this study had nucleotide mutations. Analysis of sequence homology showed that the seven PEDV isolates were closely related to PFF188, M-NIAH1005, and JS-2004. Furthermore, HB-HS and PFF188 had identical nucleotide sequences, although they were isolated in different countries at different times. This suggests that HB-HS might have been generated naturally from PFF188. However, lower amino acid identity was observed between the seven PEDV isolates and the Chinese isolates CH-IMT-06, HN/XYYYP, YM/2007. Phylogenetically, they belonged to three different clusters. Although the seven isolates and the other Table 1 . Tree topology was constructed using F84 model and bootstrap re-sampling (1,000 data sets) of the multiple alignments was used to test the statistical robustness of the trees obtained by NJ (using the program Neighbor from PHYLIP v3.69 package) d vaccine strain m epidemic strain Chinese isolates belonged to the same cluster, they were segregated into different sub-clusters. These results showed that a new genotype of PEDV prevails in China now. None of the seven PEDV isolates was phylogenetically close to the vaccine strain CV777 that is used in China. In fact, they were segregated into different sub-clusters. In addition, the seven Chinese isolates also had a change in the antigenic regions 2-14 though there has no scientific data to prove whether this change had brought about antigenic variation. All these results show that, in spite of the regular use of live vaccines in pigs, the seven field isolates of PEDV were of a different genotype able to cause outbreaks in immunized pigs in Hebei province. In conclusion, genetic variation exists in the M gene of the prevalent Chinese PEDV field isolates, and a more appropriate vaccine strain should be chosen to prevent PEDV outbreaks in China. Interferon alpha inducing property of coronavirus particles and pseudoparticles Sequence determination of the nucleocapsid protein gene of the porcine epidemic diarrhea virus confirms that this virus is a coronavirus related to human coronavirus 229E and porcine transmissible gastroenteritis virus Further analysis of the genome of porcine epidemic diarrhoea virus Identification of the epitope region capable of inducing neutralizing antibodies against the porcine epidemic diarrhea virus Molecular characterization and phylogenetic analysis of membrane protein genes of porcine epidemic diarrhea virus isolates in China Antiviral activity of quercetin 7-rhamnoside against porcine epidemic diarrhea virus Experimental infection of pigs with a new porcine enteric coronavirus, CV 777 In vivo morphogenesis of a new porcine enteric coronavirus, CV777 Assembly of spikes into coronavirus particles is mediated by the carboxy-terminal domain of the spike protein Detection of porcine epidemic diarrhea virus using polymerase chain reaction and comparison of the nucleocapsid protein genes among strains of the virus Isolation of porcine epidemic diarrhea virus (PEDV) in Korea Characterization of the coronavirus M protein and nucleocapsid interaction in infected cells Protein interactions during coronavirus assembly Research overview of porcine epidemic diarrhea virus Envelope glycoprotein interactions in coronavirus assembly Porcine epidemic diarrhea virus A new coronavirus-like particle associated with diarrhea in swine Antiviral activity of Zanthoxylum species against porcine epidemic diarrhea virus Use of an internal control in a quantitative RT-PCR assay for quantitation of porcine epidemic diarrhea virus shedding in pigs Preparation and characterization of the monoclonal antibodies against neutralizing epitope region of PEDV s protein Nucleocapsid-independent assembly of coronavirus-like particles by co-expression of viral envelope protein genes Cloning and sequence analysis of the spike gene of porcine epidemic diarrhea virus Chinju99