key: cord-0890050-bennvi2x
authors: Li, Dongliang; Feng, Hua; Liu, Yunchao; Chen, Yumei; Wei, Qiang; Wang, Juan; Liu, Dongmin; Huang, Huimin; Su, Yunfang; Wang, Dongyu; Cui, Yinglei; Zhang, Gaiping
title: Molecular evolution of porcine epidemic diarrhea virus and porcine deltacoronavirus strains in Central China
date: 2018-06-13
journal: Res Vet Sci
DOI: 10.1016/j.rvsc.2018.06.001
sha: cd876d54e28e21dd65a06f43a23a182e035d4239
doc_id: 890050
cord_uid: bennvi2x

Porcine epidemic diarrhea virus (PEDV) and porcine deltacoronavirus (PDCoV) are epizootic swine viruses. To detect and study the evolution of PEDV and PDCoV in central China (Shanxi, Henan, Hubei province), 70 clinical intestinal and fecal samples from piglets with severe watery diarrhea during August 2015 and June 2016 were collected, tested and analyzed. PEDV was more frequently detected by PCR than PDCoV. Phylogenetic analysis of S genes showed that the 10 PEDV strains from this study clustered into G2a (n = 7) and G2b (n = 3) groups. Additionally, the three G2b strains (PEDV S2△) contained the same specific 3 nt deletion in S2 as other reference strains in G2b. Interestingly, complete genome analysis indicated that CH/hubei/2016 was closer to the US INDEL strain and G2a group. CH/hubei/2016 had one recombination event in S2 gene which may have resulted from AH2012-12 (from G2b group) and CH-ZMDZY-11 (from G2a group). Furthermore, 10 purifying selection sites in S gene indicated an adaptive evolution of PEDV in central China swine herds. These results suggested that Pandemic G2a and G2b are predominant PEDV genotype circulating in central China. In addition, the deletion and recombination identified in S gene suggested PEDV strains of central exhibited an evolutionary variety. However, whether these changes affect the pathogenicity and antigenicity of wild PEDV is unknown and is worth for further investigation.

Porcine epidemic diarrhea (PED) is an acute, highly contagious swine disease caused by porcine epidemic diarrhea virus (PEDV), which leads to severe watery vomiting, diarrhea, and high mortality in newborn piglets (Pensaert and de Bouck, 1978; Song and Park, 2012) . PEDV is an enveloped single-stranded RNA Coronavirus that belongs to the family Coronaviridae, genus Alphacoronavirus. PEDV genome (28 kb in length) has seven open reading frames (ORFs) and encodes four structural proteins: spike (S), envelope (E), membrane (M) and nucleocapsid (N) as well as three non-structural proteins: replicates 1a, 1b and ORF3 (Song and Park, 2012) . Among these proteins, PEDV S glycoprotein shows an important role on host cell entry and viral virulence and could be divided structurally into the S1 and S2 regions (Song and Park, 2012) , which is also a key protein exhibiting major genetic variations among different stains (Sun et al., 2015) .

vomiting, diarrhea (Woo et al., 2012) . PDCoV has been firstly identified in Hong Kong 2012, and then successively detected in China mainland and USA Woo et al., 2012) . Like PEDV S gene, the phylogenetic relationship of PDCoV is also constructed in terms of S genes Zhai et al., 2016) .

In present study, a total of 70 samples from central China during 2015 to 2016 were simultaneously tested for PEDV and PDCoV by PCR. S genes of PEDVs and complete gene of PEDV from positive samples were sequenced, respectively. Phylogenetic, sequence, recombination and selection pressures analysis of PEDV were performed to construct evolution relationships between the strains from central China and other country, whereas one PDCoV S gene was also sequenced and analyzed. This work may provide new insights in the epidemiology of PEDV in central China.

A total of 70 clinical intestinal and fecal samples were collected from different regions of central China (Henan, Hubei, Shanxi province) between 2015 and 2016. These samples came from pigs with a variety of clinical signs including vomiting, anorexia, and severe watery diarrhea. All clinical samples were homogenized for RNA extraction and virus isolation.

Total RNA was extracted from samples using Trizol reagent (TaKaRa, Japan) according to the manufacturer's instructions. Reverse transcription was performed as previously described (Kim et al., 2001) . Briefly, using 20 ul of the reaction mixture containing 11 ul of the total RNA extraction solution, 1 ul 10 mM OligodT primer, 4 ul 5 × RTbuffer, 2 ul 10 Mm dNTP, 1 ul (40 U) RNase inhibitor, 1 ul M-MLV (100 U) Reverse Transcriptase. After 1 h kept under 42°C, the resulting cDNA was stored at −20°C. All sample's cDNA was performed by PCR method as previously described (Kim et al., 2001; Song et al., 2015) . In this study, to amplify the S gene of PEDV and PDCoV primers were designed based on PEDV CV777 strain (AF353511) and PDCoV HKU15 strain (KJ568769) ( Table 1 ). The reaction mixture consisted of 1 ul of cDNA, 10 ul ExTaq (TAKARA), 1 ul F-primer (10 pmol), 1 ul R-primer (10 poml), and RNase-free water in a total volume of 20 ul. And the amplification was carried out as follows: 95°C for 5 min, followed by 30 cycles of 95°Cfor 1 min, 53°C for 1 min and 72°Cfor 1 min, and then final elongation for 10 min at 72°C·The complete genome of PEDV was amplified as previously described . Products were visualized using 1% agarose gel under ultraviolet light.

Three pairs of primers were used to determine the PEDV S gene or PDCoV S gene (Table 1 ). The amplicons were purified and cloned into pMD-18 T, then sequenced by Sanger method. Each fragment was independently sequenced at least three times. The sequenced products were splicing to obtain the PEDV S gene. After analysis of the S gene, the positive samples were selected and sequenced for the complete genomes .

Reference PEDV and PDCoV strains are shown in Table 2 (Supplemental Material 1). PEDVs and PDCoV isolates from this study are shown in Table 3 (Supplemental Material 2)

The PEDV S genes, PEDV complete genomes, and PDCoV S genes were aligned using ClustalW in DNASTAR7.1. Phylogenetic analysis for each dataset was performed using neighbor-joining method with 1000 bootstrap replicates in MEGA7.

ClustalW program was used to align the full-length genome sequences of CH/hubei/2016 and reference sequence from GenBank (Table 2 ) (Supplemental Material 1). RDP3 recombination detection program was used to detect the potential recombination events. RDP, BootScan, MaxChi, Chimaera, SiScan, and 3Seq methods embedded in the RDP3 software package (p < 0.01) were utilized , only those recombination events supported by > 4 programs were considered to avoid dependence on a single methodology . Confirmation of the potential recombination events was carried out by BootScan analysis and fast neighbor-joining tree (Sun et al., 2016) .

PEDV S gene deduced amino acid sequences were aligned by ClustalW program in DNASTAR7.1 software, and selection pressures for this gene at individual codon sites were estimated by the ratio of nonsynonymous (dN) and synonymous (dS) mutations, which was calculated by the Single Likelihood Ancestor Counting method and Fixed Effects Likelihood method available at the Datamonkey online version of the Hy-Phy package, using a significance level of 0.1 (Mu et al., 2013) .

All the samples were simultaneously tested for PEDV and PDCoV by PCR. 84.2% (59/70samples) were positive for PEDV, but only 2.9% (2/ 70samples) samples were PDCoV positive, which were also PEDV positive. (These data suggested PEDV infection could be more commonly detected than PDCoV in central China.) These data suggested that PEDV could be more commonly detected than PDCoV in central China.

To reveal the characteristics of PEDV strains currently circulating in central China, 10 PEDV positive samples (Table 3 Supplemental Material 2) were selected and sequenced in terms of S gene. Compared to the S gene of CV777 (4152), sequencing of the field strains revealed a range of lengths: three were 6 nt longer (4158 nt), seven were 9 nt longer (4161 nt). The analysis of PEDV S gene showed that at N-terminal domain (NTD) of S gene, these 10 strains included two notable insertions at 166 nt to 177 nt and 337 nt to 339 nt, respectively, and one deletion between 477 nt to 483 nt at S gene compared with CV777. And same deletions and insertions segments were also observed in the same position on global Pandemic strains as previous study indicated (Huang et al., 2013) . Besides, a specific deletion from 3596 nt to 3598 nt at S2 gene was firstly found in three PEDV strains (CH/hubei/2015, CH/ hubei/2016, CH/huaiyang/2015) (Fig. 1) , respectively, which was named as PEDV S2△. From nucleotide identity analysis, the nucleotide of these 3 PEDV S2△strains from our study had 99.3%-99.8% identity to one another, while the identity of the remaining 7 strains was 98.3%-99.8% to one another. However, the identity between these 3 PEDV S2△strains and the remaining 7 strains were at 97.4%-98.5%. The 10 PEDV strains shared 96.9%-99.2% nucleotide identities with reference Pandemic strains, 95.2%-96.8% with reference INDEL strains and 93.5% -94.8% with reference Prototype strains, respectively. Moreover, the 3 PEDV S2△strains were closely related to reference Pandemic G2b strains with 97.6%-98.6% nucleotide identity. However, the remaining 7 strains from our study were closely related to reference Pandemic G2a strains with 97.6%-99.2% nucleotide identified.

The Fig. 2 phylogenetic tree was established based upon multiple sequences alignment of S gene according to (Huang et al., 2013) . All PEDV strains fell clearly into two genetic groups, designated as G2 (Pandemic strains) and G1 (Prototype PEDV).

Obviously, PEDV Pandemic strains (G2) consisted of two groups. In current study, these 3 PEDV S2△strains and other reference PEDV strains with 3 nt-deletion in S2 gene were grouped into G2b. The remaining 7 PEDV strains in this study were grouped into G2a. None of the strains isolated in this study were included in G1 (Prototype strain) and PEDV INDEL strain.

After culturing on vero cell, only CH/hubei/2016 was successfully isolated, and its genome was future amplified and sequenced . CH/hubei/2016 of complete genome had 28063 nt. Like PEDV Pandemic strains, CH/hubei/2016 strain contained a notable U insertion at 48 nt, and two notable nucleotide deletion between 72 nt and 73 nt (an A deletion) and between 84 nt and 85 nt (a 4-nt UUCC deletion) at the 5'untranslated region (UTR) in comparison with the prototype CV777 strain (Huang et al., 2013) .

From nucleotide identity analysis, CH/hubei/2016 strain had an identity at 98.2% with Pandemic strains G2b, 97.9%-98.0% with Pandemic strains G2a, 97.6%-97.7% with US INDEL strains, and 96.4%-97.2% with Prototype strains respectively.

As shown in Fig. 3 , all PEDV strains were clearly grouped into G1 and G2 group. PEDV Pandemic strains (G2) consisted of two groups. CH/hubei/2016 belonged to an independent subgroup which was more closely related to the US S INEDL and G2a strains than the reference strains with 3 nt deletion in S2 gene in the other subgroup.

The analysis revealed that four putatively recombinant breakpoints in CH/hubei/2016 strain were detected. One recombination event of putative breakpoint was located at position 23310 to 24012 nt on S2 gene, which was supported by 6 programs (RDP, p = 1.961E-6; GEN-ECONV, p = 1.400E-3; BootScan, p = 7.920E-5; MaxChi, p = 1.037E-3; Chimaera, p = 7.286E-5; SiScan, p = 3.351E-4; 3Seq, p = 4.564E-2). Additionally, phylogenetic trees based on the recombination region of each potential recombinant event and the non-recombinant region of each potential recombinant event displayed potential minor parent (CH-ZMDZY-11), potential major parent (AH2012-12) which had 3 ntdelection in S2 and potential recombinant strain evolutionary relationships (Fig. 4) .

A ratio of dN/dS > 0 or dN/dS < 0 indicates positive and purify in selection, respectively, and a valued dN/dS = 0 indicates neutral evolution. Significant selection is p < 0.1.The analysis revealed that PEDV strains isolated from central China from 2015 to 2016 had 7 purifying selection sites in S1; and there were 3 purifying selection sites in S2 (Fig. 5) . No significant positive selected positions were found in S gene.

The S of PDCoV CH/HN/2016 had 3480 nt which was similar with reference strains (3480 nt or 3483 nt). The S gene homology of PDCoV CH/HN/2016 strain was measured by aligning with the other S genes of PDCoV strains from USA, Korea, and Thailand. The S gene of CH/HN/ 2016 had an identity 98.3%-98.9 with most of Chinese strains, 98.2%-98.4% with USA strains, 98.0% with the CHN-AH-2004, 97.5% with Asian leopard cat coronavirus and Chinese ferret badger coronavirus and 96% with all the stains from Thailand, respectively.

The phylogenetic analysis show that all PDCoV strains fell into G1 except Thailand strains belonged to G2 (Fig. 6) . The PDCoV strains of China, USA and Korea had a closer relationship with Asian leopard cat and Chinese ferret-badger coronavirus than Thailand strains. In addition, PDCoV CH/HN/2016 from our study had a closer relationship with other PDCoV strains of China than CH/AH/2004. Further sequence analysis showed CHN/AH/2004 strain has 3-nt insertions at positions 153-155 .

Since 2012, PEDV has re-emerged and rapidly disseminated all over the world (Huang et al., 2013) . In this study, the data indicates that PEDV was more frequently found than PDCoV among the tested samples. The positive rate of PEDV in all diarrhea samples tested was 84.2%, which was similar with previous study (Li et al., 2014) . But, the positive rate of PDCoV was only 2.9% and lower than PDCoV in the south of China .

From sequence and phylogenetic analyses of PEDV S gene, the 10 PEDV strains in our study from central China were clustered into G2a and G2b. The 3 PEDV S2△strains were grouped into G2b Pandemic strains, and the remaining 7 PEDV strains were grouped into G2a Pandemic strains (Fig. 2) . In G2b, PEDV S2△strains (CH/hubei/2015, CH/hubei/2016, CH/huaiyang/2015) with a part of reference PEDVs strains possessed a same specific 3-nt deletion in S2 gene (Fig. 1) . Interestingly, phylogenetic analysis of PEDV complete genomes showed that the PEDV S2△strain (CH/hubei/2016) was not grouped into same subgroup with the other reference PEDV strains possessed S2 3 nt deletion (Fig. 3) . This phenomenon maybe explained by the fact that CH/ hubei/2016 may result from the recombination between AH2012-12 (G2b) which had 3 nt-delection in S2 and CH-ZMDZY-11(G2a) (Fig. 4) . Previous study revealed that recombination events in S1 gene between different PEDV groups may contribute to the genetic diversity of PEDV (Jarvis et al., 2016; Li et al., 2016; Sun et al., 2016) . Furthermore, Taiwan New PEDV variants have two potential recombination breakpoints in S1 and S2 (Chiou et al., 2017) . All these fact indicted that recombination, insertion and deletion in the S gene may contribute to the genetic diversity of PEDV. However, because of only one complete PEDV S2△strain genomes was included into the phylogenetic analysis, further study need to explore the evolution position of PEDV S2△strains.

Previous studies demonstrated that spike gene of Coronavirus is a determinant of virulence, and a single amino acid substitution can influence virulence (Macnamara et al., 2005; Sánchez et al., 1999) . According to previous study, the deletion and insertion in PEDV S1 gene may affect the virulence of PEDV (Lin et al., 2015) . Therefore, the specific 3 nt-deletion in PEDV S2 gene (Fig. 1 ) firstly reported in current study could also affect the pathogenicity and antigenicity of wild virus.

More effective purifying selections on RNA viruses than DNA viruses was proved and purifying selection is an important mechanism for adaption of RNA viruses (Hughes and Hughes, 2007) . Our study Fig. 4 . Detection of potential recombination events in CH hubei 2016 by RDP4. a One recombination breakpoints were located at 23310 to 24,12 nt. The analysis was performed with an RDP distance model, a window size of 1000 base pairs and a step size of 200 base pair; b Fast neighbor-joining (NJ) tree constructed using the recombinant region of each potential recombinant event; c Fast neighbor-joining (NJ) tree constructed using the non-recombinant region of each potential recombinant event. In current study, the strains with specific genetic signatures-3 nt deletion in S2 gene was firstly identified and reported (PEDV S2△strains), and one of strain (CH/hubei/2016) was clustered into different groups based on PEDV S gene or complete genomes. This fact may result from the recombination between AH2012-12 (G2b) which had 3 nt-delection in S2 and CH-ZMDZY-11(G2a). In contrast with central China of PEDV strains, PDCoV had low positive rates and had a close relationship with other PDCoV strains of China in 2012 to 2015. Future work should continue to trace the variability of PEDV strains and PDCoV strains. 

Virus like particles associated with porcine epidemic diarrhoea

Molecular epidemiology of porcine epidemic diarrhea virus in China

Two natural recombinant highly pathogenic porcine reproductive and respiratory syndrome viruses with different pathogenicities

Phylogenetic analysis of the spike (S) gene of the new variants of porcine epidemic Diarrhoea virus in Taiwan

Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States

More effective purifying selection on RNA viruses than in DNA viruses

Genomic and evolutionary inferences between American and global strains of porcine epidemic diarrhea virus

Differential detection of transmissible gastroenteritis virus and porcine epidemic diarrhea virus by duplex RT-PCR

Phylogenetic analysis of porcine epidemic diarrhea virus (PEDV) field strains in Central China based on the ORF3 gene and the main neutralization epitopes

Genome sequencing and analysis of a novel recombinant porcine epidemic diarrhea virus strain from Henan, China

Experimental infection of a US spike-insertion deletion porcine epidemic diarrhea virus in conventional nursing piglets and cross-protection to the original US PEDV infection

Origin, evolution, and virulence of porcine deltacoronaviruses in the United States

Contributions of the viral genetic background and a single amino acid substitution in an Immunodominant CD8(+) T-cell epitope to murine coronavirus Neurovirulence

Molecular evolution of porcine reproductive and respiratory syndrome virus isolates from Central China

A new coronavirus-like particle associated with diarrhea in swine

Targeted recombination demonstrates that the spike gene of transmissible gastroenteritis coronavirus is a determinant of its enteric tropism and virulence

Porcine epidemic diarrhoea virus: a comprehensive review of molecular epidemiology, diagnosis, and vaccines

Newly emerged porcine Deltacoronavirus associated with Diarrhoea in swine in China: identification, prevalence and full-length genome sequence analysis

Genomic and epidemiological characteristics provide new insights into the phylogeographical and spatiotemporal spread of porcine epidemic diarrhea virus in Asia

Epidemiology and vaccine of porcine epidemic diarrhea virus in China: a mini-review

Distinct characteristics and complex evolution of PEDV strains

Genetic properties of endemic Chinese porcine epidemic diarrhea virus strains isolated since 2010

Porcine deltacoronavirus: histological lesions and genetic characterization

Discovery of seven novel mammalian and avian coronaviruses in the genus Deltacoronavirus supports bat coronaviruses as the gene source of Alphacoronavirus and Betacoronavirus and avian coronaviruses as the gene source of Gammacoronavirus and Deltacoronavirus

An apparently new syndrome of porcine epidemic diarrhoea

Occurrence and sequence analysis of porcine deltacoronaviruses in southern China

This work was supported by grants from National Key R&D Program (2016YFD0500704 and 2017YFD0501103) and Program of Henan finance (201776-21).

The authors declare that they have no competing interests.

This article does not contain any studies with human participants or animals performed by any of the authors.

Supplementary data to this article can be found online at https:// doi.org/10.1016/j.rvsc.2018.06.001.