Selective advantage of implementing optimal contributions selection and timescales for the convergence of long-term genetic contributions.

BACKGROUND: Optimal contributions selection (OCS) provides animal breeders with a framework for maximising genetic gain for a predefined rate of inbreeding. Simulation studies have indicated that the source of the selective advantage of OCS is derived from breeding decisions being more closely aligned with estimates of Mendelian sampling terms ([Formula: see text]) of selection candidates, rather than estimated breeding values (EBV). This study represents the first attempt to assess the source of the selective advantage provided by OCS using a commercial pig population and by testing three hypotheses: (1) OCS places more emphasis on [Formula: see text] compared to EBV for determining which animals were selected as parents, (2) OCS places more emphasis on [Formula: see text] compared to EBV for determining which of those parents were selected to make a long-term genetic contribution (r), and (3) OCS places more emphasis on [Formula: see text] compared to EBV for determining the magnitude of r. The population studied also provided an opportunity to investigate the convergence of r over time.
RESULTS: Selection intensity limited the number of males available for analysis, but females provided some evidence that the selective advantage derived from applying an OCS algorithm resulted from greater weighting being placed on [Formula: see text] during the process of decision-making. Male r were found to converge initially at a faster rate than female r, with approximately 90% convergence achieved within seven generations across both sexes.
CONCLUSIONS: This study of commercial data provides some support to results from theoretical and simulation studies that the source of selective advantage from OCS comes from [Formula: see text]. The implication that genomic selection (GS) improves estimation of [Formula: see text] should allow for even greater genetic gains for a predefined rate of inbreeding, once the synergistic benefits of combining OCS and GS are realised.

Background

Selection theory [1, 2] states that sustained genetic gain () is obtained from creating a covariance between the Mendelian sampling terms () and long-term genetic contributions () of selection candidates. This in turn provides the framework for an effective solution for the management of genomic diversity, due to the relationship between and rate of inbreeding (). Various methods [3-9] have been proposed for managing genetic resources using this framework, which has become known as optimal contributions selection (OCS). OCS has been shown to generate greater for a predefined , and to generate lower for a given , compared to earlier methods that used truncation selection [10, 11].

So far, the source of the selective advantage provided by OCS has only been investigated in theoretical and simulation studies [7, 12, 13], with the conclusion that it arises from placing greater selection pressure on estimates of Mendelian sampling terms of selection candidates () than in truncation selection. This study is the first to test these predictions in a real population by assessing the weightings placed on estimated breeding value () versus , both before and during OCS. Given the findings of Avendaño et al. [12] and Grundy et al. [7], the following three hypotheses were tested: (1) OCS places more emphasis on compared to for determining the initial selection of parents; (2) OCS places more emphasis on compared to for determining which parents make a non-zero ; and (3) OCS places more emphasis on compared to for determining the magnitude of . Timescales for the convergence of of selected parents were also examined.

Methods

Population analysed

The commercial population analysed in this study consisted of 115,428 pigs comprising a single closed nucleus line that was bred continuously between 12th April 1997 and 16th August 2012, with pedigree information available for all individuals. The line was developed and also expanded over this period, leading up to the cohort studied here. that were calculated using a weighted index of traits were used at the time of selection, but the used for assessment in this study were only available at a single time point (16th August 2012). Individuals without a calculated on themselves or on their parents were excluded to ensure that the estimated Mendelian sampling term, , could be calculated for each individual. This left a total of 107,895 individuals (54,881 males and 53,014 females) for analysis, with the earliest individual being born on 23rd February 1999, which is denoted as time .

Woolliams et al. [1] defined the long-term genetic contribution, , of an ancestor to an individual as the proportion of the genes of that are expected to be derived by descent from ancestor . Using this definition of , each individual had its long-term contribution () calculated as follows:where the sum was over all individuals born in the 2012 cohort and was the total number of animals in the cohort. Generation intervals were based on the definition of Bijma and Woolliams [2], such that the total number of generations at time was defined as , where the sum is over all individuals born from time 0 up to and including time . Individuals of generation were defined as all individuals born at time for which . Using this definition, 10.4 generations were born between 23rd February 1999 and 16th August 2012, representing a generation interval of L = 1.29 years, with the base generation () consisting of individuals born between 23rd February 1999 and 21st May 2000.

Time scales for convergence of contributions

The value of for individual in a closed population will converge towards an asymptote over time [1]. To ensure the validity of the results from this study, was required to approach convergence. To test this convergence, contributions were examined as a function of time after birth for all individuals selected as parents and born during 1999. The final contributions for this 1999 cohort were those obtained from Eq. (1), i.e. their average contribution to the 2012 cohort, i.e. 13 years later. Contributions of the 1999 cohort to the cohorts born in each of the intermediate years, 2000–2011, were also calculated, denoted as , where denotes time after birth, with u = 1, 2, …, 12 years. For each intermediate year, the following linear regression model was used to determine the degree of convergence in that year:where is an error term and both the response and explanatory variables are adjusted by , which is the mean of the final contributions for all the individuals selected as parents and born during 1999. At convergence, and . This model was applied separately to contributions of males and females.

A similar methodology to that applied annually was used to assess the convergence of on a generational scale using all individuals selected as parents in the base generation , i.e. . Contributions of individuals in the base generation to each intermediate generation were calculated for generations 1 through 8 as , where , and final contributions were calculated as those made to , i.e. the final complete generation. This final generation consisted of individuals born between 6th December 2010 and 5th May 2012.

Grouping of individuals according to the selection method

An OCS algorithm was introduced for selection in this population in 2002 and was based on the approach described by Newman et al. [14] and Kinghorn [15]. An iterative algorithm was applied to calculate a response surface, with the aim of achieving maximal for a given . Prior to the introduction of this algorithm, selection had been by truncation with ad hoc restrictions aimed at controlling inbreeding. The full transition to this method of selection was completed by the end of 2005 and was used throughout the remainder of the period studied.

Individuals were assigned to either the Pre-OCS group, the OCS group, or neither of these groups using defined above. Comparisons between and within the Pre-OCS and OCS groups were conducted to test the three hypotheses described below. A total of 17,165 individuals born between 23rd February 1999 and 7th December 2001 were assigned to the Pre-OCS group, which represented two generations with . Since the OCS algorithm was introduced in 2002, only two generations were available prior to the change to OCS. A further 25,068 individuals, born between 23rd February 2006 and 7th October 2008, were assigned to the OCS group, also consisting of two generations, . Analysing two generations in the OCS period maintained consistency with the Pre-OCS period in the analyses. In addition to , the selection score () was determined for all individuals in the Pre-OCS and OCS groups from the recorded birth of offspring (, if the individual had offspring, and otherwise); it is clear that implies . Table 1 shows the numbers of individuals in the Pre-OCS and OCS groups by sex and contribution history.

Table 1

Number of males and females depending on selection group (Pre-OCS and OCS), selection score (), and whether or not the long-term contribution in 2012 () was positive

Constraint		Pre-OCS		OCS
x_i	r_i	Males	Females	Males	Females
–	–	8341	8824	13,001	12,067
> 0	–	177 (0.021)	1444 (0.164)	111 (0.009)	1839 (0.152)
> 0	> 0	118 (0.667)	219 (0.152)	35 (0.315)	179 (0.097)

The symbol ‘–’ indicates that no constraint was applied. The proportion of individuals remaining from the constraint in the previous row is in brackets

Estimates of Mendelian sampling terms for each individual , , were calculated as the deviation of the individual’s from the average of its parents. In the analyses that follow, the were calculated as deviations from the mean EBV of contemporaries, defined by all animals born within 30 days before or after birth of that individual. This on average resulted in 1497 contemporaries. Unless stated otherwise, the remaining references to refer to this adjusted . No adjustment for was required since it has an expectation of 0 for all animals.

The average reliability of the across both sexes () was approximately 0.16 in the pre-OCS period, increasing to 0.25 thereafter. Selection took place at off-test (~ 6 months of age), so both sexes had about the same volume of information available for estimating their initial breeding values. In 2012, the average across both sexes was 0.28.

Inbreeding coefficients () for individuals were obtained by subtracting 1 from the diagonal elements of the numerator relationship matrix [16]. Annual rates of inbreeding were calculated for the birth cohorts at year using two methods: (i) , calculated as:where is the average for all individuals in the cohort;

and (ii) , calculated as:where the sum is over all individuals in the cohort and is defined by Eq. (1), regardless of convergence.

Sources of selective advantages with OCS

Bivariate regressions were used to determine whether or provided the source of the selective advantages before and during OCS. The source of the selective advantage was examined for: (1) which individuals were initially selected as parents, (2) which parents went on to make a non-zero to the population, and for (3) determined the magnitude of .

Step 1: Initial selection as parents

Bivariate regression was used to determine whether or was the principal variable that determined the initial selection of an individual as a parent, separately for the Pre-OCS and OCS ancestors. A generalised linear model with a binomial distribution and a logit link function was fitted together with a parameter for over-dispersion, assuming a beta-binomial distribution for conditional on the fixed effects. Let the probability of selection for individual be , then if is the logistic link function, the following model was fitted:where is an intercept and and are the respective regression coefficients on and . Changes in deviance conditional on all other model terms were used to assess the relative importance of individual model terms. This analysis was conducted separately for each sex using the ‘quasi-binomial’ option in the ‘glm’ package of R, which estimates over-dispersion from the Pearson residuals.

Step 2: Maintenance of contributions over time for selected individuals

To examine the maintenance of contributions over time, models were fitted only to the subset of individuals that had been selected as parents, i.e. conditional on , because long-term contributions are always 0 for individuals that are not selected. For this subset, if , and 0 otherwise. The methods described for Step 1 were repeated with instead of , such that in Eq. (5).

Step 3: Magnitude of long-term contributions over time for selected individuals

An intrinsic property of OCS is the targeting of a desired magnitude of for each individual. The source of the selective advantage for determining the magnitude of was assessed using: (1) all individuals that were initially selected as parents, i.e. , regardless of their long-term contributions, and (2) restricted to parents with a non-zero long-term contribution (). The strength of the relationships of the magnitude of with both and were investigated separately for each sex in the Pre-OCS group and in the OCS group using the following linear regression model:Finally for the selected subgroup with , contour plots of against and were used to visualise and aid in interpreting the relationship of and with the magnitude of . Contour plots were produced using the ‘akima’ package [17] of R v2.15.

Results

Accumulation and convergence of contributions

A plot of the accumulation of over time is in Additional file 1: Figure S1, starting with the first individual born in the Pre-OCS group until the last selected individual, which was born on 23rd September 2011. The rate of accumulation was relatively constant over time, with a levelling off of the cumulative midway through 2011, which reflects individuals that were currently undergoing selection. The regression coefficient of accumulated on time was 2.14 × 10−3 per day, equating to an average generation interval of 468 days.

Figure 1a illustrates the male and female adjusted R-squared for the model of Eq. (2) for each year. The adjusted R-squared increased rapidly at first for both males and females, with males exceeding 75% after 3 years. For males, continued to converge but at a slower rate, reaching an adjusted R-squared of approximately 95% by 9 years. Females were slower to approach convergence, with the adjusted R-squared only exceeding 75% after 5 years, but still exceeding 95% within 9 years. For selected individuals born in 1999, the standard deviation of within the 2012 cohort was 1.05 × 10−3.

Fig. 1

a Adjusted R-squared of linear regression of the final assumed long-term genetic contribution (in 2012) of all selected males and females born in 1999 () on their contributions to individuals born in each year, from 2000 to 2011. b Adjusted R-squared of linear regression of the final assumed long-term genetic contribution (in generation 9) of all selected males and females born in generation 0 on their contributions to individuals born in each generation, from 1 to 8

The adjusted R-squared for the model of Eq. (2) on a generational basis is in Fig. 1b. This trend line had greater linearity than the yearly analysis. For this analysis, males again initially converged at a faster rate than females, with both sexes achieving an average adjusted R-squared of approximately 75% in five generations. For individuals in the base generation, the standard deviation of in generation 9 was 1.09 × 10−3.

Rates of inbreeding over time

Plots of and over time are shown in Fig. 2. The elevated value for observed in 2006 was strongly influenced by a prolific male that sired 1419 offspring, which alone contributed 0.012 to , representing 63.5% of the total contributions to for the year. The impact of this prolific individual would not be expected on a plot of , since its impact is spread across future generations beyond the endpoint in Fig. 2. is prospective in the sense that it measures the impact further down the pedigree arising from an individual’s gene flow [18, 19], whereas is retrospective in the sense that it measures the accumulated impact of changes in contributions across all previous generations.

Fig. 2

Annual rate of inbreeding based on pedigree and long-term genetic contributions )

Sources of selective advantages with optimal contributions selection

Estimates of regression coefficients from the logistic regression of the selection score on and are in Table 2. There was no evidence of extra-binomial variation among Pre-OCS females (P ≥ 0.05) but there was significant evidence in the remaining subgroups (P < 0.001): the estimated parameters for over dispersion were 1.02, 1.08, 3.7 and 2.3 for Pre-OCS and OCS females, and Pre-OCS and OCS males, respectively. In the Pre-OCS group, was negatively associated with selection in males (P < 0.001), but was positively associated with selection in females (P < 0.001). In both male and female OCS groups, both and were significant factors (P < 0.001) in promoting selection as a parent, with more positive values of either favouring selection, but with having greater regression coefficient estimates. The relative importance of these two terms differed slightly between male and female candidates. Estimates from univariate logistic regression are in Additional file 2: Table S1.

Table 2

Estimates of regression coefficients () from the bivariate logistic regression of selection score on estimated breeding values () and estimated Mendelian sampling terms ()

	Pre-OCS		OCS
	β (s.e.)	F	β (s.e.)	F
Males
EBV	− 0.12 (0.11)	1.09	0.34 (0.10)	11.88
a^	− 1.24 (0.33)	13.61	0.68 (0.22)	9.43
Females
EBV	0.01 (0.02)	0.16	0.16 (0.02)	82.24
a^	0.31 (0.08)	15.11	0.44 (0.04)	105.26

Standard errors (s.e.) are between parentheses. Approximate F-values are shown with numerator d.f. = 1 and denominator d.f. = 8338, 8821, 12,998, and 12,064 for Pre-OCS males and females, and OCS males and females, respectively

Step 2: Maintenance of contributions over time for selected individuals

Estimates of regression coefficients from the logistic regression of on and on in the OCS period are in Table 3. In the Pre-OCS group, there was no significant effect of for either sex but a positive association with was observed for the females (P < 0.05). In the OCS males, was positively associated with maintaining (P < 0.001), whereas was negatively associated with maintaining a positive contribution across generations (P < 0.05), conditional on initial selection as a parent. Among OCS females, both (P < 0.001) and (P < 0.05) were positively associated with maintaining a positive contribution. Estimates of univariate regression coefficients are in Additional file 2: Table S2 and show that, in these models, and were always positively correlated with maintaining non-zero contributions. There was no evidence of extra-binomial variation (P ≥ 0.05) for any of the four sub-groups.

Table 3

Estimates of regression coefficients () from the bivariate logistic regression of maintenance of non-zero contributions on estimated breeding values () and estimated Mendelian sampling terms ()

	Pre-OCS		OCS
	β (s.e.)	F	β (s.e.)	F
Males
EBV	0.03 (0.12)	0.07	0.96 (0.23)	17.81
a^	0.10 (0.22)	0.22	− 0.65 (0.26)	6.45
Females
EBV	− 0.01 (0.06)	0.04	0.25 (0.06)	21.31
a^	0.24 (0.11)	4.37	0.22 (0.10)	5.07

Standard errors (s.e.) are between parentheses. Approximate -values are also shown: with numerator d.f. = 1 and denominator d.f. = 174, 1441, 108 and 1836 for Pre-OCS males and females, and OCS males and females, respectively

Estimates of regression coefficients from the bivariate regression of on and conditional only on are in Table 4 and those conditional on are in Table 5. The additional restriction imposed by substantially increased the standard errors of the estimates and, as a result, estimates in Table 5 do not differ qualitatively from those in Table 4 if these uncertainties are taken into consideration. In males (conditional on ), had a significant effect on during both Pre-OCS (P < 0.05) and OCS (P < 0.01). In females (conditional on ), neither nor had a significant effect on Pre-OCS, but in the female OCS group, was a significant determinant of (P < 0.001). Estimates of univariate regression coefficients, conditional only on and conditional on are in Additional file 2: Tables S3 and S4, respectively.

Table 4

Estimates of regression coefficients () from the bivariate regression of the long-term contributions () on estimated breeding values () and estimated Mendelian sampling terms () for all selected individuals

	Pre-OCS		OCS
	β (s.e.)	F	β (s.e.)	F
Males
EBV	10.98 (5.24)	4.38	45.34 (14.37)	9.95
a^	− 6.97 (9.13)	0.58	9.78 (21.45)	0.21
Females
EBV	0.65 (0.60)	1.18	0.87 (0.51)	2.91
a^	0.77 (1.25)	0.39	3.21 (0.93)	11.92

Standard errors (s.e.) are between parentheses. The importance of the terms is assessed using F-values; numerator d.f. = 1, and denominator d.f. are 174, 1441, 108 and 1836 for Pre-OCS males and females, and OCS males and females, respectively

Table 5

Estimates of regression coefficients () from the bivariate regression of the long-term contributions () on estimated breeding values () and estimated Mendelian sampling terms () for all individuals with

	Pre-OCS		OCS
	β (s.e.)	F	β (s.e.)	F
Males
EBV	15.87 (7.19)	4.87	38.19 (53.82)	0.50
a^	− 13.92 (12.36)	1.27	53.63 (62.43)	0.74
Females
EBV	4.36 (3.07)	2.02	− 2.97 (5.00)	0.35
a^	− 4.08 (5.31)	0.59	15.97 (7.46)	4.58

Standard errors (s.e.) are between parentheses. The importance of the terms is assessed using F-values; the numerator d.f. = 1, and denominator d.f. are 115, 216, 32 and 176 for Pre-OCS males and females, and OCS males and females, respectively

A contour map illustrating the relationship of with and with is in Fig. 3 for selected males with and in Fig. 4 for selected females with . More active control of was generally observed in the OCS selection period compared to the Pre-OCS period. There was a shift in the coloured area of the plot both horizontally to the right and vertically upwards following the introduction of OCS, which shows that animals with had higher and in the OCS period than animals with in the Pre-OCS period. Examining the colour dispersal of the Pre-OCS contour plot revealed that , along the x-axis, was the more important determinant for the magnitude of in the Pre-OCS period. In the OCS period, there was a change in the distribution of , with warmer colours found predominantly in the upper-most half of the contour plot. This demonstrates that after the introduction of OCS, the , on the y-axis, was more important for determining the magnitude of . For males, the largest concentration of warmer colour in the OCS plot was in the upper right hand corner of the plot (Fig. 3), highlighting that both and were important for males with the very highest .

Fig. 3

Contour plot of the magnitude of long-term genetic contributions (), represented by warmth of colour, with regards to estimated Mendelian sampling terms () and estimated breeding values () for males, conditional on .  is plotted along the x-axis and is plotted along the y-axis with colour gradients used to indicate the magnitude of . Due to differences in the maximum value of between the sexes, a different scaling of was used between the male and female plots

Fig. 4

Contour plot of the magnitude of long-term genetic contributions (), represented by warmth of colour, with regards to estimated Mendelian sampling terms () and estimated breeding values () for females, conditional on . is plotted along the x-axis and is plotted along the y-axis with colour gradients used to indicate the magnitude of . Due to differences in the maximum value of between the sexes, a different scaling of was used between the male and female plots

Discussion

The research reported here was made possible by the availability of data from a single, large (n = 107,895) nucleus breeding line from Genus-PIC. This was a novel opportunity to use hypothesis-driven analyses to assess the source of the selective advantage under different selection strategies using longitudinal data. The particular hypotheses investigated were that adoption of OCS, to maximise gain conditional on the target rate of inbreeding, shifts the emphasis of selective advantages from towards estimates of Mendelian sampling terms, , based on both theory and previously published simulations [10]. The outcome of the analyses provides some qualified support for these hypotheses in females.

To maximise accuracy and minimise bias in an assessment of OCS, a number of conditions are required: (i) adequate periods of convergence for each method of selection; (ii) that are relevant to decision-making; (iii) stability in breeding objectives, selection pressures and rates of inbreeding; and, (iv) adequate population size over time. Not all these conditions were met in this commercial operation. In relation to (i), analysis of convergence of contributions for selected individuals (Fig. 1) suggested that eight generations were required to reach convergence, rather than the five generations found in simulations [1], and this was achieved only for the Pre-OCS group. However, the process of convergence within the OCS group was well advanced after five generations (Fig. 2) and this is supported by the evidence from the Pre-OCS group where the correlation between contributions at 5 and 8 generations after birth was 0.86. In relation to (ii), an individual’s and are dynamic and increase in accuracy as more information becomes available from descendants. Selection decisions were made using and that were available at the time of selection, but for the current study these values were only available at a single time point (16th August 2012). Therefore, the accuracies of and used in the analyses varied among individuals. Differences in accuracy are greatest between animals selected as parents, i.e. those that have offspring, and those that were not selected. Therefore, analyses conditional on initial selection as a parent (steps 2 and 3) are likely to be less affected by variation in accuracies of and than analyses for step 1. In relation to (iii), it is expected that changes in breeding objectives occur over time within a commercial population, and the approximate correlation between the selection index applied in 2000 and that applied in 2005 (calculated from the (co)variance matrix of the in the indices and the two sets of weights applied) was only moderate at 0.39. It is clear from Table 1 that selection intensity increased after initial establishment and expansion of the line, which affects quantitative comparisons of Pre-OCS and OCS groups. In relation to (iv), empirically, this population of more than 100,000 individuals born over the course of 15 years was sufficient to detect sources of selective advantages. The standard errors of the partial regression coefficients for and for selected males were relatively large compared to those for females, as a result of the higher selection intensity applied to males. This higher selection intensity in males increased the magnitude of the potential effects [1], but also greatly reduced the sample size.

One advantage of the study is that it was conducted on data collected on a real population rather than simulation, but this also means that day-to-day decision-making is inevitably driven by commercial requirements rather than strictly following selection recommendations. In summary, although the Pre-OCS group acted as a practical control for the OCS-group, the commercial setting of the data indicates that comparisons between the Pre-OCS and OCS groups should be considered qualitatively and with care. However, the comparison of outcomes from the OCS group with the theoretical expectations remains valid.

The observed linear increase in the accumulation of demonstrated a relatively uniform generation interval over time of approximately 468 days. Alternative generation intervals can be calculated for this population by either averaging the age of parents at birth of their first offspring (368 days), or after the birth of the first offspring that is selected (415 days). The latter is closer in concept to the generation interval in quantitative genetics since it includes the concept of replacement, and moves towards that provided by , which is the length of time to renew the gene pool that is destined to maintain the population in the long term.

There was evidence that the generation interval calculated from parent ages increased over time (results not shown) and this was likely associated with operational decisions concerned with expansion of the line (see Table 1). The stochastic simulation work of Bijma and Woolliams [2] and Woolliams et al. [1] found that the generation interval calculated using was shorter than that using the average age of parents of a cohort under simple mass selection with a pre-determined age structure. The primary reason for this qualitative difference between these two methods for calculating generation intervals is that there is a selective advantage to offspring born to younger parents in the mass selection scenario, since their parents have higher breeding values than the population average. With mixed-model evaluations, e.g. best linear unbiased prediction (BLUP), the genetic trend across cohorts is accounted for, which results in more effective selection among candidates of different ages.

Sources of selective advantages of optimal contributions selection

The OCS algorithm guides the optimum over generations until convergence, through consideration of both an individual’s as well as its degree of relatedness with the rest of the population. This approach is expected to reduce the of individuals from over-represented families, allowing the selection of superior animals from families with lower . As a result, individuals with greater are selected to provide comparatively larger [7, 20]. This process results in the maximum gain () conditional on the target rate of inbreeding (), since (i.e. not ) and (see Woolliams et al. [1]). In the case of perfect accuracy of (e.g. with a heritability ) and no interdependence of generations, this would take the form of a linear allocation above a cut-off, whereby selected individuals would be those with above a threshold and is expected to be linearly related to , conditional on selection [12]. Where these conditions do not hold, the relationship between and will have additional random noise and the true breeding value or will also influence . To examine these different aspects, the selection process was broken down into three consecutive steps, each more closely focused on the ultimate genetic contributors to the population, and each with its associated hypothesis.

The first hypothesis was concerned with whether or was the principal factor that determined the selection of an individual to become a parent. In practice, the selection of an individual to become a parent, unlike the development of , is only influenced by selection decisions within a narrow time-period. For females, was observed to be the most important determinant for selection of an individual to become a parent, although it was unclear if this had emerged as a result of the introduction of OCS, in spite of the influence being stronger in the OCS group.

The most convincing evidence for the importance of in determining when using OCS selection came from analyses that addressed the third hypothesis, which related the magnitude of to and , conditional on being selected initially. Of the analyses conducted in this study, these analyses most closely followed a previous simulation study [10] and directly examine the theoretical expectation that OCS will guide contributions in direct relation to among animals that are initially selected. The footprint of OCS that was expected from the simulation study was observed among females of the OCS group, with positively associated with and dominating the influence of . However, the reverse was observed for males, for which inferences are made more difficult due to the large standard errors of the bivariate coefficients. It is possible that short-term commercial needs that do not follow the outcome of the OCS algorithm are more likely to occur among males, compared to females, as male selection opportunities are greater and offer a more rapid response to the need.

Conclusions

This work represents the first evaluation in commercial practice of the impact and validation of the theory underlying OCS, which predicts a closer alignment of with than with . There was some evidence of a re-weighting in the emphasis of selection away from the towards in the final magnitude of , with females providing stronger evidence for these conclusions. The evidence from males was equivocal, in part because of higher selection intensity among males, leading to smaller numbers of males with which to conduct the regression analyses. As a result, this study provides some support to the assertion that there is an advantage in combining OCS with genomic information, which can achieve more accurate prediction of at an earlier age, to generate greater genetic progress over time, without negatively impacting the rate of inbreeding.

Additional file 1: Figure S1. Accumulation of long-term genetic contributions over time with a line of best fit.

Additional file 2: Tables S1, S2, S3, S4. Estimates of regression coefficients () from univariate regressions of selection score (x) on estimated breeding values () and estimated Mendelian sampling terms (). Standard errors (s.e.) are given in parentheses. Table S2. Estimates of regression coefficients () from univariate regressions of (conditional on individual was selected) on estimated breeding values () and estimated Mendelian sampling terms (). Standard errors (s.e.) are given in parentheses. Table S3. Estimates of regression coefficients () from simple linear regression of long-term genetic contribution (with all selected individuals, i.e. conditional on  = 1) on estimated breeding values () and estimated Mendelian sampling terms (). Standard errors (s.e.) are given in parentheses. Table S4. Estimates of regression coefficients () from simple linear regression of long-term genetic contribution (restricted to those individuals with a non-zero long-term genetic contribution) on estimated breeding values () and estimated Mendelian sampling terms (). Standard errors (s.e.) are given in parentheses.