J Gen Virol Email Content Delivery
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 1509-1519; DOI 10.1099/vir.0.81757-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary tables
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Dewals, B.
Right arrow Articles by Vanderplasschen, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Dewals, B.
Right arrow Articles by Vanderplasschen, A.
Agricola
Right arrow Articles by Dewals, B.
Right arrow Articles by Vanderplasschen, A.
© 2006 Society for General Microbiology

Evolution of Bovine herpesvirus 4: recombination and transmission between African buffalo and cattle

Benjamin Dewals1, Muriel Thirion1, Nicolas Markine-Goriaynoff1, Laurent Gillet1, Katalin de Fays1, Frédéric Minner1, Virginie Daix1, Paul M. Sharp2 and Alain Vanderplasschen1

1 Immunology–Vaccinology (B43b), Department of Infectious and Parasitic Diseases, Faculty of Veterinary Medicine, University of Liège, B-4000 Liège, Belgium
2 Institute of Genetics, University of Nottingham, Queens Medical Centre, Nottingham NG7 2UH, UK

Correspondence
Alain Vanderplasschen
A.vdplasschen{at}ulg.ac.be


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Bovine herpesvirus 4 (BoHV-4) has been isolated from cattle throughout the world, but virological and serological studies have suggested that the African buffalo is also a natural host for this virus. It has previously been found that the Bo17 gene of BoHV-4 was acquired from an ancestor of the African buffalo, probably around 1.5 million years ago. Analysis of the variation of the Bo17 gene sequence among BoHV-4 strains suggested a relatively ancient transmission of BoHV-4 from the buffalo to the Bos primigenius lineage, followed by a host-dependent split between zebu and taurine BoHV-4 strains. In the present study, the evolutionary history of BoHV-4 was investigated by analysis of five gene sequences from each of nine strains representative of the viral species: three isolated from African buffalo in Kenya and six from cattle from Europe, North America and India. No two gene sequences had the same evolutionary tree, indicating that recombination has occurred between divergent lineages; six recombination events were delineated for these sequences. Nevertheless, exchange has been infrequent enough that a clonal evolutionary history of the strains could be discerned, upon which the recombination events were superimposed. The dates of divergence among BoHV-4 lineages were estimated from synonymous nucleotide-substitution rates. The inferred evolutionary history suggests that African buffalo were the original natural reservoir of BoHV-4 and that there have been at least three independent transmissions from buffalo to cattle, probably via intermediate hosts and – at least in the case of North American strains – within the last 500 years.

The GenBank/EMBL/DDBJ accession numbers for the nucleotide sequences determined in this work are AY847305 [GenBank] –AY847311 [GenBank] and AY847322 [GenBank] –AY847348 [GenBank] .

Supplementary tables with primer details and GenBank accession numbers for the BoHV-4 sequenced regions are available in JGV Online.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Bovine herpesvirus 4 (BoHV-4) is a gammaherpesvirus that has been isolated throughout the world from healthy cattle, as well as those exhibiting a variety of diseases (Li et al., 2005Down; Markine-Goriaynoff et al., 2003bDown; Thiry et al., 1992Down). The widespread distribution of BoHV-4 in cattle (Bos taurus and Bos indicus) populations justified the conclusion that the bovine species is the natural host; hence the nomenclature of the virus. However, BoHV-4 has also been detected in a variety of other ruminants (Todd & Storz, 1983Down; Van Opdenbosch et al., 1986Down). In particular, BoHV-4 infection is very common among wild African buffalo (Syncerus caffer). An early survey performed mainly within one reserve in Kenya revealed that BoHV-4 can be isolated from the blood of healthy wild African buffaloes with a high frequency (25 %) and that almost all (94 %) animals tested had antibodies against BoHV-4 (Rossiter et al., 1989Down). More recently, our seroprevalence study of 400 sera from wild African buffaloes from numerous locations in eastern and southern Africa revealed that, independent of their geographical origin, members of this species exhibit a seroprevalence of anti-BoHV-4 antibodies of around 70 % (Dewals et al., 2005Down), a rate much higher than is generally observed in cattle populations (16–18 %) (Wellenberg et al., 1999Down). These observations suggest that African buffalo may be the original host species of BoHV-4 and that other species, including domestic cattle, may have acquired the virus more recently.

An independent observation supports the role of the African buffalo as the original host species of BoHV-4. The Bo17 gene of BoHV-4 is the only viral gene known to date that encodes a homologue of the cellular core 2 beta-1,6-N-acetylglucosaminyltransferase-mucin type (Vanderplasschen et al., 2000Down). Our phylogenetic study demonstrated that the Bo17 gene is related significantly more closely to the cellular homologue from African buffalo than to that from cattle, indicating that the viral gene was acquired from a recent ancestor of the African buffalo (Markine-Goriaynoff et al., 2003aDown). Furthermore, the gene-acquisition event was dated at around 1.5 million years ago, long after the separation of the Bos and Syncerus lineages and long before the domestication of cattle. The sequences of the Bo17 gene from nine BoHV-4 strains, three from African buffalo and six from cattle, were split into two clades according to the species of origin (Markine-Goriaynoff et al., 2003aDown), consistent with restriction-profile analyses showing a clear differentiation between BoHV-4 strains isolated from the two species (Bublot et al., 1991Down; Markine-Goriaynoff et al., 2003aDown; Rossiter et al., 1989Down). For Bo17, the divergence between the two viral lineages was estimated to have occurred about 0.7 million years ago, suggesting a relatively ancient transmission of BoHV-4 from the ancestor of the African buffalo to Bos primigenius, the wild ancestor of domesticated cattle. Among the cattle strains, the one isolated from zebu (B. indicus) formed an outgroup to those from the taurine cattle (B. taurus) lineage and the split was estimated to have occurred at about the same time as the divergence of the ancestors of B. taurus and B. indicus, around 0.2–0.3 million years ago, consistent with a host-dependent split between zebu and taurine BoHV-4 strains. Thus, whilst these data indicated that BoHV-4 has been co-evolving with the African buffalo for at least the last 1.5 million years, they also pointed to infection of cattle since long before they were domesticated, around 10 000 years ago.

These interpretations of the evolutionary history of BoHV-4 must be treated with caution because they rely on the analysis of a single gene, Bo17. Recombination has been reported in numerous species of herpesviruses (Burrows et al., 1996Down; Norberg et al., 2004Down; Poole et al., 1999Down; Thiry et al., 2005Down; Walling & Raab-Traub, 1994Down). The consequence of recombination is that different genes may have different evolutionary histories. Therefore, we set out to investigate the phylogenetic relationships among the nine BoHV-4 strains that were used in our former study (Markine-Goriaynoff et al., 2003aDown) by characterizing four more gene sequences from three different regions distributed across the genome. The results demonstrate that BoHV-4 evolution has indeed involved recombination among divergent strains and that the evolutionary history of the Bo17 gene is not typical of that of the genome as a whole. These data necessitate a new interpretation of the origins of the various lineages of BoHV-4 distributed worldwide among cattle.


   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Virus strains.
Nine strains representative of the BoHV-4 species were used (Table 1Down). The complete sequence of strain 66-p-347 has been published previously (Zimmermann et al., 2001Down). Sequences from the other eight strains were obtained in this study.


View this table:
[in this window]
[in a new window]
 
Table 1. BoHV-4 strains used in this study

 
Amplification and sequencing.
Three regions of the BoHV-4 genome – open reading frame (ORF) 16 (encoding a v-Bcl-2 homologue), Bo10 (encoding a putative glycoprotein unique to BoHV-4) and ORF 71 (encoding a v-FLIP homologue)–ORF 73 (encoding a LANA homologue) (ORF 71 and ORF 73 plus the intervening intergenic region) – were amplified by using PCR. The positions of these regions within the genome are shown in Fig. 1Down and the primers used for the amplifications are given in Supplementary Table S1 (available in JGV Online). PCR products were subsequently purified (Wizard PCR Preps DNA purification system; Promega) and sequenced by primer walking using a BigDye Terminator v3.1 cycle sequencing kit (Applied Biosystems). An ABI Prism 3730 DNA Analyser (Applied Biosystems) was used for analysis of samples.


Figure 1
View larger version (7K):
[in this window]
[in a new window]
 
Fig. 1. Regions of the BoHV-4 genome analysed in this study. The structure of the BoHV-4 genome is shown, with thepositions and orientations of the five genes within the long unique region (LUR) and the lengths of the sequences in strain 66-p-347.

 
Nucleotide sequence accession numbers.
The sequences reported in this paper have been deposited in GenBank. The accession numbers of these sequences, together with those of previously published sequences included in the analysis, are listed in Supplementary Table S2 (available in JGV Online).

Sequence analysis.
Sequences were aligned by using CLUSTAL W (Thompson et al., 1994Down), with minor manual adjustment. The numbers of non-synonymous (KA) and synonymous (KS) substitutions per site were estimated by the method of Li (1993)Down. Phylogenetic trees were estimated from nucleotide sequences by the neighbour-joining (NJ) method (Saitou & Nei, 1987Down), using distances corrected by Kimura's two-parameter method and 1000 bootstrap replicates, as implemented in CLUSTAL W. Phylogenetic trees were also estimated by the maximum-likelihood (ML) method, implemented using DNAML from the PHYLIP package (J. Felsenstein; http://evolution.genetics.washington.edu/phylip.html). No significant differences were found between the NJ and ML trees; the results from the NJ analyses are presented. For most of the sequences analysed, there was no suitable outgroup available and, so, each tree was rooted at its midpoint; for the one case where an outgroup was available (the Bo17 gene), midpoint rooting yielded the same result as outgroup rooting (Markine-Goriaynoff et al., 2003aDown). Observed substitution rates were consistent with midpoint rooting for the other genes. Recombination was inferred from topological differences between the trees for different genes and was not dependent on the position of the root (except in one case, discussed specifically). Mosaicism of sequences due to recombination was investigated by analysis of phylogenetically informative sites; the statistical significance of putative recombination breakpoints was assessed by a permutation test (Hatwell & Sharp, 2000Down).

The protein encoded by the Bo10 gene was found to contain a highly variable, Ser-rich region that was difficult to align with confidence. This region, corresponding to codons 38–94 in the 66-p-347 Bo10 sequence, was excluded from the comparative analyses.


   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
To investigate the diversity and evolution of BoHV-4 strains isolated from African buffalo and cattle, nine strains representative of the BoHV-4 species (Table 1Up) were compared for four different regions distributed across the genome (Fig. 1Up). Three strains were isolated from buffalo in Kenya, five from taurine cattle from Europe or North America and one from zebu cattle in India; these nine strains are representative of the three groups of BoHV-4 strains defined on the basis of restriction profiles. From left to right, the genomic regions selected for this study were ORF 16, the Bo10 ORF, ORF 71–ORF 73 and the previously characterized Bo17 ORF. The ORF 71–ORF 73 region includes two genes and the region between them (IG 71–73). A putative 74 codon ORF, called Bo13, overlapping ORF 71 and ORF 73 and thus spanning the IG 71–73 region, has been described in the genome of strain 66-p-347 (Zimmermann et al., 2001Down). However, we found an insertion of 1 nt within IG 71–73, resulting in a frameshift, in six of the nine viral strains. Also, in a comparison of strains 66-p-347 and M40 (the latter lacks the frameshift), there was an excess of non-synonymous over synonymous substitutions. These observations suggest that this ORF does not function as a gene.

In each region, the three European strains (V. test, MOVAR and LVR 140) were very similar, differing by no more than 2 nt and showing overall diversity (summed across regions) of <0.1 %, consistent with the similarity of the EcoRI restriction profiles of the genomes of these three strains (Markine-Goriaynoff et al., 2003aDown). Similarly, the two North American strains (DN599 and 66-p-347) differed by only 1 nt in each of two regions, again consistent with their restriction profiles. In contrast, the three buffalo strains (108, 130 and Buf.) differed from each other by up to 2.5 % overall. The most divergent strain was M40 from zebu, differing on average at 4.7 % of sites. Phylogenetic analysis of a concatenated alignment of all four regions (a total of 4093 sites compared) yielded a tree in which the European and North American strains from taurine cattle clustered, the three African buffalo strains clustered and the Indian zebu strain appeared to be the outgroup (Fig. 2aDown). This tree differed from that found previously from analysis of the Bo17 gene alone (Markine-Goriaynoff et al., 2003aDown), where the M40 strain grouped quite close to the taurine cattle virus clade. Therefore, phylogenetic analyses of each of the individual regions were compared.


Figure 2
View larger version (24K):
[in this window]
[in a new window]
 
Fig. 2. Phylogenetic analyses of BoHV-4 strains. The trees were derived from (a) a concatenated alignment of all sequences or (b) separate regions of the BoHV-4 genome. Values on internal branches indicate the percentage of bootstrap replicates in which the branch was found; only values >60 % are shown. Asterisks highlight recombination events described in the text. Horizontal branches are drawn to scale. Bar, 0.01 nucleotide substitution per site.

 
BoHV-4 strains isolated from African buffalo and from cattle have undergone inter-strain recombination events since their last common ancestor
Phylogenetic analyses of six genomic regions (the five genes plus the intergenic sequence) yielded six different trees (Fig. 2bUp). Significant topological differences among the trees must reflect recombination during the divergence of the BoHV-4 strains. In comparison to the tree from the concatenated sequences, only the IG 71–73 tree was similar, but each of the gene trees differed with respect to the position of only one or a small number of strains. This suggests that the tree from the concatenated sequences reflects the clonal evolution of the viral strains, whereas the trees from individual sequences indicate the impact of occasional inter-strain recombination events. To facilitate interpretation of the differences among the trees, the extent of synonymous nucleotide substitution (KS) between pairs of strains was calculated for each gene. Synonymous-substitution rates are expected to be similar for different genes and, so, differences among genes in KS values for the same pair of strains may reflect differences in the time since common ancestry, due to recombination.

The ORF 16 tree differed from the concatenated tree in that the topological positions of the 130 and Buf. strains switched (Fig. 2bUp). The extent of synonymous substitution between Buf. and 108 was lower than for any other gene (Table 2Down), suggesting that the ancestor of the Buf. strain acquired DNA by recombination from a recent ancestor of the 108 strain. The Bo10 tree differed from the concatenated tree in two, more subtle, ways. First, the North American strains were very close to the European strains, with KS values less than half those seen for other genes, suggesting that a recombination event occurred between the ancestors of the two groups. Second, midpoint rooting placed the 130, 108 and Buf. strains together with M40 on one side of the root. Although an alternative rooting with the same topology as the concatenated tree is possible, the KS values for Bo10 between M40 and the buffalo strains are lower than those for other genes (except Bo17) and lower than those for Bo10 between M40 and the taurine strains. This suggests the occurrence of a recombination event between the ancestor of the buffalo strains and an ancestor of the M40 strain. However, the KS values do not provide a clear indication of the direction of the exchange in this case. In the ORF 71 tree, the North American strains clustered within the African buffalo clade (Fig. 2bUp), suggesting that the ancestor of the North American strains acquired the ORF 71 sequence from an ancestor of the 108 and 130 strains. In the ORF 73 tree, the three buffalo strains did not form a monophyletic clade. This topology could be explained as the result of an acquisition of sequence by an ancestor of the 108 and 130 strains from the common ancestor of the European and North American strains, but further investigation of this region (see below) revealed a complex mosaic of evolutionary histories within this gene. Finally, the Bo17 tree differed from all others in placing the M40 strain close to the clade of taurine cattle viruses, and Bo17 KS values between M40 and the five taurine viruses were unusually low. This suggests that an ancestor of the M40 strain acquired sequence encompassing the Bo17 ORF from a common ancestor of the North American and European strains.


View this table:
[in this window]
[in a new window]
 
Table 2. Synonymous-substitution rates (KS) among BoHV-4 sequences

 
The three groups of strains defined previously from genome-restriction profiles comprised (i) the European strains, (ii) the North American strains and (iii) others, termed ‘unclassifiable’ (Markine-Goriaynoff et al., 2003aDown). It is now apparent that the unclassifiable strains are not a true group, as they are much more heterogeneous than the two other groups, are not monophyletic and have complex relationships due to interlineage recombination (Fig. 2Up).

The ORF 71–73 region of the BoHV-4 genome contains recombination breakpoints
The three contiguous genomic regions ORF 71, IG 71–73 and ORF 73 yielded three different trees (Fig. 2Up), suggesting the presence of (at least) two recombination breakpoints within this ORF 71–73 region. As recombination need not occur at the boundaries of genes, the distribution of informative sites within this region was examined in order to localize the putative breakpoints. The three European strains were considered as a single taxon, as were the two North American strains and two of the buffalo strains (108 and 130); thus, together with strains Buf. and M40, there were five taxa in the analysis, which may be connected by 15 possible unrooted tree topologies (Fig. 3Down). There were 24 phylogenetically informative sites, including a tetranucleotide at position 976–979 (GTCT in Buf. and M40; ACAA in the other strains) that we considered as a single site (Table 3Down). Each informative site was consistent with, i.e. would require a single mutation on, three of the 15 possible topologies; any of the 12 other topologies would require two mutations at the same site on different branches of the tree. Within ORF 71, there were five (of seven) informative sites consistent with topology 12, the tree found for ORF 71 (Fig. 2bUp); across the rest of the region, there was only one other such site, at the far end of ORF 73. Over the intergenic region and ORF 73, 12 (of 17) informative sites were consistent with topology 1, the tree found for ORF 73 (Fig. 2bUp). In a permutation test, this array of informative sites for topologies 1 and 12 was significantly non-random (P<0.01). These data suggest that a recombination breakpoint occurred between the informative sites at 514 and 591, close to the end of ORF 71.


Figure 3
View larger version (8K):
[in this window]
[in a new window]
 
Fig. 3. Analysis of the ORF 71–ORF 73 region for recombination breakpoints. The nine strains were distributed into five groups: the three European strains (Eur.) V. test, LVR 140 and MOVAR; the two North American strains (Am.) DN599 and 66-p-347; two of the African buffalo strains (130) 130 and 108; the third African Buf. strain (Buf.); and the Indian zebu strain (M40). The 15 possible unrooted tree topologies for the five groups are shown.

 

View this table:
[in this window]
[in a new window]
 
Table 3. The 24 phylogenetically informative sites within the ORF 71–73 region

A variable tetranucleotide at position 976–979 in the alignment is considered as a single site. The mutation types are indicated as either a transition (Ti) or a transversion (Tv), and as either synonymous (s) or non-synonymous (n). Strains with the same sequence appear to one side of the slash under ‘configuration’. Each configuration is consistent with three of the 15 tree topologies (Fig. 3Up), as indicated by black dots. Eur., the three European strains (V. test, LVR 140 and MOVAR); Am., the two North American strains (DN599 and 66-p-347); 130, two of the African buffalo strains (130 and 108); Buf., the third African Buf. strain; M40, the Indian zebu strain.

 
The tree obtained for IG 71–73, as well as for the concatenated sequences (Fig. 2bUp), corresponds to topology 3. Many of the informative sites across IG 71–73 and ORF 73 consistent with topology 1 were also consistent with topology 3 (Table 3Up). However, within the region from 930 to 1202 were four sites consistent with topology 1, but not with topology 3; these include the tetranucleotide site, which seems unlikely to have undergone similar mutations in two different branches of the tree. In this context, the important difference between the two topologies is whether the three buffalo strains cluster (as in topology 3) or not (topology 1). Phylogenetic analysis of the region from 901 to 1230 (from the middle of ORF 73) yielded a tree (Fig. 4Downb) in which strains 108 and 130 fell within the clade of taurine cattle viruses, with strong bootstrap support; this tree matches topology 4. In contrast, a tree (Fig. 4aDown) based on a concatenation of regions 550–900 and 1231–1456 (IG 71–73 plus the two ends of ORF 73) was the same as the overall concatenated tree (Fig. 2aUp). Thus, the earlier tree for ORF 73 (Fig. 2bUp) does not appear to reflect the evolutionary history of any sequence within that gene, but was a compromise between the two histories for different segments within the gene (Fig. 4Down). Simultaneous analysis for two breakpoints within the array of the informative sites supporting topology 3 or topology 4 across IG 71–73 and ORF 73 yielded significant (P=0.005) results in a permutation test. These data suggest that there are two recombination breakpoints within ORF 73, representing the introduction of an internal segment of the gene from an ancestor of the European strains into the ancestor of strains 108 and 130.


Figure 4
View larger version (10K):
[in this window]
[in a new window]
 
Fig. 4. Phylogenetic analyses of different regions of the ORF 73 gene. ORF 73 lies on the complementary strand. The trees were derived from (a) a concatenation of the ORF 71–ORF 73 intergenic region and nt 762–557 and 226–1 of ORF 73 and(b) nt 556–227 of ORF 73. Values on internal branches indicate the percentage of bootstrap replicates in which the branch was found. Horizontal branches are drawn to scale. Bar,0.01 nucleotide substitution per site.

 
Further examination of the alignments for ORF 16, Bo10 and Bo17 did not reveal any signs of mosaicism, indicating that the recombination events described above for those genes involved patches larger than the genes themselves.

Estimation of divergence time between BoHV-4 phylogeographical clades
In the concatenated sequence tree inferred above to reflect the evolutionary history of the clonal backbone of the BoHV-4 genome (Fig. 2aUp), the various strains form phylogeographical clades. To estimate the divergence times of these clades, average KS values across the five gene sequences were calculated for each pairwise comparison of strains (Table 2Up). As genes involved in recombination events would have different times of divergence from those genes inherited clonally, a corrected average was calculated, excluding pairwise comparisons involving sequences that may have been acquired by recombination. Mean corrected averages were then calculated for groups of pairwise comparisons that diverged at the same point in the tree. Assuming a rate of evolution of 5x10–8 substitutions per site per year for BoHV-4 (Markine-Goriaynoff et al., 2003aDown), estimated times of divergence between the different clades were calculated from the mean corrected KS values. These data suggest that the initial split of the M40 strain lineage from the ancestor of the other strains occurred around 850 000 years ago, whilst the last common ancestor of the taurine and buffalo strains existed around 730 000 years ago. The separation between the European and North American clades was estimated at around 260 000 years ago, whilst the last common ancestor of the buffalo strains existed around 280 000 years ago. These estimates of the divergence times between the clades, as well as the various recombination events detected in the present study, are summarized in Fig. 5Down.


Figure 5
View larger version (9K):
[in this window]
[in a new window]
 
Fig. 5. Evolutionary history of the BoHV-4 strains. The tree derived from concatenated sequences (Fig. 2aUp) is inferred to reflect the clonal evolution of the BoHV-4 genome, upon which various inter-strain recombination events (indicated by arrows) have been superimposed; double-headed arrows indicate exchanges where the direction of transfer could not be determined. The estimated dates of divergence of major lineages (from Table 2Up) are shown below the relevant nodes.

 

   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The mammalian herpesviruses appear to have generally co-evolved and co-speciated with their hosts over the last 50 million years or more (Davison, 2002Down; McGeoch, 2001Down; McGeoch et al., 1995Down; Sharp, 2002Down). Nevertheless, there are clear instances where cross-species transmission must have occurred in the past, for example to give rise to bovine and wallaby viruses that lie within the radiation of primate alphaherpesviruses. Whilst evolving with their hosts, gammaherpesviruses have captured a variety of genes from their genomes (McGeoch, 2001Down; Shackelton & Holmes, 2004Down). In an investigation of the origin of one such gene, the Bo17 gene of BoHV-4, we found that the gene was acquired from the genome of African buffalo, indicating that BoHV-4 had been transmitted from buffalo to cattle at some time after the gene-acquisition event (Markine-Goriaynoff et al., 2003aDown). Nevertheless, phylogenetic analyses of the Bo17 genes from various cattle and buffalo strains of BoHV-4 suggested that the cross-species transmission event could have occurred up to 0.7 million years ago and that the virus may have been co-evolving with cattle since then (Markine-Goriaynoff et al., 2003aDown). For the same group of strains, we have now analysed sequences of four other genes from locations across the BoHV-4 genome. These data provide evidence of a number of recombination events during the divergence of the strains, including one involving the Bo17 gene, so that it is now apparent that the evolutionary history of the BoHV-4 strains is very different from that inferred previously.

Although a number of recombination events were detected, such that none of the trees for the five genes investigated here were the same, it is clear that recombination between divergent strains has been relatively infrequent during the diversification of BoHV-4. With even a modest rate of recurrent recombination, the evolutionary histories of different genes would become quite scrambled. However, here it was possible to discern the underlying tree reflecting the clonal evolution of strains, upon which the recombination events have been superimposed (Fig. 5Up). In that tree, the cattle strains do not form a single clade, but rather the M40 strain from zebu represents the earliest-diverging lineage. From this tree alone, the most-parsimonious interpretation would be a single transmission from cattle to African buffalo. However, the estimated divergence times for the various lineages in the tree do not correspond well with the evolutionary history of domestic cattle. Although these dates rely on the use of a synonymous-substitution rate estimated for the Bo17 gene (Markine-Goriaynoff et al., 2003aDown), other BoHV-4 genes seem (as expected) to have a similar rate (Table 2Up). Estimates of the synonymous-substitution rate for herpes simplex viruses (3x10–8; Hatwell & Sharp, 2000Down) and during the divergence of Epstein–Barr virus and Cercopithecine herpesvirus 12 (2x10–8–4x10–8; Hughes, 2002Down) are a little lower than, but of the same order as, that used here for BoHV-4 (5x10–8); thus, the estimated dates are unlikely to be substantially wrong, but must be taken with some caution. Thus, the split between the M40 strain and others, at around 850 000 years ago, greatly predates the estimated time of divergence between the ancestors of zebu and taurine cattle, placed at 275 000 years ago (Bradley et al., 1996Down). Similarly, the split between the European and North American BoHV-4 strains from taurine cattle was estimated at around 260 000 years ago, whereas taurine cattle were domesticated about 10 000 years ago and the ancestors of the North American domestic cattle were taken there from Europe less than 500 years ago. These observations suggest strongly that BoHV-4 has not been co-evolving with cattle during most of the evolutionary history portrayed in the tree. Taken together with our previous conclusion that BoHV-4 infected African buffalo around 1.5 million years ago, when the Bo17 gene was acquired (Markine-Goriaynoff et al., 2003aDown), and the observation that the prevalence of BoHV-4 infection is higher in buffalo than in cattle (Dewals et al., 2005Down), these data point to the African buffalo as the original BoHV-4 reservoir, implying at least three recent cross-species transmissions to cattle to give rise to the lineages seen in zebu and in European and North American taurine cattle.

In this interpretation, the overall tree reflects divergence of, and recombination between, BoHV-4 strains infecting African buffalo over the last million years. However, the three strains isolated from buffalo form a subclade within the tree, with a much more recent coalescence time (around 280 000 years ago) than the tree as a whole (around 850 000 years ago). Buffalo are distributed widely across sub-Saharan Africa, but all three strains analysed here were sampled from Kenya and, so, it is possible that more divergent strains of BoHV-4 infect buffalo elsewhere in Africa. However, the great rinderpest epidemic of the 1890s caused mortality levels reaching 80–90 % among African ruminant species, notably causing a crash in the African buffalo population (Mack, 1970Down; Wenink et al., 1998Down). It would not be surprising if this led to a dramatic reduction in the diversity of buffalo BoHV-4 lineages. Thus, the various cattle strains of BoHV-4 would reflect more divergent lineages that survived because they were transmitted from buffalo to other species before the rinderpest epidemic. Further isolation and phylogenetic analysis of buffalo BoHV-4 strains in various sub-Saharan areas are required to test these two hypotheses.

The transmission of BoHV-4 from buffalo to cattle probably involved intermediate host species. It has been shown that BoHV-4 is able to replicate in a broad range of mammal species both in vivo and in vitro (Gillet et al., 2004Down). Certainly, any transmission prior to the domestication of cattle is unlikely to have been direct, as their ancestor, B. primigenius, was found only in Eurasia and North Africa (Troy et al., 2001Down), whereas African buffalo are restricted to sub-Saharan Africa (Simonsen et al., 1998Down; Van Hooft et al., 2002Down). Also, the relatively large divergence between North American and European BoHV-4 strains, estimated to have diverged around 260 000 years ago, is surprising if both lineages had an origin in European cattle. No North American-type BoHV-4 strains have ever been described in Europe (Bublot et al., 1990Down, 1991Down), whereas European-type BoHV-4 strains have been described in North American cattle (Henry et al., 1986Down; Shen et al., 1992Down). Thus, it seems more likely that the divergence of North American strains reflects a host jump in prehistory followed by a recent transmission to cattle. North American-like strains of BoHV-4 have been recovered from the North American bison (Bison bison) (Todd & Storz, 1983Down) and it is plausible that bison were the reservoir for the ancestors of the North American strains. Bison are believed to have first entered eastern Beringia from Asia during the middle Pleistocene (300 000–130 000 years ago) and then moved southward into central North America around 130 000–75 000 years ago (Shapiro et al., 2004Down). This scenario argues for a very recent cross-species transmission from bison to cattle, whilst the European-like strains in North America are probably due to the recent wide flow of livestock breeds associated with modern selective breeding. It would be interesting to isolate a collection of BoHV-4 strains from bison and cattle populations in North America and to establish the proportion and the distribution of North American-type and European-type BoHV-4 strains. European bison became extinct in the wild in 1919; attempts to detect evidence of BoHV-4 infection in a reintroduced population from a single area have proved unsuccessful (Borchers et al., 2002Down).

In conclusion, taken together with recent serological and viral isolation data (Dewals et al., 2005Down; Rossiter et al., 1989Down), the present study further supports the role of the African buffalo as the original natural host species of BoHV-4. It suggests that while BoHV-4 has been coevolving with the African buffalo for at least the last 1.5 million years, it has been transmitted to cattle on at least three independent occasions, and via intermediate species much more recently.


   ACKNOWLEDGEMENTS
 
A. V. is Senior Research Associate of the ‘Fonds National Belge de la Recherche Scientifique’ (FNRS). N. M.-G. and L. G. are Postdoctoral Researchers of the FNRS and B. D. is a Research Fellow of the FNRS. M. T. is a Research Fellow of the Belgian ‘Fonds pour la formation à la Recherche dans l'Industrie et dans l'Agriculture’.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Bartha, A., Juhasz, M. & Liebermann, H. (1966). Isolation of a bovine herpesvirus from calves with respiratory disease and keratoconjunctivitis. A preliminary report. Acta Vet Acad Sci Hung 16, 357–358.[Medline]

Borchers, K., Brackmann, J., Wolf, O., Rudolph, M., Glatzel, P., Krasinska, M., Krasinski, Z. A. & Frölich, K. (2002). Virologic investigations of free-living European bison (Bison bonasus) from the Bialowieza Primeval Forest, Poland. J Wildl Dis 38, 533–538.[Abstract]

Bradley, D. G., MacHugh, D. E., Cunningham, P. & Loftus, R. T. (1996). Mitochondrial diversity and the origins of African and European cattle. Proc Natl Acad Sci U S A 93, 5131–5135.[Abstract/Free Full Text]

Bublot, M., Van Bressem, M.-F., Thiry, E., Dubuisson, J. & Pastoret, P.-P. (1990). Bovine herpesvirus 4 genome: cloning, mapping and strain variation analysis. J Gen Virol 71, 133–142.[Abstract/Free Full Text]

Bublot, M., Wellemans, G., Van Bressem, M.-F., Dubuisson, J., Pastoret, P.-P. & Thiry, E. (1991). Genomic diversity among bovine herpesvirus 4 field isolates. Arch Virol 116, 1–18.[CrossRef][Medline]

Burrows, J. M., Khanna, R., Sculley, T. B., Alpers, M. P., Moss, D. J. & Burrows, S. R. (1996). Identification of a naturally occurring recombinant Epstein-Barr virus isolate from New Guinea that encodes both type 1 and type 2 nuclear antigen sequences. J Virol 70, 4829–4833.[Abstract]

Davison, A. J. (2002). Evolution of the herpesviruses. Vet Microbiol 86, 69–88.[CrossRef][Medline]

Dewals, B., Gillet, L., Gerdes, T., Taracha, E. L. N., Thiry, E. & Vanderplasschen, A. (2005). Antibodies against bovine herpesvirus 4 are highly prevalent in wild African buffaloes throughout eastern and southern Africa. Vet Microbiol 110, 209–220.[CrossRef][Medline]

Gillet, L., Minner, F., Detry, B. & 7 other authors (2004). Investigation of the susceptibility of human cell lines to bovine herpesvirus 4 infection: demonstration that human cells can support a nonpermissive persistent infection which protects them against tumor necrosis factor alpha-induced apoptosis. J Virol 78, 2336–2347.[Abstract/Free Full Text]

Hatwell, J. N. & Sharp, P. M. (2000). Evolution of human polyomavirus JC. J Gen Virol 81, 1191–1200.[Abstract/Free Full Text]

Henry, B. E., II, Ota, R. & Evermann, J. F. (1986). Genetic relatedness of disease-associated field isolates of bovid herpesvirus type 4. Am J Vet Res 47, 2242–2246.[Medline]

Hughes, A. L. (2002). Origin and evolution of viral interleukin-10 and other DNA virus genes with vertebrate homologues. J Mol Evol 54, 90–101.[CrossRef][Medline]

Li, W. H. (1993). Unbiased estimation of the rates of synonymous and nonsynonymous substitution. J Mol Evol 36, 96–99.[CrossRef][Medline]

Li, H., Gailbreath, K., Flach, E. J. & 9 other authors (2005). A novel subgroup of rhadinoviruses in ruminants. J Gen Virol 86, 3021–3026.[Abstract/Free Full Text]

Mack, R. (1970). The great African cattle plague epidemic of the 1890s. Trop Anim Health Prod 2, 210–219.

Markine-Goriaynoff, N., Georgin, J.-P., Goltz, M., Zimmermann, W., Broll, H., Wamwayi, H. M., Pastoret, P.-P., Sharp, P. M. & Vanderplasschen, A. (2003a). The core 2 beta-1,6-N-acetylglucosaminyltransferase-mucin encoded by bovine herpesvirus 4 was acquired from an ancestor of the African buffalo. J Virol 77, 1784–1792.[Abstract/Free Full Text]

Markine-Goriaynoff, N., Minner, F., de Fays, K., Gillet, L., Thiry, E., Pastoret, P. P. & Vanderplasschen, A. (2003b). L'herpèsvirus bovin 4. Ann Méd Vét 147, 215–247 (in French).

McGeoch, D. J. (2001). Molecular evolution of the {gamma}-Herpesvirinae. Philos Trans R Soc Lond B Biol Sci 356, 421–435.[Abstract/Free Full Text]

McGeoch, D. J., Cook, S., Dolan, A., Jamieson, F. E. & Telford, E. A. R. (1995). Molecular phylogeny and evolutionary timescale for the family of mammalian herpesviruses. J Mol Biol 247, 443–458.[CrossRef][Medline]

Mohanty, S. B., Hammond, R. C. & Lillie, M. G. (1971). A new bovine herpesvirus and its effect on experimentally infected calves. Brief report. Arch Gesamte Virusforsch 33, 394–395.[CrossRef][Medline]

Moreno-Lopez, J., Goltz, M., Rehbinder, C., Valsala, K. V. & Ludwig, H. (1989). A bovine herpesvirus (BHV-4) as passenger virus in ethmoidal tumours in Indian cattle. Zentralbl Veterinarmed B 36, 481–486.[Medline]

Norberg, P., Bergström, T., Rekabdar, E., Lindh, M. & Liljeqvist, J.-Å. (2004). Phylogenetic analysis of clinical herpes simplex virus type 1 isolates identified three genetic groups and recombinant viruses. J Virol 78, 10755–10764.[Abstract/Free Full Text]

Poole, L. J., Zong, J.-C., Ciufo, D. M., Alcendor, D. J., Cannon, J. S., Ambinder, R., Orenstein, J. M., Reitz, M. S. & Hayward, G. S. (1999). Comparison of genetic variability at multiple loci across the genomes of the major subtypes of Kaposi's sarcoma-associated herpesvirus reveals evidence for recombination and for two distinct types of open reading frame K15 alleles at the right-hand end. J Virol 73, 6646–6660.[Abstract/Free Full Text]

Rossiter, P. B., Gumm, I. D., Stagg, D. A., Conrad, P. A., Mukolwe, S., Davies, F. G. & White, H. (1989). Isolation of bovine herpesvirus-3 from African buffaloes (Syncerus caffer). Res Vet Sci 46, 337–343.[Medline]

Saitou, N. & Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4, 406–425.[Abstract]

Shackelton, L. A. & Holmes, E. C. (2004). The evolution of large DNA viruses: combining genomic information of viruses and their hosts. Trends Microbiol 12, 458–465.[CrossRef][Medline]

