|
|
||||||||
1 Sustainable Livestock Systems Group, Scottish Agricultural College, Bush Estate, Penicuik, Midlothian, EH26 0PH, United Kingdom
2 School of Biological Sciences, University of Edinburgh, Ashworth Laboratories, Kings Buildings, Edinburgh, EH9 3JT, United Kingdom
3 Roslin Institute (Edinburgh), Roslin, Midlothian, EH25 9PS, United Kingdom
Corresponding author: E. Wall; e-mail: eileen.wall{at}sac.ac.uk.
| ABSTRACT |
|---|
|
|
|---|
Key Words: nonadditive genetic effect inbreeding fertility
Abbreviation key: CI = calving interval, DFS = number of days from calving to first insemination, F = inbreeding coefficient, F1 = first cross breeding, INS = number of inseminations per conception, MILK = daily milk yield at d 110, NR56 = nonreturn rate after 56 d
| INTRODUCTION |
|---|
|
|
|---|
Within a breed, the nonadditive genetic effects commonly expressed are associated with either inbreeding or crossbreeding. Inbreeding arises from the mating of related individuals and results in reduced heterozygosity, and when there is dominance, this reduction in heterozygosity is expected to lead to inbreeding depression. The management of inbreeding is becoming increasingly important in domestic livestock populations as genetic evaluation procedures, such as BLUP, are more accurate in identifying elite bulls and, when used in combination with truncation selection, increase the rate of inbreeding, largely caused by a greater tendency to select related animals. This increased potential for inbreeding has been exacerbated by the use of reproductive technology, allowing an increase in the selection intensity through the extensive worldwide use of semen and increasing the reproductive rate of elite females. The crossing of breeds increases the heterozygosity in the offspring, whereby crossbred progeny have a performance advantage over the midparent mean for that trait (Shull, 1914). Interbreeding the crosses or back-crossing breaks up epistatic gene combinations present in each breed but to different amounts in different crosses (Dickerson, 1969). The impact of these phenomena can be observed in heterosis and recombination loss, respectively. Many countries, including the UK, analyze Holstein and Friesian animals together. However, in the UK they are considered distinct breeds, and, therefore, the heterosis and recombination loss between the 2 breeds is accounted for in the genetic evaluation of production (Brotherstone and Hill, 1994).
The consequences of inbreeding include inbreeding depression, which is a reduction of the mean phenotypic value, particularly for traits connected with reproduction or fitness (Falconer, 1989). Inbreeding depression has been shown to decrease milk production by approximately 9 to 26 kg of milk per lactation for each 1% of inbreeding (Thompson et al., 2000a, b). Smith et al. (1998) estimated an economic loss in relative net income of approximately $22 to $25 (losses expressed in terms of production index) for registered cows per 1% increase in inbreeding over the lifetime of a cow. Very little data exist for nonproduction related traits, such as fertility, but it is expected that inbreeding depression could be substantially greater for such fitness-related traits than for production traits. Studies that have considered the effect of inbreeding on fertility and related traits have found a nonsignificant or small effect (Smith et al., 1998; Cassell et al., 2003).
Little work has been done to estimate inbreeding in the UK dairy population. Roughsedge et al. (1999) found the average inbreeding coefficient of cows born in 1997 to be 0.4% relative to a base population born in 1960. There has been a recent move away from Friesian genes to North American Holstein genes, with 76% of the cow population in 1997 having North American Holstein founder origins (Roughsedge et al., 1999). Although inbreeding was low at 0.4%, it can be hypothesized that the rate of inbreeding will follow that of the US population but with a generation lag. Kearney et al. (2004) showed that current levels of inbreeding in the UK Holstein population were 2.6% for females and 3.1% for males relative to a base population born in 1940, using a more complete pedigree file than that available to Roughsedge et al. (1999). Kearney et al. (2004) showed that the mean level of inbreeding is currently rising at a rate of 0.17%/yr. The average inbreeding coefficient of the US Holstein population is now approximately 5% with a current average annual increase of 0.2% (AIPL, USDA, 2004).
The UK dairy population was predominately Friesian until the 1980s (Roughsedge et al., 1999; Kearney et al., 2004), but a large influx of North American Holstein genes has occurred over the past 20 yr, resulting in a steep increase in the proportion of Holstein in the population. Therefore, historical data will include cows of varying proportions of Friesian and Holstein. However, a small population of "pure" and distinct Friesians still exist. Substantial crossing and upgrading has occurred in the UK dairy population, which still retains a reasonable spread of crosses. Estimates of heterosis effects for dairy cattle performance vary across studies and traits. For example, Ahlborn-Breier and Hohenboken (1991) described a heterotic effect of 6.1% for milk yield and 7.2% for fat yield in a Holstein-Jersey cross. Akbas et al. (1993) found a lower heterotic effect of 1.8 to 2.2% for 305-d milk, fat, and protein kg between Holsteins and Friesians.
Studies suggest that poor fertility has become a major reason for involuntary culling of dairy cows in the UK (Esslemont, 1993) and worldwide (Olori et al., 2002). Research has recently led to the development of a fertility index for dairy cattle in the UK to counter this decline (e.g., Wall et al., 2003). The index is based on data from calving, insemination records, milk yield, and BCS. Published fertility proofs in the UK are based on calving interval and nonreturn rate after 56 d weighted by their relative economic weights (independent of culling). The development of genetic analyses for fertility has created an opportunity for evaluation of the possible differential impact of nonadditive genetic variation for fertility, a trait directly associated with fitness, and production traits.
The purpose of this study was to estimate the effects of inbreeding, heterosis, recombination loss, and proportion of Holstein genes on dairy cow fertility and production and to examine their potential impact on selection decisions in the UK dairy population by comparing the outcomes of genetic evaluations made with and without these effects.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Production records for first-lactation Holstein-Friesian animals with at least 3 test days were taken from 1992 until the end of 2002. Validation and editing rules were applied to these data as described in Wall et al. (2003). The pedigree of all cows in the data set was extracted from the Holstein UK database. Animals with
4 complete generations of pedigree information were extracted, leaving 408,847 records. Only animals with this degree of pedigree completeness were considered for this analysis because of the problems in estimating accurate inbreeding coefficients for animals with incomplete pedigrees (Cassell et al., 2003).
Inbreeding coefficients were calculated for all animals using the algorithm of Meuwissen and Luo (1992). Each animal was also assigned to an inbreeding class of 0, 1, 2, ..., 10; 10 to 15; or 15; where the zero class included non-inbred animals, class 1 animals had inbreeding coefficients >0% but
1%, and so on, with class 15 consisting of animals with inbreeding of
15% . Percentage Holstein was calculated for all animals in this data set, based on the average percentage Holstein of the parents. Additional assumptions on percentage Holstein based on the breed code or country of origin and date of birth of a sire or dam were also applied. For example, all registered US animals were considered to be 100% Holstein, while the older UK bulls in the pedigree were considered to be 100% Friesian. Heterosis and recombination loss were calculated for all cows as follows (Akbas et al., 1993):
![]() |
where PS and PD are the proportion of Holstein for the sire and dam, respectively. Thus, (heterosis, recombination loss) took values (1, 0), (0.5, 0.25), (0.25, 0.188), (0.125, 0.109) for the first cross (F1) and the subsequent 3 backcrosses to Holstein, respectively.
More than 67% of cows had information on CI, over 27% had BCS information, and 88% had at least one insemination record; 63% had a record for INS. A total of 11,354 sires were included, all with
3 daughters in the data set that also represented 32,264 herd-year seasons. Table 1
gives summary statistics for the data set.
|
![]() |
![]() |
![]() |
where Pijk = an observation for CI, DFS, NR56, or INS; Tijk = MILK; Vijk = BCS; hysi = fixed effect of herd-by-year-by-season of calving interaction i; hysci = fixed effect of herd-by-year-by-season of visit interaction i on BCS; monthj = fixed effect of the month of calving i; ß1 to ß10 = linear and quadratic regression coefficients of the dependent variable (P, T, or V) on age of animal at calving (Xa), DIM at test (Xd_t), DIM at BCS measurement visit (Xd_c), inbreeding coefficient (XF), percent Holstein (X%), heterosis (Xhet), and recombination loss (Xrec) expressed as deviations from their mean; sirek = the random genetic effect of sire k; and eijk = residual random error term.
The significance of each of the solutions for each of the covariates was tested using a 2-way t-test. Fitting inbreeding as a simple covariate in the model assumes that the effect of inbreeding depression on a trait is linear. The curvilinearity of inbreeding depression was also tested by including inbreeding as a classification variable (inbreeding class, as defined in Figure 1
, added as a fixed effect in the previous models).
|
| RESULTS |
|---|
|
|
|---|
|
|
Heterosis was shown to have a small but statistically significant effect on MILK, CI, and DFS (Table 2
). Heterosis had a favorable effect on yield, with an estimated difference of 0.36 kg/d between an F1 and the average of Friesian and Holstein. However, this quantity corresponded to only 1.4% of the mean. Heterosis had a favorable effect on DFS and CI, reducing them by 1.1 and 1.6 d, respectively.
Recombination Loss
The mean recombination loss was 17% (SD 11%; Table 1
). The only statistically significant effect of recombination loss was upon MILK (Table 2
), resulting in a 1.45-kg loss in milk on test nearest d 110. The loss caused by recombination seen in MILK was larger, in absolute terms, than the gain achieved from the favorable heterosis. The joint impact on the F1 and the subsequent 3 backcrosses to Holstein was 0.355, 0.186, 0.185, and 0.112 kg, respectively, suggesting that while the F1 was marginally better than expected from the additive breed composition, the subsequent crosses were marginally worse.
Inbreeding
Average inbreeding was 1.7% (SD 2.0%; Table 1
), with nearly 94% of animals having an inbreeding coefficient of <5% (Figure 1
). Over 85% of cows in this data set were inbred to some degree relative to a base population of 1940. However, a small number of cows in this population had a high inbreeding coefficient (F), up to 38%.
Figure 2
shows the trend in average F and percent Holstein by year of birth of the 408,847 cows in the data set. A steady increase in F occurred over the period shown, and the current rate of increase in inbreeding of 0.2%/yr (based on linear regression) is equal to recent estimates in the current US population and in the complete UK Holstein population (Kearney et al., 2004). The average F of cows born in 2000 is 2.5% (Figure 2
), with nearly 98% of animals being inbred to some degree. The most recent (born since 1995) sires represented in this data set have a mean F of just under 3%, with a similar rate of increase as seen in their daughters.
Table 3
shows that inbreeding had a significant and unfavorable effect on all traits. The difference between a noninbred animal and an animal of F = 10% was a 2.8-d increase in CI, a 1.7-d increase in DFS, a 0.27-unit decrease in BCS, and a 0.54-kg decrease in MILK. This pattern of change is similar, but opposite, to the predicted heterosis in the F1 between the Friesian and Holstein. As for heterosis, the impacts of F on the traits directly related to fertility were not clearly larger than those on MILK and BCS.
|
8%). The deviations of the inbreeding class solutions from that estimated from linear regressions were all within one standard error up to the inbreeding class of 8%, suggesting that the effect of inbreeding on a trait is linear up to this point. The difference between animals in F class 0 (non-inbred) and F class 10 was +5 d in CI, +1 d in DFS, 0.6 in BCS, and 0.3 kg in MILK. These results were not significantly different from the result estimated from the linear regression, with the exception of BCS. The effect of inbreeding depression on all traits for animals with inbreeding <3% was unfavorable but small, with the effect of inbreeding being most pronounced for animals in higher inbreeding classes. For example, animals in the highest 3 inbreeding classes (inbreeding >9%) had a CI that was approximately 6 d longer than a non-inbred animal. Animals with inbreeding coefficients between 15 and 30% produced 1.16 kg less MILK compared with a noninbred animal, a 5% decline in the phenotypic mean for that trait. Inbreeding had a similar (and significant) effect on DFS to that had on CI, with highly inbred animals (
15%) having 4.5 additional d to first insemination compared with noninbred animals.
|
Wall et al. (2003) estimated the genetic decline in fertility and associated traits. Over the period from 1980 to 2000, a 4.9-d increase in CI and 3.6-d increase in DFS was observed. Distribution of inbreeding for sires born in 1980, 1990, and 2000 was obtained from the results of Kearney et al. (2004). Given these distributions, the effect of inbreeding on each trait at the 3 time points was examined. Over the 20-yr period, inbreeding depression accounted for a 0.6-d increase in CI, 11% of the genetic decline observed in that trait, while 8.3% of the genetic decline in DFS was attributed to inbreeding depression.
Effect of Nonadditive Genetic Effects on Breeding Value Estimation
Inbreeding was included as a covariate in the model for each trait. Heterosis was added to the models for MILK, CI, and DFS, while recombination loss was included in the model for MILK only. The effect of inbreeding was fitted as a covariate in the model because of how similar the inbreeding class solutions were to linear regression solutions up to 8%. Percentage Holstein was significant for all traits. However, this effect is already accounted for in the prediction of fertility indices by fitting genetic groups in the sire model with pedigree. Using both of these effects would result in some confounding.
Including nonadditive effects in the model, on average, caused a very slight and unfavorable change in the index and its components (Table 4
). The mean value for CI increases (0.16 d), and NR56 decreased. The standard deviation of the PTA for each trait was slightly larger with the inclusion of the nonadditive genetic effects in the model. However, there was little change in overall rank of animals by fitting these additional effects in the model. The rank correlations between the trait PTA calculated with the nonadditive effects fitted in the model and trait PTA calculated without the non-additive genetic effects fitted in the model were consistently >0.99. Although overall ranking of bulls remained relatively unchanged, the change in position of some individual bulls was considerable. Some bulls (out of 2000) dropped in ranking more than 90 positions, and others moved up more than 60 places when ranked on the fertility index.
|
| DISCUSSION |
|---|
|
|
|---|
Miglior et al. (1995) found that 10% inbreeding caused a 4% decrease in the phenotypic mean of total lactation yield. This analysis found a smaller effect of inbreeding depression on MILK, with a 2.3% decrease in the phenotypic mean for each 10% rise in inbreeding. Smith et al. (1998) found that inbreeding depression increased CI by 2.6 d per 10% rise in inbreeding, similar to the 2.8 d seen with this analysis. Cassell et al. (2003) found a nonsignificant effect of inbreeding on days to first insemination. However, Hoeschele (1991) found that each 10% rise in inbreeding caused a 1.3-d increase in days open, much closer to the significant 1-d increase in days to first insemination seen here. Weiner et al. (1992) showed that conception at first service in sheep declined by 4.2%, with a 10% rise in inbreeding, compared with the 0.6% increase in animals returning to first service observed in this study. In the current study, the impact of inbreeding on BCS (6.1% of the mean for a 10% rise in inbreeding), when expressed as a proportion of the phenotypic standard deviation (16.3%), was similar to results described by Falconer (1989) for BW in pigs (15%), litter size in mice (23%), and milk yield in cattle (17%).
The genetic trend in MILK from 1980 to 2000 was 2.7 kg in daily milk yield at d 110 (Wall et al., 2003). Inbreeding seemed to have a slightly favorable effect (nonsignificant) on MILK for animals in inbreeding classes 1 to 3, suggesting that inbreeding has improved MILK by 0.015 kg in the past 20 yr. A similar trend was seen for test day milk after d 70 in the American Holstein population (Thompson et al., 2000a). Figure 2
shows how inbreeding increased with time, as has genetic merit for production. Therefore, in the past 2 decades, cow inbreeding has increased as well as their genetic merit, and this could be a reason for the estimated favorable effect of low inbreeding on MILK. Modeling all of the nonadditive and additive effects (and trends) precisely may be difficult, and this problem could at least partially explain this apparent artifact resulting in a favorable effect of low inbreeding on production. However, if inbreeding in the UK continues to rise in line with the US, more animals will be in the higher inbreeding classes. Thus, in 10 yr (following current trends), inbreeding will have a negative effect on MILK, reducing it by 0.12 kg.
This study examined the effects of inbreeding in a population with relatively low levels of inbreeding and a historically low rate of inbreeding. Other studies have shown that as animals move into these higher inbreeding classes and the rate of inbreeding increases, the effects on overall performance are much larger than observed in this study. For example, Weiner et al. (1994) showed an almost linear decline in overall profitability in sheep (£1.27) with each percentage rise in F in rapid inbreeding.
Interestingly, the effects of Friesian genetics were favorable for all traits (3 to 7% better fertility than for the Holstein) except MILK, where the Holstein had a 13% benefit over the Friesian. The effect of heterosis MILK, CI, and DFS was low (0.5 to 1.4%). Other studies have found much higher effects of heterosis of the order of 9 to 20% for fertility and productivity traits (e.g., McAllister et al., 1994) but generally in crosses of more diverse breeds (e.g., Holstein x Ayrshire). The effect of recombination loss was as expected for MILK, with the favorable heterosis being lost because of recombination and, therefore, breakdown of the epistatic effects between the genes with additional crossing. Brotherstone and Hill (1994) saw a similar pattern attributable to crossing for lactation yields for 5 lactations; for example, heterosis resulted in a 100-kg increase of first lactation milk, with recombination resulting in a 156-kg loss. Some studies have shown a negative effect of recombination on fertility traits (e.g., Koenen et al., 1994). Distl et al. (1998) showed a favorable recombination effect on days open, which agrees with our results.
To date, no countries include inbreeding in the genetic evaluation model for production or fitness traits (INTERBULL, 2003). This may have to change if inbreeding continues to rise in dairy populations and the effects of inbreeding depression become more pronounced, but this may require a review of the presentation of breeding values. For example, if bull A has all daughters with F = 12.5% and bull B has daughters with F = 0%, accounting for inbreeding depression by adjusting all individuals to a common F (e.g., F = 0 or the mean of the population) would improve the estimated breeding value of bull A relative to B. However, this may be misleading if A is more related to the population, as the offspring of A will be expected to display greater inbreeding depression than those of B. This presentational problem can be resolved by reference to the definition of breeding value (Falconer, 1989) and, hence, transmitting ability. Thus, the transmitting ability of an individual is the difference between the mean performance of its offspring and the population mean, assuming random mating. With this definition, it is clear that an individual males PTA should be estimated at an F equal to its coefficient of coancestry with the population of breeding females. Thus, for presentation, each individual PTA is adjusted to a different value of F.
It is important that this adjustment for presenting breeding values is not confused with the separate task of managing genetic variation and the rate of inbreeding in the population. Studies have examined the optimization of contributions of selection candidates for maximizing genetic gain at a given rate of inbreeding in dairy cattle (Weigel and Lin, 2002; Kearney et al., 2004). This would have the most effect when combined with the selection of bulls and bull dams by breeding companies, since this is a) the subpopulation that is critical for managing the long-term gene flow and b) where the contributions of the bulls and dams of greatest impact can be most easily managed. Such decisions should be based upon the group coancestry of all breeding males and females, including coancestry among the breeding males and among the breeding females, and not upon the inbreeding coefficient of the offspring. This cannot be achieved by individual breeders, as they have little control over the total contributions of individual animals to the population and can only be achieved by the cooperation of worldwide breeding organizations. Some help can be given to breeders to minimize the impact of inbreeding depression in their herd. Breeding companies can recommend matings that avoid related individuals, which can be done by providing information on the expected inbreeding of daughters.
The differences between Holsteins and Friesians in the UK are accounted for in genetic evaluations of production traits by dividing founders in the pedigree into genetic groups and including heterosis and recombination loss coefficients in the model of analysis. This accounts for the benefits and costs in crossing the 2 breeds, as occurs frequently in the UK dairy population.
Crossbreeding of dairy breeds is practiced by a small proportion of dairy farmers in many countries across the world to capture the benefits of the heterosis in F1 populations and the epistatic interaction of genes of purebred animals. Ideally, this practice should be part of a population-wide controlled program of breeding, maintaining both pure and crossbred populations. However, implementation of genetic evaluation procedures that allow for the full exploitation of the advantages of crossbreeding may prove difficult, as it would require the estimation of and accounting for relevant genetic parameters of fixed breed, heterosis, and recombination effects for all crosses (Swan and Kinghorn, 1992). The estimates of the effect of heterosis will be different depending on the 2 purebreds in question. This study showed a small heterotic advantage of a crossbred progeny over that of the average of purebred Holsteins and Friesians. Brotherstone (personal communication, 2004) also found small heterotic effects for the majority of traits of interest and only found useful heterosis in SCS for various crosses in the UK. The estimation of each of the nonadditive effects between each potential cross of breeds in a country is a large task but essential to allow for the development of crossbreed and multibreed evaluations.
The approach to the presentation of PTA accounting for F (described earlier) could be extended to the other forms of nonadditive variation discussed in this paper, such as heterosis and recombination loss; however, the advisability of this is more dependent on the breeding system in which the multiple breeds are used. In a population that is subject to breed substitution and otherwise unmanaged interbreeding, it may be appropriate to consider presentation of breeding values adjusted to expected levels of heterosis or recombination loss arising from random mating. However, where crossbreeding systems are more tightly managed, this practice would be less inappropriate, and more direct predictions of crossbred performance would be more informative.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication March 5, 2004. Accepted for publication August 27, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. Waldmann, J. Hallander, F. Hoti, and M. J. Sillanpaa Efficient Markov Chain Monte Carlo Implementation of Bayesian Analysis of Additive and Dominance Genetic Variances in Noninbred Pedigrees Genetics, June 1, 2008; 179(2): 1101 - 1112. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Parland, J. F. Kearney, M. Rath, and D. P. Berry Inbreeding Effects on Milk Production, Calving Performance, Fertility, and Conformation in Irish Holstein-Friesians J Dairy Sci, September 1, 2007; 90(9): 4411 - 4419. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. D. Dechow, G. W. Rogers, J. B. Cooper, M. I. Phelps, and A. L. Mosholder Milk, Fat, Protein, Somatic Cell Score, and Days Open Among Holstein, Brown Swiss, and Their Crosses J Dairy Sci, July 1, 2007; 90(7): 3542 - 3549. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Cerchiaro, M. Cassandro, R. Dal Zotto, P. Carnier, and L. Gallo A Field Study on Fertility and Purity of Sex-Sorted Cattle Sperm J Dairy Sci, May 1, 2007; 90(5): 2538 - 2542. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Mc Parland, J. F. Kearney, M. Rath, and D. P. Berry Inbreeding trends and pedigree analysis of Irish dairy and beef cattle populations J Anim Sci, February 1, 2007; 85(2): 322 - 331. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Gulisija, D. Gianola, and K. A. Weigel Nonparametric Analysis of the Impact of Inbreeding on Production in Jersey Cows J Dairy Sci, January 1, 2007; 90(1): 493 - 500. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |