Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

Association of bovine leptin polymorphisms with energy output and energy storage traits in progeny tested Holstein-Friesian dairy cattle sires



Leptin modulates appetite, energy expenditure and the reproductive axis by signalling via its receptor the status of body energy stores to the brain. The present study aimed to quantify the associations between 10 novel and known single nucleotide polymorphisms in genes coding for leptin and leptin receptor with performance traits in 848 Holstein-Friesian sires, estimated from performance of up to 43,117 daughter-parity records per sire.


All single nucleotide polymorphisms were segregating in this sample population and none deviated (P > 0.05) from Hardy-Weinberg equilibrium. Complete linkage disequilibrium existed between the novel polymorphism LEP-1609, and the previously identified polymorphisms LEP-1457 and LEP-580. LEP-2470 associated (P < 0.05) with milk protein concentration and calf perinatal mortality. It had a tendency to associate with milk yield (P < 0.1). The G allele of LEP-1238 was associated (P < 0.05) with reduced milk fat concentration, reduced milk protein concentration, longer gestation length and tended to associate (P < 0.1) with an increase in calving difficulty, calf perinatal mortality and somatic cells in the milk. LEP-963 exhibited an association (P < 0.05) with milk fat concentration, milk protein concentration, calving difficulty and gestation length. It also tended to associate with milk yield (P < 0.1). The R25C SNP associated (P < 0.05) with milk fat concentration, milk protein concentration, calving difficulty and length of gestation. The T allele of the Y7F SNP significantly associated with reduced angularity (P < 0.01) and reduced milk protein yield (P < 0.05). There was also a tendency (P < 0.1) for Y7F to associate with increased body condition score, reduced milk yield and shorter gestation (P < 0.1). A80V associated with reduced survival in the herd (P < 0.05).


Several leptin polymorphisms (LEP-2470, LEP-1238, LEP-963, Y7F and R25C) associated with the energetically expensive process of lactogenesis. Only SNP Y7F associated with energy storage. Associations were also observed between leptin polymorphisms and calving difficulty, gestation length and calf perinatal mortality. The lack of an association between the leptin variants investigated with calving interval in this large data set would question the potential importance of these leptin variants, or indeed leptin, in selection for improved fertility in the Holstein-Friesian dairy cow.


Leptin is a 16 kD non-glycosylated polypeptide that signals the body fat reserves and energy status to the hypothalamus. It is primarily produced by adipose tissue and acts to alter feed intake [1, 2]. Leptin binds to its receptor mainly localised on neuropeptide Y-neurons, which results in an increase in energy output and hypophagia [3, 4]. Leptin acts as a body barometer providing a critical link between energy homeostasis, appetite and reproductive function [46]. During early lactation, dairy cows are in a state of negative energy balance where fat stores and feed intake are directed towards the energetically expensive process of lactogenesis [7]. Reproduction, therefore, receives a low priority [5, 79]. With a noticeable decline in fertility in the modern dairy cow as a result of the intense selection for milk production alone [10], knowledge of combination allelic variations in the bovine leptin gene, or its receptor, may provide a mechanism to select for improved fertility without negatively affecting appetite or milk production.

The bovine leptin gene (LEP) maps to chromosome 4q32, is 16.735 kb in length, spans 3 exons and codes for a protein of 167 amino acids which includes a 21 amino acid signal sequence [11, 12]. The bovine LEP promoter stretches approximately 3 kb upstream [12]. To date, polymorphisms in the bovine LEP gene locus have been associated with serum leptin concentrations [13, 14], dry matter intake [15], feed intake [13, 1517], milk energy output [15, 1820] and energy storage [2124]. The pleiotropic effects of leptin are further evidenced by associations between LEP polymorphisms and immune function [19, 25], calf prenatal mortality [26], calving interval [27], carcass traits [13], lean meat yield [23] and growth rate [13, 28].

The bovine leptin receptor gene (LEPR) is located on chromosome 3q33 [29]. A polymorphism in exon 20 (T945M) associates with serum leptin levels in late pregnancy in Holstein-Friesian dairy cows [30]. Although no associations have been found with milk yield and feed intake in other cattle breeds [15], associations have been found with LEPR polymorphisms and subcutaneous fat deposition, carcass fat yield [31] and growth [32]. In addition, polymorphisms in LEPR in other species have been linked to reproductive performance [33] and body mass [34].

The present study aimed to quantify the associations between novel and previously reported single nucleotide polymorphisms (SNP) in the LEP and LEPR with performance traits in Holstein-Friesian sires.

Results and Discussion

SNP discovery and detection in the Holstein-Friesian study population

Sequence analysis of 1.656 kb of the regulatory region of LEP in 14 Holstein-Friesian sires, selected based on their predicted transmitting ability for body condition score (BCS) being greater than 2 standard deviation units each side of the base population, revealed a T to C substitution at position -1609 bp upstream of the transcription start site (LEP-1609).

A total of 848 Holstein-Friesian sires were genotyped for this novel and 9 previously identified SNPs (LEPR-T945M, LEP-2470, LEP-1457, LEP-1238, LEP-963, LEP-580, Y7F, R25C and A80V). The performance traits evaluated are not expressed by the sires themselves, but rather by their female progeny. Statistical analysis of the female progeny facilitates the estimation of genetic merit on individual sires. Systematic environmental effects are adjusted for and the random non-genetic variation associated with individual animal phenotypes is minimised, thereby resulting in a more accurate phenotype compared to the phenotype on the individual itself. This increased study power is particularly beneficial for low heritability traits such as fertility where the proportion of phenotypic variance attributable to additive genetic differences is low. The disadvantage of such a study design is that the performance traits included in this analysis are limited to traits routinely measured on progeny. In the present study the average number of progeny per sire was 842 daughter-parity records. The average co-ancestry among the 848 sires was 2.2%. When coupled with the mixed model methodology used and the de-regression of the predicted transmitting abilities (PTA), this implies that the associations reported herein are independent of pedigree structure.

Table 1 details the linkage disequilibrium (r2) between the LEP SNPs. Due to the complete linkage disequilibrium between phases among the SNPs LEP-1609, LEP-1457 and LEP-580, only LEP-1457 is reported in the association analysis. In addition, the strong linkage disequilibrium between SNPs LEP-1238, LEP-963 and R25C suggests that these SNPs are unlikely to act independently of each other. Indeed Banos et al. (2008) reported strong linkage disequilibrium between LEP-963 and R25C [15]. Table 2 lists the genotype and allele frequencies of LEPR and the remaining 7 LEP polymorphisms. All SNPs were segregating in this sample population and none deviated (P > 0.05) from Hardy-Weinberg equilibrium. The minor allele frequency of the SNPs varied from 0.04 to 0.48 (Table 2). Homozygous TT at position Y7F was not present in the dataset.

Table 1 Linkage disequilibrium (r2) between the SNPs in the promoter and coding region of the leptin gene
Table 2 Genotype frequencies, minor allele frequencies (MAF) and significance of deviation from Hardy-Weinberg equilibrium (HWE) for the leptin and leptin receptor SNPs

Association analysis

Allele substitution effects of the 8 SNPs (LEPR-T945M, LEP-2470, LEP-1457, LEP-1238, LEP-963, Y7F, R25C and A80V) on the energy intensive processes of reproduction and lactogenesis are detailed in Tables 3 and 4. Table 5 describes the associations between the 8 SNPs and energy storage. LEP polymorphisms exhibited associations with milk composition, angularity, calving difficulty, gestation length, calf perinatal mortality and survival in the herd. LEPR-T945M did not associate with any of the traits examined.

Table 3 Allele substitution effect (standard error in parenthesis) between the 8 SNPs and survival in the herd, reproduction and calving performance
Table 4 Allele substitution effect (standard error in parenthesis) between the 8 SNPs and milk production
Table 5 Allele substitution effect (standard error in parenthesis) between the 8 SNPs and body fatness


Leptin interacts with the reproductive axis at multiple sites [35] and is a key regulator in foetal development [36]. The G allele of LEP-1238 was associated (P < 0.05) with an increase in gestation length and tended (P < 0.10) to be associated with an increase in calving difficulty and calf perinatal mortality (Table 3). Previously, Brickell et al. reported an association in heifers between calf perinatal mortality and SNP UASMS1, a SNP in complete linkage disequilibrium with LEP-1238 [13, 26]. Leifers et al. observed that the G allele at LEP-1238 tended to associate with a quicker return to luteal activity postpartum [20]. However, in our data set calving interval is unaffected (Table 3). LEP-963 T allele was associated (P < 0.05) with easier calving and shorter gestation length (Table 3). Interestingly, Leifers et al. noted a longer interval from calving to 1st observed oestrus in cows with the T allele at LEP-963 [20].

The T allele of R25C, which associated with shorter gestation length and easier calving in our study (Table 3), was associated with lower serum leptin protein concentrations during late pregnancy and higher calf perinatal mortality in Holstein-Friesians [14, 26]. In humans, leptin is thought to play an inhibitory role in uterine contractions and has been suggested as a tocolytic agent for threatened preterm [37]. It is possible that where maternal leptin concentration is low, a preterm labour (short gestation) may occur resulting in an easier birthing event of a smaller birth weight calf with an increased risk of perinatal mortality.

The associations with gestation length remained significant or approached significance when each of the three SNPs (LEP-1238, LEP-963 and R25C) was included in the model as a class variable (additional file 1). Using this alternative analysis, calving difficulty also significantly associated with LEP-963 and tended to associate with R25C. To the authors knowledge this is the first analysis of LEP polymorphisms with either gestation length or calving difficulty. What role leptin plays in gestation length or a parturition event remains to be determined.

The T allele of LEP-2470 was associated (P < 0.05) with an increase in calf perinatal mortality (Table 3). When the genotype was included in the model as a class effect, this association approached significance (P < 0.1, see additional file 1). However, Brickell et al.[26] reported no association between this SNP and calf perinatal mortality in 381 Holstein-Friesian nulliparous heifers. Our study uses calf perinatal data from multiparious cows in addition to heifers, which may explain the discrepancy. Furthermore, despite the fewer number of records used in our study, the accuracy of the phenotype for calf mortality is likely to be greater since PTA was used here which is based on performance in daughter descendants.

Interestingly, there was no corresponding association observed between LEP-2470 and calving difficulty (maternal and direct) which suggests that fetopelvic incompatibility does not contribute to the mortality association observed with LEP-2470. Although LEP-2470 may not be the causative mutation, this SNP lies in a sequence conserved in human, rat and bovine leptin promoters and adjacent to a putative binding site for estrogen binding protein (E2BP) transcription factor. Allele specific differences in the E2BP binding matrix may influence estrogen regulation of LEP expression [38]. Previously, associations have been found with the T allele and increased serum leptin concentrations, backfat and growth rate in beef cattle [13, 28]. It is possible that the T allele is also linked to similar phenotypes in the dairy cow which may partition nutrients away from the developing fetus with a detrimental effect on calf perinatal survival [26].

None of the LEP polymorphisms selected in the present study associated with fertility as measured by calving interval i.e. the number of days between consecutive parturitions. In particular, no association was observed between SNP-1457 and calving interval (Table 3), even though AG heterozygotes of LEP-1457 commenced luteal activity 7.1 days earlier postpartum than AA homozygotes [20]. However, calving interval has a lower heritability than commencement of luteal activity postpartum [39] and is a more complex phenotype influenced not only by luteal activity but also fertilization, embryo implantation and embryo survival. In the present study calving interval and gestation length are not confounded. Gestation length in the present study is based on the direct effects of the sire on the calf itself while calving interval is based on the direct effects of the sire through its expression in his daughters (i.e., the maternal grandsire of the calf). There is a part-whole relationship between calving interval and maternal gestation length and as such PTAs for the latter are not calculated in Ireland.

In the present study, A80V associated with survival in the herd (Table 3). Survival is influenced by several factors including fertility, health, milk production and temperament. Previous work in other data sets have demonstrated associations between the T allele and increased serum leptin concentrations in late pregnancy [14], increased milk yield and protein yield [18], more rapid growth rates and heavier animals [22]. Although this SNP has no association with milk, carcass or fertility traits individually in our data set, there may be an additive effect resulting in the observed association with survival.

Milk Production

Although circulating leptin concentrations are low in early lactation [14], leptin is locally produced by the mammary gland [40] and appears to play important roles in mammary gland development [41] and lactogenesis [42, 43]. The T allele of LEP-2470 associated (P < 0.05) with reduced milk protein concentration and had a tendency to associate with increased (P < 0.10) milk yield (Table 4). In agreement with Banos et al.[15], the tendency of this SNP to associate with milk yield in Holstein-Friesian cattle does not persist in our study when the genotype is included in the model as a class effect. To the authors' knowledge, this is the first time this SNP has been tested and found to associate with milk protein.

SNPs LEP-1238, LEP-963 and R25C significantly associated with milk fat and milk protein concentration (Table 4). In addition, the G allele of LEP-1238 tended (P < 0.10) to associate with an increase in somatic cell score (SCS) in the milk and the T allele of LEP-963 tended to associate with reduced milk yield (P < 0.10). When the genotypes were included as class effects in the model, the association with milk fat concentration for each of the three SNPs approached significance (P < 0.10). However the association with protein concentration only approached significance for LEP-1238 (additional file 1). Although none of the three SNPs may be the causative mutation for the observed association to milk composition, potential mechanisms for this association may lie in the location of LEP-1238 in a putative binding site for CCAAAT/enhancer binding protein (C/EBP). The polymorphism may alter the binding affinity of C/EBP at this locus. C/EBP sites along the promoter, in other species, play an important role in regulating LEP transcription [44, 45]. In a dose dependant manner, leptin upregulates the lactogenic effect of prolactin, enhancing induction of fatty acid synthesis and milk protein gene expression [43].

Previous association analyses with LEP-1238 and LEP-963 did not find an association with milk composition or yield [15, 20]. In addition some conflicting evidence exists for the association of R25C with milk production. The T allele at R25C associated with increased milk yield, milk protein yield and decreased milk fat concentration [46]. CT heterozygous Holstein cows had greater yields of milk fat and milk protein compared to CC cows [19]. In other data sets no associations were found between R25C and milk production traits in either Polish Black and White or Holstein cattle [15, 18]. The discrepancy between our data and these studies could be due to a multitude of factors: 1) the frequency of the different alleles in the sample populations; 2) the experimental design (e.g., using estimated breeding value or actual phenotypes as the dependent variable) or the statistical models and data used (e.g., parity of the animal or period during which the phenotype was measured), 3) genotype by environment interaction or 4) the genetic background of the animals.

The rare T allele of SNP T7F associated with a decrease in protein yield (P < 0.05) and tended to associate with a decrease in milk yield (P < 0.10) (Table 4). However these associations did not hold when the genotype was analyzed as a class effect (additional file 1).

Energy Storage

Angularity in the dairy cow is a subjective measure of adiposity [47]. Only SNP T7F exhibited a significant association with angularity (P < 0.01; Table 5). It also tended to associate with BCS (P < 0.10), a trait genetically similar yet opposite to angularity (Berry et al., 2004). Indeed when the genotype was analyzed as a class effect, the significant association (P < 0.05) with angularity remained (additional file 1). This agrees with previous work which indicated an association between this SNP and BCS [15, 16, 23]. T7F results in a tyrosine to phenylalanine substitution at amino acid position 7 within the peptide signal sequence. This substitution may alter processing and export of leptin, disrupting its ability to report on levels of adiposity.

LEP-1238, LEP-963 and R25C did not associate with carcass fat, angularity or BCS (Table 5) in our study. Previously, LEP-1238 exhibited associations with energy balance, feed intake and body weight but there is some conflicting evidence about LEP-963 and these traits in various data sets [13, 15, 20]. The T allele of R25C was associated with significantly heavier calves at weaning, fatter carcasses, slower growth rates, and higher leptin mRNA levels in adipose tissue compared to the C allele, albeit in beef cattle [28, 4851, 23, 52].

Haplotype analysis

Using the SNPs in the promoter and coding region of LEP, 22 possible haplotypes were reconstructed; 16 haplotypes occurred in less than 1% in the population and were therefore grouped together for the analysis (Table 6). None of the haplotypes were associated with fertility. Haplotype C-2470G-1457G-1238C-963AY7FCR25CCA80V was associated with an increase in milk fat yield of 1.22 kg, milk protein yield of 0.97 kg and milk protein concentration of 0.0128 percentage units. No adverse associations with milk yield, calving interval, functional survival, carcass fat or calf perinatal mortality were observed. Interestingly, this haplotype contains a G allele at LEP-1238, a C allele at LEP-963 and a C allele at R25C, all of which associated with a decrease in protein concentration. The C allele at LEP-2470, which associated with an increase in protein concentration, therefore appears to overpower the influence of these other alleles on milk protein in this haplotype. The two major haplotypes C-2470G-1457C-1238T-963AY7FTR25CCA80V and C-2470A-1457G-1238C-963AY7FCR25CTA80V represent 62% of the population and only differ in genetic merit for milk protein concentration. Three minor haplotypes (12%, 11% and 3%) were associated with an increase in angularity compared to their major counterparts.

Table 6 Regression coefficients (standard errors in parenthesis) of the performance traits on the different haplotypes


Of the 10 LEP and LEPR SNPs tested, LEP-2470, LEP-1238, LEP-963, Y7F and R25C associated with the energetically expensive process of lactogenesis. Animals with the C-2470G-1457G-1238C-963AY7FCR25CCA80VLEP haplotype had a superior genetic merit for milk fat yield, protein yield and milk protein concentration. SNP Y7F associated with energy storage. Several LEP polymorphisms associated with calving difficulty, gestation length or calf perinatal mortality. However, the absence of an association between the LEP variants investigated with calving interval would question the potential importance of these variants, or indeed leptin itself, in selection for improved fertility in the Holstein-Friesian dairy cow. This information may be useful in marker assisted selection, or some derivative of such (e.g,. genomic selection placing greater prior emphasis on known QTLs), to increase the accuracy of selection, in particular milk composition, thereby increasing genetic gain.


DNA extraction, Genotyping and Study population

Genomic DNA from frozen semen straws was extracted using the following procedure: thawed semen was washed twice in phosphate buffered saline (pH 7.4). Cell pellets were harvested by centrifugation and re-suspended in 450 μL of pre-warmed extraction buffer (10 mM Tris pH 8, 10 mM EDTA pH8.0, 1% SDS, 100 mM NaCl). Fifteen μL of 2-Mercaptoethanol was added. Samples were incubated at 55°C for 15 minutes followed by the addition of 10 μL Proteinase K (20 mg/ml). Lysis occurred following an overnight incubation at 60°C. DNA was extracted using the Maxwell® instrument (Promega Corp., Madison, WT, USA), according to the manufacturer's instructions. For SNP discovery, 1.656 kb of the LEP regulatory region was amplified by polymerase chain reaction (PCR) from genomic DNA of 14 Holstein-Friesian sires using four pairs of oligonucleotide primers designed based on the published sequence (GenBank AB070368). Primer pair -1656F 5' GTGGATGCTACTGCCTCTATT 3' and -929R 5' GCCTGGTTGTTTTGCTTTTA 3', using an annealing temperature of 56°C with 2 mM MgCl2, amplified a product of 720 bp. Primer pair -1358F 5' AAGTCCCCTGTAGATGTTTTTATG 3' and -579R 5' CGGGTCCGTTTTGTTCAC 3', using an annealing temperature of 54°C with 1.5 mM MgCl2, amplified a product of 779 bp. Primer pair -630F 5' GCAACGCACGGGGCTATCAATG 3' and -446R 5' GCCTGGCCTGGAAAATCACACCT 3', using an annealing temperature of 68°C with 1.5 mM MgCl2, amplified a product of 206 bp. Primer pair -473F 5' CTTACCCCTCCACACCATCATCAA 3' and + 33R 5' GAGCCGGGCACTTACCT 3', using an annealing temperature of 60°C with 1.5 mM MgCl2, amplified a product of 522 bp. PCR amplicons were sequenced bi-directionally by Cogenics (North Carolina, U.S.A) and the data analysed using the Lazergene 6 suite of programs (DNAStar Inc., Madison WI, USA). Where sequencing reads were consistently poor, PCR amplicons were cloned in order to sequence each allele separately. For cloning, the TOPO Cloning kit (Invitrogen Ltd., Paisley, UK) was used according to the manufacturers' instructions and for each cloned amplicon, plasmid DNA from at least 4 clones was sequenced.

SNP genotyping analysis was performed by Sequenom® using the iPLEX Gold assay on a MassARRAY® Platform ( In total, 10 SNPs (LEPR-T945M, LEP-2470, LEP-1609, LEP-1457, LEP-1238, LEP-963, LEP-580, Y7F, R25C and A80V) were selected for genotyping in 848 Holstein-Friesian sires with progeny in Ireland. These sires originated either directly or indirectly (e.g., sons of sires) from international breeding programs and were representative of sires used in Ireland in recent years. Polymorphisms LEP-2470 (also referred to as UASMS-2) [13], LEP-1457 [20], LEP-1238 (also referred to as UASMS3) [13], LEP-963 [20] and LEP-580 [20] are designated according to their distance upstream from the transcription start site. SNP T7F is located in exon 2 and results in a tyrosine to phenylalanine substitution at amino acid position 7 within the peptide signal sequence [16]. The C to T SNP of R25C results in an arginine to cysteine substitution at amino acid 25 in the α-helix of the leptin protein [50]. A80V SNP in exon 3 substitutes an alanine for valine at amino acid position 80 [53]. The LEPR-T945M is a T to C polymorphism in exon 20 of LEPR which results in a threonine for methionine substitution at residue 945 [30]. Genomic DNA from twenty-five animals was genotyped in duplicate for each SNP. Concordance across the 10 SNPs and all duplicates was 100%. TESS ( was used to search for putative transcription factor binding sites.

Genotypic and phenotypic data

A measure of linkage disequilibrium r2[54] was calculated between each pair-wise leptin combination of the segregating SNPs. PHASE V2.1 [55, 56] was used to reconstruct haplotype probabilities. Daughter yield deviations (DYD) and PTA, as well as associated reliabilities for a range of performance traits evaluated by the Irish Cattle Breeding Federation in January 2009, were available for inclusion in the analysis. The range of phenotypic traits could be broadly categorised into milk energy output traits (milk production), reproduction traits (calving interval, calving difficulty, perinatal mortality and gestation length), and energy storage traits (angularity, BCS, carcass fat). Models used in genetic evaluations in Ireland, as well as variance components, are outlined by Berry et al[57]. DYDs for 305-day milk, fat and protein yield as well as geometric mean SCS (loge somatic cell count) are estimated in Ireland using a repeatability animal model across the first five lactations. DYDs expressed on the scale of PTA were used in this study. A DYD is the weighted average of the deviation of a bull's daughters from their contemporaries adjusted for merit of their dams. A PTA is the average genetic merit an animal is expected to transmit to its offspring. PTA for calving interval and survival are estimated using a multi-trait animal model, including data from the first three lactations. PTA for milk yield is used to adjust survival for differences in genetic merit of milk yield; hence, this survival trait is functional survival. PTAs for progeny carcass fat score, measured at slaughter, are estimated in a multi-trait animal model that includes carcass weight, carcass conformation, cow carcass weight, weaning weight, live-weight of the animal between 300 and 600 days of age, feed intake, and skeletal and muscular linear traits. Male progeny slaughtered between 300 and 1,200 days of age and female progeny slaughtered between 300 and 875 days of age are included in the evaluation of carcass fat. Genetic evaluations for angularity and BCS, both subjective measures of subcutaneous fat depots, are undertaken as part of a joint evaluation in the UK and Ireland. The PTAs are standardised to the mean and standard deviation of the base population. Perinatal mortality is a dichotomous variable scored by farmers as calf dead at birth or within 24 hours [58]; calving difficulty, also scored by farmers is scored on a scale of 1 (easy) to 4 (veterinary assistance required). Gestation length is the number of days from service to birth of the calf. Genetic merit for perinatal mortality and gestation length are estimated using univariate animal linear mixed models. Genetic merits for direct and maternal calving difficulty are estimated using an animal-dam linear mixed model.

PTAs were de-regressed using the procedure outlined by Berry et al[59] to remove the contribution of pedigree to the estimated genetic merit of the animal. Parental contribution to the reliability of each DYD or PTA was removed using the approach of Harris and Johnson [60] and only sires with a reliability, less parental contribution, of > 60% were retained for inclusion in the association analysis. A total of 742 sires fulfilled these criteria for inclusion in the analysis of milk, fat and protein yield as well as milk fat and protein concentration; the number of sires included in the association analysis with SCS, calving interval and survival was 703, 501, and 477 respectively. The number of sires included in the analysis of angularity, BCS and carcass fat was 521, 504 and 446, respectively. The number of sires for direct calving difficulty, maternal calving difficulty, perinatal mortality and gestation length was 575, 506, 201 and 614 respectively.


The association between each SNP and performance was quantified using weighted mixed models in ASREML [61] with individual included as a random effect, and average expected relationships among individuals accounted for through the numerator relationship matrix. Year of birth (divided into 5 yearly intervals) and percent Holstein of the individual sire were included as fixed effects in the model. In all instances the dependent variable was the DYD or deregressed PTA of interest, weighted by its respective reliability less the parental contribution. Genotype was included in the analysis as a continuous or a class variable coded as the number of copies of a given allele. In a separate analysis, the probability of each haplotype for the individual was included as fixed effect in the weighted animal mixed model with year of birth of the bull and Holstein proportion also included as fixed effects; the haplotype with the greatest frequency was not included in the model to avoid linear dependencies among haplotype effects. Similarly when genotype was included in the model as a class effect one of the genotypes was set to zero.


  1. 1.

    Henry BA, Clarke IJ: Adipose tissue hormones and the regulation of food intake. J Neuroendocrinol. 2008, 20: 842-9. 10.1111/j.1365-2826.2008.1730.x.

  2. 2.

    Konner AC, Klockener T, Bruning JC: Control of energy homeostasis by insulin and leptin: targeting the arcuate nucleus and beyond. Physiol Behav. 2009, 97: 632-8. 10.1016/j.physbeh.2009.03.027.

  3. 3.

    Fruhbeck G: Intracellular signalling pathways activated by leptin. Biochem J. 2006, 393: 7-20. 10.1042/BJ20051578.

  4. 4.

    Farooqi IS, O'Rahilly S: Leptin: a pivotal regulator of human energy homeostasis. Am J Clin Nutr. 2009, 89: 980S-984S. 10.3945/ajcn.2008.26788C.

  5. 5.

    Zieba DA, Szczesna M, Klocek-Gorka B, Williams GL: Leptin as a nutritional signal regulating appetite and reproductive processes in seasonally-breeding ruminants. J Physiol Pharmacol. 2008, 59 (Suppl 9): 7-18.

  6. 6.

    Bluher S, Mantzoros CS: Leptin in reproduction. Curr Opin Endocrinol Diabetes Obes. 2007, 14: 458-64.

  7. 7.

    Chagas LM, Bass JJ, Blache D, Burke CR, Kay JK, Lindsay DR, Lucy MC, Martin GB, Meier S, Rhodes FM, et al: Invited review: New perspectives on the roles of nutrition and metabolic priorities in the subfertility of high-producing dairy cows. J Dairy Sci. 2007, 90: 4022-32. 10.3168/jds.2006-852.

  8. 8.

    Chilliard Y, Delavaud C, Bonnet M: Leptin expression in ruminants: nutritional and physiological regulations in relation with energy metabolism. Domest Anim Endocrinol. 2005, 29: 3-22. 10.1016/j.domaniend.2005.02.026.

  9. 9.

    Peter AT, Vos PL, Ambrose DJ: Postpartum anestrus in dairy cattle. Theriogenology. 2009, 71: 1333-42. 10.1016/j.theriogenology.2008.11.012.

  10. 10.

    Dobson H, Smith R, Royal M, Knight C, Sheldon I: The high-producing dairy cow and its reproductive performance. Reprod Domest Anim. 2007, 42 (Suppl 2): 17-23. 10.1111/j.1439-0531.2007.00906.x.

  11. 11.

    Pomp D, Zou T, Clutter AC, Barendse W: Rapid communication: mapping of leptin to bovine chromosome 4 by linkage analysis of a PCR-based polymorphism. J Anim Sci. 1997, 75: 1427-

  12. 12.

    Taniguchi Y, Itoh T, Yamada T, Sasaki Y: Genomic structure and promoter analysis of the bovine leptin gene. IUBMB Life. 2002, 53: 131-5. 10.1080/15216540211465.

  13. 13.

    Nkrumah JD, Li C, Yu J, Hansen C, Keisler DH, Moore SS: Polymorphisms in the bovine leptin promoter associated with serum leptin concentration, growth, feed intake, feeding behavior, and measures of carcass merit. J Anim Sci. 2005, 83: 20-8.

  14. 14.

    Liefers SC, te Pas MF, Veerkamp RF, Chilliard Y, Delavaud C, Gerritsen R, van der Lende T: Association of leptin gene polymorphisms with serum leptin concentration in dairy cows. Mamm Genome. 2003, 14: 657-63. 10.1007/s00335-003-2275-y.

  15. 15.

    Banos G, Woolliams JA, Woodward BW, Forbes AB, Coffey MP: Impact of single nucleotide polymorphisms in leptin, leptin receptor, growth hormone receptor, and diacylglycerol acyltransferase (DGAT1) gene loci on milk production, feed, and body energy traits of UK dairy cows. J Dairy Sci. 2008, 91: 3190-200. 10.3168/jds.2007-0930.

  16. 16.

    Lagonigro R, Wiener P, Pilla F, Woolliams JA, Williams JL: A new mutation in the coding region of the bovine leptin gene associated with feed intake. Anim Genet. 2003, 34: 371-4. 10.1046/j.1365-2052.2003.01028.x.

  17. 17.

    Liefers SC, te Pas MF, Veerkamp RF, van der Lende T: Associations between leptin gene polymorphisms and production, live weight, energy balance, feed intake, and fertility in Holstein heifers. J Dairy Sci. 2002, 85: 1633-8. 10.3168/jds.S0022-0302(02)74235-5.

  18. 18.

    Madeja Z, Adamowicz T, Chmurzynska A, Jankowski T, Melonek J, Switonski M, Strabel T: Short communication: effect of leptin gene polymorphisms on breeding value for milk production traits. J Dairy Sci. 2004, 87: 3925-7. 10.3168/jds.S0022-0302(04)73531-6.

  19. 19.

    Chebel RC, Susca F, Santos JE: Leptin genotype is associated with lactation performance and health of Holstein cows. J Dairy Sci. 2008, 91: 2893-900. 10.3168/jds.2007-0891.

  20. 20.

    Liefers SC, Veerkamp RF, te Pas MF, Delavaud C, Chilliard Y, Platje M, van der Lende T: Leptin promoter mutations affect leptin levels and performance traits in dairy cows. Anim Genet. 2005, 36: 111-8. 10.1111/j.1365-2052.2005.01246.x.

  21. 21.

    Corva PM, Fernandez Macedo GV, Soria LA, Papaleo Mazzucco J, Motter M, Villarreal EL, Schor A, Mezzadra CA, Melucci LM, Miquel MC: Effect of leptin gene polymorphisms on growth, slaughter and meat quality traits of grazing Brangus steers. Genet Mol Res. 2009, 8: 105-16. 10.4238/vol8-1gmr556.

  22. 22.

    Kulig H, Kmiec M: Association between leptin gene polymorphisms and growth traits in Limousin cattle. Genetika. 2009, 45: 838-41.

  23. 23.

    Schenkel FS, Miller SP, Ye X, Moore SS, Nkrumah JD, Li C, Yu J, Mandell IB, Wilton JW, Williams JL: Association of single nucleotide polymorphisms in the leptin gene with carcass and meat quality traits of beef cattle. J Anim Sci. 2005, 83: 2009-20.

  24. 24.

    Yang D, Chen H, Wang X, Tian Z, Tang L, Zhang Z, Lei C, Zhang L, Wang Y: Association of polymorphisms of leptin gene with body weight and body sizes indexes in Chinese indigenous cattle. J Genet Genomics. 2007, 34: 400-5. 10.1016/S1673-8527(07)60043-5.

  25. 25.

    Asiamah PA, Bechtel DG, Waldner C, Buchanan FC: Effects of leptin Arg25Cys on peripheral mononuclear cell counts and antibody response to vaccination in beef cattle. Anim Genet. 2009, 40: 783-7. 10.1111/j.1365-2052.2009.01918.x.

  26. 26.

    Brickell JS, Pollott GE, Clempson AM, Otter N, Wathes DC: Polymorphisms in the bovine leptin gene associated with perinatal mortality in Holstein-Friesian heifers. J Dairy Sci. 2010, 93: 340-7. 10.3168/jds.2009-2457.

  27. 27.

    Almeida SE, Almeida EA, Moraes JCF, Weimer TA: Molecular Markers in the LEP gene and reproductive performance of beef cattle. Journal of Animal Breeding and Genetics. 2003, 120: 106-113. 10.1046/j.1439-0388.2003.00377.x.

  28. 28.

    Lusk JL: Association of single nucleotide polymorphisms in the leptin gene with body weight and backfat growth curve parameters for beef cattle. J Anim Sci. 2007, 85: 1865-72. 10.2527/jas.2006-665.

  29. 29.

    Pfister-Genskow M, Hayes H, Eggen A, Bishop MD: The leptin receptor (LEPR) gene maps to bovine chromosome 3q33. Mamm Genome. 1997, 8: 227-10.1007/s003359900398.

  30. 30.

    Liefers SC, Veerkamp RF, te Pas MF, Delavaud C, Chilliard Y, van der Lende T: A missense mutation in the bovine leptin receptor gene is associated with leptin concentrations during late pregnancy. Anim Genet. 2004, 35: 138-41. 10.1111/j.1365-2052.2004.01115.x.

  31. 31.

    Schenkel FS, Miller SP, Moore SS, Li C, Fu A, Lobo S, Mandell IB, Wilton JW: Association of SNPs in the leptin and leptin receptor genes with different fat depots in beef cattle. Proceedings of the 8th World Congress on genetics apllied to livestock production. 2006, Minas Gerais, Brazil, 03-80.

  32. 32.

    Guo Y, Chen H, Lan X, Zhang B, Pan C, Zhang L, Zhang C, Zhao M: Novel SNPs of the bovine LEPR gene and their association with growth traits. Biochem Genet. 2008, 46: 828-34. 10.1007/s10528-008-9197-z.

  33. 33.

    Chen CC, Chang T, Su HY: Characterization of porcine leptin receptor polymorphisms and their association with reproduction and production traits. Anim Biotechnol. 2004, 15: 89-102. 10.1081/ABIO-120037903.

  34. 34.

    Mutch DM, Clement K: Genetics of human obesity. Best Pract Res Clin Endocrinol Metab. 2006, 20: 647-64. 10.1016/j.beem.2006.09.006.

  35. 35.

    Moschos S, Chan JL, Mantzoros CS: Leptin and reproduction: a review. Fertil Steril. 2002, 77: 433-44. 10.1016/S0015-0282(01)03010-2.

  36. 36.

    Grisaru-Granovsky S, Samueloff A, Elstein D: The role of leptin in fetal growth: a short review from conception to delivery. Eur J Obstet Gynecol Reprod Biol. 2008, 136: 146-50. 10.1016/j.ejogrb.2007.06.021.

  37. 37.

    Wuntakal R, Hollingworth T: Leptin--a tocolytic agent for the future?. Med Hypotheses. 74: 81-2. 10.1016/j.mehy.2009.07.039.

  38. 38.

    Yi KW, Shin JH, Seo HS, Lee JK, Oh MJ, Kim T, Saw HS, Kim SH, Hur JY: Role of estrogen receptor-alpha and -beta in regulating leptin expression in 3T3-L1 adipocytes. Obesity (Silver Spring). 2008, 16: 2393-9. 10.1038/oby.2008.389.

  39. 39.

    Royal MD, Pryce JE, Woolliams JA, Flint AP: The genetic relationship between commencement of luteal activity and calving interval, body condition score, production, and linear type traits in Holstein-Friesian dairy cattle. J Dairy Sci. 2002, 85: 3071-80. 10.3168/jds.S0022-0302(02)74394-4.

  40. 40.

    Bartha T, Sayed-Ahmed A, Rudas P: Expression of leptin and its receptors in various tissues of ruminants. Domest Anim Endocrinol. 2005, 29: 193-202. 10.1016/j.domaniend.2005.03.010.

  41. 41.

    Hu X, Juneja SC, Maihle NJ, Cleary MP: Leptin--a growth factor in normal and malignant breast cells and for normal mammary gland development. J Natl Cancer Inst. 2002, 94: 1704-11.

  42. 42.

    Feuermann Y, Mabjeesh SJ, Shamay A: Mammary Fat Can Adjust Prolactin Effect on Mammary Epithelial Cells via Leptin and Estrogen. Int J Endocrinol. 2009, 2009: 427260-

  43. 43.

    Feuermann Y, Mabjeesh SJ, Shamay A: Leptin affects prolactin action on milk protein and fat synthesis in the bovine mammary gland. J Dairy Sci. 2004, 87: 2941-6. 10.3168/jds.S0022-0302(04)73425-6.

  44. 44.

    Mason MM, He Y, Chen H, Quon MJ, Reitman M: Regulation of leptin promoter function by Sp1, C/EBP, and a novel factor. Endocrinology. 1998, 139: 1013-22. 10.1210/en.139.3.1013.

  45. 45.

    Melzner I, Scott V, Dorsch K, Fischer P, Wabitsch M, Bruderlein S, Hasel C, Moller P: Leptin gene expression in human preadipocytes is switched on by maturation-induced demethylation of distinct CpGs in its proximal promoter. J Biol Chem. 2002, 277: 45420-7. 10.1074/jbc.M208511200.

  46. 46.

    Buchanan FC, Van Kessel AG, Waldner C, Christensen DA, Laarveld B, Schmutz SM: Hot topic: an association between a leptin single nucleotide polymorphism and milk and protein yield. J Dairy Sci. 2003, 86: 3164-6. 10.3168/jds.S0022-0302(03)73918-6.

  47. 47.

    Berry DP, Buckley F, Dillon P, Evans RD, Veerkamp RF: Genetic relationships among linear type traits, milk yield, body weight, fertility and somatic cell count in primiparous dairy cows. Ir J Agr Food Res. 2004, 43: 161-176.

  48. 48.

    Di Stasio L, Brugiapaglia A, Galloni M, Destefanis G, Lisa C: Effect of the leptin c.73T>C mutation on carcass traits in beef cattle. Anim Genet. 2007, 38: 316-7. 10.1111/j.1365-2052.2007.01595.x.

  49. 49.

    Nkrumah JD, Li C, Basarab JA, Guercio S, Meng Y, Murdoch B, Hansen C, Moore SS: Association of a single nucleotide polymorphism in the bovine leptin gene with feed intake, feed efficiency, growth, feeding behaviour, carcass quality and body composition. Canadian Journal of Animal Science. 2004, 84: 211-219.

  50. 50.

    Buchanan FC, Fitzsimmons CJ, Van Kessel AG, Thue TD, Winkelman-Sim DC, Schmutz SM: Association of a missense mutation in the bovine leptin gene with carcass fat content and leptin mRNA levels. Genet Sel Evol. 2002, 34: 105-16. 10.1186/1297-9686-34-1-105.

  51. 51.

    Kononoff PJ, Deobald HM, Stewart EL, Laycock AD, Marquess FL: The effect of a leptin single nucleotide polymorphism on quality grade, yield grade, and carcass weight of beef cattle. J Anim Sci. 2005, 83: 927-32.

  52. 52.

    DeVuyst EA, Bauer ML, Cheng FC, Mitchell J, Larson D: The impact of a leptin gene SNP on beef calf weaning weights. Anim Genet. 2008, 39: 284-6. 10.1111/j.1365-2052.2008.01730.x.

  53. 53.

    Haegeman A, Van Zeveren A, Peelman LJ: New mutation in exon 2 of the bovine leptin gene. Anim Genet. 2000, 31: 79-10.1111/j.1365-2052.2000.579-14.x.

  54. 54.

    Hill WG, Roberston A: Linkage Disequilibrium in finite populations. Theoretical and Applied Genetics. 1968, 38: 226-231. 10.1007/BF01245622.

  55. 55.

    Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73: 1162-9. 10.1086/379378.

  56. 56.

    Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-89. 10.1086/319501.

  57. 57.

    Berry DP, Shalloo L, Cromie AR, Veerkamp RF, Dillion P, Amer PR, Kearney JF, Evans RD, Wickham B: The economic breeding index: a generation on. I C B Federation. 2007, 29-34.

  58. 58.

    Mee JF, Berry DP, Cromie AR: Prevalence of, and risk factors associated with, perinatal calf mortality in pasture-based Holestin-Friesian cows. Animal. 2008, 2: 613-620. 10.1017/S1751731108001699.

  59. 59.

    Berry DP, Kearney F, Harris B: Genomic selection in Ireland. Interbull International Workshop. 2009, Uppsala, Sweden

  60. 60.

    Harris B, Johnson D: Approximate reliability of genetic evaluations under an animal model. J Dairy Sci. 1998, 81: 2723-8. 10.3168/jds.S0022-0302(98)75829-1.

  61. 61.

    Gilmour AR, Cullis BR, Welham SJ, Thompson R: ASREML Reference Manual. 2007, New South Wales Agriculture, Orange Agricultural Institute, Orange, NSW, Australia

Download references


The authors would like to acknowledge the Irish Cattle Breeding Federation (ICBF) for their valuable contribution to this project. The authors would also like to kindly acknowledge David Magee for his technical assistance. We gratefully acknowledge research funding provided by Teagasc under the National Development Plan, the Irish Dairy Research Trust and the Research Stimulus Fund (RSF-06-0353; RSF-06-0409) from the Department of Agriculture, Food and Fisheries, Ireland. BMK was in receipt of a Teagasc Walsh Fellowship.

Author information

Correspondence to Linda Giblin.

Additional information

Authors' contributions

LG designed and coordinated the study and participated in manuscript drafting. SB participated in the design of the study and drafting the manuscript. BMK performed sequence analysis. SMW coordinated and prepared DNA samples for genotyping. MJC reviewed sequencing data and helped draft the manuscript. DB provided statistical analysis and participated in manuscript drafting. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1:SNPs with phenotypic associations in Tables3,4and5are included in the model as a class effect. aReferent class is the homozygous genotype not presented for the respective SNP. Significance of genotype association with the performance variable: = P < 0.10; * = P < 0.05. Standard error in parenthesis. BCS: Body Condition Score. SD Units: Standard deviation Units (DOC 52 KB)

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

Giblin, L., Butler, S.T., Kearney, B.M. et al. Association of bovine leptin polymorphisms with energy output and energy storage traits in progeny tested Holstein-Friesian dairy cattle sires. BMC Genet 11, 73 (2010).

Download citation


  • Single Nucleotide Polymorphism
  • Milk Yield
  • Body Condition Score
  • Gestation Length
  • Somatic Cell Score