Skip to main content

A phylogenetic analysis of Bovine Viral Diarrhoea Virus (BVDV) isolates from six different regions of the UK and links to animal movement data

Abstract

Bovine Viral Diarrhoea Virus (BVDV) is a pestivirus which infects cattle populations worldwide and is recognised as a significant source of economic loss through its impact on health and productivity. Studies investigating the molecular epidemiology of BVDV can give invaluable information about the diversity of viral strains present in a population and this, in turn, can inform control programs, drive vaccine development and determine likely infection sources. The current study investigated 104 viral isolates from forty farms across the UK. Through phylogenetic and nucleotide sequence analysis of the 5′UTR and Npro regions of the isolates investigated, it was determined that BVDV 1a was the predominant sub-genotype. However, BVDV 1b, 1e and 1i were also identified and, for the first time in the UK, BVDV 1d. Through analysis of animal movement data alongside the phylogenetic analysis of these BVD isolates, it was possible to link animal movements to the viral isolates present on several premises and, for the first time, begin to elucidate the routes of viral transmission. With further work, this type of analysis would enable accurate determination and quantification of the true biosecurity risk factors associated with BVDV transmission.

Introduction

Bovine Viral Diarrhoea Virus (BVDV) infects cattle worldwide and is a cause of significant losses in both beef and dairy systems. The virus is a member of the pestivirus genus in the Flaviviridae family. It is a positive sense, single stranded RNA virus of approximately 12.5 kb in length[1]. Pestivirus genomes are flanked by 5′ and 3′ untranslated regions (5′UTR, 3′UTR) and encode 11–12 structural and non-structural proteins (Npro, C, Erns, E1, E2, P7, NS2/3, NS4A, NS4B, NS5A, NS4B). The genus is currently comprised of four recognised species; BVDV-1, BVDV-2, Border Disease Virus (BDV) and Classical Swine Fever Virus (CSFV) and four proposed “atypical” species Giraffe, “HoBi”, Pronghorn Antelope and Bungowannah[27].

Pestiviruses are highly variable both antigenically and genetically, hence each species is further classified. BVDV is currently divided into two genotypes[8]. More recent molecular analysis, primarily based upon sequence analysis of the 5′UTR and the non-structural Npro gene, has resulted in further sub-classification such that BVDV genotype 1 currently contains twelve sub-genotypes; BVDV 1a-k plus Deer[9, 10]. Six further sub-genotypes of BVDV1 have been proposed (1l, 1m, 1n, 1o and 1p with BVDV 1l being designated twice for diverse sets of isolates) although the literature concerning some of these members requires further clarification and is discussed in detail later[3, 1115]. BVDV 2 is currently separated into at least two sub-genotypes[1618].

Studies investigating the molecular epidemiology of BVDV can provide invaluable information about the diversity of viral strains present in a population and, in turn, inform control programs, drive vaccine development and determine likely infection sources. In a survey of bovine pestiviruses published in 1999, 62 BVD viral isolates across England and Wales were sequenced and were predominantly found to be BVDV 1a although BVDV 1b and 1i were also identified[10, 19]. A more recent study conducted between 2004 and 2009 indicated that whilst BVDV 1a was still the predominant sub-genotype, and that 1b and 1i were still circulating, the genetic diversity of BVDV in England and Wales has increased further to include BVDV 1e and 1f[20]. Here, we provide results of a further phylogenetic analysis of pestiviruses isolated from cattle across several regions of the UK including some of the most cattle dense regions in England with findings that largely support those of Strong et al.[20]. Furthermore, we have analysed cattle movement information for a number of farms participating in this study. For the first time with BVDV, we have demonstrated that it is possible to elucidate potential sources and routes of infection if animal movements are correlated with phylogenetic analyses.

Materials and methods

Farm distribution for sample collection

Samples were collected from PI animals identified on farms that were screened across six regions of the UK between 2006 and 2011 (see Table 1 and Figure 1). In summary, four farms were sampled in South-east England, fourteen in South-west England, seven in East Anglia, five in Northern England, one in Wales and nine in Eastern Scotland. Samples were labelled with “Farm number – PI number” such that the farm of origin and order of PI identification is consistent throughout. For example, Isolate 33–3 is Farm 33-third PI identified. This numbering also allows cross referencing to the farms previously described by Booth and Brownlie[21].

Table 1 Regional distribution and details of BVDV infected premises investigated in this study
Figure 1
figure 1

Regional distribution of infected farms investigated per Animal Health Divisional Office region relative to the cattle density in Great Britain in 2008. The cattle population density and Animal Health Divisional Office (AHDO) Boundaries illustrated in Great Britain are as of 1st June 2008 (“The Cattle Book 2008” © Crown copyright. Reproduced with kind permission of Ordnance Survey).

Virus RNA isolation and cDNA synthesis

Blood samples collected from PI animals were either heparinised or clotted (this was at the discretion of the clinician taking the sample). RNA was isolated from heparinised whole blood using a QIAamp RNA Blood Mini Kit (QIAGEN Ltd., Manchester, UK) and from serum using a QIAGEN QIAamp Viral RNA Kit (QIAGEN Ltd.) according to the manufacturer’s instructions.

cDNA was synthesised from isolated RNA using Invitrogen Superscript II reverse transcription kit (Invitrogen, Paisley, UK). For each cDNA reaction, 5 μL of RNA was mixed with 3 μL of pd(N)6 Random Hexamers (GE Healthcare, Hatfield, UK) and incubated at 70°C for 10 min. 8 μL of 5× Superscript II First Strand Buffer (Invitrogen), 4 μL of 0.1M DTT (Invitrogen), 2 μL of 10 mM dNTPs (Promega, Southampton, UK), 16 μL nuclease-free dH2O and 1 μL of 40 U/μL RNasin (Promega) was then added. After incubating the mixture at 42°C for 2 min, 1 μL of 200 U/μL Superscript II Reverse Transcriptase enzyme (Invitrogen) was then added to each sample and the incubation continued for a further 50 min followed by 10 min at 72°C. RNA samples were stored at −70°C and cDNA samples were stored at −20°C for future use.

RT-PCR and DNA sequencing

For each viral isolate, a 288bp and 428bp region of the 5′UTR and Npro regions respectively were amplified. At least two PCR reactions for each region, using a high fidelity Pfu DNA polymerase enzyme (Stratagene, Stockport, UK), were performed on each sample. If consensus sequence could not be obtained from sequence analysis of two PCR reactions, further reactions and sequencing were performed in order to achieve this. Primers 324 and 326 were used for amplification of a 288bp region of the 5′ UTR as described by Vilcek et al.[22]. Primers BD1[23] and a modified version of BD3[10], named BD3A, were used for amplification of a 428bp region of Npro under the same conditions described by Vilcek et al.[10] for the BD1/BD3 pair. BD3A was designed in house and the sequence of the modified primer is: 5′- CAT CCA TCT ATR CAY AYA TAA ATR TGG TAC - 3′.

For PCR reactions, the following was setup; 2 μL template cDNA, 5 μL Forward primer (10 pmol/μL), 5 μL Reverse primer (10 pmol/μL), 5 μL 10 × Pfu DNA polymerase reaction buffer (Stratagene), 1 μL 10mM dNTP, 0.5 μL Pfu DNA polymerase 2.5 U/μL (Stratagene) and 31.5 μL dH20. The PCR conditions were as follows: 94°C for 3 min initial denaturing followed by 30 cycles of 94°C for 1 min, 52°C for 1 min, 72°C for 1 min and then final elongation at 72°C for 10 min. Amplicons of the expected sizes (288bp for 5′UTR and 428bp for Npro) were gel purified or PCR purified (based on whether or not a single individual amplicon was observed) using QIAquick /PCRGel Extraction Kit (QIAGEN). Purified PCR products were sequenced by QIAGEN Genomic Services (Hilden, Germany) using Dye Terminator cycle sequencing performed with the ABI PRISM BiG DYE V3.1 Terminator Cycle Sequencing. At least two independent PCR reactions were sequenced for each sample using forward and reverse primers 324/326 or BD1/BD3A for 5′UTR and Npro respectively.

DNA sequencing and phylogenetic analysis

The sequences received were initially analysed in Vector NTI Version 11 (Invitrogen) and consensus sequences were assembled in the Contig Express feature of Vector NTI where they were checked for discrepancies. Primer sequences were clipped from each consensus prior to phylogenetic analysis conducted using MEGA version 5[24]. 5′UTR and Npro sequences in FASTA format were imported into MEGA 5 and sequence alignments performed using ClustalW. Phylogenetic analysis was performed using the Neighbour Joining Method[25] and evolutionary distances were calculated using the Kimura 2-parameter method[26]. Bootstrap analysis of the resultant tree was performed using 1000 replicates with deletion of gaps in the alignment data. Bootstrap figures >70% are reported[27].

Cattle movement analysis

Sixty-two farms gave permission to obtain animal movement records relating to their premises. This network included thirty-one farms from which PI samples were genotyped, twenty-six farms where BVDV surveillance was conducted but no PI animals identified and five farms on which PI animals were identified but genotyping not yet performed. These are listed as Farms 1–62 in Table 1. Animal movements recorded by the British Cattle Movement Service (BCMS) were sourced from the Department for Environment, Food and Rural Affairs (DEFRA) veterinary surveillance group “Rapid Analysis and Detection of Animal-related Risks” (RADAR). RADAR supplied movement data for Farms 1–62 to include cattle movements between study farms, through intermediate premises which were not in the study (anonymous), markets and showgrounds. Movement data was supplied in text file format which was imported into Microsoft Access. Movement details were extracted from the database through a series of queries producing matrices of movements between relevant premises. Movement links were scrutinised alongside the phylogenetic analysis to assess links between virus strain and animal movements.

Results

Distribution of regions sampled

One hundred and four viral isolates were collected from PI cattle on forty farms across the UK. The geographical location of the farms and the way in which the samples were submitted is summarised in Table 1. Samples from ten farms were submitted to the BVDV research group at The Royal Veterinary College (RVC) by individual practitioners or laboratories. The remaining samples came from members of organised BVDV control programmes. Six farms surveyed across East Anglia were members of an eradication programme promoted by Holstein UK (HUK) which, received part funding from the English Beef and Lamb Executive (EBLEX). Twenty-four of the farms surveyed were involved in an eradication scheme set up jointly by RVC and the Scottish Agricultural College (SAC) which recruited farms in Somerset, Southeast England, Northern England and Eastern Scotland. The Somerset/South-Eastern eradication programme is described in detail by Booth and Brownlie[21]. The distribution of farms sampled compared to UK cattle density is illustrated in Figure 1 and whilst some cattle dense regions were sampled, it should be noted that the sampling cannot be considered representative of Wales, the Welsh borders or the far southwest of England.

Phylogenetic and sequence analysis

Of the 104 viral isolates investigated, 101/104 5′UTR and 98/104 Npro sequences were obtained for phylogenetic analysis and sequence comparison. Table 2 gives details of the number of PIs investigated on each infected farm and the sub-genotypes identified. The majority of viral isolates, 85% (88/104), were typed as BVDV 1a based upon analysis of both 5′UTR and Npro and comparison to a set of reference strains detailed in Table 3. However, also identified were BVDV types 1b (Farms 18 and 71), 1d (Farm 67), 1e (Farms 1 and 68) and 1i (Farms 58 and 69). There was no noticeable geographical distribution of isolates with BVDV 1a being identified in all regions sampled except Wales where only one farm had been sampled.

Table 2 Strain allocation table
Table 3 Reference strains used in this publication

For the viral isolates investigated, the apparent sub-genotypes indicated by the 5′UTR and Npro phylogenetic trees in Figures 2 and3 respectively were largely supported by nucleotide sequence comparisons to the reference strains in Table 3. However, several issues that confused isolate sub-genotype classification became apparent when using the 5′UTR alone and this is noted in Table 2. The analysis of Npro, shown in Table 2, provided a much clearer definition of the sub-genotypes investigated. The only instance of confusion generated through Npro analysis occurred for isolates 1–1 and 1–2, which both had the highest identity (95%) to the BVDV 1e 5′UTR reference strains, yet upon Npro analysis were most identical (83%) to the BVDV 1a reference strain SD1. However, removal of SD1 from the Npro analysis left isolates 1–1 and 1–2 with closest identity (82%) to BVDV 1e, 1i, 1j and 1m reference strains and only 81% identity to NADL (1a); hence, isolates from Farm 1 clustered phylogenetically with BVDV 1e in both Figures 2 and3 (5′UTR and Npro respectively).

Figure 2
figure 2

Phylogenetic tree of a 245bp region of the 5′UTR. All field isolates are denoted by “Farm Number – PI Number”. * indicates isolates with equal 5′UTR nucleotide sequence identity to both BVDV 1a and 1j reference strains, + indicates isolates with 1% less 5′UTR nucleotide identity to BVDV 1j reference strains than to 1a reference strains, ^ indicates isolates with 1% higher 5′UTR identity to BVDV 1j reference strains than to either 1a or 1c reference strains, # indicates isolates with 1% less 5′UTR nucleotide identity to BVDV 1c reference strains than to 1a reference strains. Thirty six reference strains are included (underlined) each beginning with the genotype and sub-genotype to which they belong (1a-1m and 2) followed by their name; further details can be found in Table 3. Bootstrap figures supported by >70% of 1000 replicates are indicated.

Figure 3
figure 3

Phylogenetic tree of a 380bp region of Npro. All field isolates are denoted by “Farm Number – PI Number”. Twenty five reference strains are included (underlined) each beginning with the genotype and sub-genotype to which they belong (1a-1m and 2) followed by their name; further details can be found in Table 3. Bootstrap figures supported by >70% of 1000 replicates are indicated.

Movement data analysis

Movement data was available for Farms 1–62 in Table 1 for the period January 2005 to December 2010. Farms 63–71 had not given permission for analysis of their animal movement data. In addition to cattle movement data, the yearly infection status was available for Farms 1–42 from the study conducted by Booth and Brownlie[21], thus enabling the construction of an animal movement network linked to farm BVD status between 2005 and 2010 (Figure 4). Indirect links between premises involved, such as animal movements through other farms, markets or showgrounds were infrequent with the latter being identified as the only indirect animal movement link between farms in this study (see below). From analysis of the sequence data and animal movements, it was possible to describe three situations regarding the epidemiology of BVDV transmission.

  1. 1.

    Same farm, same isolate: This situation is demonstrated in both Figures  2 and 3 for all farms in this study where more than one PI isolate was investigated, except for Farms 2, 40 and 59 (see point 2 below). Typically, isolates from PIs born on one premise were seen to cluster phylogenetically with others from the same farm indicating that the infecting strains have very high identity. Analysis of the nucleotide identity matrices for the 5′UTR and Npro sequences where more than one PI was identified on the same farm revealed that they shared 99%–100% nucleotide identity in both regions investigated.

  2. 2.

    Same farm, different isolate: In some instances it is evident that two isolates have caused infection on the same farm; Farms 2, 40 and 59 are examples of this. For Farms 2 and 40, two clusters of PIs are present for each farm in the phylogenetic trees in Figures  2 and 3. Nucleotide sequence analysis for these isolates reveals that for both farms, isolates share 99%–100% identity in both the 5′UTR and Npro regions within each PI cluster. However, between clusters on the same farm, identity in the 5′UTR and Npro regions ranges from 95%–96% and 88%–90% respectively. For Farm 59, only two PIs were identified and these shared nucleotide sequence identities in the 5′UTR and Npro regions of only 95% and 91% respectively thus grouped separately from each other upon phylogenetic analysis (Figures  2 and 3).

  3. 3.

    Different farms, same isolate: This instance appears to occur on several farms if the 5′UTR in Figure  2 is considered alone however, the results indicate that Npro also needs to be considered to make conclusions of this nature. Viral isolates from PIs on Farms 40 and 41 (isolates share 99%–100% 5′UTR nucleotide identity), Farms 20 and 43 (isolates share 99% 5′UTR nucleotide identity), Farms 56 and 62 (isolates share 99%–100% 5′UTR nucleotide identity) and Farms 33 & 34 (isolates share 100% 5′UTR nucleotide identity) are seen to cluster as distinct groups in Figure  2. The Npro phylogenetic tree in Figure  3, confirms the links between Farms 33 & 34 and 40 & 41 with isolates clustering together and achieving 99%–100% nucleotide identity in this region. The links are less well supported by Npro phylogenetic analysis of Farms 20 & 43 and 56 & 62 with isolates from these farm groups returning 97% and 96% Npro nucleotide identity respectively. Whilst these identities are still high, the clustering of isolates from these farms is less convincing in Figure  3 when compared to that seen in Figure  2. Geographically, none of these farms share borders (data not shown); however, the possible animal/material links between farms that presented potential biosecurity pathways for virus transmission are highlighted below.

Figure 4
figure 4

Farm status and animal movement matrix between study Farms 1–42 (2005–2010): Each coloured node represents an individual farm. Farm number is stated alongside the node. The status of each farm is defined by the colour of its node such that: a “Black Node” represents a farm of unknown status, a “Green Node” represents a BVDV free farm, a “Red Node” represents a BVDV infected farm and a “Yellow Node” represents a farm that is currently eradicating BVDV. Numbers of animal movements are indicated alongside the arrow that illustrating the direction of movement. The arrow line width is proportional to the number of movements. Geographic position and relationships are not represented in the figure.

Farms 33 and 34 have been linked by direct animal movements (1 animal in 2007 and 11 animals in 2010) and also through indirect animal movements through Farm 32 throughout the study period. Figure 4 indicates that BVDV infection was detected on Farms 33 and 34 in 2006 and 2007 respectively thus the animal movements in 2010 are unlikely to be responsible for transmission of infection between these farms. It is feasible that the one movement noted in 2007 could have been responsible for the presence of the same BVDV isolate on both premises. However, there are five further premises involved in the movement network with Farms 33 and 34 in Figure 4 (Farms 28, 29, 30, 32 and 36). It can be noted that multiple movements occur between Farm 32 and 33 and Farm 32 and 34. It is quite possible that infected animals have passed between Farms 34 and 33 through Farm 32. However, it is curious that infection was not detected on Farm 32 at any point. Further to this, it can also be noted that BVDV infection was detected on Farms 29 and 36 in 2007 with the former testing bulk milk PCR positive and the latter testing positive in youngstock for BVD antibody in 2007 (data not shown). PI searches were delayed on these farms due to management changes and thus it is believed that any PIs present left these premises prior to whole herd investigations, thus viral isolates were not available from these farms sub-genotyping. The involvement of infected farms in this network illustrates the risk to the other farms involved; biosecurity advice was given to all farms involved in order to mitigate this. Direct animal movements were also noted from Farm 50 to 49 throughout 2006–2010 and from Farm 49 to 50 in 2007 (data not shown), yet Figures 2 and3 do not indicate a viral isolate that is common to both farms suggesting that they had different infection sources and that one had not infected the other.

Whilst the isolates from Farms 40 and 41 appear linked when analysing both the 5′UTR and Npro phylogenetic trees, they were not linked by direct animal movements between their premises. Despite this, the high identity noted in both the 5′UTR and Npro regions of viral isolates from these farms supports the theory that there is a common isolate in both units. It is known to the authors that the owners of Farms 40 and 41 exhibit animals at the same agricultural shows each year and it is quite feasible that animal contact may have occurred at one of these shows resulting in either concurrent infection of both farms with the same BVDV strain or one farm infecting the other at the show. As an adjunct to this, the authors are aware that Farm 40 has supplied colostrum to Farm 41 on several occasions; another potential route for transmission (this practice has now ceased consequent to our advice).

Whilst Farms 56 & 52 and 20 & 43, initially appear to have a common BVDV isolate upon 5′UTR phylogenetic and nucleotide analysis, this is less well supported by analysis of the Npro region. No direct animal movements were identified to link these farms. Upon initial investigation, there was no common link between Farms 56 and 52 since they were more than 400 miles apart. However, both use embryo transfer services which could explain a common source of infection although this would seem unlikely given the analysis of the Npro regions of BVDV isolates from both premises.

Discussion

The predominant BVDV sub-genotype identified on the forty infected farms investigated in this study was BVDV 1a. This was isolated from the four main regions sampled in this study. BVDV 1b was isolated on two premises; one in the southwest and one in the southeast of England. Whilst these findings are largely in agreement with those of Vilcek et al.[10, 19], the data in this paper also identifies a likely set of type 1d and 1e isolates supporting the conclusions of Strong et al.[20] that the genetic diversity of BVD virus in the UK has increased since 1999. There is potential that the emergence of BVDV sub-genotypes not previously reported in the UK could be a symptom of increased trade of cattle/cattle products.

No type 2 BVDV isolates were identified during the investigations in this study and this is re-assuring since the vaccines currently available in the UK only confer good cross protection against type 1 strains. Isolated occurrences of BVDV 2 have been reported in the UK[3640] but these are normally linked to cattle movements from abroad and do not yet appear to have become endemic.

Current literature suggests that there is a continuing increase in the diversity of the BVDV 1 sub-genotype[3, 1115]. When this diversity is considered alongside the data presented by Strong et al.[20] and also in this manuscript, the emergence of isolates previously unreported in the UK leads one to question whether this has consequences for UK BVDV control programmes. The vaccines currently available in the UK which offer foetal protection against BVDV were first developed 13 and 17 years ago[41, 42] and, at the time, were shown to provide good cross protection against the circulating type 1 strains. Further studies may be necessary to demonstrate that cross protection against BVDV 1 is still extensive considering the identification of increasing numbers of circulating BVDV 1 sub-genotypes. Of course, phylogenetic analysis may not provide an accurate indication of antigenic cross reactivity. Challenge studies, antigenic testing and/or investigation of the E2 and NS2-3 sequences (and predicted structures) may be required to ascertain whether the available vaccines are still cross protective against all BVDV 1 sub-genotypes. In addition to the increasing diversity of the BVDV 1 sub-genotype, the emergence of “atypical” pestivirus strains within the UK could become a particular concern. Originating from Brazil, the “HoBi” like viruses have been reported to cause severe disease and pathology in cattle in Northern Italy[43, 44]. It is thought that “HoBi” like viruses were most likely introduced into Europe via live vaccines manufactured with contaminated bovine serum[43]. Thought should be given to the threat that atypical pestiviruses such as “HoBi” represent to the UK herd since, in a national herd with limited immunity to strains outside of BVDV type 1, the effects could be severe. To this end, continued surveillance of BVDV outbreaks should be undertaken, as well as screening of products containing imported bovine serum and considerable care with the licensing and use of live vaccines.

The phylogenetic analysis performed on isolates presented in this manuscript utilises the Neighbour Joining method[25] which has been most commonly used in the published literature for analysis and classification of BVDV. The bootstrap figures for the main branches in both Figures 2 and3 appear to offer low statistical support (<70%) for the tree structure. Essentially, this is a consequence of the large number of isolates with high identity contained in the trees and the resultant difficulty in consistently defining the order within each branch. The correlation between the groupings illustrated in both phylogenetic trees and defined by nucleotide sequence analysis provides good support for the tree structures in Figures 2 and3.

From the dataset investigated, it is evident from both Figure 2 and Table 2 that investigation of the 5′UTR alone would cause some confusion regarding the sub-genotypes present since it was unclear in some instances whether isolates were BVDV 1a, 1c or 1j from sole analysis of the 5′UTR. However, analysis of the Npro region confirmed that where this confusion occurred, the isolates were actually BVDV 1a. Furthermore, one could argue that the BVDV 1i branch highlighted in Figure 2 ought to be included in the bracket demarcating BVDV 1a (along with the 1c and 1j reference strains). Again, Npro analysis was able to provide a more convincing definition of the isolates on this branch confirming them as likely BVDV 1i isolates. Becher et al.[5] has previously noted the limitations of using the 5′UTR and the increased statistical support afforded to phylogenetic trees generated based upon the Npro region. It is however worth re-iterating this point again since many of the most recent publications sub-genotyping BVDV isolates still concentrate on the 5′UTR with little additional information from other genomic regions. As a result, ideally, at least two regions of the BVDV genome should be analysed and agreement sought between them in order to define the isolates investigated with some certainty. Several studies have explored the use of three or more regions of the BVDV genome for sub-genotyping[45, 46]. The methods employed by Liu et al.[46] produced a reliable definition of the isolates investigated and these might be employed to produce, unequivocally, a set of reference strains for future phylogenetic analyses. The analyses conducted in the current manuscript indicated that the confusion in sub-genotype classification could be further enhanced depending upon which reference strains were included. As a result, it would seem sensible to suggest that an agreed database of reference strains (agreed by the pestivirus community) is constructed to enable reliable and consistent analysis that is comparable between publications. A regulated reference strain library would also help avoid the confusion caused by the assignment of BVDV 1l to multiple unique BVDV isolates by Jackova et al.[11] and Kadir et al.[12]. Peterhans et al.[3] has subsequently suggested that the Turkish BVDV 1l strains be re-classified to BVDV 1p. However, Xue et al.[14] have already described a Chinese BVDV 1p sub-genotype. These matters require some clarification to ensure clarity with future sub-genotype definition.

In addition to the more traditional use of phylogenetics to investigate and type viral isolates present in a region or country, the data presented in this manuscript links viral isolates to the farms investigated and explores the potential routes of viral transmission. The farm level data collected enabled a detailed analysis of animal movements to be performed between premises which could, in some cases, be related directly to the phylogeny of the viruses isolated. In effect we are less concerned with primary branching and definition of viral type but instead interested in the clade formation at the end of each branch. Observing clustering at branch ends on the phylogenetic trees enables rapid identification of isolates with a high identity and thus high chance of a common infection route. Three different situations were described in the results:

Same farm, same isolate: Where more than one PI is isolated on the same farm, typing of those isolates commonly shows that the virus infecting each PI will share close to 100% nucleotide identity in both the 5′UTR and Npro regions. The phenomenon of clustering of strains on the same farm is not new and has been noted previously[47]. As yet, the reason why such homology is observed within a herd, yet variation is seen between herds is unknown. Uninhibited viral replication in an infected foetus prior to the development of immunocompetance would seem an obvious place for genetic variation to occur, but it seems that this is not the case[48].

Same farm, different isolate: Where the one farm demonstrates two (or more) markedly different groupings upon phylogenetic analysis, it is likely that more than one infection event has occurred (perhaps two PIs bought in from different locations) producing clusters of PIs with two unique viral isolates.

Different farm, same isolate: Ridpath et al.[49] used sequence and phylogenetic analysis to confirm that the same strain of BVDV2 was responsible for multiple outbreaks between 1993 and 1995, however the routes of transmission were not investigated. For the first time, we show that phylogenetic and sequence analysis of BVDV isolates can be applied to determine factors influencing the epidemiology and transmission of the disease. Whilst somewhat limited, the network of farms involved in this research has enabled novel analysis of animal movement data alongside BVDV phylogeny. Where viral isolates from different farms were seen to cluster closely in both phylogenetic trees (with 99–100% nucleotide identity in both the 5'UTR and Npro regions), we identified links between the farms related to animal movements as demonstrated by Farms 33 and 34. There is a high frequency of animal movements within the network that includes these two premises and it is highly likely that either direct animal movements between these farms, or indirect movements through Farm 32 have resulted in the transmission of infection between Farms 33 and 34. There is perhaps an outstanding question regarding the status of Farm 32, especially if it is considered likely that infection has moved through this unit between Farms 33 and 34. At the point infection was detected on Farms 33 and 34, Farm 32 was being used as a facility to store heifers before they went on to the unit where they were destined to milk. If these animals were kept as separate management groups when on Farm 32, it is conceivable that the infection did not spread amongst the homebred milking herd and youngstock. This highlights the importance of taking representative samples from each management group when performing youngstock “spot tests” for BVD sero-surveillance[50]. Given the potential failure to diagnose infection on Farm 32, the other farms integrated into the movement network with Farms 33 and 34 were at an increased risk of infection before the PI animals were removed and biosecurity advice given. Had infection spread on Farm 32, a particularly high risk period would have been in 2008 when the unit was closed and animals transferred to other premises in the network. Animal movements from infected farms do not always result in infection of the destination premises but, where this does occur, we have shown that phylogenetic analysis may be used to determine the infection source. This may have implications in legal cases where typing of the viral isolates identified on the premises involved may prove or disprove that one farm infected another. Showing cattle has also been identified in this manuscript as a potential source of BVDV transmission. Expansion of an analysis of the sort described here with detailed animal movement and farm location data could allow for quantification and ranking of risk factors associated with BVDV transmission.

In conclusion, the work in this study supports the findings of Strong et al.[20] and has identified an additional isolate (BVDV 1d) that has not previously been reported in the UK. The apparently increasing diversity of strains in the UK may have consequences for vaccine efficacy and further work is required to determine whether this is the case. Whilst the strains identified in this research were all members of the BVDV 1 genotype, the increasing diversity highlights the need for continued and thorough surveillance in order to rapidly detect and react to incursions of new and atypical pestiviruses.

We have demonstrated that in depth investigations of animal movements and contacts between farms combined with phylogenetic analysis of BVDV can produce a “papertrail” of infection. This may have potential when utilised on a wider scale to quantify the risk factors associated with BVDV transmission between premises. These investigations may also have implications in legal cases where the origin of infection is under dispute.

Finally, the analysis of the data presented in this paper has demonstrated areas of pestivirus classification that require clarification. Particular issues arising from this work have been highlighted which would benefit from set guidelines in order to standardise future BVDV phylogenetic analyses. This is particularly important where these processes are used to determine the existence of new pestivirus strains/subgroups. This research, as a consequence, may also develop discussion within the pestivirus community in order to develop agreed guidelines to attain consistent investigation and classification for future epidemiological studies and newly identified isolates.

References

  1. Meyers G, Thiel HJ: Molecular characterization of pestiviruses. Adv Virus Res. 1996, 47: 53-118.

    Article  CAS  PubMed  Google Scholar 

  2. Schirrmeier H, Strebelow G, Depner K, Hoffmann B, Beer M: Genetic and antigenic characterization of an atypical pestivirus isolate, a putative member of a novel pestivirus species. J Gen Virol. 2004, 85: 3647-3652. 10.1099/vir.0.80238-0.

    Article  CAS  PubMed  Google Scholar 

  3. Peterhans E, Bachofen C, Stalder H, Schweizer M: Cytopathic bovine viral diarrhea viruses (BVDV): emerging pestiviruses doomed to extinction. Vet Res. 2010, 41: 44-10.1051/vetres/2010016.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Avalos-Ramirez R, Orlich M, Thiel HJ, Becher P: Evidence for the presence of two novel pestivirus species. Virology. 2001, 286: 456-465. 10.1006/viro.2001.1001.

    Article  CAS  PubMed  Google Scholar 

  5. Becher P, Orlich M, Shannon AD, Horner G, Konig M, Thiel HJ: Phylogenetic analysis of pestiviruses from domestic and wild ruminants. J Gen Virol. 1997, 78: 1357-1366.

    Article  CAS  PubMed  Google Scholar 

  6. Vilcek S, Ridpath JF, Van Campen H, Cavender JL, Warg J: Characterization of a novel pestivirus originating from a pronghorn antelope. Virus Res. 2005, 108: 187-193. 10.1016/j.virusres.2004.09.010.

    Article  CAS  PubMed  Google Scholar 

  7. Kirkland PD, Frost MJ, Finlaison DS, King KR, Ridpath JF, Gu X: Identification of a novel virus in pigs–Bungowannah virus: a possible new species of pestivirus. Virus Res. 2007, 129: 26-34. 10.1016/j.virusres.2007.05.002.

    Article  CAS  PubMed  Google Scholar 

  8. Ridpath JF, Bolin SR, Dubovi EJ: Segregation of bovine viral diarrhea virus into genotypes. Virology. 1994, 205: 66-74. 10.1006/viro.1994.1620.

    Article  CAS  PubMed  Google Scholar 

  9. Vilcek S, Durkovic B, Kolesarova M, Greiser-Wilke I, Paton D: Genetic diversity of international bovine viral diarrhoea virus (BVDV) isolates: identification of a new BVDV-1 genetic group. Vet Res. 2004, 35: 609-615. 10.1051/vetres:2004036.

    Article  CAS  PubMed  Google Scholar 

  10. Vilcek S, Paton DJ, Durkovic B, Strojny L, Ibata G, Moussa A, Loitsch A, Rossmanith W, Vega S, Scicluna MT, Paifi V: Bovine viral diarrhoea virus genotype 1 can be separated into at least eleven genetic groups. Arch Virol. 2001, 146: 99-115. 10.1007/s007050170194.

    Article  CAS  PubMed  Google Scholar 

  11. Jackova A, Novackova M, Pelletier C, Audeval C, Gueneau E, Haffar A, Petit E, Rehby L, Vilcek S: The extended genetic diversity of BVDV-1: typing of BVDV isolates from France. Vet Res Commun. 2008, 32: 7-11. 10.1007/s11259-007-9012-z.

    Article  CAS  PubMed  Google Scholar 

  12. Kadir Y, Christine F, Barbara BW, Zeki Y, Feray A, Aykut O, Ibrahim B, Sibilina Cedillo R, Heinz-Jurgen T, Matthias K: Genetic heterogeneity of bovine viral diarrhoea virus (BVDV) isolates from Turkey: identification of a new subgroup in BVDV-1. Vet Microbiol. 2008, 130: 258-267. 10.1016/j.vetmic.2008.01.016.

    Article  CAS  PubMed  Google Scholar 

  13. Xu X, Zhang Q, Yu X, Liang L, Xiao C, Xiang H, Tu C: Sequencing and comparative analysis of a pig bovine viral diarrhea virus genome. Virus Res. 2006, 122: 164-170. 10.1016/j.virusres.2006.05.005.

    Article  CAS  PubMed  Google Scholar 

  14. Xue F, Zhu YM, Li J, Zhu LC, Ren XG, Feng JK, Shi HF, Gao YR: Genotyping of bovine viral diarrhea viruses from cattle in China between 2005 and 2008. Vet Microbiol. 2010, 143: 379-383. 10.1016/j.vetmic.2009.11.010.

    Article  CAS  PubMed  Google Scholar 

  15. Nagai M, Hayashi M, Itou M, Fukutomi T, Akashi H, Kida H, Sakoda Y: Identification of new genetic subtypes of bovine viral diarrhea virus genotype 1 isolated in Japan. Virus Genes. 2008, 36: 135-139. 10.1007/s11262-007-0190-0.

    Article  CAS  PubMed  Google Scholar 

  16. Becher P, Orlich M, Kosmidou A, Konig M, Baroth M, Thiel HJ: Genetic diversity of pestiviruses: identification of novel groups and implications for classification. Virology. 1999, 262: 64-71. 10.1006/viro.1999.9872.

    Article  CAS  PubMed  Google Scholar 

  17. Tajima M, Frey HR, Yamato O, Maede Y, Moennig V, Scholz H, Greiser-Wilke I: Prevalence of genotypes 1 and 2 of bovine viral diarrhea virus in Lower Saxony, Germany. Virus Res. 2001, 76: 31-42. 10.1016/S0168-1702(01)00244-1.

    Article  CAS  PubMed  Google Scholar 

  18. Flores EF, Ridpath JF, Weiblen R, Vogel FS, Gil LH: Phylogenetic analysis of Brazilian bovine viral diarrhea virus type 2 (BVDV-2) isolates: evidence for a subgenotype within BVDV-2. Virus Res. 2002, 87: 51-60. 10.1016/S0168-1702(02)00080-1.

    Article  CAS  PubMed  Google Scholar 

  19. Vilcek S, Drew TW, McGoldrick A, Paton DJ: Genetic typing of bovine pestiviruses from England and Wales. Vet Microbiol. 1999, 69: 227-237. 10.1016/S0378-1135(99)00111-X.

    Article  CAS  PubMed  Google Scholar 

  20. Strong R, Errington J, Cook R, Ross-Smith N, Wakeley P, Steinbach F: Increased phylogenetic diversity of bovine viral diarrhoea virus type 1 isolates in England and Wales since 2001. Vet Microbiol. 2012, 162: 315-320.

    Article  PubMed  Google Scholar 

  21. Booth RE, Brownlie J: Establishing a pilot bovine viral diarrhoea virus eradication scheme in Somerset. Vet Rec. 2012, 170: 73-10.1136/vr.100191.

    Article  CAS  PubMed  Google Scholar 

  22. Vilcek S, Herring AJ, Herring JA, Nettleton PF, Lowings JP, Paton DJ: Pestiviruses isolated from pigs, cattle and sheep can be allocated into at least three genogroups using polymerase chain reaction and restriction endonuclease analysis. Arch Virol. 1994, 136: 309-323. 10.1007/BF01321060.

    Article  CAS  PubMed  Google Scholar 

  23. Vilcek S, Nettleton PF, Paton DJ, Belak S: Molecular characterization of ovine pestiviruses. J Gen Virol. 1997, 78: 725-735.

    Article  CAS  PubMed  Google Scholar 

  24. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.

    CAS  PubMed  Google Scholar 

  26. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.

    Article  CAS  PubMed  Google Scholar 

  27. Hillis DM, Bull JJ: An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol. 1993, 42: 182-192.

    Article  Google Scholar 

  28. Colett MS, Larson R, Gold C, Strick D, Anderson DK, Purchio AF: Molecular cloning and nucleotide sequence of the pestivirus bovine viral diarrhea virus. Virology. 1988, 165: 191-199. 10.1016/0042-6822(88)90672-1.

    Article  CAS  PubMed  Google Scholar 

  29. Deng R, Brock KV: 5' and 3' untranslated regions of pestivirus genome: primary and secondary structure analyses. Nucleic Acids Res. 1993, 21: 1949-1957. 10.1093/nar/21.8.1949.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. De Moerlooze L, Lecomte C, Brown-Shimmer S, Schmetz D, Guiot C, Vandenbergh D, Allaer D, Rossius M, Chappuis G, Dina D, Renard A, Martial JA: Nucleotide sequence of the bovine viral diarrhoea virus Osloss strain: comparison with related viruses and identification of specific DNA probes in the 5' untranslated region. J Gen Virol. 1993, 74: 1433-1438. 10.1099/0022-1317-74-7-1433.

    Article  CAS  PubMed  Google Scholar 

  31. Ridpath JF, Fulton RW, Kirkland PD, Neill JD: Prevalence and antigenic differences observed between Bovine viral diarrhea virus subgenotypes isolated from cattle in Australia and feedlots in the southwestern United States. J Vet Diagn Invest. 2010, 22: 184-191. 10.1177/104063871002200203.

    Article  PubMed  Google Scholar 

  32. Nagai M, Sakoda Y, Mori M, Hayashi M, Kida H, Akashi H: Insertion of cellular sequence and RNA recombination in the structural protein coding region of cytopathogenic bovine viral diarrhoea virus. J Gen Virol. 2003, 84: 447-452. 10.1099/vir.0.18773-0.

    Article  CAS  PubMed  Google Scholar 

  33. Baule C, van Vuuren M, Lowings JP, Belak S: Genetic heterogeneity of bovine viral diarrhoea viruses isolated in Southern Africa. Virus Res. 1997, 52: 205-220. 10.1016/S0168-1702(97)00119-6.

    Article  CAS  PubMed  Google Scholar 

  34. Schweizer M, Peterhans E: Oxidative stress in cells infected with bovine viral diarrhoea virus: a crucial step in the induction of apoptosis. J Gen Virol. 1999, 80: 1147-1155.

    Article  CAS  PubMed  Google Scholar 

  35. Ridpath JF, Bolin SR: The genomic sequence of a virulent bovine viral diarrhea virus (BVDV) from the type 2 genotype: detection of a large genomic insertion in a noncytopathic BVDV. Virology. 1995, 212: 39-46. 10.1006/viro.1995.1451.

    Article  CAS  PubMed  Google Scholar 

  36. Courtenay AE, Henderson RG, Cranwell MP, Sandvik T: BVD virus type 2 infection and severe clinical disease in a dairy herd. Vet Rec. 2007, 160: 706-707.

    Article  CAS  PubMed  Google Scholar 

  37. Drew TW, Sandvik T, Wakeley P, Jones T, Howard P: BVD virus genotype 2 detected in British cattle. Vet Rec. 2002, 151: 551-

    CAS  PubMed  Google Scholar 

  38. Wakeley PR, Turner JL, Ibata G, King DP, Sandvik T, Howard P, Drew TW: Characterisation of a type 2 bovine viral diarrhoea virus isolated from cattle in the UK. Vet Microbiol. 2004, 102: 19-24. 10.1016/j.vetmic.2004.05.005.

    Article  CAS  PubMed  Google Scholar 

  39. Cranwell MP, Jones JR, Wakeley PR: BVD virus type 2 in British cattle. Vet Rec. 2005, 156: 257-258.

    CAS  PubMed  Google Scholar 

  40. Drew T, David G, Steinbach F: Occurrence of a second strain of BVDV-2 in the UK. Vet Rec. 2007, 160: 810-

    Article  PubMed  Google Scholar 

  41. Brownlie J, Clarke MC, Hooper LB, Bell GD: Protection of the bovine fetus from bovine viral diarrhoea virus by means of a new inactivated vaccine. Vet Rec. 1995, 137: 58-62. 10.1136/vr.137.3.58.

    Article  CAS  PubMed  Google Scholar 

  42. Patel JR, Williams JK: Protection by vaccination from in utero infection by BVD virus. 1998, Sydney, Australia: Proceedings of the XX World Association for Buiatrics Congress

    Google Scholar 

  43. Decaro N, Lucente MS, Mari V, Cirone F, Cordioli P, Camero M, Sciarretta R, Losurdo M, Lorusso E, Buonavoglia C: Atypical pestivirus and severe respiratory disease in calves, Europe. Emerg Infect Dis. 2011, 17: 1549-1552.

    PubMed Central  PubMed  Google Scholar 

  44. Decaro N, Mari V, Lucente MS, Sciarretta R, Moreno A, Armenise C, Losurdo M, Camero M, Lorusso E, Cordioli P, Buonavoglia C: Experimental infection of cattle, sheep and pigs with ‘Hobi’-like pestivirus. Vet Microbiol. 2012, 155: 165-171. 10.1016/j.vetmic.2011.08.030.

    Article  PubMed  Google Scholar 

  45. Nagai M, Hayashi M, Sugita S, Sakoda Y, Mori M, Murakami T, Ozawa T, Yamada N, Akashi H: Phylogenetic analysis of bovine viral diarrhea viruses using five different genetic regions. Virus Res. 2004, 99: 103-113. 10.1016/j.virusres.2003.10.006.

    Article  CAS  PubMed  Google Scholar 

  46. Liu L, Xia H, Wahlberg N, Belak S, Baule C: Phylogeny, classification and evolutionary insights into pestiviruses. Virology. 2009, 385: 351-357. 10.1016/j.virol.2008.12.004.

    Article  CAS  PubMed  Google Scholar 

  47. Vilcek S, Alenius S, Paton DJ, Mittelholzer C, Belak S: Genetic clustering of bovine viral diarrhoea viruses in cattle farms: genetic identification and analysis of viruses directly from cattle sera. Vet J. 1999, 158: 33-38. 10.1053/tvjl.1999.0363.

    Article  CAS  PubMed  Google Scholar 

  48. Stokstad M, Brownlie J, Collins ME: Analysis of variation of bovine viral diarrhoea virus E2 sequence following transplacental infection of cattle. Vet Microbiol. 2004, 102: 141-145. 10.1016/j.vetmic.2004.06.015.

    Article  CAS  PubMed  Google Scholar 

  49. Ridpath JF, Neill JD, Vilcek S, Dubovi EJ, Carman S: Multiple outbreaks of severe acute BVDV in North America occurring between 1993 and 1995 linked to the same BVDV2 strain. Vet Microbiol. 2006, 114: 196-204. 10.1016/j.vetmic.2005.11.059.

    Article  PubMed  Google Scholar 

  50. Houe H, Baker JC, Maes RK, Ruegg PL, Lloyd JW: Application of antibody titers against bovine viral diarrhea virus (BVDV) as a measure to detect herds with cattle persistently infected with BVDV. J Vet Diagn Invest. 1995, 7: 327-332. 10.1177/104063879500700305.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank DEFRA for funding this work as part of project grant SE0777. In addition, our thanks must go to Dr Antonia Ganser (SAC Epidemiology Division, Inverness) for supplying strains collecting PI samples from herds in Scotland and Northern England, to George Caldow (SAC, St. Boswells) for the strains investigated from East Anglia and Farms 64–66, to Paul Burr (Biobest Laboratories Ltd) for submissions from Farm 56 and also to the individual practitioners who have supplied samples for this work. Our thanks must also go to Julia Ridpath for her technical advice on the PCR procedures, phylogenetic analysis, and also Angel Ortiz Pelaez for his help establishing the analysis of the animal movement data. Finally, we would like to thank Julia Ridpath, Peter Kirkland and Stefan Vilcek for not only supplying sequences for several “reference” strains utilised in this publication, but also for their support and advice about the broader issues of pestivirus taxonomy that have informed our views in this paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Richard E Booth.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RB conducted the sample collection from herds 1–42, performed RT-PCR followed by sequence and phylogenetic analysis on samples from all farms involved. CT and LA gave extensive advice on the technical aspects of the phylogenetic and sequence analysis. CT, LA, GG and JB contributed to the drafting of the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Booth, R.E., Thomas, C.J., El-Attar, L.M. et al. A phylogenetic analysis of Bovine Viral Diarrhoea Virus (BVDV) isolates from six different regions of the UK and links to animal movement data. Vet Res 44, 43 (2013). https://doi.org/10.1186/1297-9716-44-43

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1297-9716-44-43

Keywords