|
|
||||||||
ANIMAL GENETICS |
US Meat Animal Research Center, USDA, ARS, Clay Center, NE 68933-0166
| Abstract |
|---|
|
|
|---|
Key Words: carcass trait growth trait meat quality trait prion gene scrapie resistance sheep
| INTRODUCTION |
|---|
|
|
|---|
A genetic component to scrapie susceptibility has been known for many years. The 3 most important genetic markers associated with scrapie susceptibility are located at codons 136, 154, and 171 of the prion (PRNP) gene (for review see Baylis and Goldmann, 2004
). The PRNP haplotype encoding alanine, arginine, and arginine (ARR) at the respective 136, 154, and 171 positions is associated with increased resistance to scrapie, whereas the valine, arginine, and glutamine haplotype (VRQ) is associated with increased susceptibility to scrapie. The scrapie susceptibility of 3 other common PRNP haplotypes (ARQ, AHQ, and ARH) is intermediate or unknown.
Many countries currently have scrapie eradication plans utilizing these genetic markers (Byrne, 2003
). Several researchers have investigated possible antagonistic associations between the PRNP locus and milk (de Vries et al., 2005
), performance (Alexander et al., 2005
), and growth (Brandsma et al., 2004
) traits. However, consistent evidence supporting antagonistic associations has not been found. The objective of this study was to test for associations between PRNP haplotypes and a comprehensive set of growth, carcass, and meat quality traits in an F2 population of Dorset x Romanov sheep.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Carcasses were split along the dorsal midline and the length of the carcass was measured from the anterior edge of the first rib to the anterior edge of the aitchbone. Metacarpal bone length was measured as an additional indicator of skeletal development. Fat depths were measured at the midpoint of the longissimus muscle between the 12th and 13th ribs and at the midline of the fourth sacral vertebra. The left side of the carcass was retained for later use in shear force determination. The right side of the carcass was cut between the fifth and sixth ribs and between the last 2 lumbar vertebrae to produce 3 carcass subsections. Subsections from the right side corresponded to anterior (shoulder, neck, foreshank, and breast), middle (loin, rib, and flank), and posterior (sirloin, leg, and hindshank) portions of the carcass. Each subsection was individually weighed. Longissimus muscle depth, width, surface area, and subjective marbling score (on a scale from 0 = devoid to 600 = moderate) were measured between the 12th and 13th ribs.
Proximate chemical analysis was performed on each of the carcass subsections. Each subsection was frozen and ground separately to measure ash, water, protein, and ether extract (fat) content (AOAC, 1990
).
From the left side of the carcass, loin chops (2.5 cm thick) were collected for use in shear force determination. Loin chops were aged at 4°C for 14 d and stored at 20°C until shear force measurement. Storage duration ranged from 2 wk to 5 mo. Before shear force measurement, chops were thawed at 5°C, broiled to an internal temperature of 40°C, and then turned and broiled to an internal temperature of 75°C using an open-hearth electric broiler (Farberware, Bronx, NY). Loin chops were allowed to cool at 4°C for 24 h before 6 cores (1.27 cm diameter each) from 3 chops were obtained, making sure that the cores ran parallel to the muscle fibers. Each core was sheared once using an Instron Universal Testing Machine (model 1132, Instron, Canton, MA) with a Warner-Bratzler shear attachment, a 50-kg load cell, and a cross-head speed of 50 mm/min. Shear force was calculated as the mean shear force of the 6 cores.
Collection of Genotypic Data
Genomic DNA was extracted from blood, liver, or spleen tissue of 9 Dorset grandsires, 8 Romanov grand-dams, and all F1 parents and F2 lambs using Gentra Generation Capture kits (Gentra Systems Inc., Minneapolis, MN) or a standard saturated salt procedure (Miller et al., 1988
). Sheep were genotyped at 6 microsatellite loci (Figure 1
) using the standard PCR protocol of de Gortari et al. (1998)
. Five of these markers were previously linkage-mapped to ovine chromosome 13 (Maddox et al., 2001
; http://rubens.its.unimelb.edu.au/~JILLM/jill.htm). The sixth marker, \PRNP_S11, was not linkage-mapped to chromosome 13, but was identified within the ovine genome sequence (Genbank Accession No. U67922) that contains the PRNP gene (Geldermann et al., 2003
).
|
To account for the effects of callipyge genotypes on carcass and meat quality traits, all sheep were genotyped at the callipyge locus using a matrix-assisted laser desorption/ionization time-of-flight mass spectrometer assay described by Freking et al. (2002)
. Animals were genotyped as CMatCPat, CMatNPat, NMatCPat, or NMatNPat at the callipyge locus, with C representing the mutant allele and N representing the wild-type allele. Superscripts denoted the maternal (Mat) or paternal (Pat) origin of each allele. The 2 classes of heterozygous sheep (CMatNPat and NMatCPat) were differentiated using flanking marker information for ovine chromosome 18 (Freking et al., 1998a
).
Linkage distances for chromosome 13 were calculated using CRIMAP version 2.4 (Washington University, St. Louis, MO), and compared with the reported linkage distances (http://rubens.its.unimelb.edu.au/~JILLM/jill.htm). Genotyping errors and deviations from the reported linkage distances were investigated using the CHROMPIC option of CRIMAP. The lack of genotypic information and the inability to resolve pedigree inconsistencies using genotypic data resulted in deletion of data collected on 13 F2 lambs. No recombinant gametes were detected within the interval containing PRNP_136, PRNP_154, PRNP_171, and \PRNP_S11 loci (Figure 1
), and haplotypes were determined based on phase information. The probability of grandparental origin (Dorset or Romanov) of the PRNP haplotype for each gamete was estimated using Genoprob (Thallman, 2002
). Following Genoprob analysis, the \PRNP_S11 locus was no longer considered part of the PRNP haplotypes, because this locus has no known association with scrapie susceptibility.
Five PRNP haplotypes of Dorset origin (ARR, ARQ, VRQ, AHQ, and ARH) and 3 of Romanov origin (ARR, ARQ, and VRQ) were segregating in the F2 population (Table 2
). Because of the low frequency of the Dorset ARH haplotype (n = 4), these data were removed from subsequent analyses. Seven probabilities were calculated for each gamete, each corresponding to the probability of inheriting 1 of the 7 breed-specific PRNP haplotypes. Probabilities summed to 2 within each lamb, as each F2 lamb contained 2 PRNP haplotypes. Breed of origin probabilities were distributed in a bimodal manner, with the majority of probabilities greater than 0.95 or less than 0.05. The small number of intermediate probabilities was most likely due to incomplete genotypic information or low marker informativeness for some F2 chromosomes.
|
Fitting the effect of PRNP haplotype in this manner, rather than including the assumed PRNP haplotype in the statistical model as a class variable (as in other studies), allows the PRNP haplotype to be weighted using the probabilities of breed origin. Probabilities located at the upper and lower extremes of the distribution will be weighted more heavily in the analyses than will intermediate probabilities. The power to detect true effects of the resistant ARR haplotype increases as the contribution of more uncertain genotypes (intermediate probabilities) is reduced.
A single contrast among regression coefficients of PRNP haplotype covariates was constructed within each breed of origin. For Dorset, the ARR coefficient was compared with the average of ARQ, VRQ, and AHQ (1, 0.33, 0.33, 0.33), whereas for Romanov, ARR was compared with ARQ and VRQ (1, 0.5, 0.5). Contrasts tested the difference between the resistant ARR haplotype and the average of the other haplotypes, consistent with the use of haplotype information by sheep industries in many countries.
Freking et al. (1999)
discussed the nonnormal distribution of Warner-Braztler shear force residuals in these data. Dispersion of residuals increased with greater values of shear force, indicating a nonnormal distribution. Shear force values were transformed by natural logarithm before statistical analysis to account for the proportional relationship between standard deviations and means of callipyge genotypes.
Meat quality traits were analyzed using a model that included the random effect of sire, fixed effects of birth year, sex, and callipyge genotype, and 7 covariates, as described above. For analysis of transformed Warner-Bratzler shear force data, callipyge genotypes were grouped into 3 categories: NMatCPat, CMatCPat, and NN (which consisted of NMatNPat and CMatNPat sheep). For analysis of longissimus marbling data, callipyge genotypes were grouped into 4 categories (CMatCPat, CMatN-Pat, NMatCPat, and NMatNPat). For shear force, a birth year x callipyge genotype interaction was included to account for the differing distributions of shear force values in sheep born in 1994 and 1995 (Freking et al., 1999
). A covariate of chilled carcass weight and the interaction between chilled carcass weight and callipyge genotype were included in the model for both meat quality traits. To test for associations between meat quality traits and PRNP haplotypes, breed-specific contrasts were constructed as described above.
| RESULTS |
|---|
|
|
|---|
Dorset and Romanov haplotype frequencies differed greatly (Table 3
), supporting the decision to include separate covariates for Romanov- and Dorset-derived haplotypes in statistical models. Because of this situation, the precision of estimates of the contrast of the ARR regression coefficient compared with the average of the other coefficients differed between breeds of origin. Using carcass traits as an example, there were 207 gametes of Dorset origin with the ARR haplotype, and 169 Dorset gametes with ARQ, VRQ, or AHQ haplotypes. For Romanov, there were 27 ARR gametes, compared with 305 ARQ or VQR gametes.
|
|
| DISCUSSION |
|---|
|
|
|---|
A number of studies have investigated associations between PRNP genotypes and economically important traits (Prokopova et al., 2002
; de Vries et al., 2004
, 2005
; Alexander et al., 2005
). However, only a single study (Brandsma et al., 2004
) found multiple associations of PRNP genotypes with traits. In the study by Brandsma et al. (2004)
, associations were detected between the ARR haplotype and breeding values for litter size and 135-d weight in Dutch Texel sheep. Sheep homozygous for the ARR haplotype had greater breeding values for litter size and lesser breeding values for 135-d weight than sheep with zero or one copy of the ARR haplotype. In contrast, the current study did not detect an association between PRNP haplotypes and growth traits. However, the differing breed composition of the 2 studies makes direct comparison difficult. The genetic background (breed) of the population under study may influence associations between PRNP haplotypes and economically important traits.
In general, it is unlikely that detection of unfavorable associations of specific PRNP haplotypes with phenotypic variation of pertinent traits will be a common occurrence in sheep. There are several reasons for this statement. Most traits of economic importance are largely quantitative in nature, although the number of QTL reported in the scientific literature continues to increase. The PRNP locus on chromosome 13 segregates independently of loci on the other 26 ovine chromosomes. Within chromosome 13, associations between loci depend on linkage disequilibrium, a temporary phenomenon. If associations are detected, the PRNP haplotype resistant to scrapie may in fact be favorably, rather than unfavorably, associated with phenotypic variation. Consistent with goals of the US national scrapie eradication program and experimental results to date, increasing the frequency of resistant PRNP haplotypes should be considered a selection goal of breed development programs.
| IMPLICATIONS |
|---|
|
|
|---|
| Footnotes |
|---|
2 The authors thank J. K. Carnahan (US Meat Animal Research Center, Clay Center, NE) for providing genotypes at prion codons 136, 154, and 171. ![]()
3 Current address: ASC 2118, Ferris State University, Big Rapids, MI 49307. ![]()
4 Corresponding author: leymaster{at}email.marc.usda.gov
Received for publication August 22, 2005. Accepted for publication November 13, 2005.
| LITERATURE CITED |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. J. Lupton ASAS CENTENNIAL PAPER: Impacts of animal science research on United States sheep production and predictions for the future J Anim Sci, November 1, 2008; 86(11): 3252 - 3274. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Sawalha, S. Brotherstone, N. R. Lambe, and B. Villanueva Association of the prion protein gene with individual tissue weights in Scottish Blackface sheep J Anim Sci, August 1, 2008; 86(8): 1737 - 1746. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Salaris, S. Casu, and A. Carta Investigating the relationship between the prion protein locus and udder morphology traits and milk yield in Sardinian sheep J Anim Sci, November 1, 2007; 85(11): 2840 - 2845. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. D. Snowder, L. D. Van Vleck, L. V. Cundiff, G. L. Bennett, M. Koohmaraie, and M. E. Dikeman Bovine respiratory disease in feedlot cattle: Phenotypic, environmental, and genetic correlations with growth, carcass, and longissimus muscle palatability traits J Anim Sci, August 1, 2007; 85(8): 1885 - 1892. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Sawalha, S. Brotherstone, W. Y. N. Man, J. Conington, L. Bunger, G. Simm, and B. Villanueva Associations of polymorphisms of the ovine prion protein gene with growth, carcass, and computerized tomography traits in Scottish Blackface lambs J Anim Sci, March 1, 2007; 85(3): 632 - 640. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |