|
|
||||||||

* Scottish Agricultural College, West Mains Road, Edinburgh, EH9 3JG, Scotland, U.K. and
and
Roslin Institute, Roslin, Midlothian, EH25 9PS, Scotland, U.K.
| Abstract |
|---|
|
|
|---|
F) before and after the introduction of BLUP evaluations were not significantly different. The observed
F in the last generation was 1.0% for Meatlinc and 0.2% for Aberdeen Angus. The application of the dynamic selection algorithms for maximizing genetic gain at a fixed
F led to important expected increases in the rate of genetic gain (
G). When
F was restricted to the value observed in both populations, increments per year in
G of 4.6 (i.e., 17%) index units for Meatlinc and 3.5 (i.e., 30%) index units for Aberdeen Angus were found in comparison to the
G expected from conventional truncation BLUP selection. More relaxed constraints on
F allowed even higher expected increases in
G in both populations. This study demonstrates that the optimization tools constitute a potentially highly effective way of managing gain and inbreeding under a broad range of schemes in terms of scale and inbreeding level. No losses in genetic gain were associated with the use of dynamic optimization selection when schemes were compared at the same
F.
Key Words: Beef Cattle Genetic Gain Inbreeding Optimization Sheep
| Introduction |
|---|
|
|
|---|
Although inbreeding is unavoidable in closed selection programs, increases in inbreeding need to be restricted to alleviate long-term negative effects (Lamberson and Thomas, 1984
; Burrow, 1993
). Woolliams et al. (2002)
have described the rate of inbreeding as a measure of risk from the perspective of the breeding program justifying its management with arguments that go beyond avoiding inbreeding depression and loss of genetic variation in the selected trait.
Dynamic tools for maximizing genetic progress while constraining the rate of inbreeding to a predefined value are now available (Meuwissen, 1997
; Grundy et al., 1998
; Meuwissen and Sonesson, 1998
; Grundy et al., 2000
). These tools optimize the number of parents and their contributions to subsequent generations for maximizing gain for a fixed rate of inbreeding. Simulation studies showed improvements in genetic gain greater than 20% over BLUP truncation selection at the same rate of inbreeding (Meuwissen, 1997
; Grundy et al., 1998
). However, the expected benefits from optimized selection in real livestock populations remain unknown.
The main objective of this study was to evaluate the potential extra gains to be obtained by dynamic optimization algorithms in two livestock populations of sheep (Meatlinc) and beef cattle (Aberdeen Angus). This was accompanied by a description of the population structure and rates of genetic gain and inbreeding before and after the introduction of BLUP evaluations. Relationships between contributions of ancestors of the current populations and their EBV were also evaluated.
| Materials and Methods |
|---|
|
|
|---|
Data
The Aberdeen Angus is a traditional British beef breed, with a recorded pedigree extending over 50 yr. The Meatlinc is a synthetic terminal sire breed of sheep created in the United Kingdom in the early 1960s. In contrast with Aberdeen Angus, the recorded Meatlinc pedigree is relatively small, covering only 24 yr. Pedigree data and index scores for both populations were provided by the Meat and Livestock Commission (M.L.C., Milton Keynes, U.K.). The two indices provided were the BLUP-EBV for the breeding goals of U.K. terminal sire breeds of beef cattle and sheep ("beef value" and the "lean index," respectively). The beef value includes carcass weight, carcass conformation score, and carcass fat score (Amer et al., 1998
; Simm, 1998
), whereas the lean index includes carcass lean weight and carcass fat weight (Simm and Dingwall, 1989
).
The Aberdeen Angus pedigree included a total of 119,953 animals (57,431 males and 62,522 females) born from 1948 to 2000. A total of 45,472 parents (6,686 sires and 38,786 dams) were identified. The Meatlinc pedigree included a total of 12,391 animals (5,661 males and 6,730 females) born from 1974 to 2000. A total of 3,742 parents (329 rams and 3,413 ewes) were identified. Parents with unknown genealogies were considered as "base parents." This group represented 28.9% of the total number of parents (2,443 sires and 10,704 dams) in Aberdeen Angus and 7.4% of the total number of parents (35 rams and 243 ewes) in Meatlinc.
Because multitrait BLUP evaluations were introduced for both populations in 1991, the analyses performed were each applied to three periods of data. The three periods included an overall period covering all years with available information, and two periods of approximately equal length defined pre-BLUP and post-BLUP introduction.
Generation Intervals
The generation interval for each breed was computed as the average age of parents at the birth of their offspring. It was calculated for each year of birth and then averaged over years for all parents (L), for sires (Lm), and for dams (Lf).
Rates of Genetic Progress and Inbreeding
Average index scores and inbreeding coefficients of individuals born at each year were calculated. The inbreeding coefficients (F) were obtained from the additive relationship matrix that was computed using the algorithm of Meuwissen and Luo (1992)
. The rate of genetic gain (
G) and the rate of inbreeding (
F) were computed as the linear regression of the average index score and average F on the year of birth, respectively.
Both
G and
F were analyzed for the three periods in both populations. For Meatlinc, the periods were 1) overall period from 1974 to 2000, 2) pre-BLUP period from 1983 to 1991, and 3) post-BLUP period from 1992 to 2000. For Aberdeen Angus, the corresponding periods were 1) overall period from 1948 to 1999, 2) pre-BLUP period from 1983 to 1991, and 3) post-BLUP period from 1992 to 1999.
Long-Term Genetic Contributions
The effect of different cohorts of ancestors on genetic gain was investigated by studying the relationship between their long-term genetic contributions and index scores. The long-term contribution (r) of an ancestor is defined as the proportion of genes it contributes over the long term to the population (Wray and Thompson, 1990
). Over many generations, in a population thoroughly mixed, the r of an ancestor will converge to the same value for all of its descendants but will differ among ancestors (Woolliams et al., 1999
). Long-term contributions were computed following the approach used by Woolliams and Mäntysaari (1995)
. To compute r, a generation of ancestors and a generation of descendants were defined according to average generation intervals previously calculated. Thus, the ancestral and descendant generations were defined by using L. This definition ensures that r summed over all ancestors over a period of L yr equals unity (Bijma and Woolliams, 1999
). Convergence of contributions was assumed if the variance of contributions of ancestors across descendants was lower than 1.0 x 10-4. For Meatlinc (where L was about 2 yr), contributions were calculated for two generations of ancestors: a) the cohorts born between 1983 and 1984 and b) the cohorts born between 1991 and 1992 (i.e., the first generation after the introduction of BLUP evaluation). For both groups of ancestors, descendants were the cohorts born between 1999 and 2000. For Aberdeen Angus, ancestors were the cohorts born between 1976 and 1979 (L was about 4 yr from 1971 to 1988), and the descendants were the cohorts born between 1995 and 1999 (L was about 5 yr from 1988 to 2000). The regression of the long-term genetic contribution of ancestors on their index scores was calculated for each cohort of ancestors.
Optimizing Genetic Contributions for Maximizing Genetic Gain
The potential extra genetic gains expected from using selection tools based upon the algorithm described by Meuwissen (1997)
were investigated. The algorithm was used to obtain the number of individuals to be selected and the number of offspring each of them should contribute to the next cohort, to achieve the maximum
G while constraining
F to a specific value. Different restrictions on
F were considered. The algorithm maximized the following objective function (Meuwissen, 1997
):
![]() |
where ct is the solution vector of mating proportions (c) of candidates at generation t, gt is the vector of EBV of selection candidates, At is the numerator relationship matrix for selection candidates, Q is a known incidence matrix for the sex of the candidates, 1T equals [1 1], and
0 and
are Lagrangian multipliers. The restriction on the inbreeding rate was achieved each generation by setting Ct = 2[
F +(1 -
F)Ft], where Ft is the average inbreeding coefficient of selection candidates. The third term in the objective function ensures that male and female parents contribute with a half of the gene pool each. Selected candidates are those with c > 0 and will contribute to the next generation according to their c value.
The optimization described above does not take into account any constraint on the maximum contribution a particular candidate may have that may arise from reproductive limitations. This might not be a problem in males since AI techniques are often widespread in livestock populations. However, it can be unrealistic for female candidates for which high reproductive rates are less feasible, particularly in beef cattle and sheep populations. In order to obtain more-realistic results, another set of optimizations was run with an additional constraint on the female contributions. In this case, all females were selected by setting their contributions to a predefined value (i.e.,
nf, where nf is the number of female candidates). This implies that all female candidates are selected and only male mating proportions are optimized. The objective function was modified following Meuwissen (1997, Appendix)
:
![]() |
where c1t is the solution vector of mating proportions of male candidates at generation t; g1t is the vector of EBV of male candidates; A11t and A12t are submatrices of A including only male, and male by female candidates, respectively; c2t is the known vector of female mating proportions; Kt is
; A22t C2t, A22t is the relatedness matrix for female candidates; sT is a vector with constant values [0
]; and Q1 is a known incidence matrix for males analogous to Q in the unconstrained case. Software was developed in Fortran 90 to solve the objective functions described above.
Potential benefits from using optimized contributions were estimated by comparing the expected index gains obtained by using the selection algorithm after mimicking selection in 1999 to 1) the actual observed
G in 2000 and 2) the expected
G in 2000 under truncation selection (i.e., equal contributions) at the observed
F in the population being evaluated. The expected
G from truncation selection was calculated by allocating a fixed mating proportion to female candidates (i.e., equivalent to one mating) and by selecting the number of male candidates that gave the observed
F. This latter comparison allows evaluating the expected benefits from optimizing contributions independently to the benefits of selecting solely on the index. The fact that in practice selection intensity might be lower than that achievable if selection decisions include criteria other than exclusively BLUP-EBV (e.g., Lewis and Simm, 2000
) is not accounted for in the first comparison.
Candidates for the selection algorithms were defined by using both Lm and Lf. Therefore, for Meatlinc, candidates where those males born in 1999 (Lm = 1.0 yr) and those females born from 1996 to 1998 inclusive (Lf = 3.0 yr). The total number of candidates was 1,841. For Aberdeen Angus, candidates were those males and females born from 1992 to 1998 inclusive (Lm = Lf = 5.0 yr) and this gave a total number of candidates of 55,553. However, in order to reduce computing requirements, a preselection of candidates was performed by imposing a minimum index score. For Meatlinc, 1,297 candidates (395 males and 902 females) with index score equal to or greater than 179.0 were included. For Aberdeen Angus, 6,429 candidates (3,321 males and 3,108 females) with index score equal to or greater than 21.0 were included. When only male mating proportions are optimized, computer requirements were higher and, in this case, only 417 male candidates (those with index score equal to or higher than 30.0) were included in the Aberdeen Angus optimization. The index scores were those obtained from the M.L.C. genetic evaluation in 2000, and
F was constrained to a range of values including the observed inbreeding rate per generation in each breed.
| Results |
|---|
|
|
|---|
|
|
|
|
The average number of offspring per male across year in Aberdeen Angus (18.4) was very close to the upper bound of the 25% to 75% interquartile range (2 to 19; Table 1
), indicating a much more skewed distribution than for Meatlinc, where the average (37.6) falls near the mid-point of the range (27 to 43; Table 1
).
Generation Intervals
Figure 2
shows the average generation interval over years for males (Lm) and females (Lf) for Meatlinc and Aberdeen Angus. In Meatlinc, Lm was calculated from 1983 onward because ram dates of birth were not available before that year. An important increase in Lf over years was observed in the period 1976 to 1983. This increase is related to the period of establishment of this synthetic breed in which females had to be kept in the flock for more time. From 1984 onward, Lf remained unchanged and the average was 3.2 yr. In this population, Lm remained unchanged around a value of 1.0 yr over the pre-BLUP period but has slightly increased over the post-BLUP period up to 1.4 yr.
In Aberdeen Angus both Lm and Lf increased at similar rates (0.11 ± 0.02 and 0.16 ± 0.01/yr, respectively) during the period 1976 to 1987 (Figure 2b
). The average Lm and Lf in this period were 4.6 and 4.7 yr, respectively. Over the last 12 analyzed years (1988 to 2000), Lm and Lf averaged 5.2 and 5.7 yr, respectively, although since 1994 the generation intervals started to diverge. By the year 2000, Lf was around 1 yr larger than Lm. There was no evidence to link this increase in Lm and the use of BLUP-EBV.
Rates of Genetic Progress
Figure 3
shows the average index values per year of birth for Meatlinc and for Aberdeen Angus for the periods 1982 to 2000 and 1970 to 1999, respectively. Results indicate that the introduction of BLUP evaluations led to a sustained increase in the rate of genetic gain in Meatlinc from 1994. The difference between
G in the pre- and post-BLUP periods was statistically significant. For this breed,
G was 5.5 ± 1.0 (P < 0.01) index units per year in the pre-BLUP period, and 16.5 ± 0.6 (P < 0.01) index units per year in the post-BLUP period. On the other hand, the
G for Aberdeen Angus before and after the BLUP introduction were not significantly different. The pre-BLUP and post-BLUP rates of gain were 0.55 ± 0.04 (P < 0.01) and 0.46 ± 0.05 (P < 0.01) index units per year, respectively.
Rates of Inbreeding and Long-Term Contributions
The average inbreeding coefficient (F) in the Meatlinc population in 2000 was 6.3% (Figure 3a
). The
F per year for the period 1982 to 2000 was 0.19% (P < 0.001). The difference between
F in the pre- and post-BLUP periods, 0.21% ± 1.31% and 0.23% ± 0.05%, respectively, was not significant. Nevertheless, the pre-BLUP estimation of
F should be taken with caution because F fluctuated considerably in this period. Considering that the generation interval of the population in the post-BLUP was about 2.3 yr (Figure 2a
), the
F per generation in this period was 0.53%. This is equivalent to an effective size of the population (Ne) of 95 animals (i.e., Ne = 
F). On the other hand, the
F increased in the last generation up to about 1.0% (i.e., Ne = 50).
The average F in the Aberdeen Angus population in 1999 was about 0.97% (Figure 3b
). For the period 1974 to 1999,
F was 0.04% per yr (P < 0.001). As with Meatlinc, the rate of inbreeding for the pre- and post-BLUP periods were similar (0.02% ± 0.008% and 0.03% ± 0.008%, respectively). Considering that the generation interval was about 5 yr (Figure 2b
),
F per generation in the post-BLUP period was approximately 0.15% (i.e., Ne = 333). The
F in the last generation (i.e., from 1994 to 1999, inclusive) was about 0.20%.
For Meatlinc, 203 ancestors born between 1983 and 1984 (21 males and 182 females) were identified for computation of their long-term contributions to the 2,094 descendants born between 1999 and 2000. The relationship between long-term genetic contributions of these ancestors and their index values is shown in Figure 4a
. The regression coefficients of contributions on index scores were not significant for all ancestors (4.1 x 10-5; P = 0.15), or for ram ancestors (-1.6 x 10-4; P = 0.13), but were significant for ewes (6.8 x 10-5; P < 0.01). For the analysis of contributions for the first generation after BLUP evaluation, long-term genetic contributions for the 1,337 ancestors born from 1991 to 1992 (643 males and 694 females) to descendants born between 1999 and 2000 were computed. In this case, the regressions of contributions on index scores for this set of ancestors were significant for both males (3.6 x 10-5; P < 0.01) and females (2.4 x 10-5; P < 0.01).
|
Expected Increases in Genetic Merit from Applying Optimized Selection
Table 2
shows the optimization results for both populations. Three levels of
F constraints (0.5%, 1.0%, and 2.0%) were considered for Meatlinc, and four levels (0.2%, 0.5%, 1.0%, and 2.0%) were considered for Aberdeen Angus. The observed
F at the last generation were 1.0% for Meatlinc and 0.2% for Aberdeen Angus.
|
GP-2000) in 2000 in both populations (Table 2
GO-2000) from 1999 to 2000 was 16.1 index units for Meatlinc and 3.6 index units for Aberdeen Angus (see Table 2
F was restricted to the
F observed in the last generation, the
GP-2000 were 70.2 index units (i.e., 4.4-fold over
GO-2000) for Meatlinc (
F = 1.0%) and of 21.1 index units (i.e., around sixfold over
GO-2000) for Aberdeen Angus (
F = 0.2%). Further relaxation of the
F constraint led to higher increments in index gain. For instance, at the most relaxed constraint (
F = 2.0%), the expected increments over the observed
G were 65.6 index units for Meatlinc and 22.1 index units for Aberdeen Angus. However, the relaxation in the restriction on
F contributed only to relatively small increases in
G in comparison to the increases observed by optimizing contributions. It should be noted that for Meatlinc an increment in
G of 46.5 index units was expected even if
F was restricted to a value as low as 0.5%, which was the
F observed over the post-BLUP period.
At the tightest constraint in
F, the number of selected candidates was 80 (i.e., 31 males and 49 females) for Meatlinc and 149 (i.e., 68 males and 81 females) for Aberdeen Angus. As the
F restriction was less severe, the number of selected candidates decreased, as expected. For the most relaxed constraint (
F = 2.0%), it dropped down to 55 (i.e., 18 males and 37 females) for Meatlinc and to 36 (i.e., 17 males and 19 females) for Aberdeen Angus.
Optimization of Male Contributions When All Females Are Selected.
When a more realistic scenario for typical production systems, where female contributions are restricted, the algorithm still achieved significant predicted increases in index score gains (Table 2
). At the observed
F in the last generation,
GP-2000 was 32.2 index units (i.e., twofold over
GO-2000) for Meatlinc (
F = 1.0%) and 15.2 index units (i.e., around fourfold over
GO-2000) for Aberdeen Angus (
F = 0.2%). As when contributions of both sexes were optimized, the relaxation of the constraint allowed for even higher predicted increases. For the most relaxed constraint (
F = 2.0%), predicted increases in index gain were 32.6 index units for Meatlinc and 14.8 index units for Aberdeen Angus.
For Meatlinc, when
F was restricted to a value lower than that observed in the last generation (i.e., 0.5%),
GP-2000 was 8.6 index units lower than the observed index gain in 2000 as it implied a very tight constraint (Table 2
).
Because female contributions were fixed, the expected relative gains over the observed average index scores arose only from the management of the male selection intensity. The number of selected rams in Meatlinc decreased by relaxing the constraint on
F from 58 (
F = 0.5%) to 19 (
F = 2.0%), whereas for Aberdeen Angus the number of selected bulls decreased from 67 (
F = 0.2%) to 11 (
F = 2.0%) (Table 2
).
The predicted benefits over the observed
G might be overestimated because, at the current
F, breeders may be able to achieve higher selection intensities if selection were based solely on index values. The rates of gain under truncation selection (i.e., with equal mating proportions) that gave the observed
F were 27.6 index units for Meatlinc and 11.7 units for Aberdeen Angus (Table 2
). Thus, benefits of optimal selection when only male contributions were optimized over truncation selection (based exclusively on EBV) at the observed
F were 17% for Meatlinc and 30% for Aberdeen Angus. This suggests that the above twofold and fourfold expected benefits in
G for Meatlinc and Aberdeen Angus, respectively, may be overpredictions.
Relationship Between Optimized Mating Proportions and Index Scores
To achieve the restriction on
F, the more severe the
F constraint was, the higher was the number of selected candidates and the more alike were their optimized mating proportions. This behavior across
F constraints can be seen in Figure 5
for Aberdeen Angus where selection was mimicked at 1999. The same pattern of behavior was observed for Meatlinc.
|
F constraint was relaxed, the variance of optimized mating proportions among selected candidates increased from 4.3 x 10-5 (
F = 0.2%) to 4.9 x 10-4 (
F = 2.0%), whereas the variance of the index score among selected candidates decreased from 23.0 (
F = 0.2%) to 13.6 (
F= 2.0%). The highest optimal mating proportion was assigned to the individual with the highest index score (5.20 units) and ranged from 0.042 (
F= 0.2%) to 0.102 (
F = 2.0%). | Discussion |
|---|
|
|
|---|
F was increased by the introduction of BLUP. The application of dynamic selection tools for maximizing genetic gain while constraining rates of inbreeding to target levels would have led to important benefits in
G compared to what has been observed. This demonstrates the scope for this kind of optimization tool in livestock breeding programs.
The Impact of Artificial Insemination
The impact of AI on both populations can be clearly seen in the increase in mating ratios of the breeding males across years. These days, this technique is a standard reproductive technique in beef cattle; however, it is generally less widespread in sheep populations, where the AI procedures are much more complex and success rates are typically much lower. Nevertheless, AI techniques in the United Kingdom have been promoted as one of the key elements for the establishment of sire reference schemes (SRS). The SRS enabled BLUP evaluations across flocks and thereby increased the potential benefit from the use of BLUP (Simm, 1998
). The widespread use of some sires has not only led to an increase in d, but also in the variance of the number of offspring per male (results not shown).
It might be anticipated that increases in the number of offspring per male and its variance would have led to increases in
F. However,
F has remained relatively steady, particularly in the Aberdeen Angus. This has been due to the expansion of the recorded breed numbers, and in particular to an increase in the number of bulls used per year in the population over the period studied. This simple step has reduced the proportional contributions of individual males to the gene pool and so limited the expected increase in
F.
The Effect of BLUP
The introduction of BLUP evaluations seems to have led to an increase in
G for the Meatlinc, but this response was not observed in the Aberdeen Angus. One reason for this difference is the structure of the populations. Whereas the Meatlinc consists of four closely cooperating flocks, with selection policies closely defined by the selection index (H. Fell, personal communication), the Aberdeen Angus has a looser breeding pyramid, with about 200 herds and where policies of individual breeders might not be so closely determined by the society alone. Also, in beef cattle, there is likely to be a much higher use of older "proven" males via AI than in sheep. The Meatlinc may be better placed to utilize the more accurate information arising from the BLUP evaluation, and in combination with AI, the better information across flocks produced by the SRS. A further example of the impact of BLUP in the Meatlinc is the evidence of increases in Lm during the post-BLUP period where ram usage was extended for a longer period because of the better comparison across age groups made possible by BLUP.
The exclusive use of BLUP-EBV as a ranking tool for truncation selection would be expected to lead to an increase in
F (Quinton et al. 1992
), due to an increased coselection of relatives. Nevertheless, this phenomenon was not observed in either breed during the post-BLUP period, which may be particularly surprising in the Meatlinc where closer attention was paid to the index evaluations. However, the increase in breeding males used per year in both populations has proved effective in managing
F over the short term. The changes in Lm observed post-BLUP for Meatlinc would have led to a reduction in
F per year, but this benefit would be offset by the larger lifetime genetic contributions arising from rams kept more than 1 yr.
The estimated effective population sizes in the post-BLUP period are 95 for Meatlinc and 333 for Aberdeen Angus (i.e.,
F of 0.53 and 0.15% per generation, respectively). These values are above the minimum reference value of 40 animals of Goddard and Smith (1990)
for maximizing net genetic response for total economic merit in dairy cattle and fall within the critical range of 30 to 250 animals of Meuwissen and Woolliams (1994)
for balancing decreases in fitness due to inbreeding and increases in fitness due to natural selection. Nevertheless, there was a substantial decrease in effective size in Meatlinc in the last generation to a value comparable to the minimum effective size of 50 recommended by F.A.O. (1998)
. Thus, the application of methods for avoiding further future increases in
F in this population is advisable.
When Meatlinc ancestors born from 1991 to 1992 were analyzed, those with higher index values tended to have larger long-term contributions because regressions on index values were positive and significant compared to a more uniform relationship between ram usage and index scores during the early establishment phase of this synthetic population. This result clearly coincides with the higher genetic gains achieved after the implementation of BLUP-EBV in 1991. A positive association of contributions with index EBV was also observed in Aberdeen Angus for the only set of ancestors analyzed (i.e., born from 1976 to 1979). For this breed, a comparison of the distribution of contributions before and after the introduction of BLUP is difficult because at most two generations have passed since the introduction. Long-term genetic contributions require five or more generations to achieve a reasonable degree of stability. Higher regressions were observed for male ancestors in both populations, indicating higher selection intensities applied on male than on females in accordance with expectations (Woolliams et al., 1999
).
The Effect of Optimized Selection
Although the Meatlinc has increased its rate of gain using BLUP evaluations and simultaneously managed its
F by increasing the number of males selected per year (so reducing the selection intensity), further gains are possible by using the selection algorithms. The results showed that these algorithms would benefit both Aberdeen Angus and Meatlinc over a range of values of
F. The most dramatic increases in
G were obtained when selection was allowed in both males and females. However, these gains assume unrealistic reproductive rates for females. Substantial and valuable increases in
G were obtained when no selection among females was allowed. At the
F in the last generation, the benefits over the observed
G were 16.1 index units (i.e., twofold) for Meatlinc and 11.6 index units (i.e., threefold) for Aberdeen Angus.
The comparisons of the expected (from optimized mating proportions) with observed index scores, implicitly assumed that in practice selection in the two populations has been exclusively based on index values. In practical breeding schemes, however, selection decisions are based not only on EBV, but also on other factors (e.g., physical and reproductive soundness). This reduces the selection intensity and the maximum genetic merit achievable. For instance, Lewis and Simm (2000)
found that losses in selection intensity in sheep SRS lead to genetic response 0.58 to 0.69 times that obtained when strictly the best animals were selected on BLUP-EBV. Results in Meatlinc support this expectation, where the ratio
GO-2000/
GT-2000 was 0.58 (see Table 2
). This effect was more important in Aberdeen Angus, where the ratio
GO-2000/
GT-2000 was only 0.31. In this case, because of the scale of population, selection decisions might be restricted within herds or groups of breeders. In contrast, the Meatlinc SRS is managed as a single selection unit with a single selection policy and tight cooperation among flocks (H. Fell, personal communication). The predicted benefits of 17% for Meatlinc and 30% for Aberdeen Angus over the expected gain under truncation selection provides an even more realistic evaluation of the benefits of optimal selection when compared to traditional truncation selection based solely on BLUP-EBV. In addition, deterministic predictions of the rate of genetic gain at predefined rates of inbreeding (our unpublished observations) showed expected increases ranging from 20 to 40% for
F = 1.0%.
The
G from optimal selection after optimizing only male contributions of 15.2 index units (see Table 2
) at the observed
F in Aberdeen Angus was similar to the expected
G (i.e., 15.7 index units, result not shown) after optimizing the contributions of selected bulls in 1999 with observed offspring in 2000 conditional to the observed dam contribution. This indicates that the expected extra index gain from optimal selection are realistic and would be achievable by only optimizing the usage of the current set of bulls selected by the breeders using the available selection index.
In practice, about 30 rams and 700 ewes are used each year in Meatlinc, but when contributions of both sexes were optimized, the optimization algorithm implied selection of 25 males and 43 females on average (Table 2
). Similarly, in Aberdeen Angus the actual numbers of breeding animals (about 1,000 males and 10,000 females) are much larger than those obtained after applying the optimization tool on both sexes (43 males and 49 females on average, Table 2
). As pointed out before, these optimum numbers of selected candidates imply very high reproductive rates, particularly on the female side. On the other hand, much more realistic outcomes from the application of the optimization tool were obtained when only male contributions were optimized conditional to fixed female contributions. At the observed
F, the maximum contribution allocated to a male was about 0.046 and 0.048 for Meatlinc and Aberdeen Angus, respectively. This is equivalent to an expected maximum optimum number of matings per male of 83 (i.e., 2 x 0.046 x 902) for Meatlinc and 298 (i.e., 2 x 0.048 x 3,108) for Aberdeen Angus. These optimum numbers and the differential usage arising from optimizing contributions can be readily achieved through AI, which is currently a standard male reproductive technique. Moreover, they imply a maximum number of offspring per male that are well below the observed upper limit for the number of offspring per male range in each breed (see Table 1
).
The change in the slope and intercept of the regression coefficient between mating proportions and index scores with different restrictions on
F observed here agrees with the general form of the optimal solutions stated by Grundy et al. (1998)
. Basically, as a less severe constraint was imposed, fewer individuals were selected, the usage of the individuals became more unequal, and the slope of the regression was higher (Figure 5
). Although in the initial cohorts after applying the dynamic selection tools the selective advantage may be the index EBV as suggested by Figure 5
, once the use of the dynamic selection algorithm is established, the selection advantage is given by the estimated Mendelian sampling term of the index (Woolliams et al., 2002
).
Weigel and Lin (2002)
applied the algorithm of Meuwissen (1997)
in five major U.S. dairy breeds, concluding that genetic gain may be sacrificed by imposing constraints on inbreeding. Their conclusion could be somewhat misleading because they did not compare the predicted average genetic merit from the optimization to that obtained without optimization at the observed level of inbreeding. Our results clearly indicate that no losses of gain are expected when the inbreeding rate is constrained to the observed value, but additional gains are expected. The only situation where a lower relative gain with respect to the observed gain was obtained, occurred in the extreme case when a constraint equivalent to a lower than the observed
F in Meatlinc was applied, after fixing all female contributions.
The practical realization of the optimal contributions and the average expected index scores requires a coordinated policy of the use of selected candidates among the different breeding flocks or herds making up the breeding population. This could be a reasonable target in small-scale schemes with coordinated breeding policies, but clearly would be much more difficult in large schemes in which different objectives might be pursued. Hence, in breeds with large-scale breeding programs, a more reasonable approach would be to apply the optimization tool on individual herds or groups of herds with coordinated selection policies and objectives.
Two methodological aspects of this optimization should be addressed. Firstly, the use of EBV obtained in 2000 to mimic selection in 1999, instead of using the EBV obtained in 1999, is not expected to have affected the results obtained because no significant changes in candidate ranking were found when EBV from both evaluations were compared. Secondly, the preselection of candidates would not have affected the optimization outcome. Among the group of higher merit candidates, those in the bottom half were never selected, indicating that any discarded candidate would not have made a significant contribution if it had been included in the optimization. On the other hand, this allowed a significant reduction of computational requirements.
More flexible constraints may be required in breeding programs with particular features or breeding structure. These may include setting a maximum contribution per male (i.e., a minimum number of sires), a fixed contribution for a particular set of males, or a desired contribution of a group of females (e.g., in a nucleus). These constraints can be accommodated with the same tool used here (after Appendix in Meuwissen, 1997
). For instance, for Meatlinc, an additional optimization with a maximum number of 20 female mates per selected ram (i.e., a minimum of 45 selected rams) was set (not shown). Accordingly, the selection tool found a feasible solution by selecting 46 males, of which 45 were allocated a maximum fixed contribution of 0.01. Moreover, even with this highly restrictive constraint, the expected
G at the current
F was still about 6.0 index units higher than the observed
G.
Although the optimization approach used here can be realistically applied in practical livestock breeding programs, evolutionary computation strategies may provide a more flexible framework for setting a greater variety of constraints. Genetic algorithms have been used as optimization tools in livestock breeding programs (e.g., Shepherd and Kinghorn, 1998
, Meszaros et al., 1999
), and they could be extended to explicitly restrict
F (e.g., Correnti, 2002
).
After allocating optimal mating proportions to the selected candidates, the following step in a breeding program is to decide a mating policy. Sonesson and Meuwissen (2000)
found that the optimization tool used here combined with mating systems that restricted either mating pairs coancestry or offspring coancestry, achieved 22% higher response than random mating, in particular for stringent constraints on
F. Nevertheless, the extra benefit from the use of nonrandom matings was reduced as the size of the scheme increased.
Our results refer to practical livestock populations in which the main objective is to achieve the highest genetic gain for a given
F. However, the approach is also valid for conservation purposes in which the aim may be to minimize
F while achieving a predefined level of genetic gain (Villanueva et al., 2003
). Hence, genetic improvement and conservation can be taken as the extremes of a broader optimization problem with particular relative emphasis given to the gain and inbreeding.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 Correspondence: Sustainable Livestock Systems Group, Livestock Breeding, Bush Estate, Penicuik, EH26 0PH, Scotland, U.K. (phone: +44-131-535-3241; fax: +44-131-535-3121; E-mail: s.avendano{at}ed.sac.ac.uk).
Received for publication December 23, 2002. Accepted for publication August 25, 2003.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Hinrichs, M. Wetten, and T. H. E. Meuwissen An algorithm to compute optimal genetic contributions in selection programs with large numbers of candidates J Anim Sci, December 1, 2006; 84(12): 3212 - 3218. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. F. Kearney, E. Wall, B. Villanueva, and M. P. Coffey Inbreeding Trends and Application of Optimized Selection in the UK Holstein Population J Dairy Sci, October 1, 2004; 87(10): 3503 - 3509. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME |