 Research
 Open Access
 Published:
Breeding with resistant rams leads to rapid control of classical scrapie in affected sheep flocks
Veterinary Research volume 42, Article number: 5 (2011)
Abstract
Susceptibility to scrapie, a transmissible spongiform encephalopathy in sheep, is modulated by the genetic makeup of the sheep. Scrapie control policies, based on selecting animals of resistant genotype for breeding, have recently been adopted by the Netherlands and other European countries. Here we assess the effectiveness of a breeding programme based on selecting rams of resistant genotype to obtain outbreak control in classical scrapieaffected sheep flocks under field conditions. In six commerciallyrun flocks following this breeding strategy, we used genotyping to monitor the genotype distribution, and tonsil biopsies and postmortem analyses to monitor the occurrence of scrapie infection. The farmers were not informed about the monitoring results until the end of the study period of six years. We used a mathematical model of scrapie transmission to analyze the monitoring data and found that where the breeding scheme was consistently applied, outbreak control was obtained after at most four years. Our results also show that classical scrapie control can be obtained before the frequency of nonresistant animals is reduced to zero in the flock. This suggests that control at the national scale can be obtained without a loss of genetic polymorphisms from any of the sheep breeds.
Introduction
Classical scrapie in sheep is the eldest known transmissible spongiform encephalopathy (TSE) and is present in all sheepproducing countries except Australia and New Zealand [1, 2]. Infection mostly occurs at very young age and clinical signs of this fatal disease are visible after a variable incubation period of one or more years dependent of genotype [3]. Frequently observed clinical signs are uncoordinated movement (ataxia), abnormal posture and severe scratching and rubbing. During the incubation period the prion protein PrP^{Sc} slowly accumulates in the animal and can be detected for most sheep genotypes in lymphoid organs such as tonsils before clinical signs become visible [4, 5].
Scrapie control became a priority in many countries some ten years ago. This was motivated in part by the theoretical possibility that bovine spongiform encephalopathy (BSE) may in the past have been introduced into sheep through consumption of feed supplements, with potential consequences to public health [6, 7]. This possibility became apparent after experimental infection of sheep with BSE showed that sheep can be infected via the oral route and that the resulting clinical symptoms are very similar to scrapie [8, 9].
The susceptibility to scrapie is modulated by polymorphisms of the sheep prion protein (PrP) gene. For classical scrapie strains the most relevant polymorphisms occur at codons 136, 154 and 171 of the PrP gene [10–13], and for the recently discovered atypical scrapie an additional polymorphism has been identified at codon 141 [14–16]. For atypical strains the betweenanimal transmissibility in the field, if at all present, is likely to be very low [2].
In this study we focus on classical scrapie and thus on the three aforementioned codons. Four alleles (VRQ, ARQ, AHQ and ARR) were observed in this study, and each one of the corresponding 10 possible genotypes occurred in the study flocks. The VRQ allele confers high susceptibility to most strains of classical scrapie, the ARQ allele is associated with moderate susceptibility and the AHQ allele may be associated with increased resistance and longer incubation periods [17, 18]. The allele ARR is known to confer resistance to all strains of classical scrapie, with the homozygote genotype ARR/ARR being extremely resistant, and the heterozygote ARQ/ARR and AHQ/ARR genotypes being only rarely affected by classical scrapie [19, 20].
The current scrapie control regulation in the European Union consists of the following two minimum requirements: culling of all animals of susceptible genotype in infected flocks (if no genetic testing is done, all animals need to be culled), and the genetic testing of rams intended for breeding in scrapiefree flocks of "high genetic merit" followed by culling of the susceptible rams (EC Regulation No 999/2001 [21]). Several years before the EC scrapie control regulation came into force, some member states already had a national breeding programme that sheep breeders could join on a voluntary basis, notably the Netherlands (started in 1998), Great Britain (started in 2001) and France (started in 2001). In the Netherlands, the selection of resistant rams (ARR/ARR) for breeding was made compulsory for the sheep industry in November 2004.
The purpose of this study is to assess under field conditions how a breeding programme based on selecting rams of resistant genotype affects the transmission of classical scrapie. To this end we recruited six scrapieaffected commerciallyrun sheep flocks into our study that were willing to join the (initially voluntary) breeding programme. Both the genotype distribution and scrapie transmission in each of these flocks was monitored over a period of 4 to 6 years.
Materials and methods
Study population
Six commerciallyrun sheep farms (labelled from A to F) were selected in 2000 on the basis of the following criteria. We were interested in flocks with a recent history of (confirmed) scrapie cases, with individual identification and registration of the animals, with the perspective of continuation of the flock for at least five years, and with a limited sale of animals over six months of age. Furthermore, the farmer needed to be willing to take part in the research project for at least four years, and to use only ARR/ARR rams for breeding starting from, at the latest, the 2001 mating season until the end of the study period, and to purchase no animals except ARR/ARR rams during this period. The research project bought all live animals older than 12 months that were culled by the farmer. The farms kept breeding flocks with an average size of between 90 and 100 ewes. Flock A has a size of about 125 ewes (purebred Texel and Suffolks). In this flock scrapie had been confirmed in 1998, 1999 and 2000 in four sheep. Flock B and C both have a flock size of about 125 ewes (unregistered Texel). In flock B scrapie had been confirmed in 1997 and 1999 in four sheep, and flock C had about 1015 unconfirmed scrapiesuspected sheep each year since 1992. Flock D has an average size of 50 ewes (crossbreeds), in which scrapie was confirmed in one sheep in 2000. Flock E has about 150 ewes (crossbreeds). Since 1997 there had been a few scrapie cases each year in this flock, with a peak of seven cases in the 19992000 winter, and scrapie was officially confirmed in 2000 in three sheep. Flock F has 35 Swifter ewes, and scrapie had been confirmed in 1998 and 2000 in four sheep. In total the study population consisted of about 550600 sheep each year.
Monitoring approach
Each flock in the study was monitored by carrying out the following three types of analysis:

