The effects of polymorphisms in theDGAT1,leptinandgrowth hormone receptorgene loci on body energy, blood metabolic and reproductive traits of Holstein cows

Genetic analysis of milk citrate predicted by milk mid-infrared spectra of Holstein cows in early lactation

Milk citrate is regarded as an early biomarker of negative energy balance in dairy cows during early lactation and serves as a suitable candidate phenotype for genomic selection due to its wide availability across a large number of cows through milk mid-infrared spectra prediction. However, its genetic background is not well known. Therefore, the objectives of this study were to (1) analyze the genetic parameters of milk citrate; (2) identify genomic regions associated with milk citrate; and (3) analyze the functional annotation of candidate genes and quantitative trait loci (QTL) related to milk citrate in Walloon Holstein cows. In total, 134,517 test-day milk-citrate phenotypes (mmol/L) collected within the first 50 d in milk on 52,198 Holstein cows were used. These milk-citrate phenotypes, predicted by milk mid-infrared spectra, were divided into 3 traits according to the first (citrate1), second (citrate2), and third to fifth parity (citrate3+). Genomic information for 566,170 SNPs was available for 4,479 animals. A multiple-trait repeatability model was used to estimate genetic parameters. A single-step GWAS was used to identify candidate genes for citrate and post-GWAS analysis was done to investigate the relationship and function of the identified candidate genes. The heritabilities estimated for citrate1, citrate2, and citrate3+ were 0.40, 0.37, and 0.35, respectively. The genetic correlations among the 3 traits ranged from 0.98 to 0.99. The genomic correlations among the 3 traits were also close to 1.00 across the genomic regions (1 Mb) in the whole genome, which means that citrate can be considered as a single trait in the first 5 parities. In total, 603 significant SNPs located on 3 genomic regions (chromosome 7, 68.569-68.575 Mb; chromosome 14, 0.15-1.90 Mb; and chromosome 20, 54.00-64.28 Mb), were identified to be associated with milk citrate. We identified 89 candidate genes including GPT, ANKH, PPP1R16A, and 32 QTL reported in the literature related to the identified significant SNPs. These identified QTL were mainly reported associated with milk fatty acids and metabolic diseases in dairy cows. This study suggests that milk citrate in Holstein cows is highly heritable and has the potential to be used as an early proxy for the negative energy balance of Holstein cows in a breeding objective.

Candidate K232A mutation of DGAT1 gene associated with production and reproduction traits in Indian Dairy cattle

The present study was conducted on Hardhenu and Sahiwal cows to evaluate the genetic diversity and the effect of candidate K232A mutation of diacylglycerol O-acyltransferase (DGAT1) gene on performance traits. PCR product of 413 bp size targeting K232A mutation of DGAT 1 gene was genotyped by Eae1 restriction enzyme and further confirmed by sequencing. The PCR-RFLP analysis for K232A showed that the frequencies of alleles K and A were 0.80 and 0.20 in Hardhenu and 0.98 and 0.02 in Sahiwal cattle, respectively. The frequencies of genotypes for KK and KA genotype were estimated as 0.60 and 0.40 in Hardhenu and 0.96 and 0.04 in Sahiwal cattle, respectively. Chi-square test showed that K232A SNP did not meet with the Hardy-Weinberg equilibrium (p < .01) in Hardhenu cattle. Association analysis revealed that KA genotype was significantly associated with total milk yield (p < .01) and 305 days milk yield (p < .01)in Hardhenu cattle. Hardhenu cows having KA genotype comparatively revealed favorable mean values for AFS (880.31 ± 26.65), AFC (1166.36 ± 19.94), SP (157.46 ± 22.80), CI (441.41 ± 44.28) and AI/Conception (1.80 ± 0.17), indicating that KA genotype can be used as a candidate marker for simultaneous improvement of both production and reproduction traits.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Genetic structure of some candidate genes of repeat breeder syndrome in Egyptian buffaloes

Background

This study aimed to explore the association between polymorphisms in three genes: leptin (LEP), leptin receptor (LEPR), and BMP4, and incidence of repeat breeding in Egyptian buffaloes.

Methods

DNA was extracted from 160 female buffaloes, involving 108 fertile and 52 repeat breeders. Genotyping was performed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Sequence analysis and alignment were performed by employing NCBI/BLAST/blastn suite, to identify SNPs among different patterns and alleles. We utilized PredictSNP software to predict the non-synonymous SNPs influences on protein function. Moreover, the conservation score of the amino acids within the target proteins was computed by ConSurf server.

Results

The genotyping results showed that LEP and BMP4 genes were monomorphic (CC, GG) in all tested fertile and repeat breeder buffaloes. Leptin gene sequencing showed a non-synonymous C73T SNP, replacing R to C at position 25 within the leptin polypeptide (position 4 in the mature form; R4C) which is a neutral mutation, not affecting function or structure of LEP protein. For LEPR, one synonymous SNP (T102C) and two non-synonymous SNPs (A106G and C146A), triggering V967A and G954C replacements, respectively in LEPR protein. Moreover, they are neutral mutations. Sequencing results of BMP4 showed HinfI restriction site indicate fixed GG genotype (CC genotype in the anti-sense strand) in all sequenced samples. No SNPs were observed within the amplified region.

Conclusion

Genotyping and sequencing results of the surveyed three genes revealed that there is no association between these genes mutations and the incidence of repeat breeding in Egyptian buffaloes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00397-2.

Single nucleotide polymorphisms of leptin gene in five Ethiopian indigenous cattle breeds and the Korean Hanwoo breed

Considering the escalating number of scientific reports on the association between the leptin gene and diverse physiological traits and performance of cattle populations, this study was directed towards identifying SNPs in the leptin gene among five indigenous cattle breeds of Ethiopia. DNA samples were extracted from the nasal swabs of the Ethiopian indigenous cattle breeds: Arsi (n = 18), Horro (n = 20), Begait (n = 21), Boran (n = 19), and Fogera (n = 17) and the Korean Hanwoo (a representative taurine breed) (n = 20), followed by PCR amplification of exon 2 and exon 3 regions of the leptin gene and sequence analysis of the PCR products. Five SNPs, two (generating missense mutations) on exon 2 and three (generating silent mutations) on exon 3 regions, were explicated in this study. Allele frequency and genotype frequency distribution pertaining to the SNPs were recorded for the studied cattle breeds besides the minor allele frequency and deviation from the Hardy-Weinberg equilibrium. Positive F_IS index values were recorded for all the markers except SNP2, illustrative of heterozygote deficiency. MEGA X software-based evolutionary divergence analysis of the phylogenetic tree based on the SNP data revealed that the large-sized breeds, Hanwoo, Begait, Boran, and Fogera, were more closely clustered compared to the small-sized Arsi breed. Among the seven haplotypes documented from the various breeds, sequence analysis was suggestive of haplotypes 1 and 2 to be ancestral haplotypes for the leptin gene. This study is envisaged to accelerate molecular breeding programs for the genetic improvement of the Ethiopian cattle breeds.

Polymorphisms in Fatty Acid Desaturase 2 Gene are Associated with Milk Production Traits in Chinese Holstein Cows

Simple Summary

Searching for causative polymorphisms underlying the variability of milk production traits and then incorporating them into breeding programs are very effective ways to improve the efficiency and reliability of conventional dairy cattle breeding. Fatty acid desaturase 2 (FADS2) plays a pivotal role in the biosynthesis of polyunsaturated fatty acids. Previous studies provided evidence that FADS2 was one of the most downregulated genes during negative energy balance in the liver of postpartum dairy cattle. Genes involved in the energetic pathways may influence other production traits, such as protein, fat and milk yields. Therefore, in the present study, we investigated the common genetic variants of the FADS2 gene in Chinese Holstein cows. Our results provided direct evidence that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis.

Abstract

This study investigated the single nucleotide polymorphisms (SNPs) of Fatty acid desaturase 2 (FADS2) gene and further explored their genetic effects on conventionally collected milk production traits in Chinese Holstein cows using 18,264 test-day records of 841 cows. One missense mutation c. 908 C > T (SNP site in the complementary DNA sequence), which caused an amino acid change from alanine to valine (294Ala > Val), and two 3’ untranslated region (UTR) SNPs, c.1571 G > A and c.2776 A > G were finally identified. The SNP c.908 C > T was significantly associated with test-day milk yield, fat percentage and 305-day milk, fat and protein yield. In particular, the T allele of the SNP c.908 C > T showed a significant association with decreased somatic cell score (SCS) in the investigated population. Significant relationship between the SNP c.1571 G > A and 305-day milk yield showed that genotype GG was linked to the highest milk yield. Substituting the allele G for A at the c.2776 A > G locus resulted in a decrease of protein percentage. Our results demonstrated that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis.

Phenotypic relationships, genetic parameters, genome-wide associations, and identification of potential candidate genes for ketosis and fat-to-protein ratio in German Holstein cows

Energy demand for milk production in early lactation exceeds energy intake, especially in high-yielding Holstein cows. Energy deficiency causes increasing susceptibility to metabolic disorders. In addition to several blood parameters, the fat-to-protein ratio (FPR) is suggested as an indicator for ketosis, because a FPR >1.5 refers to high lipolysis. The aim of this study was to analyze phenotypic, quantitative genetic, and genomic associations between FPR and ketosis. In this regard, 8,912 first-lactation Holstein cows were phenotyped for ketosis according to a veterinarian diagnosis key. Ketosis was diagnosed if the cow showed an abnormal carbohydrate metabolism with increased content of ketone bodies in the blood or urine. At least one entry for ketosis in the first 6 wk after calving implied a score = 1 (diseased); otherwise, a score = 0 (healthy) was assigned. The FPR from the first test-day was defined as a Gaussian distributed trait (FPRgauss), and also as a binary response trait (FPRbin), considering a threshold of FPR = 1.5. After imputation and quality controls, 45,613 SNP markers from the 8,912 genotyped cows were used for genomic studies. Phenotypically, an increasing ketosis incidence was associated with significantly higher FPR, and vice versa. Hence, from a practical trait recording perspective, first test-day FPR is suggested as an indicator for ketosis. The ketosis heritability was slightly larger when modeling the pedigree-based relationship matrix (pedigree-based: 0.17; SNP-based: 0.11). For FPRbin, heritabilities were larger when modeling the genomic relationship matrix (pedigree-based: 0.09; SNP-based: 0.15). For FPRgauss, heritabilities were almost identical for both pedigree and genomic relationship matrices (pedigree-based: 0.14; SNP-based: 0.15). Genetic correlations between ketosis with FPRbin and FPRgauss using either pedigree- or genomic-based relationship matrices were in a moderate range from 0.39 to 0.71. Applying genome-wide association studies, we identified the specific SNP rs109896020 (BTA 5, position: 115,456,438 bp) significantly contributing to ketosis. The identified potential candidate gene PARVB in close chromosomal distance is associated with nonalcoholic fatty liver disease in humans. The most important SNP contributing to FPRbin was located within the DGAT1 gene. Different SNP significantly contributed to ketosis and FPRbin, indicating different mechanisms for both traits genomically.

Management of Reproductive Disease in Dairy Cows

Postpartum diseases are common in dairy cows, and their incidence contributes to reduced fertility and increased risk of culling, making their prevention and management extremely important. Reproductive efficiency has a major impact on economic success of any dairy production unit. Optimizing reproductive efficiency contributes to overall efficiency of production units, minimizing environmental impacts and contributing to sustainability of food production. Additionally, control of reproductive diseases is important for maintenance of health and welfare of dairy cows; for minimizing use of antibiotics; and ensuring a wholesome, safe, and nutritious product.

Identification of genomic regions associated with female fertility in Danish Jersey using whole genome sequence data

BACKGROUND

Female fertility is an important trait in cattle breeding programs. In the Nordic countries selection is based on a fertility index (FTI). The fertility index is a weighted combination of four female fertility traits estimated breeding values for number of inseminations per conception (AIS), 56-day non-return rate (NRR), number of days from first to last insemination (IFL), and number of days between calving and first insemination (ICF). The objective of this study was to identify associations between sequence variants and fertility traits in Jersey cattle based on 1,225 Jersey sires from Denmark with official breeding values for female fertility traits. The association analyses were carried out in two steps: first the cattle genome was scanned for quantitative trait loci using a sire model for FTI using imputed whole genome sequence variants; second the significant quantitative trait locus regions were re-analyzed using a linear mixed model (animal model) for both FTI and its component traits AIS, NRR, IFL and ICF. The underlying traits were analyzed separately for heifers (first parity cows) and cows (later parity cows) for AIS, NRR, and IFL.

RESULTS

In the first step 6 QTL were detected for FTI: one QTL on each of BTA7, BTA20, BTA23, BTA25, and two QTL on BTA9 (QTL9-1 and QTL9-2). In the second step, ICF showed association with the QTL regions on BTA7, QTL9-2 QTL2 on BTA9, and BTA25, AIS for cows on BTA20 and BTA23, AIS for heifers on QTL9-2 on BTA9, IFL for cows on BTA20, BTA23 and BTA25, IFL for heifers on BTA7 and QTL9-2 on BTA9, NRR for heifers on BTA7 and BTA23, and NRR for cows on BTA23.

CONCLUSION

The genome wide association study presented here revealed 6 genomic regions associated with FTI. Screening these 6 QTL regions for the underlying female fertility traits revealed that different female fertility traits showed associations with different subsets of the individual FTI QTL peaks. The result of this study contributed to a better insight into the genetic control of FTI in the Danish Jersey.

Detection of DGAT1 gene polymorphism and its effect on selected biochemical indicators in dairy cows after calving

The aim of the study was to detect DGAT1 K232A polymorphism in 57 dairy cows of the Slovak Spotted breed and its crossbreds, and to assess possible effect of the given polymorphism on selected metabolic indices in blood serum after calving. Using the PCR-RFLP method with improved primers enabling better differentiation of genotypes we identified 45 homozygotes for alanine variant in this locus (AA genotype), 2 homozygotes for lysine variant (KK genotype), and 10 heterozygotes (AK genotype). Genotype frequencies were 0.790 for AA genotype, 0.175 for AK genotype, and only 0.035 for KK genotype. Allele frequencies were counted as 0.877 for A allele and 0.123 for K allele. In both groups of animals (AA and AK genotype) increased mean values above the upper reference limit of lactate dehydrogenase, and total bilirubin, and decreased levels below the lower reference limit of triglycerides were detected. In the group of animals with AA genotype we also noticed decreased levels of non-esterified fatty acids. On the other hand, increased serum concentrations of total immunoglobulins were found in animals with AK and KK genotype. This is the first study concerning DGAT1 polymorphism in the Slovak Spotted breed and its association with selected biochemical indicators.

Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle

Background

Identification of single nucleotide polymorphisms (SNPs) for specific genes involved in reproduction might improve reliability of genomic estimates for these low-heritability traits. Semen from 550 Holstein bulls of high (≥ 1.7; n = 288) or low (≤ −2; n = 262) daughter pregnancy rate (DPR) was genotyped for 434 candidate SNPs using the Sequenom MassARRAY® system. Three types of SNPs were evaluated: SNPs previously reported to be associated with reproductive traits or physically close to genetic markers for reproduction, SNPs in genes that are well known to be involved in reproductive processes, and SNPs in genes that are differentially expressed between physiological conditions in a variety of tissues associated in reproductive function. Eleven reproduction and production traits were analyzed.

Results

A total of 40 SNPs were associated (P < 0.05) with DPR. Among these were genes involved in the endocrine system, cell signaling, immune function and inhibition of apoptosis. A total of 10 genes were regulated by estradiol. In addition, 22 SNPs were associated with heifer conception rate, 33 with cow conception rate, 36 with productive life, 34 with net merit, 23 with milk yield, 19 with fat yield, 13 with fat percent, 19 with protein yield, 22 with protein percent, and 13 with somatic cell score. The allele substitution effect for SNPs associated with heifer conception rate, cow conception rate, productive life and net merit were in the same direction as for DPR. Allele substitution effects for several SNPs associated with production traits were in the opposite direction as DPR. Nonetheless, there were 29 SNPs associated with DPR that were not negatively associated with production traits.

Conclusion

SNPs in a total of 40 genes associated with DPR were identified as well as SNPs for other traits. It might be feasible to include these SNPs into genomic tests of reproduction and other traits. The genes associated with DPR are likely to be important for understanding the physiology of reproduction. Given the large number of SNPs associated with DPR that were not negatively associated with production traits, it should be possible to select for DPR without compromising production.

Associations between lipid metabolism and fertility in the dairy cow

Dairy cows mobilise body tissues to support milk production and, because glucose supplies are limited, lipids are used preferentially for energy production. Lipogenic activity is switched off and lipolytic mechanisms in adipose tissue increase through changes in the expression of several key enzymes. This results in a loss of body condition, together with high circulating concentrations of non-esterified fatty acids. Changes in the synthesis, secretion and signalling pathways of somatotrophic hormones (insulin, growth hormone, insulin-like growth factor 1) and adipokines (e.g. leptin) are central to the regulation of these processes. A high reliance on fatty acids as an energy source in the peripartum period causes oxidative damage to mitochondria in metabolically active tissues, including the liver and reproductive tract. The expression of genes involved in insulin resistance (PDK4, AHSG) is increased, together with expression of TIEG1, a transcription factor that can induce apoptosis via the mitochondrial pathway. Polymorphisms in TFAM and UCP2, two autosomal mitochondrial genes, have been associated with longevity in dairy cows. Polymorphisms in many other genes that affect lipid metabolism also show some associations with fertility traits. These include DGAT1, SCD1, DECR1, CRH, CBFA2T1, GH, LEP and NPY. Excess lipid accumulation in oocytes and the regenerating endometrium reduces fertility via reductions in embryo survival and increased inflammatory changes, respectively.

Genome-wide associations for feed utilisation complex in primiparous Holstein-Friesian dairy cows from experimental research herds in four European countries

Genome-wide association studies for difficult-to-measure traits are generally limited by the sample size with accurate phenotypic data. The objective of this study was to utilise data on primiparous Holstein–Friesian cows from experimental farms in Ireland, the United Kingdom, the Netherlands and Sweden to identify genomic regions associated with the feed utilisation complex: fat and protein corrected milk yield (FPCM), dry matter intake (DMI), body condition score (BCS) and live-weight (LW). Phenotypic data and 37 590 single nucleotide polymorphisms (SNPs) were available on up to 1629 animals. Genetic parameters of the traits were estimated using a linear animal model with pedigree information, and univariate genome-wide association analyses were undertaken using Bayesian stochastic search variable selection performed using Gibbs sampling. The variation in the phenotypes explained by the SNPs on each chromosome was related to the size of the chromosome and was relatively consistent for each trait with the possible exceptions of BTA4 for BCS, BTA7, BTA13, BTA14, BTA18 for LW and BTA27 for DMI. For LW, BCS, DMI and FPCM, 266, 178, 206 and 254 SNPs had a Bayes factor .3, respectively. Olfactory genes and genes involved in the sensory smell process were overrepresented in a 500 kbp window around the significant SNPs. Potential candidate genes were involved with functions linked to insulin, epidermal growth factor and tryptophan.

The relationship between fertility and lactation characteristics in Holstein cows on United Kingdom commercial dairy farms

Higher milk production is commonly associated with poorer fertility in dairy cows. This study used a biological model of lactation to define more closely which characteristics of lactation were linked to increased calving interval, an easily recorded measure of fertility in commercial dairy herds. Large data sets from a national milk recording scheme in the United Kingdom, collected over a 10-yr period, were used to calculate the genetic and phenotypic correlations between calving interval and a range of lactation traits within the Holstein breed. A lactation curve was fitted to each lactation, and several characteristics of each lactation were calculated. These were used in a series of mixed-model bivariate analyses with calving interval to derive the genetic parameters. When heifer lactation curve trait data were used, the highest genetic correlations were found with peak yield, maximum secretion potential, and total lactation milk yield (0.59±0.06 to 0.63±0.05), reflecting the observed phenomenon of poorer fertility and higher milk production. Genetic correlations for calving interval were also calculated with the rate of increase in milk yield in early lactation (0.46±0.08), persistency (0.36±0.09), day of peak yield (0.20±0.02), and relative cell death rate (-0.12±0.09). The daily production of milk components was highly genetically correlated with calving interval, with values for fat, protein, lactose, and water being 0.73±0.15, 0.48±0.13, 0.57±0.13, and 0.50±0.13. With these results and breeding values derived from these analyses, 2 possible strategies were suggested for improving the relationship between milk yield and fertility in dairy cows. First, animals that break the correlations described above could be selected (e.g., bulls with high peak yield and low calving intervals). Second, animals with lower peak yields but better persistency could be selected to maintain total milk yield.

Leptin in farm animals: where are we and where can we go?

Fat affects meat quality, value and production efficiency as well as providing energy reserves for pregnancy and lactation in farm livestock. Leptin, the adipocyte product of the obese (ob) gene, was quickly seen as a predictor of body fat content in animals approaching slaughter and an aid to assessing reproductive readiness in females. Its participation in inflammation and immune responses that help animals survive infection and trauma has clear additional relevance to meat and milk production. Furthermore, almost a decade of discoveries of nucleotide polymorphisms in the leptin and leptin receptor genes has suggested useful applications relating to feed intake regulation, the efficiency of feed use, the composition of growth, the timing of puberty, mammogenesis and mammary gland function and fertility in cattle, pigs and poultry. The current review attempts to summarise where research has taken us in each of these aspects and speculates on where future research might lead.

Whole-genome association study for energy balance and fat/protein ratio in German Holstein bull dams

The metabolic status of cows is important to health and fertility, especially in early lactation, and energy balance (EB) and fat/protein ratio (FPR) are considered as appropriate indicators for metabolic disorders. The aim of this study was to detect SNPs (single nucleotide polymorphisms) associated with EB and FPR in German Holstein bull dams belonging to the research herd Karkendamm. Bull dams were genotyped using the Illumina Bovine SNP 50K Bead chip(®) comprising 54 001 SNPs. A total of 43 593 SNPs and 586 (EB) and 668 (FPR) bull dams passed the quality control criteria. Phenotypes were deregressed breeding values estimated via random regression animal models for lactation days 11, 20, 30, and 42 for EB and FPR. Whole-genome association analyses were carried out fitting principal components as covariates to adjust for genetic substructure. Permutation tests were applied to estimate genome-wise significance. Across all observed lactation days, 19 SNPs located in four different intervals on chromosomes 1, 14, 16, and 27 were detected. For EB, seven markers across four chromosomes were identified. There was no overlap between markers associated with FPR and EB. SNPs associated with FPR were mostly located in QTL regions for milk production traits, especially in the region of DGAT1, whereas SNPs associated with EB mainly showed positional relationships to previously described QTL regions affecting functional traits, especially fertility.

Evaluating markers in selected genes for association with functional longevity of dairy cattle

Background

Longevity expressed as the number of days between birth and death is a trait of great importance for both human and animal populations. In our analysis we use dairy cattle to demonstrate how the association of Single Nucleotide Polymorphisms (SNPs) located within selected genes with longevity can be modeled. Such an approach can be extended to any genotyped population with time to endpoint information available. Our study is focused on selected genes in order to answer the question whether genes, known to be involved into the physiological determination of milk production, also influence individual's survival.

Results

Generally, the highest risk differences among animals with different genotypes are observed for polymorphisms located within the leptin gene. The polymorphism with a highest effect on functional longevity is LEP-R25C, for which the relative risk of culling for cows with genotype CC is 3.14 times higher than for the heterozygous animals. Apart from LEP-R25C, also FF homozygotes at the LEP-Y7F substitution attribute 3.64 times higher risk of culling than the YY homozygotes and VV homozygotes at LEP-A80V have 1.83 times higher risk of culling than AA homozygotes. Differences in risks between genotypes of polymorphisms within the other genes (the butyrophilin subfamily 1 member A1 gene, BTN1A1; the acyl-CoA:diacylglycerol acyltransferase 1 gene, DGAT1; the leptin receptor gene, LEPR; the ATP-binding cassette sub-family G member 2, ABCG2) are much smaller.

Conclusions

Our results indicate association between LEP and longevity and are very well supported by results of other studies related to dairy cattle. In view of the growing importance of functional traits in dairy cattle, LEP polymorphisms should be considered as markers supporting selection decisions. Furthermore, since the relationship between both LEP polymorphism and its protein product with longevity in humans is well documented, with our result we were able to demonstrate that livestock with its detailed records of family structure, genetic, and environmental factors as well as extensive trait recording can be a good model organism for research aspects related to humans.