Genome-wide association study identifies two major loci affecting calving ease and growth-related traits in cattle

Identifying quantitative trait loci (QTL) underlying complex, low-heritability traits is notoriously difficult. Prototypical for such traits, calving ease is an important breeding objective of cattle (Bos taurus)-improving programs. To identify QTL underlying calving ease, we performed a genome-wide association study using estimated breeding values (EBVs) as highly heritable phenotypes for paternal calving ease (pCE) and related traits. The massively structured study population consisted of 1800 bulls of the German Fleckvieh (FV) breed. Two pCE-associated regions on bovine chromosomes (BTA) 14 and 21 (P = 5.72 × 10(-15) and P = 2.27 × 10(-8), respectively) were identified using principal components analysis to correct for population stratification. The two most significantly associated SNPs explain 10% of the EBV variation. Since marker alleles with negative effect on pCE have positive effects on growth-related traits, the QTL may exert their effects on the birthing process through fetal growth traits. The QTL region on BTA14 corresponds to a human chromosome (HSA) region that is associated with growth characteristics. The HSA region corresponding to the BTA21 pCE QTL is maternally imprinted and involved in the Prader-Willi and Angelman syndromes. Resequencing of positional candidate genes on BTA14 revealed a highly significantly (P = 1.96 × 10(-14)) associated polymorphism ablating a polyadenylation signal of the gene encoding ribosomal protein S20 (RPS20). Our study demonstrates the leverage potential of EBVs in unraveling the genetic architecture of lowly heritable traits.     

Fatty acid synthase effects on bovine adipose fat and milk fat

A quantitative trait locus (QTL) was identified by linkage analysis on bovine Chromosome 19 that affects the fatty acid, myristic acid (C14:0), in subcutaneous adipose tissue of pasture-fed beef cattle (99% level: experiment-wise significance). The QTL was also shown to have significant effects on ten fatty acids in the milk fat of pasture-fed dairy cattle. A positional candidate gene for this QTL was identified as fatty acid synthase (FASN), which is a multifunctional enzyme with a central role in the metabolism of lipids. Five single nucleotide polymorphisms (SNPs) were identified in the bovine FASN gene, and animals were genotyped for FASN SNPs in three different cattle resource populations. Linkage and association mapping results using these SNPs were consistent with FASN being the gene underlying the QTL. SNP substitution effects for C14:0 percentage were found to have an effect in the opposite direction in adipose fat to that in milk fat. It is concluded that SNPs in the bovine FASN gene are associated with variation in the fatty acid composition of adipose fat and milk fat.     

Single nucleotide polymorphism and haplotype effects associated with somatic cell score in German Holstein cattle

Background

To better understand the genetic determination of udder health, we performed a genome-wide association study (GWAS) on a population of 2354 German Holstein bulls for which daughter yield deviations (DYD) for somatic cell score (SCS) were available. For this study, we used genetic information of 44 576 informative single nucleotide polymorphisms (SNPs) and 11 725 inferred haplotype blocks.

Results

When accounting for the sub-structure of the analyzed population, 16 SNPs and 10 haplotypes in six genomic regions were significant at the Bonferroni threshold of P ≤ 1.14 × 10^-6. The size of the identified regions ranged from 0.05 to 5.62 Mb. Genomic regions on chromosomes 5, 6, 18 and 19 coincided with known QTL affecting SCS, while additional genomic regions were found on chromosomes 13 and X. Of particular interest is the region on chromosome 6 between 85 and 88 Mb, where QTL for mastitis traits and significant SNPs for SCS in different Holstein populations coincide with our results. In all identified regions, except for the region on chromosome X, significant SNPs were present in significant haplotypes. The minor alleles of identified SNPs on chromosomes 18 and 19, and the major alleles of SNPs on chromosomes 6 and X were favorable for a lower SCS. Differences in somatic cell count (SCC) between alternative SNP alleles reached 14 000 cells/mL.

Conclusions

The results support the polygenic nature of the genetic determination of SCS, confirm the importance of previously reported QTL, and provide evidence for the segregation of additional QTL for SCS in Holstein cattle. The small size of the regions identified here will facilitate the search for causal genetic variations that affect gene functions.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle

Background

The focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle.

Methods

Data from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (−maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models.

Results and discussion

On average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (−maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (−maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence.

Conclusions

For the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (−maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.     

A genome wide association study for QTL affecting direct and maternal effects of stillbirth and dystocia in cattle

Dystocia and stillbirth are significant causes of female and neonatal death in many species and there is evidence for a genetic component to both traits. Identifying causal mutations affecting these traits through genome wide association studies could reveal the genetic pathways involved and will be a step towards targeted interventions. Norwegian Red cattle are an ideal model breed for such studies as very large numbers of records are available. We conducted a genome wide association study for direct and maternal effects of dystocia and stillbirth using almost 1 million records of these traits. Genotyping costs were minimized by genotyping the sires of the recorded cows, and using daughter averages as phenotypes. A dense marker map containing 17 343 single nucleotide polymorphisms covering all autosomal chromosomes was utilized. The genotyped sires were assigned to one of two groups in an attempt to ensure independence between the groups. Associations were only considered validated if they occurred in both groups. Strong associations were found and validated on chromosomes 4, 5, 6, 9, 12, 20, 22 and 28. The QTL region on chromosome 6 was refined using LDLA analysis. The results showed that this chromosome most probably contains two QTL for direct effect on dystocia and one for direct effect on stillbirth. Several candidate genes may be identified close to these QTL. Of these, a cluster of genes expected to affect bone and cartilage formation (i.e. SPP1, IBSP and MEPE) are of particular interest and we suggest that these genes are screened in candidate gene studies for dystocia and stillbirth in cattle as well as other species.     

Analyzes of genome-wide association follow-up study for calving traits in dairy cattle

ABSTRACT: BACKGROUND: There is often a pronounced disagreement between results obtained from different genome-wide association studies in cattle. There are multiple reasons for this disagreement. Particularly the presence of false positives leads to a need to validate detected QTL before they are optimally incorporated or weighted in selection decisions or further studied for causal gene. In dairy cattle progeny testing scheme new data is routinely accumulated which can be used to validate previously discovered associations. However, the data is not an independent sample and the sample size may not be sufficient to have enough power to validate previous discoveries. Here we compared two strategies to validate previously detected QTL when new data is added from the same study population. We compare analyzing a combined dataset (COMB) including all data presently available to only analyzing a validation dataset (VAL) i.e. a new dataset not previously analyzed as an independent replication. Secondly, we confirm SNP detected in the Reference population (REF) (i.e. previously analyzed dataset consists of older bulls) in the VAL dataset. RESULTS: Clearly the results from the combined (COMB) dataset which had nearly twice the sample size of other two subsets allowed the detection of far more significant associations than the two smaller subsets. The number of significant SNPs in REF (older bulls) was about four times higher compare to VAL (younger bulls) though both had similar sample sizes, 2,219 and 2,039 respectively. A total of 424 SNP-trait combinations on 22 chromosomes showed genome-wide significant association involving 284 unique SNPs in the COMB dataset. In the REF data set 101 associations (73 unique SNPs) and in the VAL 24 associations (18 unique SNPs) were found genome-wide significant.Sixty-eight percent of the SNPs in the REF dataset could be confirmed in the VAL dataset. Out of 469 unique SNPs showing chromosome-wide significant association with calving traits in the REF dataset 321 could be confirmed in the VAL dataset at P < 0.05. CONCLUSIONS: The follow-up study for GWAS in cattle will depend on the aim of the study. If the aim is to discover novel QTL, analyses of the COMB dataset is recommended, while in case of identification of the causal mutation underlying a QTL, confirmation of the discovered SNPs are necessary to avoid following a false positive.     

Identification of a Short Region on Chromosome 6 Affecting Direct Calving Ease in Piedmontese Cattle Breed

Calving in cattle is affected by calf morphology and by dam characteristics. It is described by two different traits: maternal calving ease, which is the ability to generate dams with good physiological predisposition to calving, and direct calving ease, which is the ability to generate calves that are easily born. The aim of this study was to identify regions of cattle genome harboring genes possibly affecting direct calving ease in the Piedmontese cattle breed. A population of 323 bulls scored for direct calving ease (EBV) was analyzed by a medium-density SNP marker panel (54,001 SNPs) to perform a genome-wide scan. The strongest signal was detected on chromosome 6 between 37.8 and 38.7 Mb where 13 SNPs associated to direct calving ease were found. Three genes are located in this region: LAP3, encoding for a leucine aminopeptidase involved in the oxytocin hydrolysis; NCAPG, encoding for a non-SMC condensin I complex, which has been associated in cattle with fetal growth and carcass size; and LCORL, which has been associated to height in humans and cattle. To further confirm the results of the genome-wide scan we genotyped additional SNPs within these genes and analyzed their association with direct calving ease. The results of this additional analysis fully confirmed the findings of the GWAS and particularly indicated LAP3 as the most probable gene involved. Linkage Disequilibrium (LD) analysis showed high correlation between SNPs located within LAP3 and LCORL indicating a possible selection signature due either to increased fitness or breeders’ selection for the trait.     

Investigation of the genetic architecture of a bone carcass weight QTL on BTA6

A previous analysis of an F₂/Backcross Charolais × Holstein cross population identified the presence of a highly significant QTL on chromosome 6 (BTA6) affecting the proportion of bone in the carcass. Two closely linked QTL affected birth weight (BW) and body length at birth (BBL). In this report, the marker density around the QTL on BTA6 was increased, adding four additional microsatellite markers across the chromosome and 46 SNPs within the target QTL confidence interval. Of the SNPs, 26 were in positional candidate genes and the remaining 20 provided an even distribution of markers in the target QTL region. As a bone‐related trait, the sum of the bone weight for all the left fore‐ and hindquarter joints of the carcass was analysed. We also studied the BW and BBL. Analyses of the data substantially reduced the QTL confidence interval. No strong evidence was found that the QTL for the three traits studied are different, and we conclude that the results are consistent with a single pleiotropic QTL influencing the three traits, with the largest effects on the proportion of bone in the carcass. The analyses also suggest that none of the SNPs tested is the sole causative variant of the QTL effects. Specifically, the SNP in the NCAPG gene previously reported as a causal mutation for foetal growth and carcass traits in other cattle populations was excluded as the causal mutation for the QTL reported here. Polymorphisms located in other previously identified candidate genes including SPP1, ABCG2, IBSP, MEPE and PPARGC1A were also excluded. The results suggest that SNP51_BTA‐119876 is the polymorphism in strongest linkage disequilibrium with the causal mutation(s). Further research is required to identify the causal variant(s) associated with this bone‐related QTL.     

A genome scan for quantitative trait loci influencing carcass,post‐natal growth and reproductive traits in commercial Angus cattle

To gain insight into the number of loci of large effect that underlie variation in cattle, a quantitative trait locus (QTL) scan for 14 economically important traits was performed in two commercial Angus populations using 390 microsatellites, 11 single nucleotide polymorphisms (SNPs) and one duplication loci. The first population comprised 1769 registered Angus bulls born between 1955 and 2003, with Expected Progeny Differences computed by the American Angus Association. The second comprised 38 half‐sib families containing 1622 steers with six post‐natal growth and carcass phenotypes. Linkage analysis was performed by half‐sib least squares regression with gridqtl or Bayesian Markov chain Monte Carlo analysis of complex pedigrees with loki . Of the 673 detected QTL, only 118 have previously been reported, reflecting both the conservative approach to QTL reporting in the literature, and the more liberal approach taken in this study. From 33 to 71% of the genetic variance and 35 to 56% of the phenotypic variance in each trait was explained by the detected QTL. To analyse the effects of 11 SNPs and one duplication locus within candidate genes on each trait, a single marker analysis was performed by fitting an additive allele substitution model in both mapping populations. There were 53 associations detected between the SNP/duplication loci and traits with ?log₁₀P_nominal≥ 4.0, where each association explained 0.92% to 4.4% of the genetic variance and 0.01% to 1.86% of the phenotypic variance. Of these associations, only six SNP/duplication loci were located within 8 cM of a QTL peak for the trait, with two being located at the QTL peak: SST_DG156121:c.362A>G for ribeye muscle area and TG_X05380:c.422C>T for calving ease. Strong associations between several SNP/duplication loci and trait variation were obtained in the absence of any detected linked QTL. However, we reject the causality of several commercialized DNA tests, including an association between TG_X05380:c.422C>T and marbling in Angus cattle.     

DNA polymorphisms of the prion doppel gene region in four different German cattle breeds and cows tested positive for bovine spongiform encephalopathy

Polymorphisms of the prion protein gene PRNP have been shown to influence the susceptibility/resistance to prion infections in human and sheep. In addition, the T174M polymorphism within the flanking prion doppel gene (PRND) was thought to be involved in susceptibility to sporadic Creutzfeldt-Jacob disease. To study a possible influence of DNA polymorphisms of the bovine PRND gene in bovine spongiform encephalopathy (BSE), previously identified and newly isolated DNA polymorphisms were genotyped in all available German cattle that tested positive for BSE. Genotypes and calculated haplotypes were compared with breeding bulls serving as controls. Analysis of the four major breeds Schwarzbunt (Holstein Friesian), Rotbunt (Holstein Red), Fleckvieh (Simmental), and Braunvieh (Swiss Brown) resulted in the isolation of the previously known polymorphisms R50H and R132Q and two novel synonymous single nucleotide polymorphisms (SNPs) C4820T and A5063T. Comparative genotype and haplotype analysis of BSE and control animals revealed a significantly different distribution of polymorphisms C4815T and R132Q in Fleckvieh animals but not in the other breeds tested. No association to BSE susceptibility was detectable for polymorphisms R50H and A5063T.     

Genetic Variants in ADD1 Gene and their Associations with Growth Traits in Cattle

The α-adducin (ADD1) is a subunit of adducin which is a cytoskeleton heterodimeric protein. Adducin participates in oocytes chromosome meiosis of mice, prompting adducin has an effect on embryonic development. Adducin gene mutation has significantly functional change. So the present study was to identify and characterize polymorphisms within the coding region of the bovine ADD1 gene among different cattle breeds. Here, 11 novel single nucleotide polymorphisms (SNPs 1–11) were identified by DNA sequencing and polymerase chain reaction-single stranded conformational polymorphism, there were one synonymous mutation in exon 1 (SNP1); four missense mutations in exons 4, 7, and 8 (SNPs 3–6); and six mutations in introns 4, 12, 13, and 14 (SNPs 2, 7–10). The statistical analyses indicated that the some SNPs are associated with the growth traits (body length, body height, chest circumference, and hucklebone width) in Chinese Jiaxian cattle population. Our results provide evidence that polymorphisms in the ADD1 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.     

Genome‐wide association study uncovers four QTL predisposing to supernumerary teats in cattle

Supernumerary teats (hyperthelia, SNTs) are a common abnormality of the bovine udder with a medium to high heritability and a postulated oligogenic or polygenic inheritance pattern. SNTs not only negatively affect machine milking ability but also act as a reservoir for bacteria. A genome‐wide association study was carried out to identify genes involved in the development of SNTs in the dual‐purpose Fleckvieh breed. A total of 2467 progeny‐tested bulls were genotyped at 43 698 single nucleotide polymorphisms, and daughter yield deviations (DYDs) for ‘udder clearness' (UC) were used as high‐heritability phenotypes. Massive structuring of the study population was accounted for by principal components analysis‐based and mixed model‐based approaches. Four loci on BTA5, BTA6, BTA11 and BTA17 were significantly associated with the UC DYD. Three associated regions contain genes of the highly conserved Wnt signalling pathway. The four QTL together account for 10.7% of the variance of the UC DYD, whereas the major fraction of the DYD variance is attributable to chromosomes with no identified QTL. Our results support both an oligogenic and a polygenic inheritance pattern of SNTs in cattle. The identified candidate genes permit insights into the genetic architecture of teat malformations in cattle and provide clues to unravel the molecular mechanisms of mammary gland alterations in cattle and other species.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

Detection of genes influencing economic traits in three French dairy cattle breeds

A project of QTL detection was carried out in the French Holstein, Normande, and Montbéliarde dairy cattle breeds. This granddaughter design included 1 548 artificial insemination bulls distributed in 14 sire families and evaluated after a progeny-test for 24 traits (production, milk composition, persistency, type, fertility, mastitis resistance, and milking ease). These bulls were also genotyped for 169 genetic markers, mostly microsatellites. The QTL were analysed by within-sire linear regression of daughter yield deviations or deregressed proofs on the probability that the son receives one or the other paternal QTL allele, given the marker information. QTL were detected for all traits, including those with a low heritability. One hundred and twenty QTL with a chromosome-wise significance lower than 3% were tabulated. This threshold corresponded to a 15% false discovery rate. Amongst them, 32 were genome-wise significant. Estimates of their contribution to genetic variance ranged from 6 to 40%. Most substitution effects ranged from 0.6 to 1.0 genetic standard deviation. For a given QTL, only 1 to 5 families out of 14 were informative. The confidence intervals of the QTL locations were large and always greater than 20 cM. This experiment confirmed several already published QTL but most of them were original, particularly for non-production traits.     

Validation of whole genome linkage‐linkage disequilibrium and association results,and identification of markers to predict genetic merit for twinning

A previous genome‐wide search with a moderate density 10K marker set identified many marker associations with twinning rate, either through single‐marker analysis or combined linkage‐linkage disequilibrium (LLD; haplotype) analysis. The objective of the current study was to validate putative marker associations using an independent set of phenotypic data. Holstein bulls (n = 921) from 100 paternal half‐sib families were genotyped. Twinning rate predicted transmitting abilities were calculated using calving records from 1994 to 1998 (Data I) and 1999 to 2006 (Data II), and the underlying liability scores from threshold model analysis were used as the trait in marker association analyses. The previous analysis used 201 bulls with daughter records in Data I. In the current analysis, this was increased to 434, providing a revised estimate of effect and significance. Bulls with daughter records in Data II totaled 851, and analysis of this data provided the validation of results from analysis of Data I. Single nucleotide polymorphisms (SNPs) were selected to validate previously significant single‐marker associations and LLD results. Bulls were genotyped for a total of 306 markers. Nine of 13 LLD regions located on chromosomes 1, 2, 3, 6, 9, 22, 23(2) and 26 were validated, showing significant results for both Data I and II. Association analysis revealed 55 of 174 markers validated, equating to a single‐marker validation rate of 31%. Stepwise backward elimination and cross‐validation analyses identified 18 SNPs for use in a final reduced marker panel explaining 34% of the genetic variation, and to allow prediction of genetic merit for twinning rate.     

Genome Wide Analysis of Fertility and Production Traits in Italian Holstein Cattle

A genome wide scan was performed on a total of 2093 Italian Holstein proven bulls genotyped with 50K single nucleotide polymorphisms (SNPs), with the objective of identifying loci associated with fertility related traits and to test their effects on milk production traits. The analysis was carried out using estimated breeding values for the aggregate fertility index and for each trait contributing to the index: angularity, calving interval, non-return rate at 56 days, days to first service, and 305 day first parity lactation. In addition, two production traits not included in the aggregate fertility index were analysed: fat yield and protein yield. Analyses were carried out using all SNPs treated separately, further the most significant marker on BTA14 associated to milk quality located in the DGAT1 region was treated as fixed effect. Genome wide association analysis identified 61 significant SNPs and 75 significant marker-trait associations. Eight additional SNP associations were detected when SNP located near DGAT1 was included as a fixed effect. As there were no obvious common SNPs between the traits analyzed independently in this study, a network analysis was carried out to identify unforeseen relationships that may link production and fertility traits.     

Identification of single nucleotide polymorphisms in the bovine CCL2, IL8, CCR2 and IL8RA genes and their association with health and production in Canadian Holsteins

The aim of this study was to identify the presence of SNPs in the chemokine genes CCL2 and IL8 and the chemokine receptor genes IL8RA and CCR2, and assess their potential contribution to variation in estimated breeding values (EBVs) for somatic cell score (SCS) and four other traits in Canadian Holstein bulls. Pools of DNA for bulls with high (H) and low (L) EBVs for SCS were used for identification of 11 SNPs. Two unreported SNPs were found in the CCL2 gene and one SNP was found in the CCR2 gene. Previously reported SNPs (three in the IL8 gene and five in the IL8RA chemokine receptor) were also identified. Two SNPs in CCL2, three in IL8, one in IL8RA and one in CCR2 were genotyped in Canadian Holstein bulls (n = 338) using tetra primer ARMS-PCR. We investigated associations of these seven polymorphisms with three production traits (milk yield, fat yield and protein yield) and one conformation trait related to mastitis (udder depth). The allele substitution effect for the CCL2 rs41255713:T>C SNP was significant at an experimental-wise level for milk yield (247.5 +/- 79.9 kg) and protein yield (7.4 +/- 2.3 kg) EBVs (P T SNP on SCS was significant at the comparison-wise level (-0.04 +/- 0.02, P = 0.05), which might indicate a possible association in support of other published studies. Lastly, we assigned CCR2 to BTA22q24, where a previously QTL for SCS was identified.