|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ANIMAL GENETICS |


* Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg 24061;
and
Meat and Livestock Commission, Milton Keynes, UK; and
Sustainable Livestock Systems Group, Scottish Agricultural College, West Mains Road, Edinburgh EH9 3JG UK
| Abstract |
|---|
|
|
|---|
Key Words: connectedness genetic evaluation simulation bias sheep
| INTRODUCTION |
|---|
|
|
|---|
Despite these benefits, restrictions such as AI mating and commitment to use reference sires may limit producer participation in sire referencing schemes, at least for extended periods. Still, it is useful for producers to know the relative rank of their animals in the population for effective genetic improvement using outside sources of germplasm. If producers are not sufficiently connected to other management units there is risk of bias in across-unit genetic comparisons. Therefore, it would be useful to have a measure of the risk of potential bias when comparing animals across units.
Various connectedness measures have been proposed including the average prediction error of EBV differences across units (Kennedy and Trus, 1993
) and the correlation of prediction errors across units (Lewis et al., 1999
, 2005
). These measures are currently used in Canada and the United Kingdom to measure connectedness. However, the behavior of these measures over time and relative to biases due to different unit genetic means has not been well established. The objectives of this study were to use simulation to evaluate patterns of change over time in selected connectedness measures in sire referencing schemes and determine how effectively these measures quantify the level of risk due to potential bias in EBV comparisons across production units.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Simulation
Simulation Models.
Data sets were simulated using a modified version of the stochastic simulation of Lewis and Simm (2000)
. Specific parameters defining mating, reproductive, and survival characteristics were described by Kuehn et al. (2007)
. Fifteen flocks ranging in size from 40 to 140 breeding ewes were simulated, for a total of 1,030 breeding ewes per year. A single trait with a moderate within-flock heritability of 0.25, typical of growth and yield traits, was simulated for each animal. To create genetic differences among flocks, breeding values for base animals were obtained from normal distributions, with genetic means µj for flock j and a common within-flock additive variance (
2a) of 0.25. Values of µj for each flock were likewise sampled from a normal distribution with mean zero and variance 0.25. Breeding values and phenotypes for all animals were derived using the procedures of Kuehn et al. (2007)
.
Each flock underwent 10 yr of random selection followed by 15 yr of selection on the simulated performance trait. Rams selected from across the scheme as well as rams selected for use within their respective flocks were selected at random from the top one-sixth of their respective pools as ranked by BLUP EBV. Ewes still alive after their fourth lambing were culled. Enough ewe replacements were selected each year by BLUP EBV to maintain a constant flock size, resulting in an average replacement rate of 26%.
Simulation Scenarios. Connectedness was measured across flocks participating in several types of cooperative breeding schemes. Each scheme was replicated 50 times. Two schemes were used as baselines for comparison: (i) all flocks participating in an AI sire referencing scheme for all 15 yr of selection, or (ii) flocks that remained autonomous for all 15 yr of selection (no connections), with only within-flock selection. For the AI sire referencing scheme, 6 rams were chosen as a team from the 15 flocks and made available as reference sires. Half of the team was replaced each year. Each flock randomly selected 3 reference sires from the team and mated each sire to 10 ewes. Flocks were allowed to reuse 1 reference sire in consecutive years while participating in the scheme.
In addition to the baseline schemes, 4 other breeding schemes were simulated: 2 schemes involving termination of sire referencing after 5 or 10 yr, with the flocks remaining autonomous thereafter; flocks with sporadic participation in sire referencing; and natural service (as opposed to AI) sire referencing, which involved exchange of rams among flocks and a reference sire pool of 45, rather than 6, rams. Changes in connectedness and bias in estimation of breeding values across flocks were compared for these scenarios. Details on the specific assumptions of the sporadic and natural service sire referencing schemes can be found in Kuehn et al. (2007)
.
All 6 scenarios were simulated with 2 strategies for acquisition and use of rams other than the reference sires. In the first case, half of the flocks mated their excess ewes to unrelated purchased rams with no pedigree or performance data, and the rest of the flocks mated on average one-half of their excess ewes to purchased rams and one-half to their own homebred rams. Breeding values for the purchased rams were sampled from a normal distribution with mean 0.8
i and variance 0.25, where
i is the mean true breeding value of the overall scheme in year i. The value of 0.8 was chosen to reflect the genetic distance between flocks practicing across-flock vs. within-flock selection on BLUP EBV (Lewis and Simm, 2000
).
The first sire usage strategy reflects the level of purchased animal use in flocks that participate in sire referencing in the United Kingdom (Lewis and Simm, 2000
) and henceforth will be referred to as the outside sire strategy. In the second sire usage strategy, all flocks used only their own homebred sires in addition to the required reference sires. Although not likely in practice, this strategy provided a contrast to the outside sire strategy. This sire usage strategy will henceforth be referred to as the homebred only sire strategy. With these 2 sire use strategies (outside vs. homebred) and the 6 breeding schemes, 12 scenarios were explored.
The scenarios considered were the same as those used in Kuehn et al. (2007)
to evaluate the effectiveness of natural service and discontinuous sire referencing schemes at removing bias due to differing flock genetic means compared with AI sire referencing. Their results suggested that bias decreased at a decreasing rate as sire referencing schemes progressed. Bias was dramatically reduced after 5 yr with natural service and AI sire referencing, and after 10 yr with sporadic participation in the AI scheme. For this study, several connectedness measures were compared to determine whether the loss of bias could be predicted using these statistics.
EBV Prediction Model
An animal model with full relationships was used to derive EBV used for selection within the simulation scenarios. Fixed effects included flock, year, sex, and number of lambs born in a litter. In the outside sire strategies, fixed genetic groups, defined by birth year of sires introduced from outside the scheme, were also fitted. Further details on the model fitted are available in Kuehn et al. (2007)
.
Summary Statistics
Connectedness between groups of animals has traditionally been measured using functions of the inverse of the coefficient matrix (C; Foulley et al., 1992
; Kennedy and Trus, 1993
; Laloë, 1993
):
![]() |
where M is the fixed effect absorption matrix,
is the ratio of residual and additive variances (
2e/
2a), and A is the numerator relationship matrix. The inverse elements of C, multiplied by
2e,, are the prediction error (co)variances of the EBV.
For this study, 2 connectedness measures based on functions of prediction error were calculated for each year of selection in each simulated data set. The first measure was the connectedness correlation (Lewis et al., 1999
, 2005
):
![]() |
where PEC = prediction error covariance and PEV = prediction error variance; and ûi and ûj are EBV of animals i and j, respectively. The connectedness correlation was calculated and averaged each year for all pairs of ram lambs born in different flocks across all 15 flocks (lamb rij). This statistic was also calculated yearly on a flock mean basis:
![]() |
where ûi.(j.) is the mean breeding value of all animals recorded in flock i(j) since recording began. This flock connectedness correlation (flock rij) was shown by Lewis et al. (1999
, 2005)
to have a strong relationship with bias arising from different flock genetic means, with higher rij values corresponding with less bias. Flock rij is thus a measure of connectedness for all animals in the pedigree, whereas lamb rij is descriptive of the current generation.
The second connectedness statistic calculated was the scaled PEV of the difference (PEVD) in EBV between animals in different flocks:
![]() |
Like the connectedness correlation, this statistic was calculated each year both for all pairs of ram lambs from different flocks (lamb PEVD) and for the mean EBV of all animals in each flock since recording began (flock PEVD). Scaling the PEVD by the additive variance (
2a) allowed this statistic to be expressed without units. The average PEVD of EBV across management units was first suggested by Kennedy and Trus (1993)
as the most appropriate measure of connectedness because it measures the inherent risk in comparing these animals. Kennedy and Trus hypothesized that a lower average PEV corresponds to a reduced risk of comparisons across units.
The ideal connectedness statistic would measure the risk of potentially biased EBV comparisons across flocks due to different initial flock genetic means. Therefore, after the behavior of these statistics over time was assessed, candidate connectedness measures were plotted relative to the bias remaining in the system due to simulation of different initial flock genetic means. As in Lewis et al. (1999)
, to quantify the bias associated with comparisons of animals between flocks, a contrast (Lij) was calculated for all pairs of rams lambs produced in each year:
![]() |
where ui and uj are the true breeding values of animals i and j, respectively. Ram lambs were chosen for these contrasts because they were the candidates for future reference and homebred sire selection. The statistic was squared and averaged for pairs of ram lambs across flocks, resulting in a measure of the mean squared error (MSE) of prediction of differences in breeding values of ram lambs born in different flocks. The MSE is the sum of the PEVD and the squared bias in this predicted difference. In addition to the between-flock MSE (MSEB), the average squared Lij of pairs of lambs born within the same flock was also computed. This within-flock MSE (MSEW) was assumed free of bias since the lambs were born in the same flocks and, hence, would only consist of the PEV of EBV differences. Therefore, the difference between MSEB and MSEW estimates the squared bias due to different flock genetic means.
In their previous study, Kuehn et al. (2007)
found that the amount of bias in year zero of selection was dependent on whether outside sires were allowed or excluded. When both homebred and outside sires were used, MSEB averaged 0.50 at year zero; when no out-side sires were allowed, MSEB averaged 0.83. Across all years, MSEW remained at an average level of 0.27. Because the beginning level of bias in across-flock breeding value differences was strategy-dependent, a percentage decrease from the original level of bias was calculated for this study to compare with the connectedness measures as:
![]() |
where MSEB0 is the mean squared across-flock Lij in selection year 0 of the current simulation replicate and MSEBi is the mean squared Lij in selection year i.
| RESULTS |
|---|
|
|
|---|
Connectedness Correlations.
As selection and sire exchange across flocks began, the value of both of the correlation connectedness statistics (lamb rij and flock rij) increased in the sire referencing scheme scenarios (Figure 1
, Table 1
). When no sire referencing was undertaken, both lamb and flock rij necessarily remained at zero. Lamb and flock rij in continuous AI sire referencing and in AI sire referencing for 5 or 10 yr were expected to be the same after 5 yr of selection; the same was expected for continuous AI sire referencing with the AI sire referencing for 10 yr scenario after 10 yr of selection. The observed differences were within the standard error. The magnitudes of lamb and flock rij differed by a factor of approximately 10, indicating a greater level of confidence in estimates of difference in mean EBV between flocks than in estimates of difference in individual lamb EBV. This result was expected because means are less variable than individual predictions.
|
|
Both of the connectedness correlation measures increased approximately twice as fast when only homebred vs. some outside sires were used for nonreference sire matings. The introduction of nonpedigreed sires dilutes the average relationship between animals in different flocks, which results in a lower level of prediction error covariance between them and, therefore, lower connectedness correlations.
Prediction Error Differences.
As expected, both lamb and flock PEVD decreased as soon as group breeding schemes began (Figure 2
, Table 2
). If no scheme was implemented, lamb and flock PEVD slowly increased over all 15 yr, probably due to increased levels of inbreeding and average relationship within the independent flocks. Differences in lamb and flock PEVD in continuous AI sire referencing and AI sire referencing for 5 or 10 yr were expected to be the same after 5 yr of selection; the same was expected at yr 10 for the 10-yr scenario. Once again, differences between these scenarios were within the standard error.
|
|
Both lamb and flock PEVD values were higher at yr 0 if use of outside sires was not allowed. As reported by Kuehn et al. (2007)
, this discrepancy was likely due to the higher inbreeding in these scenarios (5.2 and 2.2%, respectively) after 10 yr of random mating but before selection. Inbreeding and within-flock relationships increase the additive variance between flocks; the upper limit of PEVD is proportional to the between-flock additive variance. Once breeding schemes were initiated, flock PEVD dropped more rapidly in scenarios with no outside sires, eventually converging with the flock PEVD from the other scenarios. As found when comparing flock to lamb rij, natural service sire referencing was no longer disadvantageous with regard to connectedness relative to AI sire referencing based on flock PEVD. As expected, declines in flock PEVD were slower with sporadic participation in the scheme, and flock PEVD increased slightly if flocks discontinued participation in the scheme altogether.
Connectedness Measures Relative to Bias
Kuehn et al. (2007)
illustrated that for all of these simulated breeding schemes, bias decreased at a decreasing rate as the schemes progressed. Their results show that the MSE of across-flock comparisons of lamb EBV asymptotically approach the MSE of the within-flock comparisons over time, indicating a reduction in bias in the across-flock genetic evaluations. Bias continued to decline over time and did not increase in cases where flocks discontinued sire referencing after 5 or 10 yr. Therefore, in choosing a candidate connectedness measure, it was important that the statistic was monotonic as years of selection increased. From the 4 measures described, flock rij and flock PEVD fit this criterion and were plotted against the percentage reduction in bias due to differing flock means in Figures 3
and 4
, respectively.
|
|
As shown in Figure 1
and Table 1
, if flocks disengaged from sire referencing after 5 or 10 yr, flock rij continued to increase even though new physical connections between the flocks were no longer established. This same pattern can be observed in Figure 3
, particularly when outside sires were not allowed. Bias no longer decreased after sire referencing stopped (Kuehn et al., 2007
) even though flock rij continued to increase, resulting in a disparity with the anticipated relationship between bias and flock rij.
The relationship between flock PEVD and bias is shown for all scenarios in Figure 4
. Points on the figure associated with lower values of flock PEVD coincided with more years of sire referencing. In general, as flock PEVD decreased, bias also decreased. However, unlike flock rij, the relationship between PEVD and bias is not consistent, especially between strategies where outside rams were allowed or barred. Bias was lower, whereas flock PEVD remained higher when only homebred nonreference sires were used. At the same flock PEVD, bias was higher if outside sires were allowed. The inconsistency of this pattern raises questions regarding use of flock PEVD as a diagnostic measure to manage risk in across-flock genetic evaluations.
| DISCUSSION |
|---|
|
|
|---|
Prediction error covariance between animals is zero if no connections exist. Adding connections through related animals (due to use of reference sires in this study) or through direct contemporary comparisons causes the prediction error covariance between animals in separate flocks to increase. This positive covariance reflects a correlation between errors of breeding value predictions for pairs of animals, and these errors are more likely to be in the same direction. With correlated errors, bias between flocks should decrease as prediction error covariance increases.
Evaluating connectedness from the perspective of the current generation (i.e., lamb measures) had advantages and disadvantages. Lamb measures of connectedness (lamb rij and lamb PEVD) were more sensitive to changes in the design of breeding schemes than their counterpart flock-based measures. Immediately after flocks disengaged from sire referencing, the connectedness, as measured by lamb statistics, decreased. This sudden response can be attributed to lamb connectedness being measured only on the current generation rather than using historic data captured by the flock rij and PEVD. Individual lambs are not, on average, highly related to one another across flocks and do not have enough information to produce highly accurate EBV in the year their performance was first recorded. These lamb measures therefore assess the precision of comparing lamb EBV across flocks in the current year of evaluation and could be used to rapidly detect changes in between-flock mating systems. However, they do not effectively address the rate of removal over time in bias associated with differences in initial flock genetic means.
Lamb measures are highly dependent on pedigree relationships among lambs in different flocks. When large numbers of half-sibs were born each year in different flocks through AI sire referencing, lamb rij was high relative to natural service sire referencing for which lambs in different flocks were at most cousins. However, flock rij increased as rapidly (if not faster) in natural service referencing scenarios as compared with AI sire referencing scenarios, suggesting that mean flock EBV are effectively compared with natural service schemes even though lambs are not as highly related. With natural service, individual flocks both contributed and used more rams from other member flocks in the scheme, thereby establishing strong relationships among lambs across years. As a caveat, this result may not be general. In this study, the rams selected as reference sires, both for AI and natural service, were chosen from the top one-sixth on EBV. Reference sire use was also heavy in these small flocks with a majority of ewes mated to reference sires; heterogeneity in flock genetic means declined quickly once sire referencing began. As a consequence, most flocks both used and contributed rams to the reference sire pool. This mutual exchange of rams leads to strong connectedness.
Measures of connectedness involving flock means are based on the mean EBV of all animals within each flock. Therefore, they are functions of the accuracy of prediction of flock genetic means, which are the source of bias we are concerned about in this study. Both flock rij and PEVD therefore assess how well the model predicts differences in flock genetic means. Flock rij measures correlations in errors of predicted genetic means and flock PEVD measures the precision with which differences in genetic means are estimated. Because these measures are based on information on whole flocks, they are less sensitive to changes in mating strategies between flocks and better reflect the monotonic decrease in bias observed by Kuehn et al. (2007)
and in other simulation studies (e.g., Hanocq et al., 1996
).
Flock rij was superior to flock PEVD in predicting patterns of change in percentage of bias explained over time. The general pattern in flock rij (Figure 3
) was consistent regardless of the type of breeding scheme and whether homebred rams alone were mated to ewes not bred to reference sires. Some of the noise observed when outside sires were allowed was due to the decreased bias in the system after the 10 yr of random mating that occurred before group breeding schemes were initiated. Because the outside sires were drawn from a population with a common mean, flock genetic means regressed during the 10 yr of random mating and part of the bias simulated between flocks decreased (Kuehn et al., 2007
). The amount of regression depended on the policy individual flocks adopted for use of outside rams, which varied between replicates of the simulation. This extra level of stochasticity, in addition to the sampling of flock genetic means, most likely accounts for the extra variance observed in Figure 3
for scenarios where outside sires usage was allowed.
The relationship of flock PEVD to bias, on the other hand, was highly dependent on the scenario, especially in terms of whether outside rams were allowed for nonreference sire matings. In all scenarios, bias decreased as flock PEVD decreased, but it was difficult to clearly define threshold values for PEVD to quantify the risk of comparing animals across flocks. The inconsistent pattern in PEVD relative to bias may partially be due to differences in inbreeding and average relationships within flocks in the 2 strategies used for mating surplus ewes. Because the relationship between animals within flocks must be higher when no nonpedigreed outside animals are brought into the system, the additive variance of mean comparisons across flocks should be higher, at least until connections are made between flocks. Although the values of the numerator relationship matrix were not available and could not be easily calculated with this level of replication, it may be useful to standardize flock PEVD using the average genetic relationship of animals between and within flocks as described by Laloë (1993)
.
The pattern for flock rij lends itself to possible bench-marks for management of risk due to different flock means. If 2 flocks achieve a flock rij level of 0.05, about 80% of bias due to differences in initial genetic means should be removed. Increasing flock rij to 0.10 should equate to a remaining level of bias equivalent to 10% of the initial level. These benchmarks can help producers decide whether further matings between flocks are desirable to further improve connectedness.
Establishing connections among flocks is important if the genetic means of flocks differ because biases in the genetic evaluation may otherwise occur. Connectedness measures can be increased with designed mating programs. The connectedness correlation of flock genetic means is the most consistent measure of the risk of bias in across-flock comparisons. Producers seeking to lower their risk when comparing the genetic merit of their animals to those in other flocks should enter into ram exchange programs and monitor connectedness using the flock connectedness correlation. Benchmark flock connectedness correlation values of 0.05 for good connectedness and 0.10 for superior connectedness will assist producers in deciding what level of connectedness is sufficient to permit across-flock selection decisions.
| Footnotes |
|---|
2 The authors thank G. Simm (Scottish Agricultural College, Edinburgh, UK) for helpful comments and British Sire Referencing Scheme members for their ongoing support of this research. ![]()
4 Current address: USDA-ARS, US Meat Animal Research Center, Clay Center, NE. ![]()
3 Corresponding author: larry.kuehn{at}ars.usda.gov
Received for publication May 7, 2007. Accepted for publication November 27, 2007.
| LITERATURE CITED |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. A. Kuehn, R. M. Lewis, and D. R. Notter Connectedness in Targhee and Suffolk flocks participating in the United States National Sheep Improvement Program J Anim Sci, February 1, 2009; 87(2): 507 - 515. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |