Previous Article | Next Article ![]()
Journal of Clinical Microbiology, October 2005, p. 5009-5017, Vol. 43, No. 10
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.10.5009-5017.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Veterinary Medical Diagnostic Laboratory and the Department of Veterinary Pathobiology, University of Missouri, College of Veterinary Medicine, Columbia, Missouri 65211
Received 5 April 2005/ Returned for modification 27 May 2005/ Accepted 19 July 2005
|
|
|---|
|
|
|---|
15.7-kb RNA genome and an enveloped virus particle that is 150 to 300 nm in diameter (29). The pathogenesis of CDV infection in dogs has been well characterized (2, 22, 29, 35). The genome of CDV encodes the following virion proteins: matrix (M), fusion (F), hemagglutinin (H), nucleocapsid (N), polymerase (L), and phosphoprotein (P). The H gene protein is responsible for viral attachment to the cell host (38) and may also play a role in induction of protective immunity (11). The H protein is one of the most variable morbillivirus proteins and thus has often been used to assess genetic changes between CDV isolates (8, 9, 17, 18, 20, 21, 23, 26, 28, 30). In this study, postmortem findings, immunohistochemical labeling, and reverse transcription (RT)-PCR specific for CDV genomic RNA established a definitive diagnosis of disease caused by CDV infection in four naturally occurring canine cases in Missouri that occurred between June and October 2004. Approximately 4.5 kb of genomic nucleotide sequence was determined for each of these four cases. Phylogenetic analysis performed on complete (F and H genes) and partial (P gene) gene sequences demonstrated that three genetically distinct CDV strains were detected and that these putative isolates were distant from those previously found in North America and most closely related to isolates previously described in either Asia or Europe.
|
|
|---|
Nucleic acid extraction and RT-PCR amplification. Unfixed tissues were collected at necropsy and homogenized in phosphate-buffered saline, pH 7.4. The RNA extractions were performed with the NucleoSpin RNA II kit (BD Biosciences, Inc., Palo Alto, CA) according to the manufacturer's instructions. For each RNA extraction step, strict protocols were followed to avoid cross-contamination of samples (e.g., sample extraction was performed for samples from only one case at a time, aerosol barrier pipette tips were used for all pipetting steps, gloves were changed between extractions of cases, and a negative extraction control was incorporated with every set of samples). The RNA was stored at 80°C until it was used as the template for RT-PCR amplification.
For the initial diagnosis of CDV infection, amplification of RNA by RT-PCR was performed with oligonucleotide primers previously reported for diagnosis of CDV infection in dogs (15). For nucleotide sequencing, RT-PCR amplification of RNA purified directly from tissues collected at necropsy was performed with the oligonucleotide primers shown in Table 1 and targeted to the P, F, or H gene. All RT-PCR amplifications were performed using each primer at a final concentration of 0.6 µM in a 20-µl reaction mixture with 0.8 µl of QIAGEN one-step RT-PCR enzyme mix in the manufacturer's buffer containing deoxynucleoside triphosphates (0.4 µM each, final concentration). The thermocycling conditions for amplification were 50°C (40 min) and 95°C (12 min), followed by 10 cycles of denaturation (95°C; 30 s), annealing (68°C; 20 s), and extension (72°C; 90 s), with the annealing temperature in these cycles reduced by 1°C each cycle. An additional 35 cycles of denaturation (95°C; 30 s), annealing (54°C; 20 s), and extension (72°C; 90 s) were performed, followed by a final extension (72°C; 7 min). Amplification products were separated in a 1.5% agarose, 1x Tris-acetate-EDTA gel and visualized by ethidium bromide staining and UV transillumination (33). For each step in RT-PCR amplification and analysis, strict protocols were followed to prevent cross-contamination of samples (e.g., use of dedicated pipettes for each step, use of aerosol-barrier pipette tips for all pipetting steps, changing gloves between steps, and performing all postamplification procedures in a room that was physically separate from that used for RT-PCR assembly). Negative control reactions (no template added) were demonstrated to be negative before subsequent analysis of amplification products.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide primers used for RT-PCR amplification and subsequent nucleotide sequencing
|
Alignments and phylogenetic analysis. Nucleotide and deduced amino acid sequences were aligned with CLUSTAL X (version 1.8) software. Sequences for comparison were obtained from GenBank and corresponded to nucleotide sequences of the CDV P, F, and H genes and deduced protein sequences. Phylogenetic analyses of nucleotide and amino acid alignments were performed using distance matrix methods (DNADIST or PRODIST, followed by NEIGHBOR) and maximum-parsimony methods (DNAPARS or PROTPARS) within the PHYLIP software package (13, 14). Data sets were subjected to bootstrap analysis, based on 100 resamplings of the original data set, using the SEQBOOT program within the PHYLIP software package (13, 14) to produce a majority rule consensus tree. Completed tree files were visualized using TreeView 1.5 (31).
Nucleotide sequence accession numbers. The CDV nucleotide sequences that were obtained in this study have been deposited in the National Center for Biotechnology Information database and assigned GenBank accession numbers AY964107 (case 18133; P gene), AY964108 (case 18133; F and H genes), AY964109 (case 19876; P gene), AY964110 (case 19876; F and H genes), AY964111 (case 21261; P gene), AY964112 (case 21261; F and H genes), AY964113 (case 25259; P gene), and AY964114 (case 25259; F and H genes).
|
|
|---|
RT-PCR and sequence analysis. Extracted RNA from necropsy tissue samples of all seven dogs from the four cases (which represented four independent premises) were positive for CDV, yielding an amplification product of the expected size (287 bp) using oligonucleotide primers previously described for RT-PCR detection of CDV infection in dogs (15). Samples of the same RNA from a single dog representing each of the four cases (and thus each of the four premises) were subsequently subjected to additional RT-PCR amplification using oligonucleotide primers (Table 1) designed to amplify the F and H genes, plus the intergenic space, and the central region of the P gene. Overlapping amplification products ranging in size from 562 to 1,229 bp were produced, using multiple combinations of the oligonucleotide primers shown in Table 1. The nucleotide sequence of each fragment was determined by direct sequencing of the RT-PCR products. Assembly of the nucleotide sequence data yielded full gene sequences for the F and H genes plus the intergenic region (3,968 bp) and partial sequence (540 bp) for the P gene from each of the four cases.
Alignment of the H gene nucleotide and predicted amino acid sequences demonstrated that the nucleotide sequences of the CDVs detected in cases 18133 and 21261 were 100% identical. In contrast, predicted H gene amino acid sequences from these cases were 93.4% and 95.6% identical to those from cases 19876 and 25259, respectively, while case 19876 was 94.6% identical to 25259 (Table 2). Following alignment, similar levels of nucleotide sequence identity were observed between the CDVs detected in these cases (Table 2). Alignment of sequences to CDV sequences from the public databases (Fig. 1) demonstrated identity levels ranging from 93.1% to 97.9% for the predicted amino acid sequences, with similar levels of identity for the nucleotide sequences (Table 2). A total of eight potential N-linked glycosylation sites were identified in the predicted H gene amino acid sequences with either seven (cases 19876 and 25259) or eight (cases 18133 and 21261) of these sites conserved (Fig. 1). Additionally 12 of 12 cysteine residues, 33 or 34 of 35 proline residues, and 33 or 34 of 35 glycine residues were conserved among the viruses detected in the cases in this study when compared by amino acid alignment to other CDV isolates (Fig. 1). Phylogenetic analysis of the H gene nucleotide sequences plus others from the databases using distance matrix methods (Fig. 2) demonstrated that cases 18133 and 21261 were placed into a distinct clade with phocine distemper virus 2 (PDV-2) and CDVs detected in dogs from China and Greenland. A bootstrap value of 100 for this clade suggests robust phylogenetic grouping. Phylogenetic analysis of the H gene nucleotide sequences of case 25259 placed the virus detected closest to that from a lesser panda, which was located in China. Phylogenetic analysis using maximum-parsimony analysis of H gene nucleotide sequences yielded equivalent results for CDV strains detected in cases 18133, 21261, and 25259 (data not shown). By distance matrix methods, analysis of the H gene sequence from case 19876 weakly grouped the virus detected from this case with that detected in a mink from Denmark (GenBank accession number Z47759) (Fig. 2). Maximum-parsimony analysis (data not shown) also implied a weak relationship between H gene amino acid sequences from case 19876 and the same CDV strain from a mink.
|
View this table: [in a new window] |
TABLE 2. Nucleotide and deduced amino acid sequence identities of CDV H protein genes
|
![]() View larger version (65K): [in a new window] |
FIG. 1. Alignment of deduced amino acid sequences from hemagglutinin genes of CDV strains. Only amino acids that differ from the majority sequence are shown. Identical residues are shown by dots. Potential N-linked glycosylation sites (N-X-S/T) are shaded. Cysteine ( ) and proline ( ) residues are indicated. The predicted amino acid sequences of case 21261 were identical to those shown for case 18133. Chn, China; Dnk, Denmark.
|
![]() View larger version (32K): [in a new window] |
FIG. 2. Phylogenetic tree for the complete H gene sequences of representative CDVs plus those detected in cases 18133, 21261, 19876, and 25259 (shown in boxes). The unrooted tree was generated using the distance matrix program NEIGHBOR with the jumble option invoked. Distance values were calculated by the DNADIST program within the PHYLIP software package using the Kimura two parameter. Only bootstrap values greater than 70 are shown, and the branch lengths are proportionate to genetic distances. The country of origin of each CDV is indicated by a two- or three-letter abbreviation following the isolate designation; Chn, China; Jpn, Japan; Ger, Germany; SA, South Africa; USA, United States of America; Dnk, Denmark; Trk, Turkey; Twn, Taiwan; Gnd, Greenland; Sib, Siberia. Canine distemper viruses used for comparison with the country of origin and year of isolation (if known), plus the accession numbers, are as follows: A75/17 (USA, 1975), AF164967; black leopard (USA, 1991), Z47763; black panther, strain A92-6 (USA, 1992), Z54166; Chinese leopard, strain A92-27/4 (USA, 1992), Z54156; dog (USA, 1989), Z47762; dog 4513 (Ger), Z77673; dog (Taiwan), AY378091; dog Hamamatsu (Japan, 1992 to 1994), D85754; dog Yanaka (Japan, 1992 to 1994), D85755; dog 5804 (Germany, 1990), AY386315; dog 91A (Denmark, 1991), AF478544; dog 91B (Denmark, 1991), AF478546; dog 2544 (Germany, 1995), Z77672; dog (Turkey), AY093674; dog (China), AF172411; dog GR88 (Greenland, 1988), Z47760; dog 26D (Japan, 1999), AB040766; dog 5B (Japan, 1999), AY297453; dog 98-002 (Japan, 1998), AB025270; dog HM3 (Japan, 1999), AB040767; dog KDK-1 strain (Japan, 1991), AB025271; dog Ueno (Japan, 1992 to 1994), D85753; ferret (German, 1989), X84999; giant panda (China), AF178038; javelina (USA, 1989), Z47764; lesser panda (China), AF178039; mink (Denmark, 1986), Z47759; Onderstpoort (South Africa), AF378705; PDV-2 Siberian seal (Siberia, 1988), X8499; raccoon 98-2646 (USA, 1988), AY542312 (also represents raccoon 98-2655 [USA, 1998], AY548109, and raccoon 98-2666 [USA, 1998], AY548110); raccoon 01-2689 (USA, 2001), AY465925 (also represents raccoon 01-2690 [USA, 2001]); raccoon 00-2601 (USA, 2000), AY443350; raccoon 01-2641(USA, 2001), AY526496; raccoon 01-2676 (USA, 2001), AY498692; raccoon 98-2654 (USA, 1998), AY466011; raccoon (USA, 1989), Z47765; raccoon dog (Japan, 1996), AB016776; Snyder Hill (Germany), AF259552.
|
|
View this table: [in a new window] |
TABLE 3. Nucleotide and deduced amino acid sequence identities of CDV F protein genes
|
![]() View larger version (16K): [in a new window] |
FIG. 3. Phylogenetic tree for the complete F gene sequences of representative CDVs plus those detected in cases 18133, 21261, 19876, and 25259 (shown in boxes). The unrooted tree was generated using the distance matrix program NEIGHBOR with the jumble option invoked. Distance values were calculated by the DNADIST program within the PHYLIP software package using the Kimura two parameter. Only bootstrap values greater than 70 are shown, and the branch lengths are proportionate to genetic distances. The country of origin of each CDV is indicated by a two- or three-letter abbreviation following the isolate designation: Ger, Germany; SA, South Africa; USA, United States of America; Sib, Siberia. Canine distemper viruses used for comparison with the country of origin and year of isolation (if known), plus the accession numbers, are as follows: A75/17 (USA, 1975), AF164967; dog 5804 (Germany, 1990), AY386315; dog 2544 (Germany, 1995), AJ007711; raccoon 00-2601 (USA, 2000), AY443350; raccoon 01-2689 (USA, 2001), AY649446; raccoon 98-2646 (USA, 1998), AY542312; raccoon 98-2654 (USA, 1998), AY466011; raccoon 98-2645 (USA, 1998), AY445077; Onderstepoort (South Africa), AF014953; PDV-2 (Siberia, 1988), L07075.
|
|
View this table: [in a new window] |
TABLE 4. Nucleotide and deduced amino acid sequence identities of CDV P gene fragments
|
![]() View larger version (19K): [in a new window] |
FIG. 4. Phylogenetic tree for partial P gene sequences of representative CDVs plus those detected in cases 18133, 21261, 19876, and 25259 (shown in boxes). The unrooted tree was generated using the distance matrix program NEIGHBOR with the jumble option invoked. Distance values were calculated by the DNADIST program within the PHYLIP software package using the Kimura two parameter. Only bootstrap values greater than 70 are shown, and the branch lengths are proportionate to genetic distances. The country of origin of each CDV is indicated by a two- or three-letter abbreviation following the isolate designation: Bul, Bulgaria; Jpn, Japan; Ger, Germany; SA, South Africa; USA, United States of America; Sib, Siberia; Swd, Sweden. Canine distemper viruses used for comparison with the country of origin and year of isolation (if known), plus the accession numbers, are as follows: A75/17 (USA, 1975), AF164967; dog Alaska (USA, 2003) (25); dog Hamamatsu (Japan, 1992 to 1994), AB028915; dog Jujo (Japan), AB028916; dog Yanaka (Japan, 1992 to 1994), AB028914; dog (Germany, 1993), AF259549; dog 5804 (Germany, 1990), AY386315; ferret (Germany, 1989), AF259550; PDV-2 (Siberia, 1988), AF259551; raccoon 98-2645 (USA, 1998), AY445077; raccoon 98-2646 (USA, 1998), AY542312; raccoon 98-2654 (USA, 1998), AY466011; raccoon 00-2601 (USA, 2000), AY443350; Rockborn (Sweden), AF181446.
|
|
|
|---|
Phylogenetic analyses of CDV P, F, and H gene nucleotide and amino acid sequences have been performed to study evolutionary relationships between CDV isolates and to find genetic variation among wild-type and CDV vaccine strains. The P gene is highly conserved among members of CDV strains and has been used for phylogenetic analysis by others (7, 9, 25, 32, 36). However, relatively few CDV P gene sequences of the length sequenced in this study were available in the public databases for phylogenetic comparison. The F gene is conserved within morbillivirus species (6, 24) and has also been used to determine phylogenetic relationships among these viruses (37). The H gene is more variable among CDV isolates, perhaps due to the role the protein plays in the host immune reaction (11), and thus has been widely used for phylogenetic analysis (8, 17, 18, 20, 28). Additionally, many more complete H gene sequences than complete F gene sequences are available in public databases for phylogenetic comparisons.
In the present study, phylogenetic analysis of the complete H gene sequences provided the most robust, definitive results of the three genes examined due to both the observed sequence divergence and the large number of gene sequences available for comparison relative to the F and P genes. Analysis of the H gene sequences implied an evolutionary relationship of cases 18133 and 21261 to PDV-2, as well as H gene sequences of CDV strains detected in dogs from China and Greenland. The relationship of 18133 and 21261 to PDV-2 was further supported by phylogenetic analysis of complete F gene sequences and partial P gene sequences. Analysis of H gene sequences also suggested that case 25259 may be related to CDV isolates previously detected in a lesser panda from China. Unfortunately, sequence from this panda strain was not available for F or P gene sequence analysis, and thus, the relationship could not be further supported by phylogenetic analyses. Analysis of H gene sequences by both distance and maximum-parsimony methods only weakly implied a relationship for CDV detected in case 19876 to a CDV strain detected in mink from Denmark, suggesting that the strain detected in case 19876 was the most divergent of the those studied.
Given the relationships suggested by phylogenetic analysis, it appears possible that the CDV strains infecting dogs in this study may have originated from noncanine species or, alternatively, may have been transmitted from dogs to the other species. Interspecies transmission has been widely suspected for CDV (10, 16, 19). For example, CDV isolated from captive large felids in the United States was thought to have originated from feral, nonfelid carnivores (18). In a study examining free-ranging carnivores from the Serengeti ecosystem of Africa, closely related CDV isolates were detected from four different species (9), suggesting extensive interspecies transmission. Analysis of PDV-2, which was isolated from Siberian seals, demonstrated a relationship to CDV isolates from a German dog and ferret (26), suggesting that transmission from one of these hosts may have been responsible for the widespread infection seen in these freshwater seals (27).
Recent studies have characterized CDV strains in both dogs and noncanine species in North America. The phylogenetic relationship of CDV strains isolated from raccoons concluded that these viruses were distinct from, yet strongly related to, known American-lineage CDV isolates (23). A CDV outbreak that killed several hundred sled dogs was recently described in Alaska. Sequence analysis of 540 bp of the P gene suggested a close relationship of the strain in the sled dogs to a strain of distemper virus isolated from a Siberian seal (PDV-2) (25). Pairwise alignment of CDV partial P gene sequences from the cases in the present study also demonstrated that the P gene sequences, as well as the H and F gene sequences, of 18133 and 21261 were potentially related to PDV-2 and thus may also be related to the virus detected in the Alaskan dogs in 2003. Unfortunately, only partial sequence of the P gene (but no H or F gene sequences) was available from the CDV outbreak in the Alaskan dogs, and thus, a more definitive link between these cases could not be further explored.
In dog populations with a high rate of vaccination against canine diseases, clinical cases of CDV typically occur sporadically. Speculations put forth to explain these cases are varied but often include vaccine failures, reversion of attenuated CDV vaccine strains to virulence (5), or the emergence of new strains that are sufficiently divergent to evade immune protection elicited by the vaccines used. In the cases studied here, vaccination with modified live vaccine was performed recently in three of the four cases. Importantly, the CDV strains detected in these clinical cases were clearly distinct from known vaccine strains, as has been demonstrated in previous studies of other CDV field strains (8, 17, 21, 23, 25, 28). The CDV strain Lederle, which was used as a vaccine strain in the past (19) and may still be circulating in wildlife, is phylogenetically related to the CDV strain Snyder Hill (23) and thus genetically distant from the strains detected in the present study. The widely used Galaxy-D CDV vaccine (Schering-Plough, Kenilworth, NJ) originates from the Onderstepoort strain, which is also genetically distant from the CDV strains detected in each of the four cases described here. To further explore the possibility of vaccine virus reverting to virulence, partial sequencing (979 bases; data not shown) of the H gene was performed from the same production lot of the CDV vaccine (Vanguard 5; Pfizer, Inc., Exton, PA) used in dogs from case 18133. The nucleotide sequence detected from this vaccine was 97.7% identical to that of strain A75/17 (34), which was the closest match in GenBank over the region sequenced, but only 95.0% identical to the virus detected in cases 18133 and 21261. Together, these observations suggest that a recent reversion of vaccine virus to virulence was not likely to be the cause of the clinical CDV disease seen in the dogs for which vaccine had been recently used. However, the presence of a minor though virulent strain as a contaminant of the vaccine virus stock certainly cannot be excluded.
In summary, phylogenetic analysis of H and F genes showed that the CDV strains detected in this study were genetically distinct from viruses previously detected in the continental United States and most closely matched strains from Asia or possibly Europe. In three of the four cases examined, a recent history of CDV vaccination was reported, yet none of the viruses detected were related to known vaccine strains. These data suggest the possible emergence of novel CDV strains in North America, which may require renewed efforts by vaccine manufacturers to ensure adequate protection following immunization of dogs against CDV.
We thank Susan K. Schommer and William J. Mitchell, Jr., for insightful discussion and reading of the manuscript.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»