Comparison of methods for handling censored records in beef fertility data: field data

The purpose of this study was to compare methods for handling censored days to calving records in beef cattle data, and verify results of an earlier simulation study. Data were records from natural service matings of 33,176 first-calf females in Australian Angus herds. Three methods for handling censored records were evaluated. Censored records (records on noncalving females) were assigned penalty values on a within-contemporary group basis under the first method (DCPEN). Under the second method (DCSIM), censored records were drawn from their respective predictive truncated normal distributions, whereas censored records were deleted under the third method (DCMISS). Data were analyzed using a mixed linear model that included the fixed effects of contemporary group and sex of calf, linear and quadratic covariates for age at mating, and random effects of animal and residual error. A Bayesian approach via Gibbs sampling was used to estimate variance components and predict breeding values. Posterior means (PM) (SD) of additive genetic variance for DCPEN, DCSIM, and DCMISS were 22.6d2 (4.2d2), 26.1d2 (3.6d2), and 13.5d2 (2.9d2), respectively. The PM (SD) of residual variance for DCPEN, DCSIM, and DCMISS were 431.4d2 (5.0d2), 371.4d2 (4.5d2), and 262.2d2 (3.4d2), respectively. The PM (SD) of heritability for DCPEN, DCSIM, and DCMISS were 0.05 (0.01), 0.07 (0.01), and 0.05 (0.01), respectively. Simulating trait records for noncalving females resulted in similar heritability to the penalty method but lower residual variance. Pearson correlations between posterior means of animal effects for sires with more than 20 daughters with records were 0.99 between DCPEN and DCSIM, 0.77 between DCPEN and DCMISS, and 0.81 between DCSIM and DCMISS. Of the 424 sires ranked in the top 10% and bottom 10% of sires in DCPEN, 91% and 89%, respectively, were also ranked in the top 10% and bottom 10% in DCSIM. Little difference was observed between DCPEN and DCSIM for correlations between posterior means of animal effects for sires, indicating that no major reranking of sires would be expected. This finding suggests little difference between these two censored data handling techniques for use in genetic evaluation of days to calving.     

Fertility traits in spring-calving Aberdeen Angus cattle. 2. Model comparison

The aim of this study was to investigate the possible superiority of a threshold-linear (TL) approach for calving day (CD) and calving success (CS) analysis in beef cattle over 2 multiple-trait (MT), censored models, considering CD at the first 3 calving opportunities. The CD observations on animals that failed to calve in the latter models were defined as cows being assigned a penalty value of 21 d beyond the last observed CD record within contemporary group (PEN model) or censored CD values that were randomly obtained from a truncated normal distribution (CEN-model). In the TL model, CD records were treated as missing if a cow failed to calve, and parameters were estimated in a TL analysis including CS traits (TLMISS-model). The models included the effects of contemporary group (herd x year of calving x mating management), age at calving, physiological status at mating (lactating or nonlactating cow), animal additive genetic effects, and residual. Field data included 6,763 calving records obtained from first, second, and third parities of 3,442 spring-calving Uruguayan Aberdeen Angus cows. Models were contrasted using a data splitting technique, analyzing correlations between predicted breeding values (PBV) for each pair of subsamples, by rank correlations between PBV obtained with the different models, and by inspecting percentage of sires selected in common using the different approaches at 10 and 25% hypothetical percentages of animals selected. Breeding value correlations of CD between the subsamples for the TLMISS approach were greater (0.67 to 0.68) than correlations for the censored MT models (0.49 to 0.54). Average correlations between PBV of CD in 1 subsample obtained by CEN (PEN, TLMISS) and PBV of CS in the other subsample were -0.53 (-0.55, -0.60) in the first calving opportunity (CO), -0.54 (-0.58, -0.63) in the second CO, and -0.50 (-0.49, -0.58) in the third CO. Rank correlations between PBV for CD in PEN and CEN were high (0.93 to 0.97), but correlations of either method with PBV of CD in TLMISS ranged from 0.50 to 0.71. Common identification of bulls for the top 10% of sires (25% of sires), when selected with PEN/CEN models or the TLMISS model, varied between 50 (44%) and 60 (52%). The use of the TL animal model for genetic evaluation seems attractive for genetic evaluation of fertility traits in beef cattle.     

Comparison of homogeneity and heterogeneity of residual variance using random regression test-day models for first lactation Japanese Holstein cows

Using a large‐scale data set that included first lactation test day records from 1975 to 2000 for Japanese Holsteins, genetic parameters for milk yield were estimated by using random regression (RR) test‐day models (TDM) with heterogeneous and homogeneous residual variances. It is necessary for the RR‐TDM to include a function that explains the shape of the lactation curve. The RR‐TDM with the LW curve, which combined Wilmink's curve and a Legendre polynomial, was used for fitting the model for milk yield. In recent years, increases in residual variance have been noted for Japanese dairy cattle. Thus, three kinds of heterogeneous residual variance over the calving year were considered: H1, H2 and HG. Linear and quadratic exponential functions for the calving year were used in H1 and H2, respectively. Residual variance of HG was divided into five groups according to calving year. Homogeneous residual variance was HO. All heterogeneous residual variances increased with calving year in an almost linear fashion. Residual variance increased over the study period. However, there is no need to consider heterogeneous residual variances in genetic evaluations, because the heterogeneity of residual variance over the years did not affect the ranking of top sires and cows.     

Direct and maternal (co)variance components and genetic parameters for growth and reproductive traits in the Boran cattle in Kenya

Direct and maternal (co)variance components and genetic parameters were estimated for growth and reproductive traits in the Kenya Boran cattle fitting univariate animal models. Data consisted of records on 4502 animals from 81 sires and 1010 dams collected between 1989 and 2004. The average number of progeny per sire was 56. Direct heritability estimates for growth traits were 0.34, 0.12, 0.19, 0.08 and 0.14 for birth weight (BW), weaning weight (WW), 12-month weight (12W), 18-month weight (18W) and 24-month weight (24W), respectively. Maternal heritability increased from 0.14 at weaning to 0.34 at 12 months of age but reduced to 0.11 at 24 months of age. The maternal permanent environmental effect contributed 16%, 4% and 10% of the total phenotypic variance for WW, 12W and 18W, respectively. Direct-maternal genetic correlations were negative ranging from −0.14 to −0.58. The heritability estimates for reproductive traits were 0.04, 0.00, 0.15, 0.00 and 0.00 for age at first calving (AFC), calving interval in the first, second, and third parity, and pooled calving interval. Selection for growth traits should be practiced with caution since this may lead to a reduction in reproduction efficiency, and direct selection for reproductive traits may be hampered by their low heritability.     

Partitioning of herd,year and season variation in milk production

Variances in milk and fat yields due to herd-period-season effects (period either 6 or 12 months, and season either 1, 2 or 4 months) were partitioned fitting a three nested effects model with herd, period within herd and season within period within herd as random nested effects, sire a fixed effect and linear and quadratic regressions on age at calving. The overall effects of period and season of calving were regarded as fixed effects and were removed by precorrecting records using least squares constants estimated from the same data. Environmental correlations within herd, herd-period and herd-period-season were also estimated for all period and season groupings. Records of first lactation milk and fat yields on 25 158 progeny of 69 widely used proven Friesian-Holstein sires in 832 herds in England and Wales were used.The variance components for the effects of herd, period within herd, season within period within herd and residual accounted for 31, 6, 5 and 58% of the total variance in milk yield, and 35, 8, 7 and 50% of the variance in fat yield, respectively, using a period of 12 months and a season of 4 months. Differences amongst correlations within the same herd-period-season, fitting seasons of different lengths, were small. It was therefore concluded that shorter seasons in a herd-period-season fixed effect model of sire evaluation would be of no advantage.     

Estimation of genetic parameters for functional longevity in the South African Holstein cattle using a piecewise Weibull proportional hazards model

