As feed costs make up around 70% of the running cost in a feedlot (Lima et al., 2017; Zhang et al., 2017), feedlotters require accurate predictions of expected feed intake to aid in management and planning (Illius et al., 2000; Pulina et al., 2013). While there are models available for the prediction of voluntary feed intake (VFI) in sheep (Emmans, 1997; Finlayson et al., 1995; National Research Council, 2007; Pulina et al., 2013; Vieira et al., 2013) these models are generally complex and lack breed specificity. This is a significant challenge in South Africa where a wide array of breeds (SA Stud Book, 2015; van der Merwe et al., 2020) can be found, and is further complicated by the use of crossbreeding in certain flocks.
This study therefore aimed to construct simple models for the prediction of feed intake between weaning and maturity in pure- and crossbred South African sheep.
Ethical clearance for this study was granted by the Stellenbosch University REC:ACU (ACU-2020‐14574). The project is briefly summarised below, for the full trial methodology please see Theron (2021).
Four randomly selected groups of 20 Merino and 20 Dohne Merino ewes were mated to rams of their own breed or Dorper, Dormer or Ile de France rams to create two purebred control and six crossbred trial groups of offspring, namely purebred Dohne Merino, Dohne x Dorper, Dohne x Dormer, Dohne x Ile de France, purebred Merino, Merino x Dorper, Merino x Dormer and Merino x Ile de France. At weaning (100 days) four ram and four ewe lambs from each genotype were selected for the trial. They were housed in individual pens and adapted to a feedlot diet (16.38% CP, 11.44 MJ ME/kg) for a week. After adaptation they had ad libitum access to the diet throughout the year long trial period. Individual daily dry matter intakes (DFI), intake as percentage of body weight and cumulative feed intake as well as feed conversion ratio (FCR) were calculated for the growth period.
Trial data was analysed using the Non-linear estimation procedure in Statistica 14 (TIBCO Statistica, 2020) where intake and feed conversion data were fitted to various equations (linear, quadratic, power) to find the best-fitting models. Both age and live weight were evaluated as inputs for the models.
Throughout it was found that body weight was a better predictor of feed intake than age. Daily dry matter intake followed a quadratic pattern relative to weight, concurring with results from previous studies (Lewis and Emmans, 2020, 2010; van der Merwe et al., 2022). However, these models had low explicative value (R2 < 0.249) and were deemed unsuitable for practical use. Greater success was achieved modelling intake as percentage of body weight, with a linear relationship existing between the parameters. Between 32 and 57% of the variation in intake for the various genotypic groups was explained by the linear model, thus providing low to moderate predictive accuracy.
A linear model was found to best fit the cumulative intake data. Table 1 presents the model parameters for sex and genotype groups separately, since no interaction between sex and genotype was present. The B parameter differed between sexes, where ewes had a greater value than rams, indicating a higher rate of increase in feed consumption as the B parameter refers to the slope of a linear equation. This points to ewes having a poorer FCR than rams, as more feed is required per unit weight gained. Similar results were obtained by van der Merwe et al. (2022) for purebred sheep.
Table 1
Parameter values as least squares means (± S.E.) of the linear fitting of cumulative intake to weight for pure- and crossbred lambs from weaning until maturity
Main effects | | Parameter |
| | A | B |
Sex | Ram | -268.201 ± 10.517 | 6. 631 ± 0.157 |
| Ewe | -278.652 ± 10.517 | 7.672 ± 0.157 |
| P-value | 0.486 | 0.001 |
Genotype | Dohne Merino | -284.877 ± 21.034 | 7.493 ± 0.315 |
| Dohne x Dorper | -294.480 ± 21.034 | 7.107 ± 0.315 |
| Dohne x Dormer | -275.809 ± 21.034 | 6.702 ± 0.315 |
| Dohne x Ile de France | -257.354 ± 21.034 | 6.867 ± 0.315 |
| Merino | -290.365 ± 21.034 | 7.803 ± 0.315 |
| Merino x Dorper | -287.230 ± 21.034 | 7.308 ± 0.315 |
| Merino x Dormer | -239.534 ± 21.034 | 6.958 ± 0.315 |
| Merino x Ile de France | -257.764 ± 21.034 | 6.974 ± 0.315 |
| P-value | 0.533 | 0.258 |
CI = A + B(W) where CI is cumulative intake and W refers to body weight (both in kilogram)
With more than 82% of the variation in the intake data being explained the models have a high predictive accuracy, provided that the input weights fall within the target range of 30–110 kg. The cumulative feed intake models are suitable for practical application as they have a potentially high degree of predictive accuracy and use body weight, which is regularly monitored in feedlots, as input factor. One hurdle to practical application is that the conversion of cumulative intake to DFI at any point may prove problematic. Since cumulative intake is plotted relative to weight gain, the rate of weight gain will influence DFI (van der Merwe et al., 2022). To calculate DFI based on the prediction of cumulative intake, the growth rate at that point will first have to be known. However, since growth models are available for the genotypes in this study (Theron, 2021), this is achievable. The linear model also allows for a general comparison to be made between the FCR of various groups at different time points over a given time period. A correlation of 72% exists between the slope of the cumulative intake equation and the calculated average FCR over the growth period, thus allowing for a fairly accurate estimation of FCR.
Finally, the feed conversion ratios of the various groups were compared. A power equation (\(FCR=A\times {Weight}^{B}\)) provided the best fit for the data and results from this modelling are provided in Table 2. Approximately 55% of the variation in FCR could be explained by the fitting of the model.
Table 2
Parameter values as least squares means (± S.E.) of the power equation (\(FCR=A\times {Weight}^{B}\)) fitted to body weight to predict FCR for pure- and crossbred lambs from weaning until maturity
Genotype | Sex | Parameter | RMSE | R2 |
A | B |
Dohne Merino | Ram | 0.0263e ± 0.0337 | 1.384bc ± 0.109 | 0.241 | 0.573 |
| Ewe | 0.0641bcde ± 0.0337 | 1.193bcdef ± 0.109 | 0.256 | 0.488 |
Dohne x Dorper | Ram | 0.0421de ± 0.0337 | 1.300bcde ± 0.109 | 0.284 | 0.504 |
| Ewe | 0.0293de ± 0.0337 | 1.451ab ± 0.109 | 0.239 | 0.584 |
Dohne x Dormer | Ram | 0.0327de ± 0.0337 | 1.267bcde ± 0.109 | 0.235 | 0.592 |
| Ewe | 0.2600a ± 0.0337 | 0.806g ± 0.109 | 0.203 | 0.541 |
Dohne x Ile de France | Ram | 0.1224bcd ± 0.0337 | 1.071defg ± 0.109 | 0.257 | 0.522 |
| Ewe | 0.0292de ± 0.0337 | 1.355abcd ± 0.109 | 0.232 | 0.638 |
Merino | Ram | 0.0477cde ± 0.0337 | 1.386abc ± 0.109 | 0.333 | 0.482 |
| Ewe | 0.1400cb ± 0.0337 | 1.122cdef ± 0.109 | 0.279 | 0.478 |
Merino x Dorper | Ram | 0.0332de ± 0.0337 | 1.374abcd ± 0.109 | 0.227 | 0.651 |
| Ewe | 0.0310de ± 0.0337 | 1.440ab ± 0.109 | 0.267 | 0.523 |
Merino x Dormer | Ram | 0.0982bcde ± 0.0337 | 1.038efg ± 0.109 | 0.210 | 0.608 |
| Ewe | 0.0146e ± 0.0337 | 1.628a ± 0.109 | 0.230 | 0.687 |
Merino x Ile de France | Ram | 0.1567b ± 0.0337 | 0.898fg ± 0.109 | 0.210 | 0.568 |
| Ewe | 0.0496cde ± 0.0337 | 1.291bcde ± 0.109 | 0.246 | 0.475 |
P-value | Genotype | 0.018 | 0.014 | | |
Sex | 0.663 | 0.198 | | |
Interaction | < 0.001 | < 0.001 | | |
The purebred Merino group displayed the least favourable FCR, indicating that it is the least suitable for feedlot finishing of the genotypes included in the study. Most of the other groups displayed fairly similar FCR values until 60–70 kg, whereafter more prominent variations began occurring. Feed conversion ratio differed between rams and ewes as indicated by the slope parameter of the cumulative intake equation and previous studies (Kashan et al., 2005; Kashani and Bahari, 2017; van der Merwe et al., 2022). Crossbreeding did not provide a universal improvement in FCR. The crossbred Merino genotypes displayed an improvement in FCR relative to purebred Merinos but only the Dohne x Dormer group displayed a better FCR than the purebred Dohne Merinos. This contrast is supported by literature as no universal trend is reported. Kashan et al. (2005) and Schiller et al. (2015) reported that feed conversion improved in crossbred lambs while Kiyanzad (2002) and Khaldari and Ghiasi (2018) did not find any improvement.
In conclusion it can be stated that the objective of creating practically applicable models to predict feed intake for the genotypes included in this study has been achieved with the development of linear cumulative intake models. These cumulative intake models also allowed for an estimation of FCR over the entire growth period and the relative ranking of groups for feed efficiency. It appears that no significant differences exist between these genotypic groups with regards to feed intake and that crossbreeding did not provide universal improvements in FCR.