|
|
||||||||
Sustainable Livestock Systems Group, Scottish Agricultural College, Bush Estate, Penicuik, Midlothian, EH26 0PH, United Kingdom
2 Corresponding author: gpollott{at}rvc.ac.uk
| ABSTRACT |
|---|
|
|
|---|
Key Words: genetic energy balance fertility progesterone
| INTRODUCTION |
|---|
|
|
|---|
The Langhill herd maintained by the Scottish Agricultural College at its Crichton Royal Farm (CRF) provides the opportunity to study the effects of both system and genetics on dairy traits because it is composed of 2 different levels of genetic merit with 2 contrasting systems of management. This herd is monitored for feed intake, condition scored, and weighed regularly, and feeds are sampled and analyzed routinely. In addition, milk progesterone profiles have been monitored in this herd and thus provide a more objective view of a cows reproductive physiology than do simple visible displays of estrus. This paper reports the results of a study that used both progesterone profiles and energy balance in early lactation to investigate differences in dairy cow fertility due to production system and genetic merit.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Genetic Lines
The 2 genetic lines were designated as control and selection lines. Cows were allocated to the same line as their dams and remained there throughout their productive life. Cow sires in the selection line were picked on the basis of their genetic merit for fat and protein PTA, available sires with the highest fat+protein kg of PTA being chosen at the time of AI. Sires of control line cows were selected to have the average genetic merit for fat+protein kg of UK animals at the time of AI. In both lines, matings were arranged such that the inbreeding coefficient was <6% at all matings. Cows from both lines were randomly allocated to production system groups to maintain balanced groups of cows. Cows were kept together and treated the same at all times, except where the production systems required management differences.
Production System Groups
Cows at CRF were randomly allocated to a high-concentrate or a high-forage production system at first calving in such a way as to keep the groups balanced for numbers. The high-forage system relied on homegrown feeds, including maize and other wholecrop cereals, and the cows were grazed on grass during the summer months. The winter ration consisted of grass silage, maize silage, and alkalage, at a ratio of 60:20:20 on a DM basis, plus a protein supplement. The ration was fed as a TMR. At least 75% of the DM of the ration was designed to come from forages. The target ME content was 11.5 MJ/kg of DM with a target CP content of 180 g/kg of DM. Forages were supplemented with a range of energy and protein sources to meet these targets.
The high-concentrate system cows were housed all year, with access to an exercise area during the summer months. Their ration also contained the forages mentioned above, in the same DM ratio to each other, with a supplement blend of energy and protein ingredients. The target ration ME content was 12.3 MJ/kg of DM with a target CP content of 185 g/kg of DM.
Breeding Policy
The breeding policy at CRF was designed to allow cows the maximum opportunity to express their fertility characteristics within recommended welfare guidelines. Heifers were initially mated at about 13 mo of age, depending on their live weight, to calve for the first time as close as possible to 24 mo of age. After calving, cows were bred at the first observed heat on or after d 42 of lactation. In all, a cow was permitted 7 services in a given breeding period before being removed from the herd. Cows failing to show heat after d 42 of lactation were treated by the veterinarian as appropriate for normal farm practice. No cow was artificially treated to induce estrus.
Traits and Recording
Traits Determined from the Milk Progesterone Profiles.
Fertility in dairy cows is commonly measured as a series of time-related events, or intervals between them, or both. However, the use of milk progesterone profiles allows an insight into one set of key underlying mechanisms. Milk samples were taken on Mondays, Wednesdays, and Fridays each week from cows in the herd calving from September 2003 to August 2005, during the first 140 d of the subsequent lactation. Milk progesterone analysis was carried out on whole milk by ELISA (Ridgeway Scientific, Alvington, UK) with a reading range of 1 to 15 ng/mL and milk progester-one levels recorded (ng/mL). Measurements below the lowest 1 ng/mL standard, which thus did not give a reading, were recorded as 0.1 ng/mL. The intra- and interassay coefficients of variation were 6 and 13%, respectively, and were worked out using quality control samples. At least 2 consecutive readings of
3 ng/mL were taken as evidence of luteal activity. A typical normal progesterone profile is illustrated in Figure 1
. This cow first exceeded the threshold level of progesterone of 3 ng/mL (Lamming and Bulman, 1976) on d 35, had a cycle with no service, had another cycle starting on d 53, was served, and had a positive pregnancy diagnosis on d 109. The following traits were calculated from the progesterone data using the methods outlined by Royal et al. (2000): interval from calving to first luteal activity (C-LA; d), average progesterone level (ng/mL), length of the luteal phase (d), interovulatory period (cycle length, d), and length of the interluteal period (d). In addition the number of cycles per lactation was counted.
|
Traits Recorded by, or Calculated from, Farm Observation.
Cows were milked 3 times per day and mostly fed and maintained indoors. In addition, 1 h per day was allocated for heat observation of the herd. This provided many opportunities for farm observation for signs of breeding activity and reduced the likelihood of missing signs of behavioral estrus. The days of heats, services, and calvings were recorded. This provided data from which to calculate the day of first observed heat (DFH), the day of first service, and in combination with subsequent calving information, the day of successful service. The number of heats and the number of services were counted from the data. Gestation length was calculated from the difference between the day of successful service and the next calving date, but any aborted calves were excluded from the data set. Calving interval was calculated from the dates of 2 consecutive calvings, again excluding any aborted calvings. All other interval traits were calculated from the above traits. Because of the breeding policy at CRF requiring first service to occur at the first heat on or after d 42 of lactation, several traits were recalculated from this starting point (see Tables 1 and 2![]()
). One additional trait, the incidence of silent heats, was calculated by combining the farm and cycle data. In this case, a silent heat was assumed to have occurred when a cycle, identified from the progesterone profiles, was not preceded by an observed heat. A silent heat was given a value of 1 and an observed heat the value 0 and the trait analyzed as a binary trait.
|
|
Traits Measured Once per Lactation.
Tables 1 and 2![]()
show 19 traits measured once per lactation. These were composed of basic descriptions of fertility (e.g., day of first heat), intervals between 2 fertility events (e.g., interval from first heat to first service), and traits relating to events post 42nd day of lactation. This latter group was analyzed because one of the management practices at CRF was to delay first service until the first heat on or after the 42nd day of lactation; this practice created a bias in some traits where animals were held back from naturally expressing their full potential for fertility.
The majority of these traits had skewed distributions and so they were transformed before analysis. Transformation was carried out by taking either the natural logarithm or square root of the trait, as determined by its distribution after transformation. In addition, 5 of the interval traits were also converted into binary traits. For example, the interval from first heat to successful service was coded as 0 if service occurred within 1 d of first heat and as 1 if it was greater than 1 d. This recoding to a binary trait (values 0 or 1) was carried out for intervals from first heat to first AI, first heat to successful service, first service to successful service, first heat >42 d to first AI, and first heat >42 d to successful service. In each case, the value 0 was assigned to intervals of less than 1 d. The use of transformations and binary trait analyses are shown in the results tables, as appropriate.
All traits shown in Tables 1 and 2![]()
were analyzed by REML mixed model analyses in SAS (2004). The binary transformations of the traits were analyzed using ASREML (Gilmour et al., 2002) with a logit-transformation applied to the binary data. The effects of genetic line, production system, lactation number, and year and month of calving were all fitted as fixed effects. The effect of cow within genetic line/production system group was fitted as a random effect. All first-order interactions between the fixed effects were also fitted in the model. The model used was
![]() | [1] |
where µ was the overall mean, Gi was the effect of the ith genetic line (i = control or selection line), Sj was the effect of the jth production system (j = high-forage or high-concentrate system), Lk was the effect of the kth lactation number (k = 1 to 5), Yl was the effect of the lth year of calving (l = 2003, 2004 or 2005), Mm was the effect of the mth month of calving (m = 1 to 12), An(ij) was the effect of the nth cow (n = 1 to 229), and eijklmno was the residual variance term associated with the o lactations (o = 1 to 363). The index values of the effects shown in model 1 were for C-LA, but some effects had fewer levels in some analyses where less records were available (see Tables 1 and 2![]()
).
The model was fitted using a step-down approach until only significant main effects remained in the model or main effects that were involved in a significant interaction. Results from the analysis were reported and effects considered to be significant when tested against the residual term using an F-test. The significance of the random cow effects was determined by carrying out a log-likelihood test on the likelihoods derived from 2 models, the full final model with and without the cow within genetic line/production system term included. The difference in likelihoods was tested against the chi-squared distribution with 1 df. Least squares means were computed from the full final model and converted back to real values where transformations had been used.
Traits Measured Once per Cycle.
Five traits that were characteristics of individual luteal cycles are shown in Table 3
. These traits were also logn transformed as appropriate. A similar method was used for analyzing the cycle traits as described above for the lactation traits with the addition of cycle number as a fixed effect. Also, month in these analyses referred to the month in which the cycle started. In this case, the model was
|
![]() | [2] |
and the terms were as described above with the addition of Cp, the pth cycle number (p = 1 to 6). In addition, the analysis of the incidence of silent heats included cycle length, mean progesterone level during the cycle, and days in milk as covariates.
Incidence of Cycle Abnormalities.
The incidence of the 4 cycle abnormalities were analyzed using a chi-squared test to see if there was any influence of genetic line or production system. The occurrence of an abnormality during a lactation was coded as 1 and a 2-way table of occurrence and group constructed. The table was tested against the chi-squared distribution with 1 df.
Energy Balance
Fertility is highly influenced by energy balance in early lactation (Butler and Smith, 1989; Beam and Butler, 1999; de Vries and Veerkamp, 2000). The 4 subgroups in this trial were characterized by large differences in milk yield and nutrient intake, and so the energy balance of the 4 groups would be expected to vary by design. To see if any fertility differences between groups could be explained by energy balance differences, the energy balance of each animal over the lactation was calculated from their daily live weight measurements and weekly BCS using the energy balance (EB2) methodology of Coffey et al. (2001). A typical energy content and energy balance profile is shown in Figure 2
. Several traits were calculated from the energy balance curve; day of return to positive energy balance, day of energy balance nadir, level of energy balance nadir, mean energy content of the animal over the first 25 d of lactation and the mean energy balance over the first 25 d of lactation. All these traits were added to model 1 as covariates and analyzed in a step-down approach using REML mixed model methodology in SAS. The early lactation traits C-LA and DFH were analyzed with this augmented model 1.
|
| RESULTS |
|---|
|
|
|---|
The Average Cow at CRF
Trait summaries are shown in Tables 1, 2, and 3![]()
![]()
. The 363 recorded lactations were split between the 2 genetic lines, with 191 control and 172 selection line lactations, and the 2 production systems, with 177 high-concentrate and 186 high-forage lactations. The mean lactation number was 2.25. Cows started cycling on d 32 of lactation and were observed to have their first heat on d 67. First service occurred on d 78, and cows were successfully served by d 120 with a gestation length of 282 d. Calving interval averaged 403 d. These fertility results were achieved with 4.08 luteal cycles (including the one that resulted in pregnancy), 2.63 heats, and 2.41 services. Farm policy delayed the day of first service by about 10 d on average. The average luteal cycle lasted 22.4 d and was composed of a 14.7-d luteal period and a 7.6-d interluteal period. About 61% of the completed luteal cycles were preceded by silent heats.
Analysis of Variance and Least Squares Means Summaries
A summary of the 11 analyses of variance derived from the fertility traits measured once per lactation is shown in Table 4
. Six of the traits were logn transformed, 3 were square-root transformed, and 2 were normally distributed. Five traits were unaffected by any of the effects fitted in the model: day of successful service, number of heats, number of services, calving interval, and number of heats occurring after d 42. In addition, 2 traits were only affected by the different cows used: day of first service and day of first heat after d 42. The calving to first luteal activity interval was affected by all effects except cow. Month of calving and lactation number only affected the calving to first luteal activity interval and no other traits.
|
|
An ANOVA summary for the 8 interval traits is shown in Table 5
. Lactation number did not affect any trait, and the effect of differences between cows was only found for the interval from C-LA to first heat (P < 0.05). Differences between the genetic groups was apparent for the intervals first heat to first AI (logn and binary analyses), first heat to successful service (logn and binary), and first heat after d 42 to first AI (logn). The least squares means in Table 7
show that control line cows had longer intervals than selection line cows for first heat to first AI (about 6 d), first heat to successful service (12.6 d), and first heat after d 42 to first AI (0.59 d). Control line cows also had lower conception rates than selection line cows, as shown by the binary analysis of first heat to first AI interval and first heat to successful service.
|
Three of the 5 luteal cycle traits were logn transformed (Table 6
); only mean progesterone level was normally distributed. Month of calving affected 4 of the luteal cycle traits, whereas only silent heat was affected by genetic group. In this case the selection line had a greater incidence of silent heats than the control line by about 6%.
|
Incidence of Luteal Cycle Abnormalities
The chi-squared analyses carried out on the incidence of luteal cycle abnormalities are summarized in Table 8
. There were significant differences between the genetic lines and production systems in the incidence of delayed ovulation postpartum (DOV1). The selection line and the high-concentrate system had greater than expected incidence. No difference in the incidence of subsequent delayed ovulations (DOV2) was observed between the genetic lines, but high-concentrate cows had a greater incidence of delayed ovulation between cycles than high-forage cows. No effect of genetic line or production system was observed on persistent corpora lutea in the first cycle (PCL1). Persistent corpora lutea in subsequent cycles (PCL2) were a feature of control line cows.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Hageman et al. (1991) report the fertility results from a trial using high and average production lines of Holsteins. They found differences between the 2 lines in the interval between first heat and first service and gestation length in maiden heifers, and between days open and calving interval in second-parity cows. In addition, they demonstrated significant positive linear regression coefficients for the following traits with milk yield: days to first heat, days open, and calving interval. This set of results also indicates lower levels of fertility with increasing milk yield. The high-production line had the lowest fertility. Reporting on the Langhill herd, Pryce et al. (1999) found differences between a high and average production line for the incidence of silent heats, the rate of conception to first service, calving interval, day of first heat, day of first service, and days open; they also found that the high-production line had the lowest fertility for all traits.
Other evidence for a link between genetic level and fertility can be seen from the relationship between bull index for fertility and UK production index value shown in Figure 3
. An increase in the production index of £1 was associated with a £0.05 reduction in the fertility index. The R2 of the relationship was small, however, indicating a wide range in fertility at any given production level. In addition, several authors have demonstrated both genetic and phenotypic correlations between production and fertility. For example, Wall et al. (2003) showed a range of unfavorable genetic correlations between 5 production traits and indexes with 4 fertility traits [calving interval, days to first service, nonreturn rate to 56 d (negative), and the number of inseminations per conception]. The genetic correlations ranged from 0.38 to 0.67.
|
The evidence reviewed above suggests a link between genetics and fertility in dairy cows, but several authors have highlighted the negative effect of low energy balance on fertility (Butler and Smith, 1989; Beam and Butler, 1999; de Vries and Veerkamp, 2000), and several components of energy balance are highly influenced by genetics. In the context of this study, high levels of milk production are associated with energy deficit in early lactation, which, in turn, is linked to lower levels of fertility. Thus, the relationship between genetics and fertility could be completely mediated through milk production level.
One key question addressed by the research reported here relates to whether the way in which the 2 genetic lines were formed at CRF has resulted in differences in fertility, and if so, what was the nature of this difference. Control line cows commenced luteal activity on d 25 of lactation, 5 d before selection line cows, and had their first observed heat on d 53, 14 d earlier than the selection line. This result was corroborated by the greater than expected incidence of delayed ovulation postpartum found in the luteal cycle abnormalities analysis for the selection line. Several authors have suggested that energy balance in early lactation influences the onset of cycling activity (e.g., de Vries and Veerkamp, 2000). After fitting several characteristics of early lactation energy balance to both C-LA and day of first heat in the present study, a difference was still present between the 2 lines of 5 d for C-LA and 10 d for day of first heat. This result suggests that there has been a change in the genes directly affecting fertility rather than just those affecting fertility via energy balance, milk yield, feed intake, and body fat changes indirectly.
An alternative way to look at these results may be in terms of the breed composition of the 2 lines. The current UK Holstein population is largely the result of grading up from British Friesian animals using imported Holstein semen (Simm, 1998). The breed makeup of the 2 lines may be one cause for the differences found in fertility. The mean proportion of Holstein genes in the cows used in this study was 95.7% for the control line and 96.9% for the selection line. This 1.2% is a very small difference and seems unlikely to be a major cause of differences between the 2 lines in fertility, although this difference tells us nothing about the actual genes and polymorphisms occurring in the 2 groups. Wall et al. (2005) estimated the difference between 0% (Friesian) and 100% Holstein animals as being 12.2 d for calving interval and 3.95 d for days to first service. Assuming a linear trend over the range of Holstein percentages of cows, a 1% difference would equate to 0.12 d in calving interval and 0.0395 d in days to first service, which does not explain the differences found here.
The farm at CRF has a policy of delaying first service until the first heat on or after d 42 of lactation. This policy had a greater effect on the control line than the selection line because more control line cows had their first heat before d 42 of lactation (34 vs. 15% for the selection line). The net result was a loss of 6.7 d in potential days to first service in the control line compared with the selection line. Thus, significant differences between the 2 lines in the number of cycles, interval from first heat to first AI, and interval from first heat to successful service were all due to this artificial constraint. However, the interval from first heat post d 42 to first AI was about 0.6 d shorter for selection line cows.
The genetic differences in gestation length were investigated by adding sex and birth weight to the model. When analyzed together, birth weight had a significant relationship with gestation length, but the genetic differences remained and there was no effect of sex. If sex was added to the model without birth weight, then there was a difference in gestation length between the sexes. Thus, the effect of sex on gestation length can be explained by differences in the mean weight of the sexes, but neither birth weight nor sex accounted for the differences between the 2 genetic lines.
Four of the 5 luteal cycle traits were totally unaffected by genetic line, and so these characteristics of reproduction appear to be invariant to large changes in production characteristics brought about by differing genotypes. However, there were differences between the 2 lines in the incidence of cycle abnormalities and the occurrence of silent heats. As noted above, there were differences in the incidence of delayed ovulation postpartum, and in addition, the control line cows had a higher than expected incidence of type 2 persistent corpora lutea (PCL2). The occurrence of these abnormalities results in longer intervals to observed heat and successful service, thus reducing fertility.
Production System
Energy balance is often quoted as having a key influence on fertility (Beam and Butler, 1999; Westwood et al., 2002). In early lactation, milk production and feed/energy intake are changing rapidly thus the energy balance may change rapidly too. The BCS is often quoted as an indicator of energy balance, but many CRF cows were still growing in their first and second lactation and so live weight needs to be considered as well. Thus the EB2 calculation of Coffey et al. (2001) was used to indicate energy status and model the changes in energy balance over the lactation.
Zurek et al. (1995) estimated that the number of days to the postpartum nadir of energy balance was related to C-LA such that cows resumed ovulation 15.4 ± 1.2 d after the nadir of energy balance. Butler and Smith (1989) suggested that the resumption of ovulation was related to the average energy balance during the first 20 d of lactation. There is thus a demonstrable link between energy dynamics and the resumption of fertility postpartum. However, Diskin et al. (2003) suggest that it is the return to positive energy balance that triggers ovulation.
Cows on the high-forage diet commenced luteal activity 6 d earlier than those on the high-concentrate diet, but there was no difference in the days to first heat between the 2 systems. This result was reflected in the luteal cycle abnormality analysis, which showed a higher than expected incidence of delayed ovulation postpartum for the high-concentrate cows. The differences between the 2 production systems for C-LA disappeared when characteristics of early-lactation energy balance were fitted in the model. This implies that the production system differences in C-LA were due to the effects of the 2 nutritional regimens through characteristics such as day of energy balance nadir, mean energy content, and mean energy balance over the first 25 d of lactation. The proportion of the 2 groups having their first heat before d 42 was similar (23 vs. 27% for high-concentrate and high-forage groups, respectively), and the mean loss of days to first service because of the farm service policy was 9.2 and 10.0 d, respectively.
The reproductive performance of the high-forage group was poorer than that of the high-concentrate group. The high-forage group required more luteal cycles to become pregnant and had poorer conception characteristics than the high-concentrate animals.
When birth weight and sex were added to the model analyzing gestation length, the difference between the production systems was still highly significant. Thus, as with the genetic lines, these 2 additional factors did not account for the difference in gestation length between the 2 systems, and so some other effect, or effects, must be causing this difference.
In contrast to the genetic line results, production system did have an effect on luteal cycle characteristics. High-forage cows had longer luteal phases in their cycles, by 0.9 d, but had shorter interluteal intervals, by 0.8 d. The net effect of these differences combined was a similar interovulatory period for the 2 production systems. The luteal cycle abnormality results show a higher than expected number of delayed ovulations type 2 (PCL2) for the high-concentrate cows.
| CONCLUSIONS |
|---|
|
|
|---|
Genetic selection for high productivity has been one reason for the decline in fertility in modern dairy herds, and a change in selection objectives to improve fertility is one way to reverse this trend. The incorporation of fertility characteristics in breeding value estimates and selection indices is one way to achieve this and has been implemented in several countries in recent years. Improving dairy cow diets in early lactation to encourage higher feed intakes and improve energy balance is another way forward. Alternatively, there is an argument toward less intensive systems and lower genetic merit for milk production to improve fertility. Any such approach requires a balanced economic view of the implications of these alternatives before a decision on the best way forward in any given situation is taken.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Received for publication December 4, 2007. Accepted for publication May 21, 2008.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. A. de Feu, A. C. O. Evans, P. Lonergan, and S. T. Butler The effect of dry period duration and dietary energy density on milk production, bioenergetic status, and postpartum ovarian function in Holstein-Friesian dairy cows J Dairy Sci, December 1, 2009; 92(12): 6011 - 6022. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. T. Garmo, A. D. Martin, E. Thuen, O. Havrevoll, H. Steinshamn, E. Prestlokken, A. Randby, M. Eknaes, A. Waldmann, and O. Reksen Characterization of progesterone profiles in fall-calving Norwegian Red cows J Dairy Sci, October 1, 2009; 92(10): 4919 - 4928. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |