Short communication: Genetic analysis of dairy bull fertility from field data of Brown Swiss cattle

The aim of this study was to estimate heritability and repeatability of dairy bull fertility in Italian Brown Swiss cattle. Bull fertility indicators were calving per service and nonreturn rate at 56 d after service. Data included 124,206 inseminations carried out by 86 technicians on 28,873 heifers and cows in 1,400 herds. Services were recorded from 1999 to 2008 and were performed with semen from 306 AI Brown Swiss bulls. Data were analyzed with a threshold animal model, which included the fixed effects of parity by class of days in milk of the inseminated cow (age at insemination for heifers), year-season of insemination, and status of the service bull at the time of insemination (progeny testing or proven), and the random effects of herd, technician, additive genetic, and permanent environment of inseminated heifer/cow and service bull, and residual. Also, genetic covariance between heifer/cow and service bull effects was considered in the model. Heritability and repeatability were 0.0079 and 0.0100 for nonreturn rate at 56 d after service, and 0.0153 and 0.0202 for calving per service, respectively. The low estimates obtained in the present study indicate that selection for male fertility using field data is hardly pursuable.     

Genetic evaluations for measures of the milk-flow curve in the Italian Brown Swiss population

The objective of this study was to estimate heritabilities and genetic correlations between milk-release parameters, somatic cell score, milk yield, and udder functional traits in the Italian Brown Swiss population. Data were available from 37,511 cows over a span of 12 yr (1997-2008) from 1,592 herds. Milking flows were recorded for each individual once during lactation. Three different analyses were performed to estimate variance components for all the traits of interest. The first analysis included single control data milk yield, somatic cell score, maximum milk flow, average milk flow, time of plateau, decreasing time, and total milking time, whereas the second analysis included milk-release parameters as well as total udder score, udder depth, and 305-d milk yield and somatic cell score as dependent variables. The third analysis included total milking time, 305-d milk yield and somatic cell score, total udder score, udder depth, and ratios of maximum milk flow over total milking time (R1), time of plateau (R2), and decreasing time (R3) to estimate the relationship between the shape of the milk-release curves and important milking traits. Results from the first and second analysis found similar heritabilities for milkability traits ranging from 0.05 to 0.41 with genetic correlations between production traits and flow traits ranging from low to moderate values. Positive genetic correlations were found among production, somatic cell score, and milkability traits. The third analysis showed that R1 had the greatest heritability of the ratio traits (0.37) with large genetic correlations with R2 and R3, a low correlation with 305-d somatic cell score, and no correlation with 305-d milk yield. Estimated responses to selection over 5 generations were also calculated using different indexes, which included either flow or ratio traits. The results of this study show that it is possible to use information collected through portable flowmeters to improve milkability traits. Using a set of variables or traits to describe the overall release of milk can be an advantageous selection strategy to decrease management costs while maintaining milk production.     

Genetic parameters for fertility of dairy heifers and cows at different parities and relationships with production traits in first lactation

The objectives of this study were to estimate genetic parameters for fertility of Brown Swiss cattle, considering reproductive measures in different parities as different traits, and to estimate relationships between production traits of first lactation and fertility of heifers and first-parity and second-parity cows. Reproductive indicators were interval from parturition to first service, interval from first service to conception, interval from parturition to conception, number of inseminations to conception, conception rate at first service, and nonreturn rate at 56 d after first service. Production traits were peak milk yield (pMY), lactation milk yield, and lactation length (LL). Data included 37,546 records on heifers, and 24,098 and 15,653 records on first- and second-parity cows, respectively. Cows were reared in 2,035 herds, calved from 1999 to 2007, and were progeny of 527 AI bulls. Gibbs sampling was implemented to obtain (co)variance components using both univariate and bivariate threshold and censored linear sire models. Estimates of heritability for reproductive traits in heifers (0.016 to 0.026) were lower than those in first-parity (0.017 to 0.142) and second-parity (0.026 to 0.115) cows. Genetic correlations for fertility in first- and second-parity cows were very high (>0.920), whereas those between heifers and lactating cows were moderate (0.348 to 0.709). The latter result indicates that fertility in heifers is a different trait than fertility in lactating cows, and hence it cannot be used as robust indicator of cow fertility. Heifer fertility was not related to production traits in first lactation (genetic correlations between −0.215 and 0.251). Peak milk yield exerted a moderate and unfavorable effect on the interval from parturition to first service (genetic correlations of 0.414 and 0.353 after first and second calving, respectively), and a low and unfavorable effect on other fertility traits (genetic correlations between −0.281 and 0.295). Infertility after first calving caused a strong elongation of the lactation, and LL was negatively correlated with fertility of cows after second calving, so that LL can itself be regarded as a measure of fertility. Lactation milk yield depends on both pMY and LL, and, as such, is a cause and consequence of (in)fertility.     

Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds

The objective of this study was to quantify the genetic variation in normal and atypical progesterone profiles and investigate if this information could be useful in an improved genetic evaluation for fertility for dairy cows. The phenotypes derived from normal profiles included cycle length traits, including commencement of luteal activity (C-LA), interluteal interval, luteal phase length. and interovulatory interval. In total, 44,977 progesterone test-day records were available from 1,612 lactations on 1,122 primiparous and multiparous Holstein-Friesian cows from Ireland, the Netherlands, Sweden, and the United Kingdom. The atypical progesterone profiles studied were delayed cyclicity, prolonged luteal phase, and cessation of cyclicity. Variance components for the atypical progesterone profiles were estimated using a sire linear mixed model, whereas an animal linear mixed model was used to estimate variance components for the cycle length traits. Heritability was moderate for delayed cyclicity (0.24 ± 0.05) and C-LA (0.18 ± 0.04) but low for prolonged luteal phase (0.02 ± 0.04), luteal phase length (0.08 ± 0.05), interluteal interval (0.08 ± 0.14), and interovulatory interval (0.03 ± 0.04). No genetic variation was detected for cessation of cyclicity. Commencement of luteal activity, luteal phase length, and interovulatory interval were moderately to strongly genetically correlated with days from calving to first service (0.35 ± 0.12, 0.25 ± 0.14, and 0.76 ± 0.24, respectively). Delayed cyclicity and C-LA are traits that can be important in both genetic evaluations and management of fertility to detect (earlier) cows at risk of compromised fertility. Delayed cyclicity and C-LA were both strongly genetically correlated with milk yield in early lactation (0.57 ± 0.14 and 0.45 ± 0.09, respectively), which may imply deterioration in these traits with selection for greater milk yield without cognizance of other traits.     

Genetic correlations between fertility traits and milk composition and fatty acids in Holstein-Friesian,Brown Swiss,and Simmental cattle using recursive models

This study aimed to investigate the genetic and putative causal relationships between fertility traits [i.e., days open (DO) and calving rate (CR)] and milk quality, composition, and fatty acid contents (milk composition traits) in Holstein-Friesian, Brown Swiss, and Simmental cattle, using recursive models within a Bayesian framework. Trivariate animal models were run, each including one fertility trait, one milk composition trait, and, in all models, milk yield. The DO and CR data were merged with the test days closest to the insemination date for milk composition traits. After editing, 16,468 to 23,424 records for Holstein-Friesian, 23,424 to 46,660 for Brown Swiss, and 26,105 to 35,574 for Simmental were available for the analyses. Recursive animal models were applied to investigate the possible causal influences of milk composition traits on fertility and the genetic relationships among these traits. The results suggested a potential cause-and-effect relationship between milk composition traits and fertility traits, with the first trait influencing the latter. We also found greater recursive effects of milk composition traits on DO than on CR, the latter with some putative differences among breeds in terms of sensitivity. For instance, the putative causal effects of somatic cell score on CR (on the observed scale, %) varied from ?0.96 to ?1.39%, depending on the breed. Concerning fatty acids, we found relevant putative effects of C18:0 on CR, with estimates varying from ?7.8 to ?9.9%. Protein and casein percentages, and short-chain fatty acid showed larger recursive effects on CR, whereas fat, protein, and casein percentages, C16:0, C18:0, and long-chain fatty acid had larger effects on DO. The results obtained suggested that these milk traits could be considered as effective indicators of the effects of animal metabolic and physiological status on fertility.     

Evaluation of bull fertility in Italian Brown Swiss dairy cattle using cow field data

Dairy bull fertility is traditionally evaluated using semen production and quality traits; however, these attributes explain only part of the differences observed in fertility among bulls. Alternatively, bull fertility can be directly evaluated using cow field data. The main objective of this study was to investigate bull fertility in the Italian Brown Swiss dairy cattle population using confirmed pregnancy records. The data set included a total of 397,926 breeding records from 1,228 bulls and 129,858 lactating cows between first and fifth lactation from 2000 to 2019. We first evaluated cow pregnancy success, including factors related to the bull under evaluation, such as bull age, bull inbreeding, and AI organization, and factors associated with the cow that receives the dose of semen, including herd-year-season, cow age, parity, and milk yield. We then estimated sire conception rate using only factors related to the bull. Model predictive ability was evaluated using 10-fold cross-validation with 10 replicates. Interestingly, our analyses revealed that there is a substantial variation in conception rate among Brown Swiss bulls, with more than 20% conception rate difference between high-fertility and low-fertility bulls. We also showed that the prediction of bull fertility is feasible as our cross-validation analyses achieved predictive correlations equal to 0.30 for sire conception rate. Improving reproduction performance is one of the major challenges of the dairy industry worldwide, and for this, it is essential to have accurate predictions of service sire fertility. This study represents the foundation for the development of novel tools that will allow dairy producers, breeders, and artificial insemination companies to make enhanced management and selection decisions on Brown Swiss male fertility.     

Commencement of luteal activity in three different selection lines for milk yield and fertility in Norwegian Red cows

Relationships among commencement of luteal activity (C-LA), milk yield, and energy balance (EB) were investigated in 3 selection lines of Norwegian Red cows at the Norwegian University of Life Sciences from 1994 through 2001. The cows were selected for low genetic merit for milk yield (LMP), high genetic merit for milk yield (HMP), and a combination of high indices for milk yield and fertility (HI). Breeding values for fertility were based on 56-d nonreturn rate. The material included 268 lactations from 147 cows. Milk samples for progesterone analysis were drawn 3 times weekly from 1994 through 1998, and 2 times weekly from 1999 to 2001. Commencement of luteal activity was defined as the first 2 consecutive measurements of progesterone concentration >3 ng/mL not earlier than 10 d after calving. Selection line was significantly related to C-LA, so that the least squares mean days from calving to C-LA were 22.5, 30.4, and 27.2 d for LMP, HMP, and HI cows, respectively. The HMP cows produced more milk than the LMP cows. The average milk yield in the sixth week of lactation was 24.0, 27.1, and 25.3 kg for LMP, HMP, and HI cows, respectively. The interval to C-LA decreased for the HMP and HI cows after phenotypic adjustment for EB in the model. Least squares means for the interval to C-LA were 23.2, 29.7, and 25.6 d for the LMP, HMP, and HI cows, respectively, in a model that included parity, selection lines, and EB as covariates. Cumulated EB during the first 4 wk of lactation, which itself differed between selection lines, did not fully account for differences in interval to C-LA between selection lines. Thus, the results of the present investigation indicate that selection for milk yield negatively affects C-LA over and above the effects caused by concurrent changes in EB. The increase in days to C-LA caused by selection for high yields can be reduced if selection for milk yield is combined with fertility in the breeding program.     

Genetic evaluation of fertility traits of dairy cattle using a multiple-trait animal model

A genetic evaluation system was developed for 5 fertility traits of dairy cattle: interval from first to successful insemination and nonreturn rate to 56 d of heifers, and interval from calving to first insemination, nonreturn rate to 56 d, and interval first to successful insemination of cows. Using the 2 interval traits of cows as components, breeding values for days open were derived. A multiple-trait animal model was applied to evaluate these fertility traits. Fertility traits of later lactations of cows were treated as repeated measurements. Genetic parameters were estimated by REML. Mixed model equations of the genetic evaluation model were solved with preconditioned conjugate gradients or the Gauss-Seidel algorithm and iteration on data techniques. Reliabilities of estimated breeding values were approximated with a multi-trait effective daughter contribution method. Daughter yield deviations and associated effective daughter contributions were calculated with a multiple trait approach. The genetic evaluation software was applied to the insemination data of dairy cattle breeds in Germany, Austria, and Luxembourg, and it was validated with various statistical methods. Genetic trends were validated. Small heritability estimates were obtained for all the fertility traits, ranging from 1% for nonreturn rate of heifers to 4% for interval calving to first insemination. Genetic and environmental correlations were low to moderate among the traits. Notably, unfavorable genetic trends were obtained in all the fertility traits. Moderate to high correlations were found between daughter yield-deviations and estimated breeding values (EBV) for Holstein bulls. Because of much lower heritabilities of the fertility traits, the correlations of daughter yield deviations with EBV were significantly lower than those from production traits and lower than the correlations from type traits and longevity. Fertility EBV were correlated unfavorably with EBV of milk production traits but favorably with udder health and longevity. Integrating fertility traits into a total merit selection index can halt or reverse the decline of fertility and improve the longevity of dairy cattle.     

Genetic analysis of days from calving to first insemination and days open in Danish Holsteins using different models and censoring scenarios

The objectives of this study were to estimate genetic parameters and evaluate models for genetic evaluation of days from calving to first insemination (ICF) and days open (DO). Data including 509,512 first-parity records of Danish Holstein cows were analyzed using 5 alternative sire models that dealt with censored records in different ways: 1) a conventional linear model (LM) in which a penalty of 21 d was added to censored records; 2) a bivariate threshold-linear model (TLM), which included a threshold model for censoring status (0, 1) of the observations, and a linear model for ICF or DO without any penalty on censored records; 3) a right-censored linear model (CLM); 4) a Weibull proportional hazard model (SMW); and 5) a Cox proportional hazard model (SMC) constructed with piecewise constant baseline hazard function. The variance components for ICF and DO estimated from LM and TLM were similar, whereas CLM gave higher estimates of both additive genetic and residual components. Estimates of heritability from models LM, TLM, and CLM were very similar (0.102 to 0.108 for ICF, and 0.066 to 0.069 for DO). Heritabilities estimated using model SMW were 0.213 for ICF and 0.121 for DO in logarithmic scale. Using SMC, the estimates of heritability, defined as the log-hazard proportional factor for ICF and DO, were 0.013 and 0.009, respectively. Correlations between predicted transmitting ability from different models for sires with records from at least 20 daughters were far from unity, indicating that different models could lead to different rankings. The largest reranking was found between SMW and SMC, whereas negligible reranking was found among LM, TLM, and CLM. The 5 models were evaluated by comparing correlations between predicted transmitting ability from different data sets (the whole data set and 2 subsets, each containing half of the whole data set), for sires with records from at least 20 daughters, and χ² statistics based on predicted and observed daughter frequencies using a cross validation. The model comparisons showed that SMC had the best performance in predicting breeding values of the 2 traits. No significant difference was found among models LM, TLM, and CLM. The SMW model had a relatively poor performance, probably because the data are far from a Weibull distribution. The results from the present study suggest that SMC could be a good alternative for predicting breeding values of ICF and DO in the Danish Holstein population.     

Random regression models for male and female fertility evaluation using longitudinal binary data

A longitudinal Bayesian threshold analysis of insemination outcomes was carried out using 2 random regression models with 3 (Model 1) and 5 (Model 2) parameters to model the additive genetic values at the liability scale. All insemination events of first-parity Holstein cows were used. The outcome of an insemination event was treated as a binary response of either a success (1) or a failure (0). Thus, all breeding information for a cow, including all service sires, was included, thereby allowing for a joint evaluation of male and female fertility. An edited data set of 369,353 insemination records from 210,373 first-lactation cows was used. On the liability scale, both models included the systematic effects of herd-year, month of insemination, technician, and regressions on age of service sire and milk yield during the first 100 d of lactation. The random effects in the model were the 3 or 5 random regression coefficients specific to each cow, the permanent effect of the cow, and the service sire effect. Using Model 1, the estimated heritability of an insemination outcome decreased from 0.035 at d 50 to 0.032 at d 140 and then increased continuously with DIM. The genetic correlations for insemination success at different time points ranged from 0.83 to 0.99, and their magnitude decreased with an increase in the interval between inseminations. A similar trend was observed for heritability and genetic correlations using Model 2. However, the average estimate of heritability was much higher (0.058) than those obtained using Model 1 or a repeatability model. In addition, the estimated genetic correlations followed the same trend as Model 1, but were lower and with a higher rate of decrease when the interval between inseminations increased. The posterior mean of service sire variance was 0.01 for both models, and permanent environmental variance was 0.05 and 0.02 for Models 1 and 2, respectively. Model comparison based on the Bayes factor indicated that Model 1 was more plausible, given the data.     

Genetic analysis of new progesterone-based fertility traits in dairy cows measured on-farm

The objectives of this study were (1) to analyze the agreement of a standard laboratory ELISA for progesterone (P4) with an automated on-farm ELISA kit operated under commercial conditions in 1,297 milk samples from 50 dairy cows; (2) to study the influence of the method of detection of luteal activity on genetic parameters of fertility traits based on P4 measured with an automated on-farm ELISA once weekly from wk 3 to 9 postpartum in the milk of 1,304 cows; and (3) to study the influence of sampling frequency (once or twice weekly from wk 3 to 9) on the same traits from 296 cows. Luteal activity can be detected when there is an active corpus luteum in the ovary producing P4 and indicating the onset of reproductive cyclicity after calving. The on-farm ELISA overestimated P4 contents by a mean square error of prediction of 2.76 ng/mL and had an intermediate Spearman correlation with the laboratory kit (0.54). For the second objective, the postpartum interval to the commencement of luteal activity (C-LA), proportion of luteal activity between d 15 and 63 postpartum (P-LA), calculated as the number of samples above the threshold for high P4 values divided by the number of all samples, and delay of first ovulation (DOV1), defined as C-LA occurring later than d 45 postpartum, were derived from the P4 profiles. Both C-LA and DOV1 were determined by (a) thorough qualitative visual inspection of the profile, (b) the profile's mean as threshold for the first increase in P4 postpartum, indicating commencement of luteal activity, and (c) 3 ng/mL as threshold for the first increase in P4, a value that has been used by many other studies. Similarly, P-LA was determined by using methods (b) and (c). Estimates of heritability were 0.04 to 0.13 for C-LA, 0.12 to 0.23 for P-LA, and 0.03 to 0.07 for DOV1. Genetic correlation of P-LA with C-LA and with the profile's mean P4 was ?1.00. The profile's mean had a higher estimate of heritability (0.11–0.12) than C-LA or DOV1. It can be calculated as the arithmetic mean of all P4 values of a profile, whereas C-LA, P-LA, and DOV1 need a definition of a threshold for high P4 values. We therefore suggest the profile's mean as a promising candidate for further research. For the third objective, once-weekly sampling was mimicked by neglecting every second sample, and C-LA and DOV1 shifted toward a later onset of cyclicity. Thus, a common standard for sampling regimen and detection algorithm is essential to avoid incompatibility between studies.     

Genetic parameters for digital dermatitis and correlations with locomotion, production, fertility traits, and longevity in Holstein-Friesian dairy cows

Heritability of digital dermatitis (DD) and correlations between DD and type traits related to legs and feet were estimated from a linear animal model. Data comprised 93,391 national type evaluation records of pedigreed first-lactation Holstein-Friesian cows that calved from 2002 through 2006. At the time of classification, cows were housed in different housing systems (i.e., cubicles, straw yards, slatted or loafing yards) and on pasture. The type traits evaluated were locomotion score (LOCO), rear legs side view (RLS), foot angle (FA), bone quality and leg and feet composite (L&F). In addition, cows were examined for DD lesions at classification. The relationships among these type traits, lifespan (LS), production (milk and fat), fertility (calving interval and 56-d nonreturn) and DD were examined by estimating the approximate genetic correlations from sire estimated breeding values. The study also evaluated the association between DD and the housing systems as well as the general conditions of the farm flooring where classification took place. In general, cows on pasture were less susceptible to DD than cows in other housing systems, whereas the association between DD and the flooring conditions was counterintuitive. Heritability estimate for DD was 0.011 on the 0/1 scale, which is equivalent to 0.029 on the assumed underlying normally distributed scale. Bone quality, LOCO, and L&F had moderate to high negative genetic correlations with DD, indicating that flatter, more refined bones, higher LOCO, and better L&F were associated with less incidence of DD. The genetic correlations between DD, RLS, and FA were not significantly different from zero. Digital dermatitis had moderate but negative genetic correlations with LS and milk and fat, suggesting that breeding for resistance to DD will result in an increase in both longevity and production. Cows affected with DD had a slightly shorter calving interval than healthy cows, an association found to be mediated through the reduced milk yield of these cows. Generally, the type traits included in this study had low genetic correlations with production and fertility traits whereas the associations between these traits and LS ranged from moderate to high. This indicates that good locomotion, straighter RLS, steeper FA, better L&F, and flatter, more refined bones are associated with increased longevity.     

Bayesian inference of the inbreeding load variance for fertility traits in Brown Swiss cattle

Our study investigated the inbreeding load for fertility traits in the Italian Brown Swiss dairy cattle breed. Fertility traits included continuous traits (i.e., interval from calving to first service, days open, and calving interval) and categorical traits (i.e., calving rate at first insemination and nonreturn date at d 56). We included only records of the first 3 parities of cows that calved between 2010 and 2018. We traced up the pedigree of the cows with records as far as possible, ending up with a total of 73,246 animals. The final data set consisted of 59,864 records from 34,921 cows. We analyzed all models using a Bayesian approach that included a covariate with total inbreeding in addition to systematic, permanent environment, additive genetic, and inbreeding load effects. We then evaluated the trends in heritabilities and ratios of the inbreeding load using a continuum of partial inbreeding coefficients from 0.001 to 0.100 as reference. Posterior estimates of heritabilities tended to decrease across the continuum, whereas ratios of the inbreeding load tended to increase, more noticeably in categorical traits (calving rate at first insemination and nonreturn date at d 56). From the results obtained, we confirmed the presence of heterogeneity in inbreeding depression. We then predicted the inbreeding load effects, which had a low reliability of prediction, explained by having only 513 ancestors generating inbreeding. However, reliability of prediction was high enough for some of the individuals, obtaining a favorable prediction of inbreeding load for a relevant percentage, which improved the phenotypic performance of their inbred descendants. These results make it feasible to implement breeding and management strategies that select ancestors with a favorable inbreeding load prediction. In addition, it opens the possibility to define a global index for the expected consequences of the inbreeding generated by each individual.     

Genetic analysis of male and female fertility using longitudinal binary data

A longitudinal Bayesian threshold analysis of insemination events during the first 250 d after calving of first-parity Holsteins was carried out. The outcome of an insemination event was treated as a binary response of either a success (1) or a failure (0). Thus, all breeding information for a cow, including all service sires, was included, thereby allowing for a joint evaluation of male and female fertility. An edited data set of 297,823 insemination records from 151,758 first lactation cows was used. On the liability scale, the model included the systematic effects of herd-year of insemination, technician, month of insemination, and regressions on age of service sire, 3 test days in the first 100 d of lactation (early milk yield), and days in milk at insemination. The random effects in the model were the additive breeding value, the permanent effect of the cow, and the service sire effect. Posterior mean (standard deviation) of the dispersion parameters in the model were 0.034 (0.006), 0.009 (0.001), and 0.171 (0.013) for the additive, service sire, and permanent environmental variances, respectively. The residual variance was fixed at 1, as a result of the nonidentifiability of the threshold model. The posterior mean (standard deviation) of heritability was 0.028 (0.005). This point estimate of heritability is well within the range of available estimates for the trait. Thus, these estimates suggest that some genetic variation exists that can potentially be used to improve reproductive performance or at least avoid its further deterioration. The estimate of the regression coefficient on age of service sire was 0.001, indicating better fertility among older bulls. However, this result has to be interpreted with caution given the preferential use of proven bulls on well-managed cows (as opposed to problem breeders). The estimate of the regression coefficient was negative (-0.005) for early milk yield, as expected, and positive (0.003) for days in milk at insemination. This suggests that high-producing cows are less likely to conceive at the beginning of lactation.     

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.     

Estimates of genetic parameters for Canadian Holstein female reproduction traits

Age at first insemination, days from calving to first insemination, number of services, first-service nonreturn rate to 56 d, days from first service to conception, calving ease, stillbirth, gestation length, and calf size of Canadian Holstein cows were jointly analyzed in a linear multiple-trait model. Traits covered a wide spectrum of aspects related to reproductive performance of dairy cows. Other frequently used fertility characteristics, like days open or calving intervals, could easily be derived from the analyzed traits. Data included 94,250 records in parities 1 to 6 on 53,158 cows from Ontario and Quebec, born in the years 1997 to 2002. Reproductive characteristics of heifers and cows were treated as different but genetically correlated traits that gave 16 total traits in the analysis. Repeated records for later parities were modeled with permanent environmental effects. Direct and maternal genetic effects were included in linear models for traits related to calving performance. Bayesian methods with Gibbs sampling were used to estimate covariance components of the model and respective genetic parameters. Estimates of heritabilities for fertility traits were low, from 3% for nonreturn rate in heifers to 13% for age at first service. Interval traits had higher heritabilities than binary or categorical traits. Service sire, sire of calf, and artificial insemination technician were important (relative to additive genetic) sources of variation for nonreturn rate and traits related to calving performance. Fertility traits in heifers and older cows were not the same genetically (genetic correlations in general were smaller than 0.9). Genetic correlations (both direct and maternal) among traits indicated that different traits measured different aspects of reproductive performance of a dairy cow. These traits could be used jointly in a fertility index to allow for selection for better fertility of dairy cattle.     

Genetic and genomic analysis of hyperthelia in Brown Swiss cattle

Supernumerary teats (SNT) are any abnormal teats found on a calf in addition to the usual and functional 4 teats. The presence of SNT has also been termed “hyperthelia” since the end of the 19th century. Supernumerary teats can act as an incubator for bacteria, infecting the whole udder, and can interfere with the positioning of the milking machine, and consequently, have economic relevance. Different types of SNT are observed at different positions on the udder. Caudal teats are in the rear, ramal teats are attached to another teat, and intercalary teats are found between 2 regular teats. Not all teats are equally developed; some are completely functional but most are rudimentary and not attached to any mammary gland tissue. Recently, different studies showed the poly/oligogenic character of these malformations in cattle as well as in other mammalian species. The objective of this study was to analyze the genetic architecture and incidence of hyperthelia in Swiss Brown Swiss cattle using both traditional genetic evaluation as well as imputed whole genome sequence variant information. First, phenotypes collected over the last 20 yr were used together with pedigree information for estimation of genetic variance. Second, breeding values of Brown Swiss bulls were estimated applying the BLUP algorithm. The BLUP-EBV were deregressed and used as phenotypes in genome-wide association studies. The gene LGR5 on chromosome 5 was identified as a candidate for the presence of SNT. Using alternative trait coding, genomic regions on chromosome 17 and 20 were also identified as being involved in the development of SNT with their own supernumerary mammary gland tissue. Implementing knowledge gained in this study as a routine application allows a more accurate evaluation of the trait and reduction of SNT prevalence in the Swiss Brown Swiss cattle population.     

Analysis of fertility and dystocia in Holsteins using recursive models to handle censored and categorical data

A method based on the analysis of recursive multiple-trait models was used to 1) estimate genetic and phenotypic relationships of calving ease (CE) with fertility traits and 2) analyze whether dystocia negatively affects reproductive performance in the next reproductive cycle. Data were collected from 1995 through 2002, and contained 33,532 records of CE and reproductive data of 17,558 Holstein cows distributed across 560 herds in official milk recording from the Basque Country Autonomous Community (Spain). The following fertility traits were considered: days open (DO), days to first service, number of services per pregnancy (NINS), and outcome of first insemination (OFI). Four bivariate sire and sire-maternal grandsire models were used for the analyses. Censoring existed in DO (26.49% of the data) and NINS (12.22% of the data) because of cows having been sold or culled before reaching the next parturition. To avoid bias, a data augmentation technique was applied to censored data. Threshold models were used for CE and OFI. To consider that CE affects fertility and the genetic determination of CE and fertility traits, recursive models were applied, which simultaneously considered CE as a fixed effect on fertility performance and the existence of a genetic correlation between CE and fertility traits. The effects of CE score 3 (difficult birth) with respect to score 1 (no problem) for days to first service, DO, NINS, and OFI were 8 d, 31 d, 0.5 services, and - 12% success at first insemination, respectively. These results showed poorer fertility after dystocia. Genetic correlations between genetic effects of fertility traits and CE were close to zero, except for the genetic correlations between direct effects of DO and CE, which were positive, moderate, and statistically different from 0 (0.47 ± 0.24), showing that genes associated with difficult births also reduce reproductive success.     

Atypical progesterone profiles and fertility in Swedish dairy cows

The incidence of normal and atypical progesterone profiles in Swedish dairy cows was studied. Data were collected from an experimental herd over 15 yr, and included 1,049 postpartum periods from 183 Swedish Holstein and 326 Swedish Red and White dairy cows. Milk progesterone samples were taken twice weekly until initiation of cyclical ovarian activity and less frequently thereafter. Progesterone profiles were 1) normal profile: first rise in milk progesterone above the threshold value before d 56 postpartum, followed by regular cyclical ovarian activity (70.4%); 2) delayed onset of cyclical ovarian activity: low milk progesterone the first 56 d postpartum (15.6%); 3) cessation of cyclical ovarian activity: ovarian activity resumed within 56 d postpartum, but ceased for a period of 14 d or more (6.6%); and 4) prolonged luteal phase: ovarian activity resumed within 56 d postpartum, but milk progesterone remained elevated in the nonpregnant cow for a period of 20 d or more (7.3%). Swedish Holsteins had 1.5 times higher risk of atypical profile than Swedish Red and Whites. Risk of atypical profiles was 0.5 and 0.7 times lower for older cows compared with first-parity cows; 2.3 times higher for cows in tie-stalls compared with those in loose housing; 2.6 times higher for cows calving during winter compared with summer; 0.5 times lower for cows in earlier (1994-1999) calving-year groups compared with the most recent (2000-2002); 2.5 times higher for cows with planned extended calving interval compared with conventional calving interval; and 2.2 times higher for an atypical profile in previous lactation compared with a normal profile. Cows with atypical profiles had a 15-d increase in interval from calving to first artificial insemination and an 18-d increase in interval from calving to conception. Progesterone samples taken within the first 60 d postpartum were used to calculate the percentage of samples above the threshold value of luteal activity. This measure had a significantly different mean in profiles and can be used to separate delayed onset of cyclical ovarian activity profiles and prolonged luteal phase profiles from normal. Thereby, it may be a more effective tool than measurements based only on the onset of ovarian cyclical activity in genetic evaluation of early postpartum fertility in dairy cows.     

Genetic analysis of predicted fatty acid profiles of milk from Danish Holstein and Danish Jersey cattle populations

The objective of this study was to assess the genetic variability of the detailed fatty acid (FA) profiles of Danish Holstein (DH) and Danish Jersey (DJ) cattle populations. We estimated genetic parameters for 11 FA or groups of FA in milk samples from the Danish milk control system between May 2015 and October 2016. Concentrations of different FA and FA groups in milk samples were measured by mid-infrared spectroscopy. Data used for parameter estimation were from 132,732 first-parity DH cows and 21,966 first-parity DJ cows. We found the highest heritabilities for test day measurements in both populations for short-chain FA (DH = 0.16; DJ = 0.16) and C16:0 (DH = 0.14; DJ = 0.16). In DH, the highest heritabilities were also found for saturated FA and monounsaturated FA (both populations: 0.15). Genetic correlations between the fatty acid traits showed large differences between DH and DJ for especially short-chain FA with the other FA traits measured. Furthermore, genetic correlations of total fat with monounsaturated FA, polyunsaturated FA, short-chain FA, and C16:0 differed markedly between DH and DJ populations. In conclusion, we found genetic variation in the mid-infrared spectroscopy-predicted FA and FA groups of the DH and DJ cattle populations. This finding opens the possibility of using genetic selection to change the FA profiles of dairy cattle.