key: cord-0708764-1fnvjdyn authors: Feng, Yu; Xu, Zhiwen; Zhu, Ling title: Prevalence and phylogenetic analysis of porcine deltacoronavirus in Sichuan province, China date: 2020-09-06 journal: Arch Virol DOI: 10.1007/s00705-020-04796-z sha: 6bffd7fb2a05088157d446a96ce52e34c39b067f doc_id: 708764 cord_uid: 1fnvjdyn In order to understand the prevalence and genetic diversity of porcine deltacoronavirus (PDCoV) in diarrhoeal pigs in Sichuan province, 634 clinical samples were collected from individual pigs with diarrhoea in 13 regions of Sichuan province, China, from January 2017 and June 2019. The detection results showed that the infection rate of PDCoV was relatively low in diarrhoeal pigs, 13.25% (84/634), but the infection rate of PEDV (porcine epidemic diarrhea virus) was high, 32.18% (204/634). Coinfection with PEDV was common (55.95%, 47/84) in PDCoV-infected diarrhoeal pigs. Additionally, the chance of PDCoV infection was 2.77 times higher in suckling piglets than in sows, and about 3.30 times higher in spring and winter than in summer. PDCoV/PEDV coinfection was 75% less likely in sows than in suckling piglets. The complete genomes of four Sichuan PDCoV strains were sequenced and analysed. There were some insertion-deletion signatures in the whole genome sequences of four strains, including a 6-nt deletion in the non-structural gene 2 region, a 9-nt insertion in the non-structural gene 3 region, a 3-nt deletion in the S gene region, and a distinguishing 11-nt deletion in the 3’UTR region. Phylogenetic analysis based on complete genome sequences revealed that the PDCoV Sichuan strains were closely related to other Chinese PDCoV reference strains; however, phylogenetic analysis based on S gene sequences showed that the CH/SC/2019 strain clustered in a large clade with strains from the USA, Japan, and Korea. These data advance our understanding of the genetic diversity and evolutionary characteristics of PDCoV in China and may contribute to vaccine development. Porcine deltacoronavirus (PDCoV), a member of the genus Deltacoronavirus, is a single-stranded, positive-sense enveloped RNA virus that causes enteritis and watery diarrhoea in pigs, especially nursing piglets [1] . The full genome of PDCoV is approximately 25.4 kb long, which is the smallest genome size among the porcine coronaviruses, arranged in the order 5′ untranslated region (UTR), open reading frame 1a and b (ORF1ab), spike (S), envelope (E), membrane (M), non-structural gene 6 (NS6), nucleocapsid (N), nonstructural gene 7 (NS7), and 3′UTR [2] . The PDCoV S gene is frequently used for studying genetic relationships among different PDCoV strains and for performing epidemiological investigations [3] . PDCoV was first discovered in Hong Kong, and has since been reported in many countries, including the USA, Canada, Korea, China, Thailand, Laos, and Vietnam [4] [5] [6] [7] [8] [9] . In China, PDCoV had been reported in many provinces, but little has been reported regarding the prevalence and epidemiology of PDCoV in Sichuan province, a major pig-raising province of China. Thus, in this study, we investigated the prevalence of PDCoV in diarrhoeal pigs and analysed the complete genome sequences of PDCoV isolates from Sichuan, China. In this study, a total of 634 intestinal and faecal clinical samples (405 intestinal samples and 229 faecal samples) were collected from pigs with diarrhoea outbreaks diagnosed by a pig farm veterinarian from January 2017 and June 2019. All Handling Editor: Tim Skern. * Zhiwen Xu abtcxzw@126.com 1 samples were collected from suckling piglets and sows and from pig farm in 13 regions in Sichuan province. The sampling locations are shown in Table 1 . Soon after sampling, faeces and intestinal contents were prepared as 10% suspensions in phosphate-buffered saline (PBS). The suspensions were then vortexed and centrifuged at 8000 g at 4°C for 10 min. The clarified supernatants were stored at -80°C until use for RNA extraction. Total RNAs were extracted using TRIzol Reagent (Invitrogen, USA) according to the manufacturer's instructions. RT-qPCR was carried out to assess the frequency of PDCoV monoinfections and coinfections with PEDV and transmissible gastroenteritis virus (TGEV). We used the virotype PEDV/TGEV/PDCoV RT-PCR Rgt (QIAGEN) following the manufacturer's protocol. Representative PDCoV strains from four different pig farms where severe or moderate diarrhoea had been reported were selected for sequencing of the complete genome by RT-PCR amplification of 24 regions covering the PDCoV genome as described previously [10] . These representative PDCoV strains were named CH/SC/2017, CH/SC/2018/1, CH/SC/2018, and CH/SC/2019 (GenBank accession nos. MK211169, MK005882, MK330605, and MK993519). The PCR products were purified and cloned into the vector pMD18-T. Positive clones of each amplicon was sequenced twice (Sangon Biotech, Shanghai, China). At least three independent PCR amplicons were sequenced to obtain a consensus sequence. Complete sequences of four PDCoV strains were assembled by combining overlapping contigs and trimming off primer sequences using the DNASTAR software package (DNASTAR Inc., Madison, WI, USA). Sequence alignment analysis of these four strains and other PDCoV strains with sequences in the GenBank database was conducted using the Megalign program. Phylogenetic trees based on complete genome and S gene sequences were constructed by the maximum-likelihood method using MEGA software version 6 with bootstrap analysis of 1,000 replicates. The rates of infection with PDCoV and coinfection with PDCoV and PEDV in diarrhoeal pigs were analysed according to age, year, and season. Significance was set at p < 0.05. Possible risk factors associated with virus infection were analysed by univariate and multivariate exact logistic regression models using Stata/SE 15.1 (Stata Corporation, USA). For each model the odds ratio (OR) and 95% confidence interval (CI) were calculated. Of the 634 clinical samples tested in our PDCoV prevalence survey, 84 (13.25%) were positive for PDCoV and 204 (32.18%) were PEDV positive. The positive rate for PDCoV in sows and piglets was 6.80% (14/206) and 16.36% (70/428), respectively. All of the samples were negative for TGEV, as shown in Table 2 . The results indicate that the prevalence rate of PDCoV was relatively low in diarrhoeal pigs in Sichuan and that PEDV was the main pathogen causing diarrhoea in pigs. It is also noteworthy that 47 samples that were positive for PDCoV were also positive for PEDV, showing that mixed infections with these viruses are relatively common. The PDCoV infection rates in different regions of Sichuan province are shown in Table 1 . Twelve out of the 13 regions were considered PDCoV positive, including the regions of Mianyang (14.93%), Deyang (16.67%), and Guanan (22.73%), suggesting that PDCoV is widespread on pig farms in Sichuan province. As shown in Table 3 , The results show that the prevalence of PDCoV in diarrhoeal pigs increased from 10.53% in 2017 to 17.24% in 2019; however, based on statistical analysis, we could not conclude that PDCoV infection in diarrhoeal pigs has been increasing from year to year (p > 0.05). More clinical samples are needed to determine whether the PDCoV infection rate is correlated with the The complete genomes of four novel PDCoV strains, representative of epidemic strains, were sequenced and analysed to perform a molecular characterisation of PDCoV (Fig. 1) . All four sequenced strains had a 6-nt deletion in the non-structural gene (ns) 2 region (ORF1ab), which has been observed in some other Chinese strains. Additionally, CH/SC/2019 had a 9-nt deletion in the ns 3 region (ORF1ab) when compared with most Chinese strains. This deletion is also present in strain CH/Sichuan/S27/2012, detected previously in Sichuan, and in strains from Thailand Phylogenetic analysis was carried out using the whole-genome and S gene sequences of the four sequenced strains and other reference PDCoV strains ( Fig. 2A and B . The phylogenies based on the complete genome sequences and the S gene both indicates that the nucleotide exactly matches the reference sequence. A dash (-) indicates that the nucleotide is deleted relative to the reference sequence revealed that the PDCoV strains showed obvious regional characteristics, as has been reported previously [11] . Phylogenetic analysis based on whole genome sequences indicated that the PDCoV CH/SC/2017, CH/SC/ 2018/1, and CH/SC/ 2018 strains grouped within the Chinese clade and are more closely related to strains SHJS/SL/2016 and CH/HN/2014. However, the CH/SC/2019 strain is located in another subcluster and is similar to the CH/QH/2017 and CH/GS/2017 strains, which were isolated from piglets in recent years, suggesting that the origin of these isolates may be different. A phylogenetic tree based on S gene sequences revealed that CH/SC/2017, CH/SC/2018/1, and CH/SC/2018 belong to the same group, while the CH/SC/2019 strain clustered in a large clade with strains from the USA, Japan, and Korea and demonstrated genetic differences when compared to the other Sichuan strains. In this study, RT-PCR analysis of 634 clinical samples from pigs with diarrhoea collected from 13 regions in Sichuan province from January 2017 and June 2019 revealed a positive rate of PDCoV of 13.25%, while that of PEDV was 32.18%, indicating that PEDV is still the primary pathogen causing porcine diarrhoea, which is consistent with previous studies [12] . The prevalence of PDCoV in diarrhoeal pigs in Sichuan was somewhat lower than in other provinces, such as Henan (23.49%) and Guangdong (21.8%), but a molecular-based detection study showed a similar low prevalence of PDCoV infection in Jiangsu (4.94%) and Hubei (6.45%) [13] . The prevalence of PDCoV in Qinghai and Tibet is also low, with a [14] . These data show that the prevalence of PDCoV in diarrhoeal pigs differs among the Chinese provinces, and more information regarding the molecular epidemiology of PDCoV in China is needed. In this study, the number of samples collected and regions investigated was limited, and therefore more samples need to be collected and tested, and PDCoV surveillance should be continued. Also, since the aim of this study was only to monitor the prevalence of PDCoV in diarrhoeal pigs more samples, including those from pigs without diarrhoea, are needed to evaluate the prevalence of PDCoV in Sichuan pig herds. Although the overall prevalence of PDCoV in this study was relatively low, PDCoV was detected in most regions of Sichuan, implying that PDCoV prevention measures have not been very effective, and further study is needed to prevent and control this novel enteric virus. Statistical analysis demonstrated a higher chance of PDCoV infection in winter and spring, which is consistent with previous studies [15] , suggesting that we should pay more attention to preventing infections during these seasons. Furthermore, we detected PDCoV in sows, as described previously [16] . However, unlike a previous study, our results showed that suckling piglets had a 2.77 times higher risk of PDCoV infection than sows. Previous reports have suggested that older pigs are more susceptible to PDCoV than piglets, which was not confirmed by our results [17] . Therefore, more clinical samples from older pigs are needed to test this conclusion. We plan to continue to monitor PDCoV infections in sows. Moreover, it has been reported that PDCoV/PEDV coinfection is more common in Chinese pigs than PDCoV/porcine rotavirus coinfection in US pig herds [18, 19] . Our detection data are consistent with this, also showing that suckling piglets have a higher chance of PDCoV/PEDV coinfection. The results of sequence analysis showed that the four strains had some unique deletions/insertions compared to other PDCoV strains. The CH/SC/2019 strain, like the CH/ Sichuan/S27/2012 strain, had a 6-nt deletion in the ns 2 region and a 9-nt deletion in the ns 3 region, respectively, while the others strains only had a 6-nt deletion in ns 2. Moreover, four stains had a 3-nt deletion in the deletion of the S gene, resulting in one amino acid. This deletion exists in most Chinese PDCoV isolates. The S1 region is responsible for virus receptor binding [20] , suggesting that this aa deletion may have some biological significance regarding viral replication or pathogenesis. Importantly, the CH/SC/2017 and CH/SC/2019 strains had a unique 11-nt deletion in the 3′UTR region. This is the first report of a PDCoV strain with such a deletion. However, we do not know whether these naturally occurring deletions or insertions have biological significance in PDCoV biology and pathogenesis. The four sequenced PDCoV strains were all collected from suckling piglets, but the clinical signs in these piglets were different. The piglets infected with the CH/SC/2017 and CH/SC/2019 strains suffered from severe watery diarrhoea and died after three days of diarrhoea, while the suckling piglet with the CH/SC/2018 and CH/ SC/2018/1 strains had moderate diarrhoea and then recovered. Were these different clinical signs associated with the 11-nt deletion in the 3′UTR or another deletion/insertion? We will continue to isolate PDCoV strains for further investigation. A phylogenetic tree constructed based on complete genome sequences indicated that these strains might have evolved from different Chinese strains. In the phylogenetic tree based on the S gene, the Sichuan strains clustered in a different clade. The CH/SC/2017, CH/SC/2018/1, and CH/ SC/2018 strains are in the China lineage, like most Chinese PDCoVs, while the CH/SC/2019 strain is in the USA/Japan/ South Korea lineage. This is possibly due to recombination or separate evolution Additional PDCoV sequences are needed to identify possible recombination and evolution events. In summary, we report for the first time the prevalence and genetic properties of PDCoV in diarrhoeal pigs in Sichuan, which could provide further insights into the epidemiology and evolution of PDCoV in China. Pathogenicity of 2 porcine deltacoronavirus strains in gnotobiotic pigs Isolation and characterisation of porcine deltacoronavirus from pigs with diarrhea in the United States Cryo-EM structure of porcine delta coronavirus spike protein in the pre-fusion state 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 Rapid detection, complete genome sequencing, and phylogenetic analysis of porcine deltacoronavirus Herd-level prevalence and incidence of porcine epidemic diarrhoea virus (PEDV) and porcine deltacoronavirus (PDCoV) in swine herds in Ontario Complete genome characterization of Korean porcine deltacoronavirus strain KOR/KNU14-04 Isolation and phylogenetic analysis of porcine deltacoronavirus from pigs with diarrhoea in Hebei province Retrospective investigation and evolutionary analysis of a novel porcine deltacoronavirus strain detected in Thailand from Detection and genetic characterization of deltacoronavirus in pigs Porcine Deltacoronavirus: Overview of infection dynamics, diagnostic methods, prevalence and genetic evolution Prevalence, complete genome sequencing and phylogenetic analysis of porcine deltacoronavirus in South Korea Porcine deltacoronavirus in mainland China Detection and genetic characterization of porcine deltacoronavirus in Tibetan pigs surrounding the Qinghai-Tibet Plateau of China Prevalence, phylogenetic and evolutionary analysis of porcine deltacoronavirus in Henan province The detection and phylogenetic analysis of porcine deltacoronavirus from Guangdong Province in Southern China Genetic characterization and pathogenicity of Japanese porcine deltacoronavirus Newly emerged porcine deltacoronavirus associated with diarrhoea in swine in China: identification, prevalence and full-length genome sequence analysis Rapid detection, complete genome sequencing, and phylogenetic analysis of porcine Deltacoronavirus Broad receptor engagement of an emerging global coronavirus may potentiate its diverse cross-species transmissibility Conflict of interest The authors declare that they have no conflict of interest.Ethical approval This project was approved by the appropriate ethical review committee. The "Guidelines for Experimental Animals" of the Ministry of Science and Technology (Beijing, China) were followed.