- Research article
- Open Access
First detection and genetic characterization of canine Kobuvirus in domestic dogs in Thailand
BMC Veterinary Researchvolume 15, Article number: 254 (2019)
Canine Kobuvirus (CaKoV) has been detected both in healthy and diarrheic dogs and in asymptomatic wild carnivores. In this study, we conducted a survey of CaKoV at small animal hospitals in Bangkok and vicinity of Thailand during September 2016 to September 2018.
Three hundred and seven rectal swab samples were collected from healthy dogs (n = 55) and dogs with gastroenteritis symptoms (n = 252). Of 307 swab samples tested by using one-step RT-PCR specific to 3D gene, we found CaKoV positivity at 17.59% (54/307). CaKoVs could be detected in both sick (19.44%) and healthy (9.09%) animals. In relation to age group, CaKoV could be frequently detected in younger dogs (25.45%). Our result showed no seasonal pattern of CaKoV infection in domestic dogs. In this study, we characterized CaKoVs by whole genome sequencing (n = 4) or 3D and VP1 gene sequencing (n = 8). Genetic and phylogenetic analyses showed that whole genomes of Thai CaKoVs were closely related to Chinese CaKoVs with highest 99.5% amino acid identity suggesting possible origin of CaKoVs in Thailand.
In conclusion, this study was the first to report the detection and genetic characteristics of CaKoVs in domestic dogs in Thailand. CaKoVs could be detected in both sick and healthy dogs. The virus is frequently detected in younger dogs. Thai CaKoVs were genetically closely related and grouped with Chinese CaKoVs. Our result raises the concerns to vet practitioners that diarrhea in dogs due to canine Kobuvirus infection should not be ignored.
Kobuvirus (KoV) is a single-strand positive-sense RNA virus. KoV belongs to the family Piconaviridae, genus Kobuvirus, which consists of four species Aichivirus A, B, C and D [1,2,3]. KoV has been reported in feces of several mammal species including humans, ruminants, pigs, dogs, cats, bats and rodents [3,4,5,6,7,8,9,10]. The Kobuvirus species Aichivirus A contains four types including Aichi virus 1, canine Kobuvirus 1 (CaKoV), Feline Kobuvirus 1 (FeKoV) and Murine Kobuvirus 1 (MuKoV). Canine Kobuvirus 1 (CaKoV) was first reported in dogs with acute gastroenteritis in the US in 2011 [5, 11]. CaKoV was subsequently reported in dogs in UK, Italy, Australia, Japan, Korea and China [4, 12,13,14,15]. The virus was reported in wild carnivores (Jackal and Hyena) and domestic dogs in Tanzania, Africa , in foxes in Spain  and in foxes  and wolves in Italy . Several studies have reported the detection of CaKoV infection in dogs with or without diarrhea and sometime systemic infection . To date, only 12 completed CaKoV genomes are available in the GenBank database.
During September 2016 to September 2018, the center of excellence for emerging and re-emerging infectious diseases in animals (CUEIDAs), Chulalongkorn University conducted a survey of canine Kobuvirus in domestic dogs at small animal hospitals in 5 provinces of Thailand. The survey was conducted under the Chulalongkorn University’s animal use and care protocol # 1731074. The result of this study provided the first detection and genetic characterization of CaKoV isolated from domestic dogs in Thailand.
Canine Kobuviruses in domestic dogs in Thailand
During September 2016 to September 2018, we conducted a survey of viral enteric diseases in domestic dogs in small animal hospitals in 5 provinces of Thailand (Bangkok, Nakhon Ratchasima, Ratchaburi, Suphanburi, and Tak). We tested 307 rectal swab samples for CaKoV by using one-step RT-PCR specific to 3D gene. Based on a two-year survey, we found CaKoV positivity at 17.59% (54/307). CaKoVs could be detected in both sick (19.44% (49/252)) and healthy (9.09% (5/55)) animals. Our result showed no seasonal pattern of CaKoV infection in dogs (Figs. 1 and 2). In relation to age group, CaKoV could be frequently detected in younger dogs at 25.45% (42/165) (Additional file 2: Table S2). The co-infections of CaKoV with other enteric viral pathogens were observed including CaKoV/Canine parvovirus/Canine Coronavirus (n = 6), CaKoV/Canine parvovirus (n = 20) and CaKoV/Canine Coronavirus (n = 2). In this study, 12 CaKoVs were selected and characterized by whole genome sequencing (n = 4) or 3D and VP1 gene sequencing (n = 8). The viruses were selected to represent epidemiological and demographic data such as age, date of isolation and breed. In this study, nucleotide sequences of the CaKoV were submitted to the GenBank database under the accession numbers MK201776 - MK201795 (Table 1).
Phylogeny of the Thai canine Kobuviruses
Phylogenetic analysis of whole genome of CaKoVs showed that the Thai CaKoVs were closely related to each other and clustered with Aichivirus A. The cluster Aichivirus A contains Kobuviruses from dogs, cats, rodents, bats and human. While Aichivirus B and C contain Kobuviruses from cattle and pigs, respectively. Based on whole genome sequence, Thai CaKoVs were closely related to Chinese CaKoVs sub-cluster but in separated sub-cluster from the viruses from the US, UK, Brazil and Tanzania (Fig. 3). Phylogenetic analysis of 3D and VP1 of Thai CaKoVs and reference CaKoVs from various animal species were also performed. Similarly, 3D gene of Thai CaKoVs were grouped together with Chinese CaKoVs (G1 sub-cluster) but separated from the viruses in sub-clusters G2 as well as G3 (Fig. 4). Phylogenetic analysis of VP1 gene, the viruses can be clustered into 2 major subgroups, US/EU/Africa subgroup and China/Thailand subgroup (Fig. 5).
Genetic analysis of the Thai canine Kobuviruses
We compared the nucleotide and deduced amino acid sequences of Thai CaKoVs against those of reference viruses from the US, UK, Italy, China, and Korea (Tables 2 and 3). Our results showed that whole genome of 4 Thai CaKoVs (CU-53, CU-101, CU-249 and CU-716) shared 96.7–99.3% nucleotide similarity (99.6–100% amino acid similarity) to each other and posed highest nucleotide similarity to Chinese CaKoVs including SMCD-59 (97.0% nt and 99.5% aa identity) and CH-1 (96.8% nt and 98.7% aa identity). Our analysis showed that the VP1 protein was the most diverse gene with 93.4–99.9% nucleotide similarity (96.9–100% aa similarity) among Thai CaKoVs and 82.2–96.8% with other reference CaKoVs. The most variable region of VP1 is position 201–243, especially proline rich region. Putative proline rich region at VP1–228-240 (P228XPPPPXPPXPXP240) was also observed in Thai CaKoVs as well as reference viruses (Table 4). In this study, unique amino acids were found in Thai and Chinese CaKoVs at the position, 65 V, 67D, 119L, 138T, 150P, 151M, 153D, 201S, 204Q, 205Q, 201Q, 213T and 241E (Table 4). Analysis of predicted amino acid cleavage sits of whole genome were conserved among Thai CaKoVs (Table 5).
Canine Kobuvirus (CaKoV) is an emerging pathogen in Thailand. To the best of our knowledge, the CaKoV was described in Asia in retrospective study in Korea in 2011 and have been reported in Japan, China and Australia, respectively [2, 15, 17, 21]. However, the CaKoV have never been reported in the country or South East Asia region. In this study, during the 2 year-survey program, we found CaKoV positivity at 17.59% in both sick (19.44%) and healthy (9.09%) animals. Compare to other studies, CaKoV % positivity in this study was lower than those in China (54%) and Korea (32.2%) [14, 22]. Our result showed that the CaKoV could be frequently detected in younger dogs at 27% which consistence with previous reports . Similar to other previous studies, co-infections with other enteric viral pathogens were observed such as CaKoV/Canine parvovirus and CaKoV/Canine Coronavirus [12, 14, 15]. Moreover, CaKoVs were detected in both diarrheic and non-diarrheic dogs which consistent with other studies [2, 15]. Our result supported that this virus may not be the only cause of enteric disease in dogs. Nevertheless, the CaKoV infection have still been identified in symptomatic dogs without other enteric pathogen infections . Our observation supported that the role of CaKoV as a primary pathogen of acute gastroenteritis remain unclear.
In this study, the genome size of 4 Thai CaKoVs is 7,530 bp with one ORF encoding 2,444 amino acids of a putative polyprotein, which comparable to previous reports. Genome organization of CaKoV includes leader protein (L), structural proteins (VP0, VP3, VP1), non-structural proteins (2A, 2B, 2C, 3A, 3B, 3C, 3D). Phylogenetic analyses showed that the Thai CaKoVs were closely related to each other and clustered with Aichivirus A. It is noted that Thai CaKoVs were closely related to Chinese CaKoVs sub-cluster but in separated sub-cluster from the viruses from the US, UK, Brazil and Tanzania (Fig. 3). Phylogenetic analyses of 3D gene showed similar result which Thai CaKoVs were grouped together with Chinese CaKoVs (G1 sub-cluster). This observation regarding to the sub-clusters of CaKoVs was in agreement with the previous study . On the other hand, based on VP1 gene, the viruses can be clustered into 2 major subgroups, US/EU/Africa subgroup and China/Thailand subgroup which similar to the previous reports [16, 22] (Figs. 4 and 5).
Genetic analyses of Thai CaKoVs showed that whole genome of 4 Thai CaKoVs posed highest nucleotide similarity to Chinese CaKoVs including SMCD-59 and CH-1. This observation supported phylogenetic analysis that Thai CaKoVs were closely related to Chinese CaKoVs sub-cluster but in separated sub-cluster from the viruses from the US, UK, Brazil and Tanzania. Of all viral genes, the VP1 gene was the most diverse gene among Thai CaKoVs and other reference CaKoVs. Similar observation was also reported in previous study that VP1 protein is the most variable capsid protein . It is noted that the putative proline rich region at VP1–228-240 (P228XPPPPXPPXPXP240) was observed both in Thai CaKoVs and reference viruses. Previous studies indicated that proline rich region may associate with enteric receptor binding of the viruses [14, 24]. It is noted that Thai CaKoVs posed unique PPP (VP1; 228–240), which also observed most reference viruses from China, Korea, Japan, US, UK suggesting unique characteristic. These unique amino acids were not observed in the CaKoV from the Australia (CE9), Brazil (BRA/26) and Tanzania (TZ/75, TZ82) [16, 20]. However, the association of these unique amino acids and viral pathogenesis is still need to be further investigated. Based on genetic analysis, unique amino acids at the position, 65 V, 67D, 119L, 138 T, 150P, 151M, 153D, 201S, 204Q, 205Q, 201Q, 213 T and 241E were observed. These unique amino acids of China/Thailand sub-cluster could be benefit for the detection of virus origin or diagnostic purpose in the future. Similar to previous study, analysis of predicted amino acid cleavage sits of whole genome were conserved among CaKoVs except one variation at 776/777 (VP3/VP1) which unique in wild carnivores .
In conclusion, this study is the first to report of canine Kobuvirus in dogs in Thailand. CaKoVs were mostly detected in clinical dogs of young age. However, the viruses could be detected from both healthy and sicked dogs. Genetic and phylogenetic analyses showed that whole genomes of Thai CaKoVs were closely related to Chinese CaKoVs in 2015 (SMCD-59) with high nucleotide similarity suggesting a possible origin of CaKoVs in Thailand. CaKoV is considered as an emerging viral pathogen in the domestic dogs. Since CaKoVs have never been reported in the country and SEA region, the detection and characterization of CaKoV from different parts of the regions should be extended for better understanding the epidemiology and evolution of CaKoVs. Our result raises the concerns to vet practitioners that diarrhea in dogs due to canine Kobuvirus infection should not be ignored.
Sample collection was conducted in domestic dogs at small animal hospitals in Bangkok and vicinity of Thailand During September 2016 to September 2018. 307 rectal swab samples were collected from healthy dogs (n = 55) and dogs with gastroenteritis symptoms (n = 252) including vomiting, watery diarrhea, hemorrhagic diarrhea and dehydration. The swab samples were collected from dogs of young age (< 1 year) (n = 165), adult (1–5 years) (n = 98) and older (> 5 years) (n = 44). The animal demographic data including age, sex, breed, and vaccination history were also recorded. The ethics was conducted under the Chulalongkorn University’s animal use and care protocol # 1731074. The consent to participate of the owners of the animals used in this study was obtained in writing.
Canine Kobuvirus (CaKoV) detection
All 307 samples were subjected to canine Kobuvirus identification by one step RT-PCR using primers specific to 3D gene of CaKoV . First, RNA extraction was performed using the QIAsymphony DSP Viral/Pathogen mini kit (Qiagen, Hilden, Germany) following manufacturer’s instructions. To detect CaKoV, RNA samples were screened for 3D gene of CaKoV by using one step RT-PCR assay. The primers used in this study were previously described including U1F (5′-CATGCTCCTCGGTGGTCTCA-3′) and U1R (5′-GTCCGGGTCCATCACAGGGT -3′) . Briefly, one-step RT-PCR was conducted in a total final volume of 25 μl comprising 3 μl of template RNA, 15 μl of 2xReaction Mix (Invitrogen, USA), 0.6 μl of 10 μM forward and reverse primers, 1.2 μl of SuperScript III RT (Invitrogen, USA) and distilled water to final volume 25 μl. The condition of RT-PCR assay included cDNA synthesis step at 55 °C for 30 min, next to an initial denaturation step at 94 °C for 2 min, following 40 cycles of denaturation at 94 °C for 30 s, annealing at 52 °C for 30 s and extension at 68 °C for 1 min, as well as, final extension step at 68 °C for 5 min. To confirm CaKoV, 4 μl of PCR products were run on a 1.5% agarose gel, which mixed with Red Safe at 100 V for 45 min. The expected size of CaKoV positive amplified products was 631 bp. Due to dogs showed clinical signs similar to other canine viral enteric diseases, all samples were also tested for Canine Parvovirus (n = 307), Canine Rotavirus (n = 307) and Canine Coronavirus (n = 30) [25,26,27].
Canine Kobuvirus characterization
In this study, four CaKoV positive samples (CU-53, CU-101, CU-247 and CU-716) were selected for whole genome sequencing and additional eight CaKoV positive samples were selected for 3D and VP1 gene sequencing. The CaKoVs were selected based on epidemiological and demographic data such as age, date of isolation, breed, and vaccination history. For sequencing, nucleotide sequences of each gene of the viruses were amplified by new primer sets designed by using Primer 3 plus program . List of oligonucleotide primers is provided in Additional file 1: Table S1 In brief, PCR was proceed in a final volume of 30 μl containing 2 μl of cDNA, 0.4 μM of each forward and reverse primer, 1X TopTaq Master Mix, 1X Coral Load, and distilled water. The PCR condition was set as initial denaturation at 94 °C for 3 min; 40 cycles of denaturation at 94 °C for 30 s, annealing at 50–55 °C for 45 s, extension at 72 °C for 1–1.30 min; and final extension at 72 °C for 7 min. PCR products were then purified and sequenced (1st Base Laboratories Sdn Bhd, Malaysia). Nucleotide sequences were edited, validated and assembled by using SeqMan software v.5.03 (DNASTAR Inc.; Wisconsin, USA).
Phylogenetic and genetic analyses of canine Kobuviruses
The phylogenetic and genetic analyses were performed by comparing nucleotide sequences of Thai CaKoVs with those of Kobuvirus available from the GenBank database. The reference nucleotide sequences of CaKoVs were retrived based on their different geographic locations, host species and date of isolation. Phylogenetic analysis of CaKoV was performed by using MEGA v.6.0 (Tempe, AZ, USA)  with neighbor-joining method with Kimura 2-parameter with 1,000 bootstrap replicates and Beast program with Bayesian Markov chain Monte Carlo (BMCMC) with 10,000,000 generations and an average standard deviation of split frequencies < 0.05 . For genetic analysis, the nucleotide sequences and deduced amino acids of CaKoV were aligned and compared using MegAlign software v.5.03 (DNASTAR Inc.; Wisconsin, USA). Pairwise comparison of nucleotides and amino acids of Thai CaKoV and those of reference CaKoVs were conducted. The variable and unique amino acids related to receptor binding of the viruses and host preferences of CaKoVs were monitored.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and supplement tables.
Adams MJ, Lefkowitz EJ, King AM, Harrach B, Harrison RL, Knowles NJ, Kropinski AM, Krupovic M, Kuhn JH, Mushegian AR, et al. Ratification vote on taxonomic proposals to the international committee on taxonomy of viruses (2016). Arch Virol. 2016;161(10):2921–49.
Oem JK, Lee MH, Lee KK, An DJ. Novel Kobuvirus species identified from black goat with diarrhea. Vet Microbiol. 2014;172(3–4):563–7.
Yamashita T, Ito M, Kabashima Y, Tsuzuki H, Fujiura A, Sakae K. Isolation and characterization of a new species of kobuvirus associated with cattle. J Gen Virol. 2003;84(Pt 11):3069–77.
Carmona-Vicente N, Buesa J, Brown PA, Merga JY, Darby AC, Stavisky J, Sadler L, Gaskell RM, Dawson S, Radford AD. Phylogeny and prevalence of kobuviruses in dogs and cats in the UK. Vet Microbiol. 2013;164(3–4):246–52.
Kapoor A, Simmonds P, Dubovi EJ, Qaisar N, Henriquez JA, Medina J, Shields S, Lipkin WI. Characterization of a canine homolog of human Aichivirus. J Virol. 2011;85(21):11520–5.
Khamrin P, Maneekarn N, Kongkaew A, Kongkaew S, Okitsu S, Ushijima H. Porcine kobuvirus in piglets. Emerg Infect Dis. 2009;15(12):2075–6.
Li L, Victoria JG, Wang C, Jones M, Fellers GM, Kunz TH, Delwart E. Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses. J Virol. 2010;84(14):6955–65.
Lu G, Zhang X, Luo J, Sun Y, Xu H, Huang J, Ou J, Li S. First report and genetic characterization of feline kobuvirus in diarrhoeic cats in China. Transbound Emerg Dis. 2018;65(5):1357–63.
Mohamed FF, Mansour SMG, Orabi A, El-Araby IE, Ng TFF, Mor SK, Goyal SM. Detection and genetic characterization of bovine kobuvirus from calves in Egypt. Arch Virol. 2018;163(6):1439–47.
Phan TG, Kapusinszky B, Wang C, Rose RK, Lipton HL, Delwart EL. The fecal viral flora of wild rodents. PLoS Pathog. 2011;7(9):e1002218.
Li L, Pesavento PA, Shan T, Leutenegger CM, Wang C, Delwart E. Viruses in diarrhoeic dogs include novel kobuviruses and sapoviruses. J Gen Virol. 2011;92(Pt 11):2534–41.
Di Martino B, Di Felice E, Ceci C, Di Profio F, Marsilio F. Canine kobuviruses in diarrhoeic dogs in Italy. Vet Microbiol. 2013;166(1–2):246–9.
Kong N, Zuo Y, Wang Z, Yu H, Zhou EM, Shan T, Tong G. Molecular characterization of new described kobuvirus in dogs with diarrhea in China. Springerplus. 2016;5(1):2047.
Oem JK, Choi JW, Lee MH, Lee KK, Choi KS. Canine kobuvirus infections in Korean dogs. Arch Virol. 2014;159(10):2751–5.
Soma T, Matsubayashi M, Sasai K. Detection of kobuvirus RNA in Japanese domestic dogs. J Vet Med Sci. 2016;78(11):1731–5.
Olarte-Castillo XA, Heeger F, Mazzoni CJ, Greenwood AD, Fyumagwa R, Moehlman PD, Hofer H, East ML. Molecular characterization of canine kobuvirus in wild carnivores and the domestic dog in Africa. Virology. 2015;477:89–97.
Bodewes R, Ruiz-Gonzalez A, Schapendonk CM, van den Brand JM, Osterhaus AD, Smits SL. Viral metagenomic analysis of feces of wild small carnivores. Virol J. 2014;11:89.
Di Martino B, Di Profio F, Melegari I, Robetto S, Di Felice E, Orusa R, Marsilio F. Molecular evidence of kobuviruses in free-ranging red foxes (Vulpes vulpes). Arch Virol. 2014;159(7):1803–6.
Melegari I, Sarchese V, Di Profio F, Robetto S, Carella E, Bermudez Sanchez S, Orusa R, Martella V, Marsilio F, Di Martino B. First molecular identification of kobuviruses in wolves (Canis lupus) in Italy. Arch Virol. 2018;163(2):509–13.
Ribeiro J, Headley SA, Diniz JA, Pereira AH, Lorenzetti E, Alfieri AA, Alfieri AF. Extra-intestinal detection of canine kobuvirus in a puppy from southern Brazil. Arch Virol. 2017;162(3):867–72.
Choi JW, Lee MH, Lee KK, Oem JK. Genetic characteristics of the complete feline kobuvirus genome. Virus Genes. 2015;50(1):52–7.
Li M, Yan N, Wang M, Zhang B, Yue H, Tang C. Prevalence and genomic characteristics of canine kobuvirus in Southwest China. Arch Virol. 2018;163(2):459–66.
Li CQ, Wei S, Guo DH, Wang ZH, Geng YF, Wang EY, Zhao XW, Su MJ, Wang XY, Sun DB. Prevalence and phylogenetic analysis of canine kobuviruses in diarrhoetic dogs in Northeast China. J Vet Med Sci. 2016;78(1):7–11.
Chen L, Zhu L, Zhou YC, Xu ZW, Guo WZ, Yang WY. Molecular and phylogenetic analysis of the porcine kobuvirus VP1 region using infected pigs from Sichuan Province. Virol J. 2013;10:281.
Pratelli A, Tempesta M, Greco G, Martella V, Buonavoglia C. Development of a nested PCR assay for the detection of canine coronavirus. J Virol Methods. 1999;80(1):11–5.
Buonavoglia C, Martella V, Pratelli A, Tempesta M, Cavalli A, Buonavoglia D, Bozzo G, Elia G, Decaro N, Carmichael L. Evidence for evolution of canine parvovirus type 2 in Italy. J Gen Virol. 2001;82(Pt 12):3021–5.
Ortega AF, Martinez-Castaneda JS, Bautista-Gomez LG, Munoz RF, Hernandez IQ. Identification of co-infection by rotavirus and parvovirus in dogs with gastroenteritis in Mexico. Braz J Microbiol. 2017;48(4):769–73.
Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007;23(10):1289–91.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.
We would like to thank the staffs of the Center of Excellence for Emerging and Re-emerging Infectious Diseases in Animals, Department of Veterinary Public Health for sample collection and data analysis.
This project was financial supported by the research fund under the 90th Anniversary Chulalongkorn University (Ratchadaphiseksomphot Endowment Fund) (GCUGR1125614077D). Chulalongkorn University provided financial support to the Center of Excellence for Emerging and Re-emerging Infectious Diseases in Animals for study design, sample collection, analysis and interpretation. The Thailand Research Fund supported the Royal Golden Jubilee (RGJ) Ph.D. program, for first author scholarship (RGJ-PHD/0056/2557) and TRF Senior Scholar to the corresponding author (RTA6080012).
Ethics approval and consent to participate
Ethics and consent to participate in the study was conducted under the Chulalongkorn University’s animal use and care protocol (IACUC) # 1731074.
Consent for publications
The consent to participate of the owners of the animals used in this study was obtained in writing.
All authors in this paper declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.