Analysis of the relationship between type traits and functional survival in US Holstein cattle using a Weibull proportional hazards model

Survival analysis with a Weibull proportional hazards model was used to evaluate the effects of 15 linear type traits, 5 composite traits, and final score on the functional longevity of US Holstein cows. Culling data and type classification scores (measured in first lactation) from 891,524 cows with first calving from 1993 to 2000 were used. The data were divided into 9 geographical regions to determine whether the relationship between type traits and longevity differed according to climate or management system. Functional survival was defined as the number of days from first calving until culling or censoring, after correction for 305-d mature equivalent combined fat and protein yield. The Weibull model included time-dependent effects of herd-year-season, parity-stage of lactation, and within herd-year quintile ranking for combined fat and protein yield (nested within biennium), as well as time-independent effects of age at first calving and type classification score (type traits were analyzed one at a time). Type classification scores were rounded to the nearest 5 points, and the impact of each type trait on functional survival in each region was evaluated. Mean failure time ranged from 694 d in the South to 758 d in the North East. Risk of culling differed by region for several linear type traits, and differences were greatest for regions that were most dissimilar in climate and herd management (e.g., South East, East North Central, and West). Udder depth, fore udder attachment, udder cleft, and rear legs side view were consistently associated with functional longevity, regardless of region, but, the importance of some secondary traits, such as stature or dairy form, differed by region. The survival model applied in this study easily described both linear and nonlinear relationships between type traits and longevity while accounting for important time-dependent and time-independent explanatory variables.     

Relationships Among Type Scores and Changes in Yield from First to Second Lactation

Changes between first and second lactation yields of days open adjusted milk and 3.7% FCM were correlated with classification traits. Yields within first and second lactation were corrected for herd-season, calving month within season, and age at calving. Classification scores during first lactation were available on 1341 cows appraised between 1971 and 1982 by Holstein Association of America classifiers in five North Carolina Department of Agriculture herds and a North Carolina State University research herd. All classification traits except final score were categorical. Final score, general appearance, dairy character, body capacity, mammary system, stature, rear udder, fore udder, and feet and legs were treated as continuous. Classification traits were adjusted for herd-year, age, and months in milk. First lactation milk was added as a covariate in selected analyses. Body capacity, dairy character, mammary system, back, rump, rear udder, and udder support each had significant effects on change in yield. Mammary system and final score showed a significant linear relationship with change in yield for at least one model. Dairy character and fore udder significantly predicted change in 3.7% FCM when first lactation milk was not in the model.     

Heritability of clinical mastitis incidence and relationships with sire transmitting abilities for somatic cell score, udder type traits, productive life, and protein yield

The objective of this study was to determine the relationships among daughter clinical mastitis during first and second lactations and sire transmitting abilities for somatic cell score, udder type traits, productive life, and protein yield. Data on clinical mastitis during first lactation were available for 1795 daughters (in six Pennsylvania herds, one Minnesota herd, and one Nebraska herd) of 283 Holstein sires. Data on clinical mastitis during second lactation were available for 1055 of these daughters. A total of 479 cows had 864 clinical episodes during first lactation, and 230 cows had 384 clinical episodes during second lactation. Clinical mastitis incidence and the total number of clinical episodes during each lactation were regressed on herd-season of calving (a classification variable), age at first calving, lactation length, and sire transmitting abilities taken one at a time. Linear effects, nonlinear effects, and odds ratios were estimated for sire transmitting abilities. Separate analyses were conducted on dependent variables that considered clinical mastitis from: all organisms, coagulase-negative staphylococci, coliform species, streptococci other than Streptococcus agalactiae, and the most common environmental organisms (coliform species and streptococci other than Streptococcus agalactiae). Heritability of clinical mastitis ranged from 0.01 to 0.42. Daughters of sires that transmit the lowest somatic cell score had the lowest incidence of clinical mastitis and the fewest clinical episodes during first and second lactations. Daughters of sires that transmit longer productive life, shallower udders, deeper udder cleft, and strongly attached fore udders had either fewer clinical episodes or lower clinical mastitis incidence during first and second lactations. The incidence of clinical mastitis and the number of clinical episodes per lactation may be reduced by selection for lower somatic cell score, longer productive life, shallower udders, deeper udder cleft, or strongly attached fore udders.     

Inheritance of cow durability for linear type traits

Final score and 14 linear type traits from appraisal by the Holstein Association were analyzed to assess inheritance of cow durability. Change of type scores from first to second parity (second minus first) was cow durability. Score for first parity, change of score from first to second parity, and FCM were simultaneously evaluated with a multiple-trait mixed model for each type trait. First analysis did not allow missing scores of cows for second parity, but subsequent analysis permitted missing data for second parity. When cows were required to have type scores for second parity, heritabilities for change of score were from zero (rear legs and foot angle) to 6% (final score). Genetic correlations of FCM and change of score were from -.35 (rump angle) to .24 (stature). When cows were not required to have type scores from second lactation, heritabilities for change of score were less than 8% for all traits except udder depth (17%). Except for dairy form and fore udder attachment, all genetic correlations of type scores from first lactation and change of score were positive. Seemingly, little new knowledge of linear type traits for ranking of sires was provided by appraisal of cows after first parity.     

Genetic correlations between lifetime production and linearized type in Canadian Holsteins.

Genetic and phenotypic correlations were estimated between 6 lifetime production and 28 linearized type traits using REML. The data set contained 34,322 cows, each with a record for all 34 traits. The analyses accounted for the fixed effects of herd, year-month, classifier, age at first calving, and stage of lactation. Heritabilities were low for lifetime traits and moderate for most type traits except stature, size, capacity, thurl width, and pin setting, which had high heritabilities. Most phenotypic correlations between lifetime production and type were in the range of .15 to .20 except for capacity, rump, and feet and legs, which were around .07. Genetic correlations were strong between lifetime production and angularity (.44 to .55) and dairy character (.53 to .56). Genetic correlations were low to moderate between life-time production and stature (.14 to .25), size (.07 to .18), texture (.19 to .26), style (.11 to .27), head (.15 to .23), pin setting (.10 to .16), rear udder (.19 to .25), and rear attachment (.10 to .22). The only notable negative genetic correlations were lifetime production with rear heel (-.16 to -.27), thurl width (-.18 to -.24), and fore udder (-.05 to -.11).     

Short communication: Modification of genetic evaluation of herd life from a three-trait to a five-trait model in Canadian dairy cattle

The national genetic evaluation of herd life for Canadian dairy breeds was modified from a 3-trait to a 5-trait animal model. The genetic evaluation incorporates information from daughter survival (direct herd life) and information from conformation, fertility, and udder health traits that are related to longevity (indirect herd life). Genetic evaluations for direct herd life were based on cows’ survival from first calving to 120 days in milk (DIM), from 120 to 240 DIM, from 240 DIM to second calving, survival to third calving, and survival to fourth calving, which were analyzed using a multiple-trait animal model. Sire evaluations obtained for each of the 5 survival traits were combined into an overall sire evaluation for direct herd life. Sire evaluations for indirect herd life were based on an index of sire evaluations for dairy strength, feet and legs, overall mammary, rump angle, somatic cell score, milking speed, nonreturn rate in cows, and interval from calving to first service. A multiple-trait sire model based on multiple-trait across-country evaluation methodology was used to combine direct and indirect genetic evaluations for herd life into an overall genetic evaluation for herd life. Sire evaluations for herd life were expressed as an estimated transmitting ability for the number of lactations. The transmitting ability represents expected differences among daughters for herd life; and the average herd life was set to 3 lactations.     

Correlations among linear type traits and somatic cell counts

Genetic and phenotypic correlations between linear type traits and SCC were estimated from lactation average SCC from Pennsylvania DHIA and Holstein linear type evaluations from Sire Power, Inc. and Holstein Association using REML. Correlations were estimated between linear type traits and SCC measured in first lactation and between linear type traits measured in first lactation and SCC measured in second or third lactation. Data sets ranged in size from 4294 daughters of 216 sires to 58,235 daughters of 301 sires. Phenotypic and genetic correlations between the linear traits that reflect body and locomotive characteristics and SCC were generally small and unimportant. Phenotypic correlations between udder traits and SCC were variable, but cows with higher udder depth scores (higher udders) had lower SCC. Genetic correlations between udder traits and SCC were also variable. Genetic correlations between SCC and udder depth, SCC and fore udder attachment, and SCC and teat placement were negative (favorable). Genetic correlations between teat length and SCC tended to be positive. Genetic correlations were largest in magnitude between udder depth and SCC and ranged from -.21 to -.64 (weighted mean = -.35). Selection for higher udders and closer teat placement will likely improve resistance to mastitis in dairy cattle.     

Relationship among udder and teat morphology and milking characteristics.

Teat cup liner slips, manual milking machine adjustments, milk yields, and milking times were recorded during both morning and evening milkings for 8 d on 97 Holstein cows in The Pennsylvania State University dairy herd. Fore and rear udder heights (distance from floor to udder), udder levelness, distances between teats (before and after milking), teat lengths, teat diameters, and teat end shapes were measured on the same cows. Product-moment correlations among the morphological characteristics, linear slips, manual adjustments, milk yields, and milking times were determined. Residual correlations from a model including lactation number and DIM (linear and quadratic) were also calculated. The variation among cows in machine liner slips and manual adjustments within and across lactation number and DIM can be partially explained by udder and teat morphology. Wider teats were associated with increased linear slips and increased manual adjustments. More tilted udders (rear quarters lower than front quarters) were associated with increased liner slips and tended to be associated with increased manual adjustments. In addition, larger teat diameters and longer teats tended to be associated with increased liner slips.     

Genetic parameters of conformation traits, milk yield, and herd life in Holsteins.

Genetic parameters were estimated simultaneously for 5 herd-life traits, 15 conformation (type) traits, and milk yield measured in first lactation for 128,601 Holstein cows. Heritabilities of all traits were higher in registered than in grade cows. Genetic correlations of linear type traits with first lactation yield ranged from -.48 for udder depth to .54 for dairy form. Genetic correlations among milk yield and herd-life traits were all positive except for milk-corrected herd life in grade cows. Udder traits had largest absolute genetic correlations with herd-life traits, followed by body traits and feet and leg traits. Some traits associated with body size and foot angle differed between registered and grade cows. Estimates of genetic trends from obtained parameters revealed greatest progress for milk yield from single-trait selection but also the largest changes for some type traits and milk-corrected herd life in an undesirable direction. Relative milk to type ratios of between 2:1 and 3:1 yielded 90% of the gain in milk yield with no change or slight improvement in type traits and functional herd life. Selection for type traits associated with herd life appears to be warranted to improve days of functional herd life or to decrease involuntary culling of dairy cows.     

Analysis of the relationship between type traits, inbreeding, and functional survival in Jersey cattle using a Weibull proportional hazards model

A Weibull proportional hazards model was used to analyze the effects of 13 linear type traits, final score, and inbreeding on the functional survival of 268,008 US Jersey cows in 2416 herds with first calving from 1981 to 2000. Functional survival was defined as the number of days from first calving until involuntary culling or censoring. The statistical model included the time-dependent effects of herd-year-season of calving, parity by stage of lactation interaction, and within-herd-year quintile for mature equivalent milk yield, as well as the time-independent effects of inbreeding, age at first calving, and linear type traits or final score (analyzed one at a time). Each type trait was divided into 10 classes, and the relative risk of involuntary culling was calculated for animals in each class after accounting for the aforementioned management factors. Type traits with the greatest contribution to the likelihood function were udder depth, fore udder attachment, front teat placement, and udder support. Cows with low scores for these traits had a risk of culling that was 1.3 to 1.8 times that of cows with intermediate scores. Cows with high scores for udder depth and udder support had a risk of culling only 0.7 to 0.85 as great as that of cows with intermediate scores. Intermediate scores were desirable for rear leg set, dairy form, and strength, but stature, rump angle, and rump width had negligible effects on survival. Cows with low final scores had a risk of culling that was 1.35 times that of cows with intermediate scores, whereas cows with high final scores had a risk of culling that was 0.8 times that of cows with intermediate scores. Animals with inbreeding coefficients greater than 10% had a slightly higher risk of culling than animals with inbreeding coefficients less than 5%.     

Analysis of the relationship between workability traits and functional longevity in Canadian dairy breeds

The aim of this study was to assess the effect of workability traits like milking speed and temperament on functional longevity of Canadian dairy cattle using a Weibull proportional hazards model. First-lactation data consisted of the following: 1,728,289 and 2,426,123 Holstein cows for milking temperament and milking speed, respectively, from 18,401 herds and sired by 8,248 sires; 39,618 and 60,121 Jersey cows for milking temperament and milking speed, respectively, from 1,845 herds and sired by 2,413 sires; and 54,391 and 94,847 Ayrshire cows for milking temperament and milking speed, respectively, from 1,316 herds and sired by 2,779 sires. Functional longevity was defined as the number of days from the first calving to culling, death, or censoring adjusted for production. Milking temperament and milking speed were recorded on a 1- to 5-point scale from very nervous to very calm and from very slow to very fast, respectively. The statistical model included the effects of stage of lactation; season of production; the annual change in herd size; type of milk recording supervision; age at first calving; effects of milk, fat, and protein yields calculated as within herd-year-parity deviations; herd-year-season of calving; sire; and milking temperament or milking speed class. The relative culling rate was calculated for animals in each milking temperament or milking speed class after accounting for the above-mentioned effects. The study showed that there was a statistically significant association between workability traits and functional longevity. Very nervous cows were 26, 23, and 46% more likely to be culled than very calm cows in Holstein, Ayrshire, and Jersey breeds, respectively. Similarly, very slow milkers were 36, 33, and 28% more likely to be culled than average milkers in Holstein, Ayrshire, and Jersey breeds, respectively. Additionally, very fast milkers were 11, 13, and 15% more likely to be culled than average milkers in Holstein, Ayrshire, and Jersey breeds, respectively. Producers might want to avoid consequences associated with the fast milkers such as udder health problems.     

Relationship between reproduction traits and functional longevity in canadian dairy cattle

The aim of this study was to use survival analysis to assess the relationship between reproduction traits and functional longevity of Canadian dairy cattle. Data consisted of 1,702,857; 67,470; and 33,190 Holstein, Ayrshire, and Jersey cows, respectively. Functional longevity was defined as the number of days from first calving to culling, death, or censoring; adjusted for the effect of milk yield. The reproduction traits included calving traits (calving ease, calf size, and calf survival) and female fertility traits (number of services, days from calving to first service, days from first service to conception, and days open). The statistical model was a Weibull proportional hazards model and included the fixed effects of stage of lactation, season of production, the annual change in herd size, and type of milk recording supervision, age at first calving, effects of milk, fat, and protein yields calculated as within herd-year-parity deviations for each reproduction trait. Herd-year-season of calving and sire were included as random effects. Analysis was performed separately for each reproductive trait. Significant associations between reproduction traits and longevity were observed in all breeds. Increased risk of culling was observed for cows that required hard pull, calved small calves, or dead calves. Moreover, cows that require more services per conception, a longer interval between first service to conception, an interval between calving to first service greater than 90 d, and increased days open were at greater risk of being culled.     

Relationships between herd type and genetic and environmental variances in Holsteins

Type classification data from the Holstein Association of Canada on 354,308 daughters of 13,694 sires were subdivided into 20 groups of approximately equal size according to herd average final score. Genetic and environmental variances were estimated for each group by Henderson's new method. Traits considered and their heritabilities averaged over groups were: final score .13; final class .11; dairy character .17; capacity .24; rump .16; feet and legs .10; mammary system .12; fore udder .12; rear udder .12; size .31. Regressions of genetic variances on group number ordered according to herd final score indicated significant trends only for rump and feet and legs for which higher genetic variances were found for higher type scores. Regressions of environmental variances on group number were significant for all traits. Higher environmental variances were found for lower final scores, final class, dairy character, and capacity, but the opposite was observed for size. The distribution of environmental variances across groups was U-shaped for rump, feet and legs, for udder, and rear udder; lowest variances were associated with intermediate scores for overall type.     

Phenotypic and genetic relationship between linear functional type traits and milk yield for five breeds

Lactation records of cows first calving between 18 and 35 mo were combined with linear type ratings assigned during the same lactation if before 43 mo. Phenotypic relationships were examined between final score and 13 type appraisal traits and first lactation milk yield from 2935 Ayrshire, 3154 Brown Swiss, 13,110 Guernsey, 50,422 Jersey, and 924 Milking Shorthorn records. Most phenotypic correlations between type and milk yield were low. Linear correlations of final score with first lactation milk yield were .18 to .38. Of the linear functional type traits, correlations with first lactation yield had greatest absolute values for dairy character (.19 to .53), udder depth (-.26 to -.30), and rear udder width (.20 to .31). Multiple correlations of all type appraisal traits with first lactation milk yield ranged from .41 to .59. Herd-year-season components of variance averaged 25% for type traits. Herd-year-season with sire interaction averaged 4%. Heritability estimates for final score from paternal half-sib analysis were from .11 to .21. Heritability estimates for linear traits ranged from .01 to .37. Genetic correlation between milk and final score was positive for Guernseys (.25) and Jerseys (.21). Genetic correlations between yield and most linear type traits were low to moderate except for dairy character (.53 to .77).     

Genetic analysis of herd life in Quebec Holsteins using Weibull models

Survival analysis methodologies were used to study herd life in Canadian Holstein cows. Herd life was defined as true herd life or the length of time between first calving and censoring. True herd life adjusted for 305-d milk production was defined as functional herd life. Lifetime record (censored or completed) were from 331,147 Holstein cows registered in the Programme d'Analyse des Troupeaux Laitiers du Québec (PATLQ) that calved for the first time between March 1, 1981 and March 31, 1995. The Weibull (proportional hazards) model used to analyze true herd life and functional herd life contained a Weibull baseline hazard function and the time-dependent effects of year of first calving, lactation number by stage of lactation, annual change in herd size and herd-year (random), and the time-independent effects of the milk recording option (supervised or not) and age at first calving. The model for functional herd life included also the time-dependent effect of herd-year-parity class of 305-d milk production. Genetic differences between sires with regard to the hazard function of their daughters was clearly demonstrated. The hazard rate followed a different pattern in later lactations, particularly in the first 240 d in milk. Older age at first calving was found to be associated with higher risks of culling. Changes in herd size had a small impact on the hazard function of animals. The hazard decreased as production of the cow increased. Heritability in the log scale was 0.09 for true herd life and 0.08 for functional herd life, but when heritability was expressed on the original scale, the estimates for the two traits were 0.19 and 0.15, respectively. The difference in the median survival between a bull with an estimated transmitting ability of 0.6 and another bull with an estimated transmitting ability of 1.3 was 690 d or 1.7 lactations. Rank correlations between the official estimated transmitting abilities for true herd life and functional herd life and those obtained in this study were 0.62 and 0.66, respectively.     

The genetic relationship between calving interval, body condition score and linear type and management traits in registered Holsteins

The trend to poorer fertility in dairy cattle with rising genetic merit for production over the last decade suggests that breeding goals need to be broadened to include fertility. This requires reliable estimates of genetic (co)variances for fertility and other traits of economic importance. In the United Kingdom at present, reliable information on calving dates and hence calving intervals are available for most dairy cows. Data in this study consisted of 44,672 records from first lactation heifers on condition score, linear type score, and management traits in addition to 19,042 calving interval records. Animal model REML was used to estimate (co)variance components. Genetic correlations of body condition score (BCS) and angularity with calving interval were -0.40 and 0.47, respectively, thus cows that are thinner and more angular have longer calving intervals. Genetic correlations between calving interval and milk, fat, and protein yields were between 0.56 and 0.61. Records of phenotypic calving interval were regressed on sire breeding values for BCS estimated from records taken at different months of lactation and breeding values for BCS change. Genetic correlations inferred from these regressions showed that BCS recorded 1 mo after calving had the largest genetic correlation with calving interval in first lactation cows. It may be possible to combine information on calving interval, BCS, and angularity into an index to predict genetic merit for fertility.     

Intercorrelations among milk production traits and body and udder measurements in Holstein heifers

Data from 1341 Holstein heifers of 71 sires were used to study heritabilities of and genetic and phenotypic correlations among milk production traits (308-d milk, front and rear half yields), body measurements (heart girth, withers height, body length, and rump length), udder measurements (front teat length and diameter, rear teat length and diameter, teat distance and udder height), and age at first calving. Genetic and phenotypic parameters were estimated by the multitrait restricted maximum likelihood method. Multitrait estimates of heritability ranged from .37 to .47 for first lactation yield traits, from .19 to .51 for body measurements, and from .08 to .41 for udder measurements. Age at first calving averaged 22.3 mo with a heritability estimate of .11. Milk production traits were all positively correlated with body measurements, suggesting that high producing heifers would be taller, larger, and longer than low producing heifers. Multitrait estimates of genetic and phenotypic correlations between udder height and yield traits were all negative, suggesting that high producing heifers tend to have lower udders. Of four body measurements studied, rump length showed the greatest genetic correlations with yield traits. Among six udder measurements, udder height exhibited the highest degree of associations with yield traits. Thus, rump length and udder height merit greater attention for prediction of lactational performance.