Association of length of pregnancy with other reproductive traits in dairy cattle.

The experiment involved observations of 2,514 Holstein-Friesian cows to determine the effects of environmental factors (cow's age, calving season, weight and sex of calves, housing system) and genetic factors on gestation length in dairy cattle and the correlation between gestation length and other reproductive traits (calving ease, stillbirth rates and placental expulsion). Genetic parameters were estimated based on the sires of calved cows (indirect effect) and the sires of live-born calves (direct effect). The following factors were found to contribute to prolonged gestation: increasing cow's age, male fetuses and growing fetus weight. Optimal gestation length was determined in the range of 275-277 days based on calving ease and stillbirth rates. The heritability of gestation length was estimated at 0.201-0.210 by the direct effect and 0.055-0.073 by the indirect effect. The resulting genetic correlations suggest that the efforts to optimize (prolong) gestation length could exert an adverse influence on the breeding value of bulls by increasing perinatal mortality and calving difficulty. The standard errors of the investigated parameters were relatively high, suggesting that any attempts to modify gestation length for the purpose of improving calving ease and reducing stillbirth rates should be introduced with great caution.

INTRODUCTION

Gestation length (GL), the period from effective fertilization until calving, is a reproductive trait that significantly affects cattle breeding and production. The cow’s age is the key environmental factor influencing gestation length. Gestation length is shorter in heifers than in older cows (Przysucha and Grodzki, 2009). High temperatures in the summer speed up calving and shorten gestation length (McClintock et al., 2003). High milk yield prolongs gestation, as suggested by a positive genetic correlation between GL and milk production levels (Silva et al. 1992). Both longer and shorter gestation periods contribute to a higher number of stillbirths (Norman et al., 2009). GL values should be analyzed in conjunction with calving ease (CE), stillbirth rates (SB) and placental expulsion (PE) (Hansen et al., 2004; Przysucha and Grodzki, 2009). Estimated heritability of GL based on the direct effect were determined in the range of h2 = 0.27–0.45, whereas much lower values were reported based on the indirect effect (Hansen et al., 2004; Jamrozik et al., 2005; Norman et al. 2009). The above suggests that the main genetic factors affecting GL heritability are the fetus and, to a smaller extent, it’s mother (Nogalski, 2003). The fetus’ dominant effect on GL could be related to the induction of parturition. The fetal adrenal cortex secretes cortisol, increasing placental prostaglandin levels which initiate parturition (Stabenfeldt et al., 1992). Nadarajah et al. (1989) and Cervantes et al. (2009) reported significant genetic correlations between GL vs. CE and SB, whereas Hansen et al. (2004) observed a weak correlation between GL and other reproductive traits. Owing to the high significance of CE and SB values in the production of high-yielding dairy cattle, other traits contributing to the optimization of reproductive traits should also be studied.

The objective of this study was to determine the effect of selected factors on gestation length in dairy cattle and the correlations between gestation length and other reproductive traits.

MATERIALS AND METHODS

The study investigated 2,514 Black-and-White Polish Holstein-Friesian cows kept in four cattle sheds, including two sheds in a tied-stall system and two sheds in a free-stall system. Gestation length (GL) was determined as the number of days between effective fertilization until calving. The animals were divided into five classes based on the reported GL values: up to 275 days, 275 to 277 days, 278 to 280 days, 281 to 283 days, and above 283 days. Based on the age factor, the cows were divided into primiparous, second-calf and third-calf females. Three calving seasons were identified: fall-winter (1 October to 31 January), winter-spring (1 February to 31 May) and summer (1 June to 30 September). The variability in calving ease was determined based on direct observations of calving ease (CE), stillbirth rates (SB) and placental expulsion (PE). Twin pregnancies were not taken into account. Calving ease was described in four categories: 1 - unassisted calving, 2 - calving with little intervention (assisted by one person), 3 - calving with considerable intervention (assisted by several persons), 4 - very difficult calving (hard pull) requiring veterinary intervention. Stillbirth values were assigned the following scores: 1 - born alive, survived the first 24 h after birth, 2 - dead at birth or died within the first 24 h after birth. Placental expulsion was classified as follows: 1 - spontaneous expulsion within 12 h after birth, 2 - spontaneous expulsion within 24 h after birth, 3 - placenta removed by a veterinarian.

Distribution (%) of calving ease, stillbirths and placental expulsion in respect of investigated factors is presented in Table 2. Calf sex was described as: 1 - female calf, 2 - male calf. Calves were weighed directly after birth with the accuracy of 0.1 kg. They were divided into four weight groups: up to 34 kg, 35 to 37 kg, 38 to 40 kg, and above 40 kg.

Table 2

Distribution (%) of calving ease, stillbirths and placental expulsion in respect of investigated factors

Trait	Level	Calving ease				Stillbirths		Placental expulsion
		1	2	3	4	1	2	1	2	3
Calving season	Fall-winter	66.67	21.11	8.89	3.33	95.33	4.67	70.44	19.11	10.44
	Winter-spirng	55.51	24.25	13.43	6.81	92.99	7.01	62.53	22.04	15.43
	Summer	62.66	24.03	7.14	6.17	92.53	7.47	64.29	22.73	12.99
Successive calving	First	45.25	30.28	15.85	8.63	91.90	8.10	66.73	25.00	8.27
	Second	68.88	20.82	7.09	3.20	94.51	5.49	65.90	17.39	16.70
	Third	84.13	10.71	3.17	1.98	96.43	3.57	63.49	19.05	17.46
Calf sex	1 (heifer calves)	63.34	23.99	8.95	3.72	95.61	4.39	70.10	21.28	8.61
	2 (bull calves)	59.40	22.26	11.43	6.92	92.03	7.97	61.95	21.05	16.99
Husbandry system	1 (tie-stall)	54.82	26.76	11.37	7.05	92.23	7.77	70.36	19.57	10.07
	2 (free-stall)	69.22	18.51	8.90	3.38	95.55	4.45	60.14	23.13	16.73

The effect of GL on CE was evaluated using the x2 test. We wanted to examine whether the length of pregnancy has an impact on the share of dystocia in dairy cattle.

where:

fi = numbers obtained

Fi = expected number

The observed CE, SB and PE values were not characterized by normal distribution, therefore, they were subjected to probit transformation (Żuk, 1989) prior to statistical analysis.

All investigated traits, including gestation length, were subjected to explorative analysis with the involvement of multiple analysis of variance and the GLM procedure (SAS, 2008). A general linear model was developed in view of the key investigated traits as well as first-order interactions. The following classification model was used:

where:

yijklmn = phenotypic value of trait,

μ = general average,

Hi = fixed herd effect (i = 1, .., 4),

Wj = effect of successive calving (j = 1, .., 3),

Sk = effect of calving season (k = 1, .., 3),

Ul = effect of housing system (l = 1, 2),

Pm = effect of calf sex (m = 1, 2), (H×W)ij,, .., (U×P)lm, = interactions between factors, eijklmn = random error.

The genetic parameters of gestation length and reproductive traits were estimated using Gibbs sampling method and the sire model. The components of (co-)variance were determined with the involvement of THRGIBBS1F90 (Tsuruta and Misztal, 2006) software and bivariate linear-threshold models. GL was an obligatory component in bivariate models, and the second trait was one of the investigated reproductive attributes: SB, CE and PE. The components of (co-)variance were estimated by Gibbs sampling to generate 50,000 samples, of which 10,000 were representative of the burn-in period. The number of rejected initial samples was computed with the use of POSTGIBBS1F90 (Tsuruta and Misztal, 2006) application. The results produced by neighbor samples were correlated (auto-correlated), and genetic parameters were determined based on data for every 100th sample. Standard errors of (co-)variance components were calculated as standard deviation of component values produced by successive samples.

Other statistical sources of variability were determined in the GLM procedure, subject to the analyzed dependent variable, therefore, varied linear models were applied:

y = n×1 observation vector; l = vector of liabilities;

βH, βW, βS, βU, βP = fixed effect vectors: herd (pH×1), successive calving (pW×1), calving season (pS×1), housing system (pU×1), calf sex (pP×1);

s = q×1 vector of random additive genetic sire effects;

XH, XW, XS, XU, XP = incidence matrices for fixed effects: herd (n×pH), successive calving (n×pW), calving season (n×pS), housing system (n×pU), calf sex (n×pP);

Zs = n×q incidence matrix for random additive genetic sire effects;

e = n×1 vector of random errors.

It was assumed that the liabilities conditional on all of the effects were independent and normally distributed.

The studied population’s genetic parameters were determined using the following formula:

where, σ2 = additive genetic sire variance, σ2 = error variance, where: covs = genetic sire covariance.

Genetic parameters were estimated in two replications. At the first stage, 47 half-sib groups were created based on the sires of calved cows. In the second step, 24 half-sib groups were created in view of the sires mated to cows. The minimum size of a half-sib group was 15 individuals.

RESULTS AND DISCUSSION

The average gestation length was 279.1 (±4.79) days (Table 1), and was shorter than determined for Holstein cattle by Nadarajah et al. (1989) 281.3 days, Silva et al. (1992) 280 days, and Norman et al. (2009) 281.6 days. In a study by Hansen et al. (2004), the average gestation length in Danish Holstein cows was 278.5 days. Johanson and Berger (2003) reported average GL values in American HF cattle at 277.9 days. The gestation length of Guernsey, Holstein and Jersey cows in several Florida-based farms increased by four days on average over a period of 50 years. According to the authors Silva et al. (1992), this increase could be attributed to a steep rise in production levels. Shorter gestation was noted in the summer months, and these findings are consistent with the results of previous research which validated the effect of higher temperature and dietary changes on shorter GL (Hansen et al., 2004; Przysucha and Grodzki, 2009). In this study, a longer gestation period was observed in respect of older cows (Table 3). There was a significant difference between the first and the third calving, at one day on average. In a study analyzing only primiparous cows, Hansen et al. (2004) reported a difference of 3.5 days between heifers calved at the age of 22 and 36 months. In this study, the increase in calf weight significantly contributed to a rise in GL values. In the group of the heaviest calves (>40 kg), GL was 3.9 days longer on average in comparison with the group of the lightest calves (≤34 kg). Nadarajah et al. (1989) reported a difference of 3.7 days between cows with large and small fetuses. In this study, GL values were also significantly affected by the sex of the fetus, and they were 1.8 days longer for male fetuses. According to Silva et al. (1992) and Hansen et al. (2004), GL was 1.1 days longer for male calves.

Table 1

Statistical characteristic for investigated traits

Variable	Mean	SD	Skewness
Gestation lenght (d)	279.54	4.79	−0.029
Calving ease	1.60	0.88	1.36
Stillbirths	1.06	0.24	3.61
Placental expulsion	1.47	0.72	1.17

Table 3

Least square means (LSM) of gestation length subject to calving season, calving number, calf weight and sex and housing system

Specification	Gestation length (d)
	n	LSM	SD
Calving season
Fall-winter (1 October-31 January)	900	279.1	4.65
Winter-spring (1 February–31 May)	998	279.3	4.84
Summer (1 June–30 September)	616	278.5	4.88
Successive calving
First	1,136	278.6a	4.89
Second	874	279.3	4.74
Third	504	279.6b	4.55
Calf weight (kg)
≤34	626	278.0A	4.37
35–37	670	279.2B	4.28
38–40	558	280.7C	4.60
>40	300	281.9D	4.72
Calf sex
Heifer calves	1,184	278.1A	4.39
Bull calves	1,330	279.9B	4.97
Housing system
Tie-stall	1,390	279.0	4.72
Free-stall	1,124	279.1	4.87

Values followed by different letters in columns within factors differ significantly: capital letters p≤0.01, small letters p≤0.05.

Gestation length increased with the degree of calving difficulty (Table 4). Difficult and very difficult calvings were characterized by a significant increase in GL values. A significant correlation between GL and CE was demonstrated. Difficult and very difficult calvings that required veterinary assistance had a 15.67% share of all births. The highest share of difficult and very difficult calvings was reported in the group where GL exceeded 283 days, whereas the easiest calvings were observed in the group of cows whose GL ranged between 275 and 277 days. In a study by Cervantes et al. (2009), the GL of cows experiencing calving difficulties was 2.1 longer in comparison with cows that calved without assistance. Philipsson (1976) demonstrated a correlation between difficult calving and higher GL values in Swedish Friesian heifers.

Table 4

Calving ease vs. gestation length (x2 = 35.1; p≤0.01)

Calving ease	Gestation length (d)
	<275		275–277		278–280		281–283		>283
	n	%	n	%	n	%	n	%	n	%
Unassisted	254A	55.95	408B	73.12	344B	66.70	328B	62.12	234A	51.09
Little intervention	134a	29.52	88b	15.77	90b	22.87	122a	23.11	118a	25.76
Difficult	44	9.69	38	6.81	54	10.47	46	8.71	76	16.59
Very difficult	22	4.85	24	4.30	28	5.43	32	6.06	30	6.55

Values followed by different letters in rows differ significantly: capital letters p≤0.01; small letters p≤0.05.

Both prolonged and shortened gestation contributed to a significant increase in the number of stillbirths (Figure 1). The lowest SB rate (3.23%) was observed in respect of GL values in the range of 275–277 days, whereas the share of calves that were born dead and died within 24 h after birth increased three-fold (to 10.48%) when GL exceeded 283 days. Johanson and Berger (2003) and Philipsson (1976) reported a non-linear correlation between GL and SB. In this study, optimal gestation length in Holstein cows was determined in the range of 275–277 days based on calving ease and stillbirth values. In a study evaluating more than 4,000 Friesian cows characterized by average gestation length of 278.7 days, Johanson and Berger (2003) found that the optimal GL that minimized the risk of SB was 282 days. Deviations from the optimal GL value led to a significant increase in SB rates. According to Norman et al. (2011), average GL values optimized lifetime productivity, calving ease, stillbirth rates and calving-to-first-service interval.

Figure 1

Gestation length vs. perinatal mortality and unassisted calving.

In the discussed experiment, genetic parameters were estimated with the use of Gibbs sampling and bivariate linear-threshold models. A similar approach was adopted by Lee et al. (2002), Ramirez-Valverde et al. (2001), Varona et al. (1999) who investigated the adequacy of threshold-linear models for estimating breeding value based on calving ease, birth weight and gestation length. The above authors concluded that linear-linear models and threshold-linear models produced a similar sequence of breeding values for birth weight and gestation length. A lower correlation was determined with regard to calving ease (threshold trait). Threshold-linear models delivered more accurate estimations than linear- linear models.

The heritability of gestation length that accounted for direct sire of fetus effects varied in the range of 0.201–0.210 (Table 6). Heritability values determined in view of indirect maternal effects were visibly lower at 0.054–0.073 (Table 5). GL heritability indicators reported in this study were generally lower than the values reported by other authors (Hansen et al., 2004; Crews, 2006; Norman et al., 2009; Johanson et al., 2011).

Table 6

Heritability, genetic and phenotypic (SE) correlation between gestation length and calving ease subject to direct sire effects

Parameter	Gestation length	Stillbirths	Placental expulsion	Calving ease
Heritability	0.201–0.210 (0.062)	0.195 (0.099)	0.560 (0.182)	0.552 (0.169)
Genetic correlation	-	0.179 (0.453)	0.149 (0.244)	0.242 (0.161)
Phenotypic correlation	-	0.078 (0.044)	0.046 (0.029)	0.113 (0.032)

() standard deviation is shown in brackets.

Table 5

Heritability, genetic and phenotypic (SE) correlation between gestation length and calving ease subject to indirect sire effects

Parameter	Gestation length	Stillbirths	Placental expulsion	Calving ease
Heritability	0.054–0.073 (0.032)	0.420 (0.157)	0.155 (0.067)	0.211 (0.085)
Genetic correlation	-	0.555 (0.321)	−0.340 (0.415)	0.625 (0.294)
Phenotypic correlation	-	0.088 (0.051)	0.040 (0.025)	0.108 (0.030)

() standard deviation is shown in brackets.

The results of this study as well as the cited authors’ findings were characterized by significant variance which can be attributed to numerous factors, including the classification model (linear or threshold), type of random effects, influence of breed or successive lactation (Hansen et al., 2004; Crews 2006; Cervantes et al., 2009; Norman et al., 2009; Johanson et al., 2011).

The type of random effects in the classification model had a significant influence on the estimated heritability indicators as regards calving ease (Tables 5 and 6). When the random effect of sires of calved cows was taken into account, stillbirth heritability increased two-fold, whereas the heritability of placental expulsion was reduced more than three-fold, and the heritability of calving ease decreased more than two-fold.

The model accounting for the random effect of sires of calved cows demonstrates moderate positive genetic correlations between GL vs. SB and CE (Table 5). The genetic optimization of gestation length could have a negative influence on the breeding value of bulls by increasing stillbirth frequency and calving difficulty. The use of a classification model accounting for the sire of calf model produced correlations between GL vs. stillbirth rates and calving ease that were characterized by an identical pattern but lower strength (Table 6). Similar genetic correlations (direct effect) between gestation length vs. stillbirth rates (0.18) and calving difficulty (0.38) were reported by Hansen et al. (2004) with the use of a linear-threshold model. Subject to the type of random effect, genetic correlations with a different sign were reported for correlations between GL and placental expulsion. A negative correlation was determined for the random effect of sires of calved cows, while a positive correlation was reported for the effect of sires of live-born calves.

The coefficients of all phenotypic correlations between gestation length and calving ease were positive but low. The above suggests that prolonged gestation resulted in greater calving difficulty, higher stillbirth rates and longer placenta expulsion time.

An analysis of standard deviations of the investigated parameters revealed relatively high values (Tables 5 and 6). For this reason, modifications of gestation length should be introduced with great caution during attempts to improve calving ease and reduce stillbirth rates.

CONCLUSIONS

The following factors contributed to prolonged gestation in Holstein-Friesian cows: increasing cow’s age, male fetuses and growing fetus weight. Optimal gestation length was determined in the range of 275–277 days based on calving ease and stillbirth values. In view of direct sire of fetus effects, GL heritability was found to be four-fold higher in comparison with heritability values produced by indirect maternal effects, the model accounting for the random effect of sires of calved cows demonstrated moderate positive genetic correlations between GL vs. SB and CE. The above suggests that genetic optimization of gestation length could have a negative influence on the breeding value of bulls by increasing stillbirth frequency and calving difficulty. The use of a classification model accounting for the sire of calf model produced correlations between GL vs. SB and CE that were characterized by an identical pattern but lower strength. The investigated parameters were marked by relatively high standard deviations, implying that GL modifications should be introduced with great caution during attempts to improve CE and SB values.