1.
Genotyping of all breeding animals older than 12 months at the start of the study and of all new breeding animals (own offspring and purchased rams) in subsequent years. PrP genotypes were determined (at least at codons 136, 154, and 171) by a routine TaqMan test [18] that is completely automated. A second test, based on sequencing, was used as a confirmatory test on random selected samples.

2.
Yearly tonsil biopsy of all breeding sheep older than 12 months, and testing these biopsies for presence of PrP^{Sc}, the pathogenic isoform of the prion protein [22, 23]. The biopsy technique has been described in detail by Schreuder et al. [4].

3.
Postmortem testing for the presence of PrP^{Sc} in tonsils and brains of all animals older than 12 months that were culled or died of intercurrent disease on the farm. At necrosy, the brain stem was collected together with both palatine tonsils for histopathological examination [22, 23].
The yearly sampling of blood and tonsils took place between May and September, after weaning of the lambs and before the mating season started. Of the genotyping and the test results, sheep owners only received the genotype of the rams, to be able to select resistant rams for breeding. No information about ewe genotypes or scrapie status of the animals was given during the entire project period.
Two of the six flocks (A, F) were monitored during 4 consecutive years and the remaining 4 flocks during 6 years. The monitoring length of six years was designed to be sufficient to monitor the vast majority of animals born before 2001 throughout their lives. For later cohorts the rightcensoring of the data was expected to become increasingly important to take into consideration when interpreting these data.
Assessment of scrapie cases
The presence of PrP^{Sc} was tested in the yearly tonsil biopsies and in the obex and tonsils collected at necropsy by histopathological examination [22, 23]. Brain obex: a scrapie positive diagnosis was made when the typical pattern of PrP^{Sc} (intracellular and cell membrane staining of neurons and glial cells) was found in the obex. In those cases in which the dorsal motor nucleus of the Vagus (DMNV) was not available for examination (e.g. severe autolysis or inadequate sampling) no diagnosis was made and the sample was considered to be unsuitable. A negative diagnosis was made when the DMNV was available for examination but no PrP^{Sc} could be detected. Tonsil: a diagnosis of scrapie was made when the typical pattern of PrP^{Sc} (on the cell membrane of follicular dendritic cells and in follicular macrophages) was detected in the lymphoid follicles. When a tonsil biopsy contained <3 lymphoid follicles, it was considered unsuited for scrapie diagnosis. A negative diagnosis was made when the number of lymphoid follicles was ≥3 and no PrP^{Sc} could be detected.
Because testing for the presence of PrP^{Sc} is based on the monitoring procedure and not on the presence of clinical signs of scrapie, histopathological examination can also reveal subclinal scrapie. In this study a scrapie case is defined as a sheep with detected PrP^{Sc} in brain and/or tonsils independent of the presence of clinical signs.
All animals of the flock present at the yearly sampling were tested by tonsil biopsy (and occasionally by almost coincidental postmortem examination). Scrapie prevalence was estimated for each year by the number of scrapiepositive tested animals divided by the number of tested animals. The prevalence can be underestimated as PrP^{Sc} will not be detected in tonsils for all susceptible genotypes (e.g. VRQ/ARR). If a sheep was culled within some weeks after the yearly sampling, the postmortem result was used to compute the prevalence. The incidence rate was estimated by the number of new detected scrapie cases during a year following the last yearly sampling, divided by the number of animaltime units at risk. New scrapie cases were detected by tonsil biopsy at the next yearly sampling or at post mortem examination during the current year. As VRQ/ARR scrapie cases will be tonsil negative and only detected postmortem, these cases may be incorporated relatively late in the computation of the incidence rate (in comparison of scrapie cases with tonsil positive biopsies).
Mathematical modelling
We first describe the overall approach of the modelling, before we discuss the mathematical and statistical (i.e. parameter estimation) details.
Modelling approach
The main aim of our mathematical modelling is to quantify how the scrapie transmission risk in a flock is changing in time due to the selective breeding. The appropriate measure of the transmission risk is the basic reproduction number R_{ 0 } . It is defined in our context as the expected number of secondary cases of infection produced by a single primary scrapie infection in a population in which all animals of susceptible genotype are still uninfected [24–26]. R_{ 0 } = 1 is the threshold value of the reproduction number below which the withinflock infection can not sustain itself and only small outbreaks can be expected in the flock [27, 28]. The breeding programme can thus be considered successful if the value of R_{ 0 } is brought below 1. After developing a model for scrapie transmission in a flock harbouring different genotypes, we can calculate the basic reproduction number as a function of the genotype frequencies in the flock. This enables us to evaluate the efficacy of the breeding programme in reducing the transmission risk by calculating, for each flock B to E individually, the time evolution of the basic reproduction number R_{ 0 } . For flocks A and F it was not possible to estimate R_{ 0 } , either because the number of scrapie cases was too low (flock A) or because the replacement strategy was so different from year to year that the sheep population did not have a stable age distribution (flock F). A stable host age distribution is a prerequisite for being able to define and calculate R_{ 0 } for scrapie in sheep.
In Figure 1 we schematically depict the different aspects to be taken into account in our scrapie transmission model. In our model, R_{ 0 } at a given time t is expressed in terms of a number of estimable parameters. To calculate the effect of ram selection on R_{ 0 } , the following processes need to be taken into account in the model: population genetics of the flock, agedependent culling and replacement of ewes, scrapie transmission and incubation. Therefore the estimable parameters include the agegenotype distribution at time t, the agedependent survival functions of healthy and of infected sheep, and genotypespecific, agedependent susceptibility parameters. These model parameters were all estimated from the observed data.
Modelling details
The calculation of R_{ 0 } started by formulating a transmission model, designed as a tool to analyze how the transmission risk in a flock develops under selective breeding. This model describes the probability of a change in the infection status of an animal of age a and genotype γ, together with the probability of removal of the animal (i.e. fallen stock and animals sent for slaughter), in a way defined in Table 1 and Table 2. The age a runs from 0 (corresponding to the first year of life) until a flockdependent maximum age. The model assumes lateral transmission between animals [29]. The role of the environment as a faecaloral reservoir of infectivity is implicit in this assumption, as is the possible role of infectious placentas and transmission via milk [30]. These specific possible routes of transmission were left implicit in order to avoid introducing parameters that could not be estimated from the available data (e.g. it is as yet unclear how large the contribution to transmission of an infected placenta is, relative to other possible transmission routes). In a model explicitly including the placenta route the infectiousness of the placenta would need to be dependent on the genotype of the infected animal's fetus, in view the experimental results of Andreoletti et al. [31]. In our model the infectiousness of an infected animal is fully determined by its age and genotype; for the predicted effect of a breeding programme this is a conservative approximation. Another conservative model approximation made was that the infectiousness of infected animals is independent of the genotype. The relationship between infectiousness and incubation stage (or time since infection) was assumed to be linear, in close agreement with previous estimates of this dependence [6].
In Table 2 we defined the different variables and parameters of the model. The basic reproduction number R_{ 0 } (t) can be expressed in the following way in terms of R_{ 0 } (0) (i.e. its value at t = 0, just before the start of the breeding programme) and model parameters (as listed in Table 2):
In the above expression, S is the survival function for uninfected sheep, which was assumed to be independent of time (demographic equilibrium), and ${\tilde{S}}_{\gamma}$ is the survival function for infected animals of genotype γ. The function f_{ γ } (t,a) describes the genotype distribution in time, and was calculated directly from the genotyping data. f_{ γ } (0) is the genotype distribution at t = 0. The effect of the selective breeding programme was investigated by inspecting how R_{ 0 } evolved in time as a result of the changes in f_{ γ } in time brought about by the breeding programme. We estimated the model parameters R_{ 0 } (0), $S\left(a\right),{\tilde{S}}_{\gamma}\left(a\right)$ and g_{ γ } (a) for each flock individually. The τ dependence of the infectiousness level was β(τ) taken to be linear, in close agreement with a sheeptosheep infectiousness profile estimated by Ferguson et al. [6]. We note that the estimated values for R_{ 0 } (t) are only weakly dependent on the shape of β(τ). This can be seen from Eq. (1), where β(τ) occurs in a similar way in both numerator and denominator. For a derivation of Equation 1 (for the special case of $\tilde{S}=S$, i.e. for the case that scrapie infection does not influence survival) we refer to Diekmann and Heesterbeek [32].
The calculation of R_{ 0 } (t) requires calculation of R_{ 0 } (0), i.e. the value of the basic reproduction number just before the start of the breeding programme. The approach we adopted for this calculation was tailored to an optimal use of the field data for flock C. For this flock the data were most informative due to the relatively large number of detected scrapie infections. The birth cohorts 1999 and 2000, that produced 20 cases out of the total of 31 detected in flock C, showed very similar incidence patterns: both the total scrapie incidence as well as the distribution of cases over age at first detection were almost the same for both cohorts. This motivated us to assume a constant force of infection for the time period in which the cases born in 1999 and 2000 were infected: λ(t) = λ = constant.
Clearly the force of infection should be ultimately reduced by the breeding programme if it is effective. However, due to the long incubation time of the infection, there will generally be a delay of a number of years between a reduction in R_{ 0 } and the concomitant reduction in the force of infection. Because of this effect, the assumption of a constant force of infection as experienced by the 1999 and 2000 cohorts was a consistent approximation provided that most of the infections in these cohorts were occurring in the first few years of the animal's life.
For a constant force of infection, the basic reproduction number R_{ 0 } (0) just before the start of the breeding programme can be expressed in terms of the uninfected proportion Q_{ γ } (a) of animals of age a and genotype γ as follows:
Parameter estimation
As the noninfected proportions Q_{ γ } (a) can not be estimated directly from the field data (because "negative" is not the same as "noninfected"), we relate these to a number of model parameters which subsequently will be estimated from the field data. We express the proportions Q_{ γ } (a) in terms of the conditional probability p_{ γ } (a) of being infected at age a (for details we refer to the Additional file 1). The proportions Q_{ γ } (a) can then be expressed in terms of a (small) number of parameters by choosing model parameterizations for p_{ γ } (a). The simplest model (containing only one parameter) is an ageindependent probability:
Slightly more complicated models (containing two parameters) are the following: p_{ γ } (a) decreases exponentially with age:
or p_{ γ } (a) is reduced only after the first year of life and remains the same afterwards:
We assumed that infections of ARQ/ARQ animals born in 1999 or 2000 in flock C will become first detected at most two years after the end of the year of life in which the infection took place. This assumption was motivated by the fact that 19 out of 20 cases detected in flock C were three or less years old when first detected. Based on this assumption, the delay between infection and detection could be described using two estimable parameters defined in the Additional file 1. Including also the parameters p_{γ = ARQ/ARQ} (all Models) and a (Models II and III) we thus worked with 3 (Model I) or 4 (Models II and III) estimable parameters. The model likelihood (given in the Additional file 1) was numerically optimized to obtain maximumlikelihood estimates for these parameters. For each of the models, the model fit was judged by Pearson's Chisquare statistic (and associated Pvalue from a chisquare approximation).
As for the genotype ARQ/ARR in flock C there were no cases observed, it was not possible to apply the above approach of analysis to this genotype (although the point estimate of p_{γ = ARQ/ARR} was easily calculated (p_{γ = ARQ/ARR} = 0), we would also seek to calculate an upper confidence bound of p_{γ = ARQ/ARR}). Furthermore, for all genotypes in flocks B, D and E the number of scrapie cases observed is too limited to apply the above analysis. For these flocks/genotypes we therefore followed a more pragmatic approach, in which we derived a pessimistic estimate of R_{ 0 } (0) by using an estimate of the noninfected proportions Q_{ γ } (a) based on the cohort with the highest incidence in the data. These cohorts are in all cases "early" cohorts, i.e. from before the start of the breeding programme. For this cohort we calculated the observed total fraction positive F_{ γ }^{+} as the number of detected scrapie cases divided by the total number of animals of genotype γ which were either found test positive or at least reached a given minimum age without being found positive (to avoid misclassifying animals that were infected but not yet incubating long enough to be detected). For further details we refer to the Additional file 1. We note that as the model calculations are not using the scrapie infection data for "late" cohorts (as they only use data for early cohorts to estimate the starting value R_{ 0 } (0)), they are not affected by any rightcensoring problems of detecting scrapie in the late cohorts.
Results
Genotype distribution
To monitor the genotype distribution, the results of in total 1175 sheep were used. Four alleles (VRQ, ARQ, AHQ and ARR) were observed in this study, and each one of the corresponding ten possible genotypes occurred in the study flocks. Pie charts in Figure 2 show how the observed genotype distribution of the ewe population in each of the six flocks evolved during the study period. In flocks AE we observed a genetic improvement by an increase of the ARR allele and the ARR/ARR genotype in these flocks, with even an increase up to 100% for flock D in 2006. Although in flock F ARR/ARR rams were used, unfortunately the farmer did not use the offspring of these rams to replace ewes. As a result hardly any genetic improvement was achieved. Monitoring of flock F was ceased late 2004 after the flock was sold by the owner.
More details on the genotype distribution of each flock can be found in Table S1 of the Additional file 1.
Scrapie cases
To monitor scrapie infections 3819 tonsil biopsies and postmortem samples (tonsils and obex) of 1168 sheep were analysed for the presence of PrP^{Sc}. Eleven percent of the tonsil biopsy samples were considered unsuited for diagnosis as the tonsil biopsy contained <3 lymphoid follicles. Six of the 1168 sheep (0.5%) were unsuited due to severe autolysis. For 44 postmortem sheep (3.8%) the tonsils were not available and for 3 sheep (0.3%) the obex was missing due to inadequate sampling.
In total 59 scrapie cases were detected during the entire monitoring period, of which 51 cases were sheep without the ARR allele (2 VRQ/VRQ, 16 VRQ/ARQ, 33 ARQ/ARQ) and 8 cases were VRQ/ARR (Table 3). No scrapie cases were detected in birth cohorts from 2001 onwards in any of the flocks, i.e. in animals born after the start of the control programme in 2000. The eight VRQ/ARR cases were only detected by positive obex during postmortem analysis, with all tonsil biopsies of these animals as well as the tonsils collected at post mortem being negative. For 36 of the 51 cases without the ARR allele, scrapie was detected first by tonsil biopsy when alive, and later confirmed by positive postmortem analysis. For most of the 15 remaining cases recent tonsil biopsy data were missing for different reasons. Detailed information on the 59 scrapie can be found in Table S2 of the Additional file 1.
Flocks BF were affected by a substantial scrapie infection (Table 3). Flock C is the most heavily affected flock, in which in total 31 scrapie cases were detected during the monitoring period, all of ARQ/ARQ genotype. Given the high incidence in this flock despite the absence of the VRQ allele, one might hypothesize that this outbreak involved a scrapie strain different from the strain(s) affecting the other flocks. Unfortunately the remaining flock A could not provide information on the effect of the breeding programme on scrapie transmission, as the scrapie outbreak in this flock faded out one year after the selective breeding programme started. Only a single scrapie case was detected in 2000 by tonsil biopsy and this sheep was culled in 2001. Up until (and inclusive) the last biopsy round in 2004 no further positives were found. Therefore, monitoring of flock A was ceased after that year.
Table 4 and Table 5 show the prevalence and incidence rate of scrapie in the six flocks throughout the monitoring period, irrespective of genotype. A decline of scrapie cases in time was found in the four flocks BE. This decline in both scrapie prevalence and incidence provides direct evidence for the efficacy of the breeding programme. Furthermore, none of the animals born after the start of the programme was found scrapie positive. Given the prevalence data shown in Table 4, it is unlikely (flock B) or even very unlikely (flocks CE) that a stochastic extinction of scrapie would have resulted in any of these flocks, had selective breeding not been applied. Further evidence for the efficacy of the breeding programme is given by the results for the basic reproduction number presented below. Flock F did not comply with the breeding programme (see next section); here a low scrapie prevalence remained throughout the four years of monitoring.
Parameter estimates
The main parameter of interest in our analysis is R_{ 0 } (t), and the other parameters are mainly relevant as determinants of R_{ 0 } (t) (through Eq.(1)). However, the results of the parameter estimation for the Models IIII (introduced in section 2.4.3) are also interesting by themselves as they describe possible agedependence of transmission. Table 6 displays these estimates based on the data for flock C. The best fit was obtained for Model III, as this model had the lowest value for the Chisquare statistic and the highest P value (P = 0.38). The maximumlikelihood estimates for the two parameters describing the agedependent infection probability (see Eq. (5)) are: p_{ARQ/ARQ} = 0.416 (0.273, 0.574) and α = 0.059 (0.003, 0.29). The 95% confidence intervals (in parentheses) were constructed with the score test [33]. The score test was used (and not the likelihood ratio test) to avoid any problems with the chisquare approximation of the likelihood ratio test for parameter values close or equal to the domain boundaries 0 and 1. We consider the fit of model III to be quite satisfactory, given that the data (see Additional file 1: Table S4) comprised an exceptional animal that was first found positive at the old age of six years. When this late positive was removed from the dataset, both Models II and III yielded perfect fits (P = 1). We note that the agedependency parameter α in Model III is significantly smaller than 1 (upper confidence bound equals 0.29). Therefore, Model III with α = 1 (agedependent susceptibility), is in fact rejected.
For the parameter estimation results for flock B, D and E and for p_{γ = ARQ/ARR} in flock C we refer to the Additional file 1, as their main relevance is in determining the R_{ 0 } (t) estimates presented in the next section.
Effect of the breeding programme on scrapie transmission
For each of the four flocks BE we calculated R_{ 0 } as a function of time (Figure 3). For flocks A and F it was not possible to estimate R_{ 0 } for reasons explained in section 2.4.1. The squares display R_{ 0 } calculated using upperbound estimates for the susceptibility of S/ARR genotypes and lowerbound estimates for the susceptibility of S/S genotypes, and further assuming that scrapie susceptibility is age independent in all four flocks. This combination of parameter choices yields the most conservative model prediction of the effect of the breeding programme. The corresponding R_{ 0 } results (squares in Figure 3) therefore serve as upper confidence bounds when they exceed the point estimate results (circles).
The point estimates of R_{ 0 } are obtained from point estimates of the agedependent, genotypespecific susceptibility parameters in the model. From Figure 3 we thus conclude that in flocks B, D and E R_{ 0 } is reduced to values significantly below one after four years, and in flock C this is the case already after three years. Figure 3 shows that for all three flocks BD the point estimates of R_{ 0 } (circles) are reduced below the threshold value R_{ 0 } = 1 within at most three years. For flock E the point estimate of R_{ 0 } is already below threshold at the start of the study period, however the upper confidence bound is not. It is relevant to note that the results for flocks B, D and E illustrate that scrapie control (i.e. when the upper confidence bound for R_{ 0 } declines to below 1) can be achieved before the presence of nonARR/ARR animals has reduced to zero. This can be seen from the R_{ 0 } point estimates, which are nonzero implicating the presence of nonARR/ARR animals.
Discussion
Six commerciallyrun flocks were selected to evaluate the effect of a breeding programme for scrapie control, consisting of selecting homozygote resistant rams for mating. We used genotyping to monitor the genotype distribution, and tonsil biopsies and postmortem analyses to monitor the occurrence of scrapie infection throughout the study period of six years. The results show that the breeding program leads to rapid outbreak control, typically within at most four years of selective breeding. Furthermore, our results show that to reach the situation where R_{ 0 } < 1 (i.e. outbreak control) the number of nonARR/ARR animals does not have to be reduced to zero in the flock.
For the four flocks of interest, B, C, D and E, we used mathematical modelling to estimate the basic reproduction number R_{ 0 } through time, separately for each flock. This approach of combining observational data with modelling is important to obtain more insight in the dynamics of the infection, in particular for slowprogressing diseases like scrapie. Rightcensoring of scrapie detection in the late birth cohorts of sheep may have resulted in a slight underestimation of the prevalence and incidence of the latest years of the monitoring period. However, the possible rightcensoring did not affect the calculation of R_{ 0 } , as this calculation is only using the scrapie infection data of early cohorts. At the start of the programme, the four withinflock R_{ 0 } values were estimated to be in the range of 0.88.0. This range is of the same order of magnitude as previous estimates of the withinflock R_{ 0 } for other scrapie outbreaks by Matthews et al. [34] (3.9) and Hagenaars et al. [35] (27.5).
Our results show that to reach the situation where R_{ 0 } becomes significantly smaller than 1 (as requested for outbreak control) the number of nonARR/ARR animals does not have to be reduced to zero in the flock. On a national scale this implies that it will be possible to maintain the diversity of susceptibility alleles desired in case a new scrapie (or TSE) strain targeting ARR would arise in the future. It is important to note that when a reduction of R_{ 0 } to below one is achieved, this does not mean that the breeding programme can simply be ceased. Although expected to decline with time, infection risks to susceptible genotypes, e.g. from a scrapiecontaminated environment, are expected to remain present for some time. Therefore, reintroducing the use of nonARR/ARR rams for breeding without risking to revert to a situation with R_{ 0 } > 1 requires a new breeding strategy that avoids mating such rams with nonARR/ARR ewes.
An intermediate modelling result of interest relates to flock C. The relatively high number of scrapie cases observed in this flock allowed a detailed estimation of a model parameter controlling agedependency of scrapie susceptibility. It was found that in this flock, the null hypothesis of an ageindependent susceptibility had to be rejected. The maximumlikelihood estimate of the agedependency parameter corresponds to a scenario in which scrapie infection is occurring predominantly in the first year of life. This is in line with what is often assumed to be the case, at least in flocks affected by substantial scrapie outbreaks, but has only rarely been established based on quantitative evidence [36].
Flocks B to F were affected by a substantial scrapie infection at the start of the breeding programme, and except for flock C scrapie was mainly associated with the VRQ allele. Flock C was the most heavily affected flock, in which in total 31 scrapie cases were detected, all of ARQ/ARQ genotype. Given the high incidence in this flock despite the absence of the VRQ allele in the sheep population, one might hypothesize that this outbreak involved a scrapie strain different from the strain(s) affecting the other flocks. Certainly, in this flock the estimated (relative) susceptibility of the ARQ/ARQ genotype was much higher than in the other flocks.
We expect that possible betweenstrain differences in genetic susceptibility will not limit the effectiveness of the breeding strategy found here, given that no classical scrapie strain thus far has escaped ARRassociated resistance. In line with this expectation, the breeding programme proved successful in all four flocks analyzed despite betweenflock differences in the estimated genotypespecific susceptibilities.
As a result of the control programme not only the reproduction number R_{ 0 } , but also the infection pressure (or force of infection) in the field will decrease in time. However, due to the long incubation period of a scrapie infection, a delay of a few years is expected between the reduction in R_{ 0 } and the reduction in infection pressure. This is in line with the detection of new scrapie cases during this monitoring study, which were all born before the start of the programme. Extending the breeding strategy by a removal of scrapiesusceptible ewes on the basis of their genotype would accelerate the reduction of both R_{ 0 } as well as the infection pressure. However, when considering the rapid outbreak control as observed in this study, the use of resistant rams seems sufficient and can be recommended as a control strategy in scrapieaffected countries.
References
 1.