Shapiro, B., Drummond, A. J., Rambaut, A. & 24 other authors (2004). Rise and fall of the Beringian steppe bison. Science 306, 1561–1565.[Abstract/Free Full Text]

Sharp, P. M. (2002). Origins of human virus diversity. Cell 108, 305–312.[CrossRef][Medline]

Shen, D. T., Burger, D., Tong, G. Z. & Gorham, J. R. (1992). Molecular cloning and physical mapping of bovine herpesvirus 4 strain DN 599 and comparison with two American field-isolates. Arch Virol 127, 75–87.[CrossRef][Medline]

Simonsen, B. T., Siegismund, H. R. & Arctander, P. (1998). Population structure of African buffalo inferred from mtDNA sequences and microsatellite loci: high variation but low differentiation. Mol Ecol 7, 225–237.[CrossRef][Medline]

Storz, J. (1968). Comments on malignant catarrhal fever. J Am Vet Med Assoc 152, 804–806.

Thiry, E., Pastoret, P.-P., Dessy-Doize, C., Hanzen, C., Calberg-Bacq, C. M. & other authors (1981). Reactivation of a herpesvirus in cells cultured from a testicle of a bull showing orchitis and azoospermia. Ann Méd Vét 125, 207–214 (in French).

Thiry, E., Bublot, M., Dubuisson, J., Van Bressem, M.-F., Lequarre, A.-S., Lomonte, P., Vanderplasschen, A. & Pastoret, P.-P. (1992). Molecular biology of bovine herpesvirus type 4. Vet Microbiol 33, 79–92.[CrossRef][Medline]

Thiry, E., Meurens, F., Muylkens, B. & 7 other authors (2005). Recombination in alphaherpesviruses. Rev Med Virol 15, 89–103.[CrossRef][Medline]

Thompson, J. D., Higgins, D. G. & Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22, 4673–4680.[Abstract/Free Full Text]

Todd, W. J. & Storz, J. (1983). Morphogenesis of a cytomegalovirus from an American bison affected with malignant catarrhal fever. J Gen Virol 64, 1025–1030.[Abstract/Free Full Text]

Troy, C. S., MacHugh, D. E., Bailey, J. F., Magee, D. A., Loftus, R. T., Cunningham, P., Chamberlain, A. T., Sykes, B. C. & Bradley, D. G. (2001). Genetic evidence for Near-Eastern origins of European cattle. Nature 410, 1088–1091.[CrossRef][Medline]

Vanderplasschen, A., Markine-Goriaynoff, N., Lomonte, P. & 8 other authors (2000). A multipotential beta-1,6-N-acetylglucosaminyl-transferase is encoded by bovine herpesvirus type 4. Proc Natl Acad Sci U S A 97, 5756–5761.[Abstract/Free Full Text]

Van Hooft, W. F., Groen, A. F. & Prins, H. H. T. (2002). Phylogeography of the African buffalo based on mitochondrial and Y-chromosomal loci: Pleistocene origin and population expansion of the Cape buffalo subspecies. Mol Ecol 11, 267–279.[CrossRef][Medline]

Van Opdenbosch, E., Wellemans, G. & Oudewater, J. (1986). Toevallige isolatie van het boviene herpesvirus 4 uit de long van een schaap. Vlaams Diergeneesk Tijdschr 55, 432 (in Dutch).

Walling, D. M. & Raab-Traub, N. (1994). Epstein-Barr virus intrastrain recombination in oral hairy leukoplakia. J Virol 68, 7909–7917.[Abstract/Free Full Text]

Wellemans, G., Antoine, H., Broes, A., Charlier, G. & Van Opdenbosch, E. (1983). Isolation of a herpesvirus from cows affected with post partural metritis. Ann Méd Vét 127, 481–482 (in French).

Wellenberg, G. J., van Rooij, E. M. A., Maissan, J. & Van Oirschot, J. T. (1999). Evaluation of newly developed immunoperoxidase monolayer assays for detection of antibodies against bovine herpesvirus 4. Clin Diagn Lab Immunol 6, 447–451.

Wenink, P. W., Groen, A. F., Roelke-Parker, M. E. & Prins, H. H. T. (1998). African buffalo maintain high genetic diversity in the major histocompatibility complex in spite of historically known population bottlenecks. Mol Ecol 7, 1315–1322.[CrossRef][Medline]

Zimmermann, W., Broll, H., Ehlers, B., Buhk, H.-J., Rosenthal, A. & Goltz, M. (2001). Genome sequence of bovine herpesvirus 4, a bovine Rhadinovirus, and identification of an origin of DNA replication. J Virol 75, 1186–1194.[Abstract/Free Full Text]

Received 9 December 2005; accepted 17 February 2006.


This article has been cited by other articles:


Home page
J. Virol.Home page
B. Muylkens, F. Farnir, F. Meurens, F. Schynts, A. Vanderplasschen, M. Georges, and E. Thiry
Coinfection with Two Closely Related Alphaherpesviruses Results in a Highly Diversified Recombination Mosaic Displaying Negative Genetic Interference
J. Virol., April 1, 2009; 83(7): 3127 - 3137.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary tables
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Dewals, B.
Right arrow Articles by Vanderplasschen, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Dewals, B.
Right arrow Articles by Vanderplasschen, A.
Agricola
Right arrow Articles by Dewals, B.
Right arrow Articles by Vanderplasschen, A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS