Skip to main content

Advertisement

Variation in the prion protein gene (PRNP) sequence of wild deer in Great Britain and mainland Europe

Article metrics

Abstract

Susceptibility to prion diseases is largely determined by the sequence of the prion protein gene (PRNP), which encodes the prion protein (PrP). The recent emergence of chronic wasting disease (CWD) in Europe has highlighted the need to investigate PRNP gene diversity in European deer species, to better predict their susceptibility to CWD. Here we report a large genetic survey of six British deer species, including red (Cervus elaphus), sika (Cervus nippon), roe (Capreolus capreolus), fallow (Dama dama), muntjac (Muntiacus reevesii), and Chinese water deer (Hydropotes inermis), which establishes PRNP haplotype and genotype frequencies. Two smaller data sets from red deer in Norway and the Czech Republic are also included for comparison. Overall red deer show the most PRNP variation, with non-synonymous/coding polymorphisms at codons 98, 168, 226 and 247, which vary markedly in frequency between different regions. Polymorphisms P168S and I247L were only found in Scottish and Czech populations, respectively. T98A was found in all populations except Norway and the south of England. Significant regional differences in genotype frequencies were observed within both British and European red deer populations. Other deer species showed less variation, particularly roe and fallow deer, in which identical PRNP gene sequences were found in all individuals analysed. Based on comparison with PRNP sequences of North American cervids affected by CWD and limited experimental challenge data, these results suggest that a high proportion of wild deer in Great Britain may be susceptible to CWD.

Introduction

Chronic wasting disease (CWD) is a prion disease of cervid species that is widespread in North America and has recently emerged in Europe. It belongs to a family of diseases termed transmissible spongiform encephalopathies (TSEs), which cause progressive and invariably fatal neurodegenerative disorders in humans and animals. CWD was first described in the 1960s as a wasting syndrome of captive mule deer (Odocoileus hemionus hemionus) and black-tailed deer (Odocoileus hemionus columbianus) in Colorado wildlife facilities [1]. It has since been identified in North America in wapiti (elk, Cervus canadensis) white-tailed deer (Odocoileus virginianus), moose (Alces alces) and a single captive red deer (Cervus elaphus) [2, 3]. It is currently the only animal prion disease recognised in both captive and wild populations, and efficient transmission has led to spread of the disease across North America. At the time of writing, CWD is present in at least 25 US states and three Canadian provinces [4], as well in South Korea, introduced via importation of infected wapiti [5].

Member States of the European Union (EU) and Norway were believed to be free of CWD based on a survey conducted between 2006 and 2009 on 3274 farmed and 10 049 wild deer [6], with the vast majority being red deer (10 110). However, the first European case of CWD, and the first natural case in a reindeer (Rangifer tarandus) was identified in Norway in March 2016 [7]. Subsequent culling of the affected reindeer population resulted in detection of another 18 infected animals, and increased surveillance in Scandinavia has identified further “atypical” cases in three moose and a single red deer in separated locations in Norway, and individual moose in Finland and Sweden. Although the origins of CWD in Scandinavia are uncertain, the wide geographic distribution of cases raises concern that disease may have been present and undetected for some time [8,9,10,11].

The development of prion diseases is associated with misfolding of the prion protein (PrPC) into a protease-resistant form (PrPSc). The sequence of the open reading frame (ORF) of the gene encoding the prion protein (PRNP) is strongly associated with susceptibility to prion diseases. This association has been exploited in breeding programmes such as the National Scrapie Plan in the UK, which selectively bred for disease resistant PRNP genotypes to reduce the incidence of classical scrapie in sheep [12, 13]. Cervid PRNP genetics have been widely studied in North America, where both natural and experimental infection of deer species has allowed for the identification of PRNP polymorphisms that are associated with reduced incidence of disease and/or slower disease progression. The first identified was M132L in wapiti, which encodes a methionine to leucine change at codon 132. Genotype surveys of CWD affected deer populations have produced conflicting evidence regarding the association of codon 132 variation with incidence of disease [14, 15], but CWD challenges of deer [16, 17] and transgenic mice [18] suggest that the 132L variant reduces susceptibility to infection. Other important PRNP polymorphisms associated with CWD susceptibility include S225F in mule deer, [19,20,21,22], Q95H, G96S and A116G in white-tailed deer [23,24,25,26,27,28]. Although many of these PRNP polymorphisms have been associated with reduced susceptibility to CWD, none of them appear to confer complete resistance.

Experimental transmissions of CWD provide some evidence of susceptibility in deer species found in Europe. CWD has been transmitted orally to red deer [29], reindeer [30] and Reeves’ muntjac deer [31, 32]. Fallow deer have been infected via intracerebral inoculation [33], but failed to become infected when co-grazed with infected mule deer [34]. However, due to small sample sizes in all these studies, it is difficult to draw associations between susceptibility to experimental CWD infection and host PRNP genotype. It is also not known if the PRNP genotypes of the experimental animals used in these North American studies are representative of the genetics of native wild European deer populations.

The extent of PRNP sequence variation in European deer species has not previously been studied in detail, and this genetic information will be critically important in estimating potential susceptibility of these populations to emerging CWD and informing risk assessment and control/surveillance strategies. Novel PRNP variants (not seen in North American cervids) that show evidence of association with reduced susceptibility or resistance to CWD have the potential to be used in breeding programmes for captive and farmed deer. The primary aims of this study were (i) to perform a comprehensive survey of the protein-coding sequence of the PRNP gene in the six major free-ranging or wild species of deer found in Great Britain (GB) and (ii) to compare the PRNP genotype distributions in British and European populations of red deer, as red deer are one of the most numerous wild species in both GB and mainland Europe, and are also an economically important game species due to revenue generated via hunting, tourism and venison production [35].

