Symposium review: Multiple-trait single-step genomic evaluation for hoof health.

Hoof lesions represent an important issue in modern dairy herds, with reported prevalence in different countries ranging from 40 to 70%. This high prevalence of hoof lesions has both economic and social consequences, resulting in increased labor expenses and decreasing animal production, longevity, reproduction, health, and welfare. Therefore, a key goal of dairy herds is to reduce the incidence of hoof lesions, which can be achieved both by improving management practices and through genetic selection. The Canadian dairy industry has recently released a hoof health sub-index. This national genetic evaluation program for hoof health was achieved by creating a centralized data collection system that routinely transfers data recorded by hoof trimmers into a coherent and sustainable national database. The 8 most prevalent lesions (digital dermatitis, interdigital dermatitis, interdigital hyperplasia, heel horn erosion, sole hemorrhage, sole ulcer, toe ulcer, and white line lesion) in Canada are analyzed with a multiple-trait model using a single-step genomic BLUP method. Estimated genomic breeding values for each lesion are combined into a sub-index according to their economic value and prevalence. In addition, data recorded within this system were used to create an interactive management report for dairy producers by Canadian DHI, including the prevalence of lesions on farm, their trends over time, and benchmarks with provincial and national averages.

Genetic and genomic analyses of embryo production in dairy cattle.

The Canadian dairy industry has been using invivo and invitro assisted reproductive technologies to produce embryos. Technological improvements have helped increase the number and quality of embryos produced, but genetic and genomic tools for improving these traits have yet to be assessed for the Canadian Holstein population. Genetic parameters and a genome-wide association study were performed in Canadian Holstein for the total number of embryos (NE) and the number of viable embryos (VE). Results showed potential for genetic selection for both NE and VE, with heritability estimates (± s.e.) of approximately 0.15±0.01. Genetic correlations between the number of embryos produced using different procedures (invivo and invitro) suggested that a similar number of embryos should be expected from a donor regardless of the procedure used. A region on chromosome 11 of the bovine genome was found to be significantly associated with the number of embryos, indicating a potential regulatory role of this region on embryo production. Overall, these findings are of interest for the Canadian dairy industry because they provide useful information for breeders that are interested in producing embryos from the elite donors in their herds or in the population using assisted reproductive technologies.

Symposium review: The choice and collection of new relevant phenotypes for fertility selection.

In dairy production, high fertility contributes to herd profitability by achieving greater production and maintaining short calving intervals. Improved management practices and genetic selection have contributed to reversing negative trends in dairy cow fertility, but further progress is still required. Phenotypes included in current genetic evaluations are largely interval and binary traits calculated from insemination and calving date records. Several indicator traits such as calving, health, variation in body condition score, and longevity traits also apply to genetic improvement of fertility. Several fertility traits are included in the selection indices of many countries, but for improved selection, the development of novel phenotypes that more closely describe the physiology of reproduction and limit management bias could be more effective. Progesterone-based phenotypes can be determined from milk samples to describe the heritable interval from calving to corpus luteum activity, as well as additional measures of cow cyclicity. A fundamental component of artificial insemination practices is the observation of estrus. Novel phenotypes collected on estrous activity could be used to select for cows clearly displaying heat, as those cows are more likely to be inseminated at the right time and therefore have greater fertility performance. On-farm technologies, including in-line milk testing and activity monitors, may allow for phenotyping novel traits on large numbers of animals. Additionally, selection for improved fertility using traditional traits could benefit from refined and accurate recording and implementation of parameters such as pregnancy confirmation and reproductive management strategy, to differentiate embryonic or fetal loss, and to ensure selection for reproductive capability without producer intervention. Opportunities exist to achieve genetic improvement of reproductive efficiency in cattle using novel phenotypes, which is required for long-term sustainability of the dairy cattle population and industry.

Genetic correlations of mid-infrared-predicted milk fatty acid groups with milk production traits.

The objective of this research was to estimate the genetic correlations between milk mid-infrared-predicted fatty acid groups and production traits in first-parity Canadian Holsteins. Contents of short-chain, medium-chain, long-chain, saturated, and unsaturated fatty acid groupings in milk samples can be predicted using mid-infrared spectral data for cows enrolled in milk recording programs. Predicted fatty acid group contents were obtained for 49,127 test-day milk samples from 10,029 first-parity Holstein cows in 810 herds. Milk yield, fat and protein yield, fat and protein percentage, fat-to-protein ratio, and somatic cell score were also available for these test days. Genetic parameters were estimated for the fatty acid groups and production traits using multiple-trait random regression test day models by Bayesian methods via Gibbs sampling. Three separate 8- or 9-trait analyses were performed, including the 5 fatty acid groups with different combinations of the production traits. Posterior standard deviations ranged from <0.001 to 0.01. Average daily genetic correlations were negative and similar to each other for the fatty acid groups with milk yield (-0.62 to -0.59) and with protein yield (-0.32 to -0.25). Weak and positive average daily genetic correlations were found between somatic cell score and the fatty acid groups (from 0.25 to 0.36). Stronger genetic correlations with fat yield, fat and protein percentage, and fat-to-protein ratio were found with medium-chain and saturated fatty acid groups compared with those with long-chain and unsaturated fatty acid groups. Genetic correlations were very strong between the fatty acid groups and fat percentage, ranging between 0.88 for unsaturated and 0.99 for saturated fatty acids. Daily genetic correlations from 5 to 305 d in milk with milk, protein yield and percentage, and somatic cell score traits showed similar patterns for all fatty acid groups. The daily genetic correlations with fat yield at the beginning of lactation were decreasing for long-chain and unsaturated fatty acid groups and increasing for short-chain fatty acids. Genetic correlations between fat percentage and fatty acids were increasing at the beginning of lactation for short- and medium-chain and saturated fatty acids, but slightly decreasing for long-chain and unsaturated fatty acid groups. These results can be used in defining fatty acid traits and breeding objectives.

Genome-wide association study and in silico functional analysis of the number of embryos produced by Holstein donors.

Superovulation or ovum pick-up and in vitro fertilization are technologies used to produce an increased number of embryos from elite females. Embryo production traits have been shown to be heritable, but the genes that cause this variability have not yet been assessed. The main objectives of this study were to perform a genome-wide association study (GWAS) to find single nucleotide polymorphisms (SNP) associated with embryo production traits and to identify candidate genes affecting the number of embryos produced by Holstein donors in Canada that may provide insight into the regulation of embryo production. Breeding values were estimated and de-regressed for all donors and sires using a data set of 150,971 records of superovulation or ovum pick-up and in vitro fertilization. A total of 11,607 animals were genotyped, but of that number only 5,118 were genotyped with at least a 50K SNP panel and had a de-regressed estimated breeding value reliability of at least 10%. For the GWAS, 606,406 imputed SNP on 29 autosomal chromosomes were considered after applying quality control measures. A single-SNP univariate mixed linear animal model was used to perform the GWAS, and a 5% false discovery rate was applied to adjust for multiple testing. We found 36 and 14 significant SNP associated with the total number of embryos and the number of viable embryos, respectively, with most of them located on chromosome 11. Using these significant SNP, positional genes located within 10,000 bp upstream and downstream of the SNP were retrieved. Thirteen genes were harboring or near the significant SNP for the total number of embryos, 4 of them also being near the significant SNP for viable embryos. Some of these genes (CRB2, DENND1A, MAD1L1, NDUFA8, PTGS1) could be considered as potential positional candidate genes related to the number of embryos produced by a donor. This list will need to be validated in an independent population to confirm the role of the genes for embryo production.

Invited review: Genetics and claw health: Opportunities to enhance claw health by genetic selection.

Routine recording of claw health status at claw trimming of dairy cattle has been established in several countries, providing valuable data for genetic evaluation. In this review, we examine issues related to genetic evaluation of claw health; discuss data sources, trait definitions, and data validation procedures; and present a review of genetic parameters, possible indicator traits, and status of genetic and genomic evaluations for claw disorders. Different sources of data and traits can be used to describe claw health. Severe cases of claw disorders can be identified by veterinary diagnoses. Data from lameness and locomotion scoring, activity information from sensors, and feet and leg conformation traits are used as auxiliary traits. The most reliable and comprehensive information is data from regular hoof trimming. In genetic evaluation, claw disorders are usually defined as binary traits, based on whether or not the claw disorder was present (recorded) at least once during a defined time period. The traits can be specific disorders, composite traits, or overall claw health. Data validation and editing criteria are needed to ensure reliable data at the trimmer, herd, animal, and record levels. Different strategies have been chosen, reflecting differences in herd sizes, data structures, management practices, and recording systems among countries. Heritabilities of the most commonly analyzed claw disorders based on data from routine claw trimming were generally low, with ranges of linear model estimates from 0.01 to 0.14, and threshold model estimates from 0.06 to 0.39. Estimated genetic correlations among claw disorders varied from -0.40 to 0.98. The strongest genetic correlations were found among sole hemorrhage (SH), sole ulcer (SU), and white line disease (WL), and between digital/interdigital dermatitis (DD/ID) and heel horn erosion (HHE). Genetic correlations between DD/ID and HHE on the one hand and SH, SU, or WL on the other hand were, in most cases, low. Although some of the studies were based on relatively few records and the estimated genetic parameters had large standard errors, there was, with some exceptions, consistency among studies. Various studies evaluate the potential of various data soureces for use in breeding. The use of hoof trimming data is recommended for maximization of genetic gain, although auxiliary traits, such as locomotion score and some conformation traits, may be valuable for increasing the reliability of genetic evaluations. Routine genetic evaluation of direct claw health has been implemented in the Netherlands (2010); Denmark, Finland, and Sweden (joint Nordic evaluation; 2011); and Norway (2014), and other countries plan to implement evaluations in the near future.

Short communication: Prevalence of digital dermatitis in Canadian dairy cattle classified as high, average, or low antibody- and cell-mediated immune responders.

Lameness is a major animal welfare issue affecting Canadian dairy producers, and it can lead to production, reproduction, and health problems in dairy cattle herds. Although several different lesions affect dairy cattle hooves, studies show that digital dermatitis is the most common lesion identified in Canadian dairy herds. It has also been shown that dairy cattle classified as having high immune response (IR) have lower incidence of disease compared with those animals with average and low IR; therefore, it has been hypothesized that IR plays a role in preventing infectious hoof lesions. The objective of this study was to compare the prevalence of digital dermatitis in Canadian dairy cattle that were classified for antibody-mediated (AMIR) and cell-mediated (CMIR) immune response. Cattle (n = 329) from 5 commercial dairy farms in Ontario were evaluated for IR using a patented test protocol that captures both AMIR and CMIR. Individuals were classified as high, average, or low responders based on standardized residuals for AMIR and CMIR. Residuals were calculated using a general linear model that included the effects of herd, parity, stage of lactation, and stage of pregnancy. Hoof health data were collected from 2011 to 2013 by the farm's hoof trimmer using Hoof Supervisor software (KS Dairy Consulting Inc., Dresser, WI). All trim events were included for each animal, and lesions were assessed as a binary trait at each trim event. Hoof health data were analyzed using a mixed model that included the effects of herd, stage of lactation (at trim date), parity (at trim date), IR category (high, average, and low), and the random effect of animal. All data were presented as prevalence within IR category. Results showed that cows with high AMIR had significantly lower prevalence of digital dermatitis than cattle with average and low AMIR. No significant difference in prevalence of digital dermatitis was observed between high, average, and low CMIR cows. These results indicate that having more robust AMIR is associated with lower prevalence of digital dermatitis hoof lesions.

Genetic analysis for quality of frozen embryos produced by Holstein cattle donors in Canada.

The number of embryos produced by Holstein donors has been shown to be heritable, so it could be possible to genetically select for this trait to improve the efficiency of the assisted reproductive technology (ART) in dairy cattle. Another important parameter to consider for achieving good results from ART is embryo quality because embryos of good quality have more chance of producing live offspring. The possibility of using genetic selection for increasing the quality of embryo produced from ART has yet to be assessed. The objective of this study was, therefore, to perform a genetic analysis of embryo quality of Holstein donors in Canada using data recorded by Holstein Canada. The data set used was missing quality score data for embryos transferred fresh into a recipient, so the analyses were only performed for frozen embryos. With most traits in the Canadian dairy industry being evaluated with linear models, embryo quality was also evaluated with this class of models. However, considering the categorical nature of embryo quality, a threshold model was also evaluated. Embryo quality data were analyzed with either a univariate linear animal model or a univariate binomial threshold animal model. Genetic parameters estimated from the different models were comparable. A low heritability was found for the donor (0.04 ± <0.01) and the service sire (0.02 ± <0.01), but the repeatability estimate for the donor was higher (0.17), indicating that it was worthwhile to use a repeated records model. Overall, considering the low genetic parameters estimated, slow genetic progress is expected for the quality of frozen embryos produced by Canadian Holstein donors. Rank correlations were calculated between breeding values estimated from different models. High correlations were found between all models, indicating that no substantial re-ranking of the animals is expected from the different models. So, even though a threshold model is better suited for the analysis of categorical data, a linear model could be used for the analysis of embryo quality because it is less computationally demanding.

Genetic parameters for hoof health traits estimated with linear and threshold models using alternative cohorts.

A national genetic evaluation program for hoof health could be achieved by using hoof lesion data collected directly by hoof trimmers. However, not all cows in the herds during the trimming period are always presented to the hoof trimmer. This preselection process may not be completely random, leading to erroneous estimations of the prevalence of hoof lesions in the herd and inaccuracies in the genetic evaluation. The main objective of this study was to estimate genetic parameters for individual hoof lesions in Canadian Holsteins by using an alternative cohort to consider all cows in the herd during the period of the hoof trimming sessions, including those that were not examined by the trimmer over the entire lactation. A second objective was to compare the estimated heritabilities and breeding values for resistance to hoof lesions obtained with threshold and linear models. Data were recorded by 23 hoof trimmers serving 521 herds located in Alberta, British Columbia, and Ontario. A total of 73,559 hoof-trimming records from 53,654 cows were collected between 2009 and 2012. Hoof lesions included in the analysis were digital dermatitis, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, and white line disease. All variables were analyzed as binary traits, as the presence or the absence of the lesions, using a threshold and a linear animal model. Two different cohorts were created: Cohort 1, which included only cows presented to hoof trimmers, and Cohort 2, which included all cows present in the herd at the time of hoof trimmer visit. Using a threshold model, heritabilities on the observed scale ranged from 0.01 to 0.08 for Cohort 1 and from 0.01 to 0.06 for Cohort 2. Heritabilities estimated with the linear model ranged from 0.01 to 0.07 for Cohort 1 and from 0.01 to 0.05 for Cohort 2. Despite a low heritability, the distribution of the sire breeding values showed large and exploitable variation among sires. Higher breeding values for hoof lesion resistance corresponded to sires with a higher prevalence of healthy daughters. The rank correlations between estimated breeding values ranged from 0.96 to 0.99 when predicted using either one of the 2 cohorts and from 0.94 to 0.99 when predicted using either a threshold or a linear model.

Genetic analysis of groups of mid-infrared predicted fatty acids in milk.

The objective of this study was to investigate genetic variability of mid-infrared predicted fatty acid groups in Canadian Holstein cattle. Genetic parameters were estimated for 5 groups of fatty acids: short-chain (4 to 10 carbons), medium-chain (11 to 16 carbons), long-chain (17 to 22 carbons), saturated, and unsaturated fatty acids. The data set included 49,127 test-day records from 10,029 first-lactation Holstein cows in 810 herds. The random regression animal test-day model included days in milk, herd-test date, and age-season of calving (polynomial regression) as fixed effects, herd-year of calving, animal additive genetic effect, and permanent environment effects as random polynomial regressions, and random residual effect. Legendre polynomials of the third degree were selected for the fixed regression for age-season of calving effect and Legendre polynomials of the fourth degree were selected for the random regression for animal additive genetic, permanent environment, and herd-year effect. The average daily heritability over the lactation for the medium-chain fatty acid group (0.32) was higher than for the short-chain (0.24) and long-chain (0.23) fatty acid groups. The average daily heritability for the saturated fatty acid group (0.33) was greater than for the unsaturated fatty acid group (0.21). Estimated average daily genetic correlations were positive among all fatty acid groups and ranged from moderate to high (0.63-0.96). The genetic correlations illustrated similarities and differences in their origin and the makeup of the groupings based on chain length and saturation. These results provide evidence for the existence of genetic variation in mid-infrared predicted fatty acid groups, and the possibility of improving milk fatty acid profile through genetic selection in Canadian dairy cattle.

Prediction of milk fatty acid content with mid-infrared spectroscopy in Canadian dairy cattle using differently distributed model development sets.

The fatty acid profile of milk is a prevailing issue due to the potential negative or positive effects of different fatty acids to human health and nutrition. Mid-infrared spectroscopy can be used to obtain predictions of otherwise costly fatty acid phenotypes in a widespread and rapid manner. The objective of this study was to evaluate the prediction of fatty acid content for the Canadian dairy cattle population from mid-infrared spectral data and to compare the results produced by altering the partial least squares (PLS) model development set used. The PLS model development sets used to develop the predictions were reference fatty acids expressed as (1) grams per 100 g of fatty acid, (2) grams per 100 g of milk, (3) the natural logarithmic transform of grams per 100 g of milk, and (4) subsets of samples randomly selected by removing excess records around the mean to present a more uniform distribution, repeated 10 times. Gas chromatography measured fatty acid concentration and spectral data for 2,023 milk samples of 373 cows from 4 breeds and 44 herds were used in the model development. The coefficient of determination of cross-validation (Rcv2) increased when fatty acids were expressed on a per 100 g of milk basis compared with on a per 100 g of fat basis for all examined fatty acids. The logarithmic transformation used to create a more Gaussian distribution in the development set had little effect on the prediction accuracy. The individual fatty acids C12:0, C14:0, C16:0, C18:0, C18:1n-9 cis, and saturated, monounsaturated, unsaturated, short-chain, medium-chain, and long-chain fatty acid groups had (Rcv2) greater than 0.70. When model development was performed with subsets of the original samples, slight increases in (Rcv2) values were observed for the majority of fatty acids. The difference in (Rcv2) between the top- and bottom-performing prediction equation across the different subsets for a single predicted fatty acid was on average 0.055 depending on which samples were randomly selected to be used in the PLS model development set. Predictions for fatty acids with high accuracies can be used to monitor fatty acid contents for cows in milk recording programs and possibly for genetic evaluation.

Heritabilities of measured and mid-infrared predicted milk fat globule size, milk fat and protein percentages, and their genetic correlations.

The objective of this study was to estimate the heritability of milk fat globule (MFG) size and mid-infrared (MIR) predicted MFG size in Holstein cattle. The genetic correlations between measured and predicted MFG size with milk fat and protein percentage were also investigated. Average MFG size was measured in 1,583 milk samples taken from 254 Holstein cows from 29 herds across Canada. Size was expressed as volume moment mean (D[4,3]) and surface moment mean (D[3,2]). Analyzed milk samples also had average MFG size predicted from their MIR spectral records. Fat and protein percentages were obtained for all test-day milk samples in the cow's lactation. Univariate and bivariate repeatability animal models were used to estimate heritability and genetic correlations. Moderate heritabilities of 0.364 and 0.466 were found for D[4,3] and D[3,2], respectively, and a strong genetic correlation was found between the 2 traits (0.98). The heritabilities for the MIR-predicted MFG size were lower than those estimated for the measured MFG size at 0.300 for predicted D[4,3] and 0.239 for predicted D[3,2]. The genetic correlation between measured and predicted D[4,3] was 0.685; the correlation was slightly higher between measured and predicted D[3,2] at 0.764, likely due to the better prediction accuracy of D[3,2]. Milk fat percentage had moderate genetic correlations with both D[4,3] and D[3,2] (0.538 and 0.681, respectively). The genetic correlation between predicted MFG size and fat percentage was much stronger (greater than 0.97 for both predicted D[4,3] and D[3,2]). The stronger correlation suggests a limitation for the use of the predicted values of MFG size as indicator traits for true average MFG size in milk in selection programs. Larger samples sizes are required to provide better evidence of the estimated genetic parameters. A genetic component appears to exist for the average MFG size in bovine milk, and the variation could be exploited in selection programs.

Variation in fat globule size in bovine milk and its prediction using mid-infrared spectroscopy.

The objectives of this study were to investigate the sources of variation in milk fat globule (MFG) size in bovine milk and its prediction using mid-infrared (MIR) spectroscopy. Mean MFG size was measured in 2,076 milk samples from 399 Ayrshire, Brown Swiss, Holstein, and Jersey cows, and expressed as volume moment mean (D[4,3]) and surface moment mean (D[3,2]). The mid-infrared spectra of the samples and milk performance data were also recorded during routine milk recording and testing. The effects of breed, herd nested within breed, days in milk, season, milking period, age at calving, parity, and individual animal on the variation observed in MFG size were investigated. Breed, herd nested within breed, days in milk, season, and milking period significantly affected mean MFG size. Milk fat globule size was the largest at the beginning of lactation and subsequently decreased. Milk samples with the smallest MFG on average came from Holstein cows, and those with the largest were from Jersey and Brown Swiss cows. Partial least squares regression was used to predict MFG size from MIR spectra of samples with a calibration data set containing 2,034 and 2,032 samples for D[4,3] and D[3,2], respectively. Coefficients of determination of cross validation for D[4,3] and D[3,2] prediction models were 0.51 and 0.54, respectively. The associated ratio of performance deviation values were 1.43 and 1.48 for D[4,3] and D[3,2], respectively. With these models, individual mean MFG size could not be accurately predicted, but results may be sufficient to screen samples for having either small or large MFG on average. Significant but low correlations of D[4,3] and D[3,2] with milk fat yield were estimated (0.16 and 0.21, respectively). Significant and moderate Pearson correlation coefficients for fat percent with D[4,3] and D[3,2] were assessed (0.34 and 0.36, respectively). This correlation was greater between milk fat percentage and predicted MFG size than with measured MFG size with coefficients of 0.47 and 0.49 for D[4,3] and D[3,2], respectively. The MIR prediction equations are potentially overusing the correlation between fat and MFG size and exploiting the strong relationship between the MIR spectra and total milk fat. However, the predictions of MFG size are able to determine variation in mean globule size beyond what would be achieved just by looking at the correlation with fat production.

Genetic analysis of superovulatory response of Holstein cows in Canada.

Superovulation of dairy cattle is frequently used in Canada. The cost of this protocol is high, and so is the variability of the outcome. Knowing the superovulatory potential of a donor cow could influence the breeder's decision to superovulate it or not. The main objective of this study was to perform a genetic analysis for superovulatory response of Holstein cows in Canada using data recorded by Holstein Canada, and to investigate if these data could be used for genetic evaluation. Data contained the total number of embryos and the number of viable embryos from every successful flushing performed across Canada. After editing, 137,446 records of superovulation performed between 1992 and 2014 were analyzed. A univariate repeatability animal model analysis was performed for both total number of embryos and number of viable embryos. Because both data and residuals did not follow a normal distribution, records were subject to either logarithmic or Anscombe transformation. Using logarithmic transformation, heritability estimates (SE) of 0.15 (0.01) and 0.14 (0.01) were found for total number of embryos and number of viable embryos, respectively. Using Anscombe transformation, heritability estimates (SE) of 0.17 (0.01) and 0.14 (0.01) were found for total number of embryos and number of viable embryos, respectively. The genetic correlation between the 2 traits was estimated at 0.97 using logarithmic transformation and 0.95 using Anscombe transformation. Breeding values were estimated for 54,463 cows, and 3,513 sires. Only estimated breeding values of sires having a reliability higher than 40% were considered for estimated breeding values correlations with other routinely evaluated traits. The results showed that selection for a higher response to superovulation would lead to a slight decrease in milk production, but an improvement for functional traits, including all reproduction traits. In all cases, the estimated correlations are either low or modest. We conclude that genetic selection for increased superovulatory response in donors is possible; daughters of sires with high estimated breeding values for superovulatory response will tend to yield more embryos, whereas the additive effect of service sire seems not to contribute to the variability of the 2 superovulation traits and was not significantly correlated with the additive effect of the donor.

Effect of dairy farming system, herd, season, parity, and days in milk on modeling of the coagulation, curd firming, and syneresis of bovine milk.

The objectives of this study were to characterize the variation in curd firmness model parameters obtained from coagulating bovine milk samples, and to investigate the effects of the dairy system, season, individual farm, and factors related to individual cows (days in milk and parity). Individual milk samples (n = 1,264) were collected during the evening milking of 85 farms representing different environments and farming systems in the northeastern Italian Alps. The dairy herds were classified into 4 farming system categories: traditional system with tied animals (29 herds), modern dairy systems with traditional feeding based on hay and compound feed (30 herds), modern dairy system with total mixed ration (TMR) that included silage as a large proportion of the diet (9 herds), and modern dairy system with silage-free TMR (17 herds). Milk samples were analyzed for milk composition and coagulation properties, and parameters were modeled using curd firmness measures (CFt) collected every 15 s from a lacto-dynamographic analysis of 90 min. When compared with traditional milk coagulation properties (MCP), the curd firming measures showed greater variability and yielded a more accurate description of the milk coagulation process: the model converged for 93.1% of the milk samples, allowing estimation of 4 CFt parameters and 2 derived traits [maximum CF (CF(max)) and time from rennet addition to CF(max) (t(max))] for each sample. The milk samples whose CFt equations did not converge showed longer rennet coagulation times obtained from the model (RCT(eq)) and higher somatic cell score, and came from less-productive cows. Among the sources of variation tested for the CFt parameters, dairy herd system yielded the greatest differences for the contrast between the traditional farm and the 3 modern farms, with the latter showing earlier coagulation and greater instant syneresis rate constant (k(SR)). The use of TMR yielded a greater tmax because of a higher instant curd-firming rate constant (k(CF)). Season of sampling was found to be very important, yielding higher values during winter for all traits except k(CF) and k(SR). All CFt traits were affected by individual cow factors. For parity, milk produced by first-lactation cows showed higher k(CF) and k(SR), but delays in achieving CF(max). With respect to stage of lactation, RCT(eq) and potential asymptotic CF increased during the middle of lactation and stabilized thereafter, whereas the 2 instant rate constants presented the opposite pattern, with the lowest (k(CF)) and highest (k(SR)) values occurring in mid lactation. The new challenge offered by prolonging the test interval and individual modeling of milk technological properties allowed us to study the effects of parameters related to the environment and to individual cows. This novel strategy may be useful for investigating the genetic variability of these new coagulation traits.

Modeling of coagulation, curd firming, and syneresis of milk from Sarda ewes.

This study investigated the modeling of curd-firming (CF) over time (CF(t)) of sheep milk. Milk samples from 1,121 Sarda ewes from 23 flocks were analyzed for coagulation properties. Lactodynamographic analyses were conducted for up to 60 min, and 240 CF individual observations from each sample were recorded. Individual sample CFt equation parameters (RCT(eq), rennet coagulation time; CF(P), asymptotic potential value of curd firmness; k(CF), curd-firming instant rate constant; and k(SR), curd syneresis instant rate constant) were estimated, and the derived traits (CF(max), the point at which CF(t) attained its maximum level, and tmax, the time at which CF(max) was attained) were calculated. The incidence of noncoagulating milk samples was 0.4%. The iterative estimation procedure applied to the individual coagulation data showed a small number of not-converged samples (4.4%), which had late coagulation and an almost linear pattern of the ascending part of the CF(t) curve that caused a high value of CF(P), a low value of k(CF), and a high value of k(SR). Converged samples were classified on the basis of their CF(t) curves into no-k(SR) (18.0%), low-k(SR) (72.6%), and high-k(SR) (4.5%). A CF(t) that was growing continuously because of the lack of the syneresis process characterized the no-k(SR) samples. The high-k(SR) samples had a much larger CFP, a smaller k(CF), and an anticipation of tmax, whereas the low-k(SR) samples had a fast k(CF) and a slower k(SR). The part of the average CF(t) curves that showed an increase was similar among the 3 different syneretic groups, whereas the part that decreased was different because of the expulsion of whey from the curd. The traditional milk coagulation properties recorded within 30 min were not able to detect any appreciable differences among the 4 groups of coagulating samples, which could lead to a large underestimation of the maximum CF of all samples (if predicted by a30), with the exception of the no-k(SR) samples. Large individual variability was found and was likely caused by the effects of the dairy system, such as flock size (on CF(max), t(max), and % ewes with no-k(SR) milk), flock within flock size (representing 11 to 43% of total variance for % ewes with no-k(SR) milk and CF(max), respectively), days in milk (on all model parameters and CF(max)), parity (on RCT(eq), k(SR), and CF(max)), daily milk yield (on RCT(eq) and CF(max)), and position of the individual pendulum that significantly affected model parameters and derived traits. In conclusion, the results showed that the modeling of coagulation, curd-firming, and syneresis is a suitable tool to achieve a deeper interpretation of the coagulation and curd-firming processes of sheep milk and also to study curd syneresis.

Fertility traits of purebred Holsteins and 2- and 3-breed crossbred heifers and cows obtained from Swedish Red, Montbéliarde, and Brown Swiss sires.

The objective of this study was to compare fertility traits of heifers and primiparous cows from Swedish Red × Holstein (SR × HO; n=634 and 581, respectively), Montbéliarde × Holstein (MO × HO; n=126 and 114, respectively), Brown Swiss × Holstein (BS × HO; n=59 and 50, respectively), and MO × (SR × HO) (n=241 and 139, respectively) crossbreds, versus those of Holstein heifers and cows (HO; n=3,483 and 2,549, respectively). Animals were born between 2007 and 2011, and belonged to 4 herds located in northern Italy. Heifers were compared for age at first service (AFS), age at first conception (AFC), interval between first service and conception (IFC), nonreturn rate at 56 d after first service (NR56), conception rate at first service (CR), and number of inseminations required for conception (INS). The same traits were evaluated in primiparous cows, except that AFS and AFC were replaced with days at first service (DFS) and days open (DO). The AFS, AFC, IFC, DFS, and DO traits were continuous variables and were thus analyzed under a proportional hazards Cox model that properly accounted for censoring among cows that were culled or failed to conceive. The NR56, CR, and INS traits were analyzed as binary traits using logistic regression. Our results indicated that, among heifers, SR × HO crossbreds had a better chance of having an earlier first service and conceiving earlier than HO, with hazard ratios (HR) of 1.31 for AFS and 1.34 for AFC. Similarly, MO × (SR × HO) crossbreds differed from HO heifers in this regard (HR=1.18 and 1.24, respectively). For the primiparous cows, all crossbreds showed significant differences for DFS, DO, and IFC relative to purebred HO, with the exception of the BS × HO crossbreds. The MO × HO, SR × HO, and MO × (SR × HO) crossbred cows showed increased chances of having fewer DFS (HR=1.40, 1.30, and 1.27, respectively), fewer DO (HR=1.59, 1.43, and 1.58, respectively), and fewer IFC (HR=1.52, 1.26, and 1.39, respectively) than HO cows. All crossbred genotypes, including BS × HO cows, showed higher probabilities for higher NR56, higher CR, and lower INS than purebred HO cows. Together, these findings indicate that the studied crossbred cows have higher reproductive potential than Holsteins.

Milk quality, coagulation properties, and curd firmness modeling of purebred Holsteins and first- and second-generation crossbred cows from Swedish Red, Montbéliarde, and Brown Swiss bulls.

The objective of the present study was to investigate how the crossbreeding of Holstein (HO) cows with bulls from Nordic and Alpine European breeds affect milk quality traits, traditional milk coagulation properties (MCP), and curd firmness modeling obtained from individual milk samples. A total of 506 individual milk samples were collected from evening milking at 3 commercial farms located in Northern Italy. Over the past decade, the 3 farms have followed crossbreeding programs in part of their herds, whereas the remainder of the animals consisted of purebred HO. The basic scheme was a 3-breed rotation based on the use of Swedish Red (SR) semen on HO cows (SR × HO), the use of Montbéliarde (MO) semen on first-cross cows [MO × (SR × HO)], and the use of HO semen in the third cross. In all herds, a smaller proportion of purebred HO were mated to M and Brown Swiss (BS) bulls, and these first crosses were mated to SR and MO bulls, respectively. Milk samples were analyzed for milk composition and MCP, and parameters for curd firmness were modeled. Compared with purebred HO, crossbred cows produced less milk with lower lactose content, higher fat and protein content, and a tendency for higher casein content. Crossbred cows generally produced milk with a more favorable curd-firming rate (k20) and curd firmness 30 min after rennet addition, among traditional MCP, and better trends of curd firmness measures as shown by model parameters: estimated rennet coagulation time, asymptotical potential value of curd firmness, and curd-firming instant rate constant. Among crossbred cows, SR × HO presented longer rennet coagulation time compared with MO × HO and BS × HO cows, and MO × HO showed shorter k20 compared with BS × HO cows. Among second-generation cows, those sired by SR bulls showed a lower incidence of noncoagulated samples, higher curd firmness 30 min after rennet addition and asymptotical potential value of curd firmness, and faster curd-firming instant rate constant compared with animals sired by MO bulls. Our results revealed that different sire breeds were characterized by specific technological aptitudes, but that these were not strictly related to other milk quality traits. Furthermore, the favorable characteristics (in terms of the quality and technological properties of milk) could be maintained in the third generation of 3-way crosses without negative effects on milk yield, even though the HO heritage had been reduced from 50 to 25%. Our findings, therefore, suggest that different types of sires can be chosen (depending on the intended use of the milk) to ensure the optimization of farm crossbreeding programs.

Quality traits and modeling of coagulation, curd firming, and syneresis of sheep milk of Alpine breeds fed diets supplemented with rumen-protected conjugated fatty acid.

The aim of this study was to test the modeling of curd-firming (CF) measures and to compare the sheep milk of 3 Alpine breeds supplemented with or without rumen-protected conjugated linoleic acid (rpCLA). Twenty-four ewes of the Brogna, Foza, and Lamon breeds were allotted to 6 pens (2 pens/breed) and fed a diet composed of corn grain, corn silage, dried sugar beet pulp, soybean meal, wheat bran, wheat straw, and a vitamin-mineral mixture. The rpCLA supplement (12 g/d per ewe plus 4 g/d for each lamb older than 30 d) was mixed into the diet of 1 pen per sheep breed (3 pens/treatment) to provide an average of 0.945 and 0.915 g/d per ewe of the cis-9,trans-11 C18:2 and trans-10,cis-12 C18:2 conjugated linoleic acid isomers, respectively. The trial started at 38 ± 23 d after parturition, and individual morning milk samples were collected on d 16, 23, 37, 44, and 59 of the trial. Milk samples were analyzed for composition, and duplicate samples were assessed for milk coagulation properties (MCP). A total of 180 CF measures for each sample (1 every 15s) were recorded. Model parameters were the rennet coagulation time, the asymptotic potential CF, the CF instant rate constant, the syneresis instant rate constant, the maximum CF achieved within 45 min (CFmax), and the time at achievement of CFmax. The data were analyzed using a hierarchical model that considered the fixed effects of breed, diet, lamb birth, and initial days in milk, which were tested on individual ewe (random) variance; the fixed effect of sampling day, which was tested on the within-ewe sample (random) variance; and the fixed effect of instrument or cuvette position (only for MCP), which was tested on the residual (replicates within samples) variance. The local Alpine sheep breeds displayed similar milk compositions, traditional MCP, and CF modeling parameters. Supplementation with rpCLA triggered changes in milk composition and worsened MCP (e.g., delayed rennet coagulation time, slower CF instant rate constant, and a doubling of syneresis instant rate constant), but did not influence potential CF. Overall, our results indicate that rpCLA supplementation reduced the actual maximum CF (CFmax) but did not modify the interval between rennet addition and CFmax or time to CFmax.

Milk protein composition in purebred Holsteins and in first/second-generation crossbred cows from Swedish Red, Montbeliarde and Brown Swiss bulls.

The aim of this study was to analyze milk protein composition in purebred and crossbred dairy cattle and estimate the effects of individual sources of variation on the investigated traits. Milk samples were collected from 505 cows from three commercial farms located in Northern Italy, some of which had originated from crossbreeding programs, although most were purebred Holsteins (HO). The basic crossbreeding scheme was a three-breed rotational system using Swedish Red (SR) semen on HO cows (SR×HO), Montbeliarde (MO) semen on SR×HO cows (MO×(SR×HO)) and HO semen again on MO×(SR×HO) cows. A smaller number of purebred HO from each of the herds were mated inverting the breed order (MO×HO and SR×(MO×HO)) or using Brown Swiss (BS) bulls (BS×HO) then MO bulls (MO×(BS×HO)). Milk samples were analyzed by reverse-phase HPLC to obtain protein fraction amounts (g/l) and proportions (% of total true protein). Traits were analyzed using a linear model, which included the fixed effects of herd-test-day (HTD), parity, days in milk and breed combination. Results showed that milk protein fractions were influenced by HTD, stage of lactation, parity and breed combination. The increase in protein concentration during lactation was due in particular to ß-casein (ß-CN), α S1-CN and ß-lactoglobulin (ß-LG). The higher protein content of primiparous milk was mainly due to higher concentrations of all casein fractions. The milk from crossbred cows had higher contents and proportions of κ-CN and α-lactalbumin (α-LA), lower proportions of ß-LG and greater proportion of caseins/smaller in whey proteins on milk true protein than purebred HO. The three-way crossbreds differed from two-way crossbreds only in having greater proportions of α-LA in their milk. Of the three-way crossbreds, the SR sired cows yielded milk with a smaller content and proportion of ß-LG than the MO sired cows, and, consequently, a higher proportion of caseins than whey proteins. Results from this study support the feasibility of using crossbreeding programs to alter milk protein profiles with the aim of improving milk quality and cheese-making properties.