key: cord-0992788-2gp1mqfw authors: Zhao, Xiaoya; Li, Zhili; Zeng, Xiduo; Zhang, Guanqun; Niu, Jianqiang; Sun, Baoli; Ma, Jingyun title: Sequence analysis of the spike gene of Porcine epidemic diarrhea virus isolated from South China during 2011–2015 date: 2017-06-22 journal: J Vet Sci DOI: 10.4142/jvs.2017.18.2.237 sha: 7818d1ffe60f6ea1d6b44e73402132ce41453ab2 doc_id: 992788 cord_uid: 2gp1mqfw The spike gene of porcine epidemic diarrhea virus (PEDV) was sequenced from 55 South China field strains isolated from pigs with symptoms of diarrhea. The sequences were compared within the set of field strains as well as with reference strains available in GenBank. Within the 55 South China PEDV field strains, the deduced amino acid sequence identities ranged from 93.8% to 99.9 % and ranged from 90.7% to 99.5% when compared with the foreign reference strains in GenBank. Our phylogenetic analysis showed that 10 of the 55 South China PEDV strains belonged to G1b and 45 belonged to G2b. Porcine epidemic diarrhea virus (PEDV) is a major cause of aqueous diarrhea, vomiting, and severe dehydration in pigs. It is transmitted via direct or indirect fecal-oral routes [10, 14] . As other Coronaviridae members, PEDV is an enveloped, positive-sense, single-stranded RNA virus with a small 28 kb genome. The genomic RNA has a 5′-cap and a 3′-polyadenylated tail and contains at least seven open reading frames that encode four structural proteins: spike (S), envelope (E), membrane (M), and nucleocapsid (N) [9] . Like other coronaviruses, the S protein is a type 1 transmembrane envelope glycoprotein comprising 1,383 amino acids and contains 29 predicted glycosylation sites. The S1 and S2 domains of this protein have essential roles in cellular receptor interactions that mediate viral entry and induce neutralizing antibodies in the natural host [13] . The PEDV S protein consists of three domains: a large extracellular domain, a transmembrane region, and a short cytoplasmic carboxyl terminus. There is a region of the S protein that contains the epitope(s), which is capable of inducing PEDV-neutralizing antibodies [1] . Research indicates that phage-displayed peptides, which were panned against 2C10, have antigenic similarities with the PEDV antigen epitope motif ( 1368 GPRLQPY 1374 ) and can show neutralizing activity against PEDV [6] , the CO-26K equivalent (COE) domain (aa positions 99-638), and the epitopes SS2 (aa positions 748-755) and SS6 (aa positions 764-771), which can also induce neutralizing antibodies against PEDV [4, 25] . Recently, the S gene was implicated as an important determinant of the biological properties of PEDV. These properties include genetic relationships between PEDV isolates, the epidemiological status of PEDV in the field, and associations between genetic mutations and virus function [16] . Porcine epidemic diarrhea (PED) was first reported in the United Kingdom in 1971, although the virus was not identified until 1978 when it was designated PEDV with the prototype strain CV777 [11] . Subsequently, PEDV infections have been reported in China, Japan, Korea, and Thailand and for the first time in the United States in 2013 [4, 17, 22, 24] . In China, the first occurrence of PED was in 1973, although the causative agent was not identified until 1984 [27] . Since the end of 2010, PEDV infections have been detected in more than 10 provinces in southern China. According to the previous study, there was a high prevalence rate of PED in pig herds [3] . Clinical signs of PED are seen in piglets within 7 days of birth and present as vomiting and aqueous diarrhea. Death rates during this stage range from 80% to 100% resulting in huge economic losses in the pig industry [17] . During the period from 2010 to 2013, the incidence of diarrhea among newborn piglets from 60 large-scale farms in Guangdong and some other southern provinces was 59,320, of which 45,601 piglets died, for an average mortality rate of 76.87%. These studies indicate that PED is an increasingly serious problem in China. In this study, we sequenced and analyzed the S gene of 55 PEDV strains from South China. The aim was to provide an improved understanding of the origin, evolution, and diversity of the strains in this region, and to determine the phylogenetic relationships between newly isolated and historic strains. Between February 2011 and March 2015 in Guangdong province, China, porcine intestinal and fecal samples were collected from piglets with watery diarrhea and dehydration. The intestinal and fecal samples were placed in phosphate buffer saline supplemented with penicillin G (10,000 IU/mL) and streptomycin (2 mg/mL) and centrifuged at 1,500 × g for 15 min. Viral presence in the supernatant was determined by using RT-PCR with viral gene-specific primers: (5′ to 3′) GCAACTCAAGTGTTCTCAG and GAGTCATAAAAGAAACGTCCG. RNA extraction and reverse transcription were carried out as previously described [27] . PCR products were purified by using an AxyPrep DNA gel extraction kit (Axygen, USA) according to the manufacturer's instructions, and the product was cloned into the pMD-19T vector (Takara Bio, China). The purified recombinant plasmids were sequenced by Invitrogen Trading (Invitrogen, China). Sequences of all samples detected to be positive for the complete S gene of PEDV were submitted to GenBank (National Center for Biotechnology Information, USA) under accession Nos. KP399601-KP399634 and accession Nos. KR296663-KR296683. Nucleotide and deduced amino acid sequences were aligned, edited, and analyzed with Clustal X (ver. 1.83) [5] , BioEdit (ver. 7.0.5.2), and DNASTAR software. The phylogenetic tree was constructed by using MEGA (ver. 5.05) [26] incorporating the neighbor-joining (NJ) method. Bootstrap values were estimated for 1,000 replicates. Sequence analysis revealed that the S genes of the 55 Chinese PEDV strains ranged from 4149 to 4176 nucleotides in length. Compared to the S gene of CV777, sequencing of the field strains revealed a range of lengths: one was three nt shorter (4149 nt), nine were the same length (4152 nt), four were six nt longer (4158 nt), 39 were nine nt longer (4161 nt) and two were 24 nt longer (4,176 nt). The deduced amino acid sequences of the 55 S proteins ranged from 1382 to 1391 aa in length. While the SS2 epitope was conserved in all 55 Chinese PEDV strains, most strains exhibited two or three amino acid mutations in the SS6 epitope ( Fig. 1) . Among the 55 South China field isolates, 53 strains (i.e., excluding CH-YGC-11-2013 and CH-XLC-03-2013) had the same mutation (Y/F766S) in SS6. We also found that 52 strains (i.e., excluding CH-YGC-11-2013, CH-TPC-03-2013, and CH-CCC-05-2013) had the same mutations (G/S1368S, L/F1371S, and L/X1371S) in the 2C10 region ( Fig. 1 ). Phylogenetic analysis of S gene showed that all PEDV strain sequences could be divided into four groups ( (MB-021). Forty-five of the fifty-five South China PEDV strains belonged to G2b. Sequence alignment showed nucleotide (deduced amino acid) identities in the range 94.4% to 99.9% (93.8%-99.9%) among the 55 Chinese PEDV strains and identities in the range 92.4% to 99.76% (90.7%-99.5%) with the foreign reference strains reported in GenBank. Some vaccine strains including CV777, attenuated DR13, SM98, and 83P-5_100th-passaged are used in a variety of countries. The S genes of our 55 recently isolated Chinese field strains shared identity with CV777 (93.8%-95.7%), CHGD-01 (94.6.0%-98.3%), SM98 (93.3%-95.4%), attenuated DR13 (93.3%-95.6%), and 83P-5_100thpassaged (93.7%-96.0%). The deduced amino acid sequences of these strains shared identity with CV777 (93.0%-96.1%), CHGD-01 (94.1%-98.3%), SM98 (92.1%-95.3%), attenuated DR13 (92.0%-95.7%), and 83P-5_100th-passaged (92.4%-96.1%). The S genes of the 55 recently isolated Chinese field - - - - N N - - - CH-STC-12-2014 N - - - N N - - - CH-TPC-12-2014 N - - - N N - - - CH-CLC-01-2015 N - - - N N - - - CH-GGC-01-2015 N - - N N N - - - CH-GLC-02-2015 N - - - N N - - - CH-HFEC-01-2015 - - - - N N - - - CH-HGC-01-2015 N - - - N N - - - CH-HSYKFC-01-2015 N - - - N N - - - CH-JPYC-02-2015 N - - - N N - - - CH-QTC-01-2015 N - - N N N - - - CH-SBC-03-2015 N - - - N N - - - CH-WTC-02-2015 N - - - N N - - - CH-XBC-01-2015 - - - - N N - - - CH-YGC-01-2015 - - - - N N - - - CH-ZSC-01-2015 N - N N N N - - - CH-ZWBZa-01-2015 - - - - N N - - - CH-ZWBZb-01-2015 N - - - N N - - - CH-ZWC-01-2015 - - - N - N - - - South All of the analyzed strains contained between nine and ten high-specificity N-glycosylation sites, three of which were conserved. The CV777 strain is the only vaccine source in China. Therefore, a comparative analysis of all 55 strains and CV777 was conducted. All of the Chinese field strains, except CH-SHC (Table 1 ). This study was carried out to elucidate the diversity of PEDV in South China. We determined the phylogenetic relationships between new and historic PEDV strains in this region and compared them with globally reported strains. To that end, we sequenced and analyzed the S glycoprotein gene of 55 Chinese field strains isolated from South China during 2011-2015. All 55 Chinese PEDV strains had a high sequence similarity with the previously identified Guangdong strains (CHGD-01, GDQY/2011, GDZQ/2012, GXNN/2013, and GDS01). The results further indicate the homogeneity in the S genes of PEDV field strains prevalent in South China. Our results also demonstrate that the recently identified Guangdong PEDV strains are closely related to those isolated from the USA, indicating that the USA strains may have originated from Chinese variants [12] . High morbidity and mortality were observed in pigs that were immunized with inactivated or attenuated PED vaccines based on CV777 [19] , suggesting that the CV777 vaccine is unable to provide full protection in swine. We have reported on further investigations into the S gene mutations responsible for the failure or otherwise of the existing CV777 vaccine [28] . The S genes of attenuated strains (CV777-attenuated, attenuated DR13, and 83P-5_100th-passaged) are 4149 nt in length, while the classical virulent strains, such as CV777, are 4152 nt long [16, 18] . Sequence analysis revealed length diversity among the 55 Chinese field strains in this study. The S gene of one field strain was shorter than those of the classical strains due to nucleotide deletions, while the S genes of 45 field strains were longer due to nucleotide insertion. These variations are indicative of the unique characteristics of the S genes of the 55 field strains in the present study. With the exception of the CH-YGC-11-2013 and CH-XLC-03-2013 strains, the remaining 53 South China strains contained several amino acid substitutions compared with the vaccine strains, resulting in changes in the amino acid constitution of the SS6, 2C10, and COE neutralizing epitopes. A retrospective study of the three amino acid substitutions revealed that the shift process was present not only in European strains but also in Asian PEDV strains [8] . Those results in combination with our results indicate that the three substitutions represent a marker of increased PEDV viability [8] . The domestic outbreak of PED in China was caused by new PEDV strains, the S genes of which had similar molecular characteristics [7] . According to previous studies as well as the present study, the epidemic strains of PEDV have a close genetic relationship with Asian strains including those isolated in China, Korea, and Japan after 2010 [2, 15, 20] . Therefore, a detailed description of the differences between epidemic strains and vaccine strains could aid in the development of improved vaccines against PEDV. Previous research has shown that the Severe Acute Respiratory Syndromes (SARS) corona virus (CoV) originated from Guangdong Province in southern China and has a recombinant history with lineages of type I coronavirus such as PEDV [21] . Moreover, the high frequency of recombination events in CoVs may result in the generation of novel viruses with high genetic diversity, and it could invoke unpredictable changes in virulence during human infections [23] . As such, the information we have gathered about the mutation of S genes could be informative for monitoring and predicting the potential emergence of a recombinant CoV from animals. In summary, our sequence analysis of the S genes and molecular epidemiology of 55 field strains from South China indicate that the PEDVs circulating between 2011 and 2015 in this region have evolved from historic strains within the local region. The results of our analyses will be useful when improving the PED control strategies in this region. Identification of the epitope region capable of inducing neutralizing antibodies against the porcine epidemic diarrhea virus Complete genome sequence of a variant porcine epidemic diarrhea virus strain isolated in China Detection and molecular diversity of spike gene of porcine epidemic diarrhea virus in China Molecular epidemiology of porcine epidemic diarrhea virus in China Multiple sequence alignment with the Clustal series of programs Phage-displayed peptides having antigenic similarities with porcine epidemic diarrhea virus (PEDV) neutralizing epitopes Complete genome sequence of a novel porcine epidemic diarrhea virus in South China Bioinformatics insight into the spike glycoprotein gene of field porcine epidemic diarrhea strains during 2011-2013 in Guangdong, China Porcine epidemic diarrhea virus infection: etiology, epidemiology, pathogenesis and immunoprophylaxis Genetic characterization of porcine epidemic diarrhea virus in Korea from Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea virus (PEDV) field strains in south China Pathogenesis of porcine epidemic diarrhea virus isolate (US/Iowa/18984/2013) in 3-week-old weaned pigs Immunogenicity and protective efficacy of recombinant S1 domain of the porcine epidemic diarrhea virus spike protein Novel porcine epidemic diarrhea virus variant with large genomic deletion Complete genome sequences of a Korean virulent porcine epidemic diarrhea virus and its attenuated counterpart Sequence analysis of the partial spike glycoprotein gene of porcine epidemic diarrhea viruses isolated in Korea Chinese-like strain of porcine epidemic diarrhea virus Mutations in the spike gene of porcine epidemic diarrhea virus associated with growth adaptation in vitro and attenuation of virulence in vivo Molecular characterizations of subcellular localization signals in the nucleocapsid protein of porcine epidemic diarrhea virus Porcine epidemic diarrhoea virus: a comprehensive review of molecular epidemiology, diagnosis, and vaccines Evidence from the evolutionary analysis of nucleotide sequences for a recombinant history of SARS-CoV Emergence of porcine epidemic diarrhea virus in the United States: clinical signs, lesions, and viral genomic sequences Epidemiology, genetic recombination, and pathogenesis of coronaviruses An immunohistochemical investigation of porcine epidemic diarrhea Identification of two novel B cell epitopes on porcine epidemic diarrhea virus spike protein MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods Molecular characterization and phylogenetic analysis of new variants of the porcine epidemic diarrhea virus in Gansu Genetic variation analysis of reemerging porcine epidemic diarrhea virus prevailing in central China from This research was supported by a grant from the Science and Technology Planning Project of Guangdong Province China (grant No. 2011B090400414). The authors declare no conflicts of interest.