Phenotypic effects of calving ease on the subsequent fertility and milk production of dam and calf in UK Holstein-Friesian heifers

How does a beef × dairy calving affect the dairy cow's following lactation?

For beef semen usage on dairy cows, much of the research has focused on the performance of the crossbred calves, yet little focus has been given to the subsequent performance of the cow herself. This study aimed to evaluate the performance of dairy cows for milk yield, fertility, and survival traits after giving birth to beef × dairy crossbred calves and compare this with the performance after giving birth to purebred dairy calves. Further, we aimed to study if the effect of a difficult calving was the same regardless of whether the calf was purebred dairy or beef × dairy crossbred. Phenotypic records from 587,288 calving events from 1997 to 2020 were collected from the Swedish milk recording system from cows of the dairy breeds Swedish Red (SR) and Swedish Holstein. The sire beef breeds studied were Aberdeen Angus, Hereford (combined in category LHT), Charolais, Limousin, and Simmental (category HVY). Sixteen traits were defined and grouped in 3 categories: cumulative and 305-d milk, fat, and protein yield, daily milk yield, and 75-d milk yield as yield traits; calving to first insemination interval, calving to last insemination interval, first to last insemination interval, calving interval, and number of inseminations as fertility traits; and survival to 75 d or to next calving and lactation length as measures of survival. The data were analyzed for all traits for first and second parities separately using mixed linear models, with a focus on the estimates of cow breed by service sire breed combinations. All traits in parity 2 were adjusted for previous 305-d milk yield based on the expectation that low-yielding cows would more likely to be inseminated with beef semen. Overall, milk yield was lower after beef × dairy calvings compared with the purebred dairy calvings. The largest effects were found on cumulative yields and in second parity, with lower effects for yields early in lactation and yields in first parity. The largest decrease was 13 to 14 kg (0.12 phenotypic SD) for cumulative fat yield when breeding beef breed sires with purebred SR dams. For fertility traits, for most breed combinations, the effects were not large enough to be significant. Conversely, all beef × dairy crossbred combinations showed significantly lower results for survival to the next lactation, and mostly also for lactation length. There was some indication that dairy cows with beef × dairy calvings in parity 2 that were the result of maximum 2 inseminations in parity 1, had lower survival than corresponding calvings resulting from more than 2 inseminations. This could indicate that the former cows were marked for culling already when inseminated. There was generally an unfavorable effect of a difficult calving on all traits, however, there were almost no significant interactions between calving performance and dam by sire breed combination, and these interactions were never significant in first parity.

Genetic evaluation for stillbirth and pre-weaning mortality in Australian dairy cattle

The welfare of calves is important to both farmers and consumers. Practices that increase the proportion of calves born alive and enable them to thrive through to weaning contribute to improved sustainability. Stillbirths (SB) are calvings where the calf dies at birth or within 24 h after birth. Pre-weaning mortality (PWM) refers to calves that die after the first day of life but before weaning based on termination data. Both SB and PWM are binary traits characterized by low heritability. Data collection for these traits is incomplete, compared with traits like milk yield in cows. Despite these challenges, genetic variation can be measured and used to produce breeding tools, such as estimated breeding values (EBV), to reduce calf mortality over time. The aim of this study was to compare the performance of various linear models to predict SB and PWM traits in Holstein and Jersey cattle and evaluate their applicability for industry-wide use in the Australian dairy industry. Calving records from around 2.25 million Holstein and Jersey dams were obtained from DataGene's Central Data Repository from 2000 onwards to calculate genetic parameters. About 7% of calves were recorded as stillborn in the period 2000-2021 (n = 1.48 million calvings). The prevalence of PWM was much lower than stillbirth during the same period at 2% (n = 0.89 million calves). Genetic parameters were estimated for SB direct, SB maternal and PWM using bivariate linear models with calving ease (CE) as the second trait in the model. The heritability of these calf traits was low and varied between 1 to 5% depending on the breed, trait and model. In Holstein cattle, heritabilities were 2% for PWM and SB direct and 1% for SB maternal while in Jersey cattle heritabilities were 5% for PWM, 2% for SB direct and 1% for SB maternal. The genetic trends for both SB direct and maternal in Holstein cattle indicate improvement in both traits whereas there was no apparent increase or decrease in PWM in the past 2 decades. The coefficient of genetic variation for SB direct and PWM was between 11.7 and 23.0% in Holstein and Jersey cattle demonstrating that there was considerable genetic variation in calf survival traits as a first step to using genetic selection to increase the proportion of calves born alive and calves weaned. A focus on improved calf and calving recording practices is expected to increase the reliability of genetic predictions.

Association between dam and calf measurements with overall and fetopelvic dystocia in Holstein heifers

We investigated the relationship between dam's pelvic and calf's dimensions with dystocia due to fetopelvic disproportion in the Holstein breed and estimated risk factors and dystocia probability. For this purpose, external pelvic measurements were performed in 402 heifers 15 ± 11 (1-38) days ante-partum and specific conformation measurements were obtained from their calves 1.7 ± 1.2 post-partum. Dystocia was defined as the inability of the heifer to complete parturition spontaneously within 120 min after the appearance of the amnion with normal presentation, position and posture or as having definite obstetrical obstacles within 60 min. Overall and fetopelvic disproportion dystocia incidence was 10.4% and 5.2%, respectively. Heifer measurements mainly influenced overall dystocia, whereas calf conformation was related solely with fetopelvic dystocia. Specifically, heifers with a small pelvis (hip width <49.95 cm, pelvic inlet area <333.2 cm2, pelvic volume <7799.2 cm3) had 2.8 to 3.5 times greater incidence of overall dystocia (19.0-20.8%) compared to heifers with a larger pelvis (incidence of 7.0-7.6%). Regarding calf factors, sex (male calves), body weight, chest circumference and fetlock joint circumference significantly increased the odds of experiencing dystocia due to fetopelvic disproportion compared with female, lighter or smaller calves. In a backward elimination model with independent variables treated as continuous, an area under the ROC curve of 0.66 regarding the prediction of overall dystocia based on heifer pelvic length, and of 0.64 for the prediction of fetopelvic dystocia based on fetlock joint circumference was found. The combination of the two variables in one model improved the ROC area to 0.71 regarding dystocia due to fetopelvic disproportion, reaching acceptable level of discrimination. Our findings indicate that dystocia due to fetopelvic disproportion in heifers is mainly influenced by the fetal side. Additionally, the estimation of pelvic dimensions of the dam before parturition and specific conformation characteristics of the calf during parturition, especially fetlock joint circumference, could aid obstetricians and herdsmen regarding dystocia probability and parturition surveillance.

Effects of the Breeding Strategy Beef-on-Dairy at Animal, Farm and Sector Levels

The decline in farm revenue due to volatile milk prices has led to an increase in the use of beef semen in dairy herds. While this strategy ("Beef-on-dairy" (BoD)) can have economic benefits, it can also lead to unintended consequences affecting animal welfare. Semen sale trends from breeding organizations depict increasing sales of beef semen across the globe. Calves born from such breeding strategies can perform better when compared to purebred dairy calves, especially in terms of meat quality and growth traits. The Beef-on-dairy strategy can lead to unintentional negative impacts including an increase in gestation length, and increased dystocia and stillbirth rates. Studies in this regard have found the highest gestation length for Limousin crossbred calves followed by calves from the Angus breed. This increase in gestation length can lead to economic losses ranging from 3 to 5 US$ per animal for each additional day. In terms of the growth performance of crossbred animals, literature studies are inconclusive due to the vast differences in farming structure across the regions. But almost all the studies agree regarding improvement in the meat quality in terms of color, fiber type, and intra-muscular fat content for crossbred animals. Utilization of genomic selection, and development of specialized Beef-on-dairy indexes for the sires, can be a viable strategy to make selection easier for the farmers.

Birth traits of Holstein calves compared with Holstein, Jersey, Montbéliarde, Normande, and Viking Red-sired crossbred calves

Holstein (HO) calves, 3-breed crossbred calves of Montbéliarde, Viking Red, and HO (MVH), and 3-breed crossbred calves of Normande, Jersey, and Viking Red (NJV) were compared for gestation length (GL), calf weight at birth (CW), calving difficulty (CD), and stillbirth (SB) in 2 research herds at the University of Minnesota. Calves were born from January 2009 to December 2019. For the St. Paul and Morris herds, HO calves (n = 1,121) were compared with MVH calves (n = 1,393) from primiparous and multiparous cows. For the single herd analysis at Morris, HO calves (n = 476), MVH calves (n = 922), and NJV calves (n = 405) were compared from primiparous and multiparous cows. Primiparous and multiparous births were analyzed separately because multiparous cows had multiple births, and CD and SB are likely different traits for primiparous and multiparous cows. Statistical analysis of GL, CW, CD, and SB included fixed effects of sex of calf, herd, breed group of calf, and year-season of calving. For the St. Paul and Morris herds, HO calves from primiparous (278 d) and multiparous (279 d) HO cows had shorter GL compared with MVH calves from primiparous (280 d) and multiparous (282 d) crossbred cows. The HO calves (39.4 and 43.2 kg, respectively) from primiparous and multiparous HO cows had lower CW compared with MVH calves (40.3 and 44.3 kg, respectively) from primiparous and multiparous crossbred cows. Calving difficulty and SB were not different for HO and MVH calves from primiparous and multiparous cows. For the single herd analysis at Morris, HO calves (278 and 279 d, respectively) from primiparous and multiparous HO cows had shorter GL compared with MVH calves (281 and 282 d, respectively) and NJV calves (282 and 282 d, respectively) from primiparous and multiparous crossbred cows. The CW of HO calves (38.6 and 42.0 kg, respectively) from primiparous and multiparous HO cows was lower compared with MVH calves (39.7 and 42.9 kg, respectively), but higher compared with NJV calves (35.1 and 38.0 kg, respectively) from primiparous and multiparous crossbred cows. Calving difficulty and SB did not differ for HO, MVH, and NJV calves from primiparous and multiparous cows. The longer GL for crossbred calves and higher CW for MVH calves did not increase CD and SB for primiparous and multiparous cows. Dairy producers may implement 3-breed rotational crossbreeding systems that include the HO, Jersey, Normande, Montbéliarde, and Viking Red breeds, and some breeds may increase GL and CW without an increase in CD and SB.

Birth conditions affect the longevity of Holstein offspring

Studies of dairy cow longevity usually focus on the animal life after first calving, with few studies considering early life conditions and their effects on longevity. The objective was to evaluate the effect of birth conditions routinely collected by Dairy Herd Improvement agencies on offspring longevity measured as length of life and length of productive life. Lactanet provided 712,890 records on offspring born in 5,425 Quebec dairy herds between January 1999 and November 2015 for length of life, and 506,066 records on offspring born in 5,089 Quebec dairy herds between January 1999 and December 2013 for length of productive life. Offspring birth conditions used in this study were calving ease (unassisted, pull, surgery, or malpresentation), calf size (small, medium, or large), and twinning (yes or no). Observations were considered censored if the culling reason was "exported," "sold for dairy production," or "rented out" as well as if the animals were not yet culled at the time of data extraction. If offspring were not yet culled when the data were extracted, the last test-day date was considered the censoring date. Conditional inference survival trees were used in this study to analyze the effect of offspring birth conditions on offspring longevity. The hazard ratio of culling between the groups of offspring identified by the survival trees was estimated using a Cox proportional hazard model with herd-year-season as a frailty term. Five offspring groups were identified with different length of life based on their birth condition. Offspring with the highest length of life [median = 3.61 year; median absolute deviation (MAD) = 1.86] were those classified as large or medium birth size and were also the result of an unassisted calving. Small offspring as a result of a twin birth had the lowest length of life (median = 2.20 year; MAD = 1.69) and were 1.52 times more likely to be culled early in life. Six groups were identified with different length of productive life. Offspring that resulted from an unassisted or surgery calving and classified as large or medium when they were born were in the group with the highest length of productive life (median = 2.03 year; MAD = 1.63). Offspring resulting from a malpresentation or pull in a twin birth were in the group with the lowest length of productive life (median = 1.15 year; MAD = 1.11) and were 1.70 times more likely to be culled early in life. In conclusion, birth conditions of calving ease, calf size, and twinning greatly affected offspring longevity, and such information could be used for early selection of replacement candidates.

Effect of Calving Difficulties and Calf Mortality on Functional Longevity in Polish Holstein-Friesian Cows

Simple Summary

Longevity, or the length of a cow’s productive life, is important in terms of profitability, animal welfare and environmental sustainability. In genetic evaluations, interest focuses on functional longevity, defined as a cow’s ability to avoid forced culling, an ability that increases the possibility of voluntary disposal based on economic criteria. Longevity is affected by several non-productive functional traits, among them those related to calving performance: calving ease (dystocia) and perinatal calf mortality (stillbirth). Parturition is a critical event in a cow’s life that has a number of different short- and long-term consequences. In the Polish Holstein-Friesian population, the incidence of dystocia and stillbirth is within the lower range of frequencies found in other dairy cattle populations. Our research showed that both traits affect functional longevity. Difficult calvings occur more frequently in heifers and increase the risk of involuntary culling more than in later parturitions. Additionally, a higher risk of culling is related to birth of a male calf. Moreover, the negative impact of calf mortality on longevity is also more pronounced in primiparous cows and in the case of delivery of male calves. Reducing the incidence of calving problems and perinatal mortality may improve the longevity of dairy cows.

Abstract

Longevity is one of the functional traits that considerably affect dairy herd profitability. A Weibull proportional hazards model was used to evaluate the impact of difficult calvings and calf stillbirths on cow functional longevity, defined as length of productive life corrected for milk production. The data for analysis comprised calving ease and calf mortality scores of 2,163,426 calvings, 34.4% of which came from primiparous cows. The percentage of male calves was 53.4%. Calving ease was scored as “without assistance” (34.44%), “with assistance” (62.03%), “difficult—hard pull” (3.39%), and “very difficult, including caesarean section” (0.14%). Calf mortality scores were “live born” (94.21%) and “stillborn or died within 24 h” (5.79%). The Weibull proportional hazards model included classes of calving ease or calf mortality scores × parity (1, ≥2) × sex of calf as time-dependent fixed effect. The model also included time-dependent fixed effects of year × season, parity × stage of lactation, annual change in herd size, fat yield and protein yield, time-independent fixed effect of age at first calving, and time dependent random herd × year × season. In first-parity cows, very difficult birth of a bull or heifer increased the relative risk of culling, respectively, 2.18 or 1.26 times as compared with calving without assistance. In later parities, the relative risk of culling related to very difficult calving was 2.0 times (for male calves) and 1.33 times (for female calves) higher than the relative risk of culling associated with calving without assistance. Calf mortality showed a negative impact on longevity in both heifers and cows. First-parity stillbirth increased the relative risk of culling depending on sex of calf by 18% in females and by 15% in males; in later parities the increase of the relative risk of culling was lower (by 7% for females, 9% for males). Difficult calvings and their consequences, especially in primiparous cows, may negatively influence dairy herd profitability by reducing the length of cows’ productive life.

Influence of Calving Ease on In-Line Milk Lactose and Other Milk Components

The aim of our study was to determine how the ease of calving of cows may influence changes in lactose concentration and other milk components and whether these two factors correlate with each other. To achieve this, we compared data of calving ease scores and average percentage of in-line registered milk lactose and other milk components. A total of 4723 dairy cows from nine dairy farms were studied. The cows were from the second to the fourth lactation. All cows were classified according to the calving ease: group 1 (score 1)-no problems; group 2 (score 2)-slight problems; group 3 (score 3)-needed assistance; group 4 (score 4)-considerable force or extreme difficulty. Based on the data from the milking robots, during complete lactation we recorded milk indicators: milk yield MY (kg/day), milk fat (MF), milk protein (MP), lactose (ML), milk fat/lactose ratio (MF/ML), milk protein/lactose ratio (MP/ML), milk urea (MU), and milk electrical conductivity (EC) of all quarters of the udder. According to the results, we found that cows that had no calving difficulties, also had higher milk lactose concentration. ML > 4.7% was found in 58.8% of cows without calving problems. Cows with more severe calving problems had higher risk of mastitis (SCC and EC). Our data indicates that more productive cows have more calving problems compared to less productive ones.

Classification and regression tree analysis to predict calving ease in Sussex heifers using pelvic area dimensions and morphological traits

Objectives

The objective of this study was to assess the relationship between pelvic measurements, namely pelvic height (PH), pelvic area (PA), and pelvic width, and linear body measurement traits, such as hindquarters width, rump length (RL), hip height (HH), body weight (BW), shoulder height, chest depth (CD), body length, and also calve birth weight (CBW) under different calving ease scores, namely no assistance, gently pull, and hard pull.

Material and Methods

A total of 51 first calf Sussex heifers at 24 months old, weighing approximately 440 ± 40.26 kg, were used for this study. Two-year-old bulls (n = 2), weighing approximately 800 kg, were used for mating the 51 heifers with a bull to cow ratio of 1:30 and 1:21. The fertility of bulls was assessed by a private veterinarian before the mating season. Calving ease was scored as follows: 1 = no assistance during parturition (normal), 2 = heifer assistance as gently pull, 3 = heifer assistance as hard pull, 4 = heifer cannot calf, 5 = heifer calved a dead calf, and 6 = heifer calf with the abnormal position.

Results

Descriptive statistics indicated that PH (p < 0.01) and RL (p < 0.05) were positively remarkable and different among calving ease scores. Correlation results of the no assistance Sussex heifers group revealed that CBW was positively correlated (p < 0.05) with PH. Classification and regression tree CART results demonstrated that PH, BW, CD, and HH play an important role in predicting calving ease in Sussex heifers.

Conclusion

The findings suggest that there is an association among PAs and biometric traits. The present study might help farmers select animals with higher PH, BW, CD, and HH during breeding to lower dystocia incidents during parturition. However, further studies need to be carried out in CART modeling to reduce dystocia in a larger sample size of Sussex heifers or other cattle breeds.

Calving difficulty influences rumination time and inflammatory profile in Holstein dairy cows

Difficult calving may adversely affect dairy cow health and performance. Maternal:fetal disproportion is a major cause of dystocia. Therefore, the main objective of this study was to assess the effects of dam:calf body weight ratio (D:C) on calving difficulty, rumination time, lying time, and inflammatory profile in 25 Holstein dairy cows. Using automatic monitoring systems, we monitored behavior and production in 9 primiparous and 16 pluriparous cows between dry-off and 30 d in milk. During the same period, we collected blood samples to monitor metabolism and inflammatory profile of these cows. Calvings were video recorded to assess calving difficulty and observe the duration of the expulsive stage. After parturition, the cows were separated into 3 classes according to their D:C: easy (E; D:C >17), medium (M; 14 < D:C <17), and difficult (D; D:C <14). The cows in class D showed relatively longer labor durations (108 min vs. 54 and 51 min for classes D, M, and E, respectively) and higher calving assistance rates (50% vs. 0 and 11% of calvings for classes D, M, and E, respectively) than those in the other 2 classes. Compared with the cows in classes M and E, those in class D exhibited shorter rumination times on the day of calving (176 min/d vs. 288 and 354 min/d for classes D, M, and E, respectively) and during the first week of lactation (312 min/d vs. 339 and 434 min/d for classes D, M, and E, respectively) and maintained lower rumination values until 30 DIM (399 min/d vs. 451 and 499 min/d for classes D, M, and E, respectively). Primiparous class D cows had shorter resting times during the first week after calving compared with those in class M (8 vs. 11 h/d for classes D and M, respectively). Interclass differences were found in terms of the levels of inflammation markers such as acute-phase proteins (ceruloplasmin, albumin, retinol, and paraoxonase). Moreover, cows in class D had lower plasma levels of fructosamine and creatinine after calving. Low D:C reduced postcalving rumination time and increased inflammation grade, suggesting a lower welfare of these animals at the onset of lactation. The D:C might serve as a useful index for the identification of cows at relatively higher risk of metabolic and inflammatory disease, thus helping farmers and veterinarians improve the welfare and health of these cows.

Breed- and trait-specific associations define the genetic architecture of calving performance traits in cattle

Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10⁻⁸) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10⁻⁶) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10⁻⁸). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits.

A breeding index to rank beef bulls for use on dairy females to maximize profit

The desire to increase profit on dairy farms necessitates consideration of the revenue attainable from the sale of surplus calves for meat production. However, the generation of calves that are expected to excel in efficiency of growth and carcass merit must not be achieved to the detriment of the dairy female and her ability to calve and re-establish pregnancy early postcalving without any compromise in milk production. Given the relatively high heritability of many traits associated with calving performance and carcass merit, and the tendency for many of these traits to be moderately to strongly antagonistic, a breeding index that encompasses both calving performance and meat production could be a useful tool to fill the void in supporting decisions on bull selection. The objective of the present study was to derive a dairy-beef index (DBI) framework to rank beef bulls for use on dairy females with the aim of striking a balance between the efficiency of valuable meat growth in the calf and the subsequent performance of the dam. Traits considered for inclusion in this DBI were (1) direct calving difficulty; (2) direct gestation length; (3) calf mortality; (4) feed intake; (5) carcass merit reflected by carcass weight, conformation, and fat and the ability to achieve minimum standards for each; (6) docility; and (7) whether the calf was polled. Each trait was weighted by its respective economic weight, most of which were derived from the analyses of available phenotypic data, supplemented with some assumptions on costs and prices. The genetic merit for a range of performance metrics of 3,835 artificial insemination beef bulls from 14 breeds ranked on this proposed DBI was compared with an index comprising only direct calving difficulty and gestation length (the 2 generally most important characteristics of dairy farmers when selecting beef bulls). Within the Angus breed (i.e., the beef breed most commonly used on dairy females), the correlation between the DBI and the index of genetic merit for direct calving difficulty plus gestation length was 0.74; the mean of the within-breed correlations across all other breeds was 0.87. The ranking of breeds changed considerably when ranked based on the top 20 artificial insemination bulls excelling in the DBI versus excelling in the index of calving difficulty and gestation length. Dairy breeds ranked highest on the index of calving difficulty and gestation length, whereas the Holstein and Friesian breeds were intermediate on the DBI; the Jersey breed was one of the poorest breeds on DBI, superior only to the Charolais breed. The results clearly demonstrate that superior carcass and growth performance can be achieved with the appropriate selection of beef bulls for use on dairy females with only a very modest increase in collateral effect on cow performance (i.e., 2-3% greater dystocia expected and a 6-d-longer gestation length).

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

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

Including gene networks to predict calving difficulty in Holstein, Brown Swiss and Jersey cattle

Background

Calving difficulty or dystocia has a great economic impact in the US dairy industry. Reported risk factors associated with calving difficulty are feto-pelvic disproportion, gestation length and conformation. Different dairy cattle breeds have different incidence of calving difficulty, with Holstein having the highest dystocia rates and Jersey the lowest. Genomic selection becomes important especially for complex traits with low heritability, where the accuracy of conventional selection is lower. However, for complex traits where a large number of genes influence the phenotype, genome-wide association studies showed limitations. Biological networks could overcome some of these limitations and better capture the genetic architecture of complex traits. In this paper, we characterize Holstein, Brown Swiss and Jersey breed-specific dystocia networks and employ them in genomic predictions.

Results

Marker association analysis identified single nucleotide polymorphisms explaining the largest average proportion of genetic variance on BTA18 in Holstein, BTA25 in Brown Swiss, and BTA15 in Jersey. Gene networks derived from the genome-wide association included 1272 genes in Holstein, 1454 genes in Brown Swiss, and 1455 genes in Jersey. Furthermore, 256 genes in Holstein network, 275 genes in the Brown Swiss network, and 253 genes in the Jersey network were within previously reported dystocia quantitative trait loci. The across-breed network included 80 genes, with 9 genes being within previously reported dystocia quantitative trait loci. The gene-gene interactions in this network differed in the different breeds. Gene ontology enrichment analysis of genes in the networks showed Regulation of ARF GTPase was very significant (FDR ≤ 0.0098) on Holstein. Neuron morphogenesis and differentiation was the term most enriched (FDR ≤ 0.0539) on the across-breed network. Genomic prediction models enriched with network-derived relationship matrices did not outperform regular GBLUP models.

Conclusions

Regions identified in the genome were in the proximity of previously described quantitative trait loci that would most likely affect calving difficulty by altering the feto-pelvic proportion. Inclusion of identified networks did not increase prediction accuracy. The approach used in this paper could be extended to any instance with asymmetric distribution of phenotypes, for example, resistance to disease data.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0606-y) contains supplementary material, which is available to authorized users.

Bayesian estimates of genetic relationship between calving difficulty and productive and reproductive performance in Holstein cows

The objective of present study was to estimate genetic correlations between calving difficulty and productive and reproductive traits in Iranian Holsteins. Calving records from the Animal Breeding Center of Iran, collected from 1991 to 2011 and comprising 183 203 first-calving events of Holstein cows from 1470 herds were included in the dataset. Threshold animal models included direct genetic effect (Model 1) or direct and maternal genetic effects with covariance between them (Model 2) were fitted for the genetic analysis of calving difficulty. Also, linear animal models including direct genetic effect were fitted for the genetic analysis of productive and reproductive performance traits. A set of linear-threshold bivariate models was used for obtaining genetic correlation between calving difficulty and other traits. All analyses were implemented by Bayesian approach via Gibbs sampling methodology. A single Gibbs sampling chain with 300 000 rounds was generated by the TM program. Posterior mean estimates of direct heritabilities for calving difficulty were 0.056 and 0.066, obtained from different models. Also, posterior mean estimate of maternal heritability for this trait was 0.018. Estimate of correlation between direct and maternal genetic effects for calving difficulty was negative (–0.44). Posterior mean estimates of direct heritabilities for milk yield, fat yield, protein yield, days from calving to first service, days open and first calving interval were 0.257, 0.188, 0.235, 0.034, 0.042 and 0.050 respectively. The posterior means of direct genetic correlation between calving difficulty and milk yield, fat yield, protein yield, days from calving to first service, days open and first calving interval were low and equal to –0.135, 0.030, –0.067, –0.010, –0.075 and –0.074 respectively. The results of the current study indicated that exploitable genetic variation in calving difficulty, productive and reproductive traits could be applied in designing future genetic selection plans for Iranian Holsteins.

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle

Background

Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value predictions and followed by gene identification and in silico functional analyses. The association results were further validated using five scenarios with different numbers of SNPs.

Results

Seven hundred and eighty-two SNPs were significantly associated with calving performance at a genome-wise false discovery rate (FDR) of 5%. Most of these significant SNPs were on chromosomes 18 (71.9%), 17 (7.4%), 5 (6.8%) and 7 (2.4%) and mapped to 675 genes, among which 142 included at least one significant SNP and 532 were nearby one (100 kbp). For body conformation traits, 607 SNPs were significant at a genome-wise FDR of 5% and most of them were located on chromosomes 5 (30%), 18 (27%), 20 (13%), 6 (6%), 7 (5%), 14 (5%) and 13 (3%). SNP enrichment functional analyses for calving traits at a FDR of 1% suggested potential biological processes including musculoskeletal movement, meiotic cell cycle, oocyte maturation and skeletal muscle contraction. Furthermore, pathway analyses suggested potential pathways associated with calving performance traits including tight junction, oxytocin signaling, and MAPK signaling (P < 0.10). The prediction ability of the 1206 significant SNPs was between 78 and 83% of the prediction ability of the BovineSNP50 SNPs for calving performance traits and between 35 and 79% for body conformation traits.

Conclusions

Various SNPs that are significantly associated with calving performance are located within or nearby genes with potential roles in tight junction, oxytocin signaling, and MAPK signaling. Combining the significant SNPs or SNPs within or nearby gene(s) from the HD panel with the BovineSNP50 panel yielded a marginal increase in the accuracy of prediction of genomic estimated breeding values for all traits compared to the use of the BovineSNP50 panel alone.

Electronic supplementary material

The online version of this article (10.1186/s12711-017-0356-8) contains supplementary material, which is available to authorized users.

Estimation of genetic parameters for lambing ease, birthweight and gestation length in Australian sheep

This study presents estimates of genetic parameters for lambing ease (LE), birthweight (BW) and gestation length (GL) in Australian terminal sire sheep breeds using data from the Sheep Genetics LAMBPLAN database. LE was scored on lambs on a 1–5 scale, with 1 being no assistance and 5 being other such as special veterinary assistance. The full dataset consisted of 43 448 records on LE and its two subsets (the single and twin subsets) based on the birth type of the litter were analysed. Four models with different combinations of random effects consisting of direct genetic, maternal genetic and maternal permanent environmental effects were compared. All traits were analysed using linear animal models and linear sire models with LE further analysed by threshold sire models for all datasets to evaluate the influence of datasets and models on the estimation of genetic parameters. The results showed that multiple-born lambs had shorter GL, less BW and less lambing difficulty than single-born lambs. Lambing difficulty decreased with the increase of dam age from 1.5 to 4.5 years, and then increased afterwards. Genetic parameters using linear animal models were similar to those using linear sire models for all traits. Phenotypic variance and direct heritability were higher for single-born lambs compared with twin-born lambs. No significant maternal permanent environmental effect was detected for LE. Based on results using linear animal models with the full dataset, the direct heritabilities were 0.06 ± 0.01, 0.15 ± 0.01 and 0.52 ± 0.02, the maternal heritabilities were 0.03 ± 0.01, 0.15 ± 0.01 and 0.13 ± 0.02 for LE, BW and GL, respectively. The proportions of maternal permanent environmental effects to the total variances were 0.13 ± 0.01 for BW. Low to moderate direct genetic correlations of 0.31 ± 0.09 (LE and BW), 0.24 ± 0.11 (LE and GL) and 0.08 ± 0.08 (BW and GL) were estimated using tri-variate analysis from the full dataset, indicating the trend that lambs with greater BW and longer GL would result in more lambing difficulty.

Danish holsteins favor bull offspring: biased milk production as a function of fetal sex, and calving difficulty

In a previous study from 2014 it was found that US Holstein cows that gave birth to heifer calves produced more milk than cows having bull calves. We wanted to assess whether this is also true for Danish cattle. Data from 578 Danish Holstein herds were analysed with a mixed effect model and contrary to the findings in the US, we found that cows produced higher volumes of milk if they had a bull calf compared to a heifer calf. We found a significantly higher milk production of 0.28% in the first lactation period for cows giving birth to a bull calf, compared to a heifer calf. This difference was even higher when cows gave birth to another bull calf, so having two bull calves resulted in a difference of 0.52% in milk production compared to any other combination of sex of the offspring. Furthermore, we found that farmer assisted calvings were associated with a higher milk yield. Cows with no farmer assistance or with veterinary assistance during the most recent calving produced less milk. There were also indications that dams would favor a bull fetus by decreasing milk production during the second pregnancy if the calf born in the first parity was a heifer. We hypothesize that size of calves is a confounding factor for milk production. However, calving weight was not available in the present data set to test this hypothesis.

Effect of dams' parity and age on daughters' milk yield in Norwegian Red cows

The effect of age and parity of dams on their daughters' milk yield is not well known. Lactation data from 276,000 cows were extracted from the Norwegian Dairy Herd Recording System and analyzed using a linear animal model to estimate effects of parity and age within parity of dam. The 305-d milk yield of daughters decreased as parity of dam increased. Daughters of first-parity dams produced 149 kg more milk than did daughters of seventh-parity dams. We also observed an effect of age of dam within parity on 305-d milk yield of daughters in first lactation. Dams that were young at first calving gave birth to daughters with a higher milk yield compared with older dams within the same parity. The effect of age within parity of dam was highest for second-parity dams. Extensive use of heifers would have a systematic effect, and age and parity of dam should be included in the model when planning a future strategy.

Genetic parameters for direct and maternal calving ease in Walloon dairy cattle based on linear and threshold models

Calving ease scores from Holstein dairy cattle in the Walloon Region of Belgium were analysed using univariate linear and threshold animal models. Variance components and derived genetic parameters were estimated from a data set including 33,155 calving records. Included in the models were season, herd and sex of calf × age of dam classes × group of calvings interaction as fixed effects, herd × year of calving, maternal permanent environment and animal direct and maternal additive genetic as random effects. Models were fitted with the genetic correlation between direct and maternal additive genetic effects either estimated or constrained to zero. Direct heritability for calving ease was approximately 8% with linear models and approximately 12% with threshold models. Maternal heritabilities were approximately 2 and 4%, respectively. Genetic correlation between direct and maternal additive effects was found to be not significantly different from zero. Models were compared in terms of goodness of fit and predictive ability. Criteria of comparison such as mean squared error, correlation between observed and predicted calving ease scores as well as between estimated breeding values were estimated from 85,118 calving records. The results provided few differences between linear and threshold models even though correlations between estimated breeding values from subsets of data for sires with progeny from linear model were 17 and 23% greater for direct and maternal genetic effects, respectively, than from threshold model. For the purpose of genetic evaluation for calving ease in Walloon Holstein dairy cattle, the linear animal model without covariance between direct and maternal additive effects was found to be the best choice.

Direct and maternal genetic relationships between calving ease, gestation length, milk production, fertility, type, and lifespan of Holstein-Friesian primiparous cows

As the emphasis in cattle breeding is shifting from traits that increase income toward traits that reduce costs, national breeding indices are expanding to include functional traits such as calving ease (CE). However, one issue is the lack of knowledge of genetic relationships between CE and other dairy traits. The same can be said about gestation length (GL), a potential novel selection trait with considerable heritabilities and possible genetic relationships with the calving process. This study aimed to estimate the genetic relationships between CE, GL, and other dairy traits of interest using a national data set of 31,053 primiparous cow performance records, as well as to separate direct and maternal genetic effects. Chosen dairy traits included fertility (calving interval, days to first service, nonreturn rate after 56 d, number of inseminations per conception), milk production (milk yield at d 110 in milk, accumulated 305-d milk yield, accumulated 305-d fat yield, accumulated 305-d protein yield), type (udder depth, chest width, rump width, rump angle, mammary composition, stature, body depth), and lifespan traits (functional days of productive life). To allow the separation of direct and maternal genetic effects, a random sire of the calf effect was included in the multi-trait linear trivariate sire models fitted using ASReml. Significant results showed that easily born individuals were genetically prone to high milk yield and reduced fertility in first lactation. Difficult calving primiparous cows were likely associated with being high-producing, wide and deep animals, with a reduced ability to subsequently conceive. Individuals that were born relatively early were associated with good genetic merit for milk production. Finally, individuals carrying their offspring longer were genetically associated with being wide and large animals that were themselves born relatively early. The study shows that it is feasible and valuable to separate direct and maternal effects when estimating genetic correlations between calving and other dairy traits. Furthermore, gestation length is best used as an indicator trait for lowly heritable calving traits, rather than as a novel selection trait. As estimated direct and maternal genetic correlations differ, we can conclude that genetic relationships between CE, GL, and traits of interest are present, but caution is required if these traits are implemented in national breeding indices.

Board invited review: The importance of the gestation period for welfare of calves: maternal stressors and difficult births

The prenatal period is of critical importance in defining how individuals respond to their environment throughout life. Stress experienced by pregnant females has been shown to have detrimental effects on offspring biology in humans and a variety of other species. It also is becoming increasingly apparent that prenatal events can have important consequences for the behavior, health, and productivity of offspring in farmed species. Pregnant cattle may experience many potentially important stressors, for instance, relating to their social environment, housing system and physical environment, interactions with humans and husbandry procedures, and their state of health. We examined the available literature to provide a review of the implications of prenatal stress for offspring welfare in cattle. The long-term effects of dystocia on cattle offspring also are reviewed. To ensure a transparent and repeatable selection process, a systematic review approach was adopted. The research literature clearly demonstrates that prenatal stress and difficult births in beef and dairy cattle both have implications for offspring welfare and performance. Common husbandry practices, such as transport, were shown to influence offspring biology and the importance of environmental variables, including thermal stress and drought, also were highlighted. Maternal disease during pregnancy was shown to negatively impact offspring welfare. Moreover, dystocia-affected calves suffer increased mortality and morbidity, decreased transfer of passive immunity, and important physiological and behavioral changes. This review also identified considerable gaps in our knowledge and understanding of the effects of prenatal stress in cattle.

Short communication: survival, growth to weaning, and subsequent fertility of live-born dairy heifers after a difficult birth

The experience of a difficult birth (dystocia) is traumatic and has adverse effects on the newborn in various species. Despite affecting up to 1 in 3 births in dairy cattle, studies on calves have been mostly limited to the first day of life. The objective of this study was to investigate the effects of dystocia on the survival to calving, growth to weaning, and subsequent fertility as nulliparous animals. Historical data from live-born Holstein heifer calves born from cows with various birth difficulty scores (no assistance; moderate; high difficulty) were obtained from 2 herds (Edinburgh herd: n=1,237; Crichton Royal Farm herd: n=721). Each herd was analyzed separately for birth weights, weaning weights, growth rate to weaning, number of services to conception, and age at first calving using REML and generalized linear mixed model analyses. Survival analysis (Cox proportional hazards model) was used in the Edinburgh herd to analyze the subsequent survival of live-born heifers, whereas descriptive data are presented for the Crichton Royal Farm herd. A higher mortality risk to weaning and to first service was observed in the live-born heifers that experienced moderate difficulty at birth compared with heifers born naturally. Surviving dystocial heifers had similar growth-to-weaning and fertility performance as heifers born naturally in both herds. It could be that the performance of dystocial heifers that survived to weaning was not affected or that it was compensated for by farm management. This study highlights long-term effects of the early experience of a difficult birth and thereby stresses the importance of preventing dystocia not just from the point of view of the adult cow, but also from the perspective of the calf. This would also improve farm efficiency and calf welfare.

Effect of inbreeding depression on survival of Austrian Brown Swiss calves and heifers

The aim of this study was to estimate inbreeding depression for juvenile mortality in Austrian Brown Swiss replacement heifers born in the years 2001 to 2007. After data editing, records of 69,571 animals were investigated. In total, the pedigree consisted of 203,894 animals. Mean and median inbreeding coefficients were 0.0514 and 0.0475, respectively. The following periods were defined for analyses of juvenile mortality: P1=48 h to 30 d, P2=31 to 180 d, P3=181 to the day before first calving or a maximum age of 1,200 d if no calving was reported, P4=48 h to the day before age at first calving or a maximum age of 1,200 d if no calving was reported. Mortality during the first 30 d of life was 3.2%; in the total rearing period, 9.3% of the heifers died (excluding slaughtered and exported animals). For the estimation of the effect of inbreeding and of genetic parameters a linear animal model with the random effects herd-year of birth and animal as well as the fixed effects year of birth-month of birth and dam parity, and the continuous effect of inbreeding coefficient (linear and quadratic) was applied. The model was reduced to the linear effect of inbreeding as the quadratic term was not found to have a significant effect. Inbreeding significantly affected all traits unfavorably. In the full rearing period (P4), heifer calves with an inbreeding coefficient of 0.10 had a 4.9% higher mortality rate than heifer calves with no inbreeding.

Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle

Background

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.

Methods

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.

Conclusions

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.