Printer Friendly

Evidence of recombinant strains of porcine epidemic diarrhea virus, United States, 2013.

Porcine epidemic diarrhea virus (PEDV) is an alphacoronavirus that causes enteric disease in swine. The disease, PED, is characterized by acute vomiting and watery diarrhea and causes high mortality rates in newborn piglets (1). PED was first reported in 1971 in the United Kingdom and was soon identified in many European and Asian countries (1,2). Variant PEDV strains that were fatal to young pigs, initially isolated during late 2010 in China and Southeast Asia (3,4) and in May 2013 in the United States (5,6), have posed a serious threat to the pork industry. Incidence of PEDV-associated large-scale outbreaks of diarrheal disease in China was reported at 80%-100% in suckling piglets (3,7) and outbreaks in the United States had spread to 25 states by February 2014 (http://www.aasv. org/news/story.php?id=7038), causing numerous deaths in neonatal piglets (5,6,8,9).

How the virus entered the United States remains unknown. A phylogenetic analysis based on available fulllength genomic PEDV sequences indicated that all PEDV strains were classified into 2 distinct genogroups: G1 and G2 (6). PEDV field strains isolated before 2010 and the derived vaccine strains were in the G1 genogroup, whereas all the new PEDV strains isolated since 2011 in China and the United States (US PEDV) are in G2. The US PEDV sequences were >99% identical to strains found in China in the subgroup G2a, suggesting their origin. In particular, the US PEDV are most closely related to strain AH2012, which was isolated in eastern China and was proposed to have come from multiple recombination events among G2 lineages of PEDV (6). Divergence of PEDV is driven by genetic recombination, as in other coronaviruses (10). Details of recombination events in the process are needed to investigate origins. To investigate the evolutionary process by which US PEDV strains hypothetically descended from precursors in China, we conducted a molecular epidemiologic analysis using PEDV-positive samples collected from eastern China since 2012.

The Study

A total of 169 fecal and intestinal samples were collected from pigs with typical PED symptoms on 26 farms in 4 provinces of eastern China during January 2012-July 2013. The rate of PEDV-positive samples was 56.8% (96/169) as had been determined by using reverse transcription PCR (RT-PCR) specific for the spike (S) gene (11). From the positive samples, we selected 24 representative samples (Table) to examine. Using RT-PCR, we determined the sequences of the full-length genomic cDNA for the strain CH/ZJCX-1/2012, identified the spike (S) gene for strains CH/ ZJQZ-2w/2012 and CH/ZJDX-1/2012, and identified the region encoding structural protein genes by an order of 5'-S-ORF3-E-M-N-3' (5'-spike protein-open reading frame 3-envelope-membrane-nucleoprotein-3') for the remaining 21 strains. All primers were designed based on the PEDV MN strain (GenBank accession no. KF468752). We purified and cloned PCR products into a vector using TA cloning. We used Vector NTI software (http://www. vector-nti-software.html) to assemble and analyze the sequences. We performed multiple alignments of S-ORF3-E-M-N, S, ORF3, M, N, and full-length genomes with available sequences from Asia and the United States (5,6,8,12) and performed phylogenetic analyses using the MEGA5.2 program ( with the neighbor-joining method.

Similar to most of the sequences recently documented in PEDV strains in China and the United States, the S genes of the 24 samples have a 4,161-nt sequence that, compared with the prototype CV777 strain, shows 97.9%-100% sequence identities, and contain 2 notable insertions at amino acids (aa) 56-59 (IGEN) and 139 (N) and a deletion of 2 aa (GK) at aa positions 160 and 161 at the N terminus (6). A phylogenetic analysis comparing S genes showed that, based on the complete genome (online Technical Appendix Figure 1,, all 24 strains were classified into the same group corresponding to G2. However, it is notable that the Chinese sublineage (branch) most closely related to the US PEDV strains did not include the AH2012 strain. Instead, this sublineage contained the strain CH/ZMDZY/11 and 4 other strains determined in this study. Analyses of the phylogenetic trees constructed on the basis of the S-ORF3-E-M-N genes (Figure), ORF3 or M (online Technical Appendix Figure 1) also indicated that the AH2012 strain was not closely related to the US branch, relative to the sublineage represented by strains CH/ZMDZY/11, CH/ HuBWHYQ/2012, CH/JXZS-1223L/2012, and CH/JXZS-3L/2012 (designated the ZMDZY sublineage hereafter; Figure). The exception is the N gene-based tree, in which the AH2012 was grouped more closely to the sublineage associated with the United States than the strains in the designated ZMDZY sublineage (online Technical Appendix Figure 1).

The relationship of the AH2012 strain with the 33 PEDV strains identified in China, the United States, South Korea, and Belgium (online Technical Appendix Figure 1) in nonstructural protein genes was also determined by generation of 3 phylogenetic trees based on ORF1ab, ORF1a, and ORF1b genes, respectively. In accordance with the results from the N gene and the PEDV genotyping based on the full-length genomes (6), the AH2012 strain in these trees was most closely related to the US strains (online Technical Appendix Figure 1). Therefore, the strains AH2012 and CH/ZMDZY/11 displayed different phylogenetic relationships in different genome regions. Overall, the AH2012 strain was clustered closely with the US strains in the ORF1ab and the N gene region, whereas the ZMDZY sublineage was clustered closely with the US strains in the S-ORF3-E-M region.

To accurately determine how the US strains are related to strains AH2012 and the ZMDZY sublineage, we performed a recombination analysis using the Recombination Detection Program and available PEDV sequences (13). We used a multiple-comparison-corrected p-value cutoff of 0.05 in all methods. Recombination events were identified by 6 methods (Recombination Detection Program, GENECONV, BOOTSCAN, MaxChi, CHIMAERA, and SISCAN) when the US PEDV sublineage represented by the MN strain was used as a query. By bootstrap analysis, 3 putatively major recombination breakpoints were detected at nucleotides 6699, 21840, and 26882 (online Technical Appendix Figure 2), which generated 2 regions: 1 covered the 3' half of ORF1a, complete ORF1b, and the N terminus of the S (first 1,207 nt); the other spanned partial S, ORF3, E, M, and partial N (first 504 nt) between the strain AH2012 (as the major parent) and the ZMDZY sublineage (as the minor parent). Although the second region (partial S-ORF3-E-M-partial N) of the US strains is associated with the ZMDZY sublineage, the source of genetic material in this region is not known, because none of the PEDV strains in this sublineage had a highly identical sequence to the consensus sequence of the US strains. It is possible that the other recombination breakpoints exist within the S gene, according to the bootstrap supports in this region, which may be determined by future study with available new sequence data. We showed that the emergent US PEDV strains are possibly descendent of 2 major lineages derived from the ZMDZY sublineage and AH2012 through recombination.


Our study provides further information on the origin of the US PEDV in 2013. We identified 21 S-ORF3-E-M-N genes, 2 S genes, and 1 full-length genomic cDNA of PEDV from PEDV-positive samples collected in eastern China. Comparative genomic, phylogenetic, and recombination analyses using new and known sequence data demonstrated that the AH2012 strain is likely not the direct progenitor of emergent US PEDV strains during 2013. It is possible that replacement of a region within the partial S-ORF3-E-M-partial N region of the AH2012 strain with a corresponding fragment close to the ZMDZY sublineage (including several newly identified strains) resulted in a recombinant strain related to emergence of this virus in swine in the United States. Other unidentified recombination events and accumulation of adapted mutations within the structural protein genes were also likely involved in this process.

This work was supported by the Innovation Team Projects of Zhejiang province (No. 2010R50031) and the Fundamental Research Funds for the Central Universities of China.

Mr Tian is a graduate student in the veterinary science program at Zhejiang University. His research interests include molecular epidemiology and pathology of emerging animal viruses.


(1.) Song D, Park B. Porcine epidemic diarrhoea virus: a comprehensive review of molecular epidemiology, diagnosis, and vaccines. Virus Genes. 2012;44:167-75.

(2.) Pensaert MB, de Bouck P A new coronavirus-like particle associated with diarrhea in swine. Arch Virol. 1978;58:243-7. http://dx.doi. org/10.1007/BF01317606

(3.) Pan Y, Tian X, Li W, Zhou Q, Wang D, Bi Y, et al. Isolation and characterization of a variant porcine epidemic diarrhea virus in China. Virol J. 2012;9:195.

(4.) Puranaveja S, Poolperm P, Lertwatcharasarakul P, Kesdaengsakonwut S, Boonsoongnern A, Urairong K, et al. Chinese-like strain of porcine epidemic diarrhea virus, Thailand. Emerg Infect Dis. 2009;15:1112-5.

(5.) Marthaler D, Jiang Y, Otterson T, Goyal S, Rossow K, Collins J. Complete genome sequence of porcine epidemic diarrhea virus strain USA/Colorado/2013 from the United States. Genome Announc. 2013;1:e00555.

(6.) Huang YW, Dickerman AW, Pineyro P, Li L, Fang L, Kiehne R, et al. Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States. MBio. 2013;4:e00737.

(7.) Sun RQ, Cai RJ, Chen YQ, Liang PS, Chen DK, Song CX. Outbreak of porcine epidemic diarrhea in suckling piglets, China. Emerg Infect Dis. 2012;18:161-3.

(8.) Stevenson GW, Hoang H, Schwartz KJ, Burrough ER, Sun D, Madson D, et al. Emergence of porcine epidemic diarrhea virus in the United States: clinical signs, lesions, and viral genomic sequences. J Vet Diagn Invest. 2013;25:649-54. http://dx.doi. org/10.1177/1040638713501675

(9.) Jung K, Wang Q, Scheuer KA, Lu Z, Zhang Y, Saif LJ. Pathology of US porcine epidemic diarrhea virus strain PC21A in gnotobiotic pigs. Emerg Infect Dis. 2014.

(10.) Graham RL, Baric RS. Recombination, reservoirs, and the modular spike: mechanisms of coronavirus cross-species transmission. J Virol. 2010;84:3134-46.

(11.) Song DS, Kang BK, Oh JS, Ha GW, Yang JS, Moon HJ, et al. Multiplex reverse transcription-PCR for rapid differential detection of porcine epidemic diarrhea virus, transmissible gastroenteritis virus, and porcine group A rotavirus. J Vet Diagn Invest. 2006;18:278-81.

(12.) Chen Q, Li G, Stasko J, Thomas JT, Stensland WR, Pillatzki AE, et al. Isolation and characterization of porcine epidemic diarrhea viruses associated with the 2013 disease outbreak among swine in the United States. J Clin Microbiol. 2014;52:234-43. http://dx.doi. org/10.1128/JCM.02820-13

(13.) Martin DP. Recombination detection and analysis using RDP3. Methods Mol Biol. 2009;537:185-205. http://dx.doi. org/10.1007/978-1-59745-251-9_9

Address for correspondence: Ji-Yong Zhou or Yao-Wei Huang, Key Laboratory of Animal Virology of Ministry of Agriculture, Zhejiang University, 866 Yuhangtang Road, Hangzhou 310058, PR China; email: and

Peng-Fei Tian, Yu-Lan Jin, Gang Xing, Ling-Ling Qv, Yao-Wei Huang, [1] and Ji-Yong Zhou [1]

[1] These authors were co-principal investigators.

Author affiliation: Key Laboratory of Animal Virology of Ministry of Agriculture, Zhejiang University, Hangzhou, China


Table. Summary of 24 representative porcine epidemic diarrheal
virus sequences in China determined in this study, 2012-2013

Strain/year          Province of    Sequencing       GenBank
                     collection      region *     accession no.

CH/ZJCX-1/2012        Zhejiang     Full-length      KF840537
CH/ZJHY-2/2012        Zhejiang     S-ORF3-E-M-N     KF840538
CH/ZJJS-Z/2012        Zhejiang     S-ORF3-E-M-N     KF840539
CH/ZJJS-1Z/2012       Zhejiang     S-ORF3-E-M-N     KF840540
CH/JXZS-3H/2012        Jiangxi     S-ORF3-E-M-N     KF840541
CH/JXZS-2H/2012        Jiangxi     S-ORF3-E-M-N     KF840542
CH/ZJXS212/2012       Zhejiang     S-ORF3-E-M-N     KF840543
CH/ZJHZHY-6/2013      Zhejiang     S-ORF3-E-M-N     KF840544
CH/JXJDZ-F/2012        Jiangxi     S-ORF3-E-M-N     KF840545
CH/ZJJS-2Z/2012       Zhejiang     S-ORF3-E-M-N     KF840546
CH/JXZS-3L/2012        Jiangxi     S-ORF3-E-M-N     KF840547
CH/JXZS-1223L/2012     Jiangxi     S-ORF3-E-M-N     KF840548
CH/SDZD-1/2012        Shandong     S-ORF3-E-M-N     KF840549
CH/SDZD-2/2012        Shandong     S-ORF3-E-M-N     KF840850
CH/HuBWHY Q/2012        Hubei      S-ORF3-E-M-N     KF840851
CH/ZJQZ-2/2012        Zhejiang     S-ORF3-E-M-N     KF840852
CH/ZJHZ-1C/2012       Zhejiang     S-ORF3-E-M-N     KF840853
CH/ZJHZ-2C/2012       Zhejiang     S-ORF3-E-M-N     KF840854
CH/JXJDZ-1/2012        Jiangxi     S-ORF3-E-M-N     KF840855
CH/ZJJS-4X/2012       Zhejiang     S-ORF3-E-M-N     KF840856
CH/ZJQZ-2w/2012       Zhejiang          S           KF840857
CH/ZJDX-1/2012        Zhejiang          S           KF840858
CH/JSZL-2/2013         Jiangxi     S-ORF3-E-M-N     KF840861
CH/JSZL-3/2013         Jiangxi     S-ORF3-E-M-N     KF840862

* S, spike protein; ORF, open reading frame; E, envelope;
M, membrane; N, nucleoprotein.
COPYRIGHT 2014 U.S. National Center for Infectious Diseases
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:DISPATCHES
Author:Tian, Peng-Fei; Jin, Yu-Lan; Xing, Gang; Qv, Ling-Ling; Huang, Yao-Wei; Zhou, Ji-Yong
Publication:Emerging Infectious Diseases
Geographic Code:9CHIN
Date:Oct 1, 2014
Previous Article:Treatment of giardiasis after nonresponse to nitroimidazole.
Next Article:Human Babesiosis, Maine, USA, 1995-2011.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters