Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle
© Eaglen et al.; licensee BioMed Central Ltd. 2012
Received: 26 October 2011
Accepted: 13 July 2012
Published: 28 July 2012
The focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle.
Data from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (−maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models.
Results and discussion
On average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (−maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (−maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence.
For the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (−maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.
Calving is a key event on a dairy cattle farm and successful calvings are important to financial success of the farm. Calving complications lead to an increase in veterinary and labour costs and a decrease in revenue (loss of animals and/or reduced subsequent performance) [1–3]. Furthermore, difficult calvings compromise animal welfare and thereby consumer acceptability of dairy production systems .
In recent years, dairy cattle breeders have shown an increasing interest in selection for functional traits  and gradually the focus of selection is shifting from traits that increase profit towards traits that reduce costs . For example, easy parturition and calf viability are economically important traits but are not classical production traits. Since genetic selection could improve calving performance, it is important to include calving traits in genetic evaluations, although their implementation is not straightforward.
Calving ease (CE) and stillbirth (SB) are phenotypes that are generally scored on categorical or binary scales which make them sensitive to subjectivity, especially CE . Furthermore, heritabilities of these traits are generally low  and thus much data is needed to obtain sufficiently accurate estimates that have an impact on selection indices. But above all, calving trait phenotypes are affected by two separate components, the calf’s contribution (direct effect; e.g. arising from size, hormonal balance, weight) and the dam’s contribution (maternal effect; e.g. arising from pelvic measurements, ability to respond to parturition signalling), and this complicates statistical analysis [8–10]. In quantitative genetics, the objective is to separate additive genetic variances and co-variances from other sources of variance. The statistical model fitted for calving traits should therefore allow the separation and estimation of both direct and maternal effects but there is no consensus on which is the most accurate model to achieve this objective. Various statistical models have been reported throughout the years, each aimed at improving one aspect of their predictive ability. Consequently, many different statistical models are used in routine genetic evaluations to estimate genetic parameters for calving traits , while for production traits there is a greater consensus across countries . For CE, statistical models range from sire-maternal grandsire (S-MGS) models to animal models in univariate or multi-trait form that either allow a direct-maternal genetic covariance or fix this covariance to zero .
Overall, statistical models can account for direct and maternal effects in two ways i.e. animal models that fit calf and dam effects directly and S-MGS models that fit direct and maternal effects through the sire of the calf and dam, respectively . Each of these then has a family of implementation depending on how traits are defined and modelled e.g. univariately or multi-trait, which leads to the divergence in models seen today. The genetic covariance between direct and maternal effects plays a key role in the interpretation of estimated genetic parameters and the prediction of response to selection. For example, in the case of CE, this parameter represents the genetic relationship between the animals’ genetic predisposition for ease of being born and ease of giving birth.
This genetic covariance is however suggested to be sensitive to estimation bias [12–14]. Thus, estimates of direct-maternal genetic correlations reported for calving traits, which are often negative and of moderate to high magnitude, are questioned . Theories on the source of bias in estimates of direct-maternal genetic correlations range from ignored environmental covariances [14, 16, 17] to unaccounted contemporary groups [15, 18, 19]. Estimation of the direct-maternal covariance and remaining genetic components is said to be further improved by using a multi-trait model instead of a univariate model [20, 21]. Since CE and SB are highly correlated and show low heritabilities, it has been suggested that a multi-trait model that incorporates a highly heritable and correlated indicator trait such as birth weight, calf size and/or gestation length (GL), would lead to a more optimal analysis [20–22]. In addition to models incorporating indicator traits, it has also been proposed that the extension of univariate models to a multi-trait model between parities would be useful because of the potential genetic distinctiveness of calving traits in first and later parities. This implies that models fitting first and later parities as correlated traits are theoretically more correct than models that fit parity as a non-genetic effect [23–25].
However, given the high standard errors of the estimated maternal variances and direct-maternal covariances, it appears that statistical models used to analyse calving traits can still benefit from optimization. Previous studies on the estimation of genetic parameters for calving traits have been comparing models. However, in most cases, studies limit themselves to single comparisons, such as univariate versus multi trait models within either a sire or animal model framework [21, 24]. The objective of this study was to estimate the genetic parameters of calving traits for United Kingdom (UK) Holstein-Friesian cattle with a range of present-day statistical models for national genetic evaluations of calving traits. A comparison was made between S-MGS versus animal models and univariate versus multi-trait models between traits and parities.
GL was added to the multi-trait model to study any benefits of this indicator trait to the predictive ability of the model. Specific attention was given to the potential bias of the sensitive genetic correlation between direct and maternal effects and, lastly, a ‘total heritable variance’ for calving traits was estimated, combining direct and maternal variance components. This parameter was recently reported for the analysis of traits affected by indirect genetic effects  and is suggested to represent the total amount of total genetic variance available for response to selection. Its practical use for calving traits is discussed.
Data on calving in Holstein Friesian cattle was provided by two milk recording organisations (MRO) in the UK i.e. the Cattle Information Service (CIS) and National Milk Records (NMR). GL was calculated from insemination and calving dates and restricted to 265–295 days. Only records on single born calves were used. Datasets from first and later parity records contained 43 135 and 330 753 records respectively, spanning the years 1995 to 2009. This data was checked for inconsistencies in parity, breed, calving date and age at calving. Incorrect records were discarded, as were contemporary groups with fewer than two records (herd, herd-year, sire and maternal grandsire). Age at calving was restricted to 18–48 months for heifers and 30–70 months for cows. To avoid herds where farmers unrealistically recorded the same CE score for most or all contemporary groups of calvings, the standard deviation of CE score within herd-year was calculated. Herd-year classes with standard deviations of zero were deleted when this was considered statistically improbable (97.5% confidence interval) given the herd size and CE score distribution. Sex of stillborn calves was not recorded by one of the data sources. As sex has been shown to have a considerable effect on SB  data from this specific data source was deleted for the univariate and bivariate between-parity analyses on SB. For all univariate analyses on CE and GL, the data from this source for a stillborn calf was set as missing but the remaining data obtained by the specific data source was used. The final dataset consisted of 30 640 first parity records originating from 2 098 herds representing 2 012 (service) sires and 4 783 maternal grandsires (MGS). The accompanying pedigree consisted of ~ 200 000 individuals and was 5 generations deep. Since the later parity dataset created computational problems when fitting the animal model, it was reduced by only retaining records that were of 2nd and 3rd parity and had values for each of the three traits (prior to edits on SB data described earlier). Furthermore, the thresholds for the minimum size of herd-year and sire contemporary groups were increased to 7 and 5 records, respectively. Finally, the final later parity dataset consisted of 54 744 records, originating from 2 108 herds, 1 918 sires and 5 886 maternal grandsires, with an accompanying pedigree of ~ 290 000 individuals and 4 generations deep.
Descriptive statistics of the data
Number of records
Number of dams
Number of dams with own birth record1
Number of sires
Number of maternal grandsires
GL (days), mean ± s.e.3
281.35 ± 4.89
281.17 ± 4.93
Direct and maternal genetic effects
Thus, during the life of i, the direct additive effect (Ad,i) is expressed at the start of life while, if i is a female, the maternal additive effect of i (Am,i) is expressed whenever she calves. ED,i and EM,j are the direct environmental effect, property of calf I, and maternal environmental effect, property of dam j, respectively.
In all cases, optional fixed effects and potential interaction effects were tested for significance in SAS V9.1 (P < 0.05)  and then the variance components were calculated using REML, with ASREML version 3.0 . Sex by parity and sex by age interaction effects were not significant. Prior to using multi-trait models involving GL, the relationship of GL with CE and SB was examined according to Hansen et al.  to ensure that the use of traditional bivariate models was appropriate given their assumption of linearity.
Furthermore, the limitations of analysing categorical traits, such as CE and SB, with linear models are well known. Multiple model assumptions are violated due to the fact that values of categorical data are bounded within certain limits e.g. 0 to 1 or 1 to 4. Therefore, generalized linear mixed models (GLMM) such as threshold models can be more appropriate for the analysis of categorical traits since scores are transformed by the model into values on an underlying continuous liability scale. This study aimed to evaluate statistical models with the ultimate goal of implementation in national genetic evaluations of CE, SB and GL. Although threshold models are implemented in routine national genetic evaluations in France and the USA [29, 30], in the UK and most other countries, calving traits are evaluated with a linear model. Therefore, we chose to evaluate several modelling possibilities within linear rather than threshold models. In the discussion section of this paper, we will elaborate further on this choice.
Univariate animal and S-MGS models
Bivariate models between traits
Residuals, e, and permanent environmental effects, pepe, were assumed to be MVN(0, Reσ e 2), and MVN(0, Rpeσ pe 2), where Re and Rpe denote the residual and permanent environmental 2 × 2 variance covariance matrices and σ e 2and σ pe 2 were the residual variance and permanent environmental variance.
Bivariate models between parities
To study the genetic correlation between calving traits in first and later parities, bivariate S-MGS models were fitted with first and later parities (2nd and 3rd parities combined) treated as correlated traits. The model equalled equation 4, with yi a vector representing the observations for each trait in first (y1) and later parities (y2). Random genetic effects were fitted for the sire and maternal grandsire. The fixed and random non-genetic effects were the same as in the univariate animal model. Estimates of sire and maternal grandsire variances were transformed into direct and maternal effects according to equation 3.
where σ Ad 2 and σ Am 2 are the direct and maternal additive genetic variances, σ Adm is the additive direct maternal covariance and σ e 2 is the environmental variance. To allow easy comparison with other studies, herd-year variances and permanent environmental variances were not included in the phenotypic variance but are provided in Additional file 1: Tables S1 and S2. The heritabilities and genetic direct-maternal correlations were estimated more than once by the several bivariate models and these were pooled in meta- analyses according to Corbin et al. 
Direct-maternal genetic covariance
A negative direct-maternal relationship would be worrying for the dairy cattle industry since it suggests that selecting a sire that is genetically superior for ease of birth may later cause a problem when its daughters calve. Koch  showed that, when ignored or assumed to be zero, a direct-maternal environmental covariance (cov(Ed,i, Em,i)) can cause bias in the estimated genetic parameters. Although it is possible to fit a correlation structure in the residual to avoid this problem, computational complexity is then substantially increased. Therefore, in this study, we chose to avoid this potential bias by removing from the data all individuals that appeared as both calf and dam. Residuals of these specific records would otherwise be correlated . Then to evaluate the bias, animals were reintroduced and the analyses were repeated. Throughout the paper, animal model 1 (A1) represents the animal model which was used to analyse the reduced data, whereas animal model 2 (A2) represents the animal model used to analyse the total data.
Total heritable variance
In this context, the σ TBV 2 represents the total genetic variance available for response to selection, with response predicted by R = ιρ M σ TBV 2 where ι is the selection intensity and ρ M is the accuracy of selection . This is distinct from the total heritable variance reported by Willham , Meyer  and Koch , which refers to mass selection, as explained in Eaglen et al. . Given the current selection strategies based on PTA in dairy cattle, we estimated and explored σ TBV 2 as described by Bijma et al. .
Results and discussion
Table 1 presents the descriptive statistics of the data. It shows that in the UK, approximately 20% of the calvings required assistance of some sort. Incidence of calving assistance was higher in first than in later-parity calvings, which agrees with the general consensus that calving complications are of more concern in heifers than in cows . Moreover, severe calving difficulty was experienced by approximately 4% and 2% of heifers and cows respectively. These are in line with international prevalences of calving difficulty in the Holstein breed  although comparison is not straightforward since the scoring system of CE allows for a large variety of score definitions . The incidence of SB in first and later parities (Table 1) agrees with incidences reported by Hansen et al.  and Jamrozik et al. . Table 1 also shows that there were fewer males than females in the first-parity dataset, which could indicate a bias in data recording due to the difference in value between a bull and a heifer calf in dairy cattle. Since the calving of bull calves is known to be more difficult , it is possible that CE is under-reported. However, all studies using field records for CE data in dairy cattle will likely suffer from the same problem. The frequencies of female and male calves were more equal in later parities.
Genetic parameters 1 for calving ease, stillbirth and gestation length from first parity bivariate animal models
Genetic parameters 1 for calving ease, stillbirth and gestation length from later parities bivariate animal models
Genetic parameters 1 for calving ease, stillbirth and gestation length between parities and within traits
Variances and genetic parameters for calving ease, stillbirth and gestation length from first parity univariate models
Variance components and genetic parameters3
σ P 2
σ Ad 2
σ Am 2
σ TBV 2
σ e 2
h d 2
h m 2
Accuracies of prediction (r) of average first parity PTA from 25 young 1 and older 2 sires
Average progeny group size
Bivariate with GL
Bivariate with SB
Bivariate between parities
Bivariate with CE
Bivariate with SB
Bivariate between parities
Bivariate with CE
Bivariate with GL
Bivariate between parities
Table 2 and 3 show the estimates of heritabilities and genetic correlations among traits obtained from bivariate animal models in first and later parities and Table 4 presents parameter estimates obtained by between-parity models for each trait respectively. Therefore, results for later parities in Table 4 account for selection based on first parity, whereas results in Table 3 do not. All heritabilities estimated for CE, SB and GL were within the range of previously published estimates of these traits in dairy cattle [16, 23, 33]. Heritabilities of CE were low (direct: 12% first parity and 3% later parities; maternal: 5% first parity and 2% later parities) and the direct heritability was approximately twice as large as the maternal heritability. GL appeared moderately heritable, with the direct heritability (57% first parity and 41% later parities) being considerably larger than the maternal heritability (7% first parity and 7% later parities). This supports the view that the genetic variation of this trait lies primarily in the triggering of parturition by the foetus [22, 33, 34] rather than in the maternal response to this trigger. All heritabilities were larger in first parity than in later parities, as reported elsewhere in the literature [7, 23]. This supports the general assumption that the variation in calving performance is larger in heifers than in cows [24, 35]. In addition, heritability estimates are frequency dependent when applying linear models to categorical traits.
Both direct and maternal heritabilities for SB were low, irrespectively of parity, with the direct heritability ranging from 1.8% to 2.0% (not significant in first parity) and the maternal heritability ranging from 2.0% to 3.2%. These estimates agree with values from the literature, which range from 1.6% to 10% for direct heritability and from 2.0% to 13% for maternal heritability [22, 32–34, 36, 37].
Direct and maternal genetic correlations
Within traits and within parities
The estimated genetic direct-maternal correlations for CE and GL presented in Table 23 and 4 were low to moderate (−0.52 to −0.22). For GL, the direct-maternal correlations were not significantly different from zero. For CE, a significant genetic relationship between the direct and maternal effects in first parity (−0.53) was detected. This negative direct-maternal correlation suggests that animals with a lower genetic risk of being born with difficulties are genetically prone to have more difficulty at first calving. Numerous studies confirm a negative genetic relationship between the direct and maternal effect of CE [38–40], although positive correlations also appear in the literature . The negative genetic correlation between direct and maternal effects of CE implies that dairy farmers need to base selection decisions on both the direct and the maternal PTA of a sire for CE in first parity, to avoid long-term negative consequences. An optimum index value for genetic merit in CE is therefore preferable, as discussed later.
Due to very low heritabilities and very high standard errors, the estimated direct-maternal genetic correlation of SB, obtained by the different models, were not informative. Studies estimating this covariance in large datasets (> 400 000) report correlations close to zero [20, 23, 38], although with considerable standard errors. To date, there is no clear evidence to recommend a change from the common practice of assuming this covariance as equal to zero.
Table 2 shows the estimated genetic correlations between the direct and maternal effects of CE, SB and GL in first parity heifers obtained from bivariate analyses. Table 3 shows the estimates of the same models for later parity cows (parity 2 and 3). In general, CE and SB were strongly genetically correlated, whereas the relationships of GL with CE and SB were weak to moderate. Both the direct and maternal correlations between CE and SB were positive and high in first parity (0.84; 0.85), and positive and moderate in later parities (0.37; 0.67). This suggests that both difficult birth and difficult calving are genetically associated with a higher frequency of direct and maternal stillbirth respectively, regardless of parity. The findings for the UK dataset thereby support the consensus of a strong genetic relationship between CE and SB [20, 37–41]. In this study, correlations of CE and SB were not significantly different from 1. However, Hansen et al.  and Cervantes et al.  provide evidence of genetic distinctiveness for these traits, with estimates of similar magnitude to those from this study, but with smaller standard errors. Furthermore, a simple meta-analysis  pooling estimates from this study and four other studies [20, 21, 41, 42] results in a direct and maternal genetic correlations of between CE and SB of 0.79 ± 0.02 and 0.65 ± 0.03, respectively, which suggests a genetic distinctiveness of these traits.
The genetic correlations between GL and the calving traits differed between parities (Table 2 and 3). A moderate positive genetic correlation (0.65) was found between maternal GL and maternal SB in first parity. This suggests that an individual with a longer than average gestation period is genetically more likely to give birth to a stillborn calf in first parity and vice versa. Genetic correlations between GL and SB in later parities were not significant (Table 3). A direct genetic relationship between GL and CE was detected, but only in later parities. No maternal relationship was detected in later parities. The direct effect of GL was found to be moderately correlated to the direct effect of CE. This positive correlation (0.50) between direct CE and direct GL suggests that a calf that gestates longer before birth to a multiparous dam is genetically prone to a difficult birth and vice versa. Similar positive correlations between direct GL and direct CE are reported in beef cattle , Danish Holstein cattle  and UK Holstein cattle , and support the findings from the UK dataset here. However, in this study, the genetic correlation between maternal GL and maternal CE was non-significant, although this relationship is generally reported to be low to moderate [20, 21].
All relationships between direct effects of one trait and maternal effects of the other trait (and vice versa) were non-significant, except for the genetic correlation between direct SB and maternal CE in first parity, which was high and positive. This specific relationship is difficult to estimate at the animal level and the high estimate may be due to the inaccuracy of the SB variance components. In general, studies in the literature report non-significant genetic correlations between the genetic direct effects and the genetic maternal effects between traits [20, 25, 41, 44].
Between parities within traits
Table 4 presents the genetic parameters estimated by the bivariate S-MGS models that treat first and later parity records as correlated traits. Estimated genetic correlations between first and later parities were 0.80 ± 0.12 for direct CE and 0.84 ± 0.15 for maternal CE. These estimates are similar to those estimates obtained by the threshold model reported by Wiggans et al.,  but slightly higher than those reported in general [23, 37, 41]. Among the studies estimating genetic correlations of CE between parities, there is general agreement that both direct and maternal CE are genetically distinct traits in first and later parities, which suggests that both ease of birth and ease of calving represent a different trait in heifers and in cows [23, 24, 41]. However, the standard errors reported here are too large to infer genetic distinction between first and later parities from this study alone.
Direct and maternal GL are rarely considered separately in studies that estimate between-parity correlations. Table 4 shows that different between-parity genetic correlations are found for direct GL and maternal GL. This emphasizes the fact that direct and maternal GL are separate traits, and thus must be analysed and interpreted with this in mind. For direct GL, the estimated correlation between first and later parities was near unity (0.96 ± 0.02) but the same correlation for maternal GL was lower (0.82 ± 0.13). However, in this case too, the standard error is too large to conclude that maternal GL is a distinct trait in first and later parities. Other studies on larger datasets show a correlation that is high but nevertheless significantly different from 1 [45, 46] which implies that maternally, the gestation length of a heifer and a cow are genetically distinct traits.
When fitting a between-parity S-MGS model for SB, results were difficult to obtain. Other analyses, using univariate and bivariate models, already showed the difficulty of obtaining an accurate estimate of the direct-maternal genetic correlation for SB within parities. With the between-parity model, the likelihood surface was practically flat which hampered convergence to sensible estimates.
Total heritable variance
The total heritable variance gives a holistic measure of the genetic variance affecting calving and accounts for both the maternal and direct sources of variance. Although some estimates in this study were inaccurate (in particular the direct-maternal genetic covariance for SB), the estimates of σ TBV 2 presented in Table 5 show how the maternal variance and direct-maternal genetic covariance contribute to the total genetic variance. Focusing on animal model A2, the total variance was smaller than the direct variance for CE and GL by 26% and 3%, respectively, although these differences were not significant. For SB, the very large and positive direct-maternal genetic covariance, in combination with the small direct variance, caused σ TBV 2 to be by ~ 400% larger than the additive direct variance.
When a farmer makes a selection decision based on a maternally affected trait, population mean performances change in response to both its direct and maternal breeding value. The TBV i as described by Bijma et al.  is suggested to represent the total additive value of an individual. However, it does not represent the impact of that individual on the population mean since this impact will depend on the time period and the frequency of expression of the direct and maternal effects in the population within that period. Gene flow methodology [47, 48] shows that contributions of the direct and maternal effects to genetic change in calving traits depend on several factors which determine how often the maternal effect is expressed, e.g. how many calvings, how many calves are kept as replacement heifers and the breeding system (pure breeding or crossbreeding). Therefore, while theoretically TBV i and σ TBV 2 show the importance of considering maternal effects and their interrelationship with direct effects, practically, an index value that is not the simple sum and represents the total impact of an individual would be useful to farmers. This would be in addition to the separate direct and maternal EBV that are already provided.
In this section, animal models are compared to S-MGS models, and univariate models to bivariate models. Furthermore, the benefit in treating first and later-parity calvings as correlated traits in a bivariate between-parity model is discussed. A potential bias due to an environmental direct-maternal covariance is also evaluated. All models that are discussed are linear models. Several studies have explored the advantages of threshold models over linear models for the analysis of calving traits [49, 50] given that according to the categorical nature of the traits, threshold models should theoretically be superior, as explained by Gianola . Findings show that computational requirements are greater for threshold compared to linear models and Monte Carlo methods needed to obtain the most reliable parameter estimates. However, software that can estimate variances without relying on Monte Carlo simulation methods, e.g. through the use of approximations to maximum likelihood in complex GLMM is available but limited . For calving traits, comparisons between linear and threshold models have shown very high correlations between PTA, meaning that the ranking of sires is not greatly influenced by the use of a linear model [35, 49, 53]. Threshold models have been shown to take specific interactions into account which can potentially be problematic for linear models [54, 55].
Animal model versus S-MGS model
Table 5 contains the results of the univariate analyses on first parity data and compares animal model A1 with the S-MGS model, since neither of these two models are expected to show bias due to the ignored environmental direct-maternal covariances discussed in the Methods, unlike animal model A2.
Table 5 shows that, between traits, direct heritabilities and phenotypic variances when estimated by A1 and S-MGS models were very similar but the residual variances were consistently larger for S-MGS models compared to both A1 and A2 models. The residual variance of an S-MGS model contains the default environmental variance plus a Mendelian sampling term and the remaining unexplained additive variance terms from dams totalling . Accuracies of PTA for young and older sires are presented in Table 6. Comparison of the S-MGS model with A1 and A2 models shows that there was only a small loss in accuracy when fitting the S-MGS model. In some cases, mainly for older sires for which more progeny information is available, there was an increase in accuracy of PTA when the S-MGS model was fitted as opposed to the animal model. This is probably due to the slightly higher heritabilities that were estimated by the S-MGS model (see Table 5). The computation time required with univariate animal models was 10 times greater than with univariate S-MGS models. Furthermore, when increasing the model complexity, animal models failed to converge, whereas S-MGS models were robust. The between-parity model in this study was an example where animal models failed, whereas S-MGS models performed well. Advantages in computation time and versatility of the S-MGS model therefore compensate well for the slight loss in accuracy of any resulting estimates.
Potential bias in the direct-maternal genetic correlation
Table 5 shows the estimated genetic direct-maternal correlations within traits for first parity data, when applying two univariate animal models. As described earlier, part of the data corresponding to specific dam-offspring pairs was deleted from the dataset (7.8%) to remove for a potential environmental direct-maternal covariance. In Table 5, animal model A1 represents the analysis of the edited data, whereas A2 represents the analysis of the complete data. Comparison of the results for A1 and A2 models shows that deleting records on dam-offspring pairs had only a small and non-significant impact on direct-maternal genetic correlations. The observed change implies that the environmental direct-maternal covariance was negative and small in this dataset. Since estimates of the direct-maternal genetic correlation do not differ significantly, the magnitude of the environmental covariance in this dataset is likely to be negligible and changes observed could be due to chance alone.
Univariate versus bivariate models
Table 6 demonstrates that, for calving traits, bivariate models are slightly superior to univariate traits. In particular, the maternal variance of low heritable traits (CE and SB) benefitted from including a correlated trait in the model. Accuracies of direct PTA also showed a slight improvement from bivariate analysis, in particular for the low heritable SB trait. Estimates obtained for SB heritabilities with the univariate and bivariate models (Table 5 and 2) demonstrate that these models provided nearly identical estimated direct variances, although the maternal variance showed a small but significant higher estimate with bivariate analyses. Strong genetic correlations were found between CE and SB, SB and GL, and CE and GL, which are likely to explain the increase in accuracy of PTA obtained by the bivariate model compared to the univariate model. The maternal variance of the indicator trait, GL, also benefitted from the bivariate model compared to the univariate model, although the accuracy of the maternal PTA of GL was slightly decreased with the bivariate model. Genetic direct-maternal correlations for CE and GL showed little change between univariate and bivariate models, while the estimate of this correlation for SB showed considerably more change but is too imprecise to be interpreted.
Inclusion of later parities
Calving traits in first and later parities were highly correlated, which results in a considerably greater predictive ability of PTA for all traits when later parity information is added as a correlated trait to the model (Table 6). Accuracies increased, for both direct and maternal PTA of CE and GL, when compared to the univariate model.
Heritabilities for CE, SB and GL in UK Holstein cattle were in the range of previously reported genetic parameters for these traits. Both the direct and maternal genetic variances were considerably lower in cows than in heifers. Direct and maternal effects of CE were negatively correlated but this was established only in first parity. CE and SB were genetically highly correlated traits for both direct and maternal components, especially in first parity. GL showed a moderate relationship with CE and SB, which differed between parities but implies that genetically longer gestations are associated with reduced calving performance. The three traits all had high and positive genetic correlations between parities but parities were not demonstrated as being genetically distinct for any trait with the data available. Different between-parity genetic correlations estimated for direct GL and maternal GL emphasize that these are separate traits and thus should be treated as such. Estimates of σ TBV 2 indicate that the total additive genetic variance in a calving trait may be lower than the additive direct variance when the genetic direct-maternal covariance is highly negative and the additive maternal variance is small.
Results from this study further demonstrated that estimating genetic parameters for calving traits is complex. Developing a statistical model for a maternally affected trait requires a careful balance between sufficient predictive ability and computational practicality, which in turn are affected by the size of the dataset, potential biases in data recording, the trait in question, computational facilities and the amount of time in hand. However, in general, PTA estimates for calving traits benefitted from multi-trait models. Furthermore, estimates were only slightly less accurate when a S-MGS model was fitted instead of an animal model. With the current computing facilities, S-MGS models exceeded animal models in terms of practicality, as their robustness allowed the analysis of more data and the inclusion of more traits e.g. information from later parities. In the genetic evaluation of calving traits genetic correlations between traits and between parities need to be estimated and the direct-maternal genetic correlation must be considered with caution.
The authors gratefully acknowledge funding from Defra under the Sustainable Livestock Production LINK Programme, the Scottish Government, CIS, Cogent, DairyCo, Genus, Holstein UK and NMR.
- Dekkers JC: Optimal breeding strategies for calving ease. J Dairy Sci. 1994, 77: 3441-3453. 10.3168/jds.S0022-0302(94)77287-8.View ArticlePubMedGoogle Scholar
- Dematawewa CM, Berger PJ: Effect of dystocia on yield, fertility, and cow losses and an economic evaluation of dystocia scores for Holsteins. J Dairy Sci. 1995, 80: 754-761.View ArticleGoogle Scholar
- Eaglen SA, Coffey MP, Woolliams JA, Mrode R, Wall E: Phenotypic effects of calving ease on the subsequent fertility and milk production of dam and calf in UK Holstein-Friesian heifers. J Dairy Sci. 2011, 11: 5413-5423.View ArticleGoogle Scholar
- Mee JF: Prevalence and risk factors for dystocia in dairy cattle: A review. Vet J. 2007, 176: 93-101.View ArticleGoogle Scholar
- Mark T: Applied genetic evaluations for production and functional traits in dairy cattle. J Dairy Sci. 2004, 87: 2641-2652. 10.3168/jds.S0022-0302(04)73390-1.View ArticlePubMedGoogle Scholar
- De Maturana E, Legarra A, Varona L, Ugarte E: Analysis of fertility and dystocia using recursive models to handle censored and categorical data. J Dairy Sci. 2007, 90: 2012-2024. 10.3168/jds.2005-442.View ArticleGoogle Scholar
- Meijering A: Dystocia and stillbirth in cattle – A review of causes, relations and implications. Livest Prod Sci. 1984, 11: 143-177. 10.1016/0301-6226(84)90057-5.View ArticleGoogle Scholar
- Willham RL: The covariance between relatives for characters composed of components contributed by related individuals. Biometrics. 1963, 19: 18-27. 10.2307/2527570.View ArticleGoogle Scholar
- Willham RL: The role of maternal effects in animal breeding: III. Biometrical aspects of maternal effects in animal. J Anim Sci. 1972, 35: 1288-1293.PubMedGoogle Scholar
- Falconer DS: Maternal effects and selection response. Genetics Today: Proceedings of the XIth International Congress of Genetics. 1965, Pergamon, Oxford, 763-774.Google Scholar
- Information on national evaluations providing data to the Interbull routine evaluation for calving traits, as provided by participating countries. . -[http://www-interbull.slu.se/national_ges_info2/framesida-ges.htm]
- Willham RL: Problems in estimating maternal effects. Livest Prod Sci. 1980, 7: 405-418. 10.1016/0301-6226(80)90080-9.View ArticleGoogle Scholar
- Meyer K: Variance components due to direct and maternal effects for growth traitsof Australian beef cattle. Livest Prod Sci. 1992, 31: 179-204. 10.1016/0301-6226(92)90017-X.View ArticleGoogle Scholar
- Koch RM: The role of maternal effects in animal breeding: VI. Maternal effects in beef cattle. J Anim Sci. 1972, 35: 1316-1323.PubMedGoogle Scholar
- Robinson DL: Models which might explain negative correlation between direct and maternal genetic effects. Livest Prod Sci. 1996, 45: 111-122. 10.1016/0301-6226(96)00002-4.View ArticleGoogle Scholar
- Eaglen SA, Bijma P: Genetic parameters of direct and maternal effects for calving ease in Dutch Holstein-Friesian cattle. J Dairy Sci. 2009, 92: 2229-2237. 10.3168/jds.2008-1654.View ArticlePubMedGoogle Scholar
- Koerhuis AN, Thompson R: Models to estimate maternal effects for juvenile bodyweight in broiler chickens. Genet Sel Evol. 1997, 29: 225-249. 10.1186/1297-9686-29-2-225.PubMed CentralView ArticleGoogle Scholar
- Lee C, Pollak EJ: Relationship between sire x year interactions and direct-maternal genetic correlation for weaning weight of Simmental cattle. J Anim Sci. 1997, 75: 68-75.PubMedGoogle Scholar
- Berweger Baschnagel M, Moll J, Künzi N: Comparison of models to estimate maternal effects for weaning weight of Swiss Angus cattle fitting a sire x herd interaction as an additional random effect. Livest Prod Sci. 1999, 60: 203-208. 10.1016/S0301-6226(99)00093-7.View ArticleGoogle Scholar
- Hansen M, Lund MS, Pedersen J, Christensen LG: Gestation length in Danish Holsteins has weak genetic associations with stillbirth, calving difficulty, and calf size. Livest Prod Sci. 2004, 91: 23-33. 10.1016/j.livprodsci.2004.06.007.View ArticleGoogle Scholar
- De Maturana EL, Gianola D, Rosa GJ, Wiegel KA: Predictive ability of models for calving difficulty in US Holsteins. J Anim Breed Genet. 2009, 126: 179-188. 10.1111/j.1439-0388.2008.00792.x.View ArticlePubMedGoogle Scholar
- Matilainen K, Mrode R, Strandén I, Thompson R, Mäntysaari EA: Linear-threshold animal model for birth weight, gestation length and calving ease in United Kingdom Limousin beef cattle data. Livest Prod Sci. 2009, 122: 143-148. 10.1016/j.livsci.2008.08.006.View ArticleGoogle Scholar
- Steinbock L, Näsholm A, Berglund B, Johansson K, Philipsson J: Genetic effects on stillbirth and calving difficulty in Swedish Holsteins at first and second calving. J Dairy Sci. 2003, 86: 2228-2235. 10.3168/jds.S0022-0302(03)73813-2.View ArticlePubMedGoogle Scholar
- Wiggans GR, Van Tassell CP, Cole JB, Thornton LL: Genetic correlations between first and later parity calving ease in a S-MGS model. Proceedings of the 8th World Congress on Genetics Applied to Livestock Production: 13–18 August 2008; Belo Horizonte. 2008, , , communication 01–92Google Scholar
- De Maturana EL, Wu XL, Gianola D, Wiegel KA, Rosa GJM: Exploring biological relationships between calving traits in primiparous cattle with a bayesian recursive model. Genetics. 2008, 181: 227-287.View ArticleGoogle Scholar
- Bijma P, Muir WM, Van Arendonk JAM: Multilevel selection 1: Quantitative genetics of inheritance and response to selection. Genetics. 2007, 175: 277-288.PubMed CentralView ArticlePubMedGoogle Scholar
- SAS Institute: SAS/STAT Software Release 9.1. 2006, SAS Institute Inc, CaryGoogle Scholar
- Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release 3.0. 2009, VSN International Ltd, Hemel HempsteadGoogle Scholar
- Ducrocq V: Calving ease evaluation of French dairy bulls with a heteroskedastic threshold model with direct and maternal effects. Interbull Bull. 2000, 25: 123-130.Google Scholar
- Wiggans GR, Misztal I, Van Tassell CP: Calving ease (co)variance components for aS-MGS threshold model. J Dairy Sci. 2003, 86: 1845-1848. 10.3168/jds.S0022-0302(03)73771-0.View ArticlePubMedGoogle Scholar
- Corbin LJ, Blott SC, Swinburne JE, Vaudin M, Bishop SC, Woolliams JA: Linkage disequilibrium and historical effective population size in the Thoroughbred horse. Anim Genet. 2010, 41: 8-15.View ArticlePubMedGoogle Scholar
- Jamrozik J, Fatehi J, Kistemaker GJ, Schaeffer LR: Estimates of genetic parameters for Canadian Holstein female reproduction traits. J Dairy Sci. 2005, 88: 2199-2208. 10.3168/jds.S0022-0302(05)72895-2.View ArticlePubMedGoogle Scholar
- Norman HD, Wright JR, Kuhn MT, Hubbard SM, Cole JB, VanRaden PM: Genetic and environmental factors that affect gestation length in dairy cattle. J Dairy Sci. 2009, 92: 2259-2269. 10.3168/jds.2007-0982.View ArticlePubMedGoogle Scholar
- Azzam SM, Nielsen MK: Genetic parameters for gestation length, birth date andfirst breeding date in beef cattle. J Anim Sci. 1987, 64: 348-356.Google Scholar
- Clutter AC, Berger PJ, Mattison JM: Threshold model analysis of dystocia in dairy cattle when progeny information is limited. J Dairy Sci. 1989, 72: 3264-3272. 10.3168/jds.S0022-0302(89)79487-X.View ArticlePubMedGoogle Scholar
- Luo MF, Boettcher PJ, Schaeffer LR, Dekkers JC: Estimation of genetic parameters of calving ease in first and second parities of Canadian Holsteins using Bayesian methods. Livest Prod Sci. 2002, 74: 175-184. 10.1016/S0301-6226(01)00294-9.View ArticleGoogle Scholar
- Heringstad B, Chang YM, Svendsen M, Gianola D: Genetic analysis of calving difficulty and stillbirth in Norwegian red cows. J Dairy Sci. 2007, 90: 3500-3507. 10.3168/jds.2006-792.View ArticlePubMedGoogle Scholar
- Thompson JR, Freeman AE, Berger PJ: Age of dam and maternal effects for dystocia in Holsteins. J Dairy Sci. 1980, 64: 1603-1609.View ArticleGoogle Scholar
- Carnier PA, Dal Zotto R, Groen AF, Bona M, Bittante G: Genetic parameters for direct and maternal calving ability over parities in Piedmontese cattle. J Anim Sci. 2000, 78: 2532-2539.PubMedGoogle Scholar
- Lee D, Misztal J, Bertrand JK, Rekaya R: National evaluation for calving ease, gestation length and birth weight by linear and threshold model methodologies. J Appl Genet. 2002, 43: 209-216.PubMedGoogle Scholar
- Cervantes I, Gutiérrez JP, Fernández I, Goyache F: Genetic relationships among calving ease, gestation length, and calf survival to weaning in the Asturiana de los Valley beef cattle breed. J Anim Sci. 2010, 88: 96-101. 10.2527/jas.2009-2066.View ArticlePubMedGoogle Scholar
- Thompson JR, Pollak EJ, Pelissier CL: Interrelationships of parturition problems, production of subsequent lactation, reproduction, and age at first lactation. J Dairy Sci. 1983, 66: 1119-1127. 10.3168/jds.S0022-0302(83)81909-2.View ArticlePubMedGoogle Scholar
- McGuirk BJ, Going I, Gilmour AR: The genetic evaluation of beef sires used for crossing with dairy cows in the UK: 1. Sire breed and non-genetic effects on calving survey traits. Anim Sci. 1998, 66: 35-45. 10.1017/S135772980000881X.View ArticleGoogle Scholar
- Johanson JM, Berger PJ, Tsuruta S, Misztal I: A Bayesian threshold-linear model evaluation of perinatal mortality, dystocia, birth weight, and gestation length in a Holstein herd. J Dairy Sci. 2011, 94: 450-460. 10.3168/jds.2009-2992.View ArticlePubMedGoogle Scholar
- Van Pelt ML, de Jong G, Eding H, Roelfzema JE: Analysis of calving traits with a multi-trait animal model with a correlated direct and maternal effect. Interbull Bull. 2009, 40: 138-141.Google Scholar
- Norman HD, Wright JR, Miller RH: Potential consequences of selection to change gestation length on performance of Holstein cows. J Dairy Sci. 2010, 94: 1005-1010.View ArticleGoogle Scholar
- Wolfová M, Nitter G: Relative economic weights of maternal versus direct traits in breeding schemes. Livest Prod Sci. 2004, 88: 117-127. 10.1016/j.livprodsci.2003.09.024.View ArticleGoogle Scholar
- Wall E, Coffey M, Sadeghi-Sefidmazgi A, Amer P: A model for deriving economic weights for calving ease in UK dairy cattle. Proceeding of the British Society of Animal Science Annual Conference: 4–5 April 2011; Nottingham. 2011, 015Google Scholar
- Phocas F, Laloë D: Evaluation models and genetic parameters for calving difficulty in beef cattle. J Anim Sci. 2003, 81: 933-938.PubMedGoogle Scholar
- Lee C: On the negative estimates of direct and maternal genetic correlation – A review. Asian-Aust J Anim Sci. 2002, 15: 1222-1226.View ArticleGoogle Scholar
- Gianola D: Theory and analysis of threshold characters. J Anim Sci. 1982, 54: 1079-1096.Google Scholar
- McCulloch CE: Maximum likelihood algorithms for generalized linear mixed models. J Am Stat Assoc. 1997, 92: 162-170. 10.1080/01621459.1997.10473613.View ArticleGoogle Scholar
- Weller JI, Gianola D: Models for genetic analysis of dystocia and calf mortality. J Dairy Sci. 1989, 72: 2633-2643. 10.3168/jds.S0022-0302(89)79404-2.View ArticlePubMedGoogle Scholar
- Quaas RL, Zhao Y, Pollak EJ: Describing interactions in dystocia scores with a threshold model. J Anim Sci. 1988, 66: 396-399.PubMedGoogle Scholar
- Manfredi E, Ducrocq V, Foulley JL: Genetic analysis of dystocia in dairy cattle. J Dairy Sci. 1991, 74: 1715-1723. 10.3168/jds.S0022-0302(91)78334-3.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.