key: cord-0736706-meemq3i9
authors: Song, D.; Zhou, X.; Peng, Q.; Chen, Y.; Zhang, F.; Huang, T.; Zhang, T.; Li, A.; Huang, D.; Wu, Q.; He, H.; Tang, Y.
title: Newly Emerged Porcine Deltacoronavirus Associated With Diarrhoea in Swine in China: Identification, Prevalence and Full‐Length Genome Sequence Analysis
date: 2015-08-06
journal: Transbound Emerg Dis
DOI: 10.1111/tbed.12399
sha: 5984ceec2234a1b7718a425ec9b2045148c922cc
doc_id: 736706
cord_uid: meemq3i9

To identify and characterize aetiologic agent(s) associated with an outbreak of a severe diarrhoea in piglets in Jiangxi, China, in March 2015, a nested reverse transcription–polymerase chain reaction (RT‐PCR) for the detection of porcine deltacoronavirus (PDCoV) was developed. A survey based on the nested RT‐PCR established indicated that the monoinfection of PDCoV (33.71%) and coinfection of PDCoV (19.66%) with porcine epidemic diarrhoea virus (PEDV) were common in diarrhoeal pigs in Jiangxi, China. A high prevalence of PDCoV (58.33%) in diarrhoeal samples which were PEDV negative was observed. The complete genome sequence of a representative PDCoV strain, PDCoV/CHJXNI2/2015, was determined. Phylogenetic analysis of complete genome and S protein sequences of PDCoV/CHJXNI2/2015 demonstrated that it was most closely related to Hong Kong and US PDCoVs. To our knowledge, this is the first report on the identification, prevalence, complete genome sequencing and molecular characterizations of PDCoV in diarrhoeal samples in pigs in China.

Coronaviruses are traditionally divided into three genera, Alphacoronavirus, Betacoronavirus and Gammacoronavirus (Woo et al., 2010) . Recently, a new genus, Deltacoronavirus, was first discovered in birds from Hong Kong and identified in various host species, including Asian leopard cats, and some avian species (Dong et al., 2007; Woo et al., 2009; Chan et al., 2013) . Porcine deltacoronavirus (PDCoV) was first documented as a novel coronavirus in pigs from rectal swabs in Hong Kong 2012 (Woo et al., 2012) , and then recognized and isolated in USA (Li et al., 2014; Marthaler et al., 2014a; Wang et al., 2014a,b; Hu et al., 2015) and Korea (Lee and Lee, 2014) . Clinical reports indicated that PDCoV could cause severe diarrhoea, vomiting, and dehydration in piglets, symptomatically indistinguishable from those caused by porcine epidemic diarrhoea virus (PEDV) and transmissible gastroenteritis virus (TGEV) (Wang et al., 2014a) . Experimental studies on conventional and gnotobiotic piglets showed that isolated PDCoVs caused similar clinical symptoms to those of field infections, from mild to severe diarrhoea and intestinal lesions Jung et al., 2015; Ma et al., 2015) . Molecular surveillance on diarrhoeal samples of pigs from USA indicated a 30% infection of PDCoV, which was considered as a common viral pathogen of pigs. Coinfection of PDCoV/porcine rotavirus (PoRV) was most common in pig herds in USA (Marthaler et al., 2014b) .

The full-length genome of PDCoV is about 25 kb, arranged in the order of: 5 0 untranslated region (UTR), open reading frame 1a and 1b, spike (S), envelope (E), membrane (M), non-structural protein 6 (NS 6), nucleocapsid (N), non-structural protein 7 (NS 7) and 3 0 UTR (Woo et al., 2012; Lee and Lee, 2014) . Up to date, PDCoV has been reported in Hong Kong, USA, Korea and Canada, but little is known regarding PDCoV infection in pigs in mainland China, a major pig-raising country being suffered with swine diarrhoea for years. Here, we report the identification, prevalence and full-length genome sequence analysis of newly emerged PDCoVs associated with diarrhoea in pigs in Jiangxi province, China.

Clinical signs and diagnostics In March 2015, there was an acute diarrhoea outbreak in a pig farm (Farm A), a small-scale independent commercial farm with about 100 sows, in Jiangxi, China. The clinical manifestations of the disease were characterized with profuse watery diarrhoea, vomiting, dehydration and a high mortality (>80%) in suckling piglets, resembling the diseases caused by PEDV and/or TGEV. Intestinal/faecal samples (n = 31) from diseased pigs were submitted to our laboratory for identification of aetiologic agent(s). Initially, common enteric viral and bacterial pathogens, including PEDV, TGEV, PoRV and Escherichia coli, were examined with standard assays. The results indicated that only 10 of 31 (32.25%) samples tested were positive for PEDV, but negative for the rest enteric pathogens aforementioned. Due to the features of porcine epidemic diarrhoea (PED)-like of the disease, we hypothesized that PDCoV might be a contributor for the diarrhoea outbreak in pigs and attempted to identify and characterize the potential PDCoV.

To detect PDCoVs from diarrhoeal faecal/intestinal samples of pigs, a nested reverse transcription-polymerase chain reaction (RT-PCR) was developed with the primers designed based on the nucleocapsid gene of PDCoV HKU15 strain IN2847 (GenBank accession no. KJ569769, Wang et al., 2014c) : 5 0 -outer primer (5 0 -TGCTACC TCTCCGATTCCCA-3 0 ), 3 0 -outer primer (5 0 -ATCCTGTT TGTCTGCTGGCA-3 0 ); 5 0 -inner primer (5 0 -GACACTGAG AAGACGGGTATGG-3 0 ), and 3 0 -inner primer 5 0 -TAGTT GGTTTGGTAGGTGGCTC-3 0 ), and the expected sizes of amplicons from the primers are 614 and 238 bp, respectively. Viral RNAs were extracted from faeces and/or intestinal contents by RNAplus Reagent (TaKaRa, Dalian, China) following the manufacturer's instructions, and then, the concentrations of RNAs extracted were measured by NanoDrop 2000 spectrophotometer (Thermo scientific, Wilmington, DE, USA). The first-strand cDNA was synthesized with 3 0 -outer primer according to the protocol reported previously (Song et al., 2015) . For nested PCR, initial amplification was executed with 5 0 -outer primer and 3 0 -outer primer under the following conditions: denaturation at 95°C for 4 min, 20 cycles of 94°C 9 30 s, 55°C 9 30 s, 72°C 9 45 s and consequently with a final extension at 72°C for 7 min. After that, 1 ll of initial PCR product was used as a template in the next 25 cycling PCR with 5 0 -inner primer and 3 0 -inner primer. Expected PCR products were purified, cloned and then sequenced based on the procedures (Song et al., 2015) .

To address the prevalence of PDCoV monoinfection and coinfection(s) with other enteric viral agent(s) in diarrhoeal specimens of pigs, a survey was carried out by utilizing RT-PCRs established in our laboratory for detection of PDCoV, PEDV, TGEV and PoRV. Faecal samples (n = 170: 55 from sows and 115 from piglets less than 10 days old) and intestinal samples (n = 186: all from dead piglets less than 10 days old) collected from November 2012 to March 2015 from diarrhoeal pigs on 51 pig farms in Jiangxi province, China, were examined using the nested RT-PCR established in this study.

Primers for amplification the complete genome sequence of PDCoV field strains were designed based on PDCoV HKU15 strain IN2847 (GenBank accession no. KJ569769) ( Table 1 ). Viral RNA was extracted by the method aforementioned. Fragments covering the complete genome of PDCoV were amplified on the conditions of a denaturation at 95°C for 4 min, 36 cycles (94°C 9 45 s, 53°C 9 45 s, and 72°C 9 2.5 min), and then with a final extension of 72°C for 10 min. The 5 0 -and 3 0 -RACE for the determination of the terminal sequences were performed using 5 0 /3 0 SMARTer RACE kit (Clontech, Beijing, China) according to the manufacturer 0 s instructions. Polymerase chain reaction products obtained were subject to gel purification using a gel extraction kit (TaKaRa), and afterwards cloned into pMD 18-T vectors (TaKaRa) following the manufacturer's protocol. Three to five positive clones of each amplicon were submitted to a commercial company (Sangon Biotech, Shanghai, China) for sequencing at both directions by Sanger sequencing methodology. The raw sequence fragments targeting PDCoV were imported to SeqMan in DNAStar Lasergene V 7.10 (DNAStar, Inc., Madison, WI, USA) for assembly and annotation. Phylogenetic trees based on complete genome and deduced amino acid (aa) sequences of the S protein of PDCoVs were constructed by the neighbour-joining method (Bootstrap in 1000 replicates) using MEGA 6 program (Tamura et al., 2013) .

Nested RT-PCR establishment A nested RT-PCR for detection of PDCoV was developed and firstly employed to test 31 diarrhoeal samples of pigs received by our laboratory. Of which, nine of these samples showed expected RT-PCR products, corresponding sequences via cloning and sequencing. The results of BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) search of these fragments amplified by the nested RT-PCR showed that the sequences obtained had 99% nucleotide (nt) identity with PDCoV HKU15 strain IN2847, indicating that the virus detected belonged to PDCoV.

Clinical samples collected from pigs with diarrhoea from 51 farms in Jiangxi, China, from November 2012 to March 2015 were selected for PDCoV prevalence survey. Of 356 porcine faecal/intestinal samples examined, 120 (33.71%) were PDCoV positive; 231 (64.89%) were PEDV positive; 281(78.93%) were positive either for PEDV or PDCoV and none of these samples were positive for TGEV and PRoV ( Table 2) ; 50 of 125 (40%) samples which were negative for PEDV were confirmed PDCoV positive. PDCoV/PEDV coinfections were most common; and 70 of 120 (58.33%) PDCoV-positive samples were also proved to be positive for PEDV. As for sample type, PDCoV was detected in 37.63% (70/186) of intestinal samples and 29.41% (50/170) in faecal samples, respectively. Phylogenetic trees were constructed based on the complete genome and aa sequences of the S gene of PDCoV/ CHJXNI2/2015, deltacoronavirus strains from pigs and other hosts, and coronaviruses from genera of Alphacoronavirus, Betacoronavirus and Gammacoronavirus. In the context of the full-length genome phylogenetic tree, all 23 deltacoronavirus strains were clustered into a group; and PDCoV/CHJXNI2/2015 was grouped with other 18 swineorigin deltacoronaviruses, forming a branch, which was distinct from the bird-origin deltacoronaviruses (Fig. 1, panel  a) . The topology of the S protein resembled that of fulllength genome sequences (Fig. 1, panel b ).

As a major pig-raising country in the world, pigs in China have been suffering severe diarrhoea for years, but only about 60% of the diarrhoeal samples were confirmed as infection of PEDV (Sun et al., 2012; Song et al., 2015) , which indicated other infectious/non-infectious factors might involved in the outbreaks of diarrhoea of pigs. Deltacoronavirus has been detected in pigs in several counties, and the pathogenicity and pathogenesis of this virus have been elucidated in gnotobiotic and conventional neonatal piglets Jung et al., 2015; Ma et al., 2015) . However, until now, little regarding PDCoV in China has been reported, and it is urgently needed to develop assays for detection of PDCoV, a newly emerged Of 56 samples, 31 were from Farm A (10 were PEDV positive, 9 were PDCoV positive and 6 were PEDV/PDCoV positive).

diarrhoeal virus in China, which might provide insight into the diarrhoea outbreaks in pigs.

In this study, a nested RT-PCR for detection of the newly emerged PDCoV was established and evaluated. High prevalence of PDCoV infections in diarrhoeal samples of pigs in Jiangxi, China, was confirmed. A high mortality in piglets in an outbreak with severe diarrhoea in Farm A was observed in this study. A survey based on RT-PCRs showed that the diarrhoeal samples of pigs from November 2012 to March 2015 examined revealed the infection of PDCoV was 33.71%, and coinfection of PDCoV/PEDV was 19.66%, respectively. We speculate that PDCoV was a neglected pathogen associated with swine diarrhoea in China for years and hypothesize that PDCoV/PEDV coinfection might result in the high mortality in diarrhoeal piglets. Further study is needed to elucidate the interactions between these two viruses. The PDCoV infection rate (33.71%) from the survey in Jiangxi, China, was consistent with the previous report (Marthaler et al., 2014a,b) . However, PDCoV/PEDV coinfection (19.66%) was most common in Jiangxi pigs examined, rather than PDCoV/PoRV coinfection (58%) in US pig herds tested (Marthaler et al., 2014b) .

Analysis based on complete genome sequences showed that PDCoV/CHJXNI2/2015 strain was most closely related to Hong Kong HKU 15-155 and 5 US PDCoV strains, which was similar to the previous reports (Wang et al., 2014a; Ma et al., 2015) . Like HKU 15-155 strain, PDCoV/ CHJXNI2/2015 contained a 3-nt deletion in S gene which was not present in the US and Korea PDCoV strains. More information regarding the genome and molecular epidemiology of PDCoVs is needed to extend our knowledge on the evolution of China-origin PDCoVs. Moreover, the isolation of PDCoV is urgently needed to help better understand PDCoV's virology.

In conclusion, a nested RT-PCR was developed and evaluated for detection of PDCoV in field samples from pigs. Survey on the diarrhoeal samples showed that PDCoV monoinfection and PDCoV/PEDV coinfections in diarrhoeal swine in Jiangxi, China, were frequently detected. Full-length genome of PDCoV/CHJXNI2/2015, a newly emerged China PDCoV strain, was determined, and phylogenetic analysis indicated it had a close relationship with Hong Kong and US PDCoV strains. . The trees were constructed by MEGA 6.05, using the neighbour-joining method. Bootstrap was set in 1000 replicates. Bars of 2/0.2 indicate nt or aa substitutions per site. Names of each reference strain retrieved from GenBank are followed by corresponding accession numbers of each strain.

Interspecies transmission and emergence of novel viruses: lessons from bats and birds

Pathogenicity and pathogenesis of a United States porcine deltacoronavirus cell culture isolate in 5-day-old neonatal piglets

Detection of a novel and highly divergent coronavirus from asian leopard cats and Chinese ferret badgers in Southern China

Isolation and characterization of porcine deltacoronavirus from pigs with diarrhea in the United States

Pathogenicity of 2 porcine deltacoronavirus strains in gnotobiotic pigs

Complete genome characterization of Korean porcine Deltacoronavirus strain KOR/KNU14-04/ 2014

Fulllength genome sequence of porcine deltacoronavirus strain USA/IA/2014/8734

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

Complete genome sequence of strain SDCV/USA/Illinois121/2014, a porcine deltacoronavirus from the United States

Rapid detection, complete genome sequencing, and phylogenetic analysis of porcine Deltacoronavirus

Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea viruses associated with outbreaks of severe diarrhea in piglets in Jiangxi

Outbreak of porcine epidemic diarrhea in suckling piglets

2013: MEGA6 molecular evolutionary genetics analysis version 6.0

Detection and genetic characterization of deltacoronavirus in pigs

Porcine coronavirus HKU15 detected in 9 US States

Complete genome sequence of porcine coronavirus HKU15 strain IN2847 from the United States

Comparative analysis of complete genome sequences of three avian Coronaviruses reveals a novel group 3c Coronavirus

Coronavirus genomics and bioinformatics analysis

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

This work was supported by grants from The National Natural Science Foundation of China (31260611, and 31372457), Landing Program Fund of Jiangxi Province (KJLD13029) and Natural Science Foundation of Jiangxi Province, China (20132BAB204017).