key: cord-267189-xq70rn1j
authors: Wang, Xinyu; Li, Chunqiu; Guo, Donghua; Wang, Xinyu; Wei, Shan; Geng, Yufei; Wang, Enyu; Wang, Zhihui; Zhao, Xiwen; Su, Mingjun; Liu, Qiujin; Zhang, Siyao; Feng, Li; Sun, Dongbo
title: Co-Circulation of Canine Coronavirus I and IIa/b with High Prevalence and Genetic Diversity in Heilongjiang Province, Northeast China
date: 2016-01-15
journal: PLoS One
DOI: 10.1371/journal.pone.0146975
sha: 
doc_id: 267189
cord_uid: xq70rn1j

To trace the evolution of canine coronavirus (CCoV), 201 stool samples from diarrheic dogs in northeast China were subjected to reverse transcription-polymerase chain reactions (RT-PCRs) targeting the partial M and S genes of CCoV, followed by an epidemiological analysis. M gene RT-PCRs showed that 28.36% (57/201) of the samples were positive for CCoV; of the 57 positive samples, CCoV-I and CCoV-II accounted for 15.79% (9/57) and 84.21% (48/57), respectively. A sequence comparison of the partial M gene revealed nucleotide homologies of 88.4%–100% among the 57 CCoV strains, and 88.7%–96.2% identity between the 57 CCoV strains and the Chinese reference strain HF3. The CCoV-I and CCoV-II strains exhibited genetic diversity when compared with reference strains from China and other countries. The 57 CCoV strains exhibited high co-infection rates with canine kobuvirus (CaKV) (33.33%) and canine parvovirus-2 (CPV-2) (31.58%). The CCoV prevalence in diarrheic dogs differed significantly with immunization status, regions, seasons, and ages. Moreover, 28 S genes were amplified from the 57 CCoV-positive samples, including 26 CCoV-IIa strains, one CCoV-IIb strain, and one CCoV-I strain. A sequence comparison of the partial S gene revealed 86.3%–100% nucleotide identity among the 26 CCoV-IIa strains, and 89.6%–92.2% identity between the 26 CCoV-IIa strains and the Chinese reference strain V1. The 26 CCoV-IIa strains showed genetic diversity when compared with reference strains from China and other countries. Our data provide evidence that CCoV-I, CCoV-IIa, and CCoV-IIb strains co-circulate in the diarrhoetic dogs in northeast China, high co-infection rates with CaKV and CPV-2 were observed, and the CCoV-II strains exhibited high prevalence and genetic diversity.

Canine coronavirus (CCoV) was first recognized as an enteric pathogen of dogs in 1971 [1] . CCoV is a common infection in young dogs, particularly those housed in large groups [2] [3] [4] [5] . CCoV is an enveloped, single-stranded, positive-sense RNA virus, and it belongs to the family Coronaviridae, subfamily Coronavirinae, genus Alphacoronavirus, species Alphacoronavirus-1 [6] . CCoV consists of two distinct genotypes, CCoV-I and CCoV-II; the CCoV-II viruses are further divided into two subtypes CCoV-IIa and CCoV-IIb [3, 7] . The S protein of CCoV is a glycoprotein peplomer on the viral surface, and it plays an important role in the induction of neutralizing antibodies, specific receptor binding, and cell membrane fusion. Accumulating reports attributed that the increase in the severity of CCoV infections in dogs and the emergence of CCoV variants to potential recombination events within the S gene, which occur when a host is co-infected with different CCoV types [8] [9] [10] . Therefore, CCoV has received much attention as an emerging cause of infectious disease in dogs [4, [11] [12] [13] [14] [15] [16] .

CCoV infection is a leading causes of diarrhea in dog population in China. Wang et al. (2006) reported that CCoV-II infections were very common in domestic dog, fox, and raccoon-dog populations in China [17] . Ma et al. (2008) reported the molecular characterization of the 9.36-kb 3 0 region of the CCoV 1-71 strain [18] . Gao et al. (2009) reported the isolation and identification of a CCoV strain from giant pandas in China [19] . However, information about the epidemiology of CCoVs in China is not available in the past five years. In the current study, we conducted a molecular epidemiologic investigation of CCoV in Heilongjiang province, northeast China. Moreover, the genetic evolution and co-infection of the identified CCoV strains were analyzed. Our aim was to provide insights into the epidemiology and genetic diversity of the CCoV strains circulating in northeast China.

The animal experiment, sampling, was approved by the Animal Care and Use Committee of the Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences, China. The sampling and data publication also were approved by animal's owners. The field study did not involve endangered or protected species. No specific permissions were required for locations of samples because the samples were collected from public areas or non-protection areas.

A total of 201 fecal samples were collected in the form of rectal swabs of dogs with diarrhea from animal hospitals in the Harbin, Daqing, and Mudanjiang districts of Heilongjiang province in northeast China from May 2014 to April 2015, using 3.5-ml commercial virus sampling tubes (YOCON Biological Technology Co. Ltd., Beijing, China). For all samples, animal age, animal breed, animal gender, collection date, and vaccination were recorded, respectively. Of the 201 samples, 141 were collected in Harbin, 20 were collected in Daqing, and 40 were collected in Mudanjiang. All rectal swab samples were stored at −80°C, and they were also used for etiological investigations in our other studies [20, 21] .

After 1 mL of fecal samples was centrifuged at 1,500 × g for 10 min at 4°C, the supernatant of each sample was transferred to a 1.5-ml Eppendorf tube. Viral RNA was extracted from each sample using the TIANamp Virus RNA Kit (Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer's instructions. The extracted RNA were stored at -80°C.

Molecular detection of CCoV was conducted using reverse transcription-polymerase chain reaction (RT-PCR) targeting a 409-bp fragment of the M gene of CCoV that was described by Pratelli et al. (1999) [22] . The S gene fragments used for CCoV genotyping/subtyping, including a 346-bp fragment of CCoV-I, a 694-bp fragment of CCoV-IIa, and a 370-bp fragment of CCoV-IIb, were amplified using the primers EL1F/EL1R, S5F/S6R, and CEPol-1/TGSP-2, respectively [23, 24] . Briefly, first-strand cDNA was synthesized using random primers (six nucleotides) using Moloney murine leukemia virus (RNase H-) reverse transcriptase (Novoprotein Scientific Inc., Shanghai, China) according to the manufacturer's instructions. In this study, the EmeraldAmp PCR Master Mix (2× Premix) (TaKaRa Biotechnology Co., Ltd., Dalian, China), and the Applied Biosystems GeneAmp PCR System 9700 thermal cycler (Thermo Fisher Scientific, Waltham, MA, USA) were used for PCR amplifications of all target genes. Other PCR amplification conditions were performed according to the protocols described by [23] and Erles and Brownlie (2009) [24] . The purified PCR products were directly subjected to Sanger sequencing, and each sample was sequenced three times. Sequence analysis was performed using the EditSeq program in the Lasergene DNAS-TAR™ version 5.06 software (DNASTAR Inc., Madison, WI, USA). All nucleotide sequences generated in our study were submitted to GenBank under accession numbers KT192642-KT192698 for the M gene of the 57 CCoV strains, KT222969-KT222994 for the S gene of the 26 CCoV-IIa strains, KT222995 for the S gene of the CCoV-I strain, and KT222996 for the S gene of the CCoV-IIb strain.

For the phylogenetic analysis, partial sequences of the M and S genes of CCoV strains, including CCoV-I, CCoV-IIa, and CCoV-IIb strains, were retrieved from GenBank. To construct the phylogenetic trees, a multiple alignment of all target sequences was performed using ClustalX program version 1.83. Furthermore, phylogenetic trees of all target sequences were generated from the ClustalX-generated alignments by MEGA6.06 software using the neighbor-joining method [25] . Neighbor-joining phylogenetic trees were built with the p-distance model, 1000 bootstrap replicates, and, otherwise, the default parameters in MEGA 6.

The 57 CCoV-positive samples were screened for CPV-2, CaKV, canine astrovirus (CaAstV), canine norovirus (CNoV), canine bocavirus (CBoV), and Group A-Rotavirus (CRV-A) by either PCR or RT-PCR, followed by sequencing, according to previously described protocols [26] [27] [28] [29] [30] . All nucleotide sequences generated in our study were submitted to GenBank, and the accession numbers of the target genes of co-infected viruses are shown in S1 Table. Results

The characteristics of the 57 CCoV-positive dogs, the genotyping of 57 CCoV strains, and the amino acid substitutions of M protein of 57 CCoV strains are shown in S1 Table, and a further analysis of the 57 CCoV-positive samples is shown in Table 1 . Nucleotide sequences of the partial M gene of the 57 CCoV strains identified in our study were shown in Supporting Information (S1 Fig). Of the 201 samples, 57 samples (28.36%) were positive for CCoV following RT-PCR amplification of the partial M gene, and the CCoV-positive rates of the Harbin, Daqing, and Mudanjiang districts were 26.95%, 50%, and 22.5%, respectively (Table 1 ). Of the 57 CCoV-positive samples, CCoV-I and CCoV-II accounted for 15.79% (9/57) and 84.21% (48/ 57), respectively; immunized animals and non-immunized animals accounted for 43.86% and 36.84%, respectively; 56.14% (32/57) were collected from October to December, and 66.66% (38/57) were collected from dogs aged from 2 to 4 months; the total co-infection rate of CCoV with any viruses reached 56.14%, and CaKV and CPV-2 accounted for 33.33% (19/57) and 31.58% (18/57), respectively (Table 1) . Of the 57 CCoV strains, only 28 S gene sequences were successfully amplified, of which 26 belonged to CCoV-IIa, one to CCoV-IIb, and one to CCoV-I (S1 Table) . Nucleotide sequences of the partial S gene of CCoV strains identified in our study were shown in Supporting Information (S2 Fig).

The sequence comparison of the partial M gene revealed nucleotide homologies of 88.4%-100% and amino acids homologies of 93%-100% among the 57 CCoV strains, while nucleotide and amino acid homologies of 88.7%-96.2% and 92.2%-97.4%, respectively, were observed between the 57 CCoV strains and the Chinese reference strain HF3. Of the 57 CCoV strains, the nine CCoV-I strains exhibited 94.8%-100% nucleotide identities and 94.8%-100% amino acid identities, and the 48 CCoV-II strains showed 96.8%-100% nucleotide identities and 98.3%-100% amino acid identities (Table 1) . A total of 12 amino acid substitutions were found in the partial M protein (S1 Table) ; the four substitutions, Thr178Val, Met198Ile, His206Asn, and Gln228Lys, occurred in all identified CCoV-II strains (48/48), and two substitutions, Ile198Met and Asn206His, occurred in all identified CCoV-I strains (9/9) ( Table 1 ). The sequence comparison of the partial S gene revealed nucleotide homologies of 86.3%-100% and amino acids homologies of 87.6%-100% among the the 26 CCoV-IIa strains, and nucleotide and amino acid homologies of 89.6%-92.2% and 91.1%-97.5%, respectively, between the 26 CCoV-IIa strains and the Chinese reference strain V1 (Table 1) . A phylogenetic analysis using the nucleotide sequences of the M gene revealed that the nine CCoV-I strains formed two clusters (C1 and C2), while the 48 CCoV-II strains formed six clusters (C1-C6). The 48 CCoV-II strains only exhibited a close relationship to one Chinese reference strain, NJ17, and differed genetically from most of the reference strains from China and other countries (Fig 1) . A phylogenetic analysis using partial S gene sequences demonstrated that the 26 CCoV-IIa strains were closely related to three reference strains, 5281 (Japan), TN-449 (USA), and 1086-IIa (Brazil), and differed genetically from reference strains from China and other countries (Fig 2A) . The CCoV-I strain HRB-A4 showed a close relationship to the reference strain from Italy (Fig 2B) , while the CCoV-IIb strain HRB-BB9 was closely related to reference strains from European countries (Fig 2C) .

In our study, 57 of 201 samples (28.36%) were positive for CCoV, and the CCoV-positive rate differed among the three districts of Heilongjiang province (26.95% for Harbin, 50% for Daqing, and 22.5% for Mudanjiang). The CCoV-positive rate in feces has been reported to be 43.75% [17] , the Republic of Korea [31] , Japan [32] , Albania [33] , Brazil [11] , Italy, Belgium, the Netherlands, Germany, the United Kingdom, Spain, and France, respectively [34] . These data demonstrated that CCoV infections showed clear differences among the geographical regions. The total CCoV-positive rate reported here was lower than that in a previous report in 2006 in China [17] . In our study, CCoV-II, accounting for 84.21% of the 57 CCoV strains, was the predominant CCoV type in northeast China, which is line with most reports [11, 17, 34] .

Mixed infections of canine enteric viruses frequently occur in diarrheic dogs. In our study, the total co-infection rate of CCoV with one or more CaKV, CPV-2, and CBoV strains was 56.14%; single co-infections with CaKV, CPV-2, and CBoV accounted for 33.33%, 31.58%, and 5.26%, respectively. Co-infections between CCoV and CPV-2 have been documented in dog populations in Western Europe, Japan, and Albania [32] [33] [34] . However, the high co-infection rate of CCoV with CaKV has not been reported in China and other countries. CaKV had been reported to be associated with severe enteritis in a litter of puppies [13] . It is speculated that the high prevalence of co-infections of CCoV with CaKV and CPV-2 may be associated with the occurrence of viral diarrhea in dogs in northeast China. Further studies should be conducted to understand the effect of the high co-infection rate with CaKV on the severity of clinical symptoms of CCoV infections.

In our study, CCoV-positive rate showed clear differences among seasons and ages. The high prevalence of CCoV in the feces of diarrheic dogs aged 2-4 months was partially validated in other studies [11, 32, 33] . The high positive rate (56.14%) of CCoV was found between October and December, which may be associated with seasonal variations in northeast China. At present, most dogs in northeast China are vaccinated for CPV-2, canine distemper virus, canine parainfluenza virus, canine adenovirus type 1, and canine adenovirus type 2. In our study, of the 57 CCoV-positive samples, the immunized and non-immunized animals accounted for 43.86% and 36.84%, respectively. This result demonstrated that vaccination for other canine viruses did not effect CCoV infections in the dog population in Heilongjiang province, northeast China. Given the high prevalence and co-infection rates of CCoV, although it is controversial whether CCoV vaccines provide adequate immunity [35, 36] , the immunization for CCoV is recommended in the dog population in northeast China in the future.

Fragments of the M and S genes have been used to genotype CCoV [4, 11, 24, 37] . In our study, the CCoV-I and CCoV-II genotypes of the 57 CCoV strains were successfully identified demonstrated that CCoV-I, CCoV-IIa, and CCoV-IIb co-circulated in northeast China, and that CCoV-IIa strains predominated in the dog population.

The current study is the first to reveal that CCoV-I, CCoV-IIa, and CCoV-IIb strains co-circulate in northeast China, and that there are high co-infection rate with CaKV and CPV-2. CCo-V-II strains are predominant, and they exhibit genetic diversity. Resulting data increase our understanding of the evolution of CCoV strains circulating in northeast China, and they provide valuable epidemiological information for further studies of CCoV. However, further studies are necessary to clarify the effect of co-infections and amino acid substitutions on the severity of the clinical symptoms of CCoV infections.

Supporting Information 

Recovery and characterization of a coronavirus from military dogs with diarrhea

Canine coronavirus infections in Japan: virological and epidemiological aspects

An update on canine coronaviruses: viral evolution and pathobiology

Canine coronavirus, Greece. Molecular analysis and genetic diversity characterization

Cross sectional and longitudinal surveys of canine enteric coronavirus infection in kennelled dogs: a molecular marker for biosecurity

Virus taxonomy: classification and nomenclature of viruses. Ninth report of the International Committee on Taxonomy of Viruses

Feline and canine coronaviruses: common genetic and pathobiological features

Canine enteric coronaviruses: emerging viral pathogens with distinct recombinant spike proteins

Recombinant canine coronaviruses related to transmissible gastroenteritis virus of Swine are circulating in dogs

Characterization of a recombinant canine coronavirus with a distinct receptor-binding (S1) domain

Molecular characterization of canine coronavirus strains circulating in Brazil

European surveillance for pantropic canine coronavirus

Canine kobuviruses in diarrhoeic dogs in Italy

Characterization of pantropic canine coronavirus from Brazil

Emergence of pathogenic coronaviruses in cats by homologous recombination between feline and canine coronaviruses

Full-length genome analysis of canine coronavirus type I

Detection of canine coronaviruses genotype I and II in raised Canidae animals in China

Molecular characterization of the 9.36 kb C-terminal region of canine coronavirus 1-71 strain

Isolation and identification of a canine coronavirus strain from giant pandas (Ailuropoda melanoleuca)

Prevalence and phylogenetic analysis of canine kobuviruses in diarrhoetic dogs in northeast China

Co-Circulation of the Rare CPV-2c with Unique Gln370Arg Substitution, New CPV-2b with Unique Thr440Ala Substitution, and New CPV-2a with High Prevalence and Variation in Heilongjiang Province, Northeast China

Development of a nested PCR assay for the detection of canine coronavirus

Two genotypes of canine coronavirus simultaneously detected in fecal samples of dogs with diarrhea

Sequence analysis of divergent canine coronavirus strains present in a UK dog population

MEGA6: Molecular evolutionary genetics analysis version 6.0

Antigenic analysis of canine parvovirus strains isolated in Italy

Rotavirus A genotype P[4]G2: genetic diversity and reassortment events among strains circulating in Brazil between 2005 and

Identification and characterization of bocaviruses in cats and dogs reveals a novel feline bocavirus and a novel genetic group of canine bocavirus

Molecular epidemiology of canine norovirus in dogs from Portugal

Phylogenetic analysis of astrovirus and kobuvirus in Korean dogs

M gene analysis of canine coronavirus strains detected in Korea

Detection and genotyping of canine coronavirus RNA in diarrheic dogs in Japan

Detection and genetic characterization of canine parvovirus and canine coronavirus strains circulating in district of Tirana in Albania

Western European epidemiological survey for parvovirus and coronavirus infections in dogs

Efficacy of an inactivated canine coronavirus vaccine in pups

Safety and efficacy of a modified-live canine coronavirus vaccine in dogs

Identification of coronaviruses in dogs that segregate separately from the canine coronavirus genotype

Genotype-specific fluorogenic RT-PCR assays for the detection and quantitation of canine coronavirus type I and type II RNA in faecal samples of dogs

Canine coronavirus highly pathogenic for dogs

Canine coronavirus: not only an enteric pathogen

Isolation, tissue distribution and molecular characterization of two recombinant canine coronavirus strains