Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Investigation of a Simple Model for Within-Flock Transmission of Scrapie

Abstract

Genetic control programs for scrapie in sheep build on solid knowledge of how susceptibility to scrapie is modulated by the prion protein genotype at the level of an individual sheep. In order to satisfactorily analyze the effectivity of control programs at the population level, insight is needed at the flock level, i.e., how the grouping of sheep in flocks affects the population-level transmission risk. In particular, one would like to understand how this risk is affected by between-flock differences in genotype frequency distribution. A first step is to model the scrapie transmission risk within a flock as a function of the flock genotype profile. Here we do so by estimating parameters for a model of within-flock transmission using genotyping data on Dutch flocks affected by scrapie. We show that the data are consistent with a relatively simple transmission model assuming horizontal transmission and homogeneous mixing between animals. The model expresses the basic reproduction number for within-flock scrapie as a weighted average of genotype-specific susceptibilities, multiplied by a single overall transmission parameter. The value of the overall transmission parameter may vary between flocks to account for random between-flock variation in non-genetic determinants such as management practice. Here we provide an estimate of its mean value and variation for Dutch flocks.

Citation: Hagenaars TJ, Windig JJ (2015) Investigation of a Simple Model for Within-Flock Transmission of Scrapie. PLoS ONE 10(10): e0139436. https://doi.org/10.1371/journal.pone.0139436

Editor: Hiroshi Nishiura, The University of Tokyo, JAPAN

Received: September 14, 2014; Accepted: September 14, 2015; Published: October 1, 2015

Copyright: © 2015 Hagenaars, Windig. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work was supported by the Dutch Ministry of Economic Affairs, research project codes WOT-01-002-01.06 and WOT-01-002-001.01. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Classical scrapie in sheep is a disease that potentially can be controlled by selective breeding, due to the high to full scrapie resistance of certain genotypes. Scrapie is a transmissible spongiform encephalopathy (TSE) with an incubation period of one or more years before the occurrence of clinical signs, such as uncoordinated movement, abnormal postures and severe scratching. The susceptibility to scrapie is modulated by polymorphisms of the sheep prion protein (PrP) gene [1–5]. The most important polymorphisms occur at the codons 136, 154 and 171. Five alleles (VRQ, ARQ, AHQ, ARH and ARR) are observed in The Netherlands. The VRQ allele is known to confer high susceptibility to classical scrapie, the ARQ and ARH alleles are associated with moderate susceptibility and the AHQ allele with low susceptibility. The ARR allele confers resistance, with the homozygous genotype ARR/ARR being extremely resistant. These properties make the use of exclusively ARR/ARR rams for breeding a means to breed selectively for scrapie resistance. European Union (EU) regulation since 2001 requires the selection of rams intended for breeding in scrapie-free flocks of “high genetic merit” (followed by culling of the rams with a VRQ allele). A further important scrapie control activity is the surveillance programme of testing healthy slaughtered sheep and fallen stock for scrapie by a rapid test on brainstem samples. This programme concerns animals over 18 months of age and was introduced in the EU in 2002 [6]. Since 2003, the EU requires control measures in flocks of origin of classical-scrapie positive animals in the active or passive surveillance. These measures consist of either a whole-flock cull or genotyping all animals and culling the animals of susceptible genotype and examining the brain stem of all or a sample of the culled animals of at least 12 months of age for scrapie positivity, using rapid tests.

Some member states have introduced a wider national breeding programme than requested by the EU, including The Netherlands (started in 1998), Great Britain (started in 2001) [7–11], and France (started in 2002). As these programs have to be run over many years there is a need for reliable model projections of their expected future effects [12]. In such predictive model analyses it is often desirable to quantify within-flock scrapie transmission as a function of flock genotype profile. No such quantification has been performed to date, although a statistical analysis of 30 affected flocks by McKintyre et al. [13] provided evidence that the flock genotype profile is correlated with outbreak characteristics such as mean yearly incidence, and Tongue at al. [14] identified PrP genotype and allele frequencies as flock-level risk factors for scrapie in a case-control study comprising 293 flocks. Here we use data that has accumulated due to EU statutory measures, that include the genotyping and (partial) culling and testing of flocks of origin of scrapie-positive animals. Using the Dutch culled-flocks data, we show here how this type of data can be used to calculate the basic reproduction number for within-flock transmission. We find that the data are consistent with a relatively simple transmission model assuming horizontal transmission and homogeneous mixing between animals. The model expresses the basic reproduction number for within-flock transmission of scrapie as a weighted average of genotype-specific susceptibilities, multiplied by a single parameter that is drawn from a distribution to account for random between-flock variation in non-genetic determinants such as management practice.

Materials and Methods

Culled-flocks data

The Dutch culled-flocks data (2003–2008) consist of scrapie genotyping results and scrapie infection test results in animals that were culled, as part of the mandatory scrapie control efforts, on 69 flocks of origin of scrapie index cases. The data is included as S1 Dataset. Immunohistochemistry (IHC) was used for confirmation of the positive cases detected using the rapid test. IHC and Western blotting were used to discriminate between classical and atypical scrapie. PrP genotypes were determined (at codons 136, 154, and 171) by a routine TaqMan test that is completely automated. It can detect polymorphisms 136 A to V, 154 R to H, and 171 Q to R. From 2006 onwards our TaqMan genotyping additionally distinguishes between Q and H at codon 171. When analyzing the culled-flocks data we consider total numbers of animals for each genotype across the period 2003–2008 and we therefore group the 2006–2008 ARQ and ARH results together, using the notation ARQ*. The TaqMan principle is a test in which a small part of the PrP gene is amplified. During amplification dedicated fluorescent probes are used to detect absence/presence of specific polymorphisms. A second test, based on pyro-sequencing, was used as a confirmatory test on randomly selected samples. The rapid tests used were the Prionics Check Western (2002–2006) and the Prionics Check Western SR from June 2006 onwards. Data statistics such as the overall genotype and allele frequencies across the 69 culled flocks as well as the mean detected scrapie prevalence by genotype are given in Ref. [15].

Transmission model

Our model relating within-flock basic reproduction number to the genotype distribution is a simplified version of the more general model structure described by Hagenaars et al. [16]. The most important simplification is to refrain from a stratification by age, a necessary simplification because we do not have age information on the animals tested in the culled-flocks data. The model takes the form of a genotype specific SI (susceptible-infected) model with Sγ and Iγ being respectively the proportion of animals in the flock that are susceptible and infected and have genotype γ. The change of Sγ and Iγ are modelled as follows:

Here the change in Sγ is due to the new(born) animals coming into the flock minus suspected animals becoming infected and minus suspected animals being replaced. The change in Iγ is due to suspected animals becoming infected and infected animals being replaced. The term μfγ describes the recruitment of new(born) animals of genotype γ into the flock, with μ being the replacement rate of animals and fγ = Sγ + Iγ being the frequency (proportion) of animals in the flock that has genotype γ, a frequency which is assumed to be (quasi-)stationary. This description applies to situations where no animals are bought in (closed flock), or where bought-in animals are recruited from a population with the same genotype frequency distribution. The parameter β is a transmission rate parameter; the term proportional to β represents the rate of infection transmission to animals of genotype γ, gγ the relative susceptibility of genotype γ, and hγ an infectiousness parameter. The relative susceptibility gγ is defined as the susceptibility relative to that of the reference genotype γR = ARQ*/VRQ (i.e. setting ). The infectiousness parameter hγ is introduced in order to account for between-genotype differences in how infectious an infected animal is to its flock mates. The terms μSγ and μIγ are the rates of replacement of animals. As susceptible and infected individuals are subject to the same replacement rate μ, the model neglects any scrapie-related mortality (or preferential replacement). The analysis of Matthews et al. [17] shows that vertical scrapie transmission is estimated to make only a minor contribution to the total transmission. We therefore neglect the vertical route here, as its incorporation would make the modelling considerably more complex. In this model the within-flock basic reproduction number for scrapie transmission, denoted here as , where the superscript ‘w’ is referring to “within-flock”, is expressed in terms of the genotype distribution fγ as follows:

Here the basic reproduction number is defined as the expected number of new infections in the flock caused by a single typical primary scrapie infection in the limit of negligible infection prevalence; this corresponds to the standard textbook definition [18]. The parameter ρ0 is defined as ; it serves as a transmission scale parameter and can be interpreted as a base-line value of the reproduction number corresponding to a hypothetical situation in which the flock comprises a single genotype with relative susceptibility 1. The above expression for can be obtained by employing a next-generation operator approach, and by using the observation that the operator has one-dimensional range (or equivalently: mixing is separable) as explained in Ref. [18], section 7.4.1. We set the infectiousness parameter hγ equal to one for all genotypes without ARR allele, thus assuming that these genotypes, if infected, have the same infectiousness. For those with at least one ARR allele it is set to zero, i.e. in particular we assume that the contribution of ARR/VRQ animals to R0w is negligible. This assumption is motivated by the observation that the pathogenesis in this genotype does not (or only minimally) affect the lymphoreticular system [19]. (ARR/VRQ sheep are present in 59 out of 69 culled flocks. Overall the ARR/VRQ frequency in the culled flocks was 7.1 percent (Table 6 of Ref. [15])).

In order to incorporate variation in due to causes different from the genetic content of the flock, the parameter ρ0 is taken to be a distributed quantity (i.e. randomly varying between flocks). Such causes may include farm type and management practices in particular during the lambing period (e.g. using lambing pens or not). Values for the relative susceptibilities gγ are obtained from the estimates of the detected infection prevalence in the different genotypes across all culled flocks by Hagenaars et al. [15] and listed in Table 1. We calculate the relative susceptibility from this genotype-specific prevalence, denoted by , using the following relationship that is derived in S1 Text: (1)

Here denotes the prevalence in the reference genotype. We use the genotype ARQ*/VRQ, being the most frequent genotype amongst Dutch scrapie cases, as this reference. Only the genotypes without ARR allele are relevant, as the others do not contribute to due to hγ being zero. As noted in Ref. [15], infection prevalence is approximately proportional to susceptibility when prevalence is low. The relationship (1) takes into account the non-linearity of the relationship for intermediate and high prevalence, based on approximating the dataset as one single flock. In our model we for simplicity replace the estimate gAHQ/ARQ* = 0.013 by gAHQ/ARQ* = 0. Due to the low frequency of the AHQ/ARQ* genotype and its low estimated relative susceptibility this is a good approximation.

thumbnail
Download:
Table 1. Genotype-specific scrapie risk and corresponding estimates for the relative susceptibility.

Estimates, obtained in Ref. [15] from culled-flocks data, of the genotype-specific risk of being tested scrapie positive, for genotypes (or groups of genotypes) without ARR allele, relative to the ARQ*/VRQ group of genotypes.

https://doi.org/10.1371/journal.pone.0139436.t001

Due to the absence of age structure, our model is simpler than previously published within-flock scrapie transmission models. Those publications were typically dealing with individual outbreaks for which more detailed data was available, and were reviewed in Ref. [12]. The mathematical structure of the models of most previous work (including Refs. [16,17,20–22]) reduces to that of our model when the age structure is left out, clinical onset is left implicit and disease-induced mortality is neglected.

Testing the model

The model described above assumes that, whereas the value of the transmission scale parameter ρ0 varies between flocks, the relative susceptibility parameter gγ is the same for all flocks. We seek to validate this assumption against data of culled flocks., The assumption implies that the genotype-dependent infection risk pγ,i in flock i can be expressed in terms of a flock-dependent base-line (ARQ*/VRQ) risk pARQ*/VRQ,i and a flock-independent parameter gγ as follows: (2)

Eq (2) is derived in S1 Text. The actual within-flock incidence by genotype would arise from a binomial distribution with probability pγ,j.

To test whether the model is adequate we investigated for each culled flock i with secondary cases and with at least one animal of ARQ*/VRQ genotype tested (46 flocks in total), whether it is possible to choose a value for pARQ*/VRQ,i such that for all genotypes without ARR allele the observed data are within the 95% probability range of the binomial model using gγ estimated at population level (and listed in Table 1).

Parameter estimation

For the estimation of parameters we assume that the transmission dynamics in the flocks was in an endemic phase, i.e. in (quasi-)stationary equilibrium. This assumption is necessary due to the absence of longitudinal information on prevalence. It represents an approximation, as within-flock dynamics is expected to be out of (quasi-)equilibrium for at least the early and late parts of the outbreak period. We believe that for culled flocks with at least one detected secondary scrapie case, this approximation is justifiable due to the dominance of the quasi-stationary part of outbreaks indicated by the modelling results in figure 2 of Ref. [21]. In endemic equilibrium the model equations provide a relationship between the parameter ρ0 and the proportion infected i* of animals without ARR allele as follows: (3)

Here Γ is the set of genotypes without ARR allele; this relationship is derived in S1 Text. For a subset of culled flocks, we use the above equation to estimate ρ0 for each flock separately by assuming the flock is in endemic equilibrium. In order to use these results for modelling the between-flock variation in the parameter ρ0 due to flock-specific aspects different from the genetic content of the flock, we fit a Weibull distribution to the histogram of estimates. We relate the proportion i* to the proportion i+ found positive of tested animals without ARR allele by assuming that , with Se the (unknown) test sensitivity in detecting scrapie infection in animals without ARR allele of at least 12 months of age in endemically affected flocks. Below we motivate our approximation of assuming Se to be independent of genotype. We consider different values of Se across the range [0.55–0.95] to analyze the influence of this parameter on our results, choosing Se = 0.75 as a default value. This range is motivated by the sensitivity of the test of close to 95% as evaluated on scrapie cases confirmed by Western Blot of the brainstem [23] and the notion that early on in the incubation period scrapie infection has not yet propagated to the brainstem [24]; detected scrapie prevalence in the culled flocks suggests that a sensitivity below 0.55 is unlikely. The subset of flocks is obtained by requiring that at least two positive cases were found and in addition that at least eight animals of the genotype ARQ*/VRQ (the most common genotype amongst Dutch cases [15]) in this flock were tested. The first requirement was made because with only one detected case the assumption of endemicity was deemed too crude, and the second requirement was made in order to avoid that the estimated endemic prevalence in the flock (Eq (3)) became dominated by noise.

We use a genotype-independent parameter Se for the sensitivity of the rapid test. In general the sensitivity of the scrapie test may be expected to depend on genotype, as the sensitivity depends on how far the animal has progressed towards clinical onset [19,24–27], and the incubation period is dependent on the genotype. However, as can be seen from the age-at-onset results in the Electronic Supplementary Material of the paper by Gubbins [25], in fact the incubation period distributions for the three genotypes without ARR allele (ARQ/ARQ, ARQ/VRQ, en VRQ/VRQ), are very similar. This motivates our choice to use the approximation of a genotype-independent test sensitivity. In case new data would provide evidence for a genotype dependent sensitivity of the tests used, this could of course be included in a revised analysis. We note that in the context of a rectal biopsy test evidence for such a dependence has been found [28].

FIS values and back-calculating genotype distributions

During the period 2003–2008 in which the flock cull data was gathered, selection for ARR alleles was ongoing in The Netherlands: E.g., from October 2004 onwards the use of ARR/ARR rams was obligatory for flocks with more than 10 breeding ewes (except some rare breeds), and from September 2005 until June 2007 the use of ARR/ARR rams was obligatory for all flocks (except some rare breeds) [15,29]. This forms a complication for our analysis because, if selective breeding takes place on an infected flock, genotype and allele frequencies found at the time of flock culling are not representative of the frequencies at the time that the scrapie infection became established in the flock. This latter “original” profile is the relevant one for our analysis of the basic reproduction number. To detect a possible history of selection for ARR alleles, we compare observed to expected heterozygosities for each flock. In the absence of selection genotype frequencies are expected to be in Hardy-Weinberg equilibrium. I.e. if the resistant allele R occurs in the population with frequency p and the susceptible allele S with frequency 1-p, the frequency of the homozygous genotype R/R is expected to be p2, the frequency of S/S to be (1-p)2 and the frequency of the heterozygous genotype R/S to be 2p(1-p). The deviation of Hardy-Weinberg equilibrium within a flock can be quantified by calculating the FIS value [30]: FIS = (Hobs − Hexp)/Hexp, with Hobs being the observed frequency of heterozygotes in the flock and Hexp the expected frequency of heterozygotes based on the allele frequencies in the flock assuming Hardy-Weinberg equilibrium. Negative FIS values (quantifying overrepresentation of heterozygote genotypes in comparison to expected frequencies when assuming Hardy-Weinberg equilibrium) indicate recent selective breeding.

In case of recent selective breeding we need to back-calculate from the current genotype profile the profile at the time that the current scrapie cases became infected. To this end, we assume that the original distribution is in Hardy-Weinberg equilibrium, and that the selective breeding was carried out by using ARR/ARR rams. In that case, as shown in S1 Text, we can derive the following analytical relationship between the original equilibrium value of the ARR allele frequency and the frequencies of resistant animals and of non-ARR carrying animals at the moment of culling: (4)

For flocks with negative FIS value we use Eq (4) to back-calculate. For consistency we also use equilibrium values for the genotype frequencies of flocks with a positive FIS value; in this case by calculating the Hardy-Weinberg equilibrium genotype frequencies corresponding to the allele frequencies at the moment of culling. We note that whereas these back-calculations do change the frequency of the group of non-ARR carrying (S/S) animals, the relative proportions of the different genotypes within this group remain unchanged. This implies that when testing the model assumption of flock-independent relative risks in non-ARR genotypes (as described above) we can work with the culled-flocks data before back-calculation.

Results

One example flock

To illustrate the analyses, we consider as an example a positive flock detected and culled in 2007. For this flock, 417 animals were genotyped and out of these, 200 were tested. Out of 21 ARQ*/VRQ, 55 ARQ*/ARQ*, 3 AHQ/VRQ, 19 AHQ/ARQ*, and 14 ARR/VRQ animals tested, 11 in total were found scrapie positive (including index case): 9 ARQ*/VRQ, 1 ARQ*/ARQ* and 1 AHQ/ARQ*. When we set pARQ*/VRQ, the scrapie risk for ARQ*/VRQ, equal to 9/21≈43% and multiply this with the relative risks of Table 1 to obtain expected case frequencies of the other genotypes, we find that for all these genotypes the actual case numbers are within the 95% probability range of a binomial model. For this flock the FIS value was negative (-0.08), possibly due to recent selective breeding. Based on the Hardy-Weinberg equilibrium genotype frequencies calculated using Eq (4), and setting Se to its best-fit value (see below) we estimated ρ0 = 8.3 and = 1.3 for this flock. We now describe the overall results.

Testing the model

We find that for 42 out of the 46 flocks with secondary cases and with at least one animal of ARQ*/VRQ genotype tested, it is possible to choose a value for pARQ*/VRQ,i such that for all genotypes the observed data are within the 95% probability range of the model using the relative susceptibility gγ estimated at population level. This analysis yielded values for pARQ*/VRQ ranging from 4% to 100% with a median of 43%. Amongst the four flocks outside the 95% probability range of the model, two flocks are outside the range due to having ARQ*/ARQ* positives but no ARQ*/VRQ positives, suggestive of an “ARQ” adapted scrapie strain with different relative susceptibilities. No evidence of a more widespread presence of an ARQ adapted strain was found: In only two of the 42 flocks consistent with the relative susceptibility model the detected prevalence in ARQ*/ARQ* exceeded that in ARQ*/VRQ; both had low numbers of tested ARQ*/VRQ animals. From this we conclude that the data is broadly consistent with our binomial model assuming flock-independent relative risk.

FIS values

In Fig 1 we show a histogram of the FIS values calculated for the genotype frequencies at the moment of culling for 49 flocks with at least one secondary scrapie case. Negative FIS values, representing a heterozygote frequency excess, occur for 33 flocks, indicating a recent history of selective breeding in these flocks.

thumbnail
Download:
Fig 1. FIS values histogram.

Histogram of the FIS values calculated for 49 culled flocks with at least one secondary scrapie case. Negative FIS values, representing a heterozygote frequency excess, occur for 33 flocks.

https://doi.org/10.1371/journal.pone.0139436.g001

Parameter estimation

The subset of flocks with at least two positive cases and at least eight tested animals of genotype ARQ*/VRQ contained 22 flocks, out of which 17 had a negative FIS value. In Fig 2 we plot a histogram of the estimated values for the transmission scale parameter ρ0 for these 22 flocks, based on a value for the sensitivity parameter of Se = 0.75. The corresponding histogram of values is presented as Fig 3. In Fig 2 we also plot a Weibull frequency distribution with mean and variance equal to the mean and variance of the ρ0 estimates. In Fig 4 we show how mean and variance of the distribution of ρ0 values change with the value assumed for Se. We observe that both the mean and the variance of the estimated ρ0 values are only weakly sensitive to the assumed Se, in other words that our results are robust against the uncertainty in Se.

thumbnail
Download:
Fig 2. Scale parameter distribution.

Histogram of 22 estimated values for the transmission scale parameter ρ0. The dashed line is a Weibull frequency distribution with mean and variance equal to the mean and variance of the ρ0 estimates for Se = 0.75.

https://doi.org/10.1371/journal.pone.0139436.g002

thumbnail
Download:
Fig 3. Distribution of the within-flock reproduction number.

Histogram of 22 estimated values for the reproduction number .

https://doi.org/10.1371/journal.pone.0139436.g003

thumbnail
Download:
Fig 4. Dependence on sensitivity parameter.

Mean (left-hand panel) and coefficient of variation (right-hand panel) of estimated ρ0 values as a function of the sensitivity Se.

https://doi.org/10.1371/journal.pone.0139436.g004

In Fig 5 we illustrate how this robustness arises by plotting the expected within-flock infection prevalence when Se = 1 against the “expected reproduction number” defined by . First, the results in Fig 5 display a threshold pattern close to , and this threshold is the main determinant of the estimated value. Second, we note that incorporating limitations in test sensitivity (Se<1) only affects the vertical scale and thus leaves the threshold pattern (and threshold location) unchanged. As a result, changing the value of Se has also little influence on the estimated .

thumbnail
Download:
Fig 5. Threshold pattern.

Expected within-flock infection prevalence against the “expected” within-flock reproduction number . The expected within-flock infection prevalence is computed by weighing the genotype-specific detected infection prevalence in tested animals by the frequency of the genotype in genotyped animals. As this computation uses detected prevalence it assumes Se = 1.

https://doi.org/10.1371/journal.pone.0139436.g005

Discussion

Based on genotyping and scrapie testing data from culled flocks we have developed a relatively simple transmission model for calculating the basic reproduction number for within-flock scrapie transmission in The Netherlands. The model expresses the basic reproduction number for within-flock scrapie as a weighted average of genotype-specific susceptibilities, multiplied by a single overall transmission parameter. The value of overall transmission scale parameter is allowed to vary between flocks to account for random between-flock variation in non-genetic determinants such as management practice. Indeed, our estimation for the overall transmission parameter yields a distribution that spans across a substantial range (of about one order of magnitude). Risk factor analyses in the literature have found evidence for certain aspects of flock management practice promoting scrape incidence in affected flocks, namely lambing in group pens [31], and spreading sheep compost on the land and disposing of the placenta in the compost [32]. In principle, the variation observed in the overall transmission parameter may in part be due to certain simplifying modelling assumptions made; in particular, due to the assumption of genotype-independent relative infectiousness. We further note that (unknown) age distribution differences between flocks could in reality produce between-flock differences in test sensitivity (as this sensitivity depends on age since infection, which is strongly correlated with age itself). However, given the robustness of the mean and variance of the distribution of estimated values for the overall transmission parameter, we expect that our results would be only weakly affected by any such between-flock differences.

On the basis of only the flock-level genotype profile (i.e. in absence of prevalence data) our model, due to the random between-flock variation in the overall transmission scale parameter shown in Fig 2, provides no single estimate. However, the model enables the construction of an value distribution for a large population of flocks (e.g. a national population of flocks) from a sufficiently large data set of genotype profiles. Such a distribution can serve as a building block for a population-level scrapie transmission model that incorporates both within- and between-flock transmission. We note that some of the variation in the transmission scale parameter will arise from stochastic variation in the transmission dynamics, which is expected to yield variation in observed prevalence even between flocks that had identical parameter values.

The model was parameterized to apply to the classical scrapie strain(s) dominating in the Netherlands, to which the animals with a VRQ allele are most susceptible. In a small subset of culled flocks however, the prevalence pattern indicated the presence of an ARQ-adapted scrapie strain. For such a strain, the susceptibility pattern across non-ARR alleles is different [33], leading to a different dependence of the overall susceptibility of a flock on the genotype profile. A similar model as the one developed here could be developed for countries/breeds dominated by ARQ-adapted scrapie strain if a set of relative susceptibility parameter values is available.

Our model neglected scrapie-related mortality, and assumed that animal replacement rate did not differ between genotypes. Also we used a genotype-independent parameter for the sensitivity of the rapid test, motivated by age-at-onset results by Gubbins [25] that indicate that the incubation period distributions for the three genotypes without ARR allele (ARQ/ARQ, ARQ/VRQ, and VRQ/VRQ) are very similar. Inaccuracies introduced by these approximations would cause additional widening of the estimated distribution of values for the transmission scale parameter ρ0.

Author Contributions

Conceived and designed the experiments: TJH JJW. Analyzed the data: TJH JJW. Wrote the paper: TJH JJW. Designed and carried out the analysis: TJH JJW. Drafted the manuscript: TJH. Commented on the manuscript: JJW.

References

  1. 1. Baylis M, Chihota C, Stevenson E, Goldmann W, Smith A, Sivam K, et al. (2004) Risk of scrapie in British sheep of different prion protein genotype. J Gen Virol 85: 2735–2740. pmid:15302967
  2. 2. Detwiler LA, Baylis M (2003) The epidemiology of scrapie. Revue Scientifique Et Technique-Office International Des Epizooties 22: 121–143.
  3. 3. Meydan H, Ozkan MM, Yildiz MA, Goldmann W (2013) Novel polymorphisms in ovine prion protein gene. Animal Genetics 44: 588–591. pmid:23782171
  4. 4. Vaccari G, Scavia G, Sala M, Cosseddu G, Chiappini B, Conte M, et al. (2009) Protective effect of the AT(137)RQ and ARQK(176) PrP alleles against classical scrapie in Sarda breed sheep. Veterinary Research 40.
  5. 5. Moreno CR, Cosseddu GM, Schibler L, Roig A, Moazami-Goudarzi K, Andreoletti O, et al. (2008) Identification of new quantitative trait loci (other than the PRNP gene) modulating the scrapie incubation period in sheep. Genetics 179: 723–726. pmid:18493086
  6. 6. Vilas V, Bohning D, Kuhnert R (2008) A comparison of the active surveillance of scrapie in the European Union. Veterinary Research 39: 37. pmid:18307969
  7. 7. Arnold M, Meek C, Webb CR, Hoinville LJ (2002) Assessing the efficacy of a ram-genotyping programme to reduce susceptibility to scrapie in Great Britain. Prev Vet Med 56: 227–249. pmid:12441238
  8. 8. Gubbins S, Webb CR (2005) Simulation of the options for a national control programme to eradicate scrapie from Great Britain. Preventive Veterinary Medicine 69: 175–187. pmid:15907568
  9. 9. Gubbins S, Roden JA (2006) Breeding programmes for TSE resistance in British sheep—II. Assessing the impact on the prevalence and incidence of scrapie. Preventive Veterinary Medicine 73: 17–31. pmid:16169613
  10. 10. Kao RR, Gravenor MB, McLean AR (2001) Modelling the national scrapie eradication programme in the UK. Mathematical Biosciences 174: 61–76. pmid:11730857
  11. 11. Truscott JE, Ferguson NM (2009) Control of scrapie in the UK sheep population. Epidemiology and Infection 137: 775–786. pmid:18687157
  12. 12. Gubbins S, Touzeau S, Hagenaars TJ (2010) The role of mathematical modelling in understanding the epidemiology and control of sheep transmissible spongiform encephalopathies: a review. Veterinary Research 41.
  13. 13. McIntyre KM, Gubbins S, Goldmann W, Hunter N, Baylis M (2008) Epidemiological Characteristics of Classical Scrapie Outbreaks in 30 Sheep Flocks in the United Kingdom. Plos One 3.
  14. 14. Tongue SC, Pfeiffer DU, Shearn PD, Wilesmith JW (2009) PrP genotype: A flock-level risk factor for scrapie? Preventive Veterinary Medicine 92: 309–323. pmid:19783057
  15. 15. Hagenaars TJ, Melchior MB, Bossers A, Davidse A, Engel B, van Zijderveld FG (2010) Scrapie prevalence in sheep of susceptible genotype is declining in a population subject to breeding for resistance. Bmc Veterinary Research 6.
  16. 16. Hagenaars TJ, Donnelly CA, Ferguson NM, Anderson RM (2000) The transmission dynamics of the aetiological agent of scrapie in a sheep flock. Mathematical Biosciences 168: 117–135. pmid:11121561
  17. 17. Matthews L, Woolhouse MEJ, Hunter N (1999) The basic reproduction number for scrapie. Proceedings of the Royal Society of London Series B-Biological Sciences 266: 1085–1090.
  18. 18. Diekmann O, Heesterbeek H, Britton T (2012) Mathematical Tools for Understanding Infectious Disease Dynamics. Princeton: Princeton University Press. 502 p.
    • 19. van Keulen LJM, Schreuder BEC, Meloen RH, MooijHarkes G, Vromans MEW, Langeveld JP (1996) Immunohistochemical detection of prion protein in lymphoid tissues of sheep with natural scrapie. Journal of Clinical Microbiology 34: 1228–1231. pmid:8727908
    • 20. Nodelijk G, van Roermund HJW, van Keulen LJM, Engel B, Vellema P, Hagenaars TJ (2011) Breeding with resistant rams leads to rapid control of classical scrapie in affected sheep flocks. Veterinary Research 42.
    • 21. Hagenaars TJ, Ferguson NM, Donnelly CA, Anderson NM (2001) Persistence patterns of scrapie in a sheep flock. Epidemiology and Infection 127: 157–167. pmid:11561968
    • 22. Woolhouse MEJ, Matthews L, Coen P, Stringer SM, Foster JD, Hunter N (1999) Population dynamics of scrapie in a sheep flock. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences 354: 751–756.
    • 23. Bozzetta E, Nappi R, Crudeli S, Meloni D, Varello K, Loprevite D, et al. (2011) Comparative performance of three TSE rapid tests for surveillance in healthy sheep affected by scrapie. Journal of Virological Methods 173: 161–168. pmid:21256871
    • 24. van Keulen LJM, Bossers A, van Zijderveld F (2008) TSE pathogenesis in cattle and sheep. Veterinary Research 39.
    • 25. Gubbins S (2008) Prevalence of sheep infected with classical scrapie in Great Britain: integrating multiple sources of surveillance data for 2002. Journal of the Royal Society Interface 5: 1343–1351.
    • 26. Andreoletti O, Berthon P, Marc D, Sarradin P, Grosclaude J, van Keulen L, et al. (2000) Early accumulation of PrPSc in gut-associated lymphoid and nervous tissues of susceptible sheep from a Romanov flock with natural scrapie. Journal of General Virology 81: 3115–3126. pmid:11086143
    • 27. van Keulen LJM, Schreuder BEC, Vromans MEW, Langeveld JPM, Smits MA (2000) Pathogenesis of natural scrapie in sheep. Archives of Virology: 57–71. pmid:11214935
    • 28. Gonzalez L, Dagleish MP, Martin S, Dexter G, Steele P, Finlayson J, et al. (2008) Diagnosis of preclinical scrapie in live sheep by the immunohistochemical examination of rectal biopsies. Veterinary Record 162: 397–+.
    • 29. Melchior MB, Windig JJ, Hagenaars TJ, Bossers A, Davidse A, van Zijderveld FG (2010) Eradication of scrapie with selective breeding: are we nearly there? Bmc Veterinary Research 6.
    • 30. Wright S (1969) Evolution and the genetics of populations. V. II. The theory of gene frequencies. Chicago: Univ. Chicago Press. 511 p.
      • 31. Sivam SK, Baylis M, Gravenor MB, Gubbins S (2006) Descriptive analysis of the results of an anonymous postal survey of the occurence of scrapie in Great Britain in 2002. Veterinary Record 158: 501–506. pmid:16617040
      • 32. Healy AM, Hannon D, Morgan KL, Weavers E, Collins JD, Doherty ML (2004) A paired case-control study of risk factors for scrapie in Irish sheep flocks. Preventive Veterinary Medicine 64: 73–83. pmid:15325763
      • 33. Baylis M, Goldmann W (2004) The genetics of scrapie in sheep and goats. Current Molecular Medicine 4: 385–396. pmid:15354869