Association analysis of PRNP gene region with chronic wasting disease in Rocky Mountain elk
BMC Research Notes volume 3, Article number: 314 (2010)
Chronic wasting disease (CWD) is a transmissible spongiform encephalopathy (TSE) of cervids including white-tailed (Odocoileus virginianus) and mule deer (Odocoileus hemionus), Rocky Mountain elk (Cervus elaphus nelsoni), and moose (Alces alces). A leucine variant at position 132 (132L) in prion protein of Rocky Mountain elk confers a long incubation time with CWD, but not complete resistance. However, variants in regulatory regions outside the open reading frame of PRNP have been associated with varying degrees of susceptibility to prion disease in other species, and some variants have been observed in similar regions of Rocky Mountain elk PRNP. Thus, additional genetic variants might provide increased protection, either alone or in combination with 132L.
This study provided genomic sequence of all exons for PRNP of Rocky Mountain elk. Many functional sites in and around the PRNP gene region were sequenced, and this report approximately doubled (to 75) the number of known variants in this region. A haplotype-tagging approach was used to reduce the number of genetic variants required to survey this variation in the PRNP gene region of 559 Rocky Mountain elk. Eight haplotypes were observed with frequencies over 1.0%, and one haplotype was present at 71.2% frequency, reflecting limited genetic diversity in the PRNP gene region.
The presence of 132L cut odds of CWD by more than half (Odds Ratio = 0.43; P = 0.0031), which was similar to a previous report. However after accounting for 132L, no association with CWD was found for any additional variants in the PRNP region (P > 0.05).
Chronic wasting disease (CWD) is a transmissible spongiform encephalopathy (TSE) of cervids including white-tailed and mule deer (Odocoileus virginianus and O. hemionus, respectively), Rocky Mountain elk (Cervus elaphus nelsoni), and moose (Alces alces) . This invariably fatal disease can be associated with progressive physical wasting, often accompanied by behavioral changes including listlessness, decreased interaction with other animals, decreased responsiveness to environmental stimuli, and repetitive behaviors . There is no cure for CWD or any of the TSEs, and current management practice to reduce the spread of CWD includes depopulation or permanent quarantine of infected farmed herds as well as removal of free-ranging animals showing clinical signs of disease. Despite these efforts, prevalence continues to climb and the range of CWD continues to expand [2, 3].
The TSEs are believed to be caused by a misfolded prion protein (PrPd) that is protease resistant and infectious . The prion protein is highly conserved and the normal cellular form of the protein (PrPc) is expressed in neurons of all mammals studied to date . The disease form PrPd is capable of recruiting the normal cellular form of the protein (PrPc) into the same misfolded conformation . While there is no cure for CWD, genetic variants in and around the prion gene confer varying degrees of resistance to TSEs in a range of species. In sheep, an arginine substitution at position 171 (171R) confers strong resistance to classical scrapie, the TSE of sheep and goats . Scrapie eradication programs in the U.S. and Europe include a large component of breeding animals toward resistant genotypes, and classical scrapie prevalence is in decline [7, 8].
Genetic solutions may be possible for Rocky Mountain elk, and a leucine variant at position 132 (132L) in prion protein confers a long incubation time with CWD, but not complete resistance [9–11]. However, additional genetic variants might provide increased protection, either alone or in combination with 132L. Genetic variants in regulatory regions outside the open reading frame of the prion gene have been associated with susceptibility to or incubation time with prion disease in cattle , and some genetic variants have been observed in similar regions of Rocky Mountain elk [9, 13, 14].
The many variants in any gene region are arranged in only a limited number of haplotypes, or linear arrangements of genetic variants as they exist on chromosome segments in a population. To find an important genetic variant, one need only genotype a limited set of variants that "tag" all the major haplotypes in a gene region to capture all or most of the information contained in the complete set of genetic variants in a region . This strategy is known as haplotype-tagging, and it often reduces the number of variants that need to be genotyped by 50-75% or more, greatly reducing the time, cost, and complexity of analysis required to test a region for important genetic variants [16, 17]. Importantly, because any disease-related variant will occur on one or more haplotypes in the population, haplotype tagging allows the testing of haplotypes for association with disease whether or not the underlying causative variants are in the dataset, or even known at all. This has given rise to a two-stage strategy for assessing genetic association: first, a relatively small number of individuals are genotyped to determine common haplotypes in a population, and then all individuals are genotyped with only a set of haplotype-tagging markers identified in the first round . Therefore, this study used a two-stage haplotype-tagging strategy to search for any variants in the PRNP gene region - whether previously known, discovered in this study, or as yet undiscovered - that may provide increased protection from CWD in Rocky Mountain elk.
A total of 559 captive and free-ranging Rocky Mountain elk were sampled from herds following positive CWD diagnosis in 6 states including Colorado, Montana, Minnesota, Nebraska, Oklahoma, and South Dakota. All animal procedures used were exempt by the Institutional Animal Care and Use Committee of Washington State University as no live animal use was involved; all samples were from depopulation programs approved by federal and state regulatory bodies to control the spread of CWD. In total, 120 animals tested positive and 439 tested negative for CWD by immunohistochemical methods as previously described . Briefly, formalin fixed tissues including brain, retropharyngeal lymph node, and tonsil were processed for PrPd detection using an automated monoclonal antibody immunohistochemistry assay . The standard for considering an animal CWD positive was at least one tissue positive for PrPd. Because the time between infection and appearance of detectable PrPd is unknown , elk from herds exposed to CWD but individually lacking detectable PrPd were not defined as CWD negative; these animals were defined as herd-matched controls.
Sequence and Association Analysis
Genetic variant discovery was performed on a group of 20 animals chosen for geographic diversity, including animals from Colorado, Montana, Minnesota, and Oklahoma. PCR was performed using primers as shown (additional file 1: amplification primers and conditions) and the following standard PCR conditions: an initial denaturation at 95 degrees for 5 minutes, followed by 35 cycles of 95 degrees for 30 seconds, the annealing temperature (see additional file 1) for 30 seconds, and 68 degrees for 90 seconds, and a final extension step at 68 degrees for 10 minutes. Sequencing employed BigDye chemistry (Applied Biosystems, Foster City, CA) and primers as shown (additional file 2: sequencing primers). The regions sequenced included the promoter, 5' UTR, entire open reading frame, 3' UTR, all exon splice sites, and flanking regions located approximately 20 kb upstream of the sequenced promoter region and approximately 20 kb downstream of exon 3. Deeper resequencing using a subset of these primers was performed on a group of 93 animals (46 cases, 47 controls) to identify additional variants including lower frequency alleles that may be associated with CWD.
Haplotype tagging markers (Table 1) were chosen using an r2 threshold of 0.8 in Tagger as implemented in HAPLOVIEW  to represent every variant present at greater than 5% minor allele frequency in either the geographic diversity group of 20 or in the case-control group. Genotyping assays were designed for these haplotype-tagging markers (Tables 2, 3), and genotyping was performed on the full animal set as previously described . Briefly, fluorescent TaqMan genotyping assays were performed according to manufacturer specifications, including standard cycling conditions, (Applied Biosystems, Foster City, CA) using primer and probe sets as shown in Table 2. Restriction fragment length polymorphism (RFLP) assays were performed using standard PCR conditions (see above) and restriction enzymes listed in Table 3 according to manufacturer specifications (New England Biolabs, Ipswich, MA; Fermentas Inc., Glen Burnie, MD; Roche Diagnostics, Indianapolis, IN), and visualized on 2% (mass/volume) agarose gels. One haplotype-tagging marker (g.152T > C; [dbSNP:ss115456962]) was not suitable for either fluorescent or RFLP genotyping assays, so allele-specific PCR was performed to identify genotypes for this marker. The reaction to detect the C allele of this marker used standard PCR conditions (above) with a 51 degree annealing temperature, and employed the following primers: CTAGGTGGAATCAGTCGYAC and GGACTTTGCCCAGAGGGTAG. The reaction to detect the T allele was similar except for a 61 degree annealing temperature, and it employed primers CTAGGTGGAATCAGTCTYAT and GGACTTTGCCCAGAGGGTAG. As with RFLP assays, the results were visualized on 2% (mass/volume) agarose gels. The threshold for successful marker genotyping was genotypes on at least 95% of animals. Haplotype frequencies were calculated using PHASE v2 (Table 4) [20, 21]. The threshold for successful genotypic markers was genotype identification in at least 95% of animals with at least 95% confidence in PHASE v2.
Baseline association of the 132L variant was performed using the logistic procedure of SAS 9.2 (SAS Institute, Cary, NC). Association analysis of other markers while accounting for 132L was performed on markers with greater than 5% minor allele frequency using the logistic procedure of SAS 9.2 with a model including the presence/absence of 132L in addition to presence/absence of the variant of interest. Additionally, similar logistic models were used to test full genotypes while accounting for 132L by performing exact conditional tests (to account for low minor allele frequencies) for genotypes of 132L and the variant of interest. Fisher's exact was used to perform specific comparison of CWD frequency among L132-bearing haplotypes.
Results and Discussion
While it is clearly possible to have strong genetic resistance to TSEs, Rocky Mountain elk have only been shown to have the 132L substitution which confers extended incubation time with CWD but not complete resistance [9, 10]. Additional variants in regulatory sites outside the open reading frame could provide increased protection, and this study used a haplotype-tagging approach to identify such variants. Specifically, this study expanded the available genomic DNA sequence and variants in the PRNP gene region, assessed haplotype structure in this gene region in Rocky Mountain elk, and tested common haplotype variants for association with CWD.
The current data study greatly expanded both the genomic sequence and known variants in the PRNP gene region of Rocky Mountain elk. As depicted in Figure 1, sequence data were generated on the proximal promoter region, all exons and splice sites, the 3'UTR, and flanking regions of the PRNP gene of Rocky Mountain elk. The total sequence was over 7 kb from a region spanning approximately 63 kb, of which over 40% was not previously reported to the best of our knowledge (Table 5). Further, the markers provided here approximately doubled the number of publically available variants in the PRNP gene region to 75 (Table 5). Most of these additional variants were obtained from the 20 elk chosen for geographic diversity. There was only one additional variant (dbSNP:ss119994812) observed from the deeper resequencing of additional case-control animals, and it had a low minor allele frequency (0.5%) in the case-control group.
As anticipated, these 75 total variants were found to be organized into a much smaller number of haplotypes. Out of 559 elk, only 19 haplotypes were observed 3 or more times in the sample set, and only 8 haplotypes were observed at 1.0% or greater allele frequency in the total sample, with only small variations once subdivided by captive or free-ranging state of the animals (Table 4). However, since the free-ranging animals were obtained from a smaller number of locations from only one state and since captive animals have had documented human-assisted gene flow, it is difficult to make any inferences about the genetic diversity of captive versus free-ranging animals from these data. Haplotype tagging enabled representation of all common underlying variants at r2 of 0.80 or greater using only 17 SNP markers, plus 5 derived genotypes composed of short multimarker haplotypes (Table 1). This represents a 77.3% reduction in variants that needed to be genotyped, which is comparable to other reports for haplotype tagging strategies in mammals [16, 17]. Further, the haplotype tagging procedure required coverage of every variant in each of two animal groups, even though the variants observed in each group varied somewhat. This approach was a conservative way to ensure that all markers received coverage, but at the expense of possibly overestimating the number of markers required. The small number of haplotypes reflects relatively limited diversity in the PRNP gene region among Rocky Mountain elk , which is consistent with other reports of low genetic diversity using microsatellites in Rocky Mountain elk . However, the paucity of studies based on comparable marker types precludes conclusions regarding selective sweeps in the PRNP region as compared to the rest of the genome. Overall, this group of tagging markers provided representation of the genetic diversity in the PRNP gene region for measuring local linkage disequilibrium between gene regions and for association testing with CWD.
These markers were used to investigate whether additional markers could improve the resistance provided by the previously described 132L variant (O'Rourke et al 1999). The 132L variant was underrepresented among CWD cases (P = 0.0031), occurring in cases less than half as often as the predominant genotype 132 MM (OR = 0.43, 95%CI: 0.25-0.75). However, after accounting for 132L no other variants showed significant association with CWD even on a nominal basis (P > 0.05), before any correction for multiple testing. The statistical tests specifically included comparison of CWD frequency among carriers of two haplotypes (2 and 7) that harbor L132, but no significant differences were observed (P = 0.99). Furthermore, all genotypes with any appreciable frequency showed CWD positive animals, suggesting that there is no complete resistance to CWD on the basis of common PRNP genotypes in these elk. While we are not aware of any epidemiological evidence to suggest that truly CWD resistant elk exist, future research could examine the possibility that complete genetic resistance to CWD does exist in Rocky Mountain elk, either because of a very low frequency PRNP allele or because of the influence of other genes.
Williams ES: Chronic wasting disease. Vet Pathol. 2005, 42: 530-549. 10.1354/vp.42-5-530.
Miller MW, Williams ES, McCarty CW, Spraker TR, Kreeger TJ, Larsen CT, Thorne ET: Epizootiology of chronic wasting disease in free-ranging cervids in Colorado and Wyoming. J Wildl Dis. 2000, 36: 676-690.
Argue CK, Ribble C, Lees VW, McLane J, Balachandran A: Epidemiology of an outbreak of chronic wasting disease on elk farms in Saskatchewan. Can Vet J. 2007, 48: 1241-1248.
Prusiner SB: Novel proteinaceous infectious particles cause scrapie. Science. 1982, 216: 136-144. 10.1126/science.6801762.
Roucou X, Gains M, LeBlanc AC: Neuroprotective functions of prion protein. J Neurosci Res. 2004, 75: 153-161. 10.1002/jnr.10864.
Hunter N: Scrapie: uncertainties, biology and molecular approaches. Biochim Biophys Acta. 2007, 1772: 619-628.
Melchior MB, Windig JJ, Hagenaars TJ, Bossers A, Davidse A, van Zijderveld FG: Eradication of scrapie with selective breeding: are we nearly there?. BMC Vet Res. 2010, 6: 24-10.1186/1746-6148-6-24.
Hagenaars TJ, Melchior MB, Bossers A, Davidse A, Engel B, van Zijderveld FG: Scrapie prevalence in sheep of susceptible genotype is declining in a population subject to breeding for resistance. BMC Vet Res. 2010, 6: 25-10.1186/1746-6148-6-25.
O'Rourke KI, Besser TE, Miller MW, Cline TF, Spraker TR, Jenny AL, Wild MA, Zebarth GL, Williams ES: PrP genotypes of captive and free-ranging Rocky Mountain elk (Cervus elaphus nelsoni) with chronic wasting disease. J Gen Virol. 1999, 80: 2765-2769.
Spraker TR, Balachandran A, Zhuang D, O'Rourke KI: Variable patterns of distribution of PrP(CWD) in the obex and cranial lymphoid tissues of Rocky Mountain elk (Cervus elaphus nelsoni) with subclinical chronic wasting disease. Vet Rec. 2004, 155: 295-302. 10.1136/vr.155.10.295.
Green KM, Browning SR, Seward TS, Jewell JE, Ross DL, Green MA, Williams ES, Hoover EA, Telling GC: The elk PRNP codon 132 polymorphism controls cervid and scrapie prion propagation. J Gen Virol. 2008, 89: 598-608. 10.1099/vir.0.83168-0.
Sander P, Hamann H, Pfeiffer I, Wemheuer W, Brenig B, Groschup MH, Ziegler U, Distl O, Leeb T: Analysis of sequence variability of the bovine prion protein gene (PRNP) in German cattle breeds. Neurogenetics. 2004, 5: 19-25. 10.1007/s10048-003-0171-y.
Seabury CM, Gill CA, Templeton JW, Dyar JB, Derr JN, Adelson DL, Owens E, Davis DS, Kraemer DC, Womack JE: Molecular characterization of the Rocky Mountain elk (Cervus elaphus nelsoni) PRNP putative promoter. J Hered. 2007, 98: 678-686. 10.1093/jhered/esm091.
Perucchini M, Griffin K, Miller MW, Goldmann W: PrP genotypes of free-ranging wapiti (Cervus elaphus nelsoni) with chronic wasting disease. J Gen Virol. 2008, 89: 1324-1328. 10.1099/vir.0.83424-0.
Johnson GC, Esposito L, Barratt BJ, Smith AN, Heward J, Di Genova G, Ueda H, Cordell HJ, Eaves IA, Dudbridge F, Twells RC, Payne F, Hughes W, Nutland S, Stevens H, Carr P, Tuomilehto-Wolf E, Tuomilehto J, Gough SC, Clayton DG, Todd JA: Haplotype tagging for the identification of common disease genes. Nat Genet. 2001, 29: 233-237. 10.1038/ng1001-233.
Thompson D, Stram D, Goldgar D, Witte JS: Haplotype tagging single nucleotide polymorphisms and association studies. Hum Hered. 2003, 56: 48-55. 10.1159/000073732.
Lowe CE, Cooper JD, Chapman JM, Barratt BJ, Twells RC, Green EA, Savage DA, Guja C, Ionescu-Tirgoviste C, Tuomilehto-Wolf E, Tuomilehto J, Todd JA, Clayton DG: Cost-effective analysis of candidate genes using htSNPs: a staged approach. Genes Immun. 2004, 5: 301-305. 10.1038/sj.gene.6364064.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
White SN, Mousel MR, Reynolds JO, Lewis GS, Herrmann-Hoesing LM: Common promoter deletion is associated with 3.9-fold differential transcription of ovine CCR5 and reduced proviral level of ovine progressive pneumonia virus. Anim Genet. 2009, 40: 583-589. 10.1111/j.1365-2052.2009.01882.x.
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-989. 10.1086/319501.
Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73: 1162-1169. 10.1086/379378.
Polziehn RO, Hamr J, Mallory FF, Strobeck C: Microsatellite analysis of North American wapiti (Cervus elaphus) populations. Mol Ecol. 2000, 9: 1561-1576. 10.1046/j.1365-294x.2000.01033.x.
We thank Codie Durfee for technical assistance. We also gratefully acknowledge veterinarians and technicians with the USDA Animal Plant Health Inspection Service and the state of Colorado for epidemiology and tissue collection. This work was funded by USDA-ARS grant 5348-32000-026-00 D.
Katherine I. O'Rourke is the inventor of patented monoclonal antibody 99.
SNW conceived the study, carried out the statistical analyses, and drafted the manuscript. TRS participated in sample collection and performed diagnostic pathology. JOR participated in study design, genotyping, and analysis. KIO participated in study design, sample collection, genotyping, and analysis. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Amplification primers and conditions. Primer sequences, strand directions, base positions relative to reference sequences, annealing temperatures, and included region descriptions. (XLS 47 KB)
Additional file 2: Sequencing primers. Primer sequences, strand directions, base positions relative to reference sequences, and descriptions of regions sequenced using each primer. (XLS 16 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
White, S.N., Spraker, T.R., Reynolds, J.O. et al. Association analysis of PRNP gene region with chronic wasting disease in Rocky Mountain elk. BMC Res Notes 3, 314 (2010). https://0-doi-org.brum.beds.ac.uk/10.1186/1756-0500-3-314