Genetic parameters of milk coagulation properties and their relationships with milk yield and quality traits in Italian Holstein cows

Milk coagulation properties (MCP) are an important aspect in assessing cheese-making ability. Several studies showed that favorable conditions of milk reactivity with rennet, curd formation rate, and curd strength, as well as curd syneresis, have a positive effect on the entire cheese-making process and subsequently on the ripening of cheese. Moreover, MCP were found to be heritable, but little scientific literature is available about their genetic aspects. The aims of this study were to estimate heritability of MCP and genetic correlations among MCP and milk production and quality traits. A total of 1,071 Italian Holstein cows (progeny of 54 sires) reared in 34 herds in Northern Italy were sampled from January to July 2004. Individual milk samples were collected during the morning milking and analyzed for coagulation time (RCT), curd firmness (a₃₀), pH, titratable acidity, fat, protein, and casein contents, and somatic cell count. About 10% of individual milk samples did not coagulate in 31 min, so they were removed from the analyses. Estimates of heritability for RCT and a₃₀ were 0.25 ± 0.04 and 0.15 ± 0.03, respectively. Estimates of genetic correlations between MCP traits and milk production traits were negligible except for a₃₀ with protein and casein contents (0.44 ± 0.10 and 0.53 ± 0.09, respectively). Estimates of genetic correlations between MCP traits and somatic cell score were strong and favorable, as well as those between MCP and pH and titratable acidity. Selecting for high casein content, milk acidity, and low somatic cell count might be an indirect way to improve MCP without reducing milk yield and quality traits.     

Genetic parameters for rennet- and acid-induced coagulation properties in milk from Swedish Red dairy cows

Milk coagulation is an important processing trait, being the basis for production of both cheese and fermented products. There is interest in including technological properties of these products in the breeding goal for dairy cattle. The aim of the present study was therefore to estimate genetic parameters for milk coagulation properties, including both rennet- and acid-induced coagulation, in Swedish Red dairy cattle using genomic relationships. Morning milk samples and blood samples were collected from 395 Swedish Red cows that were selected to be as genetically unrelated as possible. Using a rheometer, milk samples were analyzed for rennet- and acid-induced coagulation properties, including gel strength (G′), coagulation time, and yield stress (YS). In addition to the technological traits, milk composition was analyzed. A binary trait was created to reflect that milk samples that had not coagulated 40 min after rennet addition were considered noncoagulating milk. The cows were genotyped by using the Illumina BovineHD BeadChip (Illumina Inc., San Diego, CA). Almost 600,000 markers remained after quality control and were used to construct a matrix of genomic relationships among the cows. Multivariate models including fixed effects of herd, lactation stage, and parity were fitted using the ASReml software to obtain estimates of heritabilities and genetic and phenotypic correlations. Heritability estimates (h²) for G′ and YS in rennet and acid gels were found to be high (h² = 0.38–0.62) and the genetic correlations between rennet-induced and acid-induced coagulation properties were weak but favorable, with the exception of YS_rennet with G′_acid and YS_acid, both of which were strong. The high heritability (h² = 0.45) for milk coagulating ability expressed as a binary trait suggests that noncoagulation could be eliminated through breeding. Additionally, the results indicated that the current breeding objective could increase the frequency of noncoagulating milk and lead to deterioration of acid-induced coagulation through unfavorable genetic associations with protein content (0.38) and milk yield (−0.61 to −0.71), respectively. The outcome of this study suggests that by including more detailed compositional traits genetically associated with milk coagulation or by including milk coagulation properties directly within the breeding goal, it appears possible to breed cows that produce milk better suited for production of cheese and fermented products.     

Genetic parameters of coagulation properties, milk yield, quality, and acidity estimated using coagulating and noncoagulating milk information in Brown Swiss and Holstein-Friesian cows

The aim of this study was to estimate heritabilities of rennet coagulation time (RCT) and curd firmness (a₃₀) and their genetic correlations with test-day milk yield, composition (fat, protein, and casein content), somatic cell score, and acidity (pH and titratable acidity) using coagulating and noncoagulating (NC) milk information. Data were from 1,025 Holstein-Friesian (HF) and 1,234 Brown Swiss (BS) cows, which were progeny of 54 HF and 58 BS artificial insemination sires, respectively. Milk coagulation properties (MCP) of each cow were measured once using a computerized renneting meter and samples not exhibiting coagulation within 31 min after rennet addition were classified as NC milk. For NC samples, RCT was unobserved. Multivariate analyses, using Bayesian methodology, were performed to estimate the genetic relationships of RCT or a₃₀ with the other traits and statistical inference was based on the marginal posterior distributions of parameters of concern. For analyses involving RCT, a right-censored Gaussian linear model was used and records of NC milk samples, being censored records, were included as unknown parameters in the model implementing a data augmentation procedure. Rennet coagulation time was more heritable [heritability (h²) = 0.240 and h² = 0.210 for HF and BS, respectively] than a₃₀ (h² = 0.148 and h² = 0.168 for HF and BS, respectively). Milk coagulation properties were more heritable than a single test-day milk yield (h² = 0.103 and h² = 0.097 for HF and BS, respectively) and less heritable than milk composition traits whose heritability ranged from 0.275 to 0.275, with the only exception of fat content of BS milk (h² = 0.108). A negative genetic correlation, lower than −0.85, was estimated between RCT and a₃₀ for both breeds. Genetic relationships of MCP with yield and composition were low or moderate and favorable. The genetic correlation of somatic cell score with RCT in BS cows was large and positive and even more positive were those of RCT with pH and titratable acidity in both breeds, ranging from 0.80 to 0.94. Including NC milk information in the data affected the estimated correlations and decreased the uncertainty associated with the estimation process. On the basis of the estimated heritabilities and genetic correlations, enhancement of MCP through selective breeding with no detrimental effects on yield and composition seems feasible in both breeds. Milk acidity may play a role as an indicator trait for indirect enhancement of MCP.     

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.     

Genetic analysis of detailed milk protein composition and coagulation properties in Simmental cattle

The objective of this study was to estimate genetic parameters for milk protein fraction contents, milk protein composition, and milk coagulation properties (MCP). Contents of α_S1-, α_S2-, β-, γ-, and κ-casein (CN), β-lactoglobulin (β-LG), and α-lactalbumin (α-LA) were measured by reversed-phase HPLC in individual milk samples of 2,167 Simmental cows. Milk protein composition was measured as percentage of each CN fraction in CN (α_S1-CN%, α_S2-CN%, β-CN%, γ-CN%, and κ-CN%) and as percentage of β-LG in whey protein (β-LG%). Rennet clotting time (RCT) and curd firmness (a₃₀) were measured by a computerized renneting meter. Heritabilities for contents of milk proteins ranged from 0.11 (α-LA) to 0.52 (κ-CN). Heritabilities for α_S1-CN%, κ-CN%, and β-CN% were similar and ranged from 0.63 to 0.69, whereas heritability of α_S2-CN%, γ-CN%, and β-LG% were 0.28, 0.18, and 0.34, respectively. Effects of CSN2-CSN3 haplotype and BLG genotype accounted for more than 80% of the genetic variance of α_S1-CN%, β-CN%, and κ-CN% and 50% of the genetic variance of β-LG%. The genetic correlations among the contents of CN fractions and between CN and whey protein fractions contents were generally low. When the data were adjusted for milk protein gene effects, the magnitude of the genetic correlations among the contents of milk protein fractions markedly increased, indicating that they undergo a common regulation. The proportion of β-CN in CN correlated negatively with κ-CN% (r = −0.44). The genetic relationships between CN and whey protein composition were trivial. Low milk pH correlated with favorable MCP. Genetically, contents and proportions of α_S1- and α_S2-CN in CN were positively correlated with RCT. The relative proportion of β-CN in CN exhibited a genetic correlation with RCT of −0.26. Both the content and the relative proportion of κ-CN in CN did not correlate with RCT. Weak curds were genetically associated with increased proportions in CN of α_S1- and α_S2-CN, decreased contents of β-CN and κ-CN, and decreased proportion of κ-CN in CN. Negligible effects on the estimated correlations between a₃₀ and κ-CN contents or proportion in CN were observed when the model accounted for milk protein gene effects. Increasing β-CN and κ-CN contents and relative proportions in CN and decreasing the content and proportions of α_S1-CN and α_S2-CN and milk pH through selective breeding exert favorable effects on MCP.     

Short communication: Genetic aspects of milk differential somatic cell count in Holstein cows: A preliminary analysis

The aim of the present study was to assess genetic variation and heritability of a novel indicator of udder health, milk differential somatic cell count (DSCC), which represents the percentage of neutrophils plus lymphocytes in the total somatic cell count (SCC). Furthermore, we estimated genetic and phenotypic correlations of DSCC with other milk traits routinely measured in Italian Holstein cows. Besides DSCC, test-day data included milk yield, composition traits (i.e., fat, protein, casein, and lactose percentages), pH, milk urea nitrogen, and SCC. After editing, the final data set included 10,709 test-day records of 5,142 cows in 299 herds. Mean of DSCC was 62.07%, which means that macrophages were approximately 38% of total SCC. Comparing our results with the literature offered compelling evidence of the importance of acquiring information about the proportion of the different cell types in milk to better define the udder health status. In addition, our analysis revealed, for the first time, that DSCC is a heritable trait, and heritability (0.08 ± 0.02) was higher than that of traditional somatic cell score (0.04 ± 0.02). Nevertheless, heritability of DSCC is still low compared with milk yield and quality traits. Single-trait analysis within parity showed that DSCC was less heritable in primiparous than in multiparous cows, whereas bivariate analysis confirmed that DSCC and somatic cell score were 2 different traits, as their genetic and phenotypic correlations differed from unity. From a genetic point of view, the DSCC was positively weakly associated with milk yield, lactose percentage, and milk urea nitrogen, and negatively associated with pH. Our findings contributed to the understanding of the genetic background of DSCC and are a precursor to the potential use of DSCC in breeding programs to enhance cow resistance to mastitis. However, further research is needed to determine the weight this novel trait should receive in a selection program aimed to reduce udder health problems.     

Genetic parameters for clinical mastitis, somatic cell score, and production in the first three lactations of Swedish holstein cows

Using a mixed linear animal model, genetic parameters were estimated for clinical mastitis (MAST), lactation average somatic cell score (LSCS), and milk production traits in the first 3 lactations of more than 200,000 Swedish Holstein cows with first calving from 1995 to 2000. Heritability estimates for MAST (0.01 to 0.03) were distinctly lower than those for LSCS (0.10 to 0.14) and production traits (0.23 to 0.36). The genetic correlation between MAST and LSCS was high for all lactations (mean 0.70), implying that selection for low LSCS will reduce the incidence of mastitis. Undesirable genetic relationships with production were found for both MAST and LSCS with genetic correlations ranging from 0.01 to 0.45. This emphasizes the need for including udder health traits in the breeding goal. Genetic correlations across lactations for the same trait were positive and high for both MAST (>0.7), LSCS (>0.8), and production traits (>0.9), with the strongest correlations between second and third parity for all traits (>0.9 for udder health traits and close to unity for production traits).     

Mid-infrared spectroscopy predictions as indicator traits in breeding programs for enhanced coagulation properties of milk

The aims of this study were to investigate variation of milk coagulation property (MCP) measures and their predictions obtained by mid-infrared spectroscopy (MIR), to investigate the genetic relationship between measures of MCP and MIR predictions, and to estimate the expected response from a breeding program focusing on the enhancement of MCP using MIR predictions as indicator traits. Individual milk samples were collected from 1,200 Brown Swiss cows (progeny of 50 artificial insemination sires) reared in 30 herds located in northern Italy. Rennet coagulation time (RCT, min) and curd firmness (a₃₀, mm) were measured using a computerized renneting meter. The MIR data were recorded over the spectral range of 4,000 to 900 cm⁻¹. Prediction models for RCT and a₃₀ based on MIR spectra were developed using partial least squares regression. A cross-validation procedure was carried out. The procedure involved the partition of available data into 2 subsets: a calibration subset and a test subset. The calibration subset was used to develop a calibration equation able to predict individual MCP phenotypes using MIR spectra. The test subset was used to validate the calibration equation and to estimate heritabilities and genetic correlations for measured MCP and their predictions obtained from MIR spectra and the calibration equation. Point estimates of heritability ranged from 0.30 to 0.34 and from 0.22 to 0.24 for RCT and a₃₀, respectively. Heritability estimates for MCP predictions were larger than those obtained for measured MCP. Estimated genetic correlations between measures and predictions of RCT were very high and ranged from 0.91 to 0.96. Estimates of the genetic correlation between measures and predictions of a₃₀ were large and ranged from 0.71 to 0.87. Predictions of MCP provided by MIR techniques can be proposed as indicator traits for the genetic enhancement of MCP. The expected response of RCT and a₃₀ ensured by the selection using MIR predictions as indicator traits was equal to or slightly less than the response achievable through a single measurement of these traits. Breeding strategies for the enhancement of MCP based on MIR predictions as indicator traits could be easily and immediately implemented for dairy cattle populations where routine acquisition of spectra from individual milk samples is already performed.     

Genetic parameters for milk urea nitrogen in relation to milk production traits

The aim of this study was to estimate genetic parameters for test-day milk urea nitrogen (MUN) and its relationships with milk production traits. Three test-day morning milk samples were collected from 1,953 Holstein-Friesian heifers located on 398 commercial herds in the Netherlands. Each sample was analyzed for somatic cell count, net energy concentration, MUN, and the percentage of fat, protein, and lactose. Genetic parameters were estimated using an animal model with covariates for days in milk and age at first calving, fixed effects for season of calving and effect of test or proven bull, and random effects for herd-test day, animal, permanent environment, and error. Coefficient of variation for MUN was 33%. Estimated heritability for MUN was 0.14. Phenotypic correlation of MUN with each of the milk production traits was low. The genetic correlation was close to zero for MUN and lactose percentage (−0.09); was moderately positive for MUN and net energy concentration of milk (0.19), fat yield (0.41), protein yield (0.38), lactose yield (0.22), and milk yield (0.24), and percentage of fat (0.18), and percentage of protein (0.27); and was high for MUN and somatic cell score (0.85). Herd-test day explained 58% of the variation in MUN, which suggests that management adjustments at herd-level can reduce MUN. This study shows that it is possible to influence MUN by herd practice and by genetic selection.     

Genetic parameters of mid-infrared methane predictions and their relationships with milk production traits in Holstein cattle

Many countries have pledged to reduce greenhouse gases. In this context, the dairy sector is one of the identified sectors to adapt production circumstances to address socio-environmental constraints due to its large carbon footprint related to CH₄ emission. This study aimed mainly to estimate (1) the genetic parameters of 2 milk mid-infrared-based CH₄ proxies [predicted daily CH₄ emission (PME, g/d), and log-transformed predicted CH₄ intensity (LMI)] and (2) their genetic correlations with milk production traits [milk (MY), fat (FY), and protein (PY) yields] from first- and second-parity Holstein cows. A total of 336,126 and 231,400 mid-infrared CH₄ phenotypes were collected from 56,957 and 34,992 first- and second-parity cows, respectively. The PME increased from the first to the second lactation (433 vs. 453 g/d) and the LMI decreased (2.93 vs. 2.86). We used 20 bivariate random regression test-day models to estimate the variance components. Moderate heritability values were observed for both CH₄ traits, and those values decreased slightly from the first to the second lactation (0.25 ± 0.01 and 0.22 ± 0.01 for PME; 0.18 ± 0.01 and 0.17 ± 0.02 for LMI). Lactation phenotypic and genetic correlations were negative between PME and MY in both first and second lactations (?0.07 vs. ?0.07 and ?0.19 vs. ?0.24, respectively). More close scrutiny revealed that relative increase of PME was lower with high MY levels even reverting to decrease, and therefore explaining the negative correlations, indicating that higher producing cows could be a mitigation option for CH₄ emission. The PME phenotypic correlations were almost equal to 0 with FY and PY for both lactations. However, the genetic correlations between PME and FY were slightly positive (0.11 and 0.12), whereas with PY the correlations were slightly negative (?0.05 and ?0.04). Both phenotypic and genetic correlations between LMI and MY or PY or FY were always relatively highly negative (from ?0.21 to ?0.88). As the genetic correlations between PME and LMI were strong (0.71 and 0.72 in first and second lactation), the selection of one trait would also strongly influence the other trait. However, in animal breeding context, PME, as a direct quantity CH₄ proxy, would be preferred to LMI, which is a ratio trait of PME with a trait already in the index. The range of PME sire estimated breeding values were 22.1 and 29.41 kg per lactation in first and second parity, respectively. Further studies must be conducted to evaluate the effect of the introduction of PME in a selection index on the other traits already included in this index, such as, for instance, fertility or longevity.     

Genetic parameters for handling and milking temperament in Danish first-parity Holstein cows

The human-animal relationship in dairy cattle is reflected in the trait “temperament” in breeding programs and is mainly based on observations by farmers. However, farmers' knowledge of an individual cow's temperament decreases with an increased herd size, and this has been the case in many countries during the last decades. The aim of this study was to investigate if temperament recorded by classifiers and automatic milking systems is heritable, and estimate the genetic relationship with farmer-assessed temperament. Farmer-assessed temperament is defined as the overall assessment of the individual cows' temperament at milking and handling. Data on handling temperament were recorded by Danish classifiers from October 2016 to April 2017 on a 1 to 9 scale specially designed for this purpose. Data from automatic milking systems were recorded from January 2010 until April 2017, where connection time and number of attachments per teat were classified as milking temperament traits. Estimated heritabilities were relatively low for handling temperament (0.13) and farmer-assessed temperament (0.10). For milking temperament traits, connection time showed higher heritability than number of attachments per teat (0.36 and 0.26, respectively). The genetic correlation between farmer-assessed temperament and handling temperament was highly favorable (0.84). The genetic correlations between handling temperament and the 2 milking temperament traits, connection time and number of attachments per teat, were low (?0.02 and ?0.10, respectively). Moderate genetic correlations were estimated between farmer-assessed temperament and connection time (?0.29) and between farmer-assessed temperament and number of attachments per teat (?0.37). The genetic correlations and heritabilities suggest a basis for further investigations of the possibility of including handling or milking temperament traits (or both) in the breeding program for temperament in dairy cattle.     

Short communication: estimates of genetic variation of milk fatty acids in US Holstein cows

Interest in changing the milk fatty acid profile is growing. However, little is known about the genetic variability of milk fatty acids in the US Holstein population. Therefore, genetic parameters for milk fatty acids were estimated using a single-trait, mixed, linear animal model on 592 individual milk samples from 233 daughters of 53 sires in a cow herd genetically representative of the US Holstein population. Heritability (h²) and repeatability (r) estimates ± standard errors for yields of individual fatty acids ranged from 0.00 ± 0.08 (C4:0) to 0.43 ± 0.13 (C12:0) for heritabilities and from 0.21 ± 0.05 (C18:1) to 0.43 ± 0.05 (C12:0) for repeatabilities. Saturated (h² = 0.23 ± 0.12; r = 0.36 ± 0.05) and de novo synthesized fatty acids (C6:0 to C14:0; h² = 0.30 ± 0.13; r = 0.40 ± 0.05) had numerically higher estimates than did monounsaturated (h² = 0.09 ± 0.09; r = 0.22 ± 0.05) and polyunsaturated fatty acids (h² = 0.08 ± 0.09; r = 0.27 ± 0.05). For relative proportions of individual fatty acids, the greatest heritability and repeatability estimates were obtained for C8:0 (h² = 0.18 ± 0.12; r = 0.36 ± 0.05), C10:0 (h² = 0.22 ± 0.13; r = 0.46 ± 0.05), C12:0 (h² = 0.18 ± 0.12; r = 0.46 ± 0.05), C16:0 (h² = 0.09 ± 0.12; r = 0.48 ± 0.05), C16:1 (h² = 0.49 ± 0.13; r = 0.49 ± 0.05), and C18:0 (h² = 0.24 ± 0.11; r = 0.39 ± 0.05). Our results suggest the existence of genetic variability of milk fatty acids, in particular of medium-and long-chain fatty acids (C8:0 to C18:0), which could be used to improve the nutritional and textural properties of milk fat by selective breeding.     

Survival analysis as a statistical methodology for analyzing factors that affect milk coagulation time in Holstein-Friesian and Brown Swiss cows

The aim of this study was to propose the survival analysis technique as a statistical approach for the analysis of rennet coagulation time (RCT) able to make use of coagulating and noncoagulating (NC) milk information in order to estimate potential sources of variations that affects RCT. A total of 1,025 Italian Holstein-Friesian cows (HF; progeny of 54 sires) and 1,234 Brown Swiss cows (BS; progeny of 58 sires) reared in 34 and 38 herds, respectively, were milk-sampled once. Rennet coagulation time was analyzed with a semiparametric proportional hazard model (i.e., a Cox model), with the NC samples considered as censored records. Furthermore, a different censoring scenario, with a new end point at 18 min, was considered after the rearrangement of the time space originally used for the observation of RCT. The percentage of NC samples was almost 10% for HF and 3.5% for BS cows in in the 31 min set, whereas it increased to 44 and 24.9%, respectively, in the 18 min set. Estimated hazard ratios indicated that the most important factors affecting the coagulation process were herd, days in milk, casein number, and milk acidity (expressed in terms of titratable acidity) for both HF and BS, whereas the SCS was relevant only for BS. The survival model seems to be particularly suitable for this analysis, as it can properly account for censored and uncensored records and appropriately use all available information. Moreover, this methodology allows us to rearrange the time space used for the observation of RCT and to define alternative traits (i.e., RCT with an end point at 18 min). Our restriction of the time space and the increased percentage of censored records did not highlight any substantial differences in terms of the risk of coagulating with respect to the traditional 31 min testing time. Although further research is needed to investigate the effect of these sources of variation on cheese yield, our results indicate that casein number, acidity, and SCS may be used as indicator traits for enhancing the technological properties of bovine milk.     

Genetic parameters for major milk fatty acids and milk production traits of Dutch Holstein-Friesians

The objective of this study was to estimate genetic parameters for major milk fatty acids and milk production traits. One morning milk sample was collected from 1,918 Holstein-Friesian heifers located in 398 commercial herds in the Netherlands. Each sample was analyzed for total percentages of fat and protein, and for detailed fatty acid percentages (computed as fatty acid weight as a proportion of total fat weight). Intraherd heritabilities were high for C4:0 to C16:0, ranging from 0.42 for C4:0 to 0.71 for C10:0. Saturated and unsaturated C18 fatty acids had intraherd heritability estimates of approximately 0.25, except for C18:2 cis-9, trans-11, which was 0.42. Standard errors of the heritabilities were between 0.07 and 0.12. Genetic correlations were high and positive among C4:0 to C14:0, as well as among unsaturated C18, but correlations of C4:0 to C14:0 with unsaturated C18 were generally weak. The genetic correlation of C16:0 with fat percentage was positive (0.65), implying that selection for fat percentage should result in a correlated increase of C16:0, whereas unsaturated C18 fatty acids decreased with increasing fat percentage (−0.74). Milk fat composition can be changed by means of selective breeding, which offers opportunities to meet consumer demands regarding health and technological aspects.     

Short communication: Genetic correlation and heritability of milk coagulation traits within and across lactations in Holstein cows using multiple-lactation random regression animal models

Genetic parameters of milk rennet coagulation time (RCT) and curd firmness (a₃₀) among the first 3 lactations in Holstein cows were estimated. The data set included 39,960 test-day records from 5,216 Estonian Holstein cows (the progeny of 306 sires), which were recorded from April 2005 to May 2010 in 98 herds across the country. A multiple-lactation random regression animal model was used. Individual milk samples from each cow were collected during routine milk recording. These samples were analyzed for milk composition and coagulation traits with intervals of 2 to 3 mo in each lactation (7 to 305 DIM) and from first to third lactation. Mean heritabilities were 0.36, 0.32, and 0.28 for log-transformed RCT [ln(RCT)] and 0.47, 0.40, and 0.62 for a₃₀ for parities 1, 2, and 3, respectively. Mean repeatabilities for ln(RCT) were 0.53, 0.55, and 0.56, but 0.59, 0.61, and 0.68 for a₃₀ for parities 1, 2 and 3, respectively. Mean genetic correlations between ln(RCT) and a₃₀ were −0.19, −0.14, and 0.02 for parities 1, 2, and 3, respectively. Mean genetic correlations were 0.91, 0.79, and 0.99 for ln(RCT), and 0.95, 0.94, and 0.94 for a₃₀ between parities 1 and 2, 1 and 3, and 2 and 3, respectively. Due to these high genetic correlations, we concluded that for a proper genetic evaluation of milk coagulation properties it is sufficient to record RCT and a₃₀ only in the first lactation.     

Short communication: Association analysis of diacylglycerol acyltransferase (DGAT1) mutation on chromosome 14 for milk yield and composition traits,somatic cell score,and coagulation properties in Holstein bulls

The aim of the present study was to determine the allele frequencies of the diacylglycerol acyltransferase (DGAT1) K232A mutation in Italian Holstein bulls and to estimate the effect of the mutation on milk yield, composition, somatic cell score, and coagulation traits (rennet coagulation time and curd firmness). For this purpose, 349 Italian Holstein bulls were genotyped for the DGAT1 mutation on chromosome 14. Association analysis was performed by regressing the number of copies for the K allele on the deregressed estimated breeding value of the individual. Breeding values were calculated using field data routinely collected in Northeast Italy. The frequencies of the AA, KA, and KK genotypes were 59.6, 32.1, and 8.3%, respectively, and the minor allele frequency (K variant) was 24.7%. The K allele was significantly associated with greater fat yield and fat, protein, and casein percentages and with reduced protein:fat ratio. The association between the DGAT1 mutation and somatic cell score was not significant, whereas a favorable association between presence of the K allele and milk coagulation properties was found. Results from the present study confirmed the effect of the diallelic DGAT1 polymorphism K232A on milk production traits and, for the first time, provided evidence that this mutation also affects milk coagulation properties in the Italian Holstein breed.     

Genetic and phenotypic correlations between milk coagulation properties, milk production traits, somatic cell count, casein content, and pH of milk

Genetic and phenotypic correlations between milk coagulation properties (MCP: coagulation time and curd firmness), milk yield, fat content, protein content, ln(somatic cell count) (SCS), casein content, and pH of milk and heritability of these traits were estimated from data consisting of milk samples of 4664 Finnish Ayrshire cows sired by 91 bulls. In addition, differences in average estimated breeding values (EBV) for the above traits between the cows with noncoagulating (NC) milk and those with milk that coagulated (CO samples) were examined. The estimations were carried out to study the possibilities of indirect genetic improvement of MCP by use of the above characteristics. The genetic and phenotypic correlations between MCP and the milk production traits were low or negligible. The genetic associations between desirable MCP and low SCS were rather strong (-0.45 to 0.29). Desirable MCP correlated both genetically and phenotypically with low pH of milk (-0.51 to 0.50). The rather high heritability estimates for curd firmness in different forms (0.22 to 0.39), and the wide variation in the proportion of daughters producing NC milk between the sires (0 to 47%) suggested that noncoagulation of milk is partly caused by additive genetic factors. Based on the genetic correlations between curd firmness and SCS and the high EBV for SCS obtained for the cows with NC-milk, it is possible that the loci causing noncoagulation of milk and increasing somatic cell count of milk are closely linked or partly the same. One means to genetically improve MCP and to reduce the occurrence of NC milk could thus be selection for low somatic cell count of milk.     

Genetic parameters for claw disorders and the effect of preselecting cows for trimming

Claw disorders are important traits relevant to dairy cattle breeding from an economical and welfare point of view. Selection for reduced claw disorders can be based on hoof trimmer records. Typically, not all cows in a herd are trimmed. Our objectives were to estimate heritabilities and genetic correlations for claw disorders and investigate the effect of selecting cows for trimming. The data set contained 50,238 cows, of which 20,474 cows had at least one claw trimming record, with a total of 29,994 records. Six claw trimmers scored 14 different claw disorders: abscess (AB), corkscrew claw (CC), (inter-)digital dermatitis or heel erosion (DER), double sole (DS), hardship groove (HG), interdigital hyperplasia (IH), interdigital phlegmon (IP), sand crack (SC), super-foul (SF), sole hemorrhage (SH), sole injury (SI), sole ulcer (SU), white line separation (WLS), yellow discoloration of the sole (YD), and a combined claw disorder trait. Frequencies of the claw disorders for trimmed cows ranged from 0.1% (CC, YD, HG) to 23.8% (DER). More than half of the cows scored had at least one claw disorder. Heritability on the observed scale ranged from 0.02 (DS, SH) to 0.14 (IH) and on the underlying scale from 0.05 to 0.43 in trimmed cows. Genetic correlations between laminitis-related claw disorders were moderate to high, and the same was found for hygiene-related claw disorders. The effect of selecting cows for trimming was first investigated by including untrimmed cows in the analyses and assuming they were not affected by claw disorders. Heritabilities on the underlying scale showed only minor changes. Second, different subsets of the data were created based on the percentage of trimmed cows in the herd. Heritabilities for IH, DER, and SU tended to decrease when a higher percentage of cows in the herd were trimmed. Finally, a bivariate model with a claw disorder and the trait “trimming status” was used, but heritabilities were similar. Heritability for trimming status was relatively high (0.09). Genetic correlations of trimming status with claw disorders were generally moderate to high. To conclude, the effect of selecting cows for trimming on the heritability for claw disorders is negligible. Selecting herds with a high fraction of cows being trimmed tended to decrease heritability. Trimming status, as such, is a heritable trait and correlated with claw disorders and is therefore an interesting trait to include in the genetic evaluation.     

Genetic parameters for atypical reproductive patterns in dairy cows estimated from in-line milk progesterone profiles

Our aim was to estimate genetic parameters of atypical reproductive patterns and estimate their genetic correlation with milk production and classical fertility traits for commercial dairy cows. In contrast with classical fertility traits, atypical reproductive patterns based on in-line milk progesterone profiles might have higher heritability and lower genetic correlation with milk production. We had in-line milk progesterone profiles available for 12,046 cycles in 4,170 lactations of 2,589 primiparous and multiparous cows (mainly Holstein Friesian) from 14 herds. Based on progesterone profiles, 5 types of atypical reproductive patterns in a lactation were defined: delayed ovulation types I and II, persistent corpus luteum types I and II, and late embryo mortality. These atypical patterns were detected in 14% (persistent corpus luteum type II) to 21% (persistent corpus luteum type I) of lactations. In 47% of lactations, at least 1 atypical pattern was detected. Threshold model heritabilities for atypical reproduction patterns ranged between 0.03 and 0.14 and for most traits were slightly higher compared with classical fertility traits. The genetic correlation between milk yield and calving interval was 0.56, whereas genetic correlations between milk yield and atypical reproductive patterns ranged between ?0.02 and 0.33. Although most of these correlations between milk yield and atypical reproductive patterns are still unfavorable, they are lower compared with the correlations between classical fertility traits and milk yield. Therefore selection against atypical reproductive patterns may relax some constraints in current dairy breeding programs, to enhance genetic progress in both fertility and milk yield at a steady pace. However, as long as the target trait for fertility is calving interval, atypical reproductive patterns will not add additional value to the breeding goal in the near future due to the low number of available records.     

Short communication: Genetic parameters of individual fatty acids in milk of Canadian Holsteins

The objective of the present study was to estimate heritabilities of milk fatty acids (FA) and genetic and phenotypic correlations among milk FA and milk production traits in Canadian Holsteins. One morning milk sample was collected from each of 3,185 dairy cows between February and June 2010 from 52 commercial herds enrolled in Valacta (Ste-Anne-de-Bellevue, Quebec, Canada). Individual FA percentages (g/100 g of total FA) were determined for each sample by gas chromatography. After editing the data, 2,573 cows from 46 herds remained. Genetic parameters were estimated using multitrait animal models fitted under REML. The model included fixed effects of age at calving and stage of lactation each nested within parity and random effects of herd-year-season of calving, animal, and residual. The pedigree of animals with data was traced back 5 generations on both the male and female sides to account for relationships among animals. The estimates of heritability for individual FA ranged from 0.01 to 0.39, with standard errors ranging from 0.01 to 0.06. Generally, monounsaturated FA (MUFA) and saturated FA (SFA) showed higher heritability estimates than polyunsaturated FA (PUFA). Overall, SFA were negatively genetically correlated with MUFA and PUFA, whereas genetic correlations between MUFA and PUFA were positive. The SFA showed positive associations with fat yield and fat percentage, whereas unsaturated FA were negatively associated with fat yield and fat percentage. Bovine milk FA composition could be improved through genetic selection, and selection for MUFA or against SFA could alter the bovine milk fat profile in a desirable direction.