Detwiler LA, Baylis M: The epidemiology of scrapie. Rev Sci Tech. 2003, 22: 121143.
 2.
Fediaevsky A, Tongue SC, Noremark M, Calavas D, Ru G, Hopp P: A descriptive study of the prevalence of atypical and classical scrapie in sheep in 20 European countries. BMC Vet Res. 2008, 4: 1910.1186/17466148419.
 3.
Hadlow WJ, Kennedy RC, Race RE: Natural infection of Suffolk sheep with scrapie virus. J Infect Dis. 1982, 146: 657664.
 4.
Schreuder BE, van Keulen LJ, Vromans ME, Langeveld JP, Smits MA: Preclinical test for prion diseases. Nature. 1996, 381: 56310.1038/381563a0.
 5.
Schreuder BE, van Keulen LJ, Vromans ME, Langeveld JP, Smits MA: Tonsillar biopsy and PrP^{Sc} detection in the preclinical diagnosis of scrapie. Vet Rec. 1998, 142: 564568. 10.1136/vr.142.21.564.
 6.
Ferguson NM, Ghani AC, Donnelly CA, Hagenaars TJ, Anderson RM: Estimating the human health risk from possible BSE infection of the British sheep flock. Nature. 2002, 415: 420424. 10.1038/nature709.
 7.
Kao RR, Gravenor MB, Baylis M, Bostock CJ, Chihota CM, Evans JC, Goldmann W, Smith AJA, McLean AR: The potential size and duration of an epidemic of bovine spongiform encephalopathy in British sheep. Science. 2002, 295: 332335. 10.1126/science.1067475.
 8.
Foster JD, Parnham D, Chong A, Goldmann W, Hunter N: Clinical signs, histopathology and genetics of experimental transmission of BSE and natural scrapie to sheep and goats. Vet Rec. 2001, 148: 165171. 10.1136/vr.148.6.165.
 9.
Houston EF, Gravenor MB: Clinical signs in sheep experimentally infected with scrapie and BSE. Vet Rec. 2003, 152: 333334. 10.1136/vr.152.11.333.
 10.
Goldmann W, Hunter N, Smith G, Foster J, Hope J: Prp genotype and agent effects in scrapie: change in allelic interaction with different isolates of agent in sheep, a natural host of scrapie. J Gen Virol. 1994, 75: 989995. 10.1099/00221317755989.
 11.
Hunter N, Foster JD, Goldmann W, Stear MJ, Hope J, Bostock C: Natural scrapie in a closed flock of Cheviot sheep occurs only in specific PrP genotypes. Arch Virol. 1996, 141: 809824. 10.1007/BF01718157.
 12.
Hunter N: PrP genetics in sheep and the implications for scrapie and BSE. Trends Microbiol. 1997, 5: 331334. 10.1016/S0966842X(97)010810.
 13.
Spiropoulos J, Casalone C, Caramelli M, Simmons MM: Immunohistochemistry for PrP^{Sc} in natural scrapie reveals patterns which are associated with PrP genotype. Neuropathol Appl Neurobiol. 2007, 33: 398409. 10.1111/j.13652990.2007.00800.x.
 14.
Benestad SL, Sarradin P, Thu B, Schonheit J, Tranulis MA, Bratberg B: Cases of scrapie with unusual features in Norway and designation of a new type, Nor98. Vet Rec. 2003, 153: 202208. 10.1136/vr.153.7.202.
 15.
Moum T, Olsaker I, Hopp P, Moldal T, Valheim M, Moum T, Benestad SL: Polymorphisms at codons 141 and 154 in the ovine prion protein gene are associated with scrapie Nor98 cases. J Gen Virol. 2005, 86: 231235. 10.1099/vir.0.804370.
 16.
Saunders GC, Cawthraw S, Mountjoy SJ, Hope J, Windl O: PrP genotypes of atypical scrapie cases in Great Britain. J Gen Virol. 2006, 87: 31413149. 10.1099/vir.0.817790.
 17.
Belt PB, Meuleman IH, Schreuder BE, Bosde Ruijter J, Gielkens AL, Smits MA: Identification of five allelic variants of the sheep PrP gene ansd their association with natural scrapie. J Gen Virol. 1995, 76: 509517. 10.1099/00221317763509.
 18.
Bossers A, Schreuder BEC, Muileman IH, Belt P, Smits MA: PrP genotype contributes to determining survival times of sheep with natural scrapie. J Gen Virol. 1996, 77: 26692673. 10.1099/0022131777102669.
 19.
Baylis M, Chihota C, Stevenson E, Goldmann W, Smith A, Sivam K, Tongue S, Gravenor MB: Risk of scrapie in British sheep of different prion protein genotype. J Gen Virol. 2004, 85: 27352740. 10.1099/vir.0.798760.
 20.
Hunter N, Goldmann W, Smith G, Hope J: The association of a codon136 Prp gene variant with the occurrence of natural scrapie. Arch Virol. 1994, 137: 171177. 10.1007/BF01311184.
 21.
Anonymous, European Commission: Regulation (EC) No 999/2001 of the European Parliament and of the Council of 22 May 2001 laying down rules for the prevention, control and eradication of certain transmissible spongiform encephalopathies. Off J Eur Communities. 2001, 147: 140.
 22.
van Keulen LJ, Schreuder BE, Meloen RH, MooijHarkes G, Vromans ME, Langeveld JP: Immunohistochemical detection of prion protein in lymphoid tissues of sheep with natural scrapie. J Clin Microbiol. 1996, 34: 12281231.
 23.
van Keulen LJM, Schreuder BEC, Meloen RH, Poelenvan den Berg M, MooijHarkes G, Vromans MEW, Langeveld JPM: Immunohistochemical detection and localization of prion protein in braintissue of sheep with natural scrapie. Vet Pathol. 1995, 32: 299308. 10.1177/030098589503200312.
 24.
Gubbins S, Touzeau S, Hagenaars TJ: The role of mathematical modelling in understanding the epidemiology and control of sheep transmissible spongiform encephalopathies: a review. Vet Res. 2010, 41: 4210.1051/vetres/2010014.
 25.
Hagenaars TJ, Donnelly CA, Ferguson NM, Anderson RM: The transmission dynamics of the aetiological agent of scrapie in a sheep flock. Math Biosci. 2000, 168: 117135. 10.1016/S00255564(00)000481.
 26.
Hagenaars TJ, Ferguson NM, Donnelly CA, Anderson NM: Persistence patterns of scrapie in a sheep flock. Epidemiol Infect. 2001, 127: 157167. 10.1017/S0950268801005738.
 27.
Anderson RM, May RM: Infectious diseases of humans: dynamics and control. 1991, Oxford and New York, Oxford University Press, 757:
 28.
Diekmann O, Heesterbeek JAP, Metz JAJ: On the definition and the computation of the basic reproduction ratio R_{ 0 } in models for infectiousdiseases in heterogeneous populations. J Math Biol. 1990, 28: 365382. 10.1007/BF00178324.
 29.
Ryder S, Dexter G, Bellworthy S, Tongue S: Demonstration of lateral transmission of scrapie between sheep kept under natural conditions using lymphoid tissue biopsy. Res Vet Sci. 2004, 76: 211217. 10.1016/j.rvsc.2003.11.007.
 30.
Konold T, Moore SJ, Bellworthy SJ, Simmons HA: Evidence of scrapie transmission via milk. BMC Vet Res. 2008, 4: 1410.1186/17466148414.
 31.
Andreoletti O, Lacroux C, Chabert A, Monnereau L, Tabouret G, Lantier F, Berthon P, Eychenne F, LafondBenestad S, Elsen JM, Schelcher F: PrP(Sc) accumulation in placentas of ewes exposed to natural scrapie: influence of foetal PrP genotype and effect on ewetolamb transmission. J Gen Virol. 2002, 83: 26072616.
 32.
Diekmann O, Heesterbeek JAP: Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. Wiley Series in Mathematical and Computational Biology, Chichester, Wiley. 2000
 33.
Cox DR, Hinkley DV: Theoretical Statistics. 1974, London, Chapman and Hall
 34.
Matthews L, Woolhouse MEJ, Hunter N: The basic reproduction number for scrapie. Proc R Soc London B Biol Sci. 1999, 266: 10851090. 10.1098/rspb.1999.0747.
 35.
Hagenaars TJ, Donnelly CA, Ferguson NM, Anderson RM: Dynamics of a scrapie outbreak in a flock of Romanov sheep  estimation of transmission parameters. Epidemiol Infect. 2003, 131: 10151022. 10.1017/S0950268803001055.
 36.
Matthews L, Coen PG, Foster JD, Hunter N, Woolhouse MEJ: Population dynamics of a scrapie outbreak. Arch Virol. 2001, 146: 11731186. 10.1007/s007050170113.
Acknowledgements
This research was funded by the Dutch Ministry of Agriculture, Nature and Food Quality, the Dutch Commodity Boards of Meat, Milk and Eggs, the Dutch Organisation for Agriculture and Horticulture (LTONederland) and EU project ScrapieFreeSheep (QLK5CT200101733). We thank MCM de Jong and MA Smits for discussions; W Buist for statistical analysis; A Bossers, J Briaire, and F Verburg for genotyping; G Werkhoven, P Roskam, P Riedstra, D Schoenmakers, A Spithoven, G Visser, and J van Zandbergen for technical assistance, K de Haan for support and technical assistance; J Vroman for immunohistochemistry; W Wouda and his colleagues for pathology. We acknowledge the participation of the sheep farmers to the observational field study.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
GN, HvR, and PV participated in the design of the study; GN, TH and LvK performed research; GN, TH and BE analyzed data; GN, HvR, and TH drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Supplementary data
Additional file 1: . Tables with details of genotype distribution of the six flocks and detailed information of all scrapie cases are shown. Details of quantitative analyses are provided: the calculation of the noninfected proportions Q_{γ} (a), the estimation of parameters for γ = ARQ/ARQ in flock C, the estimation of parameters for flock B, D and E and of upper confidence bound for p_{ARQ/ARR} in flock C, the calculation of upper confidence bound for R_{0}, and the distribution of survival times. (DOC 342 KB)
Authors’ original submitted files for images
Below are the links to the 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.
About this article
Cite this article
Nodelijk, G., van Roermund, H.J., van Keulen, L.J. et al. Breeding with resistant rams leads to rapid control of classical scrapie in affected sheep flocks. Vet Res 42, 5 (2011). https://doi.org/10.1186/12979716425
Received:
Accepted:
Published:
Keywords
 Breeding Programme
 Bovine Spongiform Encephalopathy
 Scrapie
 Basic Reproduction Number
 Sheep Flock