Next Article in Journal
The Effect of In Ovo Administration of Rosemary Essential Oil on Hatchability, Relative Hatching Weight, and Embryo Mortality Rate in Japanese Quail (Coturnix coturnix japonica)
Previous Article in Journal
Phylogenetic and Spatiotemporal Analyses of Porcine Epidemic Diarrhea Virus in Guangxi, China during 2017–2022
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementation of Feed Efficiency in Iranian Holstein Breeding Program

1
Department of Animal Sciences, College of Agriculture, Isfahan University of Technology, Isfahan 83111-84156, Iran
2
Department of Animal Science, University of Tehran, Karaj P.O. Box 3158711167-4111, Iran
3
Department of Animal Science, Faculty of Agriculture, Bu-Ali Sina University, Hamedan 65176-58978, Iran
4
Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, Beltsville, MD 20705-2350, USA
*
Author to whom correspondence should be addressed.
Animals 2023, 13(7), 1216; https://doi.org/10.3390/ani13071216
Submission received: 24 February 2023 / Revised: 25 March 2023 / Accepted: 27 March 2023 / Published: 31 March 2023
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:

Simple Summary

Including feed efficiency in breeding objectives for dairy cattle is desirable due to the potential benefit for increased sustainability and herd profitability. An investigation of selection indices that included the direct and indirect selection of feed efficiency traits showed that the best total response was observed when residual feed intake was directly included as a selection criterion. Using indicator traits (e.g., type) could be a useful proxy in the absence of direct genetic evaluations for feed efficiency. Breeding objectives of the current Iranian selection index, which currently focuses on future production and marketing environment, could include feed efficiency in future evaluations.

Abstract

This study aimed to evaluate the economic impact of improving feed efficiency on breeding objectives for Iranian Holsteins. Production and economic data from seven dairy herds were used to estimate the economic values of different traits, and a meta-analysis was conducted to analyze the genetic relationships between feed efficiency and other traits. Economic weights were calculated for various traits, with mean values per cow and per year across herds estimated at USD 0.34/kg for milk yield, USD 6.93/kg for fat yield, USD 5.53/kg for protein yield, USD −1.68/kg for dry matter intake, USD −1.70/kg for residual feed intake, USD 0.47/month for productive life, and USD −2.71/day for days open. The Iranian selection index was revised to improve feed efficiency, and the feed efficiency sub-index (FE$) introduced by the Holstein Association of the United States of America was adopted to reflect Iran’s economic and production systems. However, there were discrepancies between Iranian and US genetic coefficients in the sub-index, which could be attributed to differences in genetic and phenotypic parameters, as well as the economic value of each trait. More accurate estimates of economic values for each trait in FE$ could be obtained by collecting dry matter intake from Iranian herds and conducting genetic evaluations for residual feed intake.

1. Introduction

Feed costs account for 40 to 60% of the total costs associated with milk production [1]. For the dairy industry, collecting sufficient information on feed intake to enable the development and implementation of genetic evaluations for feed efficiency was and remains an expensive and difficult task [2]. However, efforts to include feed efficiency as part of dairy breeding objectives began in the 1990s because feed efficiency has a significant impact on the profitability of dairy production. Feed efficiency traits have been studied from various perspectives in the USA, Australia, Canada, Netherlands, Denmark, Sweden, and Finland, including estimation of economic value [3], breeding value [2], genetic evaluation [4,5], and finally inclusion in breeding objectives [6]. Improved feed efficiency is also desirable due to its potential benefits toward a sustainable environment by reducing greenhouse gas emissions per animal as well as farm profitability by reducing the input of natural resources for equivalent output production [7,8].
Different traits have been used to study genetic diversity for feed efficiency in dairy cattle. Feed intake measured as DMI is a key component for calculating feed efficiency in dairy cattle, and its heritability has been reported as being in a range from 0.08 to 0.52 [9,10,11]. Residual feed intake (RFI) is a valuable metric and preferred trait for feed efficiency; it is estimated from DMI and is defined as the difference between observed and predicted feed intake after accounting for body size, body weight change, and milk production [12]. With the concept of RFI [13], the definitions of residual methane (adjusting methane for production, feed, and weight) have been suggested [14,15,16] to explain express emission in terms of CH4 per kg of output, as CH4 intensity (MeI = g CH4/kg milk), or input such as methane yield (MeY = g CH4/kg feed). As RFI has been genetically adjusted for the main energy sinks (e.g., milk yield, BW), it is genetically uncorrelated with those energy sink traits. Milk production and BW-related traits are often part of the total merit index traits in the genetic evaluation of dairy cattle. Thus, after eliminating the genetic correlations of RFI with milk production and BW, RFI becomes a more independent trait for representing FE in the selection index [17]. The estimated heritability of RFI ranges from 0.14 to 0.38 in dairy cows [18,19,20]. The heritability of RFI depends largely on the environmental correlation between feed intake and milk production. RFI is dependent and independent of production, based on the phenotypic and genotypic regression of feed intake on milk production, respectively [21].
Breeding objectives and selection indices define the direction of selection in a livestock population. One of the most important tasks in developing an economic selection index is to estimate the economic values (EVs) of specific traits because these EVs depend on variety of revenue and expense variables [22]. Therefore, selection index coefficients for different traits must be quantified to assess an economic selection index. Calculation of RFI EVs in dairy cattle is complicated due to the dynamic changes in body weight (BW) and body condition score (BCS) during the lactation cycle [23]. The commonly used RFI model employs constant partial regressions of feed intake on energy sinks. Considering the complexity of the metabolic changes in cows across lactation stages, it is possible that the partial regression coefficients of feed intake on milk production, body maintenance, and BW change may differ across lactation stages. Therefore, the general RFI model with constant partial regression coefficients of feed intake on energy sinks may not always hold for the entire lactation period, influencing RFI estimation [17]. In addition, the lack of a cost-efficient method for measuring RFI on commercial dairy farms has limited its inclusion in breeding goals for dairy production systems. The main obstacle to integrating feed efficiency traits in breeding programs for dairy cattle is the restricted phenotypic database, which results in low genomic reliability. Thus, optimizing the collection of feed intake data from multiple sources is needed to increase the reference population size so that genomic prediction accuracy for feed efficiency can be improved [20].
The direct measurement of DMI and RFI currently is not possible in the Iranian dairy production system. Therefore, using economically important traits such as type traits with a moderate-to-high correlation with DMI and RFI is an alternative option. In addition, to improve genetic gain, Iran imports dairy bull semen from countries that include feed efficiency traits in their selection indices such as the USA, Canada, Germany, and the Netherlands.
Although feed efficiency traits are not documented in Iran, the genetic difference between cows for feed intake provides an opportunity to reduce costs through genetic selection. However, the estimation of economic coefficients for various traits of Iranian dairy cows is needed to form an economic selection index. The objectives of this study were (1) to investigate the genetic analysis of feed efficiency traits (DMI and RFI) of Holstein cows using meta-analysis, (2) to develop a bioeconomic model to estimate the EVs of DMI and RFI, and (3) to revise the current Iranian selection index for Holstein dairy cows to improve feed efficiency.

2. Materials and Methods

2.1. Meta-Analysis

Quantitative findings of numerous studies have been collected and integrated using meta-analysis techniques to learn more about the statistical significance and relevance of feed efficiency traits. This quantitative technique uses specific measures to indicate the strength of variable relationships for studies included in the analysis [24]. The genetic parameters used in this study (122 heritabilities and 95 genetic correlations) originated from a comprehensive review of 37 scientific studies that were published between 1991 and 2021. Those studies estimated heritability and genetic and phenotypic correlations for feed intake (DMI and RFI), production (milk, fat, and protein yields; MY, FY, and PY, respectively), reproduction (days open; DO), longevity (productive life; PL), type (stature, chest width, body depth, angularity, rump angle, and rump width), BW, and BCS traits of Holstein dairy cows in different countries.
Pooled Heritability. Koots’ method was used to calculate the weighted mean of heritability [25], and pooled heritability ( h pooled 2 ) was
h pooled 2 = i = 1 n ( h i 2 / SE h i 2 2 ) i = 1 n ( 1 / SE h i 2 2 )
where n is the number of published studies, h i 2 is the heritability estimate from published study i , and SE is standard error calculated as the weighted mean of standard error divided by the square root of the number of records; weighted mean of standard error was calculated as follows [26]:
SE = i = 1 n σ i 2 x i 2 i = 1 n x i
where σ i 2 is sampling variance and x i is a number of records.
Pooled Correlations. Meta-analysis of genetic correlations among traits of reported values in the literature was determined by transforming those values to an approximate normal scale using Fisher’s Z transformation [27]:
Z = 0.5 l o g ( 1 + r i j 1 r i j )
SE Z = 1 / ( n 3 )
Z = i = 1 x ( Z i / SE Z i 2 )
Z pooled = Z i = 1 x ( 1 / SE Z i 2 )
and
r g = ( e 2 Z 1 ) / ( e 2 Z + 1 )
where r i j is the genetic correlation between traits i and j; x is the number of records observed per trait; and Z is the weighted mean for the Z-transformed genetic correlations, which is used to estimate the pooled genetic correlation ( Z pooled ) and back-transformed to obtain the weighted mean genetic correlation ( r g ). This methodology assumes a genetic correlation of 1 between traits in similar target groups, such as production, reproduction, longevity, and feed intake.

2.2. Economic Weights

This study focused on seven Holstein dairy herds that were selected as typical examples of Iranian herds with respect to milk production, reproduction, culling, and productive life cows. These farms implemented a management, feeding, and housing system that is commonly used across the majority of Iranian dairy herds. All of the participating farms were systematically recorded by the Animal Breeding Center of Iran and were located in four Iranian provinces (Isfahan, Tehran, Qazvin, and Razavi Khorasan), which are known for housing the majority of industrial dairy farms. Table 4 provides descriptive statistics on production data and variables utilized to compute EVs. Market conditions in 2020 were the basis for generating the input parameters to derive EV. All costs and prices were converted from Iranian rial currency (IRR) to US currency (USD) using a conversion rate of 1 USD = 42,000 IRR.
The EVs were estimated based on the Iranian data from the 7 large Holstein dairy herds for each trait group (production, longevity, reproduction, and feed intake) using the bioeconomic model described in [28]. The economic weights (EWs) were calculated by multiplying the EV of each trait by cumulative discounted expressions (CDEs). The term CDEs refers to a correction for the accumulation of values over time, where “cumulative” refers to this accumulation. The term “discounted” is used to describe the process of adjusting future returns to reflect their present-day value [29].
The economic value of DMI was expressed in kg of dry matter intake in each lactation based on production and body size. The DMI was estimated as a function of DIM, BW, and corrected 3.2% milk fat [30]. The cost of 1 kg of DMI was calculated and compared using 4 different scenarios based on parity, production level, forage-to-concentrate (F:C) ratio, and concentrate source. Diets were formulated using Cornell Net Carbohydrate and Protein System software (CNCPS, version 5.0; Cornell University, Ithaca, NY, USA) [31]. The diet for the first parity was designed with a single F:C ratio of 45:55. For the second parity and later, the F:C ratio was 40:60, and the F:C ratio for high- (above 40 kg), moderate- (around 35 kg), and low-production (below 25 kg) cows was 35:65, 48:52, and 50:50, respectively. Two different concentrate sources were investigated for each of the cow groupings.
The definition of RFI and predicted daily DMI for dairy cows was determined using nutrient requirements of the Iranian production system with respect to net energy and metabolizable protein based on the equations found in [30]. The required amounts of net energy, metabolizable protein, and DMI were 4.55 (Mcal/kg of DMI), 121 (Mcal/kg of DMI), and 22.6 (kg) for a high-production cow with a BW of 680 kg and an average milk yield of 40 kg/day, respectively. Mean RFI was calculated from the difference between actual and predicted daily DMI. The EV for feed efficiency traits was defined as the partial derivative of the bioeconomic model with respect to the DMI and RFI traits included in the breeding objective.
The following relationship was used to calculate the economic value.
EV = T P high T P l o w T V high T V l o w
where T V h i g h and T V l o w represent increases and decreases in the value of trait i, respectively; T P high and T P l o w represent the total profit in the evaluated production system that was considered if the average attribute i increased or decreased, respectively.

2.3. Breeding Objectives

Traits of economic importance for MY, FY, PY, PL, DO, DMI, and RFI were defined in the Iranian dairy breeding program to explore 2 different breeding objectives ( H 1 includes DMI while H 2 includes RFI). Both H 1 and H 2 breeding objectives were considered to investigate the impact of including feed efficiency traits (DMI and RFI) on the genetic and economic gains.
Multiplying the EW of each trait by the EBV of its corresponding trait gives the aggregate breeding value (breeding objective):
H 1 = EW MY ( EBV MY ) + EW FY ( EBV FY ) + EW PY ( EBV PY ) + EW PL ( EBV PL ) + EW DO ( EBV DO ) + EW DMI ( EBV DMI )
H 2 = EW MY ( EBV MY ) + EW FY ( EBV FY ) + EW PY ( EBV PY ) + EW PL ( EBV PL ) + EW DO ( EBV DO ) + EW RFI ( EBV RFI )

2.4. Selection Index Scenarios

Several indices were established to examine the genetic and economic responses to selection. Effect of breeding objective on response per trait per generation was investigated using the selection index theory [32]. The optimal selection index coefficients (b) for a set of phenotypic traits maximize the correlation of phenotypic index I with breeding objectives H. To establish distinct selection criteria for each index scenario, the phenotypic values of MY, FY, PL, DO, DMI, and RFI for bulls with at least 60 daughters with a record for each trait were examined. The calculation of selection index with several traits required variability, heritability, relative EVs, and correlations among their phenotypes and genotypes.
A total of 9 different selection indices (I) have been chosen to investigate the direct and indirect response to selection based on heritabilities and their genetic correlations used from the meta-analysis for study feed efficiency traits. For the first scenario, 7 different selection indices based on the direct selection of DMI and RFI were developed and compared for simplicity and effectiveness in determining genetic and economic trends:
I 1 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b DO ( DO ) I 2 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b DO ( DO ) + b DMI ( DMI ) I 3 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b DO ( DO ) + b RFI ( RFI ) I 4 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b BW ( BW ) + b DMI ( DMI ) I 5 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b BCS ( BCS ) + b DMI ( DMI ) I 6 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b BW ( BW ) + b RFI ( RFI ) I 7 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b BCS ( BCS ) + b RFI ( RFI )
The b were calculated as b = P 1 G a , where P is a t × t phenotypic (co)variance matrix of correlated traits, G is a t × m genetic (co)variance matrix for m traits in the breeding objective with t correlated traits in the selection index incorporating the additive genetic relationships between sources of information, and a is an m × 1 vector of relative EWs in the breeding objective. The predicted response to selection for trait m in the breeding objective was calculated based on the method outlined by [33]
s ( b G ) / b G a
where s is selection intensity, which is assumed to be 1.4 and equivalent to the top 20% of sires selected, and b is a transposed vector of selection index weights. Accuracy of a selection index was also calculated as follows [33]:
b G a / a C a
where a is a vector of EW for m traits in a breeding objective and C is a m × m genetic covariance matrix among the traits in the breeding objective.
For the second scenario, indirect or indicator traits were included as selection criteria to improve feed efficiency. In commercial dairy herds, quantification of direct feed efficiency traits is challenging. Therefore, type traits that are more heritable and simpler to measure than DMI and RFI were investigated using meta-analysis to determine the genetic correlation between the existing type traits and feed efficiency traits. The 2 type traits with the highest genetic correlations with feed efficiency traits were chest width (CW) and angularity (ANG) with DMI and CW and stature (ST) with RFI. Those correlated type traits were included as indirect indicators of DMI and RFI in 2 additional selection indices:
I 8 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b DO ( DO ) + b CW ( CW ) + b ANG ( ANG ) I 9 = b MY ( MY ) + b FY ( FY ) + b PY ( PY ) + b PL ( PL ) + b DO ( DO ) + b CW ( CW ) + b ST ( ST )

2.5. Feed Efficiency Sub-Index

The methodology and equations applied by [34] for its feed efficiency sub-index (FE$) were adapted for Iranian economics and production systems. The FE$ index is defined as the net profit that a producer receives from increased production relative to PTA for MY, FY, PY, and feed saved (FSAV). The PTA for FSAV represents the expected weight of feed saved per lactation based on BW composite and RFI evaluations. For this study, FE$ was calculated as the sum of PTA for MY, FY, PY, and FSAV weighted for EV ( w ):
FE $ = w MY ( PTA MY ) + w FY ( PTA FY ) + w PY ( PTA PY ) + w FSAV ( PTA FSAV )
The w i for traits included in the Iranian FE$ formula were estimated as b i , FE $ ( EW i ) , where b i , FE $ is the genetic regression of FE$ on trait i and EW i is the EW for trait i under the Iranian production system. The b i , FE $ shows the change in the correlated trait as FE$ increases by 1 unit genetically and is calculated as follows:
  b i , FE $ = r i , FE $ h i 2 σ P i 2 h FE $ 2 σ P FE $ 2
where r i , FE $ is the genetic correlation between trait i and FE$; h i 2 and h FE $ 2 are heritabilities of trait i and FE$, respectively; and σ P i 2 and σ P FE $ 2 are phenotypic variance of trait i and FE$, respectively.
The US genetic evaluations for FE$, MY, FY, PY, and FSAV released in August 2021 for 10,852 Holstein AI bulls were obtained from the Alta Genetics (Watertown, WI) website (https://bullsearch.altagenetics.com/us/BS/Advanced/HO) to calculate Iranian FE$. To obtain Iran’s FE$ at the scale of the United States, the FE$ formula for both countries (USA and Iran) was multiplied by the genetic value of that trait ( PTA i ) in proportion to the estimated coefficients of the countries themselves. Finally, a regression was used to convert the PTA for FE$ of US bulls ( PTA FE $ USA ) to PTA for FE$ estimated under the Iranian production system ( PTA FE $ I r a n ) :
PTA FE $ Iran = β 0 + β 1 ( PTA FE $ USA )
where β 0 is the intercept and β 1 is the regression coefficient.

3. Results

3.1. Meta-Analysis

Table 1 shows the weighted mean heritabilities and summary statistics for feed efficiency, production, reproduction, longevity, type, BW, and BCS traits calculated by meta-analysis across dairy cattle breeds and countries from the published literature. For feed efficiency, the weighted mean heritability estimate for DMI (0.21) was higher than that for RFI (0.19). Weighted mean heritabilities estimated for production traits were moderate (0.22) and higher than for DO (0.03) and PL (0.10). Weighted mean heritabilities for type traits ranged from 0.23 for angularity to 0.46 for stature. The weighted mean heritabilities for body size traits were estimated to be 0.23 for BCS, and 0.58 for BW. The standard errors for weighted mean heritabilities of all traits were relatively low (<0.1).
Table 2 shows the weighted mean genetic correlations and summary statistics between feed efficiency traits (DMI and RFI) and production, reproduction, longevity, type, BW, and BCS traits calculated by meta-analysis for Holsteins across countries from the published literature. Estimated weighted mean genetic correlations between DMI and production traits were strongly positive and ranged from 0.51 (FY) to 0.68 (MY). A strong positive weighted mean genetic correlation (0.49) was found between DMI and PL, whereas a weakly negative weighted mean genetic correlation (−0.14) was found between DMI and DO. Estimated weighted mean genetic correlations of DMI with type traits were all positive (0.12 to 0.58), with the largest correlations for angularity (0.58) and chest width (0.55). For RFI, weighted mean genetic correlations with production traits were weak and ranged from −0.04 for PY to 0.08 for MY. Estimated weighted mean genetic correlation for RFI was strongly negative with DO (−0.50) and moderately negative with PL (−0.23). Weighted mean genetic correlations between RFI and type traits ranged from weakly negative (−0.06) to moderately positive (0.28), with the highest correlations for stature (0.28) and chest width (0.17).
Table 3 shows the heritabilities and phenotypic and genetic correlations of feed efficiency, production, reproduction, longevity, type, BW, and BCS traits used from meta-analysis to calculate optimal weights for selection criteria in each index.

3.2. Economic Weights

Mean values of production data, prices, and costs of each unit for computing EVs are in Table 4. EVs, genetic expressions (genetic improvement of genes in various traits at different times on different numbers of animals), and estimated EWs for each breeding objective trait are in Table 5. An increase in MY, FY, and PY by 1 kg/cow/year results in an average EWs of USD 0.34, USD 6.93, and USD 5.53, respectively (Table 5). Mean EWs for DO and PL were USD −2.71 day/cow/year and USD 0.47 month/cow/year, respectively. The EW for 1 kg DMI/cow/year varied between herds and ranged from USD −1.43 to USD −1.88 with a mean EW of USD −1.68/cow/year. For a reduction of 1 kg DMI/day during lactation, the mean EW for RFI was estimated to be USD −1.70/cow/year.

3.3. Selection Index Scenarios

The expected genetic and economic gains for each breeding objective trait using different selection indices are presented in Table 6. For the first scenario when a feed efficiency trait (DMI or RFI) was included in the selection index (I2 to I7), genetic and economic gains were generally higher for FY and DO compared with the baseline index (I1), whereas both gains were generally lower for PL. The total selection index response when DMI or RFI was included was equal to or higher than for I1 except for I2. Including DMI or RFI as a feed efficiency trait in I2 to I7 increased index accuracy by 1 to 8 percentage points compared to the baseline accuracy of I1.
For the second scenario in which type traits that were genetically correlated with feed efficiency were included in the selection index rather than DMI or RFI (I8 and I9), genetic and economic gains for MY, FY, PY, PL, and DO (Table 6) followed the same pattern as found for the first scenario when compared with I1. For DMI, higher genetic and economic gains for MY, FY, PL, and DO; lower gains for PY; higher index accuracy (2 percentage points); and higher total selection index response (USD 10/feed unit) were found for direct selection (I2) compared with indirect selection (I8). For RFI, higher genetic and economic gains for MY, FY, and DO; lower gains for PY and PL; higher index accuracy (1 percentage point); and equal total selection index response were found for direct selection (I3) compared to indirect selection (I9). Regardless of direct or indirect selection for feed efficiency, the total selection index response was generally higher if RFI rather than DMI was used in the selection index.

3.4. Feed Efficiency Sub-Index

The genetic regression coefficients for MY, FY, PY, and FSAV included in FE$ were estimated to be 0.015, 0.451, 0.518, and 1.000, respectively. Multiplying those coefficients by the EVs of the traits resulted in EW estimates of 0.002 for MY, 1.42 for FY, 1.27 for PY, and 0.762 for FSAV, as well as a conversion Equation (1) of PTA FE $ Iran = 44.24 + 1.262 ( PTA FE $ USA ) .

4. Discussion

4.1. Meta-Analysis

The meta-analysis in this study evaluated the extent of utilizable genetic variation in feed intake, production, reproduction, type, BW, and BCS traits of Holstein dairy cows. The heritability estimates reported in the studies were 0.10 for DMI [35], from 0.15 to 0.18 for 305-day feed intake traits [36], from 0.21 to 0.40 for DMI across lactations [10], 0.27 for lactation total DMI [37], and from 0.09 to 0.38 for RFI [14,38]. Those heritability estimates indicated that the regulation of feed intake has a substantial genetic component and can be successfully altered through selection [39]. Difficulty in accurately estimating genetic parameters for feed intake traits could result from the limited number of observations recorded per year or the lack of a suitable statistical model. As a reproductive trait, DO was generally reported with low heritability and strongly influenced by herd management, climate differences, and level of milk production [40,41,42,43]. As a longevity trait, PL is defined as the length of time from the first day of calving to culling or death of the cow. Although PL is primarily affected by environmental effects such as year and age at first calving, season, and herd, it also has a genetic component. Most studies reported low heritability for PL, but estimates depended on the model used for estimating genetic variation [44,45,46]. Due to the fact that heritability estimates for type traits generally ranged from moderate to high, indirect selection could be useful for economic traits of interest that have a strong genetic correlation with type traits [14,18,47,48,49]. Heritability estimates for BCS ranged from 0.17 to 0.34 [11,14,48,49]. Several studies reported that the highest heritability for BCS occurred when mean BCS was lowest; BCS tends to be lowest in early lactation when it is less influenced by differences in management practices [50,51,52]. For BW, heritability estimates ranged from 0.35 to 0.71 [11,14,36,38] and varied at different stages of a dairy cow’s life [53].
Meta-analysis revealed a weak genetic correlation of RFI with production traits, which is due to the fact that RFI is genetically independent of production traits, [13,14,18,54], indicating that the indirect selection for production traits will not improve RFI but will enhance DMI due to a strong genetic correlation of DMI with production traits [9,14,36]. However, the strong positive genetic correlation between RFI and DMI [14,38] suggests that selection of low RFI can improve the feed efficiency of dairy cattle. For reproductive traits, the low negative genetic correlation between DMI and DO (−0.14 ± 0.29) reported by [55] indicated that higher-yield cows had higher DMI and longer DO or lower reproductive performance. The same study reported a high negative genetic correlation between DO and RFI (−0.51 ± 0.76), indicating that cows with low RFI have longer DO and consequently poorer reproductive performance. However, standard errors for this genetic correlation were high, and additional research with a larger dataset is needed to confirm this interpretation. In contrast, Coleman et al. [56] reported no association between fertility traits and RFI. The genetic correlation between feed efficiency traits (DMI and RFI) and PL indicates that the selection of DMI against RFI can lead to longer stays of dairy cows in the herd. Although low heritability, voluntary culling, and post-deletion records can affect genetic gain for PL, this longevity trait has received much more attention in the past decades due to its high economic value and direct impact on profitability [57]. In addition, selecting dairy cattle for longevity traits in herds can reduce the cost of heifer replacements and increase the PL of high-yielding cows [58].
In most dairy cattle breeding programs, type traits are recorded because they are easy and affordable to measure and are heritable [59,60]. Several studies reported medium-to-high heritability for type traits [14,37,49,61], which agreed with the meta-analysis heritability estimate. Some of the type traits in this study, such as chest width, angularity, and stature, were used as indicator (or indirect) traits to select DMI and RFI in dairy cows due to the medium-to-high genetic correlations between those traits. The magnitude and direction of genetic correlation between body size traits (BCS and BW) and feed efficiency traits (DMI and RFI) showed that the genetic variation of DMI and RFI in dairy cows is explained better by BW than by BCS. Therefore, BW can be used as an indirect trait to select for DMI and RFI.

4.2. Economic Weights

Despite the same trait definition and pricing system for milk yield, the estimated EVs for production traits in this study differed from previous studies [28,62]. Sadeghi-Sefidmazgi et al. [28] reported EVs of USD 0.15/kg/cow/year for MY, USD 1.36/kg/cow/year for FY, and USD −1.02/kg/cow/year for PY in contrast to USD 0.34, USD 6.93, and USD 5.53, respectively (Table 5); Ghiasi et al. [62] reported an EW for MY of USD 0.193/kg/cow/year. The variation in the EVs of MY is the result of differences in the sale price of milk and production costs, including feed and labor costs. The EV of DO for all milk producing abilities was reported to be positive (from USD 0.21 to USD 0.40) when calving interval was increased from 12 to 13 month and negative (USD −0.04 to USD −0.23) when it was increased from 13 to 15 month [63]. Recently, the EV of DO was evaluated for Japanese Holstein cows and ranged from USD 3.59 to 7.72 [64], which is not similar to the estimated EV for DO of USD −2.71 in this study (Table 5). Sadeghi-Sefidmazgi et al. [28] reported an EV of USD 10.13/month for PL, which is higher than the EV reported in this study (USD 2.77/month/cow/year) (Table 5). Differences in maintenance costs as well as profits obtained through PL vary among different economic conditions around the world, which is the primary reason for differentiating EV from EW among studies.
To assess the economic importance of feed efficiency in the breeding objective for dairy cattle, the EWs of feed efficiency traits should be determined under current Iranian dairy production conditions. In this study, the different EWs for DMI estimated in different scenarios showed that the EW was highly sensitive to feed and animal factors influencing the feed intake of dairy cows as well as the difference between concentrate and roughage price. This sensitivity requires a consideration of production circumstances in defining breeding objectives.
The calculation of EVs for RFI was based on a daily DMI increase or decrease of 1 kg/day/animal. Therefore, the EV for RFI shows how the profit of a production system expressed per cow per year would decrease if the daily DMI of cows increased by 1 kg/day/cow over the whole rearing period, and the unit for RFI is monetary units/cow/year. Hietala et al. [65] calculated EWs for additional feed efficiency traits of Finnish Ayrshires, and the marginal EVs for RFI of breeding heifers and cows were estimated to be USD −30.6 and USD −66.96/kg DMI/day/cow/year, respectively. Krupová et al. [3] estimated the marginal EVs for RFI to be USD −66.18 and USD −65.568/kg DMI/day/cow/year for cows and breeding heifers in the dairy system, respectively, and USD −24.54, USD −13.56, and USD −7.248/kg DMI/day/cow/year for cows, breeding heifers, and fattened animals in the cow-calf system, respectively. In the 2021 US net merit index [6], the EV for RFI was estimated to be USD −0.14/kg of PTA RFI. The estimated EV for feed intake in this study was in the range of the estimated EV of other studies. The variation in the EVs for RFI can be attributed to differences in trait definition and expression of unit traits as well as assumptions related to management systems.

4.3. Selection Index Scenarios

In most countries, production traits were included as the initial economic traits in selection indices for dairy cattle. During the last 2 decades, breeding objectives have been extended to longevity, reproduction, health, and functionality traits for dairy cattle breeding programs worldwide [60,66,67]. The importance of integrating feed efficiency traits in Iranian breeding goals for genetic improvement has been determined due to the fact that studies in several countries have shown the substantial impact on the profitability of dairy production. Although data on feed intake in the dairy production system of Iran are missing and not recorded, forming a selection index using a meta-analysis may be valuable. Consequently, the importance of including two feed efficiency traits (DMI and RFI) in breeding objectives was investigated through meta-analysis in this study.
Based on the selection criteria defined for the first scenario, the direct selection for DMI (I2) resulted in higher economic and genetic gains for MY and FY than did the direct selection for RFI (I3), whereas the direct selection for RFI resulted in higher economic and genetic gains for PY, PL, and DO (Table 6). Due to the fact that RFI is genetically independent of production, its use as a selection criterion is equivalent to using restricted selection criteria to stabilize the production constant [21]. An unfavorable genetic and economic trend was observed for DO if either DMI or RFI were included in the selection index. The presence of antagonistic genetic correlations among traits included in an index is crucial for the response to selection. Additionally, the genetic improvement for DO relies on the heritability and economic weight of all traits incorporated in the index. Therefore, anticipating genetic gains based on a single affecting factor is not sufficient. Cows with low RFI can easily mobilize body reserves if changes in BCS are not included in the prediction of RFI or if BCS does not show all changes in body tissue mobilization. Due to a modest unfavorable genetic correlation between RFI and fertility, mobilization of body reserves has a negative effect on fertility [68,69,70,71]. The results of this study are consistent with the results of previous studies that reported unfavorable genetic correlations between RFI and fertility measures [55] and between DMI and DO [72]. Profit per feed unit (total selection index response) was higher when BW or BCS was substituted for DO in selection indices I4 through I7 than for the baseline index I1 (Table 6). Selection criteria of MY, FY, PY, PL, BW, and RFI were determined to result in the most profitable index (I6) when compared to others in the first scenario and used a suitable profit function to estimate the relative EW of the traits. Including the direct selection for RFI in I6 resulted in the highest total selection index response, which has also been reported by [73].
The direct selection for feed intake is not commonly used around the world because collecting feed intake data in most commercial herds is costly [74]. In several studies, phenotypic [75,76] and genetic [36,47,77,78] correlations between feed intake and type traits have been reported. An alternative to direct selection on feed intake data is to use indirect traits such as type traits to select for feed efficiency in dairy cows. Indirect selection on DMI or RFI was examined using I8 and I9, respectively. Genetic and economic gains for traits in the breeding objectives were generally higher for I8 than for I9 (Table 6). However, the total selection index response for I8 was USD 20/feed unit less than that for I9, and neither index had a total response higher than that for the baseline index.
The inclusion of RFI or DMI in breeding objectives in the absence of a robust genetic correlation between other breeding objective traits results in unfavorable selection responses. In a study of Australian dairy systems across breeds, increasing survival and reducing milk volume (i.e., increasing fat and protein content), live weight, DMI, SCC, and poor fertility were the only combination of traits that could increase net income and reduce greenhouse gas emissions per cow and per unit of milk [79]. In addition, investigation of the economic effect of genomic selection for RFI indicated that integrating genomic breeding values for RFI in the Australian Profit Ranking selection index would improve the annual rate of gain in profitability for Australian dairy cattle by 3.8% [80] FSAV has been included in two Australian selection indices (the Balance Performance Index, BPI; and the Health Weighted Index, HWI) since 2015 and Bolormaa et al. [81] reported that the reliability of lifetime RFI using trivariate analysis increased from 11% to 20% compared with the 2015 model, indicating that the genetic trend for FSAV is expected to shift in a favorable direction when FSAV is included in those indices. In our study, both the direct and indirect selection on RFI in a genetic–economic index were economically more profitable than direct or indirect selection on DMI. The use of RFI has been preferred because it accounts for milk production, body size traits, and parity [82] and is a potential target for genetic selection due to its connection to herd profitability and decreased environmental impact of dairy cattle. A study by Marinho et al. [83] showed that the most feed-efficient cows ate approximately 4 kg/day less dry matter and produced the same amount of energy-corrected milk (ECM) with improved reproductive performance compared with the least feed-efficient cows.

4.4. Feed Efficiency Sub-Index

In April 2021, the EWs for MY, FY, PY, and FSAV (the four main components in FE$) were estimated to be USD 0.0018/kg, USD 3.42/kg, USD 3.82/kg, and USD 0.24/kg, respectively [34]. The discrepancy between the Iranian and US EW for each trait in FE$ can be attributed to differences in genetic and phenotypic parameters and EV for each trait. To account for differences in Iranian and US dairy production and economic conditions, the PTA FE$ of US bulls should be multiplied by the conversion coefficients found for (1) to estimate Iranian PTA FE$. A more precise Iranian EV for each trait in FE$ can be determined by collecting DMI data from Iranian herds and initiating genetic evaluation for RFI.

5. Conclusions

Including feed efficiency in breeding objectives is desirable due to the potential benefit for increased sustainability and herd profitability. An investigation of different selection indices that included direct and indirect selection for feed efficiency traits showed that the best total selection index response was observed when RFI was directly included as a selection criterion. Although the inclusion of RFI or DMI in selection criteria is recommended, using indicator traits (e.g., type traits) could be a useful proxy to select for feed efficiency in the absence of direct genetic evaluation for feed efficiency traits. Breeding objectives of the current Iranian selection index, which currently focus on future production and marketing environment, should include feed efficiency in future evaluations.

Author Contributions

Conceptualization, A.S.-S. and S.T.; methodology, A.S.-S. and S.N.; software, A.S.-S. and S.N.; validation, A.S.-S., S.T., P.Z. and G.R.G.; formal analysis, S.N.; investigation, A.S.-S. and S.T.; resources, A.S.-S. and G.R.G.; data curation, S.N.; writing—original draft preparation, S.N.; writing—review and editing, A.S.-S. and S.T.; visualization, S.N.; supervision, A.S.-S., P.Z., G.R.G. and S.T.; project administration, A.S.-S.; funding acquisition, A.S.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Animal Care and Use Committee approval was not required for this study since the data used were obtained from an existing database.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are available upon request from the author Ali Sadeghi-Sefidmazgi: [email protected] and with permission from the Animal Breeding Center of Iran (ABCI).

Acknowledgments

The authors of the article would like to thank the herd managers for collaborating and generously providing the information needed. Special thanks are extended to Suzanne Hubbard for assistance in proofreading the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Connor, E.E. Invited review: Improving feed efficiency in dairy production: Challenges and possibilities. Animal 2015, 9, 395–408. [Google Scholar] [CrossRef]
  2. Pryce, J.E.; Gonzalez-Recio, O.; Nieuwhof, G.; Wales, W.J.; Coffey, M.P.; Hayes, B.J.; Goddard, M.E. Hot topic: Definition and implementation of a breeding value for feed efficiency in dairy cows. J. Dairy Sci. 2015, 98, 7340–7350. [Google Scholar] [CrossRef]
  3. Krupová, Z.; Krupa, E.; Michaličková, M.; Wolfová, M.; Kasarda, R. Economic values for health and feed efficiency traits of dual-purpose cattle in marginal areas. J. Dairy Sci. 2016, 99, 644–656. [Google Scholar] [CrossRef] [PubMed]
  4. de Jong, G.; van der Linde, R.; de Haas, Y.; Schopen, G.; Veerkamp, R. Genetic evaluation for feed intake in the Netherlands and Flanders, impact on efficiency and responses. In Proceedings of the 2016 Interbull Meeting, Puerto Varas, Chile, 24–28 October 2016. [Google Scholar]
  5. de Jong, G.; de Haas, Y.; Veerkamp, R.; Schopen, G.; Bouwmeester-Vosman, J.; van der Linde, R. Feed intake genetic evaluation: Progress and an index for saved feed cost. In Proceedings of the 2019 Interbull Meeting, Cincinnati, OH, USA, 23–26 June 2019; pp. 1–4. [Google Scholar]
  6. VanRaden, P.M.; Cole, J.B.; Neupane, M.; Toghiani, S.; Gaddis, K.P.; Tempelman, R.J. Net merit as a measure of lifetime profit: 2021 revision. In AIPL Research Reports NM$8 (05–21); USDA: Beltsville, MD, USA, 2021. [Google Scholar]
  7. Knapp, J.R.; Laur, G.L.; Vadas, P.A.; Weiss, W.P.; Tricarico, J.M. Invited review: Enteric methane in dairy cattle production: Quantifying the opportunities and impact of reducing emissions. J. Dairy Sci. 2014, 97, 3231–3261. [Google Scholar] [CrossRef] [PubMed]
  8. Hurley, A.M.; López-Villalobos, N.; McParland, S.; Lewis, E.; Kennedy, E.; O’Donovan, M.; Burke, J.L.; Berry, D.P. Genetics of alternative definitions of feed efficiency in grazing lactating dairy cows. J. Dairy Sci. 2017, 100, 5501–5514. [Google Scholar] [CrossRef] [PubMed]
  9. Persaud, P.; Simm, G.; Hlll, W.G. Genetic and phenotypic parameters for yield, food intake and efficiency of dairy cows fed ad libitum 1. Estimates for ‘total’lactation measures and their relationship with live-weight traits. Anim. Sci. 1991, 52, 435–444. [Google Scholar] [CrossRef]
  10. Manzanilla-Pech, C.I.V.; Veerkamp, R.F.; Calus, M.P.L.; Zom, R.; Van Knegsel, A.; Pryce, J.E.; De Haas, Y. Genetic parameters across lactation for feed intake, fat-and protein-corrected milk, and liveweight in first-parity Holstein cattle. J. Dairy Sci. 2014, 97, 5851–5862. [Google Scholar] [CrossRef]
  11. Köck, A.; Ledinek, M.; Gruber, L.; Steininger, F.; Fuerst-Waltl, B.; Egger-Danner, C. Genetic analysis of efficiency traits in Austrian dairy cattle and their relationships with body condition score and lameness. J. Dairy Sci. 2018, 101, 445–455. [Google Scholar] [CrossRef]
  12. Koch, R.M.; Swiger, L.A.; Chambers, D.; Gregory, K.E. Efficiency of feed use in beef cattle. J. Anim. Sci. 1963, 22, 486–494. [Google Scholar] [CrossRef]
  13. Berry, D.P.; Crowley, J.J. Cell Biology Symposium: Genetics of feed efficiency in dairy and beef cattle. J. Anim. Sci. 2013, 91, 1594–1613. [Google Scholar] [CrossRef]
  14. Manzanilla-Pech, C.I.V.; Veerkamp, R.F.; Tempelman, R.J.; Van Pelt, M.L.; Weigel, K.A.; VandeHaar, M.; Lawlor, T.J.; Spurlock, D.M.; Armentano, L.E.; Staples, C.R. Genetic parameters between feed-intake-related traits and conformation in 2 separate dairy populations—The Netherlands and United States. J. Dairy Sci. 2016, 99, 443–457. [Google Scholar] [CrossRef]
  15. Manzanilla-Pech, C.I.V.; Stephansen, R.B.; Difford, G.F.; Løvendahl, P.; Lassen, J. Selecting for feed efficient cows will help to reduce methane gas emissions. Front. Genet. 2022, 13, 885932. [Google Scholar] [CrossRef] [PubMed]
  16. Richardson, C.M.; Nguyen, T.T.T.; Abdelsayed, M.; Moate, P.J.; Williams, S.R.O.; Chud, T.C.S.; Schenkel, F.S.; Goddard, M.E.; van den Berg, I.; Cocks, B.G.; et al. Genetic parameters for methane emission traits in Australian dairy cows. J. Dairy Sci. 2021, 104, 539–549. [Google Scholar] [CrossRef] [PubMed]
  17. Li, B. Genetic Properties of Feed Efficiency and Related Traits in Dairy Cattle. Ph.D. Thesis, Swedish University of Agricultural Sciences, Uppsala, Sweden, 2018. [Google Scholar]
  18. Veerkamp, R.F.; Emmans, G.C. Sources of genetic variation in energetic efficiency of dairy cows. Livest. Prod. Sci. 1995, 44, 87–97. [Google Scholar] [CrossRef]
  19. Connor, E.E.; Hutchison, J.L.; Norman, H.D.; Olson, K.M.; Van Tassell, C.P.; Leith, J.M.; Baldwin, R. Use of residual feed intake in Holsteins during early lactation shows potential to improve feed efficiency through genetic selection. J. Anim. Sci. 2013, 91, 3978–3988. [Google Scholar] [CrossRef]
  20. Li, B.; VanRaden, P.M.; Guduk, E.; O’Connell, J.R.; Null, D.J.; Connor, E.E.; VandeHaar, M.J.; Tempelman, R.J.; Weigel, K.A.; Cole, J.B. Genomic prediction of residual feed intake in US Holstein dairy cattle. J. Dairy Sci. 2020, 103, 2477–2486. [Google Scholar] [CrossRef]
  21. Kennedy, B.W.; Van der Werf, J.H.; Meuwissen, T.H. Genetic and statistical properties of residual feed intake. J. Anim. Sci. 1993, 71, 3239–3250. [Google Scholar] [CrossRef]
  22. McArthur, A.T.G. Weighting breeding objectives—An economic approach. In Proceedings of the VI Annual Conference of Australian Association of Animal Breeding and Genetics, Perth, Western Australia, 19–21 July 1987; pp. 170–187. [Google Scholar]
  23. Williams, Y.J.; Pryce, J.E.; Grainger, C.; Wales, W.J.; Linden, N.; Porker, M.; Hayes, B.J. Variation in residual feed intake in Holstein-Friesian dairy heifers in southern Australia. J. Dairy Sci. 2011, 94, 4715–4725. [Google Scholar] [CrossRef]
  24. Shelby, L.B.; Vaske, J.J. Understanding meta-analysis: A review of the methodological literature. Leis. Sci. 2008, 30, 96–110. [Google Scholar] [CrossRef]
  25. Koots, K.R.; Gibson, J.P.; Smith, C.; Wilton, J.W. Analyses of published genetic parameter estimates for beef production traits. 1. Heritability. Anim. Breed. Abstr. 1994, 62, 309–338. [Google Scholar]
  26. Reckhow, D. Chapter XXIII: Parametric Statistics; Elsevier: Amsterdam, The Netherlands, 1999. [Google Scholar]
  27. Steel, R.G.D.; Torrie, J.H. Principles and Procedures of Statistics; McGraw-Hill: New York, NY, USA, 1960. [Google Scholar]
  28. Sadeghi-Sefidmazgi, A.; Moradi-Shahrbabak, M.; Nejati-Javaremi, A.; Miraei-Ashtiani, S.R.; Amer, P.R. Breeding objectives for Holstein dairy cattle in Iran. J. Dairy Sci. 2012, 95, 3406–3418. [Google Scholar] [CrossRef] [PubMed]
  29. Jiang, X.; Groen, A.F.; Brascamp, E.W. Discounted Expressions of Traits in Broiler Breeding Programs. Poult. Sci. 1999, 78, 307–316. [Google Scholar] [CrossRef]
  30. NRC. Nutrient Requirements of Dairy Cattle; National Academies Press: Washington, DC, USA, 2001. [Google Scholar]
  31. Fox, D.G.; Tedeschi, L.; Tylutki, T.; Russell, J.; Van Amburgh, M.; Chase, L.; Pell, A.; Overton, T. The Cornell Net Carbohydrate and Protein System model for evaluating herd nutrition and nutrient excretion. Anim. Feed Sci. Technol. 2004, 112, 29–78. [Google Scholar] [CrossRef]
  32. Hazel, L.N.; Lush, J.L. The efficiency of three methods of selection. J. Hered. 1942, 33, 393–399. [Google Scholar] [CrossRef]
  33. Dekkers, J.C. Design and Economics of Animal Breeding Strategies; Notes for Summer Short Course, 14–18 July; Iowa State University: Ames, IA, USA, 2003. [Google Scholar]
  34. Holstein Association USA. TPI Formula—April 2021. 2021. Available online: https://www.holsteinusa.com/genetic_evaluations/ss_tpi_formula.html (accessed on 9 September 2022).
  35. Hüttmann, H.; Stamer, E.; Junge, W.; Thaller, G.; Kalm, E. Analysis of feed intake and energy balance of high-yielding first lactating Holstein cows with fixed and random regression models. Animal 2009, 3, 181–188. [Google Scholar] [CrossRef]
  36. Vallimont, J.E.; Dechow, C.D.; Daubert, J.M.; Dekleva, M.W.; Blum, J.W.; Barlieb, C.M.; Liu, W.; Varga, G.A.; Heinrichs, A.J.; Baumrucker, C.R. Genetic parameters of feed intake, production, body weight, body condition score, and selected type traits of Holstein cows in commercial tie-stall barns. J. Dairy Sci. 2010, 93, 4892–4901. [Google Scholar] [CrossRef]
  37. Manafiazar, G.; Zimmerman, S.; Basarab, J.A. Repeatability and variability of short-term spot measurement of methane and carbon dioxide emissions from beef cattle using GreenFeed emissions monitoring system. Can. J. Anim. Sci. 2016, 97, 118–126. [Google Scholar] [CrossRef]
  38. Difford, G.F.; Løvendahl, P.; Veerkamp, R.F.; Bovenhuis, H.; Visker, M.; Lassen, J.; De Haas, Y. Can greenhouse gases in breath be used to genetically improve feed efficiency of dairy cows? J. Dairy Sci. 2020, 103, 2442–2459. [Google Scholar] [CrossRef] [PubMed]
  39. Shonka, B.N.; Tao, S.; Dahl, G.E.; Spurlock, D.M. Genetic regulation of prepartum dry matter intake in Holstein cows. J. Dairy Sci. 2015, 98, 8195–8200. [Google Scholar] [CrossRef]
  40. Haile-Mariam, M.; Morton, J.M.; Goddard, M.E. Estimates of genetic parameters for fertility traits of Australian Holstein-Friesian cattle. Anim. Sci. 2003, 76, 35–42. [Google Scholar] [CrossRef]
  41. El-Bayoumi, K.M.; El-Tarabany, M.S.; Abdel-Hamid, T.M.; Mikaeil, O.M. Heritability, genetic correlation and breeding value for some productive and reproductive traits in Holstein cows. Res. Opin. Anim. Vet. Sci. 2015, 5, 65–70. [Google Scholar]
  42. Ayalew, W.; Aliy, M.; Negussie, E. Estimation of genetic parameters of the productive and reproductive traits in Ethiopian Holstein using multi-trait models. Asian-Australas. J. Anim. Sci. 2017, 30, 1550–1556. [Google Scholar] [CrossRef]
  43. Brzáková, M.; Zavadilová, L.; Přibyl, J.; Pešek, P.; Kašná, E.; Kranjčevičová, A. Estimation of genetic parameters for female fertility traits in the Czech Holstein population. Czech J. Anim. Sci. 2019, 64, 199–206. [Google Scholar] [CrossRef]
  44. Roxström, A.; Ducrocq, V.; Strandberg, E. Survival analysis of longevity in dairy cattle on a lactation basis. Genet. Sel. Evol. 2003, 35, 305. [Google Scholar] [CrossRef] [PubMed]
  45. de Mello, F.; Kern, E.L.; Bretas, A. Longevity in Dairy Cattle. J. Adv. Dairy Res. 2014, 2, 126. [Google Scholar] [CrossRef]
  46. Imbayarwo-Chikosi, V.E.; Dzama, K.; Halimani, T.E.; Van Wyk, J.B.; Maiwashe, A.; Banga, C.B. Genetic prediction models and heritability estimates for functional longevity in dairy cattle. South Afr. J. Anim. Sci. 2015, 45, 105–121. [Google Scholar] [CrossRef]
  47. Veerkamp, R.F.; Brotherstone, S. Genetic correlations between linear type traits, food intake, live weight and condition score in Holstein Friesian dairy cattle. Anim. Sci. 1997, 64, 385–392. [Google Scholar] [CrossRef]
  48. Tsuruta, S.; Misztal, I.; Lawlor, T.J. Changing definition of productive life in US Holsteins: Effect on genetic correlations. J. Dairy Sci. 2005, 88, 1156–1165. [Google Scholar] [CrossRef]
  49. Bilal, G.; Cue, R.I.; Hayes, J.F. Genetic and phenotypic associations of type traits and body condition score with dry matter intake, milk yield, and number of breedings in first lactation Canadian Holstein cows. Can. J. Anim. Sci. 2016, 96, 434–447. [Google Scholar] [CrossRef]
  50. Koenen, E.P.C.; Veerkamp, R.F. Genetic covariance functions for live weight, condition score, and dry-matter intake measured at different lactation stages of Holstein Friesian heifers. Livest. Prod. Sci. 1998, 57, 67–77. [Google Scholar] [CrossRef]
  51. Jonest, H.E.; White, I.M.S.; Brotherstone, S. Genetic evaluation of Holstein Friesian sires for daughter condition-score changes using a random regression model. Anim. Sci. 1999, 68, 467–475. [Google Scholar] [CrossRef]
  52. Koenen, E.P.C.; Veerkamp, R.F.; Dobbelaar, P.; De Jong, G. Genetic analysis of body condition score of lactating Dutch Holstein and Red-and-White heifers. J. Dairy Sci. 2001, 84, 1265–1270. [Google Scholar] [CrossRef]
  53. Berry, D.; Buckley, F.; Dillon, P.; Evans, R.; Rath, M.; Veerkamp, R. Genetic parameters for level and change of body condition score and body weight in dairy cows. J. Dairy Sci. 2002, 85, 2030–2039. [Google Scholar] [CrossRef] [PubMed]
  54. Veerkamp, R.F.; Emmans, G.C.; Cromie, A.R.; Simm, G. Variance components for residual feed intake in dairy cows. Livest. Prod. Sci. 1995, 41, 111–120. [Google Scholar] [CrossRef]
  55. Vallimont, J.E.; Dechow, C.D.; Daubert, J.M.; Dekleva, M.W.; Blum, J.W.; Liu, W.; Varga, G.A.; Heinrichs, A.J.; Baumrucker, C.R. Feed utilization and its associations with fertility and productive life in 11 commercial Pennsylvania tie-stall herds. J. Dairy Sci. 2013, 96, 1251–1254. [Google Scholar] [CrossRef]
  56. Coleman, J.; Berry, D.P.; Pierce, K.M.; Brennan, A.; Horan, B. Dry matter intake and feed efficiency profiles of 3 genotypes of Holstein-Friesian within pasture-based systems of milk production. J. Dairy Sci. 2010, 93, 4318–4331. [Google Scholar] [CrossRef]
  57. Krupova, Z.; Wolfová, M.; Wolf, J.; Oravcová, M.; Margetín, M.; Peškovičová, D.; Krupa, E.; Daňo, J. Economic values for dairy sheep breeds in Slovakia. Asian-Australas. J. Anim. Sci. 2009, 22, 1693–1702. [Google Scholar] [CrossRef]
  58. VanRaden, P.M.; Wiggans, G.R. Productive life evaluations: Calculation, accuracy, and economic value. J. Dairy Sci. 1995, 78, 631–638. [Google Scholar] [CrossRef]
  59. Schaeffer, L.R. Effectiveness of model for cow evaluation intraherd. J. Dairy Sci. 1983, 66, 874–880. [Google Scholar] [CrossRef]
  60. Miglior, F.; Muir, B.L.; Van Doormaal, B.J. Selection indices in Holstein cattle of various countries. J. Dairy Sci. 2005, 88, 1255–1263. [Google Scholar] [CrossRef]
  61. Pryce, J.E.; Coffey, M.P.; Brotherstone, S. The genetic relationship between calving interval, body condition score and linear type and management traits in registered Holsteins. J. Dairy Sci. 2000, 83, 2664–2671. [Google Scholar] [CrossRef] [PubMed]
  62. Ghiasi, H.; Nejati-Javaremi, A.; Pakdel, A.; González-Recio, O. Selection strategies for fertility traits of Holstein cows in Iran. Livest. Sci. 2013, 152, 11–15. [Google Scholar] [CrossRef]
  63. Holmann, F.J.; Shumway, C.R.; Blake, R.W.; Schwart, R.B.; Sudweeks, E.M. Economic value of days open for Holstein cows of alternative milk yields with varying calving intervals. J. Dairy Sci. 1984, 67, 636–643. [Google Scholar] [CrossRef]
  64. Sasaki, O.; Takeda, H.; Nishiura, A. The economic value of days open in Holstein cows in Japan based on simulated changes in conception rate. Anim. Sci. J. 2020, 91, e13342. [Google Scholar] [CrossRef] [PubMed]
  65. Hietala, P.; Wolfová, M.; Wolf, J.; Kantanen, J.; Juga, J. Economic values of production and functional traits, including residual feed intake, in Finnish milk production. J. Dairy Sci. 2014, 97, 1092–1106. [Google Scholar] [CrossRef]
  66. Leitch, H.W. Comparison of international selection indices for dairy cattle breeding. In Proceedings of the Open Session of the Interbull Annual Meeting, Ottawa, ON, Canada, 5–6 August 1994. [Google Scholar]
  67. Kargo, M.; Hjortø, L.; Toivonen, M.; Eriksson, J.A.; Aamand, G.P.; Pedersen, J. Economic basis for the Nordic Total Merit index. J. Dairy Sci. 2014, 97, 7879–7888. [Google Scholar] [CrossRef]
  68. Pryce, J.E.; Coffey, M.P.; Simm, G. The relationship between body condition score and reproductive performance. J. Dairy Sci. 2001, 84, 1508–1515. [Google Scholar] [CrossRef]
  69. Berry, D.P.; Buckley, F.; Dillon, P.; Evans, R.D.; Rath, M.; Veerkamp, R.F. Genetic relationships among body condition score, body weight, milk yield, and fertility in dairy cows. J. Dairy Sci. 2003, 86, 2193–2204. [Google Scholar] [CrossRef]
  70. Dechow, C.D.; Rogers, G.W.; Klei, L.; Lawlor, T.J.; VanRaden, P.M. Body condition scores and dairy form evaluations as indicators of days open in US Holsteins. J. Dairy Sci. 2004, 87, 3534–3541. [Google Scholar] [CrossRef]
  71. Bastin, C.; Loker, S.; Gengler, N.; Sewalem, A.; Miglior, F. Genetic relationships between body condition score and reproduction traits in Canadian Holstein and Ayrshire first-parity cows. J. Dairy Sci. 2010, 93, 2215–2228. [Google Scholar] [CrossRef]
  72. Brito, L.F.; Oliveira, H.R.; Houlahan, K.; Fonseca, P.A.S.; Lam, S.; Butty, A.M.; Seymour, D.J.; Vargas, G.; Chud, T.C.S.; Silva, F.F. Genetic mechanisms underlying feed utilization and implementation of genomic selection for improved feed efficiency in dairy cattle. Can. J. Anim. Sci. 2020, 100, 587–604. [Google Scholar] [CrossRef]
  73. Houlahan, K.; Schenkel, F.S.; Hailemariam, D.; Lassen, J.; Kargo, M.; Cole, J.B.; Connor, E.E.; Wegmann, S.; Junior, O.; Miglior, F. Effects of incorporating dry matter intake and residual feed intake into a selection index for dairy cattle using deterministic modeling. Animals 2021, 11, 1157. [Google Scholar] [CrossRef] [PubMed]
  74. Løvendahl, P.; Ridder, C.; Friggens, N.C. Limits to prediction of energy balance from milk composition measures at individual cow level. J. Dairy Sci. 2010, 93, 1998–2006. [Google Scholar] [CrossRef] [PubMed]
  75. Gravert, H.O. Genetic factors controlling feed efficiency in dairy cows. Livest. Prod. Sci. 1985, 13, 87–99. [Google Scholar] [CrossRef]
  76. Sieber, M.; Freeman, A.E.; Kelley, D.H. Relationships between body measurements, body weight, and productivity in Holstein dairy cows. J. Dairy Sci. 1988, 71, 3437–3445. [Google Scholar] [CrossRef]
  77. Brotherstone, S. Genetic and phenotypic correlations between linear type traits and production traits in Holstein-Friesian dairy cattle. Anim. Sci. 1994, 59, 183–187. [Google Scholar] [CrossRef]
  78. Kennedy, B.W.; Dekkers, J.C.M.; Moore, R.K.; Jairath, L. Genetic and phenotypic parameter estimates between production, feed intake, feed efficiency, body weight and linear type traits in first lactation Holsteins. Can. J. Anim. Sci. 1999, 79, 425–431. [Google Scholar] [CrossRef]
  79. Bell, M.J.; Eckard, R.J.; Haile-Mariam, M.; Pryce, J.E. The effect of changing cow production and fitness traits on net income and greenhouse gas emissions from Australian dairy systems. J. Dairy Sci. 2013, 96, 7918–7931. [Google Scholar] [CrossRef]
  80. Hayes, B.J.; Van Der Werf, J.H.J.; Pryce, J.E. Economic benefit of genomic selection for residual feed intake (as a measure of feed conversion efficiency) in Australian dairy cattle. Recent Adv. Anim. Nutr. 2011, 18, 31–36. [Google Scholar]
  81. Bolormaa, S.; MacLeod, I.M.; Khansefid, M.; Marett, L.C.; Wales, W.J.; Miglior, F.; Baes, C.F.; Schenkel, F.S.; Connor, E.E.; Manzanilla-Pech, C.I.V.; et al. Sharing of either phenotypes or genetic variants can increase the accuracy of genomic prediction of feed efficiency. Genet. Sel. Evol. 2022, 54, 60. [Google Scholar] [CrossRef]
  82. Bach, A.; Terré, M.; Vidal, M. Symposium review: Decomposing efficiency of milk production and maximizing profit. J. Dairy Sci. 2020, 103, 5709–5725. [Google Scholar] [CrossRef] [PubMed]
  83. Marinho, M.N.; Zimpel, R.; Peñagaricano, F.; Santos, J.E.P. Assessing feed efficiency in early and mid lactation and its associations with performance and health in Holstein cows. J. Dairy Sci. 2021, 104, 5493–5507. [Google Scholar] [CrossRef] [PubMed]
Table 1. Weighted mean heritabilities and summary statistics for feed efficiency, production, reproduction, longevity, type, body weight, and body condition score traits calculated by meta-analysis across dairy cattle breeds and countries.
Table 1. Weighted mean heritabilities and summary statistics for feed efficiency, production, reproduction, longevity, type, body weight, and body condition score traits calculated by meta-analysis across dairy cattle breeds and countries.
TraitN 1Weighted Mean h2SERange of h2 in the Literature
DMI, kg/day10 (11)0.210.070.08 to 0.52
Residual feed intake, kg/day5 (7)0.190.090.09 to 0.38
Milk, kg21 (24)0.220.030.12 to 0.39
Fat, kg15 (16)0.220.050.06 to 0.25
Protein, kg15 (16)0.220.040.13 to 0.36
Productive life, month2 (2)0.100.030.01 to 0.11
Days open, day5 (4)0.030.010.002 to 0.03
Stature, cm/scale5 (7)0.460.010.40 to 0.60
Chest width, scale3 (5)0.260.010.19 to 0.31
Body depth, scale5 (7)0.390.010.28 to 0.41
Angularity, scale4 (6)0.230.010.21 to 0.29
Rump angle, scale3 (5)0.370.060.20 to 0.42
Rump width, scale2 (4)0.270.050.18 to 0.40
Body weight, kg5 (7)0.580.040.35 to 0.71
Body condition, score4 (5)0.230.020.17 to 0.34
1 Number of articles with number of databases in parentheses.
Table 2. Weighted mean genetic correlations and summary statistics between feed efficiency traits and production, reproduction, longevity, type, body weight, and body condition score traits calculated by meta-analysis for Holsteins across countries.
Table 2. Weighted mean genetic correlations and summary statistics between feed efficiency traits and production, reproduction, longevity, type, body weight, and body condition score traits calculated by meta-analysis for Holsteins across countries.
Feed Efficiency TraitCorrelated TraitN 1Weighted Mean Genetic CorrelationSERange of Genetic Correlations in the Literature
Dry matter intakeResidual feed intake3 (5)0.530.070.38 to 0.89
Milk yield8 (9)0.680.080.36 to 0.78
Fat yield4 (4)0.510.070.15 to 0.53
Protein yield4 (4)0.550.080.25 to 0.56
Productive life1 (2)0.490.080.48 to 0.51
Days open1 (2)−0.140.09−0.14 to −0.15
Stature5 (6)0.440.010.32 to 0.57
Chest width2 (3)0.550.090.45 to 0.68
Body depth4 (5)0.400.090.26 to 0.49
Angularity3 (4)0.580.09−0.02 to 0.60
Rump angle1 (2)0.150.040.10 to 0.21
Rump width2 (3)0.120.070.04 to 0.18
Body weight4 (6)0.520.010.35 to 0.71
Body condition score4 (5)0.360.08−0.04 to 0.71
Residual feed intakeMilk yield4 (5)0.080.05−0.05 to 0.35
Fat yield3 (3)−0.010.07−0.07 to 0.20
Protein yield3 (3)−0.040.09−0.03 to −0.06
Productive life1 (2)−0.230.09−0.22 to −0.24
Days open1 (2)−0.500.04−0.48 to −0.51
Stature2 (3)0.280.090.12 to 0.43
Chest width2 (3)0.170.090.06 to 0.39
Body depth2 (3)0.080.090.05 to 0.11
Angularity2 (3)0.140.090.08 to 0.41
Rump angle1 (2)−0.060.04−0.06 to −0.07
Rump width1 (2)−0.030.01−0.02 to −0.06
Body weight2 (3)0.150.070.03 to −0.26
Body condition score2 (3)−0.130.04−0.14 to 0.46
1 Number of articles with number of databases in parentheses.
Table 3. Heritabilities (diagonal, bold) and phenotypic (above diagonal), and genetic (below diagonal) correlations of feed efficiency, production, reproduction, longevity, type, body weight, and body condition score traits used from meta-analysis to calculate optimal selection index weights for selection criteria formulas.
Table 3. Heritabilities (diagonal, bold) and phenotypic (above diagonal), and genetic (below diagonal) correlations of feed efficiency, production, reproduction, longevity, type, body weight, and body condition score traits used from meta-analysis to calculate optimal selection index weights for selection criteria formulas.
Trait 1DMIRFIMYFYPYPLDOSTCWANGBWBCS
DMI0.210.750.660.600.670.030.090.220.350.240.35−0.02
RFI0.530.19−0.07−0.05−0.09−0.06−0.020.120.170.120.010.04
MY0.680.080.220.690.850.150.290.100.130.100.09−0.23
FY0.51−0.010.410.220.600.230.140.010.080.190.14−0.15
PY0.55−0.040.960.560.220.260.130.130.020.210.20−0.16
PL0.49−0.230.640.540.630.10−0.10−0.03−0.030.01−0.020.00
DO−0.14−0.500.420.320.35−0.600.030.08−0.030.01−0.07−0.17
ST0.440.280.120.100.12−0.190.170.46−0.010.060.380.08
CW0.550.170.120.100.12−0.190.170.050.26−0.110.450.24
ANG0.580.140.260.130.14−0.130.420.10−0.030.23−0.07−0.20
BW0.520.15−0.29−0.030.03−0.22−0.230.940.84−0.180.580.31
BCS0.36−0.13−0.30−0.27−0.31−0.48−0.23−0.130.72−0.650.850.23
1 DMI= dry matter intake, RFI = residual feed intake, MY = milk yield, FY = fat yield, PY = protein yield, PL = productive life, DO = days open, ST = stature, CW = chest width, ANG = angularity, BW = body weight, BCS = body condition score.
Table 4. Mean values of production data, prices, and costs of each unit of the variables considered for computing economic values.
Table 4. Mean values of production data, prices, and costs of each unit of the variables considered for computing economic values.
VariableMean
Production data
 305-day milk, kg11,750
 305-day fat, kg399.5
 305-day protein, kg364.3
 Longevity, year3.86
 Age at first calving, month24
 Calving interval, day420
 Stillbirth rate, %0.04
 Calf mortality, %0.01
 Productive cow mortality, %0.01
 Heat detection rate0.50
 Estrous cycle, day21
 Conception rate0.47
 Pregnancy rate0.23
 Live weight of culled calf, kg240
 Live weight of culled heifer, kg500
 Live weight of culled cow, kg680
 Proportion of male calves sold to a feedlot0.41
 Male and female calves reared until 3 month of age/cow/year, n0.83
 Female calves reared/cow/year, n0.40
 Involuntary culling rate of calves from birth until 3 month of age0.005
 Survival rate of heifers from 3 month of age until calving0.98
 Involuntary culling rate of heifers as a proportion of surplus heifers0.003
Prices
 Base milk, USD/kg1.07
 Accessory payment for milk fat, USD/kg0.14
 Accessory payment for milk protein, USD/kg0.14
 Male calf price, USD/calf1071
 Female calf price, USD/calf1429
 Replacement heifer, USD/heifer6667
 Culled calves, USD/kg7.14
 Culled heifers, USD/kg7.14
 Culled cows, USD/kg5.95
 Conventional domestic semen dose, USD12.14
 Conventional imported semen dose, USD24.05
 Sexed imported semen dose, USD80.95
 Insemination, USD26.89
Costs
 Base milk, USD/kg0.75
 Milk fat accessory, USD/kg7.36
 Milk protein accessory, USD/kg7.37
 NEL, USD/Mcal0.79
 MP, USD/kg0.002
 Calf rearing, USD/calf1000
 Rearing from 3 month of age until calving, USD/heifer4000
 Rearing from 3 to 21 month of age, USD/heifer3600
Table 5. Economic values, genetic expressions, and estimated economic weights for each breeding objective trait.
Table 5. Economic values, genetic expressions, and estimated economic weights for each breeding objective trait.
TraitEconomic Value
(USD/Unit of Trait/Cow/year)
Genetic ExpressionEconomic Weight
(USD/Unit of Trait)
Milk, kg0.341.000.34
Fat, kg6.931.006.93
Protein, kg5.531.005.53
Days open, day−2.711.00−2.71
Productive life, month2.770.170.47
DMI, kg−1.681.00−1.68
Residual feed intake, kg−1.701.00−1.70
Table 6. Expected genetic and economic gains for breeding objective traits based on 9 selection indices, index accuracy, and total index response.
Table 6. Expected genetic and economic gains for breeding objective traits based on 9 selection indices, index accuracy, and total index response.
TraitSelection Index 1
I1I2I3I4I5I6I7I8I9
Genetic gain (USD/unit gain)
Milk yield, kg658.88668.30664.40649.61655.34708.77657.67659.64656.76
Fat yield, kg12.3513.0312.8411.8012.4513.5712.4712.4512.53
Protein yield, kg16.9416.1616.6617.1917.0115.7816.9716.9216.86
Productive life, month3.072.582.613.133.062.903.132.912.90
Days open, day2.27−0.430.28−1.800.721.800.312.121.90
DMI, kg/day10.30−2.70−7.019.77
Residual feed intake, kg/day−0.77−5.26−7.03−2.10
Economic gain (USD/unit gain)
Milk yield, kg224.00227.20225.90220.90222.08241.00223.60224.30223.30
Fat yield, kg85.6490.3089.0081.7886.3094.0686.4786.2886.89
Protein yield, kg93.7089.3492.1595.0894.0787.3093.8793.6093.28
Productive life, month8.507.137.288.688.498.038.678.078.05
Days open, day−6.141.16−0.784.88−1.96−4.88−0.85−5.75−5.17
DMI, kg/day−17.314.5311.78−16.42
Residual feed intake, kg/day1.318.9511.973.58
Index accuracy0.890.920.910.960.970.960.930.900.90
Total selection index response, (USD/unit gain)410400410420420430420390410
1 Traits included in indices: I1 (milk yield, fat yield, protein yield, productive life, and days open), I2 (milk yield, fat yield, protein yield, productive life, days open, and DMI), I3 (milk yield, fat yield, protein yield, productive life, days open, and residual feed intake), I4 (milk yield, fat yield, protein yield, productive life, body weight, and DMI), I5 (milk yield, fat yield, protein yield, productive life, body condition score, and DMI), I6 (milk yield, fat yield, protein yield, productive life, body weight, and residual feed intake), I7 (milk yield, fat yield, protein yield, productive life, body condition score, and residual feed intake), I8 (milk yield, fat yield, protein yield, productive life, days open, chest width, and angularity), and I9 (milk yield, fat yield, protein yield, productive life, days open, chest width, and stature).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nadri, S.; Sadeghi-Sefidmazgi, A.; Zamani, P.; Ghorbani, G.R.; Toghiani, S. Implementation of Feed Efficiency in Iranian Holstein Breeding Program. Animals 2023, 13, 1216. https://doi.org/10.3390/ani13071216

AMA Style

Nadri S, Sadeghi-Sefidmazgi A, Zamani P, Ghorbani GR, Toghiani S. Implementation of Feed Efficiency in Iranian Holstein Breeding Program. Animals. 2023; 13(7):1216. https://doi.org/10.3390/ani13071216

Chicago/Turabian Style

Nadri, Sara, Ali Sadeghi-Sefidmazgi, Pouya Zamani, Gholam Reza Ghorbani, and Sajjad Toghiani. 2023. "Implementation of Feed Efficiency in Iranian Holstein Breeding Program" Animals 13, no. 7: 1216. https://doi.org/10.3390/ani13071216

APA Style

Nadri, S., Sadeghi-Sefidmazgi, A., Zamani, P., Ghorbani, G. R., & Toghiani, S. (2023). Implementation of Feed Efficiency in Iranian Holstein Breeding Program. Animals, 13(7), 1216. https://doi.org/10.3390/ani13071216

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop