Next Article in Journal
Investigating the Current Status of SARS-CoV-2 Antibodies in Hospital Staff
Previous Article in Journal
Divergent Hantavirus in Somali Shrews (Crocidura somalica) in the Semi-Arid North Rift, Kenya
Previous Article in Special Issue
Magnolol, a Neolignan-like Drug, Inhibits Porcine Epidemic Diarrhea Virus Replication in Cultured Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 12
Issue 5
10.3390/pathogens12050686
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Genetic Characterization and Evolution of Porcine Deltacoronavirus Isolated in the Republic of Korea in 2022

by
Hye-Ryung Kim
1,2,
Jonghyun Park
1,2,
Kyoung-Ki Lee
3,
Hye-Young Jeoung
3,
Young S. Lyoo
4,
Seung-Chun Park
1,5 and
Choi-Kyu Park
1,*
1
College of Veterinary Medicine & Animal Disease Intervention Center, Kyungpook National University, Daegu 41566, Republic of Korea
2
DIVA Bio Incorporation, Daegu 41519, Republic of Korea
3
Animal and Plant Quarantine Agency, Gimcheon 39660, Republic of Korea
4
College of Veterinary Medicine, Konkuk University, Seoul 05029, Republic of Korea
5
Laboratory of Veterinary Pharmacokinetics and Pharmacodynamics, College of Veterinary Medicine and Cardiovascular Research Institute, Kyungpook National University, Daegu 41566, Republic of Korea
*
Author to whom correspondence should be addressed.
Pathogens 2023, 12(5), 686; https://doi.org/10.3390/pathogens12050686
Submission received: 31 March 2023 / Revised: 4 May 2023 / Accepted: 5 May 2023 / Published: 7 May 2023
(This article belongs to the Special Issue Porcine Coronaviruses Infection)
Browse Figures
Versions Notes

Abstract

:
Porcine deltacoronavirus (PDCoV) is an emerging coronavirus that causes diarrhea in nursing piglets. Since its first outbreak in the United States in 2014, this novel porcine coronavirus has been detected worldwide, including in Korea. However, no PDCoV case has been reported since the last report in 2016 in Korea. In June 2022, the Korean PDCoV strain KPDCoV-2201 was detected on a farm where sows and piglets had black tarry and watery diarrhea, respectively. We isolated the KPDCoV-2201 strain from the intestinal samples of piglets and sequenced the viral genome. Genetically, the full-length genome and spike gene of KPDCoV-2201 shared 96.9–99.2% and 95.8–98.8% nucleotide identity with other global PDCoV strains, respectively. Phylogenetic analysis suggested that KPDCoV-2201 belongs to G1b. Notably, the molecular evolutionary analysis indicated that KPDCoV-2201 evolved from a clade different from that of previously reported Korean PDCoV strains and is closely related to the emergent Peruvian and Taiwanese PDCoV strains. Furthermore, KPDCoV-2201 had one unique and two Taiwanese strain-like amino acid substitutions in the receptor-binding domain of the S1 region. Our findings suggest the possibility of transboundary transmission of the virus and expand our knowledge about the genetic diversity and evolution of PDCoV in Korea.

1. Introduction

Porcine deltacoronavirus (PDCoV) is an emerging enteropathogenic coronavirus that causes typical clinical symptoms, including severe diarrhea, vomiting, dehydration, and mortality in seronegative neonatal piglets [1,2,3,4,5]. PDCoV is an enveloped, single-stranded, positive-sense RNA virus in the genus Deltacoronavirus, family Coronaviridae, of the order Nidovirales, with a genome of approximately 25.4 kb in length [4,6]. The genome arrangement of PDCoV is as follows: 5′-untranslated region (UTR), open reading frame 1a/1b (ORF1a/1b), spike (S), envelope (E), membrane (M), nonstructural protein 6 (NS6), nucleocapsid (N), NS7, NS7a, and 3′-UTR [6]. PDCoV was first detected in 2012 in pig rectal swabs collected from Hong Kong as a part of molecular surveillance studies to investigate the diversity of coronaviruses in mammals and birds [6]. In early 2014, an outbreak of PDCoV-associated diarrhea broke out in pig farms in Ohio, United States, which led to the identification of the first clinical case of PDCoV infection in pigs [2,4,7,8]. Since then, PDCoV has been detected worldwide, including in Korea, Japan, China, Thailand, Vietnam, Laos, Peru, and Mexico [3,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]. According to previous phylogenetic analyses, PDCoV can be divided into three or four groups based on their full-length genome sequences [15,24]. In Korea, PDCoV (KNU14-04 strain; GenBank accession no. KM820765) was first identified in April 2014 from the fecal samples of pigs with diarrhea [17]. Subsequent studies have shown that all Korean PDCoV strains identified during 2014–2016 have a high sequence identity with strains reported in the United States [10,18,21]. However, since 2016, there have been no reports of PDCoV outbreaks in Korea. In June 2022, we detected PDCoV RNAs from diarrheal pigs in a Korean domestic pig farm, successfully isolated the virus (designated KPDCoV-2201), and analyzed its full-length genome sequences to elucidate the genetic characteristics and evolution of KPDCoV-2201 currently detected in Korea.

2. Materials and Methods

2.1. Sample Collection and RNA Extraction

In June 2022, eight fecal and one intestinal tissue samples were submitted to the Animal Disease Intervention Center of Kyungpook National University from a pig farm located in Gyeongsangbuk-do. Each sample was separately suspended in phosphate-buffered saline, homogenized, and centrifuged. The supernatant of each sample was used to extract total DNA/RNA using Patho Gene-spin™ DNA/RNA Extraction Kit (iNtRON Biotechnology, Sungnam, Korea).

2.2. Detection of Enteric Pathogens

To detect swine enteric pathogens, real-time reverse transcription polymerase chain reaction (qRT–PCR) was performed for 10 viral pathogens [PDCoV, porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis virus (TGEV), porcine enteric alphacoronavirus (PEAV), swine acute diarrhea syndrome coronavirus (SADS-CoV), porcine torovirus (PToV), porcine sapelovirus (PSV), group A rotavirus (RVA), group B rotavirus (RVB), group C rotavirus (RVC)] using a commercial one-step qRT–PCR kit (THUNDERBIRD™ Probe One-step qRT–PCR kit, TOYOBO, Osaka, Japan). The previously reported 10 qRT-PCR assays were performed using the target gene-specific primer and probe sets of each viral pathogen, as described in their studies [7,25,26,27,28,29,30]. The reaction composition and conditions were performed according to the manufacturer’s instructions.

2.3. Isolation and Genome Sequencing of PDCoV

Sequence analysis of the viral genome was performed on seven samples with a cycle threshold (Ct) of ≤25 among the nine PDCoV-positive samples. The spike (S) protein genes were sequenced by RT-PCR using PDCoV-specific primer sets (Table S1). We attempted to isolate PDCoV from the PDCoV-positive intestinal samples using swine testis (ST) cells (ATCC CRL-1746), as described previously [31]. We designed 11 pairs of primers for full-length genome sequencing of KPDCoV-2201 based on the conserved regions of the 127 complete genome sequences available in GenBank (Table S1). The synthesis of each cDNA fragment, PCR amplification, purification of the amplified PCR products, and determination of the sequences of the 5′- and 3′-termini of the viral genome were performed as previously described by Park et al. 2022 [32]. The sequences of each product were analyzed using the Sanger’s method by a commercial company (BIONICS, Daejeon, Republic of Korea). One complete genome sequence (strain KPDCoV-2201) was deposited to GenBank (accession numbers OQ718500).

2.4. Phylogenetic and Molecular Evolutionary Analyses

To determine the genetic relationship between KPDCoV-2201 and other strains, nucleotide (nt) and amino acid (aa) sequences were aligned using MAFFT [33], which is available in Geneious Prime (https://www.geneious.com, accessed on 21 July 2022), together with 127 other full-length PDCoV genome sequences and 273 complete spike gene sequences available in GenBank. Phylogenetic trees were constructed using the RAxML method available in Geneious Prime (https://www.geneious.com, accessed on 21 July 2022) based on the general time-reversible (GTR) nucleotide substitution with a gamma distribution model [34]. These trees were subjected to bootstrap analysis with 1000 replicates to determine the percentage reliability values for each internal node of the tree [35]. All phylogenetic trees were visualized using the iTOL phylogenetic tree viewer [36]. To estimate the temporal distributions of KPDCoV-2201 and the phylogenetic relationships in South Korea, we constructed a time-scaled phylogenetic tree of G1b PDCoV strains using BEAST v.1.10.4 [37] with a GTR + gamma (GTR + γ) substitution model. Time-scaled maximum clade credibility tree was constructed using TreeAnnotator v.1.10.4 (https://beast.community/treeannotator, accessed on 22 July 2022) in BEAST and was visualized using FigTree v.1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/, accessed on 22 September 2022).

3. Results

3.1. Detection and Isolation of KPDCoV-2201

In June 2022, pigs of all ages on a farm located in Gyeongsangbuk-do exhibited the clinical symptoms of watery diarrhea. Approximately 10% of the sows presented with diarrhea symptoms, whereas 90% of the piglets showed diarrhea symptoms with a mortality rate of 40%, lower than the typical mortality rate of PEDV infection (Figure 1). Etiological diagnostic results showed that all nine samples collected from the farm were positive for PDCoV but negative for the remaining nine viral antigens (Table 1). Subsequently, seven PDCoV-positive samples with high viral titers (Ct value of <25) were selected for the sequencing of S gene. The results of the sequencing analysis showed that all S genes from all seven samples shared 100% homology at the nucleotide level. Further, we selected the intestinal sample with the highest viral load and inoculated it into ST cells for virus isolation. After one passage of the cell culture, we observed a remarkable cytopathic effect in the infected cells, indicating that the PDCoV Korean isolate was successfully isolated in this study (Figure S1).

3.2. Full-Length Genomic Characterization of KPDCoV-2201

The full-length genome sequence of KPDCoV-2201 was analyzed to determine the genetic characteristics of PDCoV, which was recently detected in Korea. The genome of KPDCoV-2201 is 25,422 nt in length, excluding the 3′-poly-(A) tail. It is arranged as follows: 5′-UTR comprising 539 nt, ORF1a and ORF1b encoding a replicase of 18,803 nt, S gene comprising 3483 nt, E gene comprising 252 nt, M gene comprising 654 nt, NS6 gene comprising 285 nt, N gene comprising 1326 nt, NS7 gene comprising 603 nt, and 3′-UTR comprising 392 nt. The full-length genome of KPDCoV-2201 was submitted to the NCBI GenBank database (GenBank accession no. OQ718500). A comparison of the nucleotide identity of the full-length genome of KPDCoV-2201 showed 96.9–99.2% identity with global PDCoV strains, with 98.7–99.2% identity with Korean PDCoV strains. When the genome of Korean PDCoV strains and that of KPDCoV-2201 were compared, it was found that the M gene was the most conserved, with 99.7–99.8% nucleotide homology. In contrast, the NS7 gene was the most mutated, with 97.5–98.8% nucleotide homology. Interestingly, the E and M genes of KPDCoV-2201 were more similar to those of the Peruvian PDCoV strain (GenBank accession no. MT227371) than to those of other Korean PDCoV strains (Table 2). According to a previously reported structural analysis, the S protein PDCoV contains the S1-NTD (N-terminal domain of S1), S1-CTD (C-terminal domain of S1), CH-N (N-terminal central helix), CH-C (C-terminal central helix), FP (fusion peptide), HR-N (N-terminal heptad repeat), and HR-C (C-terminal heptad repeat), and the PDCoV receptor binding domain (RBD) is located at the C-terminal of the S1 subunit (S1-CTD) of the S protein [38]. Based on the S gene of KPDCoV-2201, the aa sequences of two PDCoV strains, which are closely related to KPDCoV-2201 according to phylogenetic analysis, were aligned and compared with those of the previously reported eight Korean PDCoV strains (Figure S2A). The alignment results revealed that KPDCoV-2201 has 18–31 aa substitutions with Korean PDCoV strains, of which 16 aa substitutions (K96R, V277L, V326I, T387S, R342K, E/K485G, K534N, I550L, A624T, A630T, L632M, N639S, R/Q642K, A706T, K815T, and E854D) were not identical to those of all Korean PDCoV sequences, and 2 aa substitutions (I190V and D284N) were identical to those of only one strain each. KPDCoV-2201 had three aa differences from the previously reported eight Korean PDCoV strains in the S1-CTD region, but two of three were identical to those of the recently emerged PDCoV in the Taiwan strain (GenBank accession no. MZ712038) (Figure S2B).

3.3. Phylogenetic and Molecular Evolutionary Analyses

For phylogenetic analysis, a RAxML tree was constructed using the genomic sequence of the KPDCoV-2201 strain along with previously reported 127 complete genomic sequences or 273 spike gene sequences, respectively (Figure 2). According to a previous report [15], PDCoV strains were classified into two genogroups (G)s, G1 and G2. G1 is divided into G1a (China) and G1b (United States). G2 includes PDCoV originating from Southeast Asia. Based on phylogenetic analysis of the full-length genome, the PDCoV-2201 strain was classified as G1b clustered with the US-like PDCoVs from the United States, Japan, Taiwan, Peru, Haiti, and the previously reported Korean strains (Figure 2A). Moreover, in the phylogenetic tree constructed using the S gene, KPDCoV-2201 was more closely related to PDCoV strains from Taiwan and Peru and was not clustered with other Korean PDCoV strains (Figure 2B).
To elucidate the evolutionary characteristics of Korean PDCoV strains, times to the most recent common ancestor (tMRCA) were analyzed using the S gene of strains belonging to the G1b group. The tMRCA of G1b was estimated as 2010.856, with a 95% highest posterior density (HPD) ranging from 2008.49 to 2012.77, which diverged to clades 1 and 2. These clades individually had tMRCA between 2011 and 2012. The Korean PDCoV strains that caused outbreaks from 2014 to 2016 belonged to clade 1, which also included most of the American, Japanese, and Mexican strains, and one Haitian strain. The tMRCA of this clade was 2012.29, with a 95% HPD interval ranging from 2011.18 to 2013.31. However, KPDCoV-2201 belonged to clade 2, which included the American, Japanese, and Mexican strains and the Peruvian and Taiwanese strains, which caused the recent outbreaks, with a tMRCA of 2011.78 (95% HPD, 2010.88–2012.57) (Figure 3).

4. Discussion

Recent studies have reported recombinant PDCoV strains in Southeast Asia and highly pathogenic strains in China, suggesting that PDCoV exhibits continuous antigenic changes [14,15,39]. Moreover, a case of PDCoV infection in humans was reported in Haiti [40], highlighting the risk of cross-species transmission and zoonosis. Therefore, it is necessary to conduct further research to expand insights into the PDCoV infection status and its genetic epidemiology for the global distribution of the virus.
To our knowledge, there has been no study on PDCoV in Korea for approximately the past 7 years. In 2022, KPDCoV-2201, a Korean PDCoV strain, was isolated from ST cells, and its full-length genome was sequenced in this study. Phylogenetic analysis based on the 128 full-length genome sequences showed that global PDCoV strains were classified into three distinct genogroups: G1a (China group), G1b (United States group), and G2 (Southeast Asia group) as previously reported [14,15] (Figure 2). All Korean PDCoV strains, including KPDCoV-2201 isolated in this study, were classified into the G1b genogroup, but the KPDCoV-2201 was more closely related to Taiwan PDCoV strains than to previously reported Korean PDCoV strains.
Due to research gaps, it remains unclear whether the origin of KPDCoV-2201 is the result of evolutions within Korean PDCoV strains or whether it was recently introduced as part of the Asian epidemic. Therefore, to further elucidate the evolutionary characteristics of the KPDCoV-2201 strain, tMRCA analysis was conducted using the S gene sequences of strains belonging to the G1b group. As a result, the KPDCoV-2201 was classified as Clade 2 and shared the same tMRCA with recently reported Taiwanese and Peruvian strains, which differed from previously reported Korean PDCoV strains belonging to Clade 1 (Figure 3). It is noteworthy that the occurrences of Clade 2 PDCoV strains increased compared to their minor occurrence in 2014-2016. The increasing occurrence of Clade 2 strains may indeed complicate PDCoV epidemiology and make the prevention and control of the disease more challenging.
A previous study reported that the titers of the recently emerged Taiwanese PDCoV strains were higher than those of the strain detected in 2015, indicating severe symptoms due to single infections with these PDCoV strains [41]. In our study, we confirmed that severe clinical symptoms were observed as a consequence of a single infection with PDCoV without coinfection by the other nine types of enteric viral pathogens in sows and piglets. These results suggested that PDCoV could potentially be the primary causative agent of severe diarrhea in some pig farms. However, diagnosing for PDCoV infection may not have been routinely conducted in diagnostic laboratories since its clinical symptoms are similar to those of other major viral pathogens causing severe enteritis, such as PEDV and TGEV. Therefore, it is essential to diagnose PDCoV in domestic pig populations exhibiting the symptoms of diarrhea, and further research is warranted on the pathogenicity of the current circulating PDCoV.

5. Conclusions

We isolated the KPDCoV-2201 strain from a Korean domestic pig farm where sows and piglets had severe watery diarrhea. Genetic and evolutionary analysis of KPDCoV-2201 revealed that the newly identified Korean strain is distinct from previously reported Korean strains but is closely related to non-Korean strains derived from Taiwan and Peru. Our findings suggest the possibility of transboundary transmission of the virus and expand our knowledge about the genetic diversity and evolution of PDCoV in South Korea.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens12050686/s1, Table S1: Primer sets for the complete genome sequencing of PDCoV in this study; Figure S1: KPDCoV-2201 strain causes cytopathic changes in ST cells; Figure S2: Alignment of the amino acid sequences of spike protein in KPDCoV-2201 with Korean PDCoV strains (A). Details of amino acid substitutions of KPDCoV-2201 that differ from KNU14-04 (B).

Author Contributions

Conceptualization, C.-K.P.; methodology, H.-R.K., J.P., C.-K.P.; software, H.-R.K., J.P.; validation, H.-R.K.; formal analysis, H.-R.K., J.P.; investigation, H.-R.K.; resources, K.-K.L., H.-Y.J.; data curation, H.-R.K.; writing—original draft preparation, H.-R.K.; writing—review and editing, Y.S.L., S.-C.P., C.-K.P., K.-K.L.; visualization, H.-R.K., J.P.; supervision, C.-K.P.; project administration, C.-K.P.; funding acquisition, K.-K.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Animal and plant quarantine agency [grant number B-1543069-2023-25-01].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jung, K.; Hu, H.; Eyerly, B.; Lu, Z.; Chepngeno, J.; Saif, L.J. Pathogenicity of 2 Porcine Deltacoronavirus Strains in Gnotobiotic Pigs. Emerg. Infect. Dis. 2015, 21, 650–654. [Google Scholar] [CrossRef]
  2. Li, G.; Chen, Q.; Harmon, K.M.; Yoon, K.J.; Schwartz, K.J.; Hoogland, M.J.; Gauger, P.C.; Main, R.G.; Zhang, J. Full-length genome sequence of porcine deltacoronavirus strain USA/IA/2014/8734. Genome Announc. 2014, 2, e00278-14. [Google Scholar] [CrossRef]
  3. Song, D.; Zhou, X.; Peng, Q.; Chen, Y.; Zhang, F.; Huang, T.; Zhang, T.; Li, A.; Huang, D.; Wu, Q.; et al. Newly emerged porcine deltacoronavirus associated with diarrhoea in swine in China: Identification, prevalence and full-length genome sequence analysis. Transbound. Emerg. Dis. 2015, 62, 575–580. [Google Scholar] [CrossRef]
  4. Wang, L.; Byrum, B.; Zhang, Y. Detection and Genetic Characterization of Deltacoronavirus in Pigs, Ohio, USA, 2014. Emerg. Infect. Dis. 2014, 20, 1227–1230. [Google Scholar] [CrossRef]
  5. Wang, L.; Hayes, J.; Sarver, C.; Byrum, B.; Zhang, Y. Porcine deltacoronavirus: Histological lesions and genetic characterization. Arch. Virol. 2015, 161, 171–175. [Google Scholar] [CrossRef]
  6. Woo, P.C.; Lau, S.K.; Lam, C.S.; Lau, C.C.; Tsang, A.K.; Lau, J.H.; Bai, R.; Teng, J.L.; Tsang, C.C.; Wang, M.; et al. 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. J. Virol. 2012, 86, 3995–4008. [Google Scholar] [PubMed]
  7. Marthaler, D.; Raymond, L.; Jiang, Y.; Collins, J.; Rossow, K.; Rovira, A. Rapid Detection, Complete Genome Sequencing, and Phylogenetic Analysis of Porcine Deltacoronavirus. Emerg. Infect. Dis. 2014, 20, 1347–1350. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, L.; Zhang, Y.; Byrum, B. Complete Genome Sequence of Porcine Coronavirus HKU15 Strain IN2847 from the United States. Genome Announc. 2014, 2, e00291-14. [Google Scholar] [CrossRef]
  9. Chen, F.; Zhu, Y.; Wu, M.; Ku, X.; Yao, L.; He, Q. Full-length genome characterization of Chinese porcine deltacoronavirus strain CH/SXD1/2015. Genome Announc. 2015, 3, e01284-15. [Google Scholar] [CrossRef] [PubMed]
  10. Chung, H.C.; Nguyen, V.G.; Oh, W.T.; My Le, H.T.; Moon, H.J.; Lee, J.H.; Kim, H.K.; Park, S.J.; Park, B.K. Complete genome sequences of porcine deltacoronavirus strains DH1/2016 and DH2/2016 isolated in South Korea. Genome Announc. 2017, 5, e01706-16. [Google Scholar] [CrossRef]
  11. Dong, N.; Fang, L.; Zeng, S.; Sun, Q.; Chen, H.; Xiao, S. Porcine Deltacoronavirus in Mainland China. Emerg. Infect. Dis. 2015, 21, 2254–2255. [Google Scholar] [CrossRef]
  12. Saeng-Chuto, K.; Stott, C.J.; Wegner, M.; Senasuthum, R.; Tantituvanont, A.; Nilubol, D. Retrospective investigation and evolutionary analysis of a novel porcine deltacoronavirus strain detected in Thailand from 2008 to 2015. Arch. Virol. 2017, 162, 2103–2108. [Google Scholar] [CrossRef]
  13. Saeng-Chuto, K.; Lorsirigool, A.; Temeeyasen, G.; Vui, D.T.; Stott, C.J.; Madapong, A.; Tripipat, T.; Wegner, M.; Intrakamhaeng, M.; Chongcharoen, W.; et al. Different Lineage of Porcine Deltacoronavirus in Thailand, Vietnam and Lao PDR in 2015. Transbound. Emerg. Dis. 2016, 64, 3–10. [Google Scholar] [CrossRef]
  14. Saeng-Chuto, K.; Jermsutjarit, P.; Stott, C.J.; Vui, D.T.; Tantituvanont, A.; Nilubol, D. Retrospective study, full-length genome characterization and evaluation of viral infectivity and pathogenicity of chimeric porcine deltacoronavirus detected in Vietnam. Transbound. Emerg. Dis. 2020, 67, 183–198. [Google Scholar] [CrossRef]
  15. Stott, C.J.; Sawattrakool, K.; Saeng-Chuto, K.; Tantituvanont, A.; Nilubol, D. The phylodynamics of emerging porcine deltacoronavirus in Southeast Asia. Transbound. Emerg. Dis. 2022, 69, 2816–2827. [Google Scholar] [CrossRef]
  16. Suzuki, T.; Shibahara, T.; Imai, N.; Yamamoto, T.; Ohashi, S. Genetic characterization and pathogenicity of Japanese porcine deltacoronavirus. Infect. Genet. Evol. 2018, 61, 176–182. [Google Scholar] [CrossRef] [PubMed]
  17. Lee, S.; Lee, C. Complete genome characterization of Korean porcine deltacoronavirus strain KOR/KNU14-04/2014. Genome Announc. 2014, 2, e01191-14. [Google Scholar] [CrossRef] [PubMed]
  18. Lee, J.H.; Chung, H.C.; Nguyen, V.G.; Moon, H.J.; Kim, H.K.; Park, S.J.; Lee, C.H.; Lee, G.E.; Park, B.K. Detection and Phylogenetic Analysis of Porcine Deltacoronavirus in Korean Swine Farms, 2015. Transbound. Emerg. Dis. 2016, 63, 248–252. [Google Scholar] [CrossRef]
  19. Lorsirigool, A.; Saeng-Chuto, K.; Madapong, A.; Temeeyasen, G.; Tripipat, T.; Kaewprommal, P.; Tantituvanont, A.; Piriyapongsa, J.; Nilubol, D. The genetic diversity and complete genome analysis of two novel porcine deltacoronavirus isolates in Thailand in 2015. Virus Genes 2017, 53, 240–248. [Google Scholar] [CrossRef] [PubMed]
  20. Janetanakit, T.; Lumyai, M.; Bunpapong, N.; Boonyapisitsopa, S.; Chaiyawong, S.; Nonthabenjawan, N.; Kesdaengsakonwut, S.; Amonsin, A. Porcine deltacoronavirus, Thailand, 2015. Emerg. Infect. Dis. 2016, 22, 757–759. [Google Scholar] [CrossRef]
  21. Jang, G.; Lee, K.-K.; Kim, S.-H.; Lee, C. Prevalence, complete genome sequencing and phylogenetic analysis of porcine deltacoronavirus in South Korea, 2014–2016. Transbound. Emerg. Dis. 2017, 64, 1364–1370. [Google Scholar] [CrossRef]
  22. Pérez-Rivera, C.; Ramírez-Mendoza, H.; Mendoza-Elvira, S.; Segura-Velázquez, R.; Sánchez-Betancourt, J.I. First report and phylogenetic analysis of porcine deltacoronavirus in Mexico. Transbound. Emerg. Dis. 2019, 66, 1436–1441. [Google Scholar] [CrossRef] [PubMed]
  23. Vicente-Huaman, J.; Gómez-Quispe, O.E. Evaluation of a Porcine deltacoronavirus eradication program in a full-cycle pig farm in Peru. J. Adv. Vet. Anim. Res. 2021, 8, 300–306. [Google Scholar] [CrossRef] [PubMed]
  24. He, W.-T.; Ji, X.; He, W.; Dellicour, S.; Wang, S.; Li, G.; Zhang, L.; Gilbert, M.; Zhu, H.; Xing, G.; et al. Genomic Epidemiology, Evolution, and Transmission Dynamics of Porcine Deltacoronavirus. Mol. Biol. Evol. 2020, 37, 2641–2654. [Google Scholar] [CrossRef]
  25. Chun, Y.-H.; Jeong, Y.-J.; Park, S.-I.; Hosmillo, M.; Shin, D.-J.; Kwon, H.-J.; Kang, S.-Y.; Woo, S.-K.; Kang, M.-I.; Park, S.-J.; et al. Development of One-Step Real-Time Reverse Transcription Polymerase Chain Reaction Assays for Rapid Detection of Porcine Group C Rotaviruses. J. Vet. Diagn. Investig. 2010, 22, 74–77. [Google Scholar] [CrossRef]
  26. Huang, X.; Chen, J.; Yao, G.; Guo, Q.; Wang, J.; Liu, G. A TaqMan-probe-based multiplex real-time RT-qPCR for simultaneous detection of porcine enteric coronaviruses. Appl. Microbiol. Biotechnol. 2019, 103, 4943–4952. [Google Scholar] [CrossRef]
  27. Kim, J.K.; Kim, H.R.; Kim, D.Y.; Kim, J.M.; Kwon, N.Y.; Park, J.H.; Park, J.Y.; Kim, S.H.; Lee, K.K.; Lee, C.; et al. A simple colorimetric detection of porcine epidemic diarrhea virus by reverse transcription loop-mediated isothermal amplification assay using hydroxynaphthol blue metal indicator. J. Virol. Methods 2021, 298, 114289. [Google Scholar] [CrossRef] [PubMed]
  28. Pan, Z.; Lu, J.; Wang, N.; He, W.T.; Zhang, L.; Zhao, W.; Su, S. Development of a TaqMan-probe-based multiplex real-time PCR for the simultaneous detection of emerging and reemerging swine coronaviruses. Virulence 2020, 11, 707–718. [Google Scholar] [CrossRef]
  29. Kumari, S.; Ray, P.K.; Singh, R.; Desingu, P.A.; Sharma, G.T.; Saikumar, G. Development of a Taqman-based real-time PCR assay for detection of porcine sapelovirus infection in pigs. Anim. Biotechnol. 2020, 31, 264–267. [Google Scholar] [CrossRef]
  30. Masuda, T.; Nagai, M.; Yamasato, H.; Tsuchiaka, S.; Okazaki, S.; Katayama, Y.; Oba, M.; Nishiura, N.; Sassa, Y.; Omatsu, T.; et al. Identification of novel bovine group A rotavirus G15P[14] strain from epizootic diarrhea of adult cows by de novo sequencing using a next-generation sequencer. Vet. Microbiol. 2014, 171, 66–73. [Google Scholar] [CrossRef]
  31. Hu, H.; Jung, K.; Vlasova, A.N.; Chepngeno, J.; Lu, Z.; Wang, Q.; Saif, L.J. Isolation and Characterization of Porcine Deltacoronavirus from Pigs with Diarrhea in the United States. J. Clin. Microbiol. 2015, 53, 1537–1548. [Google Scholar] [CrossRef] [PubMed]
  32. Park, J.; Kim, H.; Kim, J.; Lee, K.; Kim, W.; Lyoo, Y.S.; Kwon, O.; Park, C.; Park, S. First report of Porcine respirovirus 1 in South Korea. Transbound. Emerg. Dis. 2022, 69, 4041–4047. [Google Scholar] [CrossRef] [PubMed]
  33. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef]
  34. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  35. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  36. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  37. Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 2018, 4, vey016. [Google Scholar] [CrossRef]
  38. Ji, W.; Peng, Q.; Fang, X.; Li, Z.; Li, Y.; Xu, C.; Zhao, S.; Li, J.; Chen, R.; Mo, G.; et al. Structures of a deltacoronavirus spike protein bound to porcine and human receptors. Nat. Commun. 2022, 13, 1467. [Google Scholar] [CrossRef]
  39. Xu, Z.; Zhong, H.; Zhou, Q.; Du, Y.; Chen, L.; Zhang, Y.; Xue, C.; Cao, Y. A Highly Pathogenic Strain of Porcine Deltacoronavirus Caused Watery Diarrhea in Newborn Piglets. Virol. Sin. 2018, 33, 131–141. [Google Scholar] [CrossRef]
  40. Lednicky, J.A.; Tagliamonte, M.S.; White, S.K.; Elbadry, M.A.; Alam, M.M.; Stephenson, C.J.; Bonny, T.S.; Loeb, J.C.; Telisma, T.; Chavannes, S.; et al. Independent infections of porcine deltacoronavirus among Haitian children. Nature 2021, 600, 133–137. [Google Scholar] [CrossRef]
  41. Hsueh, F.-C.; Wu, C.-N.; Lin, M.Y.-C.; Hsu, F.-Y.; Lin, C.-F.; Chang, H.-W.; Lin, J.-H.; Liu, H.-F.; Chiou, M.-T.; Chan, K.R.; et al. Phylodynamic analysis and spike protein mutations in porcine deltacoronavirus with a new variant introduction in Taiwan. Virus Evol. 2021, 7, veab096. [Google Scholar] [CrossRef] [PubMed]
Pathogens 12 00686 g001 550
Figure 1. PDCoV affected pigs on the farm of isolate KPDCoV-2201 strain. (A) Black tarry diarrhea in sows; (B) watery diarrhea in piglets.
Figure 1. PDCoV affected pigs on the farm of isolate KPDCoV-2201 strain. (A) Black tarry diarrhea in sows; (B) watery diarrhea in piglets.
Pathogens 12 00686 g001
Pathogens 12 00686 g002 550
Figure 2. Phylogenetic trees of porcine deltacoronavirus (PDCoV) constructed using 128 global PDCoV full-length genomes (A) and 274 S genes (B). Red circles indicate the KPDCoV-2201 strain isolated in this study. Black circles indicate the previously reported Korean PDCoV strains. Genogroups (G) of PDCoV are indicated by different background colors of red (G1a), blue (G1b), and green (G2). The countries of the strains are listed on the right side of each strain with the indicated color. The right figure panel is enlarged to depict G1b, including Korean PDCoV strains. Scale bars indicate nt substitutions per site.
Figure 2. Phylogenetic trees of porcine deltacoronavirus (PDCoV) constructed using 128 global PDCoV full-length genomes (A) and 274 S genes (B). Red circles indicate the KPDCoV-2201 strain isolated in this study. Black circles indicate the previously reported Korean PDCoV strains. Genogroups (G) of PDCoV are indicated by different background colors of red (G1a), blue (G1b), and green (G2). The countries of the strains are listed on the right side of each strain with the indicated color. The right figure panel is enlarged to depict G1b, including Korean PDCoV strains. Scale bars indicate nt substitutions per site.
Pathogens 12 00686 g002
Pathogens 12 00686 g003 550
Figure 3. Time-scaled phylogenetic tree of G1b PDCoV strains constructed using BEAST. G1b PDCoV strains could be clustered into two clades based on the divergence of tMRCA. The black triangles represent the tMRCA of each clade. The value of each tMRCA is indicated on the node. Evolutionary clades of G1b PDCoV strains are indicated by different background colors of green (clade 1) and red (clade 2). The countries of the strains are presented on the tips of each strain with the displayed color. A red star indicates the KPDCoV-2201 strain isolated in this study.
Figure 3. Time-scaled phylogenetic tree of G1b PDCoV strains constructed using BEAST. G1b PDCoV strains could be clustered into two clades based on the divergence of tMRCA. The black triangles represent the tMRCA of each clade. The value of each tMRCA is indicated on the node. Evolutionary clades of G1b PDCoV strains are indicated by different background colors of green (clade 1) and red (clade 2). The countries of the strains are presented on the tips of each strain with the displayed color. A red star indicates the KPDCoV-2201 strain isolated in this study.
Pathogens 12 00686 g003
Table
Table 1. Details and diagnostic results of clinical samples.
Table 1. Details and diagnostic results of clinical samples.
SampleAgeResults of Real-Time Reverse Transcription PCR (qRT-PCR)
PEDVPEAVTGEVPToVPDCoVSADSPSVGARVGBRVGCRV
Feces4 daysN/AN/AN/AN/A24.75N/AN/AN/AN/AN/A
Feces4 daysN/AN/AN/AN/A21.08N/AN/AN/AN/AN/A
Feces4 daysN/AN/AN/AN/A19.04N/AN/AN/AN/AN/A
Feces4 daysN/AN/AN/AN/A19.87N/AN/AN/AN/AN/A
Feces4 daysN/AN/AN/AN/A21.84N/AN/AN/AN/AN/A
Feces4 daysN/AN/AN/AN/A21.96N/AN/AN/AN/AN/A
FecesSowN/AN/AN/AN/A30.32N/AN/AN/AN/AN/A
FecesSowN/AN/AN/AN/A32.29N/AN/AN/AN/AN/A
Intestine4 daysN/AN/AN/AN/A14.45N/AN/AN/AN/AN/A
Abbreviations: PEDV, porcine epidemic diarrhea virus; PEAV, porcine enteric alphacoronavirus; TGEV, transmissible gastroenteritis virus; PToV, porcine torovirus; PDCoV, porcine deltacoronavirus; SADS, swine acute diarrhea syndrome coronavirus; PSV, porcine sapelovirus; GARV, porcine group A rotavirus; GBRV, porcine group B rotavirus; GCRV, porcine group C rotavirus; N/A, not amplified.
Table
Table 2. Comparison of the full-length genome of Korean and global PDCoV strains with that of KPDCoV-2201.
Table 2. Comparison of the full-length genome of Korean and global PDCoV strains with that of KPDCoV-2201.
Strain/YearGenbank Accession NumberCountryFull-Length *5′-UTRORF1aORF1bSEMNS6NNS73′-UTR
% Identity to KPDCoV-2201
KNU14-04/2014KM820765Korea99.299.399.299.398.898.899.798.998.898.898.0
KNU16-07/2016KY364365Korea99.299.399.299.398.798.899.798.998.598.398.2
KNU16-11/2016KY926512Korea98.798.598.998.997.998.499.898.697.997.597.4
DH1/2016KY354363Korea99.199.199.299.298.698.899.799.698.898.798.2
DH2/2016KY354364Korea99.199.199.299.298.698.899.799.698.898.798.2
AH2004/2004KP757890China98.499.198.498.897.497.698.998.997.797.897.5
HKU15-155/2010NC_039208China98.698.898.798.997.698.499.198.298.398.298.0
Sichuan/2019MK993519China97.798.397.598.297.498.498.598.697.998.094.9
IL2768/2014KJ584355United States99.299.399.399.398.898.899.899.398.698.398.2
YMG/2014LC260044Japan99.2†99.399.498.898.899.898.998.898.898.2
HaNoi6/2015KX834351Vietnam97.898.096.997.896.198.898.998.698.297.798.0
S5015L/2015KU051649Thailand96.998.596.797.595.898.498.398.296.596.795.9
P1_16_BTL_0115/2016KX118627Laos97.098.396.797.696.098.498.398.297.897.295.9
Peru_isolate_2019/2019MT227371Peru99.2†99.399.398.699.210099.698.998.8†
NTU_C253_21/2021OM777140Taiwan99.0†99.199.198.598.899.898.998.898.598.5
* To avoid errors caused by UTR, the full-length genome identity was evaluated except for the UTR sequence. † UTR sequence was not fully sequenced.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, H.-R.; Park, J.; Lee, K.-K.; Jeoung, H.-Y.; Lyoo, Y.S.; Park, S.-C.; Park, C.-K. Genetic Characterization and Evolution of Porcine Deltacoronavirus Isolated in the Republic of Korea in 2022. Pathogens 2023, 12, 686. https://doi.org/10.3390/pathogens12050686

AMA Style

Kim H-R, Park J, Lee K-K, Jeoung H-Y, Lyoo YS, Park S-C, Park C-K. Genetic Characterization and Evolution of Porcine Deltacoronavirus Isolated in the Republic of Korea in 2022. Pathogens. 2023; 12(5):686. https://doi.org/10.3390/pathogens12050686

Chicago/Turabian Style

Kim, Hye-Ryung, Jonghyun Park, Kyoung-Ki Lee, Hye-Young Jeoung, Young S. Lyoo, Seung-Chun Park, and Choi-Kyu Park. 2023. "Genetic Characterization and Evolution of Porcine Deltacoronavirus Isolated in the Republic of Korea in 2022" Pathogens 12, no. 5: 686. https://doi.org/10.3390/pathogens12050686

APA Style

Kim, H. -R., Park, J., Lee, K. -K., Jeoung, H. -Y., Lyoo, Y. S., Park, S. -C., & Park, C. -K. (2023). Genetic Characterization and Evolution of Porcine Deltacoronavirus Isolated in the Republic of Korea in 2022. Pathogens, 12(5), 686. https://doi.org/10.3390/pathogens12050686

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop