Modelling Marek's Disease Virus (MDV) infection: parameter estimates for mortality rate and infectiousness
- Katherine E Atkins^{1, 4}Email author,
- Andrew F Read^{2},
- Nicholas J Savill^{1},
- Katrin G Renz^{3},
- Stephen W Walkden-Brown^{3} and
- Mark EJ Woolhouse^{1}
DOI: 10.1186/1746-6148-7-70
© Atkins et al; licensee BioMed Central Ltd. 2011
Received: 17 November 2010
Accepted: 11 November 2011
Published: 11 November 2011
Abstract
Background
Marek's disease virus (MDV) is an economically important oncogenic herpesvirus of poultry. Since the 1960s, increasingly virulent strains have caused continued poultry industry production losses worldwide. To understand the mechanisms of this virulence evolution and to evaluate the epidemiological consequences of putative control strategies, it is imperative to understand how virulence is defined and how this correlates with host mortality and infectiousness during MDV infection. We present a mathematical approach to quantify key epidemiological parameters. Host lifespan, virus latent periods and host viral shedding rates were estimated for unvaccinated and vaccinated birds, infected with one of three MDV strains. The strains had previously been pathotyped to assign virulence scores according to pathogenicity of strains in hosts.
Results
Our analyses show that strains of higher virulence have a higher viral shedding rate, and more rapidly kill hosts. Vaccination enhances host life expectancy but does not significantly reduce the shedding rate of the virus. While the primary latent period of the virus does not vary with challenge strain nor vaccine treatment of host, the time until the maximum viral shedding rate is increased with vaccination.
Conclusions
Our approach provides the tools necessary for a formal analysis of the evolution of virulence in MDV, and potentially simpler and cheaper approaches to comparing the virulence of MDV strains.
Background
Marek's Disease Virus (MDV) is an oncogenic poultry herpesvirus of considerable economic importance to the poultry industry. Virus strains have become increasingly virulent since the 1960s [1, 2]. The causes of this viral evolution are unclear [3]. Quantification of patterns of viral shedding and virus-induced host mortality are necessary for a rigorous understanding of the epidemiology of a disease, not least to identify increases in virulence. Here we develop methods to do this.
MDV is an airborne pathogen with infection occurring via inhalation [4]. Virus shedding occurs by infected feather follicle epithelium [5]. The resulting dust and dander from dead stratified cells and moulted feathers can then remain in the environment and act as a reservoir for chicken infection. Clinical signs are varied and result in significant morbidity and mortality depending on host genetic susceptibility and virulence of the MDV strain [6]. Symptoms include polyneuritis (an enlargement of multiple peripheral nerves), visceral lymphoma (tumours affecting organs such as the heart, liver, spleen etc.), acute transient paralysis, immunosuppression, brain oedema and acute rash. There has been a change in the types of clinical signs since the disease was first noted [7, 8], when chronic polyneuritis was the only sign. Since then, the list of clinical signs described above expanded gradually over the decades [2].
Almost all industrialised countries have experienced MD losses in their poultry industry and a crude estimate of the cost of Marek's disease is said to be in the range of US $1-2 billion annually [7]. Control of Marek's disease is predominantly via vaccination of chickens. Crucially the MDV vaccine was the first vaccine to be developed against any cancer [9]. Three types of vaccine have been developed for use against MD. These are herpesvirus of turkeys (HVT), non-pathogenic serotype 2 MDV and non-pathogenic serotype 1 [10]. These vaccines have been used in different sequences in different countries and the vaccine types have also been commonly combined in bivalent or trivalent vaccines.
Definitions of virulence are numerous and varied [11]. Microbiologists equate virulence with the notions of both infectivity and severity of disease, whereas evolutionary biologists focus on evolutionary fitness of either the pathogen or the host [12]. Zoologists tend to focus on host fitness [13, 14] with most mathematical models describing this specifically as host mortality [14, 15].
There have been various attempts to define virulence in the context of MDV [16]. For example, Witter [1] relates the percentage of HVT or Bivalent (HVT + serotype 2) vaccinated chickens, that when infected with a particular strain develop gross lesions or die of Marek's disease within eight weeks to the figure in unvaccinated chickens challenged with the same strain. This percentage score uses the protective effect of vaccines as its metric for virulence. In addition, [17] use neurovirulence as a tool for pathotyping MDV strains, noting that many of the very virulent strains are synonymous with high levels of neuropathology. This method circumvents the need for both vaccine-based definitions and therefore high numbers of birds, and also may cut the experiment time. In the first pathotyping regime, there are four recognised MDV serotype 1 pathotypes, each occupying part of the 'continuum of virulence' [1] defined as: mild, m; virulent, v; very virulent, vv; very virulent plus, vv+. This grouping correlates very well with the second pathotyping regime [18].
Control strategies for MDV require an understanding of the epidemiology of the disease, in particular how virulence relates to key parameters such as viral shedding rates and duration of infectious period. There have been efforts to find correlates with virulence of MDV isolates, most notably with viral load (virus within a bird tissue). This has been achieved by cell culture techniques [19] and intra-cellular detection during the first 10 days of infection [20]. However, the relationship between viral load (or replication) and virulence has not been convincingly tested [16]. PCR testing methods have been developed in order to quantify viral loads [20–26], which have enabled the viral loads in shed dust to be directly measured [27, 28]. The infectious period is defined as the total time in which an infected bird sheds virus [29]. This is determined by the clearance rate and the lifespan of the individual host. Since there is no recovery from MDV infection, the infectious period is defined by a total of four parameters: the disease-induced mortality, the latent period (the time between infection and infectiousness of a host) and two other non-disease associated parameters, namely non-disease induced mortality and the maximum lifespan of an individual (the farm slaughter time in the case of broilers).
Here we present methods for parameter estimation for MDV isolates, allowing a formal quantification of infected host lifespan, viral shedding rates and viral latent periods. The parameters determined from the analysis allow comparison between isolates of different virulence in both unvaccinated and vaccinated birds. It is anticipated that the parameter estimates will be used for a better understanding of the pathology and aetiology of the disease itself and as a platform for investigation into the causes of virulence evolution.
Results
In the first section, we used survival analysis with a Weibull mortality function to build a statistical model for the lifespan of birds infected with MDV. In the second section, we developed a dynamic model to simulate the shedding of virus by a group of birds and used Bayesian techniques to estimate the four key infectiousness parameters: primary and secondary latent periods and viral shedding rates.
Mortality
Survival analysis: Estimated coefficients of the covariates in the survival analysis.
Symbol | Coefficient | Value (95% CI via McMC) | z | p |
---|---|---|---|---|
r | Shape Parameter | 4.18 (3.38,4.99) | -11.60 | < 0.001 |
β _{0} | Intercept | 4.54 (4.23,4.87) | 26.08 | < 0.001 |
β _{1} | Virulence Score | -0.53 (-1.01,-0.06) | -2.04 | 0.040 |
β _{2} | HVT vaccine | 0.44 (0.24,0.66) | 3.95 | < 0.001 |
β _{3} | Bivalent vaccine | 0.35 (0.19,0.52) | 4.28 | < 0.001 |
Shedding
Viral shedding unvaccinated: median and 95% credible intervals for the posterior distribution of the four parameters for both replicates of the sham vaccinated factorial study for different virulence scores, v.
VIRULENCE SCORE | ||||||
---|---|---|---|---|---|---|
PARAMETER | 16.5 | 36 | 46 | |||
Rep 1 | Rep 2 | Rep 1 | Rep 2 | Rep 1 | Rep 2 | |
1 Latent Period (days) | 6 (2,6) | 4 (0,6) | 4 (0,6) | 4 (0,6) | 3 (0,6) | 3 (0,6) |
2 Latent Period (days) | 10 (7,14) | 9 (6,13) | 9 (6,13) | 9 (6,13) | 9 (6,13) | 9 (5,13) |
1 Shedding Rate (logVCN/mg dust) | 4.87 (3.79,5.25) | 3.37 (2.79,4.05) | 3.86 (3.30,4.56) | 3.28 (2.71,3.95) | 3.20 (2.64,3.90) | 2.89 (2.34,3.59) |
2 Shedding Rate (logVCN/mg dust) | 6.99 (6.72,7.21) | 6.89 (6.72,7.07) | 7.34 (7.17,7.51) | 7.58 (7.41,7.76) | 7.69 (7.51,7.85) | 7.52 (7.32,7.72) |
Viral shedding HVT vaccination: median and 95% credible intervals for the posterior distribution of the four parameters for both replicates of the HVT vaccinated factorial study for different virulence scores, v.
VIRULENCE SCORE | ||||||
---|---|---|---|---|---|---|
PARAMETER | 16.5 | 36 | 46 | |||
Rep 1 | Rep 2 | Rep 1 | Rep 2 | Rep 1 | Rep 2 | |
1 Latent Period (days) | 5 (0,6) | 2 (0,6) | 6 (4,6) | 6 (2,6) | 3 (0,6) | 3 (0,6) |
2 Latent Period (days) | 12 (8,17) | 10 (7,13) | 11 (8,14) | 12 (8,16) | 10 (6,13) | 9 (6,13) |
1 Shedding Rate (logVCN/mg dust) | 3.79 (3.34,4.20) | 2.44 (1.93,3.02) | 4.58 (4.19,4.92) | 4.56 (4.13,4.91) | 2.70 (2.16,3.42) | 4.15 (3.58,4.82) |
2 Shedding Rate (logVCN/mg dust) | 6.83 (6.61,7.04) | 6.35 (6.18,6.52) | 7.49 (7.28,7.70) | 7.53 (7.31,7.75) | 7.50 (7.32,7.67) | 7.29 (7.11,7.46) |
Viral shedding bivalent vaccination: median and 95% credible intervals for the posterior distribution of the four parameters for both replicates of the Bivalent vaccinated factorial study for different virulence scores, v.
VIRULENCE SCORE | ||||||
---|---|---|---|---|---|---|
PARAMETER | 16.5 | 36 | 46 | |||
Rep 1 | Rep 2 | Rep 1 | Rep 2 | Rep 1 | Rep 2 | |
1 Latent Period (days) | 4 (0,6) | 2 (0,5) | 3 (0,6) | 3 (0,6) | 3 (0,6) | 3 (0,6) |
2 Latent Period (days) | 14 (9,19) | 23 (18,27) | 15 (10,19) | 9 (6,13) | 10 (6,13) | 10 (6,13) |
1 Shedding Rate (logVCN/mg dust) | 4.87 (3.27,5.19) | 3.37 (3.27,3.83) | 3.86 (2.92,4.70) | 3.28 (2.93,4.14) | 3.20 (2.94,4.18) | 2.89 (2.49,4.69) |
2 Shedding Rate (logVCN/mg dust) | 6.23 (6.02,6.45) | 5.85 (5.61,6.09) | 7.16 (6.94,7.37) | 6.95 (6.78,7.13) | 6.96 (6.78,7.12) | 7.28 (7.11,7.45) |
Secondary latent period: linear regression results for effect of transformed virulence score and vaccine treatment on the secondary latent period (days).
Estimate | Std. Error | t value | Pr(> |t|) | |
---|---|---|---|---|
Intercept | 10.16 | 4.30 | 2.36 | 0.04 |
Virulence Score | -1.64 | 6.98 | -0.24 | 0:82 |
VaccHVT | 2.60 | 6.10 | 0.43 | 0:68 |
VaccBiv | 19.29 | 6.10 | 3.16 | < 0.01 |
Virulence Score:VaccBiv | -24.82 | 9.88 | -2.51 | 0.03 |
Virulence Score:VaccHVT | -1.82 | 9.88 | -0.18 | 0.86 |
Secondary viral shedding: linear regression results for effect of transformed virulence score and vaccine treatment on the second viral shedding rate (VCN/mg dust).
Estimate | Std. Error | t value | Pr(> |t|) | |
---|---|---|---|---|
Intercept | -3.24 × 10^{7} | 1.35 × 10^{7} | -0.24 | 0.03 |
Virulence Score | 9.79 × 10^{7} | 2.19 × 10^{7} | 4.47 | < 0.01 |
VaccHVT | 8.51 × 10^{6} | 1.91 × 10^{7} | 0.45 | 0.66 |
VaccBiv | 1.70 × 10^{7} | 1.91 × 10^{7} | 0.89 | 0.39 |
Virulence Score:VaccBiv | -5.73 × 10^{7} | 3.19 × 10^{7} | -0.19 | 0.09 |
Virulence Score:VaccHVT | -2.37 × 10^{7} | 3.10 × 10^{7} | -0.77 | 0.46 |
Discussion
The survival model fits the nine survival curves well on the whole, with both replicates combined into one fit. The delay which characterises the mortality profile for Marek's disease is captured well by the Weibull model. In most cases, the viral shedding model captures the data well.
The viral shedding model assumes that all birds become initially infected with MDV. While it is difficult to model individual birds based on the resolution of the data available, we do point out that the assumption that all birds become infected precipitates the estimated per bird viral shedding rates being perhaps an over-estimate of the actual rates. On the other hand, the qPCR methods used estimate total viral copies which is not necessarily equivalent to total infectious material and therefore serves as an upper limit of viable infectious DNA.
Our estimate of the primary latent period of between 2 and 6 days is much shorter than earlier estimates of around 13 days [31] but only slightly shorter than more recent estimates based on new PCR techniques which have detected significant quantities of virus in feather tips [32] and feather dander [27, 28] at 7 days post infection. This discrepancy may be due to a number of factors: detectability of virus at low quantities, the sampling time of the experiments and heterogeneities between the experiments. The estimates here also have confidence intervals between 0-6 days which suggest the data are not able to give a great deal of certainty on the exact time of onset of viral shedding. Further analysis might take into account biological knowledge through prior distributions to get a more exact value, and past experimental studies would suggest shedding begins towards the end of the estimated interval [27, 28].
The pathotyping of MDV isolates, similar to that pioneered by Witter and colleagues [1, 16], is a useful tool for comparing the pathogenicity of strains with respect to others in vaccinated birds. However several of the more recent isolated, highly virulent strains have ranks above 80, in some cases 90 on Witter's ranking system. Should more virulent strains emerge, it will be difficult to categorise them via this method, as the scale is truncated at 100. The method of fitting a Weibull distribution proposed in this paper allows for a greater flexibility of pathotyping, in that an isolate may be categorised according to its mean time to kill the host. The subsequent effect of vaccination of the host on this metric can then be established.
Biologically, it may seem more sustainable to pathotype according to the lifespan of an unvaccinated bird, but experimentally it may also be a solution to the costly procedure of having two or three groups of birds, each vaccinated with a different treatment. Even a simple 'time to death' metric can be performed in relatively few birds, as was the case for the myxoma virus in rabbits in Australia. Myxoma strains were pathotyped effectively when the case fatality rate proved experimentally intensive and ineffective in discriminating between similarly virulent strains [33, 17] noted this restriction of current pathotyping regimes, instead favouring a system based on neuropathology to circumvent these methodological and experimental limitations. The analysis presented here differs from both of these past approaches due to its epidemiological focus on virulence. However, we point out that the methodology adopted here only analyses three distinct strains. Therefore, a wider range of strain virulences and places of origin will be required to sufficiently extrapolate the current results to the whole virulence spectrum.
Conclusions
This paper uses a mathematical approach to estimate key epidemiological parameters for Marek's Disease Virus (MDV). We used survival analysis with a Weibull mortality function to build a statistical model for the lifespan of birds infected with MDV. We were able to compare the relative impact of virulence score and vaccine covariates by calculating credible intervals via Bayesian methods. We developed a dynamic model to simulate the shedding of virus by a group of birds and used Bayesian techniques to estimate the four key infectiousness parameters: primary and secondary latent periods and viral shedding rates.
We applied these methods to the data for three MDV strains varying moderately in virulence. From the survival analysis, we found that increasing the virulence score of a virus strain decreases the lifespan of the host (Figure 2). Vaccinating the host with either the first or second generation vaccine increases its lifespan, but neither one to any greater degree (Table 1). From the viral shedding analysis, we concluded primary latent period is not affected by vaccination or virulence score, however, secondary latent period increases with Bivalent vaccine treatment (Table 5). Secondary (long term) viral shedding rate (measured in VCN per mg dust) is higher for strains of higher virulence (Table 5 and Figure 6). Vaccination does not reduce long term viral shedding rate significantly. Together, these results imply that more virulent viruses shed more infectious material into the environment and Bivalent vaccines reduce the infectious period, especially for less virulent strain infections.
The results described here are qualitatively in line with current thinking on MDV pathogenesis, but the new methods presented allow a formal comparison of strains in an epidemiological meaningful way. It is anticipated that this analysis will allow a more quantitative understanding of the aetiology of MDV and the mechanisms governing a functional rise in virulence for a strain. With this analysis it is hoped that a more epidemiologically significant classification of MDV strains will emerge which will eventually provide not only insight into the reasons for evolution to highly virulent strains but also a metric to evaluate future risks and the efficacy of possible control strategies.
Methods
Experimental Data
We look at a single experiment [30] with two response variables measured. Groups of maternal antibody positive IsaBrown layer chicks were inoculated with one of three vaccine treatments: sham (diluent only), HVT (first generation industry vaccine) or Bivalent (second generation industry vaccine) at one day old. At 5 days post vaccination, the birds were challenged with 500 pfu (plaque forming units) of one of three MDV strains: 04CRE, MPF57 and 02LAR. These independently sampled isolates were Australian in origin and have been pathotyped with virulence scores (the percentage of HVT or Bivalent vaccinated chickens infected with a particular strain which develop gross lesions or die of MDV within eight weeks [34]) of 16.5, 36 and 46 respectively. The viruses were all isolated from field outbreaks and used at passage level 4 to 7 in chick kidney cells. The experiment was therefore a 3 × 3 factorial design study, with vaccine type fully cross-factored against virus strain. Each treatment combination was replicated in two separate isolators resulting in 26 or 27 birds per combination for the pathotyping component of the study. The air in each isolator was changed 8-20 times per hour. Further details can be found in [30]. The experiment was approved by the University of New England Animal Ethics Committee with approval number AEC05/076. Due to Animal Ethics regulations, birds thought to be nearing death were euthanised and classified as mortality in the experimental results. The experiment lasted 56 days post infection (dpi) and all birds were humanely killed and examined for MD lesions on the last day. Therefore for the purposes of this analysis, the data were censored after day 55. All birds were given a large dose of MDV (500 pfu) via intra-abdominal injection and it was therefore assumed that all the birds were infected and would shed the virus. Differences between unvaccinated and vaccinated birds in detection of MD have been reported [35], but these have been in studies where infection has been by injection with 50 pfu. Increasing the MDV dose fourfold (from 50 to 200 pfu) significantly increased the detected incidence of MD in unvaccinated birds from 64% to 81% [27].
The two variables recorded were the time to death of a bird (measured in days) and the density of virus within each isolator every week (measured in logVCN per mg of dust in the environment). The virus titres in the dust were determined via quantitative real-time PCR of dust collected from the isolator exhaust ducts, the methods of which have been discussed in separate studies [25]. The data are shown in Figures 1 (mortality), 3, 4 and 5 (shedding).
Parameter Estimates
We assumed that each infected bird starts shedding virus after an initial delay, then after a further delay, the rate at which it sheds virus rises to a long-term stable value. We estimated five key epidemiological parameters: host lifespan, primary latent period (the time until an infected bird becomes infectious), secondary latent period (the delay between an infectious bird first shedding virus and when it sheds virus at its secondary long term rate) and viral shedding rates (both primary and secondary, measured in VCN/mg dust). A survival analysis was undertaken for the mortality data because multiple covariates were needed to find the host lifespan. These covariates were vaccine treatment and virus virulence score. Latent period and rate of viral shedding could not be estimated directly from the data since the response variable (measured weekly in VCN per mg dust) could vary with the number of birds housed together, which changed as birds died. Therefore a dynamic model capturing this variability was used to estimate the two infectiousness parameters.
For ease of notation, log is written to represent log_{10} throughout the paper.
Mortality
where δ_{ i } is zero when the ith observation is censored and unity elsewhere. This function was maximized via the Newton-Raphson algorithm such that $\mathcal{L}\left(\stackrel{\u0303}{\lambda},\stackrel{\u0303}{r}\right)=max\mathcal{L}\left(\lambda ,r\right)$where $\stackrel{\u0303}{\lambda}$ and $\stackrel{\u0303}{r}$ are the maximum likelihood estimates. Note that the virulence score, v, of an isolate is a percentage measure and to be used as an explanatory variable in a regression analysis, it should be transformed such that ${v}_{T}=arcsin\sqrt{0.01v}$. The regression fitted the maximum likelihood estimates for β, however further Bayesian analysis was required to estimate the associated credible intervals when the covariate covariance matrix does not approximate the identity matrix. An McMC framework was set up in WinBUGS [36] to calculate the posterior distributions and the credible intervals of the Weibull survival function.
The prior distributions of each parameter were assumed to be uninformative and were thus taken as uniform, to ensure equivalence to maximum likelihood estimation. The burn-in was set to 22,000 iterations with every 10th sample taken from the subsequent 100,000 iterations [37, 38].
Viral Shedding
The dynamic model tracked the density of MDV over time in each isolator. In this model it is assumed:
1) After MDV infection, there is a delay before virus is first shed [31]. This is termed the primary latent period.
2) After this primary latent period, virus is shed at a constant rate (the primary viral shedding rate, measured in VCN per mg dust) for a set period of time, termed the secondary latent period.
3) Once this secondary latent period is over, virus is shed at a constant rate (the secondary viral shedding rate) until the termination of the experiment.
Additionally, the density of virus (VCN per mg dust) was calculated at the end of each day and any removal of birds was assumed to occur at the start of the day. While all parameters are assumed to be constant across all the birds in a single isolator, they can vary between isolators.
There are therefore four parameters estimated per isolator: the primary/secondary latent periods (in days) and the primary/secondary rates of viral shedding (VCN/mg dust). The primary and secondary latent periods are denoted by T^{1} and T^{2} respectively. Birds shed virus at a rate, a^{1} (VCN per mg dust) from T^{1}+1 to T^{1}+T^{2} days post infection; and at a rate a^{2} (VCN per mg dust) from T^{2}+1 until the end of the experiment (at 8 weeks post infection). The density of MDV in the dynamic model could be compared to the amount of virus recorded in the data (sampled every 7 days). We therefore estimated ${a}_{j}^{1}\left(v\right)$, ${a}_{j}^{2}\left(v\right)$, ${T}_{j}^{1}\left(v\right)$ and ${T}_{j}^{2}\left(v\right)$ for each virulence score v ∈ {16.5, 36, 46} and each vaccine status j ∈ {Sham, HVT, Bivalent}.
The model accounts for the removal of birds given in the mortality data and the number of complete hourly air changes (varied as a parameter α). The daily quantity of dust shed per bird was calculated by fit ting a cubic spline to dust MDV data provided from the same experiment [30].
To find the time point after which the data is assumed to have plateaued, four sets of points were examined: weeks 5, 6, 7 and 8. The differences between the logged data of weeks 7-8 provided the highest probability of both being drawn from a normal distribution (p = 0.8, Anderson-Darling n = 18) and that the distribution has a mean of zero (p = 0.40). The resulting standard deviation of the difference between log-transformed data was estimated to be 0.33 logVCN/mg dust.
Posterior distributions were found by means of McMC realisations by calculating the likelihood of the data given the parameter values, with lognormal errors of σ^{2} = 0.33^{2}/2. The burn-in time was set to 30,000 iterations and the posterior distribution was taken as every 10th sample from the following 90,000 iterations [37, 38].
Declarations
Acknowledgements
We thank BBSRC, Pfizer, the Australian Poultry CRC, R. Howey for supplying Bayesian sampling code and M. Chase-Topping for statistical advice.
Authors’ Affiliations
References
- Witter RL: Increased virulence of Marek's disease virus field isolates. Avian Diseases. 1997, 41: 149-163. 10.2307/1592455.PubMedView ArticleGoogle Scholar
- Osterrieder N, Kamil JP, Schumacher D, Tischer BK, Trapp S: Marek's disease virus: from miasma to model. Nature Reviews Microbiology. 2006, 4: 753-761.View ArticleGoogle Scholar
- Witter RL: The changing landscape of Marek's disease. Avian Pathology. 1998, 27: S46-S53. 10.1080/03079459808419292.View ArticleGoogle Scholar
- Calnek BW, Adldinger HK, Kahn DE: Feather Follicle Epithelium: a Source of Enveloped and Infectious Cell-Free Herpesvirus from Marek's Disease. Avian Diseases. 1970, 14 (2): 219-233. 10.2307/1588466.PubMedView ArticleGoogle Scholar
- Carrozza JH, Fredrickson TN, Prince RP, Luginbuhl RE: Role of desquamated epithelial cells in transmission of Marek's disease. Avian Diseases. 1973, 17: 767-781. 10.2307/1589045.PubMedView ArticleGoogle Scholar
- Baigent SJ, Smith LP, Nair VK, Currie RJW: Vaccinal control of Marek's disease: current challenges, and future strategies to maximize protection. Veterinary Immunology and Immunopathology. 2006, 112: 78-86. 10.1016/j.vetimm.2006.03.014.PubMedView ArticleGoogle Scholar
- Morrow C, Fehler F: Marek's disease: a worldwide problem. Marek's disease: an Evolving Problem. Elsevier Academic Press 2004.Google Scholar
- Nair V, Kung HJ: Marek's disease virus oncogenicity: molecular mechanisms. Marek's Disease: an Evolving Problem. Elsevier Academic Press 2004.Google Scholar
- Churchill A, Payne L, Chubb R: Immunization against Marek's Disease using a live attenuated virus. Nature. 1969, 221 (5182): 744-747. 10.1038/221744a0.PubMedView ArticleGoogle Scholar
- Bublot M, Sharma J: Vaccination against Marek's disease. Marek's Disease: an Evolving Problem. Edited by: Nair V, Davison F. Elsevier Academic Press 2004.Google Scholar
- Dieckmann U: Adaptive dynamics of pathogen-host interactions. Adaptive Dynamics Of Infectious Diseases: In Pursuit of Virulence Management. Cambridge University Press 2002.View ArticleGoogle Scholar
- Read AF: The Evolution of Virulence. Trends in Microbiology. 1994, 2: 73-76. 10.1016/0966-842X(94)90537-1.PubMedView ArticleGoogle Scholar
- Levin S, Pimentel D: Selection of Intermediate Rates of Increase in Parasite-Host Systems. American Naturalist. 1981, 117: 308-315. 10.1086/283708.View ArticleGoogle Scholar
- Anderson RM, May RM: Directly transmitted diseases: Control by vaccination. Science. 1982, 215 (4536): 1053-1060. 10.1126/science.7063839.PubMedView ArticleGoogle Scholar
- Tompkins DM, Dobson AP, Arneberg P, Begon ME, Cattadori IM, Greenman JV, Heesterbeeck JAP, Hudson PJ, Newborn D, Pugliese A, Rizzoli AP, Rosa R, Rosso F, Wilson K: Parasites and Host Population Dynamics. The Ecology of Wildlife Diseases. Oxford University Press 2002.Google Scholar
- Witter RL, Calnek BW, Buscaglia C, Gimeno IM, Schat KA: Classification of Marek's disease viruses according to pathotype: philosophy and methodology. Avian Pathology. 2005, 34 (2): 75-90. 10.1080/03079450500059255.PubMedView ArticleGoogle Scholar
- Gimeno IM, Witter RL, Neumann U: Neuropathotyping: A New System to Classify Marek's Disease Virus. Avian Diseases. 2002, 46 (4): 909-918. 10.1637/0005-2086(2002)046[0909:NANSTC]2.0.CO;2. [http://www.jstor.org/stable/1593315]PubMedView ArticleGoogle Scholar
- Gimeno IM, Witter RL, Hunt HD, Reddy SM, Lee LF, Silva RF: The pp38 gene of Marek's disease virus (MDV) is necessary for cytolytic infection of B cells and maintenance of the transformed state but not for cytolytic infection of the feather follicle epithelium and horizontal spread of MDV. Journal of Virology. 2005, 79 (7): 4545-4549. 10.1128/JVI.79.7.4545-4549.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Calnek BW, Witter RL: Marek's disease - a model for herpesvirus oncology. CRC Critical Reviews of Microbiology. 1986, 12: 293-319.View ArticleGoogle Scholar
- Yunis R, Jarosinski KW, Schat KA: Association between rate of viral genome replication and virulence of Marek's disease herpesvirus strains. Virology. 2004, 328: 142-150. 10.1016/j.virol.2004.07.017.PubMedView ArticleGoogle Scholar
- Baigent SJ, Petherbridge LJ, Howes K, Smith LP, Currie RJW, Nair VK: Absolute quantification of Marek's disease virus genome copy number in chicken feather and lymphocyte samples using real-time PCR. Journal of Virological Methods. 2005, 123: 53-64. 10.1016/j.jviromet.2004.08.019.PubMedView ArticleGoogle Scholar
- Bumstead N, Sillibourne J, Rennie M, Ross N, Davison F: Quantification of Marek's disease virus in chicken lymphocytes using the polymerase chain reaction with fluorescence detection. Journal of Virological Methods. 1997, 65: 75-81. 10.1016/S0166-0934(96)02172-6.PubMedView ArticleGoogle Scholar
- Burgess SC, Davison TF: A quantitative duplex PCR technique for measuring amounts of cell-associated Marek's disease virus: difference s in two populations of lymphomas cells. Journal of Virological Methods. 1999, 82: 27-37. 10.1016/S0166-0934(99)00075-0.PubMedView ArticleGoogle Scholar
- Islam A, Harrison B, Cheetham BF, Mahony TJ, Young PL, Walkden-Brown SW: Differential amplification and quantification of Marek's disease viruses using real-time polymerase chain reaction. Journal of Virological Methods. 2004, 119: 103-113. 10.1016/j.jviromet.2004.03.006.PubMedView ArticleGoogle Scholar
- Islam A, Cheetham BF, Mahony TJ, Young PL, Walkden-Brown SW: Absolute quantification of Marek's disease virus and herpesvirus of turkeys in chicken lymohocyte, feather tip and dust samples using real-time PCR. Journal of Virological Methods. 2006, 132: 127-134. 10.1016/j.jviromet.2005.10.009.PubMedView ArticleGoogle Scholar
- Reddy SM, Witter RL, Gimeno I: Development of quantitative-competitive polymerase chain reaction assay for serotype 1 Marek's disease virus. Avian Diseases. 2000, 44: 770-775. 10.2307/1593048.PubMedView ArticleGoogle Scholar
- Islam AFMF, Walkden-Brown SW, Groves PJ, Underwood GJ: Effects of vaccine dose, virus challenge dose and interval from vaccination to challenge on protection of broiler chickens against Marek's disease virus challenge. Australian Veterinary Journal. 2007, 85 (9): 348-355. 10.1111/j.1751-0813.2007.00195.x.PubMedView ArticleGoogle Scholar
- Islam AFMF, Walkden-Brown SW, Groves PJ, Underwood GJ: Kinetics of Marek's disease virus (MDV) infection in broiler chickens 1: effect of varying vaccination to challenge interval on vaccinal protection and load of MDV and herpesvirus of turkey in the spleen and feather dander over time. Avian Pathology. 2008, 37 (3): 225-235. 10.1080/03079450701802230.View ArticleGoogle Scholar
- Witter R, Solomon JJ, Nazerian Champion: Long-term studies of Marek's disease infection in individual chickens. Avian Diseases. 1971, 15 (2): 346-65. 10.2307/1588707.PubMedView ArticleGoogle Scholar
- Renz KG: In vitro and in vivo characterisation of selected Australian isolates of Marek's disease virus. 2008, Ph.D Thesis, University of New England, Australia, [https://e-publications.une.edu.au/vital/access/services/Download/une:1909/SOURCE04?view=true]Google Scholar
- Baigent SJ, Davison F: Marek's disease virus: Biology and life cycle. Marek's Disease: an Evolving Problem. Elsevier Academic Press 2004.Google Scholar
- Baigent SJ, Smith LP, Currie RJW, Nair VK: Replication kinetics of Marek's disease vaccine virus in feathers and lymphoid tissus using PCR and virus isolation. Journal of General Virology. 2005, 86: 2989-2998. 10.1099/vir.0.81299-0.PubMedView ArticleGoogle Scholar
- Fenner F, Fantini B: Biological Control of Vertebrate Pests: The History of Myxomatosis - an Experiment in Evolution. CABI Publishing 1999.Google Scholar
- Walkden-Brown S, Cooke J, Islam A, Renz K, Hussain Z, Islam A, Tannock G, Groves P: Pathotyping of Australian isolates of MDV. Proceedings of the 8th International Marekís Disease Symposium: 6-10th July, Townsville, Australia. 2008Google Scholar
- Islam AFMF, Walkden-Brown SW, Wong CW, Groves PJ, Burgess SC, Arzey KE, Young PL: Influence of vaccine deposition site on post-vaccinal viraemia and vaccine efficacy in broiler chickens following in ovo vaccination against Marek's disease. Avian Pathology. 2001, 30: 525-533. 10.1080/03079450120078725.PubMedView ArticleGoogle Scholar
- Lunn D, Thomas A, Best N, Spiegelhalter D: WinBUGS Bayesian modelling framework: concepts, structure, and extensibility. Statistics and Computing. 2000, 10: 325-337. 10.1023/A:1008929526011.View ArticleGoogle Scholar
- Raftery A, Lewis S: How many iterations in the Gibbs sampler?. Bayesian Statistics. Oxford University Press, 4: 1992.Google Scholar
- Carlin B, Lewis T: Bayesian Methods for Data Analysis. Chapman & Hall/CRC 2008.Google Scholar
Copyright
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.