|
|
||||||||
1 Southeastern Cooperative Wildlife Disease Study, Department of Population Health, College of Veterinary Medicine, University of Georgia, Athens, Georgia, 30602, USA
2 Department of Pathology, College of Veterinary Medicine, University of Georgia, Athens, Georgia 30602, USA
3 Corresponding author (email: mmurphy{at}vet.uga.edu)
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Prior studies of the molecular epidemiology of EHDV-2 have shown that viruses within this serotype are evolving very slowly and likely exist as a contiguous population throughout the southeastern and western US (Murphy, 2003). Although EHDV-1 was first isolated from New Jersey white-tailed deer in 1955 (Shope et al., 1960), because of the paucity of isolates, no comparable studies have been performed with this serotype. During the autumn of 1999 (mid-Augustlate September), a large outbreak of hemorrhagic disease in white-tailed deer caused by EHDV-1 occurred in a wide band in the eastern US, extending from Georgia to New Jersey (Fig. 1
). This outbreak not only afforded a unique opportunity to explore the genetic relatedness among EHDV-1 isolates within a regional outbreak, but also enabled us to compare recent and older isolates to assess long-term genetic change.
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
All viruses (Table 1
) were isolated in cow pulmonary artery endothelial (CPAE) or baby hamster kidney (BHK21) cells as previously described (Murphy et al., 2005) and passaged fewer than four times. Viruses from Texas were isolated from the blood of clinically normal white-tailed deer fawns; all other viruses were isolated from moribund and deceased white-tailed deer or moribund mule deer (Odocoileus hemionus). Infected BHK21 cells in 75 cm2 tissue-culture flasks were used for total RNA extraction using RNAzolB (Tel-Test Incorporated, Friendswood, Texas, USA), according to the manufacturers protocol.
|
|
Gel-purified S10 and L2 gene PCR products (729 base pairs and 594 base pairs in length, respectively) were sequenced on the forward and reverse strands at the Molecular Genetics Instrumentation Facility, University of Georgia (Athens, Georgia, USA) using dye-deoxy terminator nucleotides (Big Dye, Applied Biosystems, Foster City, California, USA) and an ABI 3100 Capillary Sequencer (Applied Biosystems). A panel of internal and flanking primers (Table 2
) was used in multiple sequencing reactions to generate overlapping contiguous fragments and a subsequent consensus sequence, corresponding to bases 21765 (S10) and 5311124 (L2) of the published sequences (S10, GenBank accession number L29023 [Jensen and Wilson, 1995]; L2, GenBank accession number D10767 [Iwata et al., 1992]). Nucleotide sequences determined in this study are available in GenBank, accession numbers DQ899833-DQ899874.
Neighbor-joining sequence analysis
Chromatograms were viewed and contiguous fragments (contigs) built using the Se-quencher program (Gene Codes Corporation, Ann Arbor, Michigan, USA). Contigs were formatted with the ToPir program of the Wisconsin package (GCG, version 8.0 software, Accelrys, San Diego, California, USA) and aligned with the Clustal X program (Thompson et al., 1994) using default parameters. The Modeltest program, version 3.06 (Posada and Crandall, 1998), was utilized to estimate models of nucleotide substitution for each locus. Pairwise distances between isolates were calculated using a maximum likelihood approach as implemented in the PAUP* program, version 4.0, beta 10 (Swofford, 2002), to correct for multiple substitutions. The Neighbor-Joining search algorithm in PAUP* was utilized for generating a phylogram and bootstrapped consensus tree (100 replicates) for each gene segment. Trees were outgroup rooted using the cognate gene region of BTV-13 (S10 gene, Genbank accession number AFO44712 [Pierce et al., 1998]) and the Ibaraki (Japanese) strain of EHDV-2 (L2 gene, GenBank accession number ABO30735 [Ohashi et al., 2002]). The overall number of nucleotide differences among isolates at each locus was calculated using the MEGA program (Kumar et al., 2001). MEGA was also utilized to translate the open reading frame of each isolate, as well as compare the number of amino acid substitutions between viruses.
| RESULTS |
|---|
|
|
|---|
For the analyses of the S10 locus dataset, a maximum likelihood model of nucleotide substitution was selected, and model parameters estimated utilizing a hierarchical likelihood ratio test method as implemented in the Modeltest 3.06 and PAUP* software programs. The outgroup taxon (BTV-13) was removed from the model selection and parameter estimation process because its inclusion in the estimation decreased the log likelihood score of the best-fit model selected for the ingroup and thus is assumed to evolve via a different mechanism from the ingroup taxa. The estimated parameters include a transition-to-transversion ratio of 20.49, and frequencies of each nucleotide (A = 0.3231, C = 0.1866, G = 0.2602, and T = 0.2301). The rate of substitution among sites was estimated to be equal, and no sites were considered invariant.
In the S10 phylogram (Fig. 2
), isolates segregate into two distinct clades. The eastern clade contains the 1955 prototype EHDV-1 strain, originally isolated from a New Jersey white-tailed deer in 1955, as well as a relatively conserved (maximum number of nucleotide changes between samples = 7) grouping of isolates collected in the eastern US between 1991 and 1999. The western clade contains isolates collected in the western US between 2000 and 2001 and has a maximum of 10 nucleotide substitutions between isolates. Within the eastern clade, there is further division of isolates. The largest subclade contains the majority of the 1999 isolates, most of which are identical. Also within the eastern clade are three smaller clades, containing a 1996 TN isolate, a 1991 TN isolate, and two 1999 Georgia isolates. None of the isolates in the S10 western clade are identical, even though four of the isolates (TX(A)(D)) are from the same herd, county, and year.
|
The L2 phylogram (Fig. 3
) displays nearly identical region-specific groupings as found in the S10 phylogram. Several of the 1999 isolates are identical at this locus, as are two of the western isolates. In both phylograms the branch lengths among the western isolates are more heterogeneous than the branch lengths connecting the eastern isolates. This indicates eastern and western regionspecific differences in the level of genetic variability.
|
| DISCUSSION |
|---|
|
|
|---|
As with the 1999 EHDV-1 viruses, a high degree of genetic relatedness among EHDV-2 isolates from white-tailed deer in the US has also been reported (Murphy et al., 2005). With EHDV-2, there is an apparent lack of a regional segregation (topotyping) in phylograms based on both the S10 and L2 gene sequences (Murphy, 2003). White-tailed deer in the southeastern US are exposed to EHDV-2 in two-to-three-year cycles (Stallknecht et al., 2002), and it is possible that these frequent HD events in a nearly contiguous deer host population would create a situation in which EHDV-2 could "travel" between the eastern and western US over a virtual animal "bridge." This may result in the formation of a single virus population.
In contrast, EHDV-1 appears to exist in distinct regional populations representing the eastern and western US, as evidenced by the topotyping effect seen in both the S10 and L2 phylograms. Unlike EHDV-2, outbreaks of HD caused by EHDV-1 have been described only twice, in 1955 (Shope et al., 1960) and 1975 (McConnell et al., 1977), both occurring in New Jersey. However, serologic and virus isolation evidence indicates that white-tailed deer in the southeastern US are occasionally exposed to this serotype, especially in the coastal plain areas, and in Texas, where annual enzootic activity involving multiple EHDV and BTV occurs (Stallknecht et al., 1996, 2002). It is possible that the regional variation detected among the EHDV-1 isolates is related to the lack of frequent epizootic activity involving this serotype; this is in direct contrast to the regular cycles of transmission observed with EHDV-2. Therefore, unlike EHDV-2, EHDV-1 may be more restricted to localized enzootic areas and in the absence of frequent epizootic events is unable to traverse over long geographic spans.
The phylograms of both the S10 and L2 loci also demonstrate a variation in branch length among the eastern versus western EHDV-1 isolates. This suggests a different level of nucleotide sequence heterogeneity among the viruses from these two regions. Nearly all of the 1999 isolates were identical at both loci, even though they originated from different geographic locations at different time points. The paucity of nucleotide sequence variation during an epizootic is not surprising, in that peak viremia (and therefore the peak transmission period) in the white-tailed deer host likely would occur several days before the emergence of an antibody response that could exert selective pressure upon these viruses (Quist et al., 1997; Gaydos et al., 2002a). This high level of genetic conservation also has been described among isolates collected from diseased deer during epizootics of EHDV-2related HD in the eastern and western US (Murphy, 2003; Murphy et al., 2005). The observation of a higher level of genetic variation among viruses cocirculating in an enzootic transmission pattern (Texas isolates) in contrast to the genetic conservation observed among isolates from an epizootic transmission pattern (1999 eastern isolates) is not unique to EHDV-1. A similar phenomenon has been described for other arboviruses including Venezuelan equine encephalitis (family Togaviridae, genus Alphavirus), wherein the substitution rate among isolates associated with epizootics up to 35 years apart was 20-fold lower than that demonstrated for a closely related enzootic subtype (Powers et al., 1997).
In contrast to the conservation observed among viruses in the eastern subclades, the branches connecting the isolates in the western subclades of both the S10 and L2 phylograms are of variable length. This suggests a higher level of genetic variation among the western isolates. Although based on only a small number of western samples (n = 6), this finding may be related to the broader geographic area represented by the western isolates, and the fact that these isolates represent areas in which both enzootic (Texas) and epizootic (Ida-ho, Wyoming) transmission occur. All of the 2000 Texas isolates originated from white-tailed deer at the same research facility in a hyperenzootic region (Stallknecht et al., 1996), and it is evident that several strains of EHDV-1 were cocirculating among these animals. Additionally, the 2000 Texas isolates were made from clinically normal fawns, and three of these fawns were viremic in the presence of serum-neutralizing antibodies to this serotype (Gaydos et al., 2002b). The simultaneous presence of virus and antibodies could provide an opportunity for antibody selection and continuous generation of viral variants in white-tailed deer from enzootic foci. It is interesting that the two 1999 EHDV-1 isolates from Georgia (1999 GA(B) and 1999 GA(C)), which did not group with the other 1999 isolates in either the S10 or L2 phylograms, also originated from a location immediately adjacent to the coastal plain, an area of recognized enzootic activity (Stallknecht et al., 1995).
The 1999 outbreak provides an example of epizootic transmission, with rapid dissemination, host mortality, and very little viral genetic variation. The genetic similarity between the 1999 eastern viruses and other eastern isolates from previous years is consistent with the presence of known enzootic foci of EHDV-1 as exist in the coastal plain of the eastern US (Stallknecht et al., 1995). Viruses within these geographic foci may be genetically distinct from those in western endemic foci, as is present in Texas (Stallknecht et al., 1996), since a strong topotyping effect was seen among eastern and western EHDV-1 isolates. Although a topotype effect also was observed among EHDV-2 viruses before 1990, it was not demonstrable among more recent (19902001) isolates (Murphy, 2003). This difference in phylogenetic pattern demonstrates the need for additional studies to further delineate the molecular epidemiology of EHDV-1.
| ACKNOWLEDGMENTS |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
CHENEY, I. W., M. D. LARSON, J. O. MECHAM, AND W. C. WILSON. 1995. Geographical genetic variation in the gene encoding VP3 from the Alberta isolate of epizootic hemorrhagic disease virus. Virus Research 36: 279286.[Medline]
, M. YAMAKAWA, P. ROY, J. O. MECHAM, AND W. C. WILSON. 1996. Molecular characterization of the segment 2 gene of epizootic hemorrhagic disease virus serotype 2: Gene sequence and genetic diversity. Virology 224: 555560.[Medline]
DEMAULA, C. D., K. R. BONNEAU, AND N. J. MACLACHLAN. 2000. Changes in the outer capsid proteins of bluetongue virus serotype ten that abrogate neutralization by monoclonal antibodies. Virus Research 67: 5966.[Medline]
GAYDOS, J. K., W. R. DAVIDSON, F. ELVINGER, D. G. MEAD, E. W. HOWERTH, AND D. E. STALLKNECHT. 2002a. Innate resistance to epizootic hemorrhagic disease in white-tailed deer. Journal of Wildlife Diseases 38: 713719.[Abstract]
, D. E. STALLKNECHT, D. KAVANAUGH, R. J. OLSON, AND E. FUCHS. 2002b. Dynamics of maternal antibodies to hemorrhagic disease viruses (Reoviridae: Orbivirus) in white-tailed deer. Journal of Wildlife Diseases 38: 253257.[Abstract]
GOULD, A. R., AND B. T. EATON. 1990. The amino acid sequence of the outer coat protein VP2 of neutralizing monoclonal antibody-resistant, virulent and attenuated bluetongue viruses. Virus Research 17: 161172.[Medline]
HOLBROOK, F. R. 1996. Biting midges and the agents they transmit. In The biology of disease vectors, B. J. Beaty and W. C. Marquardt (eds.). University Press of Colorado, Niwot, Colorado, pp. 110116.
HUISMANS, H., AND B. J. ERASMUS. 1981. Identification of the serotype-specific and group-specific antigens of bluetongue virus. Onderstepoort Journal of Veterinary Research 48: 5158.[Medline]
IWATA, H., T. CHUMA, AND P. ROY. 1992. Characterization of the genes encoding two of the major capsid proteins of epizootic hemorrhagic disease virus indicates a close genetic relationship to bluetongue virus. Journal of General Virology 73: 915924.
JENSEN, M. J., AND W. C. WILSON. 1995. A model for the membrane topology of the NS3 protein as predicted from the sequence of segment 10 of epizootic hemorrhagic disease virus serotype 1. Archives of Virology 140: 799805.[Medline]
KNUDSON, D. L., AND T. P. MONATH. 1990. Orbiviruses. In Fields virology, 2nd Edition, B. N. Fields, D. M. Knipe, R. M. Chanock, M. S. Hirsch, J. L. Melnick, T. P. Monath, and B. Roizman (eds.). Raven Press, New York, New York, pp. 14051433.
KUMAR, S., K. TAMURA, I. B. JAKOBSEN, AND M. NEI. 2001. MEGA2: Molecular evolutionary genetics analysis software. Arizona State University, Tempe, Arizona, www.megasoftware.net.
MCCONNELL, P. A., R. C. LUND, AND N. R. BOSS. 1977. The 1975 outbreak of hemorrhagic disease among white-tailed deer in northwestern New Jersey. Transactions of Northeast Section of the Wildlife Society 33: 3544.
MURPHY, M. D. 2003. Molecular epidemiology of epizootic hemorrhagic disease viruses in the United States. PhD Thesis, University of Georgia, Athens, Georgia, 115 pp.
, E. W. HOWERTH, N. J. MACLACHLAN, AND D. E. STALLKNECHT. 2005. Genetic variation among epizootic hemorrhagic disease viruses in the southeastern United States: 19782001. Infection Genetics and Evolution 5: 157165.
NETTLES, V. F., W. R. DAVIDSON, AND D. E. STALLKNECHT. 1992. Surveillance for hemorrhagic disease in white-tailed deer and other wild ruminants, 19801989. In Proceedings of the annual conference of the Southeastern Association of Fish and Wildlife Agencies 46: 138146.
OHASHI, S., K. YOSHIDA, T. YANASE, AND T. TSUDA. 2002. Analysis of intratypic variation evident in an Ibaraki virus strain and its epizootic hemorrhagic disease virus serogroup. Journal of Clinical Microbiology 40: 36843688.
PEARSON, J. E., G. A. GUSTAFSON, A. L. SHAFER, AND A. D. ALSTAD. 1992. Distribution of bluetongue in the United States. In Bluetongue, African horsesickness and related orbiviruses, T. E. Walton and B. I. Osburn (eds.). CRC Press, Boca Raton, Florida, pp. 128139.
PIERCE, C. M., U. B. R. BALASURIYA, AND N. J. MACLACHLAN. 1998. Phylogenetic analysis of the S10 gene of field and laboratory strains of bluetongue virus from the United States. Virus Research 55: 1527.[Medline]
POSADA, D., AND K. A. CRANDALL. 1998. Modeltest: Testing the model of DNA substitution. Bioinformatics 14: 817818.
POWERS, A. M., M. S. OBERSTE, A. C. BRAULT, R. RICO-HESSE, S. M. SCHMURA, J. F. SMITH, W. KANG, W. P. SWEENEY, AND S. C. WEAVER. 1997. Repeated emergence of epidemic/epizootic Venezuelan equine encephalitis from a single genotype of enzootic subtype ID virus. Journal of Virology 71: 66976705.[Abstract]
QUIST, C. F., E. W. HOWERTH, D. E. STALLKNECHT, J. BROWN, T. PISELL, AND V. F. NETTLES. 1997. Host defense responses associated with experimental hemorrhagic disease in white-tailed deer. Journal of Wildlife Diseases 33: 584599.[Abstract]
ROY, P. 1996. Orbivirus structure and assembly. Virology 216: 111.[Medline]
SHOPE, R. E., L. G. MACNAMARA, AND R. MANGOLD. 1960. A virus-induced epizootic hemorrhagic disease of Virginia white-tailed deer (Odocoileus virginianus). Journal of Experimental Medicine 111: 155170.[Abstract]
STALLKNECHT, D. E., E. W. HOWERTH, AND J. K. GAYDOS. 2002. Hemorrhagic disease in white-tailed deer: Our current understanding of risk. Transactions of the North American Wildlife and Natural Resources Conference 67: 7586.
, M. P. LUTTRELL, K. E. SMITH, AND V. F. NETTLES. 1996. Hemorrhagic disease in white-tailed deer in Texas: A case for enzootic stability. Journal of Wildlife Diseases 32: 695700.[Abstract]
, V. F. NETTLES, E. A. ROLLOR III, AND E. W. HOWERTH. 1995. Epizootic hemorrhagic disease and bluetongue virus serotype distribution in white-tailed deer in Georgia. Journal of Wildlife Diseases 31: 331338.[Abstract]
SWOFFORD, D. L. 2002. PAUP*: Phylogenetic analysis using parsimony (*and other methods). Version 4.0 beta 10. Sinauer Associates, Sunderland, Massachusetts,
THOMPSON, J. D., D. G. HIGGINS, AND T. J. GIBSON. 1994. Clustal W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position specific gap penalties and weight matrix choice. Nucleic Acids Research 22: 46734680.
Received for publication 21 July 2003.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |