Skip to main content

Prion gene haplotypes of U.S. cattle



Bovine spongiform encephalopathy (BSE) is a fatal neurological disorder characterized by abnormal deposits of a protease-resistant isoform of the prion protein. Characterizing linkage disequilibrium (LD) and haplotype networks within the bovine prion gene (PRNP) is important for 1) testing rare or common PRNP variation for an association with BSE and 2) interpreting any association of PRNP alleles with BSE susceptibility. The objective of this study was to identify polymorphisms and haplotypes within PRNP from the promoter region through the 3'UTR in a diverse sample of U.S. cattle genomes.


A 25.2-kb genomic region containing PRNP was sequenced from 192 diverse U.S. beef and dairy cattle. Sequence analyses identified 388 total polymorphisms, of which 287 have not previously been reported. The polymorphism alleles define PRNP by regions of high and low LD. High LD is present between alleles in the promoter region through exon 2 (6.7 kb). PRNP alleles within the majority of intron 2, the entire coding sequence and the untranslated region of exon 3 are in low LD (18.0 kb). Two haplotype networks, one representing the region of high LD and the other the region of low LD yielded nineteen different combinations that represent haplotypes spanning PRNP. The haplotype combinations are tagged by 19 polymorphisms (htSNPS) which characterize variation within and across PRNP.


The number of polymorphisms in the prion gene region of U.S. cattle is nearly four times greater than previously described. These polymorphisms define PRNP haplotypes that may influence BSE susceptibility in cattle.


Transmissible spongiform encephalopathies (TSEs) have been identified in humans, sheep, goats, deer, elk, moose, cattle, cats, and mink [1]. A cattle TSE, bovine spongiform encephalopathy (BSE), was first diagnosed among Holstein/Friesian cattle in the United Kingdom [2] and has since been detected in at least twenty five countries including the United States. The BSE agent is the probable cause of the human TSE, variant Creutzfeldt-Jakob Disease (vCJD) [3, 4], transmitted from cattle to people via the food chain.

Variation in the prion gene (PRNP) correlates with TSE progression in humans [5, 6], sheep [7], and mice [8]. In cattle, a 23-bp insertion/deletion (indel) polymorphism in the putative promoter region and a 12-bp indel within intron I have been associated with German BSE-affected animals [9]. These polymorphisms are present in U.S. cattle [10]. However, most of PRNP has not been characterized in a population as diverse as U.S. cattle outside of the coding region and 3'UTR of exon III, and portions of the promoter and intron I [1012]. Consequently, the extent of PRNP polymorphisms, linkage between PRNP alleles, recombination events, and haplotype diversity within PRNP is not known.

Public health concerns associated with vCJD and economic impacts of BSE on the cattle industry worldwide compel a thorough characterization of the genetic variation of bovine PRNP. Single nucleotide polymorphism (SNP) discovery in small populations introduce ascertainment biases of SNP properties [13, 14], and partial sequencing of genes in deep populations may characterize haplotype networks that extend past the sequenced region [15], yet still miss significant variation within the gene. The aim of this study was to characterize the extent of linkage disequilibrium (LD) and haplotype networks within PRNP ranging from the promoter past the 3'UTR (25.2 kb) in 192 U.S. cattle (16 beef and five dairy breeds). Reported here are 287 newly identified PRNP polymorphisms, the frequencies of 388 PRNP polymorphisms in U.S. beef and dairy cattle, a reference map of LD and haplotypes throughout PRNP, and the identification of 19 haplotype tagging SNPs (htSNPs) that are effective in U.S. populations of cattle. These results provide a reference framework for accurate and comprehensive evaluation of PRNP variation and its relationship to BSE.


Amplification and sequence coverage of the bovine PRNP gene from the promoter region through the 3'UTR (25.2 kb) in U.S. cattle

The PRNP gene was sequenced in 192 beef and dairy cattle; 16 beef and five dairy breeds from 24 overlapping amplicons (Figure 1 and see additional file 1). The amplification primers do not hybridize with PRNP regions containing any of the polymorphisms observed in this study, nor do any of the 150 sequencing primers used for redundant coverage of PRNP nucleotides (see additional file 1). Two or more high quality or unambiguous heterozygous reads were obtained for each PRNP nucleotide throughout 24.8 kb of the 25.2 kb region for approximately 95% of the cattle. Regions of PRNP that correspond with ambiguous sequence from more than five percent of the 192 animals were identified in the promoter region (95 bp), intron 1 (96 bp), and intron 2 (225 bp). These regions are attributable to closely positioned indels with heterozygous genotypes of high frequency in the cattle populations or stretches of mononucleotide repeats, both of which interfere with collection of high quality sequence. The positioning of these problematic loci was such that it was not possible to design amplification/sequencing reactions to cover the areas with high quality sequence.

Figure 1

Sequence coverage of PRNP in U.S. beef and dairy cattle. A. Physical map of PRNP that was sequenced in 192 beef and dairy cattle (25.2 kb) on 24 overlapping amplicons (blue bars). Orange and yellow color arrows represent untranslated and coding regions, respectively. B. Overall quality of PRNP sequence coverage. PRNP nucleotide sequence. with a phred score greater than 20 from at least two sequencing reads from the same animal was mapped to the corresponding nucleotide on a reference sequence [GenBank:AJ298878]. The blue line depicts the sequence quality of the 192 haploid cattle genomes comprising MBCDP2.1, the red line represents the 192 haploid cattle genomes comprising MDCP1.5, and the black line represents all 384 genomes combined. C. Close-up of PRNP sequence quality in an arbitrarily selected 2-kb window. Vertical arrows represent the location of polymorphisms (black for SNPs and white for indels). Polymorphisms with high minor allele frequencies are reflected by dips in sequence quality that correspond with heterozygous animals. The low phred score locus approximately at 48.5 kb is the result of animals heterozygous for two indels.

PRNP polymorphisms in U.S. cattle

A total of 388 polymorphisms (351 SNPs and 37 indels) were observed in PRNP gene sequences from the 384 chromosomes present in all 192 cattle (Figure 2 and see additional file 2). Two hundred and eighty-seven of the polymorphisms were not described in GenBank and literature searches as of July 10, 2006, and all were identified in non-coding regions of PRNP. The majority of polymorphisms (382/388) were observed in the multi-breed beef diversity panel (17 breeds, 192 chromosomes). In contrast, 158 polymorphisms were observed in the multi-breed dairy diversity panel (five breeds, 192 chromosomes), of which six were unique to the panel. Polymorphisms were observed in subgroups of cattle as follows (Figure 3): Bos taurus (240/388), British (161/388), Continental (216/388), Composites of U.S. Brahman (331/388), and Holstein (137/388). A total of 261 polymorphisms were used for haplotype inference across the five subgroups of cattle. One hundred and twenty polymorphisms were excluded from haplotype inference in all five subgroups of cattle due to low minor allele frequencies and an additional seven were excluded by Hardy-Weinberg testing. Sixty three of the 388 polymorphisms were only observed in one animal of the multi-breed beef diversity panel, a composite of U.S. Brahman.

Figure 2

PRNP polymorphism minor allele frequencies in U.S. beef (MBCDP2.1, n = 96) and dairy cattle (MDCP1.5, n = 96). Hatched bars represent the proportion of polymorphisms where the minor allele was observed in only one animal (singleton). Green and red colors represent the proportion of polymorphisms that have or have not been reported, respectively. Polymorphism numbers are show above bars. Multi-allelic polymorphisms are not represented.

Figure 3

PRNP polymorphism frequencies in U.S. cattle populations. PRNP polymorphism frequencies are shown for populations of Bos taurus (n = 94), Holstein (n = 86), Composite of U.S. Brahman (n = 20), British (n = 39), and Continental (n = 51). Hatched bars represent the proportion of polymorphisms where the minor allele was observed in only one animal (singleton). Polymorphism numbers are shown above the bars.

Linkage disequilibrium of PRNP alleles

A 6.7-kb region of high LD was identified in U.S. beef and dairy cattle from the 5' promoter region through exon 1, intron 1, exon 2, and part of intron 2 (Figure 4). Adjacent to the high LD region, an 18.0-kb portion of PRNP, containing the majority of intron 2 and all of exon 3, displayed markedly less LD (region of low LD). The region of high LD was not restricted to a cattle subgroup and was observed in the beef cattle diversity panel; dairy cattle diversity panel, B. taurus, British, Continental, and Holstein cattle subgroups (minor allele frequency ≥ 0.05, Hardy Weinberg p > 0.01). Including the polymorphisms with minor alleles of low frequency, 115 polymorphisms were identified in the 6.7 region of high LD, of which 45 have alleles in LD.

Figure 4

Pairwise plots of LD between PRNP alleles in U.S. cattle populations. Numbers refer to the physical location of polymorphisms in a PRNP consensus sequence file [GenBank:DQ457195] and the LD plots. Dark and light shading indicate high and low r2 values, respectively. Horizontal green bars highlight the region of high LD within PRNP that is flanked by polymorphisms at positions 449 and 7176.

Haplotype inference and networks

Nineteen PRNP haplotypes were inferred by the Expected Maximization (EM) algorithm on at least four chromosomes in one or more of the following subgroups: B. taurus, British, Continental, Holstein, and U.S. Brahman composite. Six of the haplotypes were observed in a homozygous state in one or more cattle subgroup. The 19 PRNP haplotypes represented 62% or more of the haploid genomes in all subgroups of cattle except U.S. Brahman composite (42%) and could be defined by the alleles of 19 htSNPs (Figure 5, Table 1). Because the haplotypes span two distinct PRNP regions defined by high and low LD, a Median-Joining network for each region was constructed (Figure 6). The network within the 6.7-kb region of high LD contains "sub-haplotypes" phased from nine of the 19 htSNPs and shows a linear stepwise relationship of alleles (Figure 6; network 1, Table 1). The network within the PRNP region of low LD contains sub-haplotypes phased from the remaining ten htSNPs and has a distinctly looped structure, indicating multiple unresolved allele relationships (Figure 6; network 2, Table 1). Sub-haplotype combinations from the two networks effectively account for the regions of high and low LD and yield haplotypes that span PRNP (Figure 6).

Figure 5

Haplotype tagging SNPs in five U.S. cattle populations. Numbers indicate the nucleotide position of htSNP alleles in a PRNP consensus sequence file [GenBank:DQ457195]. Black boxes indicate the cattle populations in which an SNP defines a haplotype.

Table 1 Bovine PRNP haplotype frequencies in populations of U.S. cattle
Figure 6

Bovine PRNP haplotypes defined by two Median-Joining networks. Two networks were constructed from 19 htSNPs. Network one was generated from htSNPs 449–6811 (n = 9), and network two from htSNPs 8631–21680 (n = 10). The physical positions of all 19 htSNPs are shown with solid vertical lines along the PRNP physical map. The dashed vertical line shows the position of a non-htSNP at nucleotide 7176 of the PRNP consensus sequence file [GenBank:DQ457195]. Circles represent extant PRNP haplotypes and are not scaled to frequency. Grey rectangles represent PRNP haplotypes not observed in our sample. Sequence for each network combination is presented in Table 1.


Sequencing the PRNP gene in 192 cattle representing 21 breeds resulted in the identification of 388 polymorphisms and detection of a region of high LD in the 5' non-coding region of PRNP. We identified 19 common PRNP haplotypes in U.S. cattle and characterized 19 htSNPs with the power to monitor these haplotypes. These results provide the means and a context for testing PRNP variation for an association with BSE.

Cattle present a challenge for LD analysis due to their complex history. Current haplotype patterns in domestic cattle are influenced by an evolutionary split of B. taurus and Bos indicus subspecies over 100,000 years ago [16], stratification of B. taurus germplasm within British and Continental populations, multiple population bottlenecks, gene migration, and inter-subspecies breeding. The cattle samples used in this study were classified into subgroups that addressed their natural history, an approach applied in other studies [15, 17]. Analyses of cattle subgroups by B. taurus, British, Continental, and Holstein classifications yielded similar LD regions within PRNP, indicating that the LD is consistent within and between B. taurus populations. A similar trend was observed within the small sample of B. indicus germplasm that was represented by the U.S. Brahman composite subgroup. The conservation of LD in PRNP indicates that some htSNPs should apply across populations.

Two regions of PRNP reflect different population histories (recombination, gene conversion, allelic fixation) by both pairwise r2 measurements of LD and haplotype Median-Joining Networks. Pairwise r2 values are lowered by historical recombination and the emergence of alleles on different haplotype lineages [17]. Consequently, alleles with high pairwise r2 values predict linear haplotype networks, and alleles with low pairwise r2 values predict complex or looped haplotype networks. The region of high LD that encompasses a 6.7-kb portion of PRNP extends to the 5' boundary of the PRNP locus sequenced in this study. Additional high LD may extend upstream of the sequenced locus. PRNP SNP alleles within this region are highly correlated with each other, and haplotypes phased within the region yield a linear network. In contrast, 18.0 kb of PRNP, which includes the entire protein coding and untranslated region of exon 3, has alleles in low LD and the associated haplotype network is complex. Alleles in this region may not be correlated with alleles elsewhere on the chromosome.

A previous study implicated a 23-bp indel in the promoter region and a 12-bp indel within intron 1 of PRNP with an association with susceptibility to BSE [9]. Both of these polymorphisms lie within the region of high LD described here, and their alleles are strongly correlated with the alleles of 43 other polymorphisms detected in this study. Until both chromosomal boundaries flanking the region of LD are determined, the number of polymorphisms with alleles in LD with those associated with BSE is unknown.

Although the diversity panels used to sequence PRNP represent a broad sample of U.S. cattle, it is likely that additional diversity within PRNP is present at low frequency in U.S. herds. This hypothesis is supported by the PRNP sequence from a single animal that accounted for 16.2 % of all observed polymorphisms. Some countries, including the U.S., have detected atypical BSE cases at exceedingly low frequencies with increased surveillance [1820]. Atypical BSE can differ from typical BSE by brain distribution and plaque morphology of the protease-resistant prion isoform (PrPres), or by western immunoblot profile of PrPres following proteinase K digestion [1820]. Diverse chromosomes in cattle populations could confound interpretations of PRNP variation identified from individual cases of either typical or atypical BSE.


The number of polymorphisms in the prion gene region of U.S. cattle is nearly four times greater than previously described. PRNP is divided into regions of high and low LD. The 19 htSNPs identified in this study define haplotype combinations from the two PRNP regions that may influence BSE susceptibility in cattle.


DNA panels used for PRNP DNA polymorphism discovery

The U.S. Meat Animal Research Center (USMARC) Beef Cattle Discovery Panel 2.1 (MBCDP2.1) was used to sample the breadth of PRNP genetic diversity in popular U.S. beef breeds. The selection and assembly of this panel has been previously described [21]. MBCDP2.1 consists of 96 bulls from the following breeds; Angus (n = 8), Hereford (n = 8), Limousin (n = 8), Simmental (n = 7), Charolais (n = 6), Beefmaster (n = 5), Red Angus (n = 6), Gelbvieh (n = 6), Brangus (n = 5), Salers (n = 5), Brahman (n = 6), Shorthorn (n = 5), Maine-Anjou (n = 5), Longhorn (n = 4), St. Gertrudis (n = 4), Chianina (n = 4), and Holstein (n = 4). The USMARC Dairy Cattle Panel (MDCP1.5) was used to sample the breadth of PRNP genetic diversity in dairy cattle and was a dairy subset of the USMARC-FSIS Random Market Cattle Panel version 1.1 described elsewhere [22]. MDCP1.5 consists of 96 dairy cows of the following breeds; Holstein (n = 82), Jersey (n = 7), Guernsey (n = 3), Aryshire (n = 2), and Brown Swiss (n = 2).

Primer design, PCR, and cycle sequencing

Reference sequence for the bovine PRNP gene was used as a template for primer design [GenBank:AJ298878], [23, 24]. Primers were designed to amplify 24 overlapping amplicons that collectively span a 25.2-kb region of PRNP (Oligo 6.61). Nested sequencing primers were designed for each amplicon to provide nucleotide coverage in both directions. Following preliminary experiments of amplification primer performance, the 192 cattle genomes comprising MBCDP2.1 and MDCP1.5 were subjected to 40 rounds of PCR with conditions as described [25] (see additional file 1). Following an Exonuclease I digestion [26], the PRNP amplicons were sequenced with BigDye terminator chemistry on an ABI 3730 capillary sequencer (PE Applied Biosystems, Foster City, CA).

Polymorphism detection, sequence quality, and coverage

Sequences from the 192 animals of the multi-breed beef and dairy panels were processed for polymorphism detection with Phred and Phrap [27, 28], Polyphred 3.5 [29], and Consed software [30]. A physical map linked to the PRNP consensus sequence and the location of polymorphisms was constructed in Vector NTI (v7.1). The map was annotated with all amplification and sequencing primers connected with the PRNP sequence. Replacement primers for those that hybridized to genomic loci containing polymorphisms were designed and used for additional amplification and sequencing of PRNP regions. PRNP nucleotide sequence with a phred score greater than 20 from at least two sequencing reads from the same animal was mapped to the corresponding nucleotide on reference sequence [GenBank:AJ298878]. Sequence compromised by SNP loci under associated amplification or sequencing primers was not analyzed for sequence coverage or the determination of genotypes. Regions reflecting poor sequence quality (<95% animal coverage) were identified, additional amplicons and sequencing primers were designed, and additional sequencing was performed. PRNP allele genotypes were mapped to reference sequence [GenBank:AJ298878] and stored in a relational database. A file of PRNP sequence annotated with all polymorphisms observed in this study and their frequencies in the beef diversity panel (MBCDP2.1) and dairy diversity panel (MDCP1.5) has been deposited in GenBank [GenBank:DQ457195].

Definitions of animal subgroups

The 192 beef and dairy animals whose genomic DNA comprise diversity panels MBCDP2.1 and MDCP1.5 were sorted into five subgroups based on breed composition. The B. taurus subgroup (n = 94) consisted of Angus (n = 8), Hereford (n = 8), Limousin (n = 8), Simmental (n = 7), Charolais (n = 6), Red Angus (n = 6), Gelbvieh (n = 6), Salers (n = 5), Shorthorn (n = 5), Maine-Anjou (n = 5), Texas Longhorn (n = 4), Chianina (n = 4), Aryshire (n = 2), Brown Swiss (n = 2), Guernsey (n = 3), Jersey (n = 7), and Holstein (n = 8; randomly selected from 86 to avoid over-representation of the Holstein breed). The British subgroup (n = 39) consisted of Angus (n = 8), Hereford (n = 8), Red Angus (n = 6), Shorthorn (n = 5), Aryshire (n = 2), Jersey (n = 7), and Guernsey (n = 3). The Continental subgroup (n = 51) consisted of Charolais (n = 6), Chianina (n = 4), Gelbvieh (n = 6), Limousin (n = 8), Maine-Anjou (n = 5), Salers (n = 5), Simmental (n = 7), Brown Swiss (n = 2), and Holstein (n = 8; randomly selected from 86). The Holstein subgroup consisted of 86 Holsteins. The final subgroup, Composites of U.S. Brahman (n = 20) consisted of Beefmaster (n = 5), Brahman (n = 6), Brangus (n = 5), and Santa Gertrudis (n = 4).

LD estimation, haplotype inference, and median-joining network analyses

Unphased PRNP genotypes were assembled for each animal in datasets of the B. taurus, British, Continental, Holstein, and Composite of U.S. Brahman subgroups. Polymorphisms with more than two alleles, a minor allele frequency <0.05, or those not in Hardy-Weinberg equilibrium (Chi-square p < 0.01) were excluded from further analyses. Our cattle populations were not the result of random mating, violating an assumption of Hardy-Weinberg equilibrium. However, the Hardy-Weinberg test facilitated the identification of common haplotypes within the subgroups by excluding polymorphisms where the minor allele was amplified in a particular breed, yet had a low overall frequency.

The extent of LD between the PRNP alleles of each dataset was calculated with pairwise r2 values (Haploview v3.2 [31]). Regions of LD were determined through visual inspection of LD graphs. Haplotypes were inferred in Haploview using the EM algorithm and a minimal set of polymorphisms was identified that collectively tagged all observed haplotypes predicted on four or more chromosomes in one or more of the five subgroups. The htSNPs identified across the five subgroup datasets were combined into a single set of 19 htSNPS. The 19 htSNPs were used to infer haplotypes within the five subgroup datasets. Median-Joining networks of PRNP haplotypes were constructed in Network (v4.111)[32].


  1. 1.

    Aguzzi A, Sigurdson CJ: Antiprion immunotherapy: to suppress or to stimulate?. Nat Rev Immunol. 2004, 4 (9): 725-736. 10.1038/nri1437.

    PubMed  CAS  Article  Google Scholar 

  2. 2.

    Wells GA, Scott AC, Johnson CT, Gunning RF, Hancock RD, Jeffrey M, Dawson M, Bradley R: A novel progressive spongiform encephalopathy in cattle. Vet Rec. 1987, 121: 419-420.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Bruce ME, Will RG, Ironside JW, McConnell I, Drummond D, Suttie A, McCardle L, Chree A, Hope J, Birkett C, Cousens S, Fraser H, Bostock CJ: Transmissions to mice indicate that 'new variant' CJD is caused by the BSE agent. Nature. 1997, 389 (6650): 498-501. 10.1038/39057.

    PubMed  CAS  Article  Google Scholar 

  4. 4.

    Hill AF, Desbruslais M, Joiner S, Sidle KC, Gowland I, Collinge J, Doey LJ, Lantos P: The same prion strain causes vCJD and BSE. Nature. 1997, 389 (6650): 448-450. 10.1038/38925.

    PubMed  CAS  Article  Google Scholar 

  5. 5.

    Wadsworth JD, Asante EA, Desbruslais M, Linehan JM, Joiner S, Gowland I, Welch J, Stone L, Lloyd SE, Hill AF, Brandner S, Collinge J: Human prion protein with valine 129 prevents expression of variant CJD phenotype. Science. 2004, 306 (5702): 1793-1796. 10.1126/science.1103932.

    PubMed  CAS  Article  Google Scholar 

  6. 6.

    Prusiner SB, Scott MR: Genetics of prions. Annu Rev Genet. 1997, 31: 139-175. 10.1146/annurev.genet.31.1.139.

    PubMed  CAS  Article  Google Scholar 

  7. 7.

    Baylis M, Goldmann W: The genetics of scrapie in sheep and goats. Curr Mol Med. 2004, 4 (4): 385-396. 10.2174/1566524043360672.

    PubMed  CAS  Article  Google Scholar 

  8. 8.

    Barron RM, Baybutt H, Tuzi NL, McCormack J, King D, Moore RC, Melton DW, Manson JC: Polymorphisms at codons 108 and 189 in murine PrP play distinct roles in the control of scrapie incubation time. J Gen Virol. 2005, 86 (Pt 3): 859-868. 10.1099/vir.0.80525-0.

    PubMed  CAS  Article  Google Scholar 

  9. 9.

    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 protein gene (PRNP) in German cattle breeds. Neurogenetics. 2004, 5: 19-25. 10.1007/s10048-003-0171-y.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Seabury CM, Womack JE, Piedrahita J, Derr JN: Comparative PRNPgenotyping of U.S. cattle sires for potential association with BSE. Mamm Genome. 2004, 15: 828-833. 10.1007/s00335-004-2400-6.

    PubMed  CAS  Article  Google Scholar 

  11. 11.

    Heaton MP, Leymaster KA, Freking BA, Hawk DA, Smith TP, Keele JW, Snelling WM, Fox JM, Chitko-McKown CG, Laegreid WW: Prion gene sequence variation within diverse groups of U.S. sheep, beef cattle, and deer. Mamm Genome. 2003, 14 (11): 765-777. 10.1007/s00335-003-2283-y.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Seabury CM, Honeycutt RL, Rooney AP, Halbert ND, Derr JN: Prion protein gene (PRNP) variants and evidence for strong purifying selection in functionally important regions of bovine exon 3. Proc Natl Acad Sci USA. 2004, 101 (42): 15142-15147. 10.1073/pnas.0406403101.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  13. 13.

    Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R: Ascertainmant bias in studies of huamn genome-wide polymorphism. Genome Res. 2005, 15: 1496-1502. 10.1101/gr.4107905.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  14. 14.

    Soldevila M, Andres AM, Ramirez-Soriano A, Marques-Bonet T, Calafell F, Navarro A, Bertranpetit J: The prion protein gene in humans revisited: Lessons from a worldwide resequencing study. Genome Res. 2006, 16: 231-239. 10.1101/gr.4345506.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  15. 15.

    Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.

    PubMed  CAS  Article  Google Scholar 

  16. 16.

    Troy CS, MacHugh DE, Bailey JF, Magee DA, Loftus RT, Cunningham P, Chamberlain AT, Sykes BC, Bradley DG: Genetic evidence for Near--Eastern origins of European cattle. Nature. 2001, 410: 1088-1091. 10.1038/35074088.

    PubMed  CAS  Article  Google Scholar 

  17. 17.

    Consortium HM: A haplotype map of the human genome. Nature. 2005, 437: 1299-1320. 10.1038/nature04226.

    Article  Google Scholar 

  18. 18.

    Biacabe AG, Laplanche JL, Ryder S, Baron T: Distinct molecular phenotypes in bovine prion diseases. EMBO Rep. 2004, 5 (1): 110-115. 10.1038/sj.embor.7400054.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  19. 19.

    Casalone C, Zanusso G, Acutis P, Ferrari S, Capucci L, Tagliavini F, Monaco S, Caramelli M: Identification of a second bovine amyloidotic spongiform encephalopathy: molecular similarities with sporadic Creutzfeldt-Jakob disease. Proc Natl Acad Sci U S A. 2004, 101 (9): 3065-3070. 10.1073/pnas.0305777101.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  20. 20.

    Yamakawa Y, Hagiwara K, Nohtomi K, Nakamura Y, Nishijima M, Higuchi Y, Sato Y, Sata T: Atypical proteinase K-resistant prion protein (PrPres) observed in an apparently healthy 23-month-old Holstein steer. Jpn J Infect Dis. 2003, 56 (5-6): 221-222.

    PubMed  Google Scholar 

  21. 21.

    Heaton MP, Chitko-McKown CG, Grosse WM, Keele JW, Keen JE: Interleukin-8 haplotype structure from nucleotide sequence variation in commercial populations of U.S. beef cattle. Mamm Genome. 2001, 12: 219-226. 10.1007/s003350010269.

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Heaton MP, Keen JE, Clawson ML, Harhay GP, Bauer N, Shultz C, Green BT, Durso L, Chitko-McKown CG, Laegreid WW: Use of bovine single nucleotide polymorphism markers to verify sample tracking in beef processing. J Am Vet Med Assoc. 2005, 226 (8): 1311-1314. 10.2460/javma.2005.226.1311.

    PubMed  CAS  Article  Google Scholar 

  23. 23.

    Hills D, Schlaepfer J, Comincini S, MacLean I, Dolf G, Ferretti L, Olsaker I, Williams JL: Sequence variation in the bovine and ovine PRNP genes. Animal Genetics. 2003, 34: 183-190. 10.1046/j.1365-2052.2003.00977.x.

    PubMed  CAS  Article  Google Scholar 

  24. 24.

    Hills D, Comincini S, Schlaepfer J, Dolf G, Ferretti L, Williams JL: Complete genomic sequence of the bovine prion gene (PRNP) and polymorphisms in its promoter region. Animal Genetics. 2001, 32: 231-233. 10.1046/j.1365-2052.2001.0769a.x.

    PubMed  CAS  Article  Google Scholar 

  25. 25.

    Clawson ML, Heaton MP, Chitko-McKown CG, Fox JM, Smith TPL, Snelling WM, Keele JW, Laegreid WW: Beta-2-microglobulin haplotypes in U.S. beef cattle and association with failure of passive transfer in newborn calves. Mamm Genome. 2004, 15: 227-236. 10.1007/s00335-003-2320-x.

    PubMed  CAS  Article  Google Scholar 

  26. 26.

    Smith TP, Godtel RA, Lee RT: PCR-based setup for high-throughput cDNA library sequencing on the ABI 3700 automated DNA sequencer. Biotechniques. 2000, 29 (4): 698-700.

    PubMed  CAS  Google Scholar 

  27. 27.

    Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.

    PubMed  CAS  Article  Google Scholar 

  28. 28.

    Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.

    PubMed  CAS  Article  Google Scholar 

  29. 29.

    Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25 (14): 2745-2751. 10.1093/nar/25.14.2745.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  30. 30.

    Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8 (3): 195-202.

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.

    PubMed  CAS  Article  Google Scholar 

  32. 32.

    Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16 (1): 37-48.

    PubMed  CAS  Article  Google Scholar 

Download references


We thank Gennie Schuller-Chavez, Tammy Sorensen, and Jacky Carnahan for sequence analysis; Bucky Herman, Phil Anderson, Jim Wray, and Randy Bradley for database support; and Joan Rosch for manuscript preparation. This research was supported by the USDA National Research Initiative, Competitive Grant No. 2005-35212-15890 (W.W.L., M.L.C., M.P.H.) and the Agricultural Research Service.

The use of product and company names is necessary to accurately report the methods and results; however, the USDA neither guarantees nor warrants the standard of the products, and the use of names by the USDA implies no approval of the product to the exclusion of others that may also be suitable.

Author information



Corresponding author

Correspondence to Michael L Clawson.

Additional information

Authors' contributions

MLC participated in the project conception and design, data generation and management, sequence analyses, polymorphism genotyping, LD and haplotype analyses, and drafted the manuscript. MPH participated in the project conception and design, data generation, sequence analyses, script development regarding sequence coverage and polymorphism genotype mapping, and suggested improvements regarding the data pipeline. JWK participated in sequence analyses and wrote scripts for a semi-automated pipeline that stored DNA trace files and polymorphism genotypes within a relational database and mapped sequence phred scores and polymorphism genotypes onto a reference sequence. TPLS performed PRNP sequencing, analyzed sequence, and suggested improvements regarding the data pipeline. GPH participated in LD and haplotype analyses. WWL participated in the project conception and design, data generation, LD and haplotype analyses, and suggested improvements regarding the data pipeline.

Electronic supplementary material

Oligonucleotides for

Additional file 1: PRNP amplification and sequencing. Annealing location, USMARC number, DNA sequence, orientation, function, and annealing temperature for 198 oligonucleotides used for PRNP amplification and sequencing. (XLS 1 MB)

Allele and genotype frequencies of

Additional file 2: PRNP polymorphisms in the USMARC beef and dairy cattle panels. Physical positions and PRNP regions of all 388 polymorphisms observed in this study; their allele and genotype frequencies in the beef and dairy cattle panels, and the GenBank accession numbers for 101 polymorphisms previously reported. (XLS 158 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Clawson, M.L., Heaton, M.P., Keele, J.W. et al. Prion gene haplotypes of U.S. cattle. BMC Genet 7, 51 (2006).

Download citation


  • Linkage Disequilibrium
  • Bovine Spongiform Encephalopathy
  • Haplotype Network
  • High Linkage Disequilibrium
  • Diversity Panel