|
|
||||||||



* Department of Animal Science, Agricultural University of Norway, Box 5025, N-1432 Ås, Norway
GENO Breeding and A.I. Association, Box 5025, N-1432 Ås, Norway
| ABSTRACT |
|---|
|
|
|---|
Key Words: female fertility nonreturn rate variance components genetic change
Abbreviation key: MSE = mean squared error, NRF = Norwegian Dairy Cattle, NR56D0 = 56-d nonreturn rate in virgin heifers
| INTRODUCTION |
|---|
|
|
|---|
The female fertility trait, nonreturn rate in virgin heifers, has been selected for in Norwegian Dairy Cattle (NRF) since 1972. The trait has received considerable weight in the total merit index for sires (8% to 15%) over the entire period. Documentation of the effect of long-term selection for improved female fertility is therefore possible in NRF.
The outcome of an insemination depends on both female and male fertility. Maternal effects include the female ability to produce fertile eggs and further to develop the embryo into a viable calf. Direct effects can be partitioned into the effect of transmitting genes from service sire and cow to the embryo and quality of the semen. In France, Boichard and Manfredi (1994) found that the service sire contributed 0.8% of the total phenotypic variance for conception rate. Weigel and Rekaya (2000) reported linear model repeatability estimates for service sire effect on 60-d nonreturn rate between 0.3% and 0.5%.
The aims of this study were to investigate the quality of Norwegian fertility data and to validate alternative models to be used for estimating variance components and genetic change of 56-d nonreturn rate (NR56D0) in NRF virgin heifers.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Records of a total of 2,183,618 first inseminations in the period from September 1, 1978, to December 1, 2000, were available for this study. All first and second crop daughters (with first insemination in the period) of sires that were progeny tested after 1980 were included. Both the sire and service sire had to be present in the sire pedigree file, and a record was accepted if age at first insemination was between 10 and 31 mo. Records were discarded if the heifer was treated by a veterinarian for reproductive disorders or was culled within 56 d after first insemination. Double insemination was defined as a new insemination occurring within 0 to 5 d, and the records were excluded if different sires were used in a double insemination. The data editing process is described in Table 1
. Data used in the analysis contained 1,632,961 records from daughters of 2945 sires. A total of 3163 service sires were represented in the data. A sire pedigree file was built by tracing the pedigree of sires of heifers and service sires back as many generations as possible. A total of 3532 males, the oldest were born in 1940, were included in the pedigree file.
|
|
|
|
Model 1:
![]() | ([1]) |
where
| Yijklmno | = | observation of 56-d nonreturn (0 = return, 1 = nonreturn) of heifer 0;
| Ai | = | fixed effect of age i at first insemination, in weeks, with 92 classes ranging from 10 to 31 mo of age;
| Mj | = | fixed effect of calendar month j at first insemination in 12 classes;
| DIk | = | fixed effect of double insemination k in 2 classes;
| hyl | = | random effect of herd-year l;
| sm | = | random effect of sire m;
| ssn | = | random effect of service sire n; and
| eijklmno | = | random residual term.
|
Model 2:
![]() | ([2]) |
where
| HYl | = | fixed effect of herd-year class l (replacing hyl in [1]); and other terms defined in [1].
|
Model 3:
![]() | ([3]) |
where
| MYj | = | fixed effect of month-year class in 266 classes (replacing Mj in [1]); and other terms as defined in [1].
|
Models 4, 5, and 6 were similar to Models 1, 2, and 3, respectively, but without the effect of service sire.
The following (co)variance structures were assumed for the random effects:
![]() |
where
,
,
, and
were variance components for effects of herd-year, sire of heifer, service sire, and residual, respectively. A was the additive relationship matrix, and I was the identity matrix. In the models including service sire, the random effects of sire and service sire were assumed to be correlated with the common additive relationship matrix A. The following expectations were assumed for the (co)variances of sire of heifer and service sire effects (Willham, 1963; Van Vleck, 1978):
![]() |
These were used to solve additive direct (
), additive maternal (
) variance and their covariance (
d,m) components. REML estimates of the (co)variance components were obtained using the VCE4 program (Neumaier and Groeneveld, 1998). The PTAs, given the estimated (co)variance components, were calculated using iteration on data routines in the DMU package (Jensen and Madsen, 1994).
Model Validation
Validating estimates of genetic trend.
A method described by Boichard et al. (1995) was used to validate the six models. In this method, one assumes that successive evaluations of a sire have the same expectation, equal to the true breeding value, and display only random variation associated with new information. Five data sets were created by restricting year of first insemination from 1978 to 1981, 1985, 1989, 1993, and 1997, respectively. Successive sire evaluations were carried out with increasing amounts of information.
To test for bias in genetic change, the following regression model was used (Boichard et al., 1995):
![]() |
where
| u, v | = | vectors of sire evaluations from two releases, u is estimated with data of daughters in the first period and v depends both on the data used for u and on subsequent data;
| g | = | effect of birth year of sire with incidence matrix X;
| b | = | coefficient, variability changes between methods, here b = 1;
| t | = | known vector based on the number of additional daughters each year;
| ![]() | = | bias in the estimate of annual genetic trend; and
| e | = | random error term, with and W as a weighting matrix.
|
Weights and t were calculated as described by Boichard et al. (1995). A total of 303 sires in 1981, 801 sires in 1985, 1335 sires in 1989 and 1861 sires in 1993 were included in the weighted regression analyses.
Goodness of fit and predictive ability.
The six models were compared according to their ability to predict randomly excluded observations. From herd-year groups containing more than 2 heifers, 10,000 records were excluded at random. New (co)variance components and new solutions of effects for each of six models were re-estimated using the dataset without the 10,000 records. The predicted (
) values for the 10,000 records were estimated by using solutions for the different effects in the models. The observed (y) had the value 1 or 0 (NR56D0). The statistics of mean square error (MSE) and the correlation (r
,y) between y and
were used to evaluate the models:
![]() |
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The REML estimates of variance components from the six models were all small, and heritabilities varied from 1.2% to 1.4%. This is in agreement with Pedersen and Jensen (1996) and Hodel et al. (1995), who found heritability estimates of 0.8% and 1.1% for 56 and 90-d nonreturn rate in heifers, respectively. The three models including service sire allowed estimation of direct additive heritabilities (0.91.0%) and maternal additive hertitabilities (1.2%). The direct effect includes viability of the embryo and the effect of sperm quality, while the uterine environment, egg quality, and embryo vitality are effects in the maternal component of variance. Using a model including the effect of service sire and a similar (co)variance structure for 60 d nonreturn in heifers, Jansen (1986) found that the additive direct heritabilities were slightly larger (0.021) and the additive maternal heritabilities slightly smaller (0.007) than in the present study. The genetic correlation between sire and service sire (0.203 to 0.257) deviated from those of Jansen (1986), where the correlation was -0.290, and also from the genetic correlation between sire and service sire effects of -0.11 reported by Boichard and Manfredi (1994). This lack of consistency may be explained by standardization of frozen semen, which reduces service sire variance and also affects the covariance between sire and service sire effects. As discussed in Jansen (1986), different semen processing could influence the result of nonreturn rates in service sires. Kommisrud et al. (1996) found a difference of 1.2% (P < 0.01) in 60-d nonreturn between two different semen extenders.
In all six models, genetic change was positive with the steepest gradient obtained from models 1 and 4. Although the models resulted in different trends (Figure 4
), the model validation (Table 3
) did not reveal significant bias for any period. The method used was not able to detect bias in the evaluated models. This could be due to very small heritabilities, lack of bias, or both. The differences between breeding values estimated with a large group of first crop daughters and breeding values estimated with additional information from second crop daughters on the same sire were quite small.
Models 1, 3, 4, and 6 containing random herd-year effect, showed smaller MSE and higher correlation between observed and predicted values (Table 5
) than Model 2 and 5, which treated herd-year as fixed. This is due to small herd sizes in Norway, and because treating herd-year as random leads to improved utilization of the information compared to models with a fixed herd-year effect. Henderson (1973) and Ugarte (1992) both recommend modeling herd-year as fixed when sires are used nonrandomly over herds. However, semen of each NRF test sire is distributed and used at random throughout Norway. The use of elite sires in each herd is regulated, based on herd size and semen availability of the sire. Only a few doses of the highest ranking sires can be used in each herd. Therefore no bias from nonrandom use of sires across herds in this data should exist.
Similar needs as those reported by Henderson (1973) for separating additive breeding values from environmental effects exist when data are available over a long time span. Heringstad et al. (1999b) compared models with fixed versus random contemporary group effects in an analysis of clinical mastitis and concluded that the model must include a fixed structure to be able to separate environmental trend from genetic change. In this study, solutions of month-year with Model 6 resulted in different values for the same month between year (Figure 5
).
Although, models with service sire are nearer the true biological situation, the goodness-of-fit and predictive ability were not improved in models with service sire included compared to corresponding models without service sire (Table 5
). Models 4 and 6 showed almost equal validation results for goodness-of-fit and predictive ability. Because Model 6 includes effects for month-year of insemination while Model 4 includes only month of insemination, Model 6 is preferred irrespective of equal validation results for goodness-of-fit and predictive ability.
NR56D0 was analyzed with a linear sire model, although a threshold model would have been more appropriate for categorical data (Gianola, 1982). However, in view of the large daughter groups, making averages of 01 records close to normally distributed. Linear models have been shown to work well when binomal probabilities fall within the range of 20 to 80% (Van Vleck, 1971). The rank correlation between EBVs from linear and threshold models should therefore be high, as estimated by Weller and Ron (1992) for conception rate (> 0.99). Estimates of genetic change have been found slightly larger for threshold models (Heringstad et al., 2002), making results from a linear model more conservative.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Corresponding author:
I. M. A.-Ranberg; e-mail:
ina.ranberg{at}ihf.nlh.no.
Received for publication November 29, 2002. Accepted for publication February 23, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Heringstad, X.-L. Wu, and D. Gianola Inferring relationships between health and fertility in Norwegian Red cows using recursive models J Dairy Sci, April 1, 2009; 92(4): 1778 - 1784. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Konig, Y. M. Chang, U. U. v. Borstel, D. Gianola, and H. Simianer Genetic and Phenotypic Relationships Among Milk Urea Nitrogen, Fertility, and Milk Yield in Holstein Cows J Dairy Sci, November 1, 2008; 91(11): 4372 - 4382. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Walsh, F. Buckley, K. Pierce, N. Byrne, J. Patton, and P. Dillon Effects of Breed and Feeding System on Milk Production, Body Weight, Body Condition Score, Reproductive Performance, and Postpartum Ovarian Function J Dairy Sci, November 1, 2008; 91(11): 4401 - 4413. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. E. Pollott and M. P. Coffey The Effect of Genetic Merit and Production System on Dairy Cow Fertility, Measured Using Progesterone Profiles and On-Farm Recording J Dairy Sci, September 1, 2008; 91(9): 3649 - 3660. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Andonov, J. Odegard, I. A. Boman, M. Svendsen, I. J. Holme, T. Adnoy, V. Vukovic, and G. Klemetsdal Validation of Test-Day Models for Genetic Evaluation of Dairy Goats in Norway J Dairy Sci, October 1, 2007; 90(10): 4863 - 4871. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. David, L. Bodin, G. Lagriffoul, C. Leymarie, E. Manfredi, and C. Robert-Granie Genetic Analysis of Male and Female Fertility After Artificial Insemination in Sheep: Comparison of Single-Trait and Joint Models J Dairy Sci, August 1, 2007; 90(8): 3917 - 3923. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Banos, S. Brotherstone, and M. P. Coffey Prenatal Maternal Effects on Body Condition Score, Female Fertility, and Milk Yield of Dairy Cows J Dairy Sci, July 1, 2007; 90(7): 3490 - 3499. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T. Kuhn, J. L. Hutchison, and G. R. Wiggans Characterization of Holstein Heifer Fertility in the United States J Dairy Sci, December 1, 2006; 89(12): 4907 - 4920. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Heringstad, I. M. Andersen-Ranberg, Y. M. Chang, and D. Gianola Short communication: Genetic analysis of nonreturn rate and mastitis in first-lactation Norwegian Red cows. J Dairy Sci, November 1, 2006; 89(11): 4420 - 4423. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Heringstad, Y. M. Chang, I. M. Andersen-Ranberg, and D. Gianola Genetic analysis of number of mastitis cases and number of services to conception using a censored threshold model. J Dairy Sci, October 1, 2006; 89(10): 4042 - 4048. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Averill, R. Rekaya, and K. Weigel Random regression models for male and female fertility evaluation using longitudinal binary data. J Dairy Sci, September 1, 2006; 89(9): 3681 - 3689. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Yan, C. S. Mayne, T. W. J. Keady, and R. E. Agnew Effects of dairy cow genotype with two planes of nutrition on energy partitioning between milk and body tissue. J Dairy Sci, March 1, 2006; 89(3): 1031 - 1042. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. M. Chang, I. M. Andersen-Ranberg, B. Heringstad, D. Gianola, and G. Klemetsdal Bivariate Analysis of Number of Services to Conception and Days Open in Norwegian Red Using a Censored Threshold-Linear Model J Dairy Sci, February 1, 2006; 89(2): 772 - 778. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Gonzalez-Recio, Y. M. Chang, D. Gianola, and K. A. Weigel Number of Inseminations to Conception in Holstein Cows Using Censored Records and Time-Dependent Covariates J Dairy Sci, October 1, 2005; 88(10): 3655 - 3662. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Jamrozik, J. Fatehi, G. J. Kistemaker, and L. R. Schaeffer Estimates of Genetic Parameters for Canadian Holstein Female Reproduction Traits J Dairy Sci, June 1, 2005; 88(6): 2199 - 2208. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. A. Averill, R. Rekaya, and K. Weigel Genetic Analysis of Male and Female Fertility Using Longitudinal Binary Data J Dairy Sci, November 1, 2004; 87(11): 3947 - 3952. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. L. Muir, J. Fatehi, and L. R. Schaeffer Genetic Relationships Between Persistency and Reproductive Performance in First-Lactation Canadian Holsteins J Dairy Sci, September 1, 2004; 87(9): 3029 - 3037. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. M. Andersen-Ranberg, G. Klemetsdal, and B. Heringstad Short Communication: Validation of Two Animal Models for Estimation of Genetic Trends for Female Fertility in Norwegian Dairy Cattle J Dairy Sci, December 1, 2003; 86(12): 4090 - 4092. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |