key: cord-0857468-276r1lh7
authors: Wang, Peng; Zhu, Jinyan; Liu, Xinze; Guo, Jiaojiao; Gu, Xuejia; Ruan, Wenke
title: Isolation and recombinant analysis of variants of porcine epidemic diarrhea virus strains from Beijing, China
date: 2019-02-09
journal: Virusdisease
DOI: 10.1007/s13337-019-00513-w
sha: dba692d157c4a975bbd9b246d54e39471746713e
doc_id: 857468
cord_uid: 276r1lh7

Porcine epidemic diarrhea virus (PEDV) is a highly infectious virus infecting pigs with high morbidity, especially for newborn piglets. Several PEDV strains were isolated from the intestinal tracts of diarrheic piglets from the Beijing area, China. Sequencing of the whole-genome of the PEDV isolates (GenBank numbers MG546687-MG546690) yielded sequences of 28033–28038 nt. The phylogenetic tree revealed that these strains from the Beijing area belonged to group II, while the vaccine strain, CV777, belonged to group I. We also determined the genetic correlation between these strains and CV777 strain. However, it showed that these strains in the Beijing area had unique mutations. The sequence identity of PEDV strains showed that these strains are most similar to these strains LZW, CH/JX-1/2013, USAIllinois972013, USAKansas1252014, CH/GDZQ/2014, SHQPYM2013, AJ1102, CHZMDZY11, KoreaK14JB01, and CHYJ130330, respectively. The possible recombination events indicate that PEDV in this studies were possibly recombinant strain formed by parent strains USAIllinois972013, KoreaK14JB01, CHYJ130330, and CHZMDZY11. These PEDV strains has been genetic recombination and mutations. The variant strains characterized in this study help to the evolutionary analysis of PEDV.

Porcine epidemic diarrhea (PED) is a highly infectious disease of swine. It causes severe watery diarrhea, vomiting, dehydration and death in suckling and nursing piglets that are typically 1-2 weeks of age [1, 2] . Porcine epidemic diarrhea virus (PEDV) is a member of the Coronaviridae family [3, 4] . The virus strain was first isolated in Belgium in 1973 and designated coronavirus CV777 [1] . There were only a few reports in the following decades [5] [6] [7] . However, a new PEDV strain has emerged in Asia since 2010 [7] [8] [9] [10] . In China, the positive rate of PEDV in some provinces was over 50% in 2010-2011 [7] . In USA, classical PEDV first identified in the state of Ohio in 2013, and a variant PEDV was reported in 2014 [11] . In Japan, more than 490,000 pigs died of PED from 2013 to 2015 [12] . In Korea, PED outbreaks reoccurred in 2013 [13, 14] . In Europe, isolates were reported from the Ukraine, Belgium, Holland, France, Germany, Portugal and Italy from 2013 to 2015 [15] . These isolates were closely related to those strains reported from North America [11, 15] .

Diarrhea disease broke out in some large-scale pig farms in the Beijing area of China from 2015 to 2016. The newborn piglets had high mortality rates, resulting in heavy loss to pig farms. The CV777 vaccines used failed. PEDV strains were isolated and identified from the small intestine of diarrhea samples obtained at these pig farms. Whole genome sequencing and recombinant analysis was performed, with the aim of examining the genomes for viral mutations and recombination. The molecular epidemiological information of PEDV can help develop efficient diagnostic reagent and vaccine.

Intestine samples of dead piglets with diarrhea were collected from several piglets infected with PEDV from 2015 to 2016. These infected piglets were chosen from 4 different large-scale swine farms located in Beijing and Hebei provinces in China. PEDV detection were positive via RT-PCR. For total RNA extraction, samples were frozen with liquid nitrogen. All procedures were approved by the Animal Care and Use Committee of Beijing University of Agriculture (Beijing, China).

The fresh intestine samples were collected under sterile conditions, made then into slurry, filter sterilization. Vero E6 cells were washed with PBS twice. After adsorption for 60 min at 37°C, PEDV growth medium was added [16] . The Vero E6 cell cultures were observed for 4 days for cytopathic effects (CPE).

Total RNA was obtained from intestine samples using TRIzol (Invitrogen, USA), per the manufacturer's instructions (Promega, USA), and was used for cDNA synthesis. PCR was performed to amplify the target gene using specific primers (reference sequence number: NC003436.1), as shown in Table 1 . PCR reaction contained primers, cDNA, dNTP mix and Hifi polymerase (Transgen, China) in 19 Hifi buffer. PCR products were sequenced from both ends by Sangon Biotech co. ltd. The coding sequences for PEDV were deposited in the Gen-Bank database under accession numbers MG546687-MG546690.

Nucleotide and amino acid sequences of PEDV were aligned by the ClustalW method. Sequence alignments of the whole genome and construction of trees were also performed by the neighbor-joining method in MEGA7 software. Information of the other strains that were chosen in this study, including years and locations of isolation and GenBank accession numbers proposed, are shown in Table 2 .

Comparison and analysis of the whole PEDV genome sequences and the reference PEDV strains were conducted using MegAlign programs (DNAStar, Madison, WI).

The putative recombinant sequence of these strains in this study and its parental strains were conducted using recombination detection program (RDP) version 4.97 with seven recombination detection methods (RDP, GENE-CONV, Chimaera, MaxChi, BootScan, SiScan, and 3Seq). In order to improve the reliability of detection, the recombination breakpoints through at least 5 of these methods was taken as confirmatory for any putative recombination event (p value is 0.01). The phylogenetic tree of recombination and non-recombination were constructed to verify the reliability of these results.

Vero E6 cells were used for viral multiplication. Cytopathic effects (CPE) appeared 48 h after inoculating virus, and more than 70% CPE with cell clustering, detachment, and syncytium formation were observed at 72 h. These results indicated that the isolated PEDV could cause pathological changes in Vero E6 cells (Fig. 1 ).

The organization of the genome of PEDV was characteristic of the gene order 5 0 -ORF1a/1b-S-ORF3-E-M-N-3 0 [2, 4, 17, 18] . Sequencing of the complete gene of the four isolates (GenBank numbers: MG546687-MG546690) yielded sequences of 28033-28038 bp, which compared with 18 typical strains available in GenBank ( Table 2 ). The phylogenetic tree showed that these 22 PEDV strains divided into two groups (Fig. 2 ).

The sequence identity between PEDV strains in this studies and typical strains in terms of their according genes or regions were calculated. The highest percentage of nucleotide identities of 5 0 -UTR, ORF1a/1b, S, M, N, E, ORF3, and 3 0 -UTR of PEDV compared with other PEDV strains were show in Table 3 . For ORF1a/1b gene, these strains are most similar to these strains, LZW, SHQPYM2013, AJ1102, CH/YNKM-8/2013, CH/JX-1/ 2013, USAKansas1252014, and CHZMDZY11, respectively. For S gene, these strains are most similar to these strains LZW, KoreaK14JB01, USAIllinois972013, USA-Kansas1252014, and CH/GDZQ/2014, respectively. The percentage of nucleotide identities of S of PEDV strains in this study compared with recombined strains SeCoV is 90.4-90.7%. For M gene, these strains are most similar to these strains USAIllinois972013 and USAKansas1252014. (Table 3) . The names of the strains, years and places of isolation, as well as GenBank accession numbers proposed, are shown in Table 2 . The strains in this study are indicated by underline Isolation and recombinant analysis of variants of porcine epidemic diarrhea virus strains from… 297 Frequent mutation and recombination of genes may occur in epidemics of PEDV. In order to detect these possible recombination events, 19 reference strains were used as parent strains for RDP, GENECONV, Chimaera, MaxChi, BootScan, SiScan, and 3Seq program analysis. The results showed that PEDV in this study were possibly recombinant strain formed by parent strains USAIlli-nois972013, KoreaK14JB01, CHYJ130330 and CHZMDZY11 (Fig. 3) . USAIllinois972013 is major parent of PEDV strain CH/BJ9/2015, the recombination region locates on 4669-8176nt; KoreaK14JB01 is minor parent of PEDV strains CH/BJ2/2016 and CH/BJ11/2016, the recombination region locates on 21744-27721nt; CHYJ130330 is major parent of PEDV strain CH/HB9/ 2016 and CHZMDZY11 is minor parent of CH/HB9/2016, their recombination region locates on 2218-15942nt. However, there may be a deviation in the prediction of software. It requires further confirms these recombination events.

Since 2011, a large number of piglets died as a result of PEDV outbreaks, causing much economic loss in China [2-4, 7-10, 18-21] . In this study, the deceased piglets that provided intestine samples had been immunized with the CV777 vaccine previously. This showed that the classic vaccine strain exhibited lower protective efficacy against the variant PEDV [22, 23] . The phylogenetic trees showed that, these PEDV isolates from the Beijing area belonged to group II, while the traditional vaccine strain, CV777, belonged to group I based on the comparison of several important genes in PEDV, such as ORF1a/1b, S, M, N, E, and ORF3. However, N, as a conserved gene, was inconsistent with the above results, which also suggested that the N sequence or protein is suitable for universal detection. The ORF1a/1b gene were two long ORFs which encoding the non-structural replicase polyproteins of PEDV. The strain from 2015 in the Beijing area shared the greatest homology in ORF1a/1b with LZW, which was also isolated in Beijing; it showed the continuity of heredity. The S sequence of isolates in 2015 in the Beijing area showed a closer relationship with a previous isolate, LZW. Other isolates in 2016 in the Beijing area showed a closer relationship with previous isolates in Guangdong. This revealed that S genes accumulated variations faster based on inherent genetic characteristics or recombinant genes from other isolates [10, 20, 24] .

The 1-1180 bp section in the S1 region of PEDV genes is prone to mutate [24] . Sequence analysis revealed that isolated PEDV strains and these variant PEDV strains were reported had insertion in 167-178 sites (GTGAAAAC-CAGG) and 421-423 sites (TGA), also deletion in 218 (A), 471-474 (ATAT) and 481-482 (TG) comparing with CV777.

These insertion, deletion, or introducing stop codon (TGA) of sequence resulted in frameshift mutation and abnormally transcripts of S1 gene. It is a reason for failed vaccination with CV777 vaccine. However, different from other recent isolates, these strains in the Beijing area had shared sites with CV777 in the S1 region, such as A5G (this specific nucleotide change led to the amino acid change K-R), C70T, and T1043C (which led to the amino acid change F-S). This determined the genetic correlation between these strains in the Beijing area and the CV777 strain. It was remarkable that E sequences of isolates in the Beijing area had distinct differences from other isolates. There were unique variations in the E gene in these strains, such as T27G. This showed that these isolates from the Beijing area had unique mutations. This might provide a gene reference for other mutations of PEDV in the future.

Viral gene recombination and mutation results in evolution [25] . Gene recombination and mutation often occurs in coronavirus, such as PEDV, SARSCoV, and infectious bronchitis virus (IBV) [26] [27] [28] [29] . There is evidence that the S gene of epidemic PEDV strains and genome of transmissible gastroenteritis virus (TEGV) strains have recombined when both were present in the same host to form a new recombinant porcine coronavirus, SeCoV [30, 31] . This S sequence of SeCoV strain share above 90% homology with S sequence of PEDV strains in this study. However, analysis of the recombination of PEDV strains showed that, no recombination event was occurred between SeCoV strain and PEDV strains in this study. These PEDV strains isolated in Beijing area shows recombination events with PEDV isolated in USA, PEDV isolated in Korea, and PEDV isolated in China in sequence. It suggests the accelerated spread of PEDV virus.

Viral recombination and mutation causes immune failure. The new and effective vaccine is critical to control the occurrence of PEDV. The study of the effective vaccine must be based on the characteristic analysis of variant viruses [32] . The new data on the evolutionary analysis of PEDV is provided in this study. These data will be helpful for research on PEDV vaccine and control.

Several PEDV strains were isolated from the intestinal tracts of diarrheic piglets from the Beijing area, China. These PEDV strains has been genetic recombination and mutations. The variant strains characterized in this study help to the evolutionary analysis of PEDV.

A new coronavirus-like particle associated with diarrhea in swine

Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea virus (PEDV) field strains in south China

Sequence and phylogenetic analysis of nucleocapsid genes of porcine epidemic diarrhea virus (PEDV) strains in China

Whole-genome analysis of porcine epidemic diarrhea virus (PEDV) from eastern China

An outbreak of swine diarrhea of a new-type associated with coronavirus-like particles in Japan

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

Porcine epidemic diarrhea in China

Isolation and characterization of a variant porcine epidemic diarrhea virus in China

Genetic variability and phylogeny of current chinese porcine epidemic diarrhea virus strains based on spike, ORF3, and membrane genes

Bioinformatics insight into the spike glycoprotein gene of field porcine epidemic diarrhea strains during 2011-2013 in Guangdong, China

New variant of porcine epidemic diarrhea virus

Isolation and molecular characterization of porcine epidemic diarrhea viruses collected in Japan in 2014

New emergence pattern with variant porcine epidemic diarrhea viruses

Outbreak-related porcine epidemic diarrhea virus strains similar to US strains

Porcine epidemic diarrhea: a retrospect from Europe and matters of debate

Cell culture isolation and sequence analysis of genetically diverse US porcine epidemic diarrhea virus strains including a novel strain with a large deletion in the spike gene

Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea virus (PEDV) field isolates in Korea

Isolation, molecular characterization and an artificial infection model for a variant porcine epidemic diarrhea virus strain from Jiangsu Province, China

Evolutionary characterization of the emerging porcine epidemic diarrhea virus worldwide and 2014 epidemic in Taiwan

Detection and phylogenetic analysis of porcine epidemic diarrhea virus in central China based on the ORF3 gene and the S1 gene

Characterization of Chinese porcine epidemic diarrhea virus with novel insertions and deletions in genome

Porcine epidemic diarrhea: a review of current epidemiology and available vaccines

Cross protective immune responses in nursing piglets infected with a US spike-insertion deletion porcine epidemic diarrhea virus strain and challenged with an original US PEDV strain

S1 domain of the porcine epidemic diarrhea virus spike protein as a vaccine antigen

Epidemiology, genetic recombination, and pathogenesis of coronaviruses

Recombination in avian gamma-coronavirus infectious bronchitis virus

Evolution of infectious bronchitis virus in China over the past two decades

Evolution and variation of the SARS-CoV genome

Viral evolution and the emergence of SARS coronavirus

Porcine epidemic diarrhea virus and discovery of a recombinant swine enteric coronavirus

Characterization of a novel chimeric swine enteric coronavirus from diseased pigs in central eastern Europe in 2016

Evaluation of a recombination-resistant coronavirus as a broadly applicable, rapidly implementable vaccine platform

Acknowledgements The authors would like to thank the pig farm owners their collaboration.Funding This work was supported by the National Natural Science Foundation of China (No. 31572499). The funders did not play any role in the design, conclusions or interpretation of the study.

The authors declare that they have no competing interests.Ethics approval and consent to participate Experimental animals used in the study were cared for in accordance with the internationally accepted standards for the care and handling of laboratory animals. All procedures were approved by the Animal Care and Use Committee of Beijing University of Agriculture (Beijing, China). The Animal Ethics Committee approval number was BUAACUC2016-1. All samples were collected after informed written consent was obtained from the owners.