Materials and methods

Samples from British deer species

A total of 1003 samples collected from wild deer in GB (England, Scotland and Wales) were analysed. Archived DNA samples were provided by JP, SPE and SLS for red deer, sika deer, and red/sika hybrids, and by KB for roe deer. The procedures for sample collection and genomic DNA preparation for these archives have been described previously [36,37,38]. Additional samples were collected by stalkers in the British Deer Society to give representation of geographical areas and species not covered by the existing archives. These included samples collected during routine culling, and samples collected from culled deer at a meat processing facility. Samples consisted of small sections of ear tips, which were preserved in 20% DMSO/saturated NaCl solution during transportation/storage [39] For each sample, further information including the sex, estimated age, and postcode/grid coordinates of the location of the animal, were recorded. The distribution of sampling sites for each deer species are shown in Figure 1.

Figure 1
figure1

Maps to show the sampling locations of deer samples from Great Britain. Each point represents one sampling location with ≤ 28 animals sampled. A Red deer (red circles), B roe deer (green circles), C sika (yellow squares) and fallow (blue circles), D muntjac (orange circles) and Chinese water deer (dark blue squares), or both muntjac and Chinese water deer (pink triangles). E Red deer samples in Scotland separated into regions: South Scotland (1), Argyll (2), Central Scotland (3), North Highlands (4) and Hebrides (5)

In total, samples from 480 red deer were analysed, of which 388 were from the archived collection, and 92 samples were collected by BDS deer stalkers. Samples from 83 sika deer were analysed, 73 samples were archived DNA, and 10 stalker collected samples from the South of England. All samples collected by members of the BDS were identified by animal phenotype only, whereas archived DNA samples of red and sika deer were selected both on phenotypical appearance and microsatellite genotyping data [38, 40], whereby a Q value of Q = 1 indicated a pure red deer, Q = 0 indicated a pure sika deer. Nine sika-red hybrids were also included in this study and were classified as having 0.091 ≤ Q ≤ 0.772.

For analysis red deer populations were separated into the following regions: Southern England (mainly Exmoor National Park) (n = 78), Northern England (Lake District) (n = 32), Southern Scotland (n = 23), Argyll (n = 47), Central Highlands (n = 142), Northern Highlands (n = 114) and Hebrides (n = 44) (Figure 1).

A total of 297 samples from roe deer were analysed, 171 from a DNA archive [37], and 126 collected by BDS deer stalkers. Samples from 66 fallow deer and 41 muntjac deer were collected by BDS deer stalkers. A total of 27 Chinese Water deer samples were analysed, 3 from a DNA archive [41] and 24 supplied by BDS deer stalkers.

Samples from other European red deer

Genotype analysis was also performed on 50 wild red deer collected from six counties in Norway, and 66 archived DNA samples from wild red deer (n = 46) and wild sika (n = 20) collected in five regions of the Czech Republic [42, 43]. In central Europe, 2 distinct lineages of red deer are present. Of the 46 Czech red deer, 24 were of Western lineage, and 22 of Eastern lineage.

Extraction and purification of genomic DNA

Genomic DNA was extracted from small sections of ear-tip tissue, digested with proteinase K (PK) at 37 °C overnight, followed by standard phenol/chloroform purification, ethanol precipitation and re-suspension in 1XTE buffer, as described previously [44]. Norwegian samples were prepared from brain tissue using DNeasy® Blood & Tissue Kit (Qiagen, Oslo, Norway) according to manufacturer’s instructions.

PRNP gene amplification and sequencing: British and Czech samples

The ORF of cervid PRNP (771 bp) was amplified by PCR using AmpliTaq Gold 360 master mix (Thermo Fisher Scientific) and primer pairs − 143d (ATGGAATGTGAAGAACATTTATGACCTA) or − 213d (AGGTCAACTTTGTCCTTGGAGGAG) with + 139u (TAAGCGCCAAGGGTATTAGCAT) or AR2 (GCAAGAAATGAGACACCACCAC) for British and Czech red, roe and sika deer, − 213d and AR3 (ACCACTACAGGGCTGCAGGTA) for fallow and Chinese water deer, and − 213d and AR2 (GCAAGAAATGAGACACCACCAC) for muntjac deer. PCR conditions were 95 °C for 5 min followed by 40 cycles of 95 °C for 30 s, 61 °C for 30 s and 72 °C for 1 min, and a final 72 °C for 10 min. PCR fragments were sequenced by sanger sequencing using primer + 70u (GCTGCAGGTAGATACTCCCTC) and BigDye® reagents (Life Technologies, Paisley, UK) as recommended by the manufacturer. Sequence data were analysed using Chromas and DNAStar Lasergene version 14.

PRNP gene amplification and sequencing: Norwegian samples

The ORF was amplified using PCR primers Ce19_F (ATTTTGCAGATAAGTCATC) and Ce 778_R (AGAAGATAATGAAAACAGGAAG) [14]. PCR conditions were 95 °C for 2 min, followed by 36 cycles at 95 °C for 30 s, 51 °C for 30 s and 72 °C for 45 s, then a final cycle at 72 °C for 10 min. PCR fragments were directly sequenced using Ce19_F and Ce 778_R and BigDye® reagents (Life Technologies, Paisley, UK) as recommended by the manufacturer. Sequence data were analysed using Seqscape v3.0 software, Sequence Scanner v2.0 (Applied Biosystems) and MEGA7 version 7.0.26 [45].

Statistical analysis

The Chi square test was used to compare haplotype frequencies between different regions and countries.

Accession numbers

Sequences were deposited in Genbank with the following accession numbers: Capreolus capreolus MK103016; Dama dama MK103017; Cervus nippon MK103018, MK103019; Muntiacus reeveesi MK103020–MK103023; Hydropotes inermis MK103024–MK103026; Cervus elaphus with 247L MK103027.

Results

The ORF encoded by exon 3 of the PRNP gene is highly conserved among cervid species, and non-synonymous polymorphisms identified in individual species are usually reported as variations from a consensus sequence [25]. Following this convention, the amino acid substitutions resulting from PRNP sequence polymorphisms identified in the deer species examined in our survey are summarized in Table 1.

Table 1 Amino acid variation within the ORF of cervid PRNP in the surveyed deer species

British red deer

Four single nucleotide polymorphisms (SNPs) were identified in 480 British red deer samples: a synonymous SNP at codon 136 (t/c, nucleotide position 408), and three non-synonymous polymorphisms at codons 98 (a/g, nucleotide position 292), 168 (c/t, nucleotide position 502) and 226 (g/c, nucleotide position 676) giving rise to amino acid changes threonine to alanine (T98A), proline to serine (P168S) and glutamine to glutamic acid (Q226E), respectively. Linkage was found between positions 408 and 676 (codons 136 and 226), such that the haplotypes were either 408t-676c or 408c-676g. The SNPs resulting in substitutions at positions 98 and 168 occurred exclusively on a Q226 background and the following four haplotypes were inferred: T98-P168-E226 (TPE), T98-P168-Q226 (TPQ), A98-P168-Q226 (APQ) and A98-S168-Q226 (ASQ).

Norwegian and Czech red deer

Only the TPE and TPQ haplotypes were identified in Norwegian red deer. The synonymous codon 136 polymorphism was observed in the same haplotype linkage as described above. Czech red deer had TPE, TPQ and APQ haplotypes, as well as a novel non-synonymous SNP at codon 247 (a/c nucleotide position 739) leading to a conservative isoleucine to leucine substitution which occurred in linkage with APQ (APQ-L247).

British sika, Czech sika and sika-red deer hybrids

Scottish pure sika (n = 73), Czech sika (n = 20), English sika (n = 8/10) were all homozygous for the TPQ haplotype. Both British and Czech sika had a synonymous polymorphism at codon 133 (g/c, position 399) at frequencies of 0.08% and 0.1% respectively. Of the nine genetically determined sika-red hybrids from Argyll (0.091 ≤ Q ≤ 0.772) three were homozygous TPQ/TPQ, four were heterozygous TPE/TPQ, and two were heterozygous for the APQ haplotype (APQ/TPQ and APQ/TPE). The two animals heterozygous for the APQ haplotype were characterised as red deer-like (Q = 0.71, 0.74), as were two of four TPE/TPQ heterozygotes (Q = 0.77, 0.62). Of the 10 English animals identified as sika by phenotype only, two were TPE/TPQ heterozygous and were therefore possibly sika-like hybrids.

Roe and fallow deer

The PRNP ORF of British roe deer (n = 297) all encoded the haplotype described above as TPQ. Roe deer were monomorphic in PRNP sequence, with no variation in nucleotide sequence in any of the animals analysed. Likewise, fallow deer (n = 66) showed no variation in PRNP sequence and encoded the TPE haplotype with an additional amino acid substitution of asparagine for serine at codon 138 (TPE-N138).

Muntjac and Chinese water deer

The PRNP ORF of all 41 muntjac deer analysed encoded serine at position 98, allowing their PRNP haplotype to be described as S98P168Q226 (SPQ). In addition, synonymous polymorphisms were identified at nucleotide positions 15 (c/t), 126 (g/a) and 606 (c/t) at varying frequencies. The PRNP coding region of 27 Chinese water deer revealed two amino acid sequence polymorphisms; a serine to asparagine substitution at codon 100 (g/a, nucleotide position 299) and a 24 bp deletion of the octapeptide repeat (ccccatggaggtggctggggtcag) which corresponds to the in-frame loss of the amino acid sequence PHGGGWGQ. In most species, a glycine-rich peptide sequence (octapeptide repeat) is present in five consecutive copies in the N terminal region of wild type PrPC, whilst Chinese water deer had a four octapeptide repeat haplotype. From the genotypes it was inferred that the deletion was linked to 100S (100S-Δ71–78). The haplotype frequencies were 0.21 for 100S, 0.25 for 100SΔ71–78, and 0.54 for 100 N. Seven animals were heterozygous for 100SΔ71–78, whilst three were homozygous S100Δ71–78/S100Δ71–78. The most common genotype in this sampling was N100/N100 (0.35).

Red deer genotype and haplotype frequencies and distributions

Red deer showed the greatest variation in PRNP sequence amongst the six deer species surveyed, with four haplotypes identified in the British population. Regional variation was evident within Great Britain in the distribution and frequency of the four PRNP haplotypes and resulting genotypes (Tables 2 and 3). For the TPQ and TPE haplotypes, frequencies in different regions of Scotland and England ranged from 0.11 to 0.62 and 0.31 to 0.88 respectively. The difference was most striking in Southern England, which had a very high frequency (0.88) of the TPE haplotype, whereas in Northern England, the TPQ haplotype was predominant (0.67). The APQ and ASQ haplotypes were found almost exclusively in Scottish red deer, with only one heterozygous TPE/APQ individual identified in Northern England. The frequencies of APQ and ASQ haplotypes showed regional variation within Scotland, with the highest frequencies seen in the Central Highlands, Argyll and Southern Scotland (Table 3). When compared using the Chi squared test, the Central Highlands, Lowlands and Argyll were significantly different in haplotype frequency when compared to the Northern Highlands and Hebrides (P < 0.05), as were the overall haplotype frequencies of England and Scotland (P < 0.05).

Table 2 Genotype frequencies of PRNP polymorphisms in British red deer populations
Table 3 Haplotype frequencies of PRNP polymorphisms in British red deer populations

The haplotype frequencies in GB were also compared by Chi squared test with those in red deer from Norway and the Czech Republic, as well as previously published surveys of Scotland, Italy and Spain [46, 47] (Table 4). The number of samples was too low to assess the extent of regional variation within Norway and the Czech Republic. In Norway, where only the TPQ and TPE haplotypes were identified, haplotype frequencies were very similar to those found in Southern England, i.e. 0.89 for TPE and 0.11 for TPQ. In the Czech Republic, frequencies of the TPQ and TPE haplotypes were approximately equal (0.35 and 0.44, respectively), while the APQ and APQ-L247 haplotypes were present at frequencies of 0.16 and 0.05, respectively. The APQ-L147 haplotype was only seen in the Czech red deer belonging to the Western lineage (n = 24) [43]. The haplotype frequencies presented in Table 4 were compared across countries, and whilst most differences between countries were statistically significant (P < 0.05), those in Spain and the Czech Republic were not.

Table 4 Comparison of haplotype frequencies across European red deer populations

Discussion

This study represents the largest survey to date of PRNP genetic variation in the wild (free-ranging) deer population of any European country and is the first to include free-ranging muntjac and Chinese water deer. The analysis included over 1003 animals from six different species of deer present in Britain, as well as red deer and sika deer from mainland Europe. Ten haplotypes were identified, five of which were present in red deer. The majority of polymorphisms identified in the ORF (exon 3) of the PRNP gene have been reported previously [32, 46,47,48,49], apart from the I247L polymorphism in red deer from the Czech Republic, and the octapeptide repeat deletion in Chinese water deer. Deletion of the third glycine-rich octapeptide repeat is not an uncommon variation, having previously been described in chamois [47], lions [50], and some primates [51]. This is the first time a functional octapeptide deletion variant has been described in deer.

Red deer showed the greatest variation in PRNP sequence, with four haplotypes identified in Great Britain, designated TPQ, TPE, APQ and ASQ. The TPQ, TPE and APQ haplotypes have been previously reported in surveys of Italian, Scottish (as haplotypes 1, 8 and 10) [46] and Spanish red deer [47]. The ASQ haplotype was not identified in Italy and Spain, and only reported previously in a single Scottish red deer [46]. Non-synonymous polymorphisms previously identified in Italian red deer at codons 59 (G59S) and 208 (M208I) in single individuals [46, 52], as well as low frequency synonymous polymorphisms at codons 15, 21, 59, 78 and 79 [46] were not present in our surveyed populations.

Taking our results together with previously published surveys, the total number of European red deer analysed for PRNP sequence variation stands at 1108 individuals. In every region or country surveyed, TPQ and TPE haplotypes were present at the highest frequencies (Tables 3 and 4). In comparison to other European red deer populations, there was a high frequency of the TPE haplotype in South England (predominantly the Exmoor National Park) and in Norway, which is possibly due to historical population declines of red deer with refugial populations in these areas [53]. The APQ haplotype was absent from both Southern England and Norway but was found at similar frequencies in Scottish (0.16), Czech (0.16), Spanish (0.21) and Italian (0.10) red deer [46, 47]. In our study, it was evident that the frequency of the APQ haplotype showed marked variation between different Scottish regions (Table 3). This may explain why our frequency of the APQ haplotype was higher (0.16) than the previous Scottish survey (0.05); in the previous study, samples were collected from the Isle of Rum and an unspecified location on mainland Scotland, which may not have coincided with regions where we identified the highest APQ haplotype frequencies. The ASQ and APQ-L247 haplotypes were present at frequencies of 0.04 and 0.05 in Scottish and Czech red deer populations respectively and may be unique to these populations. Our results emphasize that there are large regional variations in PRNP genotype frequencies in British red deer populations, and similar variation in genotype frequencies between the red deer populations surveyed in European countries to date. Regional variations in the genetics of Scottish red deer populations have been previously reported, in part due to natural barriers [36] and in part due to introductions of non-native red deer into Scotland from England and further afield [54].

In contrast to red deer, our roe deer samples showed no variation in their PRNP sequence. Whilst the number of analysed animals across studies stands now at 541, only one synonymous change at codon 24 has previously been reported in two individuals [48]. Considering that roe deer are now the most widespread species in the UK, it is possible that this lack of diversity may be due to a population bottleneck following the near-extinction of British roe deer prior to the 1800 s [37]. However, a similar lack of PRNP diversity was also evident in European roe deer populations from Italy [46], Scandinavia [48] and Spain [47]. Similarly, we found little evidence of sequence variation in sika, fallow and muntjac deer, with the result that each of these species is predicted to express a single specific variant of the PrP protein. However, with a relatively small number of animals analysed (103, 92 and 41 respectively), less frequent polymorphisms may still be found in larger surveys. The PRNP sequence we identified in British fallow deer (described here as TPE-N138), was identical to previously published sequences from Scandinavia [48] and Spain [47], as well as fallow deer previously used in experimental challenge experiments [34, 55]. Our data support the conclusion that the N138 change is fixed in the fallow deer sequence and not polymorphic. The three haplotypes identified in 27 wild Chinese water deer is more than might be expected considering that the UK population was established relatively recently from small numbers of escapees from captive collections.

Comparison of the PRNP genotypes in European deer with those in North American cervid species with European data can be used to help predict their potential susceptibility to CWD, although confirmation of these predictions under natural or experimental challenge is essential. The TPQ haplotype, which we identified at high frequencies in British red (0.29), Roe (1.0), Sika (1.0) and Chinese Water deer (0.21), is identical to PRNP sequences found in North American cervids that are highly susceptible to CWD. Our analysis, and published data to date for European deer species, has not identified any PRNP polymorphisms that are associated with reduced CWD susceptibility in North American cervids, e.g. M132L, G96S [25]. This suggests that a relatively high proportion of free-ranging/wild deer in Britain and Europe may be susceptible to CWD, or at least to the predominant strains of CWD in North America.

The association of non-synonymous PRNP polymorphisms found only in British/European deer with susceptibility to CWD has yet to be determined. Amino acid substitutions at codon 98 are found in red (T98A) and muntjac (T98S) deer, and nearby codons 95 and 96 play a role in CWD susceptibility in white-tailed deer [24]. However, muntjac deer with identical PRNP sequences can be experimentally infected with North American CWD isolates [31] and recently camels, which encode alanine at the equivalent position, were reported to develop natural prion disease [56].

The substitution of serine for proline at codon 168 in red deer is of particular interest, as this position is located in the β2-α2 loop of the prion protein, an area critical for protein structure and stability [57]. The murine prion protein with a serine substitution at the equivalent residue showed no difference to wildtype in efficiency of conversion to PrPSc in an in vitro cell-free conversion assay [57]. However, in sheep the substitution of leucine for proline at codon 168 was associated with very prolonged incubation periods for experimental BSE infection [58], and a reduction of cell-free conversion efficiencies [59]. The ASQ haplotype seen in red deer has not been tested in vitro or by prion challenge, so any potential association with prion disease is still unresolved.

The novel PRNP haplotype TPQ-L247 identified in Czech red deer may be less likely to influence susceptibility as it represents a conservative amino acid substitution within the C-terminal signal peptide which is cleaved from the mature prion protein [60], however it is possible that it may have unknown effects that influence pathogenesis.

In North America codon 226 encodes either E (e.g. wapiti) or Q (e.g. mule deer and white-tailed deer), whereas in European red deer both variants are present and heterozygous animals are common. Codon 226 is in the third α-helix of PrPC, a part of the protein involved in pathogenic processes during prion disease. Nearby polymorphisms include S225F, which is associated with reduced CWD susceptibility in mule deer [20], and Q222K, which is strongly associated with prion disease resistance in goats [61, 62]. Studies in transgenic mice expressing cervid PrP with Q226 or E226 suggest that variation at 226 may contribute to differences in disease progression and susceptibility to certain prion strains [63, 64]. Following experimental oral transmission of CWD to four deer of EE226, QQ226 and EQ 226 genotypes, all animals became infected [65]. The small sample size of this study does not allow a definite conclusion regarding the association of codon 226 with CWD susceptibility, or the effect of heterozygosity at this position.

The TPQ-N100 haplotype identified in Chinese water deer has not been reported in any other cervid but N100 is common to other mammalian prion proteins including bank vole, a species highly susceptible to CWD challenges [66]. The fact that fallow deer failed to become infected with CWD when exposed to infected mule deer for up to 7 years suggests that the substitution of asparagine at codon 138 may be associated with relative resistance to infection, but it is possible that other species-specific genetic factors also play a role [34].

In summary, this study has provided an analysis of PRNP sequence variation in the six major deer species found in the wild in Great Britain and adds to our knowledge of PRNP variation in European deer species in general. Based on a comparison of the sequences with those found in North American cervids affected by CWD, it appears likely that a large proportion of the British deer population would be susceptible to CWD. The effect of specific PRNP haplotypes only found in European species, such as ASQ in Scottish red deer, on CWD susceptibility are still to be determined.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CWD:

chronic wasting disease

PRNP :

prion protein gene

SNP:

single nucleotide polymorphism

EU:

European Union

GB:

Great Britain

ORF:

open reading frame

References

  1. 1.

    Williams ES, Young S (1980) Chronic wasting disease of captive mule deer: a spongiform encephalopathy. J Wildl Dis 16:89–98

  2. 2.

    Williams ES, Miller MW (2002) Chronic wasting disease in deer and elk in North America. Rev Sci Tech 21:305–316

  3. 3.

    Schwabenlander MD, Culhane MR, Hall SM, Goyal SM, Anderson PL, Carstensen M, Wells SJ, Slade WB, Armién AG (2013) A case of chronic wasting disease in a captive red deer (Cervus elaphus). J Vet Diagn Invest 25:573–576

  4. 4.

    USGS (2018) Distribution of chronic wasting disease in North America Sept 2018. https://www.usgs.gov/media/images/distribution-chronic-wasting-disease-north-america-sept-2018. Accessed 22 Oct 2018

  5. 5.

    Lee Y-H, Sohn H-J, Kim M-J, Kim H-J, Lee W-Y, Yun E-I, Tark D-S, Cho I-S, Balachandran A (2013) Strain characterization of the Korean CWD cases in 2001 and 2004. J Vet Med Sci 75:95–98

  6. 6.

    EFSA (2010) Scientific opinion on the results of the EU survey for Chronic Wasting Disease (CWD) in cervids. EFSA J 8:1861

  7. 7.

    Benestad SL, Mitchell G, Simmons M, Ytrehus B, Vikoren T (2016) First case of chronic wasting disease in Europe in a Norwegian free-ranging reindeer. Vet Res 47:88

  8. 8.

    Pirisinu L, Tran L, Chiappini B, Vanni I, Di Bari MA, Vaccari G, Vikøren T, Madslien KI, Våge J, Spraker T, Mitchell G, Balachandran A, Baron T, Casalone C, Rolandsen CM, Røed KH, Agrimi U, Nonno R, Benestad SL (2018) Novel type of chronic wasting disease detected in moose (Alces alces), Norway. Emerg Infect Dis 24:2210–2218

  9. 9.

    Våge J, Hopp P, Vikoren T, Madslien K, Tarpai A, Moldal T, Benestad S (2018) The surveillance programme for Chronic Wasting Disease (CWD) in free ranging and captive cervids in Norway 2017 https://www.vetinst.no/overvaking/chronic-wasting-disease-vilt-cwd/_/attachment/. Accessed 20 Mar 2019

  10. 10.

    Evira (2018) Moose Found Dead in Forest with Chronic Wasting Disease https://mmm.fi/en/article/-/asset_publisher/evira-metsaan-kuolleella-hirvella-naivetystauti. Accessed 9 July 2019

  11. 11.

    SVA (2019) Wasting Disease (CWD) discovered on moose in Norrbotten County. https://www.sva.se/om-sva/pressrum/nyheter-fran-sva/avmagringssjuka-cwd-upptackt-pa-alg-i-norrbottens-lan. Accessed 1 Apr 2019 (in Swedish)

  12. 12.

    Ortiz-Pelaez A, Thompson CE, Dawson M, Hunter N (2014) The impact of the National Scrapie Plan on the PRNP genotype distribution of the British national flock, 2002–2012. Vet Rec 174:530

  13. 13.

    Dawson M, Del Rio Vilas V (2008) Control of classical scrapie in Great Britain. Practice 30:330–333

  14. 14.

    O’Rourke KI, Besser TE, Miller MW, Cline TF, Spraker TR, Jenny AL, Wild MA, Zebarth GL, Williams ES (1999) PrP genotypes of captive and free-ranging Rocky Mountain elk (Cervus elaphus nelsoni) with chronic wasting disease. J Gen Virol 80:2765–2769

  15. 15.

    Perucchini M, Griffin K, Miller MW, Goldmann W (2008) PrP genotypes of free-ranging wapiti (Cervus elaphus nelsoni) with chronic wasting disease. J Gen Virol 89:1324–1328

  16. 16.

    Hamir AN, Gidlewski T, Spraker TR, Miller JM, Creekmore L, Crocheck M, Cline T, O’Rourke KI (2006) Preliminary observations of genetic susceptibility of elk (Cervus elaphus nelsoni) to chronic wasting disease by experimental oral inoculation. J Vet Diagn Invest 18:110–114

  17. 17.

    O’Rourke KI, Spraker TR, Zhuang D, Greenlee JJ, Gidlewski TE, Hamir AN (2007) Elk with a long incubation prion disease phenotype have a unique PrP d profile. NeuroReport 18:1935–1938

  18. 18.

    Green KM, Browning SR, Seward TS, Jewell JE, Ross DL, Green MA, Williams ES, Hoover EA, Telling GC (2008) The elk PRNP codon 132 polymorphism controls cervid and scrapie prion propagation. J Gen Virol 89:598–608

  19. 19.

    Brayton KA, O’Rourke KI, Lyda AK, Miller MW, Knowles DP (2003) A processed pseudogene contributes to apparent mule deer prion gene heterogeneity. Gene 326:167–173

  20. 20.

    Jewell JE, Conner MM, Wolfe LL, Miller MW, Williams ES (2005) Low frequency of PrP genotype 225SF among free-ranging mule deer (Odocoileus hemionus) with chronic wasting disease. J Gen Virol 86:2127–2134

  21. 21.

    Fox KA, Jewell JE, Williams ES, Miller MW (2006) Patterns of PrPCWD accumulation during the course of chronic wasting disease infection in orally inoculated mule deer (Odocoileus hemionus). J Gen Virol 87:3451–3461

  22. 22.

    Wolfe LL, Fox KA, Miller MW (2014) “Atypical” chronic wasting disease in PRNP genotype 225FF mule deer. J Wildl Dis 50:660–665

  23. 23.

    O’Rourke KI, Spraker TR, Hamburg LK, Besser TE, Brayton KA, Knowles DP (2004) Polymorphisms in the prion precursor functional gene but not the pseudogene are associated with susceptibility to chronic wasting disease in white-tailed deer. J Gen Virol 85:1339–1346

  24. 24.

    Johnson C, Johnson J, Vanderloo JP, Keane D, Aiken JM, McKenzie D (2006) Prion protein polymorphisms in white-tailed deer influence susceptibility to chronic wasting disease. J Gen Virol 87:2109–2114

  25. 25.

    Robinson SJ, Samuel MD, O’Rourke KI, Johnson CJ (2012) The role of genetics in chronic wasting disease of North American cervids. Prion 6:153–162

  26. 26.

    Johnson CJ, Herbst A, Duque-Velasquez C, Vanderloo JP, Bochsler P, Chappell R, McKenzie D (2011) Prion protein polymorphisms affect chronic wasting disease progression. PLoS One 6:e17450

  27. 27.

    Keane DP, Barr DJ, Bochsler PN, Hall SM, Gidlewski T, O’Rourke KI, Spraker TR, Samuel MD (2008) Chronic wasting disease in a Wisconsin white-tailed deer farm. J Vet Diagn Invest 20:698–703

  28. 28.

    Otero A, Duque Velásquez C, Johnson C, Herbst A, Bolea R, Badiola JJ, Aiken J, McKenzie D (2019) Prion protein polymorphisms associated with reduced CWD susceptibility limit peripheral PrPCWD deposition in orally infected white-tailed deer. BMC Vet Res 15:50

  29. 29.

    Martin S, Jeffrey M, González L, Sisó S, Reid HW, Steele P, Dagleish MP, Stack MJ, Chaplin MJ, Balachandran A, Gonzalez L, Siso S, Reid HW, Steele P, Dagleish MP, Stack MJ, Chaplin MJ, Balachandran A (2009) Immunohistochemical and biochemical characteristics of BSE and CWD in experimentally infected European red deer (Cervus elaphus elaphus). BMC Vet Res 5:26

  30. 30.

    Mitchell GB, Sigurdson CJ, O’Rourke KI, Algire J, Harrington NP, Walther I, Spraker TR, Balachandran A (2012) Experimental oral transmission of chronic wasting disease to reindeer (Rangifer tarandus tarandus). PLoS One 7:e39055

  31. 31.

    Nalls AV, McNulty E, Powers J, Seelig DM, Hoover C, Haley NJ, Hayes-Klug J, Anderson K, Stewart P, Goldmann W, Hoover EA, Mathiason CK (2013) Mother to offspring transmission of chronic wasting disease in reeves’ muntjac deer. PLoS One 8:e71844

  32. 32.

    Nalls AV, McNulty E, Hoover CE, Pulscher LA, Hoover EA, Mathiason CK (2017) Infectious prions in the pregnancy microenvironment of chronic wasting disease-infected reeves’ muntjac deer. J Virol 91:e00501–e00517

  33. 33.

    Hamir AN, Greenlee JJ, Nicholson EM, Kunkle RA, Richt JA, Miller JM, Hall M (2011) Experimental transmission of chronic wasting disease (CWD) from elk and white-tailed deer to fallow deer by intracerebral route: final report. Can J Vet Res 75:152–156

  34. 34.

    Rhyan JC, Miller MW, Spraker TR, McCollum M, Nol P, Wolfe LL, Davis TR, Creekmore L, O’Rourke KI (2011) Failure of fallow deer (Dama dama) to develop chronic wasting disease when exposed to a contaminated environment and infected mule deer (Odocoileus hemionus). J Wildl Dis 47:739–744

  35. 35.

    Pérez-Espona S, Pemberton JM, Putman R (2009) Red and sika deer in the British Isles, current management issues and management policy. Mamm Biol 74:247–262

  36. 36.

    Pérez-Espona S, Pérez-Barbería FJ, Mcleod JE, Jiggins CD, Gordon IJ, Pemberton JM (2008) Landscape features affect gene flow of Scottish Highland red deer (Cervus elaphus). Mol Ecol 17:981–996

  37. 37.

    Baker KH, Hoelzel A (2013) Evolution of population genetic structure of the British roe deer by natural and anthropogenic processes (Capreolus capreolus). Ecol Evol 3:89–102

  38. 38.

    Smith SL, Senn HV, Pérez-Espona S, Wyman MT, Heap E, Pemberton JM (2018) Introgression of exotic Cervus (nippon and canadensis) into red deer (Cervus elaphus) populations in Scotland and the English Lake District. Ecol Evol 8:2122–2134

  39. 39.

    Amos B, Hoelzel AR (1991) Long-term preservation of whale skin for DNA analysis. In: Genetic ecology of whales and dolphins. Report of the International Whaling Commission. pp 99–103

  40. 40.

    Senn HV, Pemberton JM (2009) Variable extent of hybridization between invasive sika (Cervus nippon) and native red deer (C. elaphus) in a small geographical area. Mol Ecol 18:862–876

  41. 41.

    Esin A, Bergendahl LT, Savolainen V, Marsh JA, Warnecke T (2018) The genetic basis and evolution of red blood cell sickling in deer. Nat Ecol Evol 2:367–376

  42. 42.

    Krojerova-Prokesova J, Barančeková M, Voloshina I, Myslenkov A, Lamka J, Koubek P (2013) Dybowski’s sika deer (Cervus nippon hortulorum): genetic divergence between natural primorian and introduced Czech populations. J Hered 104:312–326

  43. 43.

    Krojerova-Prokesova J, Barančeková M, Koubek P (2015) Admixture of eastern and western European red deer lineages as a result of postglacial recolonization of the Czech Republic (Central Europe). J Hered 106:375–385

  44. 44.

    Tan BC, Blanco ARA, Fiona Houston E, Stewart P, Goldmann W, Gill AC, de Wolf C, Manson JC, McCutcheon S (2012) Significant differences in incubation times in sheep infected with bovine spongiform encephalopathy result from variation at codon 141 in the PRNP gene. J Gen Virol 93:2749–2756

  45. 45.

    Kumar S, Stecher G, Tamura K (2016) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol 33:1870–1874

  46. 46.

    Peletto S, Perucchini M, Acín C, Dalgleish MP, Reid HW, Rasero R, Sacchi P, Stewart P, Caramelli M, Ferroglio E, Bozzetta E, Meloni D, Orusa R, Robetto S, Gennero S, Goldmann W, Acutis PL (2009) Genetic variability of the prion protein gene (PRNP) in wild ruminants from Italy and Scotland. J Vet Sci 10:115–120

  47. 47.

    Pitarch JL, Raksa HC, Arnal MC, Revilla M, Martínez D, Fernández De Luco D, Badiola JJ, Goldmann W, Acín C (2018) Low sequence diversity of the prion protein gene (PRNP) in wild deer and goat species from Spain. Vet Res 49:33

  48. 48.

    Wik L, Mikko S, Klingeborn M, Stéen M, Simonsson M, Linné T (2012) Polymorphisms and variants in the prion protein sequence of European moose (Alces alces), reindeer (Rangifer tarandus), roe deer (Capreolus capreolus) and fallow deer (Dama dama) in Scandinavia. Prion 6:256–260

  49. 49.

    Jeong H-J, Lee J-B, Park S-Y, Song C-S, Kim B-S, Rho J-R, Yoo M-H, Jeong B-H, Kim Y-S, Choi I-S (2009) Single-nucleotide polymorphisms in prion protein gene of the Korean subspecies of Chinese water deer (Hydropotes inermis argyropus). Korean J Vet Res 49:59–62

  50. 50.

    Stewart P, Campbell L, Skogtvedt S, Griffin KA, Arnemo JM, Tryland M, Girling S, Miller MW, Tranulis MA, Goldmann W (2012) Genetic predictions of prion disease susceptibility in carnivore species based on variability of the prion gene coding region. PLoS One 7:e50623

  51. 51.

    Schätzl HM, Da Costa M, Taylor L, Cohen FE, Prusiner SB (1995) Prion protein gene variation among primates. J Mol Biol 245:362–374

  52. 52.

    Kaluz S, Kaluzova M, Flint APF (1997) Sequencing analysis of prion genes from red deer and camel. Gene 199:283–286

  53. 53.

    Haanes H, Røed KH, Flagstad Rosef O (2010) Genetic structure in an expanding cervid population after population reduction. Conserv Genet 11:11–20

  54. 54.

    Pérez-Espona S, Hall RJ, Pérez-Barbería FJ, Glass BC, Ward JF, Pemberton JM (2013) The impact of past introductions on an iconic and economically important species, the red deer of Scotland. J Hered 104:14–22

  55. 55.

    Hamir AN, Kunkle RA, Nicholson EM, Miller JM, Hall SM, Schoenenbruecher H, Brunelle BW, Richt JA (2008) Preliminary observations on the experimental transmission of chronic wasting disease (CWD) from elk and white-tailed deer to fallow deer. J Comp Pathol 138:121–130

  56. 56.

    Babelhadj B, Di Bari MA, Pirisinu L, Chiappini B, Gaouar SBS, Riccardi G, Marcon S, Agrimi U, Nonno R, Vaccari G (2018) Prion disease in dromedary camels, Algeria. Emerg Infect Dis 24:1029–1036

  57. 57.

    Kirby L, Agarwal S, Graham JF, Goldmann W, Gill AC (2010) Inverse correlation of thermal lability and conversion efficiency for five prion protein polymorphic variants. Biochemistry 49:1448–1459

  58. 58.

    Goldmann W, Houston F, Stewart P, Perucchini M, Foster J, Hunter N (2006) Ovine prion protein variant A136R154L168Q171 increases resistance to experimental challenge with bovine spongiform encephalopathy agent. J Gen Virol 87:3741–3745

  59. 59.

    Kirby L, Goldmann W, Houston F, Gill AC, Manson JC (2006) A novel, resistance-linked ovine PrP variant and its equivalent mouse variant modulate the in vitro cell-free conversion of rPrP to PrPres. J Gen Virol 87:3747–3751

  60. 60.

    Gill AC, Castle AR (2018) The cellular and pathologic prion protein. Handb Clin Neurol 153:21–44

  61. 61.

    Lacroux C, Perrin-Chauvineau C, Corbière F, Aron N, Aguilar-Calvo P, Torres JM, Costes P, Brémaud I, Lugan S, Schelcher F, Barillet F, Andréoletti O (2014) Genetic resistance to scrapie infection in experimentally challenged goats. J Virol 88:2406–2413

  62. 62.

    Aguilar-Calvo P, Fast C, Tauscher K, Espinosa J-C, Groschup MH, Nadeem M, Goldmann W, Langeveld J, Bossers A, Andreoletti O, Torres J-M (2015) Effect of Q 211 and K 222 PRNP polymorphic variants in the susceptibility of goats to oral infection with goat bovine spongiform encephalopathy. J Infect Dis 212:664–672

  63. 63.

    Angers RC, Seward TS, Napier D, Green M, Hoover E, Spraker T, O’Rourke K, Balachandran A, Telling GC (2009) Chronic wasting disease prions in elk antler velvet. Emerg Infect Dis 15:696–703

  64. 64.

    Angers RC, Kang H-E, Napier D, Browning S, Seward T, Mathiason C, Balachandran A, McKenzie D, Castilla J, Soto C, Jewell J, Graham C, Hoover EA, Telling GC (2010) Prion strain mutation determined by prion protein conformational compatibility and primary structure. Science 328:1154–1158

  65. 65.

    Balachandran A, Harrington NP, Algire J, Soutyrine A, Spraker TR, Jeffrey M, González L, O’Rourke KI (2010) Experimental oral transmission of chronic wasting disease to red deer (Census elaphus elaphus): early detection and late stage distribution of protease-resistant prion protein. Can Vet J 51:169–178

  66. 66.

    Di Bari MA, Nonno R, Castilla J, D’Agostino C, Pirisinu L, Riccardi G, Conte M, Richt J, Kunkle R, Langeveld J, Vaccari G, Agrimi U (2013) Chronic wasting disease in bank voles: characterisation of the shortest incubation time model for prion diseases. PLoS Pathog 9:e1003219

Download references

Acknowledgements

The authors would like to thank the British Deer Society and their members who were involved with collecting samples, especially John Bruce, who played an integral role in generating funding for the project, and in coordinating sample collection. We would also like to thank Helen Senn, Kirsty Macdonald of Highland game, Andy Foster of the Scottish Deer Centre, and Tobias Warnecke of Imperial College London for their help collecting and providing samples. We would like to thank Paula Stewart for her technical assistance. We are also grateful to Michael Tranulis for his critical review of the manuscript.

Funding

The work described in this paper was funded by Grants from Defra (Project SE2032), the British Deer Society and BBSRC (BBS/E/D/20002172) to ALR, WG and FH.

Author information

Study design and manuscript writing: ALR, JMP, WG, FH; sample collection and DNA preparation: ALR, HW, MEG, HT, KB, SLS, SP-E, JK-P; sample analysis: ALR, HW, MEG, HT, WG. All authors read and approved the final manuscript.

Correspondence to Amy L. Robinson.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark