J. Anim Sci.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Schenkel, F. S.
Right arrow Articles by Williams, J. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schenkel, F. S.
Right arrow Articles by Williams, J. L.
J. Anim. Sci. 2005. 83:2009-2020
© 2005 American Society of Animal Science


ANIMAL GENETICS

Association of single nucleotide polymorphisms in the leptin gene with carcass and meat quality traits of beef cattle1

F. S. Schenkel*,2, S. P. Miller*, X. Ye*, S. S. Moore{dagger},3, J. D. Nkrumah{dagger}, C. Li{dagger}, J. Yu{dagger}, I. B. Mandell*, J. W. Wilton* and J. L. Williams{ddagger}

* University of Guelph, Guelph, Canada N1G-2W1; and {dagger} University of Alberta, Edmonton, Canada T6G-2P5; and and {ddagger} Roslin Institute, Roslin, United Kingdom EH25 9PS


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Studies with different populations are required to properly characterize the robustness of associations of polymorphisms in candidate genes with economically important traits across beef cattle populations before this sort of genetic information can be used efficiently in breeding and management decisions. The objective of this study was to evaluate the association of previously reported SNP in the bovine leptin gene with carcass and meat quality traits from a large sample of crossbred beef cattle. Five SNP (UASMS1, UASMS2, UASMS3, E2JW, and E2FB) were genotyped on 1,111 crossbred bulls, heifers, and steers. The measured traits included fat, lean, and bone yield (%) by partial rib dissection, grade fat, LM area, HCW, quality grade, LM i.m. fat, and tenderness evaluation of LM and semitendinosus muscle. Only four SNP were analyzed (UASMS1, UASMS2, E2JW, and E2FB), because UASMS1 and UASMS3 were completely linked. A uni-variate mixed-inheritance animal model was used to evaluate the association of either genotypes or haplo-types with the traits. The two leptin exon 2 SNP were associated with fat and lean yield and grade fat (E2JW, P < 0.01; E2FB, P < 0.05), and they interacted in their effect on LM tenderness (P < 0.01). The leptin promoter SNP were either not associated with any of the traits (UASMS2) or with fat yield only (UASMS1). Three hap-lotypes (TCAC, CCAT, TTAC) were at high frequency in the population (88%) and had similar effects on all the traits. Compared with the common haplotypes, one haplotype (CCTT) showed a significantly different effect on fat and lean yield and grade fat (P < 0.01), and one haplotype (TTTT) had a different effect on LM tenderness (P < 0.03). Therefore, important associations between SNP within the leptin gene with lean yield, fatness (fat yield and subcutaneous fat), and tenderness were detected. Results confirm some of the previously reported associations, but diverge with respect to others, showing that further efforts are required to validate some prospective associations.

Key Words: Beef Cattle • Carcass Traits • Haplotype Analysis • Leptin Gene • Meat Quality • Single Nucleotide Polymorphism


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Leptin, the hormone product of the obese (leptin) gene located on BTA 4 (Stone et al., 1996Go; Pomp et al., 1997Go), is considered to play a role in the regulation of appetite, energy partition, and body composition (Houseknecht et al., 1998Go; Baile et al., 2000Go). Leptin is synthesized and expressed predominantly by adipocytes (Houseknecht et al., 1998Go) and relates to the feedback system that regulates long-term body fat weight and composition (Hossner, 1998Go). Studies have found associations of serum leptin concentration with carcass adipose depots and carcass characteristics of beef cattle (Minton et al., 1998Go; Geary et al., 2003Go).

Leptin has been considered a candidate gene for performance, carcass, and meat quality traits in beef (Fitzsimmons et al., 1998Go; Buchanan et al., 2002Go; Lagonigro et al., 2003Go). Several SNP have been reported in the leptin gene (Buchanan et al., 2002Go; Lagonigro et al., 2003Go; Nkrumah et al., 2005Go). Associations of molecular polymorphisms within exon 2 (Buchanan et al., 2002Go; Nkrumah et al., 2004Go) or the promoter region (Crews et al., 2004Go; Nkrumah et al., 2005Go) of the leptin gene with carcass and meat quality traits recently were reported in beef cattle, with some associations not being consistently verified across studies.

Before this sort of genetic information can be used efficiently in breeding and management decisions, studies with different populations are required to properly characterize the robustness of the associations of leptin polymorphisms with economically important traits across beef cattle populations. The objective of this study was to evaluate the association of five previously reported SNP in the bovine leptin gene with carcass and meat quality traits from a large sample of crossbred beef cattle.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Cattle
The animals were commercially fed heifers (165), steers (231), and bulls (61) from industry sires, heifers (40), steers (375), and bulls (48) from the University of Guelph breeding project, and steers (191) from a University of Guelph feeding trial carried out at a feedlot in Rockwood, Ontario, Canada. The three sources of cattle were identified as Commercial, Elora, and Rockwood, respectively. Animals were crossbred, with breed composition formed by several breeds. The major contributing breeds were Angus (AN), Charolais (CH), Limousin (LI), and Simmental (SM). The average contribution of these four breeds to the breed composition of animals having any fraction of the mentioned breeds were 0.46, 0.50, 0.50, and 0.50 for AN, CH, LI, and SM, respectively, for Commercial cattle; 0.24, 0.36, 0.38, and 0.41 for Elora cattle; and 0.51, 0.53, 0.59, and 0.41 for Rockwood cattle. Thus, with the primary use of pure-bred sires, the Elora cattle represented more advanced crosses in dams, whereas the Commercial and Rock-wood cattle were more representative of first-generation crossbreds.

Rockwood and Commercial cattle represented AI breeding, as well as some natural service. In both cases, sires were known and had extended pedigree information. Elora and Rockwood animals were a result of ongoing research at the University of Guelph. Elora cattle originated from three breeding herds with exclusive AI breeding and coordinated sire use across herds and extended pedigree information. Breeding herds were located at the University’s Elora Beef Cattle Research Centre, New Liskeard Agricultural Research Station, and the Agriculture and Agri-food Canada Research Station in Kapuskasing. Cattle from the three herds were fed postweaning at the University of Guelph Elora Beef Cattle Research Centre feedlot and were involved in various postweaning trials. These trials and treatments were accommodated into the statistical analyses. Animals were slaughtered on the basis of a target commercial finishing endpoint of 8 mm of backfat thickness.

DNA Isolation, Polymorphism Detection, and Genotyping
The sources of DNA were frozen steaks, with the exception of 48 bulls from Elora, from which venous blood was taken. The DNA was isolated using the standard phenol/chloroform method (Sambrook et al., 1989Go; Hoelzel, 1992Go). The cell lysis buffer was modified as a 0.8 M concentration of urea, including 0.8 M urea, 2% SDS, 100 mM Tris·HCl, 200 mM NaCl, and 2 mM EDTA, pH 7.5.

Five SNP (E2FB, E2JW, UASMS1, UASMS2, and UASMS3) in the leptin gene were investigated. Two SNP, E2FB (Buchanan et al., 2002Go) and E2JW (Lagonigro et al., 2003Go, originally referred to as 252-SNP), were located within the leptin exon 2, whereas UASMS1, UASMS2, and UASMS3 (Nkrumah et al., 2005Go) were within the leptin promoter region.

The genotyping of each SNP was carried out using the 5' nuclease allelic discrimination assay on an ABI Prism 7700 sequence detector (Applied Biosystems, Inc., Foster City, CA). Details of the procedures were described by Nkrumah et al. (2005)Go. The DNA from a subset of the genotyped animals was sequenced across each polymorphism, and the sequence results were used to confirm the genotypes obtained by discrimination assays.

Phenotypic Information
Information on LM tenderness at 2 (SFL2), 7 (SFL7), 14 (SFL14), and 21 d (SFL21) postmortem and of semitendinosus muscle tenderness at 7 d (SFS7) postmortem, chemical fat (CF), grade fat (GFAT), quality grade (QG), LM area (LMA), lean (LEANYL), fat (FATYL), and bone (BONEYL) yield and HCW were available on most of the 1,111 genotyped animals (Table 1Go).


View this table:
[in this window]
[in a new window]
 
Table 1. Number of phenotypic records on the carcass and meat quality traits with corresponding means, standard deviations, and coefficients of variation
 
Warner-Bratzler shear force measurements (kg) were used as an objective method of assessing tenderness (Shackelford et al., 1999Go). Grade fat is the backfat thickness measurement taken at the 12th- and 13th-rib interface. Longissimus muscle area is the measure of the LM area at the 12th- and 13th-rib interface using a tracing of the muscle. Chemical fat was determined by the chemical analysis of a core meat sample that determined the percentage of i.m. fat. Lean, fat, and bone yield were determined by dissection of a four-bone rib section. Quality grade is the marbling grade used for grading in Canada, with most carcasses falling in one of three grades (A, AA, AAA). Because only a few carcasses were classified as Prime, those animals were combined with AAA carcasses for the analyses. A complete description of carcass measures available and methods used are given in Laborde et al. (2001)Go.

Commercial and Rockwood cattle were slaughtered at Better Beef Ltd. slaughter plant in Guelph, whereas the Elora cattle were slaughtered at the University of Guelph Meat Science Laboratory. All standardized carcass and meat quality measures were made at the University of Guelph Meat Science Laboratory, except for GFAT, HCW, LMA, and QG for Commercial and Rockwood cattle, which had these measures taken at Better Beef by University of Guelph meat scientists.

A total of 1,104, 1,111, 1,106, 1,068, and 1,109 animals had genotypes available for UASMS1, UASMS2, UASMS3, E2FB, and E2JW SNP, respectively. The genotypes for UASMS1 and UASMS3 were almost perfectly linked. Only three out of 1,104 genotypes for UASMS1 and UASMS3 did not match each other (that is, the C and T alleles in UASMS1 were not associated with the C and G alleles in UASMS3, respectively, in only three animals). Thus, UASMS3 was dropped from the analyses, and allele frequencies and association with the traits for UASMS1 would be extended to UASMS3.

The genotypes of all animals were used to determine the allelic frequencies. For the study of association between SNP in the leptin gene and carcass and meat quality traits, only animals with required phenotypic information and with genotypes available for all four SNP (UASMS1, UASMS2, E2JW, and E2FB) were used. The resulting number of records ranged from 720 for SFL2 to 924 for QG. Table 1Go gives the number of records, mean, SD, and CV of the analyzed traits.

Statistical Analyses
All analyses were performed using the statistical software SAS (SAS Inst., Inc., Cary, NC) and ASREML (Gilmour et al., 2000Go). Descriptive characteristics of quantitative traits were obtained using SAS PROC MEANS. Allele frequencies were tabulated and compared by {chi}2 analysis using SAS PROC FREQ.

Genotype Analyses.
Association of the genotypes with the traits was evaluated by genetic analysis using AS-REML, fitting a mixed inheritance model (SNP genotypes plus polygenic effects). The model included SNP genotypes as fixed effects:


[1]

where Yijklm is the trait measured in the mth animal of kth sex and lth slaughter group; u is the overall mean for the trait; Geni(j) is the effect of the ith genotype for jth SNP (UASMS1, UASMS2, E2JW, and E2FB) in the leptin gene; Sexk is the fixed effect of the kth sex (bull, heifer, and steer); Slgl is the fixed effect of the lth slaughter group (94 levels); ß1, ß2, ß3, and ß4 are the regression coefficients on breed composition of AN, CH, LI, and SM, respectively; Polm is the random additive genetic (polygenic) effect of the mth animal; and eijklm is the residual random effect associated with the mth animal.

Following Fernando et al. (1998)Go, as genotypes were known, the mixed-model equations of Henderson (1984)Go for Model [1] were used in the analyses. The additive relationship matrix based on the general pedigree was used to model the covariances among polygenic effects. Animals originated from 125 sires and all sires were known. With respect to the dams, 43% of the animals had dams identified. Average size of paternal half-sib families was 8.9. Percentages of sires with fewer than five, 6 to 10, 11 to 15, and more than 15 offspring were 36, 26.4, 27.2, and 10.4%, respectively.

Slaughter groups were defined as animals from the same source (Commercial or Rockwood), and with the same slaughter date or animals from Elora coming from the same trial and feed treatment, and killed in the same season (December to February, March to May, June to August, and September to November).

The repeated shear force measurements of LM across postmortem periods were analyzed individually within each period, as the average shear force over periods (SFLavg), and as the intercept and slope of the individual linear regression of shear force measurements on postmortem days. The effect of the four SNP in the leptin gene on quality grade was analyzed by {chi}2 analysis (PROC FREQ), as well as a linear trait using AS-REML, applying Model [1]. In this case, scores of 1, 2, and 3 were assigned to quality grades A, AA, and AAA, respectively.

To keep reasonable probability values for Type I error, two levels of tests were performed. For initial assessment of the results, an overall value of P < 0.05 ({alpha}) was used. For a more detailed review of the results, a modified Bonferroni correction was used (; Mantel, 1980Go) to account for the number of tests. The value of n was determined using a SNP-wise approach combined with grouping traits according to type (Ye, 2003Go). Traits were grouped into two groups as follows: carcass yield traits (LEANYL, FATYL, BONEYL, GFAT, LMA, and HCW) and meat quality traits (CF, QG, SFL2, SFL7, SFL14, SFL21, SFLavg, and SFS7). Because there were four SNP, n was equal to 24 (4 x 6) and 32 (4 x 8) for carcass and meat quality traits, respectively, with the corresponding modified Bonferroni-corrected significance levels of 0.010 and 0.009.

Initially, two-way interactions between SNP were fit into the model, but there was a significant interaction only between E2JW and E2FB SNP for shear force of LM. For all other traits, the interactions were dropped from the model. For LM shear force, the joint genotype effect of E2JW and E2FB were included in the model.

Variances were estimated from the data and assumed known for estimation and testing purposes. Probabilities associated with the Wald F-statistics output by ASREML were obtained using error degrees of freedom that accounted for the estimated fixed effects, but ignored the fact that variances were estimated. This, however, should not be a problem, because the number of records on all traits was relatively high and variances were estimated by translation invariant functions of the data by REML (Henderson, 1984Go).

Haplotype Analyses.
Association of the haplotypes for the SNP in the leptin gene and carcass and meat quality traits was evaluated by genetic analysis using AS-REML, applying Model [1] replacing genotype effects by regressions on haplotype probabilities. The haplotype probabilities were reconstructed using the algorithm and HAPROB software developed by Boettcher et al. (2004)Go. This software estimates probabilities of haplo-type combinations for members of half-sib families, given that genotypes are known for all siblings, but unknown for all parents. The accuracy of reconstruction of the halfsibs’ haplotypes by the HAPROB software is high. For instance, the accuracy varies from 64 to 94% for reconstruction of haplotypes of individuals from halfsib families of two to 10 offspring, when three loci with three alleles are considered (Boettcher et al., 2004Go).

Table 2Go shows the 16 possible haplotypes and their corresponding probabilities. The three most frequent haplotypes had a summed probability of 0.88. Therefore, there were many rare haplotypes, and the least probable ones were joined into one group, which was referred to as haplotype 10.


View this table:
[in this window]
[in a new window]
 
Table 2. Haplotype probabilities in the beef population
 

    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Allele Frequencies
The {chi}2 test for differences in allele frequency among breeds (animals with breed composition ≥ 5/8 for a given breed) was not significant for any SNP in the leptin gene (Table 3Go). However, for UASMS1, Simmental tended to have a lower frequency of the C allele than the other breeds (P = 0.11) and, for E2FB, Angus tended to have a lower frequency of the C allele than the other breeds (P = 0.15). The latter result is in line with that reported by Buchanan et al. (2002)Go, where Angus had a lower frequency of the C allele than Charolais and Simmental, and with that reported by Nkrumah et al. (2004)Go, where cattle from an Angus foundation line had a lower frequency of the C allele than cattle from other lines with different foundation breeds.


View this table:
[in this window]
[in a new window]
 
Table 3. Allele frequencies within breeds and in the entire beef population for the UASMS1, UASMS2, E2JW, and E2FB SNP in the leptin gene
 
Considering all genotyped cattle, the T allele was predominant over the C allele for UASMS1 and E2FB (57.6 vs. 42.4%, and 61.1 vs. 38.9%, respectively), whereas for UASMS2, the C allele was much more common than T (73.8 vs. 26.1%). The E2JW SNP showed the largest difference in allele frequencies in the population. For this SNP, the T allele was rare compared with the A allele (4.0 vs. 96.0%). These overall frequencies are in agreement to those reported by Nkrumah et al. (2005)Go for UASMS1, Crews et al. (2004)Go and Nkrumah et al. (2005)Go for UASMS2, Buchanan et al. (2002)Go and Nkrumah et al. (2004)Go for E2FB, and Lagonigro et al. (2003)Go for E2JW.

The frequencies of genotypes were in agreement with Hardy-Weinberg equilibrium (Falconer and Mackay, 1996Go) within all SNP (the probabilities of the {chi}2 tests for deviation from the equilibrium were equal to 0.75, 0.17, 0.98, and 1.00 for UASMS1, UASMS2, E2JW, and E2FB, respectively.

Equilibrium in genotypic frequencies when considering jointly two SNP was tested by a {chi}2 test of expected and observed frequencies of gametic types (Falconer and Mackay, 1996Go). The only two SNP whose genotypes were jointly in equilibrium were UASMS1 and E2JW. All other pairwise tests showed a significant disequilibrium (P < 0.01).

Genotype Analyses
Genotypes did not significantly influence LMA, BON-EYL, CF, HCW, SFS7, and QG (Table 4Go). Genotypes for E2JW and E2FB significantly influenced LEANYL, FATYL, and GFAT, whereas genotypes for UASMS1 (or, alternatively, UASMS3) significantly influenced FATYL. The analyses of LM shear force in each particular postmortem day, and as an average shear force over the postmortem days, showed a significant effect of the E2JW.E2FB genotypes (the interaction E2JW by E2FB was significant; Table 4Go).


View this table:
[in this window]
[in a new window]
 
Table 4. Tests for the association of single nucleotide polymorphisms in leptin gene with lean (LEANYL), fat (FATYL), and bone (BONEYL) yield; grade fat (GFAT); chemical fat (CF), LM area (LMA); HCW; LM shear force (SFL) at 2, 7, 14 and 21 d postmortem; average shear force (SFLavg); semitendinosus muscle shear force at 7 d postmortem (SFS7); and quality grade (QG) in the beef cattle population
 
Table 5Go shows the least squares means for UASMS1, E2JW, and E2FB genotypes and breeds with the corresponding significance levels of the Wald F-tests for LEA-NYL, FATYL, and GFAT. The same table also presents the estimated polygenic heritability for the traits. The heritabilities for FATYL, LEANYL, and GFAT were moderate to high (0.62, 0.52, and 0.45, respectively), which agree with expected values from literature (Mrode et al., 1990Go; Burrow et al., 2001Go; Utrera et al., 2004Go).


View this table:
[in this window]
[in a new window]
 
Table 5. Association of single nucleotide polymorphisms in the leptin gene with lean yield (LEANYL), fat yield (FATYL), and grade fat (GFAT) in the beef cattle population
 
For E2JW, there were only two animals with genotype TT, which were excluded from the analyses. Therefore, only solutions for genotypes AA and AT were obtained. The T allele was associated with less FATYL and GFAT and more LEANYL compared with the A allele. The estimated differences between the heterozygous and homozygous genotypes were –1.5%, –1.2 mm, and 1.9% for FATYL, GFAT, and LEANYL, respectively (P < 0.05 for all differences), corresponding to 0.29, 0.36, and 0.38 phenotypic SD of the corresponding traits, respectively.

For E2FB, the C allele was associated with less FA-TYL and GFAT and more LEANYL than the T allele. The estimated differences between the homozygous genotypes CC and TT were –1.9%, –1.0 mm, and 2.3% for FATYL, GFAT, and LEANYL (P = 0.09, 0.19, and 0.05, respectively). The heterozygous genotype had, however, similar FATYL, LEANYL, and GFAT to the homozygous TT genotype, indicating a large degree of dominance of T over C. Differences of the CC genotype and the heterozygous genotype were all significant (P < 0.05) and correspond to 0.43, 0.44, and 0.40 phenotypic SD of the corresponding traits, respectively.

For UASMS1, the C allele was associated with less FATYL, with the estimated difference between the homozygous genotypes CC and TT equal to –1.5% (P < 0.05). The heterozygous genotype had similar FATYL as the homozygous CC genotype, indicating a large degree of dominance of C over T. There was a trend (P < 0.15) that the C allele might be associated with less GFAT and more LEANYL. The estimated differences between the homozygous genotypes CC and TT were –0.6 mm and 1.4% for GFAT and LEANYL, respectively.

Differences between genotypes for E2JW were also significant considering the modified Bonferroni correction for multiple tests, which was not the case for E2FB and UASMS1, indicating stronger evidence for the association of E2JW genotypes with FATYL, GFAT, and LEANYL than for E2FB and UASMS1.

Table 5Go shows that there was a significant effect of breed on FATYL, GFAT, and LEANYL. As expected (Burrow et al., 2001Go), Angus was the fattest breed with the least LEANYL. Simmental had the least FATYL and GFAT, followed by Limousin and Charolais; however, Limousin showed the greatest LEANYL. There was no significant effect of sex on FATYL, GFAT, and LEANYL, likely because this effect was partially confounded with slaughter group, which had a highly significant effect (Table 5Go).

The analyses of LM shear force in each particular postmortem day showed that the joint E2JW.E2FB genotypes had a significant effect on tenderness. Table 6Go presents the least squares means for the E2JW.E2FB genotypes, with the corresponding significance levels for the Wald F-tests.


View this table:
[in this window]
[in a new window]
 
Table 6. Association of single nucleotide polymorphisms in the leptin gene with tenderness (LM shear force; SFL) in the beef cattle population measured at different postmortem days (2, 7, 14, and 21 d)
 
Genotype AT.TT was significantly associated with tougher LM. Differences in shear force between genotypes AT.TT and AT.CT were substantial (1.98, 1.12, 1.05, and 0.69 kg for SFL2, SFL7, SFL14, and SFL21, respectively). Estimates for genotype AT.CC were not obtained because there was only one animal with this genotype. Differences in shear force between genotypes AA.TT and AA.CT were smaller and mostly nonsignificant. The magnitude of the differences between AT.TT vs. AT.CT and AA.TT vs. AA.CT illustrates the interaction between the E2JW and E2FB SNP, where a larger difference exists for the AT E2JW SNP genotype. Estimates for genotypes AA.CT and AA.CC were not significantly different.

Table 6Go also gives the estimated polygenic heritabilities, which for SFL7 and SFL21 were lower than for other postmortem days, but were within the expected literature range (Burrow et al., 2001Go).

Results for average shear force over the four postmortem measures (SFLavg) are shown in Table 7Go. In addition to the genotypes’ least squares means, the means for breeds also are presented. Results for E2JW.E2FB genotypes were in line with those found within the different postmortem days, with the genotype AT.TT having the toughest LM over the entire postmortem period.


View this table:
[in this window]
[in a new window]
 
Table 7. Association of single nucleotide polymorphisms in the leptin gene with average tenderness of LM across different postmortem days (SFLavg) and with the intercept (SFLin) and slope (SFLsl) of the regression of tenderness measurements on days postmortem
 
To further investigate the effects of SNP genotypes on tenderness, a linear regression of the repeated measures of shear force on postmortem days for each animal was estimated and the individual intercepts (SFLin) and slopes (SFLsl) were analyzed by Model [1]. Results presented in Table 7Go show a significant effect of E2J-W.E2FB genotypes on the intercept, but not on the slope, indicating that E2JW.E2FB genotypes did not influence the effect of aging on beef tenderization. The heritabilities of the slope and intercept of LM tenderness on aging times (Table 7Go) indicate that the intercept is moderately heritable, whereas the slope has low heritability.

There was a trend for a breed effect on tenderness (Table 7Go), with Simmental having the toughest LM. Slaughter group had a highly significant effect on tenderness, and sex was not significant, likely due to the partial confounding of sex with slaughter group.

Assuming the estimated allele frequencies for E2FB, E2JW, and UASMS1, and using the estimated additive (a = 1/2 homozygous genotype 1 – 1/2 homozygous genotype 2) and dominance deviation (d = heterozygous genotype – [1/2 homozygous genotype 1 + 1/2 homozygous genotype 2]) effects for the alleles, the percentage of phenotypic variation explained by each polymorphism was calculated using standard formula (Falconer and Mackay, 1996Go): %V = 100 x (2pq [a + d(q – p)]2 + [2pqd]2)/{sigma} 2p), where %V is the percentage of phenotypic variation explained by the polymorphism, and is the phenotypic variance of the trait. For E2JW, as the TT genotype effect was not estimated, it was assumed either that the T allele shows complete dominance over C, or that the T allele has an additive effect only.

The %V for tenderness explained by genotypes E2JW and E2FB was determined considering SNP main effects only. The least squares means for E2JW genotypes (AA and AT) were 4.33 ± 0.21 and 4.70 ± 0.24 kg, respectively, and for E2FB genotypes (CC, CT, and TT) were 4.19 ± 0.27, 4.16 ± 0.23, and 4.95 ± 0.26 kg, respectively.

Genotypes for E2FB explained 3.8, 3.9, 3.7, and 8.4% of the phenotypic variance for FATYL, GFAT, LEA-NYL, and SFLavg, respectively. Genotypes for E2JW explained either 0.7, 1.0, 1.0, and 1.0%, or 0.6, 0.9, 0.9, and 0.9%, when either additive or complete dominance effects of the T allele were assumed. The UASMS1 SNP explained 3.6% of phenotypic variance for FATYL. Thus, the expected percentage of the phenotypic variation explained by E2JW was much smaller than that explained by E2FB for all traits, but this was mainly because the T allele in E2JW is rare. If, however, frequency of the T allele increased, for instance by selection, to 0.50 and assuming an additive effect and no change in phenotypic variance, then E2JW genotypes would explain 4.3, 6.3, 6.4, and 6.5% of the phenotypic variation for FATYL, GFAT, LEANYL, and SFLavg, respectively.

Haplotype Analyses
The linear effect of each of 10 haplotypes was estimated. There were three highly frequent haplotypes in the beef population (88% of all haplotypes) whose effects did not differ for any trait analyzed, even though they differ with respect to the alleles in all SNP, except E2JW. This may indicate an effect of another SNP linked to the four SNP or some degree of epistasis among the SNP within the same chromosome. The average effect of the three common haplotypes was used as a control and all other haplotypes were contrasted against this average.

Differences in allele effects within each SNP were obtained through a linear contrast of haplotype solutions, which differed by only one allele at a given SNP. Haplotype 10 was not used in these contrasts because it comprised the joint effect of several rare haplotypes.

None of the other haplotypes significantly differed from the three most frequent haplotypes in their effect on LMA, BONEYL, CF, HCW, SFS7, and QG (data not shown). For FATYL, GFAT, and LEANYL, the effect of Haplotype 7 (CCTT) was significantly different from the three most frequent haplotypes in the population as shown in Table 8Go. Replacing the three most frequent haplotypes by Haplotype 7 would significantly decrease FATYL and GFAT by –2.26% and –1.84 mm, respectively, and increase LEANYL by +2.42%, corresponding to 0.44, 0.55, and 0.48 phenotypic SD of the corresponding traits, respectively.


View this table:
[in this window]
[in a new window]
 
Table 8. Association of haplotypes for single nucleotide polymorphisms in the leptin gene with lean yield (LEANYL), fat yield (FATYL), and grade fat (GFAT) in the beef cattle population
 
Table 8Go also shows the estimates of differences in the effects of the alleles within each SNP. In agreement with the genotype analyses, the T allele for E2JW decreased FATYL and GFAT and increased LEANYL compared with the A allele. With respect to E2FB, the results also agreed with the genotype analyses, where the C allele decreased FATYL and GFAT and increased LEANYL compared with the T allele.

Differences in the effects of UASMS1 alleles on FA-TYL, GFAT, and LEANYL were all nonsignificant. Nevertheless, the genotype analyses showed that UASMS1 genotypes significantly influenced FATYL. The contrast between haplotype effects is, however, estimating only the additive linear effect of the alleles.

For LM shear force at 2 and 14 d postmortem and for the average shear force over the 21 d postmortem, the effect of Haplotype 8 (TTTT) was significantly different from the three most frequent haplotypes in the population as shown in Table 9Go. Replacing the three most frequent haplotypes with Haplotype 8 would significantly decrease tenderness, increasing SFL2, SFL14, and SFLavg by 1.06, 0.58, and 0.55 kg, corresponding to 0.62, 0.46, and 0.53 phenotypic SD of the corresponding measurements, respectively.


View this table:
[in this window]
[in a new window]
 
Table 9. Association of haplotypes for single nucleotide polymorphisms in the leptin gene with LM shear force in the beef cattle population
 
Table 9Go also shows the estimates of differences in the effects of the alleles within each SNP. In agreement with the genotype analyses, the T allele for E2JW increased toughness compared with the A allele. There were no significant differences between alleles for E2FB. The genotype analysis showed a significant interaction between the E2JW and E2FB SNP, which is not accounted for when estimating allele effects by contrasts between haplotype effects. With respect to the other two SNP, the allele differences also were non-significant.

Analysis of the intercept and slope of individual regressions of shear force measurements on days postmortem (data not shown) revealed that the effect of Haplotype 8 (TTTT) was significantly different from the effect of the most frequent haplotypes for the intercept (0.96 ± 0.39 kg), and that there was a significant difference between the alleles T and A (0.57 ± 0.30 kg) for E2JW, but not for the other SNP. There were no significant differences in the haplotype effects on the slope of the regressions.

Adjustment to Different End Points
Analyses also were carried out adjusting records for either a common HCW or a common slaughter age through the inclusion of a fixed linear regression on either HCW or slaughter age in Model [1] (data not shown). Results were, however, very similar to those from the analyses without adjustment (Tables 4Go to 7GoGoGo). For instance, the significance levels for the Wald F-tests for the effects of E2JW, E2FB, and UASMS1 genotypes on FATYL were equal to 0.010, 0.012, and 0.014 adjusting for HCW, and equal to 0.009, 0.019, and 0.020 adjusting for slaughter age. The least squares means for E2JW (AT and TT), E2FB (CC, CT, and TT), and UASMS1 (CC, CT, and TT) genotypes were 22.4 and 23.9%; 21.8, 24.0, and 23.7%; and 22.8, 22.3, and 24.4% respectively, adjusting for HCW. Adjusting for slaughter age, the same features were 22.9 and 24.4%; 22.4, 24.5, and 24.2%; and 23.4, 22.8, and 24.8%, respectively. Therefore, differences between genotypes were similar when two alternative endpoints were considered.


    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Results reported are statistical associations whose biological and genetic validity cannot be proven within the scope of this study. Although two different analyses (genotype and haplotype) were carried out, strengthening the validity of some associations found, the nature of the quantitative data analyzed, and the intrinsic difficulties of this type of analyses preclude strong inferences on the biological and genetic meaning of the associations. Nonetheless, the findings constitute important evidence that, together with other literature reports, could lead to new investigations to better understand these associations, which would include the exclusion of possible effects of other nearby polymorphisms in linkage disequilibrium with the SNP investigated in this study.

The association of E2FB with FATYL, GFAT, and LEANYL found in the current study corroborates the results reported by Nkrumah et al. (2004)Go, where E2FB was genotyped in 144 commercial cattle from five genetic lines with different foundation breeds. These authors concluded that animals carrying the T allele vs. the C allele produce carcasses with poorer grades and lower lean meat yields, but which do not differ in carcass marbling. Buchanan et al. (2002)Go also reported a significant E2FB genotype effect on grade fat and average fat (mean value of three measures of backfat thickness along the 12th rib), with the T allele associated with higher fat, but with no significant association with carcass marbling score. Crews et al. (2004)Go, however, did not find association of E2FB with carcass traits of 433 Charolais and Charolais-cross steers, which included backfat thickness and carcass marbling.

As in the current study, Crews et al. (2004)Go did not find a significant association of E2FB with HCW and LMA. These authors also did not find an association of E2FB with either DMI or residual feed intake, which agrees with the results of Nkrumah el al. (2004)Go, who reported only a trend (P > 0.10) of TT animals having positive residual feed intake.

Strong evidence of an association of E2JW with FA-TYL, GFAT, and LEANYL was found in the current study. In the original study that described this SNP, however, Lagonigro et al. (2003)Go reported a nonsignificant association of E2JW with the percentage of s.c. and ultrasound backfat thickness (at 10 mo of age) from 169 Holstein-Charolais F2 bull calves. The same authors also reported no significant association of E2JW with carcass i.m. fat and marbling score, which agrees with the results of the present study, showing no association of E2JW with either CF or QG.

Lagonigro et al. (2003)Go reported a significant association of E2JW with the average feed intake of bull calves from 6 to 12 mo of age, with the AT genotype having higher daily intake than the AA genotype. The findings of the current investigation showed, however, that AT animals had lower FATYL and GFAT and higher LEA-NYL than AA animals, which might create the expectation that AT animals would have lower feed intake, if growth rate was unchanged, which is contrary to the findings of Lagonigro et al. (2003)Go.

A significant association of E2FB and E2JW with tenderness of LM across aging times was found in the current investigation. In fact, these two SNP interacted in their effect on tenderness. No other results were found in the available literature on the effect of these two SNP on beef tenderness for comparison.

Phenotypic and genetic relationships between marbling and tenderness (measured by either shear force evaluation or taste panel) are not especially high, but show favorable direction (Bertrand et al., 2001Go), indicating that higher marbling is slightly associated with higher tenderness. Analyses of LM tenderness also were carried out adjusting records for either CF or slaughter age through the inclusion of a fixed linear regression on either CF or slaughter age in Model [1] (data not shown). Results were, however, very similar to those from the analyses without adjustment. For instance, the probabilities for Wald F-tests for the effects of the joint E2JW.E2FB genotypes on SFLavg were equal to 0.001, adjusting for either CF or slaughter age. The least squares means for E2JW.E2FB genotypes (AA.CC, AA.CT, AATT, AT.CT, and AT.TT) were 4.12, 4.23, 4.50, 3.99, and 5.28 kg, adjusting for CF, respectively. The same features adjusting for slaughter age were 4.14, 4.23, 4.49, 3.99, and 5.31 kg, respectively.

Results indicate that the estimated genotype effects on tenderness were not due to differences in marbling or age of the animals. An association of the leptin receptor gene with pork meat tenderness was reported by Choi et al. (2003)Go. These authors found an association between microsatellite polymorphisms within the leptin receptor gene and LM shear force of 354 Korean Native x Landrace F2 boars at 12 wk of age.

In the present study, UASMS2 was not significantly associated with any of the carcass and meat quality traits considered. Nevertheless, Crews et al. (2004)Go reported a significant association of UASMS2 with carcass marbling score, LM area, and HCW of 433 Charolais and Charolais-cross steers, and Nkrumah et al. (2005)Go reported significant association of UASMS2 with ultrasound backfat thickness and marbling score of 150 crossbred animals (131 steers and 19 bulls). Contrasting results, such as these across studies, illustrate the need to validate associations across different populations before adoption is practical in widespread industry breeding programs. Crews et al. (2004)Go and Nkrumah et al. (2005)Go also reported a significant association of UASMS2 with daily DM intake and residual feed intake, indicating that UASMS2 might be associated with feed efficiency.

The UASMS1 SNP (or, alternatively, UASMS3) was significantly associated with FATYL and tended to have relationships with GFAT and LEANYL. In agreement with the present results, Nkrumah et al. (2005)Go reported a significant association of UASMS3 SNP with ultrasound backfat thickness. The same authors also reported an association of UASMS3 with daily DMI and BW.

Our findings confirm the association previously reported in the literature of the E2FB leptin exon 2 SNP with carcass lean meat yield and fatness (fat yield and grade fat). This SNP alone explains approximately 4% of the phenotypic variation for these traits. The T allele is associated with lower lean meat yield and higher fatness; however, the increased fatness does not translate into higher i.m. fat (marbling).

The other SNP in the leptin exon 2, E2JW, also is associated with carcass lean meat yield and fatness (fat yield and grade fat), with no association with i.m. fat. The T allele is implicated in higher lean meat yield and lower fatness. This SNP alone explains approximately 1% of the phenotypic variation for these traits. This allele is, however, quite rare, and the variance associated with this SNP increases to 4 to 6% of phenotypic variance if the frequency of the T allele was to increase to 50%, with an assumed additive effect.

The E2JW and E2FB polymorphisms are associated with tenderness of LM, and they interact in their effect. Individually, these two SNP explain approximately 1 and 8% of the phenotypic variation on tenderness, respectively. As mentioned before, the E2JW T allele is quite rare. The variance associated with this SNP would increase substantially as the T allele moves to a more intermediate frequency. Nevertheless, no other studies on the effect of these SNP on tenderness of beef cattle have been published, and no QTL for tenderness in BTA 4 have been reported. Clearly, more investigation is needed to confirm the results for tenderness found in this investigation to exclude, for instance, possible effects of other genes in linkage disequilibrium with the leptin SNP.

Two SNP in the leptin promoter, UASMS1 and UASMS3, are completely linked in the population and are significantly associated with fat yield. Another leptin promoter polymorphism, UASMS2, is not significantly associated in this population with any carcass and meat quality traits analyzed, which disagrees with two previously reported studies on this polymorphism.

Three particular haplotypes (TCAC, CCAT, and TTAC) within the leptin gene are highly frequent in the population and do not differ in their effects on carcass and meat quality traits, even though they carry different alleles. This might indicate the effect of other SNP linked to the four SNP considered in this study or some degree of epistasis among the SNP within the same chromosome.


    Implications
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Important associations between two single nucleotide polymorphisms within the leptin exon 2 and lean yield were detected. Current genetic improvement schemes that use estimated breeding values from ultrasound carcass data could potentially benefit through the use of these associations. The significant associations of those polymorphisms with tenderness are even more interesting, as current improvement programs for this potentially important trait are quite limited. Results confirm some of the previously reported associations, but diverge with respect to others, showing that further efforts are required to validate some prospective associations.


    Footnotes
 
1 Researchers from Univ. of Guelph acknowledge Beef Improvement Ontario (BIO) and Ontario Cattlemen’s Association for providing data; BIO, Natural Sciences and Engineering Research Council of Canada, and Ontario Ministry of Agriculture and Food for financial support; and the Canadian Foundation for Innovation and Ontario Innovation Trust for support of required computing infrastructure. Researchers from Univ. of Alberta acknowledge the financial support from the Canada/Alberta Beef Industry Development Fund through Grant No. 2000AB364 awarded to S. S. Moore. We thank P. J. Boettcher for his contributions to the haplotype analysis. Back

3 The author holds a U.S. patent application (currently under review) on the association of the three SNP in the leptin gene promoter (UASMS1, UASMS2, and UASMS3) used in this study with fatness traits in cattle. Back

2 Correspondence: Dept. of Anim. and Poultry Sci., Room 018 (phone: 519-824-4120, ext. 58650; fax: 519-767-0573; e-mail: schenkel{at}uoguelph.ca).

Received for publication January 5, 2005. Accepted for publication June 3, 2005.


    Literature Cited
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 


Baile, C. A., M. A. Della-Fera, and R. J. Martin. 2000. Regulation of metabolism and body fat mass by leptin. Annu. Rev. Nutr. 20:105–127.[Medline]

Bertrand, J. K., R. D. Green, W. O. Herring, and D. W. Moser. 2001. Genetic evaluation of beef carcass traits. J. Anim. Sci. 79(E. Suppl.):E190–E200.[Abstract/Free Full Text]

Boettcher, P. J., G. Pagnacco, and A. Stella. 2004. A Monte Carlo approach for estimation of haplotype probabilities in half-sib families. J. Dairy Sci. 87:4303–4310.[Abstract/Free Full Text]

Buchanan, F. C., C. J. Fitzsimmons, A. G. Van Kessel, T. D. Thue, D. C. W. Sim, and S. M. Schmutz. 2002. Association of a missense mutation in the bovine leptin gene with carcass fat content and leptin mRNA levels. Genet. Sel. Evol. 34:105–116.[Medline]

Burrow, H. M., S. S. Moore, D. J. Johnston, W. Barendse, and B. M. Bindon. 2001. Quantitative and molecular genetic influences on properties of beef: A review. Aust. J. Exp. Agric. 41:893–919.

Choi, B. H., T. H. Kim, Y. M. Cho, H. Y. Lee, J. T. Jeon, and I. C. Cheong. 2003. Association study between porcine LEPR-derived microsatellite polymorphisms and economic traits. J. Anim. Sci. Tech. 45:679–688.

Crews, D. H., J. D. Nkrumah, J. Yu, and S. S. Moore. 2004. Association of single nucleotide polymorphisms in the bovine leptin gene with feedlot and carcass characteristics of crossbred steers. Can. J. Anim. Sci. 84:749–750. (Abstr.)

Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to Quantitative Genetics. Longman Group Ltd., Essex, U.K.

Fernando, R. L., C. Stricker, and T. Wang. 1998. Detection and utilization of single genes without DNA assays. J. Dairy Sci. 81:64–75.[Abstract]

Fitzsimmons, C. J., S. M. Schmutz, R. D. Bergen, and J. J. McKinnon. 1998. A potential association between the BM1500 microsatellite and fat deposition in beef cattle. Mamm. Genome 9:432–434.[Medline]

Geary, T. W., E. L. McFadin, D. M. MacNeil, E. E. Grings, R. E. Short, R. N. Funston, and D. H. Keisler. 2003. Leptin as a predictor of carcass composition in beef cattle. J. Anim. Sci. 81:1–8.[Abstract/Free Full Text]

Gilmour, A. R., B. R. Cullis, S. J. Welham, and R. Thompson. 2000. ASREML Reference Manual. IACR-Rothamsted Experimental Station, Harpenden, U.K.

Henderson, C. 1984. Applications of Linear Models in Animal Breeding. Univ. of Guelph, Guelph, Ontario, Canada.

Hoelzel, A. R. 1992. Molecular Genetic Analysis of Populations: A Practical Approach. Oxford Univ. Press, Oxford, U.K.

Hossner, K. L. 1998. Cellular, molecular and physiological aspects of leptin: Potential applications in animal production. Can. J. Anim. Sci. 78:463–472.

Houseknecht, K. L., C. A. Baile, R. L. Matteri, and M. E. Spurlock. 1998. The biology of leptin: A review. J. Anim. Sci. 76:1405–1420.[Abstract/Free Full Text]

Laborde, F. L., I. B. Mandell, J. J. Tosh, J. W. Wilton, and J. G. Buchanan-Smith. 2001. Breed effects on growth performance, carcass characteristics, fatty acid composition, and palatability attributes in finishing steers. J. Anim. Sci. 79:355–365.[Abstract/Free Full Text]

Lagonigro, R., P. Wiener, F. Pilla, J. A. Woolliams, and J. L. Williams. 2003. A new mutation in the coding region of the bovine leptin gene associated with feed intake. Anim. Genet. 34:371–374.[Medline]

Mantel, N. 1980. Assessing laboratory evidence for neoplastic activity. Biometrics 36:381–399.[Medline]

Minton, J. E., D. J. Bindel, J. S. Droullard, E. C. Titgemeyer, D. M. Grieger, and C. M. Hill. 1998. Serum leptin is associated with carcass traits in finishing cattle. J. Anim. Sci. 76(Suppl.):231. (Abstr.)

Mrode, R. A., C. Smith, and R. Thompson. 1990. Selection for rate and efficiency of lean gain in Hereford cattle. II. Evaluation of correlated responses. Anim. Prod. 51:35–46.

Nkrumah, J. D., C. Li, J. B. Basarab, S. Guercio, Y. Meng, B. Murdoch, C. Hansen, and S. S. Moore. 2004. Association of a single nucleotide polymorphism in the bovine leptin gene with feed intake, feed efficiency, growth, feed behaviour, carcass quality and body composition. Can. J. Anim. Sci. 84:211–219.

Nkrumah, J. D., C. Li, J. Yu, C. Hansen, D. H. Keisler, and S. S. Moore. 2005. Polymorphisms in the bovine leptin promoter associated with serum leptin concentration, growth, feed intake, feeding behavior, and measures of carcass merit. J. Anim. Sci. 83:20–28.[Abstract/Free Full Text]

Pomp, D., T. Zou, A. C. Clutter, and W. Barendse. 1997. Rapid communication: mapping of leptin to bovine chromosome 4 by linkage analysis of a PCR-based polymorphism. J. Anim. Sci. 75:1427.[Free Full Text]

Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular Cloning, a Laboratory Manual. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY.

Shackelford, S. D., T. L. Wheeler, and M. Koohmaraie. 1999. Evaluation of slice shear force as an objective method of assessing beef longissimus tenderness. J. Anim. Sci. 77:2693–2699.[Abstract/Free Full Text]

Stone, R. T., S. M. Kappes, and C. W. Beattie. 1996. The bovine homologue of the obese gene maps to chromosome 4. Mamm. Genome 7:399–400.[Medline]

Utrera, A. R., and L. D. Van Vleck. 2004. Heritability estimates for carcass traits of cattle: A review. Genet. Mol. Res. 3:380–394.[Medline]

Ye, X. 2003. Identification of quantitative trait loci in swine using a candidate gene approach. Ph.D. Diss., Univ. of Guelph, Guelph, Ontario, Canada.


This article has been cited by other articles:


Home page
J ANIM SCIHome page
J. L. Lusk
Association of single nucleotide polymorphisms in the leptin gene with body weight and backfat growth curve parameters for beef cattle
J Anim Sci, August 1, 2007; 85(8): 1865 - 1872.
[Abstract] [Full Text] [PDF]


Home page
Poult. Sci.Home page
M. Lei, C. Luo, X. Peng, M. Fang, Q. Nie, D. Zhang, G. Yang, and X. Zhang
Polymorphism of Growth-Correlated Genes Associated with Fatness and Muscle Fiber Traits in Chickens
Poult. Sci., May 1, 2007; 86(5): 835 - 842.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
B. S. Sharma, I. Leyva, F. Schenkel, and N. A. Karrow
Association of toll-like receptor 4 polymorphisms with somatic cell score and lactation persistency in holstein bulls.
J Dairy Sci, September 1, 2006; 89(9): 3626 - 3635.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
R. Bergen, S. P. Miller, J. W. Wilton, and I. B. Mandell
Genetic correlations between live yearling bull and steer carcass traits adjusted to different slaughter end points. 2. Carcass fat partitioning.
J Anim Sci, March 1, 2006; 84(3): 558 - 566.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Schenkel, F. S.
Right arrow Articles by Williams, J. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schenkel, F. S.
Right arrow Articles by Williams, J. L.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS