Fine mapping of milk production QTL on BTA6 by combined linkage and linkage disequilibrium analysis

Combined linkage and linkage disequilibrium analysis were used to refine the position of a previously detected QTL affecting milk production traits on bovine chromosome 6. Through a series of single- and multitrait and single- and multipoint QTL analyses, the QTL could be positioned to a 7.5-cM interval surrounded by the markers BMS2508 and FBN12. The most significant results were found for fat percentage and protein percentage. This effect seemed to be caused by a QTL allele embedded in one specific marker haplotype that caused a reduction in fat and protein yields and a concomitant increase of milk yield, thus resulting in a marked reduction of fat and protein percentages.     

Characterization of the DGAT1 mutations and the CSN1S1 promoter in the German Angeln dairy cattle population

The identification of quantitative trait loci (QTL) and genes with influence on milk production traits has been the objective of various mapping studies in the last decade. In the centromeric region of Bos taurus autosome (BTA) 14, the acyl-CoA:diacylglycerol acyltransferase1 gene (DGAT1) has been identified as the most likely causative gene underlying a QTL for milk fat yield and content. Recently, a second polymorphism in the promoter of DGAT1 emerged as an additional source of variation. In this study, the frequencies and the effects of alleles at the DGAT1 K232A and at the DGAT1 promoter variable number of tandem repeat (VNTR) locus on BTA14, and of alleles at the CSN1S1 (α_S1-casein-encoding gene) promoter on BTA6 in the German Angeln dairy cattle population were investigated. Analyzed traits were milk, fat, protein, lactose, and milk energy yield, fat, protein, lactose, and milk energy content and somatic cell score. The lysine variant of the DGAT1 K232A mutation showed significant effects for most of the milk production traits. A specific allele of the DGAT1 promoter VNTR showed significant effects on the traits lactose yield and content, milk energy content, and SCS compared with the other alleles. Additionally, a regulation mechanism between the DGAT1 K232A mutation and the DGAT1 promoter VNTR was found for fat yield and content, which could be caused by an upper physiological bound for the effects of the DGAT1 gene. At the CSN1S1 promoter, 2 of 4 alleles showed significant allele substitution effects on the milk yield traits.     

Caprine CSN1S1 haplotype effect on gene expression and milk composition measured by Fourier transform infrared spectroscopy

The Norwegian dairy goat population has a high frequency of a CSN1S1 (α_S1-casein) haplotype with negative effects on protein and fat content. It is characterized by a single point deletion in exon 12 of CSN1S1, leading to a truncated protein and hence a low content of α_S1-casein in the milk. This haplotype together with another haplotype with a deletion in exon 9 are called “weak” haplotypes. “Strong” haplotypes, on the other hand, have positive effects on important milk production traits. We show that expression of CSN1S1 in the mammary gland of lactating goats is significantly lower in animals with 2 weak haplotypes. Moreover, the effects of defective alleles were not detected in animals having 1 strong and 1 weak haplotype. Expression levels of other genes in the casein cluster were not affected by the CSN1S1 haplotypes investigated. Milk samples from goats with 2 weak haplotypes could be distinguished from the other milk samples using Fourier transform infrared (FTIR) spectroscopy and partial least squares discriminant analysis (PLS-DA). The PLS-DA components were related to spectra of pure caseins and whey proteins, hence FTIR has a potential for identifying milk samples with low α_S1-casein content and different protein composition. The results indicate that FTIR-based measurements can be incorporated into breeding plans, or for selection of milk samples with high casein content, which in turn may improve cheese-making properties of the milk.     

Effects of casein haplotypes on milk production traits in Italian Holstein and Brown Swiss cattle

The objective of this study was to estimate the effects of different haplotypes of the casein genes on milk production traits in Italian dairy cattle. Traits of interest were yields of milk, fat, and protein, and percentages of fat and protein in milk. The data included 728 multiparous records from 347 Holsteins and 773 records from 298 Brown Swiss cows. Records were preadjusted for effects of age and parity, season of calving, and region, and expressed as deviations from herdmate averages. Twenty half-sib families were represented in each breed. Haplotype probabilities were estimated for each animal and phenotypes were regressed on these probabilities. Nine haplotypes were observed in Holsteins and 17 were identified among the Brown Swiss. For Holsteins, significant effects were observed for protein percentage, with some indication of an effect for fat percentage. For the Brown Swiss, effects of haplotypes were significant for milk yield and fat and protein percentages. Effects were strongest for protein percentage. Correlation coefficients of solutions across breeds tended to be strong and positive, indicating that the same haplotypes had similar estimated effects in the 2 breeds. Although the data were limited (<350 cows in each study), this latter result may suggest that genes in the casein complex itself are responsible for the effects observed, rather than loci that are physically linked on either side of the casein cluster.     

Statistical modeling of candidate gene effects on milk production traits in dairy cattle

A major objective of dairy cattle genomic research is to identify genes underlying the variability of milk production traits that could be useful in breeding programs. The candidate gene approach provides tools for searching for causative polymorphisms affecting quantitative traits. Genes with a possible effect on milk traits in cattle can be involved in different physiological pathways, such as triglyceride synthesis [acyl-CoA:diacylglycerol acyltransferase 1 gene (DGAT1)], fat secretion from the mammary epithelial tissue (butyrophilin), or entire-body energy homeostasis regulation (leptin and leptin receptor). In this study, based on data from 252 Black and White bulls from the active Polish dairy population, effects and potential interactions of 9 single nucleotide polymorphisms in the butyrophilin, DGAT1, leptin, and leptin receptor genes were investigated. Additionally, the effect of the number of additive, dominance, and epistatic genetic effects fitted into the model on the estimates of model parameters and model selection was illustrated. Phenotypic records were daughter yield deviations for milk, fat, and protein yields, obtained from a routine national genetic evaluation. Out of all the analyzed polymorphisms, DGAT1 K232A had a much larger effect on milk traits than the other single nucleotide polymorphisms considered. Estimates of the additive genetic effect of K232A expressed as half of the difference between Lys- and Ala-encoding variants were −107.4 kg of milk, 5.4 kg of fat, and −1.6 kg of protein at first parity, as well as −120 kg of milk and 6.8 kg of fat at second parity. In terms of model selection, it was demonstrated that the modified version of Bayesian information criterion selects models with the parameterization reflecting the genetic background of the analyzed trait, while the Bayesian information criterion chooses models that are too highly parameterized.     

Testing candidate gene effects on milk production traits in dairy cattle under various parameterizations and modes of inheritance

The major objectives of this study were 1) to assess the statistical properties of models commonly used for the estimation of single nucleotide polymorphism (SNP) effects under the assumption of various modes of inheritance and various parameterizations of SNP genotypes using simulated data, and 2) to compare effects of the selected polymorphisms located within butyrophilin (BTN1A1), diacylglycerol acyltransferase 1 (DGAT1), leptin (LEP), and leptin receptor (LEPR) candidate genes on milk production traits using data from 2 dairy cattle breeds (190 Jersey cows and 475 Polish Holstein-Friesian cows). Simulation results showed that type I error and power were not dependent on the assumed parameterization, but differences were observed regarding confidence intervals of estimated SNP effects. In the presence of epistasis, correct confidence intervals for all (epistatic and nonepistatic) SNP and all modes of inheritance were provided only by the parameterization proposed by C. H. Kao and Z. B. Zeng in 2002. However, if no dominance effect was included in the model, confidence intervals for SNP effects were correct for all parameterizations. Results based on real data showed that for both breeds the additive effects of polymorphisms were generally similar, except for LEPR, which had a different allele associated with increased fat content in Holstein-Friesians than in Jerseys. In both breeds, DGAT1 had the largest additive effect of the polymorphisms considered, but its effect on most milk traits was more pronounced in Jerseys than in Holstein-Friesians. Evidence of epistasis was found between LEPR and DGAT1, as well as between LEPR and BTN1A1, but only for milk content traits and only in the Holstein-Friesian breed. There was also more evidence for dominance in the Holstein-Friesian breed than in the Jersey breed.     

Identification of quantitative trait loci underlying milk traits in Spanish dairy sheep using linkage plus combined linkage disequilibrium and linkage analysis approaches

In this study, 2 procedures were used to analyze a data set from a whole-genome scan, one based on linkage analysis information and the other combing linkage disequilibrium and linkage analysis (LDLA), to determine the quantitative trait loci (QTL) influencing milk production traits in sheep. A total of 1,696 animals from 16 half-sib families were genotyped using the OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) and analysis was performed using a daughter design. Moreover, the same data set has been previously investigated through a genome-wide association (GWA) analysis and a comparison of results from the 3 methods has been possible. The linkage analysis and LDLA methodologies yielded different results, although some significantly associated regions were common to both procedures. The linkage analysis detected 3 overlapping genome-wise significant QTL on sheep chromosome (OAR) 2 influencing milk yield, protein yield, and fat yield, whereas 34 genome-wise significant QTL regions were detected using the LDLA approach. The most significant QTL for protein and fat percentages was detected on OAR3, which was reported in a previous GWA analysis. Both the linkage analysis and LDLA identified many other chromosome-wise significant associations across different sheep autosomes. Additional analyses were performed on OAR2 and OAR3 to determine the possible causality of the most significant polymorphisms identified for these genetic effects by the previously reported GWA analysis. For OAR3, the analyses demonstrated additional genetic proof of the causality previously suggested by our group for a single nucleotide polymorphism located in the α-lactalbumin gene (LALBA). In summary, although the results shown here suggest that in commercial dairy populations, the LDLA method exhibits a higher efficiency to map QTL than the simple linkage analysis or linkage disequilibrium methods, we believe that comparing the 3 analysis methods is the best approach to obtain a global picture of all identifiable QTL segregating in the population at both family-based and population-based levels.     

Short communication: Influence of composite casein genotypes on additive genetic variation of milk production traits and coagulation properties in Holstein-Friesian cows

The aim of the study was to quantify the effects of composite β- and κ-casein (CN) genotypes on genetic variation of milk coagulation properties (MCP); milk yield; fat, protein, and CN contents; somatic cell score; pH; and titratable acidity (TA) in 1,042 Italian Holstein-Friesian cows. Milk coagulation properties were defined as rennet coagulation time (RCT) and curd firmness (a₃₀). Variance components were estimated using 2 animal models: model 1 included herd, days in milk, and parity as fixed effects and animal and residual as random effects, and model 2 was model 1 with the addition of composite β- and κ-CN genotype as a fixed effect. Genetic correlations between RCT and a₃₀ and between these traits and milk production traits were obtained with bivariate analyses, based on the same models. The inclusion of casein genotypes led to a decrease of 47, 68, 18, and 23% in the genetic variance for RCT, a₃₀, pH, and TA, respectively, and less than 6% for other traits. Heritability of RCT and a₃₀ decreased from 0.248 to 0.143 and from 0.123 to 0.043, respectively. A moderate reduction was found for pH and TA, whereas negligible changes were detected for other milk traits. Estimates of genetic correlations were comparable between the 2 models. Results show that composite β- and κ-CN genotypes are important for RCT and a₃₀ but cannot replace the recording of MCP themselves.     

Short communication: influence of the muscle hypertrophy mutation of the myostatin gene on milk production traits and milk fatty acid composition in dual-purpose Belgian Blue dairy cattle

The objective of this study was to estimate genetic effects for the muscle hypertrophy mutation (mh) of the myostatin gene for conventional milk production traits and for milk fatty acid composition in dual-purpose Belgian Blue dairy cows. For the present study, only cows from a single herd, in which genotype frequencies were as balanced as possible (0.266 for +/+, 0.523 for mh/+, and 0.211 for mh/mh), were chosen to avoid confounding between herd and genotype effects. A total of 109 cows with 3,190 test-day records for fat, protein, and milk yields and 1,064 test-day records for saturated and monounsaturated fatty acids were used for the calculations. Variance component and gene effect estimations were performed via expectation-maximization REML and BLUP methods, respectively, using a multi-trait mixed test-day model with an additional fixed regression on the muscle hypertrophy genotype. Results showed that one copy of the wild-type “+” allele led to a significant additive effect of 26.35 g/d for fat yield. Significant dominance effects of 23.22 g/d for protein yield and 30.28 g/d for fat yield were also observed. In contrast, a nonsignificant trend was observed in favor of lower saturated fatty acid contents in milk for one copy of the mutant “mh” allele. Concerning milk, fat, and protein yields, our results confirmed literature results indicating a superior effect of the “+” allele compared with the mutant allele. Therefore, the selection of the “+” allele has the potential to increase conventional milk production traits in the dual-purpose Belgian Blue breed. However, when focus is given to milk fatty acid profile, a possible antagonistic effect between the benefit of the “+” allele for higher milk production and the “mh” allele for reduced saturated fatty acid content in milk should be confirmed in further studies.     

Milk composition traits in Jamunapari goats: Genetic parameters and effect of CSN1S1 gene variants

Milk protein polymorphism was analysed to improve the protein content in milk. The present study characterises the CSN1S1 gene and the effect of allelic combinations on milk composition traits in Jamunapari goats. The allelic variants obtained from sodium dodecyl sulfate polyacrylamide gel electrophoresis and polymerase chain reaction–restriction fragment length polymorphism were confirmed by cloning and sequencing. Genetic parameters were obtained from 518 records from 48 sires and 131 dams. The A, B and F alleles were observed in the population and the protein percentage in milk was significantly (P < 0.01) affected by allelic variants. The frequencies of A, B and F alleles were 0.456, 0.503 and 0.041, respectively. The protein content in milk was highest in the goats with AB genotype followed by AA > BB > BF > AF > FF. The goats with AB genotype had a significantly (P < 0.01) higher protein percentage in milk than goats with BF (t = 5.311, df = 113), AF (t = 8.13, df = 123) and FF (t = 9.55, df = 115) genotypes. The direct heritability for protein percentage was 0.441. Parity and season of birth had significant effects (P < 0.05) on the solids‐not‐fat percentage and lactose concentrations. The CSN1S1 AA, AB and BB genotypes should be selected for higher protein content and to improve milk quality and processing traits in Indian goats.     

Effect of goat milk composition on cheesemaking traits and daily cheese production

Cheese yield is strongly influenced by the composition of milk, especially fat and protein contents, and by the efficiency of the recovery of each milk component in the curd. The real effect of milk composition on cheesemaking ability of goat milk is still unknown. The aims of this study were to quantify the effects of milk composition; namely, fat, protein, and casein contents, on milk nutrient recovery in the curd, cheese yield, and average daily yield. Individual milk samples were collected from 560 goats of 6 different breeds. Each sample was analyzed in duplicate using the 9-laboratory milk cheesemaking assessment, a laboratory method that mimicked cheesemaking procedures, with milk heating, rennet addition, coagulation, curd cutting, and draining. Data were submitted to statistical analysis; results showed that the increase of milk fat content was associated with a large improvement of cheese yield because of the higher recovery of all milk nutrients in the curd, and thus a higher individual daily cheese yield. The increase of milk protein content affected the recovery of fat, total solids, and energy in the curd. Casein number, calculated as casein-to-protein ratio, did not affect protein recovery but strongly influenced the recovery of fat, showing a curvilinear pattern and the most favorable data for the intermediate values of casein number. In conclusion, increased fat and protein contents in the milk had an effect on cheese yield not only for the greater quantity of nutrients available but also for the improved efficiency of the recovery in the curd of all nutrients. These results are useful to improve knowledge on cheesemaking processes in the caprine dairy industry.     

Search for pleiotropic QTL on chromosome BTA6 affecting yield traits of milk production

The primary aim of this study was to investigate whether previous findings of similar quantitative trait loci (QTL) positions for correlated yield traits are due to a pleiotropic QTL. We applied a multitrait variance component based QTL mapping method to a dataset involving five granddaughter families from the German Holstein dairy cattle population. The marker map contained 16 microsatellite markers, distributed across chromosome BTA6. A chromosomewise significance threshold was used, because BTA6 is known to harbor QTL for several milk traits. To evaluate the results from the multivariate, across-family analysis, we also conducted single-family analyses using the least squares method of QTL estimation. The results provided two significant QTL findings at 49 and 64 cM for milk yield in different families and putative QTL at 68 cM for fat yield and at 71 cM for protein yield in another family. The results for fat and protein yield were confirmed by a univariate, across-family variance components analysis. The multivariate analysis of three bivariate trait combinations resulted in a significant pleiotropic QTL finding at 68 cM for fat yield and protein yield, bracketed by markers TGLA37 and FBN13. The estimates of variance contribution due to this QTL were 23% and 25%, respectively.     

Effects of bovine casein hydrolysate as a dry cow therapy on prevention and cure of bovine intramammary infection,milk production,and somatic cell count in the subsequent lactation

The primary objectives were to investigate the efficacy of bovine casein hydrolysate (bCNH) as a dry cow therapy at (1) preventing new intramammary infection (IMI) postpartum of all bacteria and coagulase-negative staphylococci (CNS), and (2) curing existing subclinical infections, mainly of CNS. The secondary objective was to measure the effects of bCNH on milk yield, composition, and somatic cell count (SCC) during the lactation postcalving. The trial was conducted as a randomized, blinded controlled experiment. Israeli Holstein dairy cows (n = 170) in first or higher lactations were recruited from 4 large commercial dairy herds. Cows were enrolled following clinical examination and bacteriological sampling of each quarter, which was the experimental unit. Random allocation was implemented at the cow level. All quarters of 100 cows were treated with 1 dose of bCNH (60 mg diluted in 20 mL of sterile solution) and those of 70 control cows were treated with saline solution. Clinical assessment of each cow's general appearance, teat-end leakage, and teat morphology was performed for 0, 1, 2, 3, 7, and 14 d after treatment, together with follow-up clinical observation and clinical examination of udder quarters. Quarter aseptic milk samples were obtained for bacteriological culture 48 h pretreatment, at time of treatment, and 3 and 5 d postcalving. Multivariable analyses were conducted to study the effects of bCNH on cure and prevention of IMI, adjusting for parity, farm, average of daily milk yield for 305 d, and average of monthly SCC values for 305 d of previous lactation. The odds of preventing IMI in cows treated with bCNH at dry-off were 2.15 times higher [95% confidence interval (CI): 1.15 to 4.00] than in cows treated with saline. Prevention was mostly of CNS. The odds of preventing CNS in cows treated with bCNH at dry-off were 2.20 times higher (95% CI: 1.58 to 3.07) than in control cows. The odds of curing IMI caused by CNS in cows treated with bCNH at dry-off were 4.80 times higher (95% CI: 0.75 to 30.75) than in saline-treated cows. Log SCC, adjusted to that of the previous lactation, was lower in the bCNH group compared with controls for 305 d in milk postcalving. The average milk yield per day for 305 d, adjusted to average daily milk yield of previous lactation, was higher by 2.1 kg in the bCNH group compared with controls (95% CI: 1.21 to 3.20). Clinical assessment of udders and cows posttreatment showed no negative effects of bCNH. Incidence of stillbirth, clinical mastitis, retained placenta, endometritis (5 to 12 d postcalving), ketosis, abortions, and reproduction did not differ between the 2 groups. Results suggest that a single intramammary administration of bCNH at dry-off effectively increases milk yield and lowers SCC, prevents new IMI during the dry period, and may be a beneficial alternative for curing existing IMI at dry-off, mainly by CNS.     

Investigating the genetic polymorphism of sheep milk proteins: a useful tool for dairy production

Sheep is the second most important dairy species after cow worldwide, and especially in the Mediterranean and Middle East regions. In some countries, the difficult environmental conditions require a peculiar adaptation and, in these contexts, sheep are able to provide higher quality protein than cattle. In the least‐developed countries, the amount of dairy sheep and ovine milk production is progressively increasing. In order to improve dairy productions, in particular those with local connotations, it is necessary to obtain in‐depth information regarding milk quality and rheological properties. The genetic polymorphisms of milk proteins are often associated with quantitative and qualitative parameters in milk and are potential candidate markers that should be included in breeding strategies similar to those already available for cattle. Due to the current and growing interest in this topic and considering the large amount of new information, the aim of this study was to review the literature on sheep milk protein polymorphisms with a particular emphasis on recent findings in order to give scientists useful support. Moreover, the effects of different protein variants on milk yield and composition are discussed. © 2014 Society of Chemical Industry     

Molecular characterization of bovine CSN1S2*B and extensive distribution of zebu-specific milk protein alleles in European cattle

The B allele of the bovine α_S2-casein gene (CSN1S2) was characterized at the molecular level and the distribution of zebu-specific milk protein alleles was determined in 26 cattle breeds originating from 3 continents. The CSN1S2 *B allele is characterized by a C → T transition affecting nucleotide 17 of exon 3, which leads to a change in the eighth amino acid of the mature protein, from Ser to Phe (i.e., TCC →TCC). DNA-based methods were developed to identify carriers of CSN1S2*B and the other alleles (CSN1S2 *A, C, and D) at the same locus. CSN1S2*B and other zebu-specific milk protein alleles and casein haplotypes are widely distributed in European cattle breeds, particularly those of southeastern origin. Alleles CSN1S2 *B and CSN3*H are important in searching for zebu imprints in European cattle breeds. Diversity estimates at the milk protein loci were highest in the zebus followed by southeastern European taurines. Anatolian Black had the highest number of zebu alleles among European taurines. Common, group, and intergroup haplotypes occurred in the breeds and demonstrated relationships that concurred with developmental histories, genetic makeup, and, in particular, exposed the extent of zebu influence on southeastern European cattle.     

Detection of quantitative trait loci affecting milk production traits on bovine chromosome 6 in a Chinese Holstein population by the daughter design

Fourteen microsatellite markers with a coverage of 63.5 cM on bovine chromosome 6 were selected, and 26 sire families with 2,260 daughters were analyzed for mapping quantitative trait loci (QTL) affecting 5 milk production traits in a Chinese Holstein population. In the analyses across 26 families and within the largest significant families with a one-QTL model fitted, a QTL near BMS470 was detected that affected fat yield at the 5% experiment-wide significance level. When a 2-QTL model was fitted in the across-family analysis, it was found that there might exist 2 QTL affecting the 3 yield traits, although the exact or empirical thresholds for the significance testing were unknown. In all analyses, the results for milk yield and protein yield were generally consistent, which might have resulted from the same genetic background for milk and protein yield.     

Association of candidate gene polymorphisms with milk technological traits,yield, composition,and somatic cell score in Italian Holstein-Friesian sires

Advances in DNA-based marker technology have enabled the identification of genomic regions underlying complex phenotypic traits in livestock species. The incorporation of detected quantitative trait loci into genetic evaluation provides great potential to enhance selection accuracies, hence expediting the genetic improvement of economically important traits. The objective of the present study was to investigate 96 single nucleotide polymorphisms (SNP) located in 53 candidate genes previously reported to have effects on milk production and quality traits in a population of highly selected Holstein-Friesian bulls. A total of 423 semen samples were used to genotype the bulls through a custom oligo pool assay. Forty-five SNP in 32 genes were found to be associated with at least 1 of the tested traits. Most significant and favorable SNP trait associations were observed for polymorphisms located in CCL3 and AGPAT6 genes for fat yield (0.037 and 0.033 kg/d, respectively), DGKG gene for milk yield (0.698 kg/d), PPARGC1A, CSN1S1, and AGPAT6 genes for fat percentage (0.127, 0.113, and 0.093%, respectively), GHR gene for protein (0.064%) and casein percentage (0.053%), and TLR4 gene for fat (0.090%), protein (0.066%), and casein percentage (0.050%). Somatic cell score was favorably affected by GHR (?0.095) and POU1F1 (?0.137), and interesting SNP-trait associations were observed for polymorphisms located in CSN2, POU1F1, and AGPAT6 genes for rennet coagulation time (?0.592, ?0.558, and ?0.462 min, respectively), and GHR and CSN2 genes for curd firmness 30 min after rennet addition (1.264 and 1.183 mm, respectively). In addition to the influence of individual SNP, the effects of composite genotypes constructed by grouping SNP according to their individual effects on traits considered in the analysis were also examined. Favorable and significant effects on milk traits were observed for 2 composite genotypes, one including 10 SNP and the other 4 SNP. The former was associated with an increase of milk (0.075 kg/d), fat (0.097 kg/d), protein (0.083 kg/d), and casein yields (0.065 kg/d), and the latter was associated with an increase of fat (0.244%), protein (0.071%), and casein percentage (0.047%). Although further research is required to validate the identified SNP loci in other populations and breeds, our results can be considered as a preliminary foundation for further replication studies on gene-assisted selection programs.     

Effects of composite beta- and kappa-casein genotypes on milk coagulation, quality, and yield traits in Italian Holstein cows

The aim of the study was to estimate the effect of the composite CSN2 and CSN3 genotypes on milk coagulation, quality, and yield traits in Italian Holstein cows. A total of 1,042 multiparous Holstein cows reared on 34 commercial dairy herds were sampled once, concurrently with monthly herd milk recording. The data included the following traits: milk coagulation time; curd firmness; pH and titratable acidity; fat, protein, and casein contents; somatic cell score; and daily milk, fat, and protein yields. A single-trait animal model was assumed with fixed effects of herd, days in milk, parity, composite casein genotype of CSN2 and CSN3 (CSN2-CSN3), and random additive genetic effect of an animal. The composite genotype of CSN2-CSN3 showed a strong effect on both milk coagulation traits and milk and protein yields, but not on fat and protein contents and other milk quality traits. For coagulation time, the best CSN2-CSN3 genotypes were those with at least one B allele in both the CSN2 and CSN3 loci. The CSN3 locus was associated more strongly with milk coagulation traits, whereas the CSN2 locus was associated more with milk and protein yields. However, because of the tight linkage between the 2 loci, the composite genotypes, or haplotypes, are more appropriate than the single-locus genotypes if they were considered for use in selection.     

Fine mapping and association analysis of a quantitative trait locus for milk production traits on Bos taurus autosome 4

To fine map a quantitative trait locus (QTL) affecting milk production traits previously associated with microsatellite RM188, we implemented an interval mapping analysis by using microsatellite markers in a large Israeli Holstein half-sib sire family, and linkage disequilibrium (LD) mapping in a large set of US Holstein bulls. Interval mapping located the target QTL to the near vicinity of RM188. For the LD mapping, we identified 42 single nucleotide polymorphisms (SNP) in 15 genes in a 12-Mb region on bovine chromosome 4. A total of 24 tag SNP were genotyped in 882 bulls belonging to the University of California Davis archival collection of Holstein bull DNA samples with predicted transmitted ability phenotypes. Marker-to-marker LD analysis revealed 2 LD blocks, with intrablock r² values of 0.10 and 0.46, respectively; outside the blocks, r² values ranged from 0.002 to 0.23. A standard additive/dominance model using the generalized linear model procedure of SAS and the regression module of HelixTree software were used to test marker-trait associations. Single nucleotide polymorphism 9 on ARL4A, SNP10 on XR_027435.1, SNP12 on ETV1, SNP21 on SNX13, and SNP24 were significantly associated with milk production traits. We propose the interval encompassing ARL4A and SNX13 genes as a candidate region in bovine chromosome 4 for a concordant QTL related to milk protein traits in dairy cattle. Functional studies are needed to confirm this result.