Emerg Infect DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention20031041287435909-072610.3201/eid1601.090726ResearchRecombinant Canine Coronaviruses in Dogs, EuropeRecombinant Canine Coronaviruses in Dogs, EuropeDecaroNicolaMariVivianaEliaGabriellaAddieDiane D.CameroMicheleLucenteMaria StellaMartellaVitoBuonavogliaCanioFaculty of Veterinary Medicine of Bari, Valenzano, Italy (N. Decaro, V. Mari, G. Elia, M. Camero, M.S. Lucente, V. Martella, C. Buonavoglia)Feline Institute Pyrenees, Etchebar, France (D.D. Addie)Address for correspondence: Nicola Decaro, Department of Veterinary Public Health, Faculty of Veterinary Medicine of Bari, Strada per Casamassima km 3, 70010 Valenzano, Bari, Italy; e-mail: n.decaro@veterinaria.uniba.it120101614147

Subtype IIb originates from recombination with porcine transmissible gastroenteritis virus.

Coronaviruses of potential recombinant origin with porcine transmissible gastroenteritis virus (TGEV), referred to as a new subtype (IIb) of canine coronavirus (CCoV), were recently identified in dogs in Europe. To assess the distribution of the TGEV-like CCoV subtype, during 2001–2008 we tested fecal samples from dogs with gastroenteritis. Of 1,172 samples, 493 (42.06%) were positive for CCoV. CCoV-II was found in 218 samples, and CCoV-I and CCoV-II genotypes were found in 182. Approximately 20% of the samples with CCoV-II had the TGEV-like subtype; detection rates varied according to geographic origin. The highest and lowest rates of prevalence for CCoV-II infection were found in samples from Hungary and Greece (96.87% and 3.45%, respectively). Sequence and phylogenetic analyses showed that the CCoV-IIb strains were related to prototype TGEV-like strains in the 5′ and the 3′ ends of the spike protein gene.

Keywords: Canine coronavirusrecombinant strainsEuropegenetic analysisresearch

Coronaviruses (CoVs) (order Nidovirales, family Coronaviridae) are exceptionally prone to genetic evolution through accumulation of point mutations in genes encoding for structural and nonstructural proteins and homologous recombination among members of the same antigenic group (1). CoVs are organized by antigenic group. The first group is subdivided into subgroups 1a and 1b. Subgroup 1a includes highly related viruses (i.e., porcine transmissible gastroenteritis virus [TGEV] and its derivative porcine respiratory coronavirus [PRCoV], feline coronaviruses [FCoVs], and canine coronaviruses [CCoVs]) (2). According to a proposal by the Coronavirus Study Group of the International Committee of Taxonomy of Viruses, and given the virus’ close genetic relatedness (i.e., >96% aa identity in the key replicase 1ab domains), TGEV, FCoV, and CCoV should not be considered separate viruses. Instead, they should be considered host range variants of the same species (3).

CCoVs exemplify the genetic evolution and complexity of CoVs. To date, 2 CCoV genotypes are known, CCoV-I and CCoV-II (4); they share up to 96% of nucleotide identity in the viral genome (Lorusso et al., unpub. data) but are highly divergent in the spike protein gene (5). In addition, CCoV-I displays a novel open reading frame (ORF3) that encodes for a putative glycosylated protein, which is likely secreted from the infected cells (6). The 2 CCoV genotypes are commonly detected in the feces of dogs with diarrhea and often simultaneously infect the same dog (7). Both CCoV genotypes have been associated with mild clinical signs in pups, although hypervirulent strains have been reported to cause severe, fatal enteritis (810); a pantropic variant was responsible for systemic disease in natural and experimental conditions (1113).

It has been postulated that TGEV originated from CCoV-II through cross-species transmission, which is supported by the high genetic relatedness between the 2 viruses and by the presence of ORF3 remnants in CCoV-II and TGEV genomes (6). More recently, novel CCoV-II strains, which likely originated from a double recombination event with TGEV, occurring in the 5′ end of the spike protein gene, have been isolated (14,15). Accordingly, genotype II has been further divided into 2 subtypes, CCoV-IIa and CCoV-IIb, including extant and TGEV-like CCoVs, respectively (14).

To assess the distribution of the TGEV-like CCoVs in canine populations of different geographic areas of Europe, we used a subtype IIb-specific reverse transcription–PCR (RT-PCR) assay. In addition, we evaluated the genetic relationship among the identified strains by using sequence and phylogenetic analyses.

MethodsSample Origin

During 2001–2008, a total of 1,172 fecal samples were collected from dogs with acute enteritis in 14 European countries: Italy (n = 760), United Kingdom (n = 199), Greece (n = 81), Hungary (n = 41), Portugal (n = 33), Spain (n = 32), Belgium (n = 7), Romania (n = 6), Bulgaria (n = 3), Germany (n = 3), Sweden (n = 3), Slovakia (n = 2), Poland (n = 1), and Slovenia (n = 1). Some samples were retrieved from previous studies (7,1618), whereas the remaining samples were sent to our laboratory (from veterinarians, breeders, or other researchers) for routine diagnostic investigations.

RNA Extraction

Specimens were homogenized (10% wt/vol) in Dulbecco modified Eagle medium and subsequently clarified by centrifuging at 2,500 × g for 10 min. For RNA extraction, 140 μL of the supernatants were then used by means of QIAamp Viral RNA Mini Kit (QIAGEN S.p.A., Milan, Italy); according to the manufacturer’s protocol, RNA templates were stored at –70°C until use.

CCoV RNA Detection, Quantification, and Genotyping

For rapid detection and quantification of CCoV RNA, all RNA extracts were subjected to a previously established TaqMan-based real-time RT-PCR (16) with minor modifications. Briefly, a 1-step method was adopted by using the Platinum Quantitative PCR ThermoScript One-Step System (Invitrogen S.R.L., Milan, Italy) and the following 50-µL mixture: 25 µL of master mix, 300 nM of primers CCoV-forward (5′-TTGATCGTTTTTATAACGGTTCTACAA-3′) and CCoV-reverse (5′-AATGGGCCATAATAGCCACATAAT-3′), 200 nM of probe CCoV-Pb (5′-FAM-ACCTCAATTTAGCTGGTTCGTGTATGGCATT-BHQ1-3′), and 10 μL of template RNA. To obtain a standard curve for absolute quantification, we simultaneously analyzed duplicates of log10 dilutions of standard RNA (16). The thermal profile consisted of reverse transcription at 50°C for 20 min, activation of Platinum Taq DNA polymerase at 95°C for 2 min, 45 cycles of denaturation at 95°C for 15 s, annealing at 48°C for 30 s, and extension at 60°C for 30 s.

The positive samples were characterized by 2 distinct genotype-specific assays (17) performed by using the Platinum Quantitative PCR ThermoScript One-Step System (Invitrogen S.R.L.) and the following oligonucleotide sets (final concentrations were 600 and 200 nM for primers and probes, respectively): primer pair CCoVI-F (5′-CGTTAGTGCACTTGGAAGAAGCT-3′)/CCoVI-R (5′-ACCAGCCATTTTAAATCCTTCA-3′) and probe CCoVI-Pb (5′-FAM -CCTCTTGAAGGTACACCAA-TAMRA-3′) for CCoV-I; primer pair CCoVII-F (5′-TAGTGCATTAGGAAGAAGCT-3′)/CCoVII-R (5′-AGCAATTTTGAACCCTTC-3′) and probe CCoVII-Pb (5′-FAM-CCTCTTGAAGGTGTGCC-TAMRA-3′) for CCoV-II. The thermal protocol was as described for CCoV detection except for different annealing temperatures (i.e., 53°C and 48°C for CCoV-I and CCoV-II, respectively).

Development of RT-PCRs Specific for Classical and TGEV-like CCoVs

Considering the high divergence observed in the 5′ end of the spike gene between classical (subtype IIa) and TGEV-like (subtype IIb) CCoVs, specific CCoV-IIa and CCoV-IIb gel-based RT-PCRs were developed. Primer 20179 (sense, 5′-GGCTCTATCACATAACTCAGTCCTAG-3′) binds a conserved region at the 3′ end of ORF1b and was recruited from a previous study (13), whereas antisense primers INS-R-dg (5′-GCTGTAACATAKTCRTCATTCCAC-3′) and 174-268 (5′-CAACATGTAACCTTTGTCTGTGATCTGC-3′) target regions at the 5′ end of the spike protein gene of feline CoV-II (FCoV-II)/classical CCoV-II and TGEV/TGEV-like CCoV, respectively. Separate RT-PCRs with primers 20179/INS-R (CCoV-IIa) or 20179/174-268 (CCoV-IIb) were conducted by using SuperScript One-Step RT-PCR for Long Templates (Invitrogen S.R.L.), according to the manufacturer’s instructions. The following thermal protocol was used: reverse transcription at 50°C for 30 min, inactivation of Superscript II RT at 94°C for 2 min, 40 cycles of 94°C for 30 s, 55°C for 30 s, 68°C for 30 s, and final extension at 68°C for 10 min. The PCR products were detected by using electrophoresis through a 1.5% agarose gel and examination under UV light after ethidium bromide staining.

RT-PCR for Amplification of the 3′ End of the Spike Protein Gene of CCoV-II

To rule out any potential infection by true TGEV strains and to confirm the recombinant origin of the TGEV-like CCoVs, we submitted 20 samples that were positive for CCoV-IIb and negative for CCoV-I to RT-PCR amplification of the 3′ end of the spike protein gene of CCoV-II (7). Primers S5 (5′-TGCATTTGTGTCTCAGACTT-3′) and S6 (5′-CCAAGGCCATTTTACATAAG-3′) were used in the RT-PCR, performed according to the protocol described for CCoV subtyping.

RT-PCR of the ORF7a/7b Region

To rule out the presence of true TGEV strains in the dog feces that were positive by CCoV-IIb–specific assay, we used an RT-PCR that had been proven to discriminate between TGEV and CCoV on the basis of amplicon size (19). In fact, primers N3SN (5′-GTGTTTGATGACACACAGGTTGAG-3′) and R3AS (5′-GCTTACCATTCTGTACAAGAGGTAG-3′) target the 3′ end of the viral genome, where CCoV/FCoV and TGEV display 2 (ORFs 7a and 7b) and 1 (ORF7) accessory genes, respectively. As controls, the following reference group-1a CoVs were used: TGEV-Purdue (kindly provided by P. Cordioli, Istituto Zooprofilattico Sperimentale di Lombardia ed Emilia Romagna, Brescia, Italy), FCoV-I-249/04 (20), FCoV-II-29/92 (21), CCoV-I-Elmo/02 (5), CCoV-IIa-CB/05 (11), CCoV-IIb-341/05, CCoV-IIb-174/06, CCoV-IIb-430/07, and CCoV-IIb-119/08 (14).

Sequence and Phylogenic Analyses

The RT-PCR products obtained with primer pairs 20179/174-268 and S5/S6 from 26 samples having positive CCoV-IIb–specific assay results and being representative of the different geographic areas were subjected to direct sequencing at the BaseClear B.V. (Leiden, the Netherlands). The sequences were manually edited and analyzed by using BioEdit software (22) and National Center for Biotechnology Information (www.ncbi.nlm.nih.gov) and European Molecular Biology Laboratory (www.ebi.ac.uk) analysis tools. Phylogenetic and molecular evolutionary analyses were conducted by using Mega 4.1, beta (23). Phylogenetic trees on the basis of partial 5′ (339-nt) and 3′ (520-nt) ends of the spike protein gene were elaborated by using parsimony and neighbor-joining methods, which supplied statistical support with bootstrapping >1,000 replicates. Group-2 CoV canine respiratory CoV-240/05 (24) was used as an outgroup. The nucleotide sequences of the analyzed CCoV-IIb strains were deposited in GenBank under accession nos. GQ130243–GQ130268 and GQ148749–GQ148774 for 5′ and 3′ ends of the spike gene, respectively.

ResultsCCoV Detection and Quantification

CCoV RNA was detected in 493 (42.06%) of 1,172 fecal samples from dogs with diarrhea. Viral RNA titers ranged from 1.25 × 101 to 7.56 × 107 copies/μL of template. For each geographic region of origin, detection rates of the CCoV real-time RT-PCR were Italy 330/760 (43.42%), United Kingdom 54/199 (27.13%), Greece 45/81 (55.5%), Hungary 32/42 (78.05%), Portugal 12/33 (36.36%), Spain 2/32 (6.25%), Belgium 4/7 (57.14), Romania 4/6 (66.66%), Bulgaria 1/3 (33.33%), Germany 3/3 (100%), Sweden 3/3 (100%), Slovakia 1/2 (50%), Poland 1/1 (100%), and Slovenia 1/1 (100%) (Table 1).

Distribution of CCoV genotypes and TGEV-like strains (CCoV-IIb) in dogs in European countries, 2001–2008*
CountryNo. samplesNo. or % positive
CCoV, no.CCoV-I no.†CCoV-II, no.†CCoV-I+II, no.‡CCoV-IIb, no.CCoV-IIb/II, %
Italy760330541381383412.32
United Kingdom199541324171024.39
Greece814516151413.45
Hungary413203023196.87
Portugal331263300
Spain32201100
Belgium7413000
Romania64112133.33
Bulgaria3101000
Germany3300300
Sweden33111150
Slovakia2100100
Poland1101000
Slovenia
1
1
1
0
0
0
0
Total1,172493932181827819.50

*CCoV, canine coronavirus; TGEV, porcine transmissible gastroenteritis virus.
†Infections with a single CCoV genotype.
‡Infections with both CCoV genotypes.

CCoV Genotype and Subtype Distribution

The geographic distribution of the CCoV types and subtypes is reported in Table 1. Genotype-specific amplification assays showed that 93 (18.86%) of 493 CCoV-positive samples were positive for CCoV-I, and 218 (44.22%) were positive for CCoV-II. In addition, 182 samples (36.92%) were positive for both genotypes. Both genotypes were found to circulate in most European countries that had been sampled, and an overall prevalence of CCoV-II was found in all countries except Greece, where the 2 genotypes were detected approximately to the same extent.

By using the developed TGEV-like RT-PCR, we found that 78 (19.50%) of 400 samples containing CCoV-II strains, alone or in combination with CCoV-I, were positive for CCoV-IIb. The remaining 322 CCoV-II strains were positive for subtype IIa, whereas mixed infections caused by both CCoV-II subtypes were not detected in any samples. Almost all CCoV–IIb-positive samples were from Italy (34/276; 12.32%), the United Kingdom (10/41; 24.39%), and Hungary (31/32; 96.87%). Single CCoV-IIb strains were detected in samples from Greece (1/29; 3.45%), Romania (1/3; 33.33%), and Sweden (1/2; 50%), whereas no samples from the other countries had TGEV-like CCoVs.

RT-PCR of ORF7a/7b Region

After RT-PCR with primer pair N3SN/R3AS, CCoV and FCoV reference strains yielded an amplicon >1,000 bp, with the exception of TGEV-like CCoV 341/05, which gave a 929-bp product as a consequence of a 154-nt deletion in ORF7b (14). In contrast, a 367-bp product was obtained from TGEV-Purdue, as previously described (19). This pattern of amplification agreed with the absence of ORF7b in TGEV (25). All 20 samples that were positive according to the CCoV-IIb–specific assay and negative according to the CCoV-I TaqMan assay were confirmed to contain true CCoV strains because they yielded RT-PCR products considerably larger than the 367-bp amplicon obtained from TGEV.

Sequence and Phylogenetic Analyses

We selected the following TGEV-like strains for sequence analysis of the 5′ and 3′ ends of the spike protein gene: 12 strains from Italy, 8 from the United Kingdom, and 5 from Hungary, plus the single strains from Greece, Romania, and Sweden. All RT-PCR products were sequenced except those obtained from the samples from Greece and Sweden, which yielded weak bands despite the considerable viral RNA titers in the original fecal samples (5.02 × 104 and 5.76 × 106 RNA copies/μL of template, respectively).

The obtained sequences were compared with each other (Table 2) and with 3 CCoV-IIa, 4 CCoV-IIb, 2 FCoV-II, 2 CCoV-I, 3 FCoV-I, and 3 TGEV reference sequences (Table 3). Sequence comparison of the TGEV-like strains showed overall nucleotide identity of 83.6%–99.6% and 92.7%–100% in the 5′ and 3′ ends of the spike gene, respectively. By analyzing the TGEV-like strains by country of origin, we found the highest genetic variability among the viruses from Italy (86.8%–99.3% and 94.5%–99.8% of nucleotide identity in the gene 5′ and 3′ ends, respectively), whereas the strains from Hungary showed the highest relatedness (92.7%–99.3% and 98.4%–99.6% of nucleotide identity in the gene 5′ and 3′ ends, respectively) (Table 2). The TGEV-like strains exhibited the best identity to prototype strains (14) in both the 5′ (85.9%–99.3%) and the 3′ end (92.7%–99.8%) of the spike gene, whereas a slightly lower identity was found to the old strain UCD1 (19) in the 5′ end, which is the only sequence available in the GenBank database. When the 5′ end was analyzed, the identified TGEV-like CCoVs were more related to classical TGEVs (76.6%–84.4%) than to type IIa CCoVs (34.2%–38.3%). In contrast, analysis of the 3′ end of the spike gene showed nucleotide identities of 87.1%–93.1% to TGEV and of 90.8%–99.0% to CCoV-IIa. The best identities among CCoV-IIa isolates were to strain Insavc-1, which has been proposed as intermediate virus between CCoV and TGEV (26) and to the more recent pantropic strain CB/05 (11) in the 5′ and 3′ ends, respectively.

Percent nucleotide identities in the 5′ (top right) and 3′ (bottom left) ends of the spike gene of CCoV-IIb (TGEV-like) strains from dogs in Europe, 2001–2008*
CountryItalyUnited KingdomHungaryRomaniaOverall
Italy86.8–99.3
94.5–99.883.6–92.786.5–99.089.1–99.383.6–99.3
United Kingdom92.7–97.088.5–99.6
97.6–99.887.8–97.087.8–93.183.6–99.6
Hungary94.7–10092.9–94.992.7–99.3
98.4–99.694.7–99.386.5–99.3
Romania95.3–99.693.5–94.598.6–99.6ND87.8–99.3
Overall92.7–10092.7–99.892.9–10093.5–99.683.6–99.6
92.7–100

*CCoV, canine coronavirus; TGEV, porcine transmissible gastroenteritis virus; ND, not done because only a single strain from Romania was analyzed.

Percent nucleotide identity of representative TGEV-like CCoVs to group-1a CoV reference strains in the 5′ and 3′ ends of the spike gene from dogs from Europe, 2001–2008*
CoV strainGenBank accession no.% Identity
5′ end†3′ end‡
CCoV-IIb-341/05EU85636188.1–99.394.7–99.4
CCoV-IIb-174/06EU85636286.8–99.093.3–97.0
CCoV-IIb-430/07EU92479085.9–97.392.7–97.0
CCoV-IIb-119/08EU92479189.1–99.093.7–99.8
CCoV-IIb-UCD1AF11624876.6–85.7NA
TGEV-PurdueNC_00230676.6–84.491.4–93.1
TGEV-TSDQ20144776.2–84.490.8–92.5
TGEV-96-1933AF10442074.0–81.887.1–88.8
PRCoV-RM4Z24675NA90.2–91.6
PRCoV-86-137004X60089NA90.2–91.6
CCoV-IIa-Insavc-1D1309635.4–38.392.1–95.3
CCoV-IIa-CB/05DQ11222635.1–37.793.7–99.0
CCoV-IIa-BGF10AY34216034.2–37.590.8–94.3
FCoV-II-79-1146NC_00702534.0–36.690.2–91.6
FCoV-II-79-1683X8079934.6–37.590.4–91.9
CCoV-I-Elmo/02AY30702037.6–40.461.8–62.8
CCoV-I-23/03AY30702131.8–34.860.3–61.2
FCoV-I-KU-2D3204437.4–40.859.1–60.3
FCoV-I-BlackEU18607236.2–39.358.9–61.0
FCoV-I-UCD1AB08822238.3–42.757.9–58.9

*TGEV, porcine transmissible gastroenteritis virus; CCoV, canine coronavirus; NA, sequence not available.
†Nucleotides 1–305 of the spike gene of reference strain CCoV-IIb-341/05.
‡Nucleotides 3574–4085 of the spike gene of reference strain CCoV-IIb-341/05.

At the phylogenetic level, the sequenced strains were grouped in the same cluster with TGEV and prototype CCoV-IIb strains at the 5′ end of the spike protein gene, displaying an obvious distance to both CCoV-IIa/FCoV-II and CCoV-I/FCoV-I (Figure, panel A). The TGEV-like CCoVs from the United Kingdom formed a unique clade, whereas the strains detected in eastern Europe were mixed with CCoV-IIb viruses from Italy. The old TGEV-like strain, UCD1, clustered with TGEV isolates. At the 3′ end of the same gene, subtype IIa and IIb strains segregated together and were separated from the FCoVs and the porcine CoVs TGEV and PRCoV (Figure, panel B). In addition, the strains from the United Kingdom were again grouped in a separate subcluster.

Phylogenetic analysis of canine coronavirus (CCoV) type IIb. Maximum parsimony trees based on partial 5´ (A) and 3´ (B) ends of the spike protein gene of group-1a coronaviruses (CoVs). For phylogenetic tree construction, the following reference strains were used (GenBank accession numbers are in parentheses): porcine transmissible gastroenteritis virus (TGEV) Purdue (NC_002306), TS (DQ201447), 96-1933 (AF104420), porcine respiratory coronavirus (PRCoV, only 3´ end) RM-4 (Z24675), 86-137004 (X60089), CCoV-IIa CB/05 (DQ112226), Insavc-1 (D13096), BGF10 (AY342160), CCoV-IIb 341/05 (EU856361), 174/06 (EU856362), 430/07 (EU924790), 118/08 (EU924791), UCD-1 (AF116248, only 5´ end), CCoV-I Elmo/02 (AY307020), 23/03 (AY307021), feline coronavirus (FCoV) type I Black (EU186072), KU-1 (D32044), UCD-1 (AB088222), FCoV-II 79-1146 (NC_007025), 79-1183 (X80799). The tree is rooted on the group-2 CoV canine respiratory coronavirus (CRCoV) 240/05 (EU999954). Statistical support was provided by bootstrapping >1,000 replicates. Scale bars indicate estimated numbers of nucleotide substitutions per site.

Discussion

CoVs are exceptionally prone to variability through accumulation of point mutations and recombination events. A CCoV strain displaying close relatedness to porcine CoVs in the N-terminus of the spike protein, which is related to porcine CoVs, was isolated ≈20 years ago, but the subsequent molecular characterization was restricted to the 5′ end of the spike gene (19). Additional TGEV-like CCoVs were reported more recently in Italy (14) and the United Kingdom (15). However, full biological and molecular characterization was carried out only for the isolates from Italy, showing that the TGEV-like strains are likely recombinant with TGEV at the level of the 5′ end of the spike gene (14). Experimental infection of CCoV-seronegative beagle pups showed that TGEV-like (i.e., subtype IIb) CCoV induces clinical signs resembling those of classical (i.e., subtype IIa) CCoVs, that is, mild diarrhea for a few days. (14) Unlike pantropic CCoV (1113), TGEV-like CCoV was not able to spread systemically. Of the 4 recombinant strains detected in Italy, 2 had originated from eastern Europe, but at the phylogenetic level they were mixed with strains from Italy. The prevalence of this CCoV subtype in the canine population has not been determined in previous studies.

Our epidemiologic investigation assessed the distribution of CCoV-IIb in the dog population of Europe. Approximately 50% of the analyzed samples were positive for CCoV, showing the presence of CCoV-I or CCoV-II. Mixed infections caused by both genotypes were detected in <40% of the CCoV-positive samples, considerably lower than previously reported percentages (7,17). Approximately 20% of the CCoV–II-positive samples contained TGEV-like strains (Table 1). However, the prevalence of this CCoV subtype differed by geographic origin of the samples; the highest detection rates (96.87%) were observed in Hungary. The recombinant origin of all strains characterized by RT-PCR was confirmed by sequence analysis of 5′ and 3′ ends of the spike gene and by RT-PCR of the ORF7a/7b region. The selected 26 TGEV-like strains were related to prototype strains from Italy (14) in both the 5′ and the 3′ ends of the spike gene. A comparison with the prototype UK strains reported by Erles and Brownlie (15) was not possible, however, because the unique spike sequence deposited in the GenBank is located more downstream of the gene with respect to the sequences that we obtained.

On the basis of the spike gene sequences, the strains from Italy and eastern Europe were closely related, whereas the strains from the United Kingdom were more genetically distant (Table 3). This pattern of segregation was confirmed by phylogenic analysis, which showed that viruses detected in the United Kingdom formed a separate cluster with respect to the samples from Italy, Hungary, and Romania (Figure). The genetic relatedness between the TGEV-like strains from Italy and those from eastern Europe may be accounted for by extensive dog importation to Italy. In addition, dog exchange between eastern Europe and Italy has been associated with the reemergence of canine infectious hepatitis (27) and the spread of the arctic lineage of canine distemper virus (28).

In the 5′ end of the spike gene, the old TGEV-like strain, UCD1, was found to be genetically more related to true TGEV isolates than to recent CCoV-IIb strains, thereby accounting for recombination events occurring at different times. With the exception of strain UCD1, analysis of archival samples found the oldest TGEV-like strain in 2001, about 4 years before this CCoV subtype was reported in Italy (14).

Our study confirms that recombinant CCoVs are effectively circulating in dogs in different European countries. Considering the genetic distance in the spike protein, this circulation questions the efficacy of vaccines, which are based on classical (CCoV-IIa) strains, against the emerging TGEV-like (CCoV-IIb) viruses. Only vaccination trials and subsequent challenges by TGEV-like strains might assess whether the poor cross-reactivity between CCoV-IIa and CCoV-IIb observed in a previous study (14) might affect the immune response of dogs against the recombinant viruses.

Suggested citation for this article: Decaro N, Mari V, Elia G, Addie DD, Camero M, Lucente MS, et al. Recombinant canine coronaviruses in dogs, Europe. Emerg Infect Dis [serial on the Internet] 2010 Jan [date cited]. Available from http://www.cdc.gov/EID/content/16/1/41.htm

Acknowledgments

We are grateful to Lucrezia Favia for her excellent assistance with part of the experimental work.

Grants from the Italian Ministry of Health, Ricerca Finalizzata 2007 project, “Mammalian coronaviruses: molecular epidemiology, vaccine development, and implications for animal and human health,” supported this work.

Dr Decaro is an associate professor at the Faculty of Veterinary Medicine of Bari, Italy. His research interests include viral pathogens in dogs, with particular emphasis on canine coronavirus and parvovirus.

ReferencesEnjuanes L, Brian D, Cavanagh D, Holmes K, Lai MMC, Laude H, Coronaviridae. In: van Regenmortel MHV, Fauquet CM, Bishop DHL, Carstens EB, Estes MK, Lemon SM, et al., editors. Virus taxonomy. Classification and nomenclature of viruses. New York: Academic Press; 2000 p. 835–9.Gorbalenya AE Genomics and evolution of the Nidovirales. In: Perlman S, Gallagher T, Snijder EJ, editors. Nidoviruses. Washington: ASM Press; 2008 p. 15–28.de Groot RJ, Ziebuhr J, Poon LL, Woo PC, Talbot P, Rottier PJ M, et al. Revision of the family Coronaviridae Taxonomic proposal of the Coronavirus Study Group to the ICTV Executive Committee; 2008 [cited 2009 Oct 20]. Available from http://talk.ictvonline.org/media/p/1230.aspxDecaro N, Buonavoglia C An update on canine coronaviruses: viral evolution and pathobiology.Vet Microbiol 2008;132:22134 10.1016/j.vetmic.2008.06.00718635322Pratelli A, Martella V, Decaro N, Tinelli A, Camero M, Cirone F, Genetic diversity of a canine coronavirus detected in pups with diarrhoea in Italy.J Virol Methods 2003;110:917 10.1016/S0166-0934(03)00081-812757915Lorusso A, Decaro N, Schellen P, Rottier PJ, Buonavoglia C, Haijema BJ, Gain, preservation and loss of a group 1a coronavirus accessory glycoprotein.J Virol 2008;82:103127 10.1128/JVI.01031-0818667517Pratelli A, Decaro N, Tinelli A, Martella V, Elia G, Tempesta M, Two genotypes of canine coronavirus simultaneously detected in fecal samples of dogs with diarrhea.J Clin Microbiol 2004;42:17979 10.1128/JCM.42.4.1797-1799.200415071054Sanchez-Morgado JM, Poynter S, Morris TH Molecular characterization of a virulent canine coronavirus BGF strain.Virus Res 2004;104:2731 10.1016/j.virusres.2004.02.03815177889Evermann JF, Abbott JR, Han S Canine coronavirus–associated puppy mortality without evidence of concurrent canine parvovirus infection.J Vet Diagn Invest 2005;17:610416475526Escutenaire S, Isaksson M, Renström LH, Klingeborn B, Buonavoglia C, Berg M, Characterization of divergent and atypical canine coronaviruses from Sweden.Arch Virol 2007;152:150714 10.1007/s00705-007-0986-117533554Buonavoglia C, Decaro N, Martella V, Elia G, Campolo M, Desario C, Canine coronavirus highly pathogenic for dogs.Emerg Infect Dis 2006;12:492416704791Decaro N, Martella V, Elia G, Campolo M, Desario C, Cirone F, Molecular characterisation of the virulent canine coronavirus CB/05 strain.Virus Res 2007;125:5460 10.1016/j.virusres.2006.12.00617275120Decaro N, Campolo M, Lorusso A, Desario C, Mari V, Colaianni ML, Experimental infection of dogs with a novel strain of canine coronavirus causing systemic disease and lymphopenia.Vet Microbiol 2008;128:25360 10.1016/j.vetmic.2007.10.00818031959Decaro N, Mari V, Campolo M, Lorusso A, Camero M, Elia G, Recombinant canine coronaviruses related to transmissible gastroenteritis virus of swine are circulating in dogs.J Virol 2009;83:15327 10.1128/JVI.01937-0819036814Erles K, Brownlie J Sequence analysis of divergent canine coronavirus strains present in a UK dog population.Virus Res 2009;141:215 10.1016/j.virusres.2008.12.00919162099Decaro N, Pratelli A, Campolo M, Elia G, Martella V, Tempesta M, Quantitation of canine coronavirus RNA in the faeces of dogs by TaqMan RT-PCR.J Virol Methods 2004;119:14550 10.1016/j.jviromet.2004.03.01215158596Decaro N, Martella V, Ricci D, Elia G, Desario C, Campolo M, 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.J Virol Methods 2005;130:728 10.1016/j.jviromet.2005.06.00516024100Decaro N, Desario C, Addie DD, Martella V, Vieira MJ, Elia G, Molecular epidemiology of canine parvovirus in Europe.Emerg Infect Dis 2007;13:1222417953097Wesley RD The S gene of canine coronavirus, strain UCD-1, is more closely related to the S gene of transmissible gastroenteritis virus than to that of feline infectious peritonitis virus.Virus Res 1999;61:14552 10.1016/S0168-1702(99)00032-510475084Campolo M, Desario C, Ricci D, Elia G, Decaro N, Martella V, Identification of a feline coronavirus type I strain from a cat with feline infectious peritonitis by RT-PCR and phylogenetic analysis.New Microbiol 2005;28:1273316035257Buonavoglia C, Sagazio P, Cirone F, Tempesta M, Marsilio F Isolamento e caratterizzazione di uno stipite di virus della peritonite infettiva felina.Veterinaria.1995;9:914Hall TA BioEdit: a user-friendly biological sequence alignment and analysis program for Windows 95/98/NT.Nucleic Acids Symposium Series.1999;41:958Kumar S, Dudley J, Nei M, Tamura K MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences.Brief Bioinform 2008;9:299306 10.1093/bib/bbn01718417537Lorusso A, Desario A, Mari V, Campolo M, Lorusso E, Elia G, Molecular characterization of a canine respiratory coronavirus strain detected in Italy.Virus Res 2009;141:96100 10.1016/j.virusres.2008.12.01119162098de Groot RJ, Andeweg AC, Horzinek MC, Spaan WJ Sequence analysis of the 3′-end of the feline coronavirus FIPV 79-1146 genome: comparison with the genome of porcine coronavirus TGEV reveals large insertions.Virology 1988;167:37063201747Horsburgh BC, Brierley I, Brown TD Analysis of a 9.6-kb sequence from the 3′ end of canine coronavirus genomic RNA.J Gen Virol 1992;73:284962 10.1099/0022-1317-73-11-28491431811Decaro N, Campolo M, Elia G, Buonavoglia D, Colaianni ML, Lorusso A, Infectious canine hepatitis: an “old” disease reemerging in Italy.Res Vet Sci 2007;83:26973 10.1016/j.rvsc.2006.11.00917197003Martella V, Elia G, Lucente MS, Decaro N, Lorusso E, Banyai K, Genotyping canine distemper virus (CDV) by a hemi-nested multiplex PCR provides a rapid approach for investigation of CDV outbreaks.Vet Microbiol 2007;122:3242 10.1016/j.vetmic.2007.01.00517275219