A genome-wide association study for quantitative trait loci of show-jumping in Hanoverian warmblood horses

Multi-omics reveals changed energy metabolism of liver and muscle by caffeine after mice swimming

In recent years, numerous studies have investigated the effects of caffeine on exercise, and provide convincing evidence for its ergogenic effects on exercise performance. However, the precise mechanisms underlying these ergogenic effects remain unclear. In this study, an exercise swimming model was conducted to investigate the effects of orally administered with caffeine before swimming on the alterations of proteome and energy metabolome of liver and muscle after swimming. We found proteins in liver, such as S100a8, S100a9, Gabpa, Igfbp1 and Sdc4, were significantly up-regulated, while Rbp4 and Tf decreased after swimming were further down-regulated in caffeine group. The glycolysis and pentose phosphate pathways in liver and muscle were both significantly down-regulated in caffeine group. The pyruvate carboxylase and amino acid levels in liver, including cysteine, serine and tyrosine, were markedly up-regulated in caffeine group, exhibiting a strong correlation with the increased pyruvic acid and oxaloacetate levels in muscle. Moreover, caffeine significantly decreased the lactate levels in both liver and muscle after swimming, potentially benefiting exercise performance.

Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach

BACKGROUND

The Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds.

RESULTS

The k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne).

CONCLUSIONS

Distinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.

Population structure and genomic footprints of selection in five major Iranian horse breeds

The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log₁₀ rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log₁₀ rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.

Genomic Correlations Between the Gaits of Young Horses Measured by Accelerometry and Functional Longevity in Jumping Competition

Functional longevity is essential for the well-being of horses and the satisfaction of riders. Conventional selection using longevity breeding values calculated from competition results is not efficient because it takes too long to obtain reliable information. Therefore, the objective was to identify early criteria for selection. We assessed two types of early criteria: gait traits of young horses and QTLs. Thus, our aim was to estimate the genetic correlation between gait traits and longevity and to perform a genome-wide association study (GWAS) for longevity. Measurements of gaits by accelerometry were recorded on 1,477 show jumping horses that were 4 to 5 years old. Gait analysis provided 9 principal components describing trot, canter, and walk. Longevity estimated breeding values (EBVs) for stallions were calculated using a survival analysis of more than 900,000 years of performances by 179,448 show jumping horses born from 1981 onwards. Longevity was measured as the number of years spent in competition. Model included region and month of birth, age at first competition, year, and performance level. Longevity EBVs were deregressed to obtain weighted pseudo-performances for 1,968 stallions. Genomic data were available for 3,658 jumping horses. Seventy-eight percent of the horses measured for gaits and twenty-five percent of those measured for longevity were genotyped. A GWAS of longevity revealed no significant QTLs. Genetic parameters between each of the 9 principal components of the gait variables and longevity were evaluated with a bi-trait animal linear mixed model using single-step GBLUP analysis with the relationship matrix constructed from genomic data and genealogy (24,448 ancestors over four generations). The heritability of the gait traits varied from 0.11 to 0.44. The third principal component for trot (high lateral activity) and the first principal component for canter (high dorsoventral activity and low stride frequency) were moderately genetically correlated with higher longevity: r_g = 0.38 (0.15) and 0.28 (0.13), respectively. Our study revealed that functional longevity is a polygenic trait with no major genes. We found new correlations between longevity and gait traits. Before using gait characteristics in a selection plan, these correlations need to be understood better at the biomechanical level.

Detection of Selection Signatures Underlying Production and Adaptive Traits Based on Whole-Genome Sequencing of Six Donkey Populations

Simple Summary

After a long period of artificial selection, the donkey now presents a variety of body types and production performance values. In this experiment, we performed selective signal scanning on the second-generation resequencing data of six different varieties. The regions and candidate genes related to artificial selection were identified to provide reference for future breeding.

Abstract

Donkeys (Equus asinus) are an important farm animal. After long-term natural and artificial selection, donkeys now exhibit a variety of body sizes and production performance values. In this study, six donkey breeds, representing different regions and phenotypes, were used for second-generation resequencing. The sequencing results revealed more than seven million single nucleotide variants (SNVs), with an average of more than four million SNVs per species. We combined two methods, Z-transformed heterozygosity (ZHp) and unbiased estimates of pairwise fixation index (di) values, to analyze the signatures of selection. We mapped 11 selected regions and identified genes associated with coat color, body size, motion capacity, and high-altitude adaptation. These candidate genes included staining (ASIP and KITLG), body type (ACSL4, BCOR, CDKL5, LCOR, NCAPG, and TBX3), exercise (GABPA), and adaptation to low-oxygen environments (GLDC and HBB). We also analyzed the SNVs of the breed-specific genes for their potential functions and found that there are three varieties in the conserved regions with breed-specific mutation sites. Our results provide data to support the establishment of the donkey SNV chip and reference information for the utilization of the genetic resources of Chinese domestic donkeys.

Genomic Divergence in Swedish Warmblood Horses Selected for Equestrian Disciplines

The equestrian sport horse Swedish Warmblood (SWB) originates from versatile cavalry horses. Most modern SWB breeders have specialized their breeding either towards show jumping or dressage disciplines. The aim of this study was to explore the genomic structure of SWB horses to evaluate the presence of genomic subpopulations, and to search for signatures of selection in subgroups of SWB with high or low breeding values (EBVs) for show jumping. We analyzed high density genotype information from 380 SWB horses born in the period 2010-2011, and used Principal Coordinates Analysis and Discriminant Analysis of Principal Components to detect population stratification. Fixation index and Cross Population Extended Haplotype Homozygosity scores were used to scan the genome for potential signatures of selection. In accordance with current breeding practice, this study highlights the development of two separate breed subpopulations with putative signatures of selection in eleven chromosomes. These regions involve genes with known function in, e.g., mentality, endogenous reward system, development of connective tissues and muscles, motor control, body growth and development. This study shows genetic divergence, due to specialization towards different disciplines in SWB horses. This latter evidence can be of interest for SWB and other horse studbooks encountering specialized breeding.

Genome-Wide Association Analyses of Equine Metabolic Syndrome Phenotypes in Welsh Ponies and Morgan Horses

Equine metabolic syndrome (EMS) is a complex trait for which few genetic studies have been published. Our study objectives were to perform within breed genome-wide association analyses (GWA) to identify associated loci in two high-risk breeds, coupled with meta-analysis to identify shared and unique loci between breeds. GWA for 12 EMS traits identified 303 and 142 associated genomic regions in 264 Welsh ponies and 286 Morgan horses, respectively. Meta-analysis demonstrated that 65 GWA regions were shared across breeds. Region boundaries were defined based on a fixed-size or the breakdown of linkage disequilibrium, and prioritized if they were: shared between breeds or across traits (high priority), identified in a single GWA cohort (medium priority), or shared across traits with no SNPs reaching genome-wide significance (low priority), resulting in 56 high, 26 medium, and seven low priority regions including 1853 candidate genes in the Welsh ponies; and 39 high, eight medium, and nine low priority regions including 1167 candidate genes in the Morgans. The prioritized regions contained protein-coding genes which were functionally enriched for pathways associated with inflammation, glucose metabolism, or lipid metabolism. These data demonstrate that EMS is a polygenic trait with breed-specific risk alleles as well as those shared across breeds.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Comparative population genomics unveils candidate genes for athletic performance in Hanoverians

Equine athletes have a genetic heritage that has been evolved for millions of years, which provides an opportunity to study the genetics of locomotion pattern and performance in mammals. The Hanoverian, a breed originating in Germany, is arguably among the most athletic of horse breeds, as well as possessing a balanced character and beautiful appearance. Here, we compared the whole genomes of Hanoverian with three other horse breeds (Akhal-Teke, Franches-Montagnes, and Standardbred), using the fixation index (Fst) and cross-population composite likelihood ratio (XP-CLR) methods for testing the multi-locus allele frequency differentiation between populations. We identified 299 and 485 positively selected genes using the Fst and XP-CLR methods, respectively. Further functional analyses showed that the ACTA1 gene is potentially involved in athletic performance in the Hanoverian breed, consistent with its role observed in human population. In addition, three other loci on chromosomes 1 and 20 were identified to be potentially involved in equine physical performance. The selected candidate genes identified in this study may be useful in current breeding efforts to develop improved breeds in regard to athletic performance.

A genome-wide scan for diversifying selection signatures in selected horse breeds

The genetic differentiation of the current horse population was evolutionarily created by natural or artificial selection which shaped the genomes of individual breeds in several unique ways. The availability of high throughput genotyping methods created the opportunity to study this genetic variation on a genome-wide level allowing detection of genome regions divergently selected between separate breeds as well as among different horse types sharing similar phenotypic features. In this study, we used the population differentiation index (FST) that is generally used for measuring locus-specific allele frequencies variation between populations, to detect selection signatures among six horse breeds maintained in Poland. These breeds can be classified into three major categories, including light, draft and primitive horses, selected mainly in terms of type (utility), exterior, performance, size, coat color and appearance. The analysis of the most pronounced selection signals found in this study allowed us to detect several genomic regions and genes connected with processes potentially important for breed phenotypic differentiation and associated with energy homeostasis during physical effort, heart functioning, fertility, disease resistance and motor coordination. Our results also confirmed previously described association of loci on ECA3 (spanning LCORL and NCAPG genes) and ECA11 (spanning LASP1 gene) with the regulation of body size in our draft and primitive (small size) horses. The efficiency of the applied FST-based approach was also confirmed by the identification of a robust selection signal in the blue dun colored Polish Konik horses at the locus of TBX3 gene, which was previously shown to be responsible for dun coat color dilution in other horse breeds. FST-based method showed to be efficient in detection of diversifying selection signatures in the analyzed horse breeds. Especially pronounced signals were observed at the loci responsible for fixed breed-specific features. Several candidate genes under selection were proposed in this study for traits selected in separate breeds and horse types, however, further functional and comparative studies are needed to confirm and explain their effect on the observed genetic diversity of the horse breeds.

A genome-wide association study for harness racing success in the Norwegian-Swedish coldblooded trotter reveals genes for learning and energy metabolism

Background

Although harness racing is of high economic importance to the global equine industry, significant genomic resources have yet to be applied to mapping harness racing success. To identify genomic regions associated with harness racing success, the current study performs genome-wide association analyses with three racing performance traits in the Norwegian-Swedish Coldblooded Trotter using the 670 K Axiom Equine Genotyping Array.

Results

Following quality control, 613 horses and 359,635 SNPs were retained for further analysis. After strict Bonferroni correction, nine genome-wide significant SNPs were identified for career earnings. No genome-wide significant SNPs were identified for number of gallops or best km time. However, four suggestive genome-wide significant SNPs were identified for number of gallops, while 19 were identified for best km time. Multiple genes related to intelligence, energy metabolism, and immune function were identified as potential candidate genes for harness racing success.

Conclusions

Apart from the physiological requirements needed for a harness racing horse to be successful, the results of the current study also advocate learning ability and memory as important elements for harness racing success. Further exploration into the mental capacity required for a horse to achieve racing success is likely warranted.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0670-3) contains supplementary material, which is available to authorized users.

Foundations of performance – factors that contribute to excellence in equine exercise

Horses are renowned for their incredible capacity for a range of athletic activities, and participation in athletic events arguably represents the most critical strut of the equine industry. Successful performance is typically a primary focus during participation in competitive athletic events, and relies upon a variety of innate physiological and structural factors of the athlete. However, a wide range of external factors also influence performance, and many of these can be readily manipulated. Therefore, thorough assessment of the individual’s inherent capacity for a specific athletic discipline must be combined with optimisation of external factors including nutrition and training to promote excellent performance. Recent progress in methods of athlete selection and monitoring of training responses are assisting continued improvements in equine performance.

Quantitative analysis of short- and long-distance racing performance in young and adult horses and association analysis with functional candidate genes in Spanish Trotter horses

The association of five candidate genes with sporting performance in young and adult Spanish Trotter horses (STHs) was performed according to a previous selection based on quantitative analysis of the trait time per kilometre (TPK). A total of 334 516 records of TPK from 5958 STHs were used to estimate the estimated breeding values (EBVs) at different age groups (young and adults horses) throughout the range of distances (1600-2700 m) using a bicharacter random regression model. The heritability estimated by distance ranged from 0.16 to 0.40, with a different range for the two age groups. Considering the animals with the best and the worst deregressed EBV, 321 STHs were selected for SNP genotyping in MSTN, COX4I2, PDK4, DMRT3 and CKM genes. An association analysis based on ridge and logistic regression revealed that the young trotters with genotype GG in PDK4 (p < 0.05) and AA of DMRT3 (p < 0.001) SNPs show the best potential in short-distance races, while those carrying the genotype AA in DMRT3 (p < 0.001) and CC in CKM (p < 0.05) genes seem to be the best in long-distance races. Adult trotters with genotype AA in DMRT3 also display greater speed (p < 0.05) and endurance (p < 0.001).

Genome-wide association study for jumping performances in French sport horses

A genome-wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo-Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54,602 SNPs available on this chip, 44,424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single-marker mixed model and a haplotype-based mixed model (significant: P < 1E-05; suggestive: P < 1E-04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.

Genome-wide association study for wool production traits in a Chinese Merino sheep population

Genome-wide association studies (GWAS) provide a powerful approach for identifying quantitative trait loci without prior knowledge of location or function. To identify loci associated with wool production traits, we performed a genome-wide association study on a total of 765 Chinese Merino sheep (JunKen type) genotyped with 50 K single nucleotide polymorphisms (SNPs). In the present study, five wool production traits were examined: fiber diameter, fiber diameter coefficient of variation, fineness dispersion, staple length and crimp. We detected 28 genome-wide significant SNPs for fiber diameter, fiber diameter coefficient of variation, fineness dispersion, and crimp trait in the Chinese Merino sheep. About 43% of the significant SNP markers were located within known or predicted genes, including YWHAZ, KRTCAP3, TSPEAR, PIK3R4, KIF16B, PTPN3, GPRC5A, DDX47, TCF9, TPTE2, EPHA5 and NBEA genes. Our results not only confirm the results of previous reports, but also provide a suite of novel SNP markers and candidate genes associated with wool traits. Our findings will be useful for exploring the genetic control of wool traits in sheep.

Genome-wide analysis of DNA methylation patterns in horse

Background

DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known.

Results

We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density.

Conclusions

We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-598) contains supplementary material, which is available to authorized users.

Applied equine genetics

Genome sequencing of the domestic horse and subsequent advancements in the field of equine genomics have led to an explosion in the development of tools for mapping traits and diseases and evaluating gene expression. The objective of this review is to discuss the current progress in the field of equine genomics, with specific emphasis on assembly and analysis of the reference sequence and subsequent sequencing of a Quarter Horse mare; the genomic tools currently available to researchers and their implications in genomic investigations in the horse; the genomics of Mendelian and non-Mendelian traits; the genomics of performance traits and considerations regarding genetic testing in the horse. The whole-genome sequencing of a Quarter Horse mare has provided additional variants within the equine genome that extend past single nucleotide polymorphisms to include insertions/deletions and copy number variants. Equine single nucleotide polymorphism arrays have allowed for the investigation of both simple and complex genetic traits while DNA microarrays have provided a tool for examining gene expression across various tissues and with certain disease conditions. Recently, next-generation sequencing has become more affordable and both whole-genome DNA sequencing and transcriptome-wide RNA sequencing are methodologies that are being applied to equine genomic research. Research in the field of equine genomics continues to expand rapidly as the cost of genotyping and sequencing decreases, resulting in a need for quality bioinformatics software and expertise to appropriately handle both the size and complexity of these data.

Genomic selection: Status in different species and challenges for breeding

Technical advances and development in the market for genomic tools have facilitated access to whole-genome data across species. Building-up on the acquired knowledge of the genome sequences, large-scale genotyping has been optimized for broad use, so genotype information can be routinely used to predict genetic merit. Genomic selection (GS) refers to the use of aggregates of estimated marker effects as predictors which allow improved individual differentiation at young age. Realizable benefits of GS are influenced by several factors and vary in quantity and quality between species. General characteristics and challenges of GS in implementation and routine application are described, followed by an overview over the current status of its use, prospects and challenges in important animal species. Genetic gain for a particular trait can be enhanced by shortening of the generation interval, increased selection accuracy and increased selection intensity, with species- and breed-specific relevance of the determinants. Reliable predictions based on genetic marker effects require assembly of a reference for linking of phenotype and genotype data to allow estimation and regular re-estimation. Experiences from dairy breeding have shown that international collaboration can set the course for fast and successful implementation of innovative selection tools, so genomics may significantly impact the structures of future breeding and breeding programmes. Traits of great and increasing importance, which were difficult to improve in the conventional systems, could be emphasized, if continuous availability of high-quality phenotype data can be assured. Equally elaborate strategies for genotyping and phenotyping will allow tailored approaches to balance efficient animal production, sustainability, animal health and welfare in future.

A post-GWAS replication study confirming the PTK2 gene associated with milk production traits in Chinese Holstein

Our initial genome-wide association study (GWAS) demonstrated that two SNPs (ARS-BFGL-NGS-33248, UA-IFASA-9288) within the protein tyrosine kinase 2 (PTK2) gene were significantly associated with milk production traits in Chinese Holstein dairy cattle. To further validate if the statistical evidence provided in GWAS were true-positive findings, a replication study was performed herein through genotype-phenotype associations. The two tested SNPs were found to show significant associations with milk production traits, which confirmed the associations observed in the original study. Specifically, SNPs lying in the PTK2 gene were also detected by sequencing 14 unrelated sires in Chinese Holsteins and a total of thirty-three novel SNPs were identified. Thirteen out of these identified SNPs were genotyped and tested for association with milk production traits in an independent resource population. After Bonferroni correction for multiple testing, twelve SNPs were statistically significant for more than two milk production traits. Analyses of pairwise D' measures of linkage disequilibrium (LD) between all SNPs were also explored. Two haplotype blocks were inferred and the association study at haplotype level revealed similar effects on milk production traits. In addition, the RNA expression analyses revealed that a non-synonymous coding SNP (g.4061098T>G) was involved in the regulation of gene expression. Thus the findings presented here provide strong evidence for associations of PTK2 variants with dairy production traits and may be applied in Chinese Holstein breeding program.