Non‐genetic factors influencing functional longevity and the heritability of the trait were estimated in South African Holsteins using a piecewise Weibull proportional hazards model. Data consisted of records of 161,222 of daughters of 2,051 sires calving between 1995 and 2013. The reference model included fixed time‐independent age at first calving and time‐dependent interactions involving lactation number, region, season and age of calving, within‐herd class of milk production, fat and protein content, class of annual variation in herd size and the random herd–year effect. Random sire and maternal grandsire effects were added to the model to estimate genetic parameters. The within‐lactation Weibull baseline hazards were assumed to change at 0, 270, 380 days and at drying date. Within‐herd milk production class had the largest contribution to the relative risk of culling. Relative culling risk increased with lower protein and fat per cent production classes and late age at first calving. Cows in large shrinking herds also had high relative risk of culling. The estimate of the sire genetic variance was 0.0472 ± 0.0017 giving a theoretical heritability estimate of 0.11 in the complete absence of censoring. Genetic trends indicated an overall decrease in functional longevity of 0.014 standard deviation from 1995 to 2007. There are opportunities for including the trait in the breeding objective for South African Holstein cattle.     

A prototype national cattle evaluation for sustained reproductive success in Hereford cattle

The objective of this research was to develop a prototype system for national cattle evaluation that would facilitate selection for improved fertility of daughters from Hereford sires. Raw data for this analysis were the birth dates of calves as reported by breeders to the American Hereford Association. Records from females entered this analysis with the reporting of a birth date for their first calf. At that time, females were required to be in contemporary groups of at least 3 animals and to have at least 2 additional paternal half-sibs also represented in the data. To explicitly define "sustained reproductive success," the philosophy taken was that a female that maintained a calving interval of 425 d or less would be considered successful. Females failing to meet this criterion were considered to be at the end of their successful lifetime. Data were analyzed using methodology for survival analysis with grouped data. Fixed contemporary groups were modeled as being time dependent, reflecting the females exposed for breeding in the same herd-year-season. Sire effects were time independent and considered random. Also included in the analysis were time-independent covariates for maternal weaning weight and total maternal calving ease from the national cattle evaluation of the American Hereford Association. Records from females still successfully in production at the time of this analysis, those that were transferred, those with calving intervals less than 280 d, and those that were successful until becoming donor dams for embryo transfer were considered censored. A total of 36,866 females contributed to this analysis, with 14,143 of these having censored records. The median number of females in a contemporary group was 6. A total of 3,323 sires had daughters with records. The median number of daughters per sire was 7. Heritability of sustained reproductive success on the underlying scale estimated from these data was approximately 0.05. Additional data accumulated over time will improve this genetic evaluation. Sustained reproductive success is important to the commercial beef industry, and results from this evaluation are expected to enhance the assessment of economic value of Hereford seedstock.     

Estimation of genetic parameters for milk and fertility traits within and between low,medium and high dairy production systems in Kenya to account for genotype-by-environment interaction

Dairy records from the Dairy Recording Service of Kenya were classified into low, medium and high production systems based on mean 305-day milk yield using the K-means clustering method. Milk and fertility records were then analysed to develop genetic evaluation systems accounting for genotype-by-environment interaction between the production systems. Data comprised 26,638 lactation yield, 3,505 fat yield, 9,235 age at first calving and 17,870 calving interval records from 12,631 cows which were descendants of 2,554 sires and 8,433 dams. An animal model was used to estimate variance components, genetic correlations and breeding values for the production systems. Variance components increased with production means, apart from genetic group variances, which decreased from the low to the high production system. Moderate heritabilities were estimated for milk traits (0.21–0.27) and fat traits (0.11–0.38). Low heritabilities were estimated for lactation length (0.04–0.10) and calving interval (0.03–0.06). Moderate heritabilities (0.25–0.26) were estimated for age at first calving, except under the high production system (0.05). Within production systems, lactation milk yield, 305-day milk yield and lactation length had high positive genetic correlations (0.52–0.96), while lactation milk yield and lactation length with age at first calving had negative genetic correlations. Milk yield and calving interval were positively correlated except under the low production system. The genetic correlations for lactation milk yield and 305-day milk yield between low and medium (0.48 ± 0.20 and 0.46 ± 0.21) and low and high production systems’ (0.74 ± 0.15 and 0.62 ± 0.17) were significantly lower than one. Milk yield in the low production system is, therefore, a genetically different trait. The low genetic correlations between the three production systems for most milk production and fertility traits suggested that sires should be selected based on progeny performance in the targeted production system.     

Genetic evaluation of the longevity of the Holstein population in Japan using a Weibull proportional hazard model

The objective of this study was to confirm the stability of the genetic estimation of longevity of the Holstein population in Japan. Data on the first 10 lactation periods were obtained from the Livestock Improvement Association of Japan. Longevity was defined as the number of days from first calving until culling or censoring. DATA1 and DATA2 included the survival records for the periods 1991–2003 and 1991–2005, respectively. The proportional hazard model included the effects of the region-parity-lactation stage-milk yield class, age at first calving, the herd–year–season, and sire. The heritabilities on an original scale of DATA1 and DATA2 were 0.119 and 0.123, respectively. The estimated transmitting abilities (ETAs) of young sires in DATA1 may have been underestimated, but coefficient δ, which indicated the bias of genetic trend between DATA1 and DATA2, was not significant. The regression coefficient of ETAs between DATA1 and DATA2 was very close to 1. The proportional hazard model could steadily estimate the ETA for longevity of the sires in Japan.     

Study of codes of disposal at different parities of Large White sows using a linear censored model

To study the genetic relationship between three grouped reasons for sow removal (SR) in consecutive parities, accounting for censoring, 13,838 records from Large White sows were analyzed. Data were from seven pure-line farms having, on average, 5.9% unknown SR. Three traits were subjectively defined, each corresponding to a classification of SR (reproductive [RR], nonreproductive [RN], and others [RO]). Records for each trait could take one of five categories, according to parity at removal (0 to 4 or later). A multivariate linear censored model was implemented. The model to estimate (co)variance components and parameters included the effects of year-season, region, contemporary group, and additive genetic effects. The most common SR was related to reproduction (48.5%). Diseases of different origin and cause, old age/parity, and sow death or loss accounted for about 18, 7, and 4% of total culls, respectively. Estimates of variance components showed heterogeneity of additive genetic and residual variances for the three traits. Estimates of heritability were 0.18, 0.13, and 0.15 for RR, RN, and RO, respectively. Genetic correlations between removal codes were high (> or =0.90). Results suggest sizeable additive genetic variances exist for parity at removal and different codes of removal. Different SR reasons seem to operate similarly or as a closely related genetic trait associated with fitness. In particular, RN and RO seem to be genetically indistinguishable. Data structure, definition, and volume are major limitations in studies of sow survival. A multiple-trait censored model is preferred to evaluate reasons of sow disposal. Grouped removal causes seem to be strongly genetically correlated but with heterogeneous variances, suggesting that combining all removal causes and treating the trait as parity at disposal is an alternative approach.     

Genotype by environment interaction for growth due to altitude in United States Angus cattle

The objectives of this study were to determine if sires perform consistently across altitude and to quantify the genetic relationship between growth and survival at differing altitudes. Data from the American Angus Association included weaning weight (WW) adjusted to 205 (n = 77,771) and yearling weight adjusted to 365 (n = 39,450) d of age from 77,771 purebred Angus cattle born in Colorado between 1972 and 2007. Postweaning gain (PWG) was calculated by subtracting adjusted WW from adjusted yearling weight. Altitude was assigned to each record based upon the zip code of each herd in the database. Records for WW and PWG were each split into 2 traits measured at low and high altitude, with the records from medium altitude removed from the data due to inconsistencies between growth performance and apparent culling rate. A binary trait, survival (SV), was defined to account for censored records at yearling for each altitude. It was assumed that, at high altitude, individuals missing a yearling weight either died or required relocation to a lower altitude predominantly due to brisket disease, a condition common at high altitude. Model 1 considered each WW and PWG measured at 2 altitudes as separate traits. Model 2 treated PWG and SV measured as separate traits due to altitude. Models included the effects of weaning contemporary group, age of dam, animal additive genetic effects, and residual. Maternal genetic and maternal permanent environmental effects were included for WW. Heritability estimates for WW in Model 1 were 0.28 and 0.26 and for PWG were 0.26 and 0.19 with greater values in low altitude. Genetic correlations between growth traits measured at different altitude were moderate in magnitude: 0.74 for WW and 0.76 for PWG and indicate possibility of reranking of sires across altitude. Maternal genetic correlation between WW at varying altitude of 0.75 also indicates these may be different traits. In Model 2, heritabilities were 0.14 and 0.27 for PWG and 0.36 and 0.47 for SV. Genetic correlation between PWG measured at low and high altitude was 0.68. Favorable genetic correlations were estimated between SV and PWG within and between altitudes, suggesting that calves with genetics for increased growth from weaning to yearling also have increased genetic potential for SV. Genetic evaluations of PWG in different altitudes should consider preselection of the data, by using a censoring trait, like survivability to yearling.     

Three statistical approaches for genetic analysis of disease resistance to vibriosis in Atlantic cod (Gadus morhua L.)

The objective of this study was to estimate and compare variance components and sire breeding values for disease resistance to vibriosis in Atlantic cod (Gadus morhua L.) using 3 statistical approaches. A total of 3,576 individually tagged juvenile cod from 50 full-sib families were infected with Vibrio anguillarum, which causes vibriosis, a frequently reported disease in cod aquaculture. The experimental fish were progeny of captured wild cod from populations of southern coastal cod (POP1), and northern coastal cod and northeast Arctic cod (combined as POP2 in the genetic analyses). Fish were randomly assigned to 1 of 3 test tanks, and daily mortality was recorded until the termination of the experiment at d 31 postinfection. Variance components were estimated separately for the 2 populations using a Cox regression model, univariate linear model, and a linear model that accounted for censoring. With all approaches, the additive genetic sire variance estimated from POP1 was greater than for POP2. Heritability estimates across populations varied from 0.08 to 0.17 depending on the method used. The Cox regression model and univariate linear model resulted in greater heritability estimates for POP1 (0.10 and 0.16) than for POP2 (0.08 and 0.13), whereas the contrary was true with a linear model that accounted for censoring (0.17 vs. 0.14). The predicted breeding values for the sires from the 3 approaches were highly correlated (0.97 to 0.99). This is likely due to the fact that censoring only occurred at the end of the test; i.e., observations of the most resistant fish were censored. The considerable genetic variation found in this study suggests that vibriosis resistance may be improved through selective breeding. The univariate linear model, even without censoring of the data, was robust for the estimation of breeding values using the present data. Therefore, inclusion of vibriosis resistance in the multivariate linear estimation of breeding values for the traits of economic importance in Atlantic cod seems appropriate.     

Bayesian inference for genetic parameter estimation on growth traits for Nelore cattle in Brazil,using the Gibbs sampler

This data set consisted of over 29 245 field records from 24 herds of registered Nelore cattle born between 1980 and 1993, with calves sires by 657 sires and 12 151 dams. The records were collected in south‐eastern and midwestern Brazil and animals were raised on pasture in a tropical climate. Three growth traits were included in these analyses: 205‐ (W205), 365‐ (W365) and 550‐day (W550) weight. The linear model included fixed effects for contemporary groups (herd‐year‐season‐sex) and age of dam at calving. The model also included random effects for direct genetic, maternal genetic and maternal permanent environmental (MPE) contributions to observations. The analyses were conducted using single‐trait and multiple‐trait animal models. Variance and covariance components were estimated by restricted maximum likelihood (REML) using a derivative‐free algorithm (DFREML) for multiple traits (MTDFREML). Bayesian inference was obtained by a multiple trait Gibbs sampling algorithm (GS) for (co)variance component inference in animal models (MTGSAM). Three different sets of prior distributions for the (co)variance components were used: flat, symmetric, and sharp. The shape parameters (ν) were 0, 5 and 9, respectively. The results suggested that the shape of the prior distributions did not affect the estimates of (co)variance components. From the REML analyses, for all traits, direct heritabilities obtained from single trait analyses were smaller than those obtained from bivariate analyses and by the GS method. Estimates of genetic correlations between direct and maternal effects obtained using REML were positive but very low, indicating that genetic selection programs should consider both components jointly. GS produced similar but slightly higher estimates of genetic parameters than REML, however, the greater robustness of GS makes it the method of choice for many applications.     

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

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

Genetic analysis of twinning rate and milk yield using a threshold‐linear model in Japanese Holsteins

The objective of this study was to estimate genetic parameters and breeding values for the twinning rate of the first three parities (T1, T2 and T3) and 305‐day milk yield in first lactation (MY), using a four‐trait threshold‐linear animal model in Japanese Holsteins. Data contained 1 323 946 cows calving between 1990 and 2007. Twinning was treated as a binary character: ‘single’ or ‘twin or more’. Reported T1, T2 and T3 were 0.70%, 2.87%, and 3.73%, respectively. Individual 305‐day milk yield was computed using a multiple trait prediction for cows with at least eight test‐day records. (Co)variance components were estimated via Gibbs sampling for randomly sampled subsets. Posterior means of heritabilities for T1, T2 and T3 were 0.11, 0.16 and 0.14, respectively. Genetic correlations between parities were 0.92 or greater. Genetic correlations of MY with twinning rate were not ‘significant’ (i.e. their 95% highest probability density intervals contained zeros). Multiple births at different parities were considered as the same genetic trait. The average evaluations of T1 (T2) for sires born before 1991 was 0.48% (2.25%) compared with a mean of 0.76% (3.37%) for sires born after 1992. A recent increase in the reported twinning rate reflects the positive genetic trend for sires in Japanese Holsteins.     

Restricted maximum likelihood estimation of additive genetic variance when selected base animals are considered fixed.

A method to estimate genetic parameters with a model that considers selected base animals as fixed was investigated. The model estimates genetic variance as a conditional variance based on the Mendelian sampling of gametes from the base parents. In a simulation study, 20 sires were selected and each was mated to 20 dams to create 400 animals for the next generation. Selection was for five generations, but only animals of Generations 4 and 5 were assumed to have performance records and known parents. Simulated values for additive genetic and residual variance were 10. Estimated genetic variance was 8.58 when base animals were assumed random and 6.03 when they were assumed fixed. Residual variance was overestimated in the latter case. When males of Generation 4 were not selected to have progeny, estimated genetic variance was 9.91. It was concluded that estimates for genetic parameters in a model with base animals assumed as fixed were not biased by selection of base animals, but a new bias was introduced if descendants of fixed base animals were selected. Estimation of genetic variance from dairy records of daughters of AI test bulls gave differences of up to 8% when the model removed bias from selected base animals.     

Genetic (co)variances for calving difficulty score in composite and parental populations of beef cattle: I. Calving difficulty score, birth weight, weaning weight, and postweaning gain

Heritability of 2-yr-old heifer calving difficulty score was estimated in nine purebred and three composite populations with a total of 5,986 calving difficulty scores from 520 sires and 388 maternal grandsires. Estimates were 0.43 for direct (calf) genetic effects and 0.23 for maternal (heifer) genetic effects. The correlation between direct and maternal effects was -0.26. Direct effects were strongly positively correlated with birth weight and moderately correlated with 200-d weight and postweaning gain. Smaller negative correlations of maternal calving difficulty with direct effects of birth weight, weaning weight, and postweaning gain were estimated. Calving difficulty was scored from 1 to 7. Predicted heritabilities using seven optimal scores were similar to those using four scores. The predicted heritability using only two categories was reduced 23%. Phenotypic and direct genetic variance increased with increasing average population calving difficulty score. The estimated direct and maternal heritabilities for 2-yr-old calving difficulty score were larger than many literature estimates. These estimates suggested substantial variance for direct and maternal genetic effects. The direct effects of 2-yr-old calving difficulty score seemed to be much more closely tied to birth weight than were maternal effects.     

Genetics of milk yield and fertility traits in Holstein-Friesian cattle on large-scale Kenyan farms

Purebred Holstein-Friesian cows are the main exotic breed used for milk production on large, medium, and small farms in Kenya. A study was undertaken on seven large-scale farms to investigate the genetic trends for milk production and fertility traits between 1986 and 1997 and the genetic relationships between the traits. This involved 3,185 records from 1,614 cows, the daughters of 253 sires. There was a positive trend in breeding value for 305-d milk yield of 12.9 kg/ yr and a drop in calving interval of 0.9 d/yr over the 11-yr period. Bulls from the United States (U.S.) had an average total milk yield breeding value 230 kg higher than the mean of all bulls used; Canada (+121 kg), Holland (+15 kg), the United Kingdom (U.K., 0 kg), and Kenya (-71 kg) were the other major suppliers of bulls. Average breeding values of bulls for calving interval by country of origin were -1.31 (Canada), -1.27 (Holland), -0.83 (U.S.), -0.63 (Kenya), and 0.68 d (U.K.). The genetic parameters for 305-d milk yield were 0.29 (heritability), 0.05 (permanent environment effect as proportion of phenotypic variance) resulting in an estimated repeatability of 0.34. Using complete lactation data rather than 305-d milk yield resulted in similar estimates of the genetic parameters. However, when lactation length was used as a covariate heritability was reduced to 0.25 and the permanent environment effect proportion increased to 0.09. There was little genetic control of either lactation length (heritability, 0.09) or calving interval (heritability, 0.05); however, there were strong genetic correlations between first lactation milk yield, calving interval, and age at first calving.     

Estimates of genetic parameters for carcass measures of body composition and growth in swine

Records for pigs included in an experiment on reciprocal recurrent selection conducted from 1956 through 1971 at the USDA Beltsville Agriculture Research Center were analyzed to obtain estimates of heritabilities and genetic correlations and to derive prediction equations for estimating weight of lean cuts (WTLC) and percentage of lean cuts of shrunk slaughter weight (LCPC). Lean cuts growth rate (LCGR) was then estimated as WTLC/age of pig at slaughter. The base population consisted of two unrelated crossbred strains. A total of 1,294 records of F1 and F2 crossbred pigs were analyzed with one barrow and one gilt from each litter. Estimates of heritabilities and genetic correlations were computed with sire components of variance and covariance from a nested analysis of variance with an assumed model of years, strain-lines within years, sire within strain-lines, dams within sires and residual. Degrees of freedom were 307 for sires in strain-lines, 270 for dams in sires and 646 for residual. Heritability (h2) estimates were .42 +/- .13, .41 +/- .13 and .27 +/- .18 for WTLC, LCPC and LCGR, respectively, and .71 +/- .16, .38 +/- .13, .31 +/- .13 and .25 +/- .15 for carcass length, average backfat thickness, longissimus muscle area and ADG in BW, respectively. These estimates were apparently the first published genetic estimates involving LCGR based on carcass data. It was recommended that prediction equations to estimate WTLC, LCPC and LCGR for use in swine testing programs be derived from current meat-type pigs.     

Analysis of gestation length in American Simmental cattle

Records of gestation length (71,461) for Simmental cattle were distributed with mean 284.3 d and standard deviation 5.52 d. Gestation length was found to increase with percent Simmental and was 1.9 d longer for calves born to mature dams than for those born to heifer dams. Bull calves experienced gestation lengths 1.5 d longer than heifer calves. Sire, maternal grandsire, residual and total variances were estimated to be 2.42, .58, 22.78 and 25.78 d2, respectively, by Henderson's Method III. Heritability of gestation length was calculated to be .374 from the sire variance and .09 from the maternal grandsire variance. Direct additive genetic variance was considered to be of greater importance than maternal additive genetic variance. Correlations between the evaluations of sires for gestation length and heifer calving ease, birth weight and weaning weight were .26, .26 and .13, respectively.