The Comparison of Four Economical Selection Indices for Improving the Performance of Kermani Sheep under Rural Production System


1 Department of Animal Science, Kerman Agricultural and Natural Resources Research and Education Center, Agricultural Research, Education and Extension Organization (AREEO), Kerman, Iran

2 Department of Animal Science, Faculty of Agriculture, University of Jiroft, Jiroft, Iran


Developing effective selection programs for improving the performance of livestock requires taking the economical selection indices into account. The present study was performed for comparing four selection indices developed in Kermani sheep through simulated data. The relative importance of traits was determined based on the estimated economic values and consequently, the most beneficial traits were applied to develop the selection indices. The considered selection indices were different combinations of traits including ewe body weight (EBW), annual wool weight (AWW) and total weaning weight of lambs per exposed ewe (TWWEE). The first index included all three mentioned traits (EBW, AWW, and TWWEE). The second, third and fourth indices were included (AWW, TWWEE), (EBW, TWWEE) and (TWWEE), respectively. In general, the obtained results indicated that under all the developed selection indices the aggregate genotype, selection index and economic progress increased by increasing in the population size and decreasing in ram ratio, but resulted in reduced inbreeding average. In more cases, the aggregate genotype decreased by removing a trait from the index. The comparisons among indices indicated that the most suitable selection index for this breed under rural production system is the first index which included EBW, AWW, and TWWEE. This index had a maximum selection index average, aggregate genotype, and economic progress.



Kermani is a meat and wool type, fat-tailed, medium-sized and white-wool sheep breed, indigenous to the south-eastern areas of Iran and appropriately adapted to the unfavorable climatic conditions and pastures of low quantity and quality which is prevalent in the region (Mohammadabadi et al. 2017). The wool of Kermani sheep is fairly thick with a carpet grade.Mutton is the main source of red meat in Iran and now meat production from the sheep does not sufficient for increasing consumer demand. The primary objective of a breeding program for any livestock species is to maximize the rate of genetic progress for economically important traits (Rashidi et al. 2008) such as body weight and reproductive performance.Increasing the profitability of any sheep production system is of great importance (Harris, 1970). Hossein-Zadeh and Ardalan (2008) pointed out that meat production is one of the primary goals in sheep breeding of Iran. Ghiasi et al. (2016) reported that the number of lambs produced by each ewe constitutes the main part of income while milk production is of secondary importance and wool has no income or is little. Falconer and Mackay (1996) recommended the application of selection indices for multi-trait selection in animal populations. Selection indices incorporate data on biological and also economic aspects of production simultaneously into a single value known as aggregate genotype (Parish et al. 2011; Paswan et al. 2016). Various breeding objectives and selection indices have been proposed in several livestock species such as indigenous chickens in Kenya (Okeno et al. 2012), Rayeni Cashmere goat (Kargar Borzi et al. 2017), wormsilk (Seidavi et al. 2008), buffalo (Agudelo-Gómez et al. 2016) and Aberdeen Angus cattle (Campos et al. 2014). In Iran, under rural production systems, most sheep breeding enterprises may be unfavorably affected by the factors such as limited numbers of breeding males which increase the level of inbreeding in the flock. In order to increase the genetic progress and decrease the rate of inbreeding, the establishment of a suitable selection index considering factors such as flock size, ram ratio and traits of interest is necessary. As Kermani sheep breed mainly reared under the rural system, there is no information on the suitable selection indices for improving the performance of flocks under such systems. Up to our knowledge, there is no information on such selection indices. Therefore, the present study was performed for comparing different selection indices in Kermani sheep under rural production systems.



Breeding objective and selection criteria

Different areas of the Kerman province, which located in the south-eastern part of Iran, were chosen for gathering the required information. Five farmers were considered finally to take part in the present investigation. All farmers expected to increase their economic profit from their sheep flocks through higher meat and wool output and also reduced costs. Therefore, higher reproduction rate, higher growth rate and higher wool weight of the sheep are of primary interest. To select animals for the increased quantity of meat and wool produced, ewe body weight, annual wool weight and total weaning weights of lambs per exposed ewe, were used as selection criteria. The breeding objectives of rural production systems in Kermani sheep were acquired from economic analysis and the conversation with flock holders. So, breeding goals (meat production and wool production) and selection criteria (body weight (kg), total weaning weight for each exposed ewe (kg) and wool production (kg)) for Kermani sheep were considered.


Structure of (co)variance matrix

The phenotypic and genetic (co)variance matrix (Table 1) was derived from the data of experimental breeding station of Kermani sheep, located in Shahrbabk, southeast of Iran, which collected over a period of 18 years (1993-2011). (co)variance components for the traits were estimated by the Average Information Restricted Maximum Likelihood (AI-REML) method fitting an animal model. For this purpose, the WOMBAT computer program was used (Meyer, 2006). The required information which was accessible from those records also was reported in some research (Bahreini Behzadi et al. 2007; Mokhtari et al. 2008; Mokhtari et al. 2010).


Economic analysis

Microsoft excel spreadsheets were used to estimate productive and reproductive efficiency as well as of the costs and incomes. The systems, production costs and revenues, and profit equation were estimated as explained by Kosgey et al. (2003). In the present study, a deterministic static model was used for calculating the economic values (EVs) of decisive traits of Kermani sheep. The total annual profit of the flock was derived as the difference between costs and revenues of the system. The average prices in 2011 were used and all costs and prices were expressed in Dollar. The productive and time units were the ewe and year, respectively. The inputs for the production system were feed, management, and fixed costs. The outputs were revenues from the sale of cull ewes and rams, excess lambs and wool. The input parameters were taken from the five Kermani sheep flocks with a total size of 635 heads. The economic value of each trait (Vi) was obtained by Eq. (1)

Vi= Pi – P                                 (1)


P and Pi: profits before an increase and after increase of the trait by 1%, keeping all other traits at their mean value simultaneously (Lobo et al. 2011).


Data simulation and flock structure

To determine the appropriate selection index for Kermani sheep, the productive population was simulated. For this purpose, the base population was randomly established based on genetic and phenotypic (co)variance matrices (Table 1) of Kermani sheep and the means of the studied flocks (Table 2) by applying the Visual Basic 6.0 programming language. The assumptions and the statistical distribution of the simulated parameters are presented in Table 3. After simulating the phenotypic records of the base population, for obtaining the phenotypic records of the offspring random mating was done between top adult rams and ewes (the best parents) as follows:

y0= µ + 0.5gb + 0.5gd + m + e

m= 0.5 (2-Fb-Fd) 0.5ZLG


y0: vector of offspring phenotypic records.

µ: vector of mean of the trait considered.

gb and gd: vectors of breeding values for rams and ewes respectively.

m: vector of Mendelian sampling.

e: vector of residuals.

Z: vector of random numbers normally distributed with a mean of zero and variance equal to one.

LG: lower triangular matrix obtained from Chalsky analysis G.

Fb and Fd: inbreeding coefficients for rams and ewes, respectively.

Top adult rams and ewes were selected based on economic selection index. To choose the best parents, the selection index for each animal was established. Flock size, ram ratio, numbers of generations, numbers of traits, litter size, percentage of ram and ewe in the flock, longevity of rams and ewes in the flock, age of the animal at maturity and mortality rate were defined similarly to the considered Kermani sheep flocks (Table 2). The overlap between generations was considered. Selection indices were derived by changing the ram ratio from 2 to 8 percent.


Developing the selection indices

Four selection indices with different herd size and ram ratio were proposed (I1–I4) to compare genetic and economic gain and inbreeding rate with computer simulation after selected 10 years. The traits included in indices were: ewe body weight (EBW), annual wool weight (AWW) and total weaning weight of lambs per exposed ewe (TWWEE). The first index (I1) was considered as a full selection index that included all the considered traits. To decrease the cost of trait measurement, other selection indices were proposed to detect the appropriate selection index. The traits contained in each index were as follows: EBW, AWW, and TWWEE (I1); AWW and TWWEE (I2); EBW and TWWEE (I3) and TWWEE (I4).

Economic gain under each index was defined as:

E= ∑(EVi×BVi)     (2)


EVi: absolute economic value of the ith trait involved in the index.

BVi: breeding value for the ith trait included in the index.

Aggregate genotype (H) was defined as:

H= ∑vi × gi     (3)


vi: relative economic value of the ith trait in the aggregate genotype.

gi: breeding value of the ith trait in the aggregate genotype.



Economic analysis

The absolute and relative economic values of studied traits are shown in Table 1. Relative economic values of EBW, AWW, and TWWEE traits were - 1.95, 1.00 and 14.48, respectively. The calculated economic values, exceptEBW, were positive. Positive economic value for a trait indicated that genetic progress in that trait would positively influence on the profitability of the breeding system. Negative economic values for body weight traits were alsoreported in literature (Bett et al. 2011; Kargar Borzi et al. 2017). Economic increase comes from higher meat and wool output as well as reduction of costs in sheep flocks. Hence, higher reproduction and growth rate and wool weight are economically important. The economic value for total weaning weight of lambs per exposed ewe was larger than the other traits (Table 1), implied that TWWEE has higher influence on the costs of production system than other studied traits.


The comparison of selection indices

The values of inbreeding under the considered selection indices (I1 to I4), different flock size and ram ratio are shown in Table 4. The obtained results showed that inbreeding was reduced by increased flock size and ram ratio, accordingly. The lowest value of inbreeding was obtained under I4 which included TWWEE in flock size of 300 and ram ratio of 0.08. Kosgey et al. (2003) reported that the increase in flock size would be associated with decreased genetic gain and inbreeding. Therefore, for preventing the increase of inbreeding in small-sized flocks, the ram ratio should be increased (Kosgey et al. 2003). As shown in Table 4, the maximum value for selection index (183.54) was obtained under the first index, I1, with flock size of 200 and ram ration of 0.02 while the minimum one (135.10) was obtained under I2, flock size of 100 and ram ratio of 0.08. Inbreeding ranged from 0.052 (under I4, flock size of 300 and ram ratio of 0.08) to 0.288 (under I1, flock size of 100 and ram ratio of 0.02). The values of aggregate genotype and economic gain under different flock size and ram ratio across four selection indices (I1 to I4) are shown in Table 5.


Table 1 Estimates of mean, economic values and structure of (co)variance matrix for Kermani sheep1


1 Genetic variance in diagonal and residual variance in diagonal (in parentheses); residual (below diagonal) and genetic (above diagonal) covariances.

EBW: ewe body weight; AWW: annual wool weight and TWWEE: total weaning weight for each exposed ewe.

SD: standard deviation.


Table 2 Overview of the assumed values for the input variables of the study



Table 3 The assumptions and statistical distribution of the simulated parameters


y: vector of phenotypic records; µ: vector of mean of traits; g: vectors of additive genetic value for traits; e: vector of residuals; G: matrix of additive genetic (co)variance; R: residuals (co)variance; LG and LR: lower triangular matrix obtained from Chalsky analysis G and R respectively and Z1 and Z2: vectors of random numbers normally distributed with a mean of zero and variance equal to one.


Table 4 The values of inbreeding and selection index under different flock size and ram ratios of Kermani sheep


I1: included ewe body weight (EBW), annual wool weight (AWW), and total weaning weight of lambs per exposed ewe (TWWEE); I2: included annual wool weight (AWW) and total weaning weight of lambs per exposed ewe (TWWEE); I3: included ewe body weight (EBW) and total weaning weight of lambs per exposed ewe (TWWEE); I4: included total weaning weight of lambs per exposed ewe (TWWEE).


Under the first selection index (I1),the highest amount of aggregate genotype value (183.64) and economic gain (123.65 $) was obtained. By considering the second selection index (I2), aggregate genotype value and economic gain were varied between 136.10 to 169.12 and 93.79 to 115.57 $, respectively. Under the third selection index (I3), a minimum amount of aggregate genotype value (136.04), and economic gain (86.39 $) was obtained. Under the fourth selection index (I4), aggregate genotype value and economic gain were varied between 145.45 to 175.20 and 79.78 to 114.73 $, respectively. Under all developed selection indices, in flock size of 300 and ram ratio of 0.02 the highest amount of aggregate genotype value and economic gain was achieved. A reduction in ram percent will precipitate genetic progress because of the high selection intensity (Falconer and Mackay, 1996). Economic values were positive for AWW and TWWEE but not for EBW, implying that by increasing body weight of ewes, the profit would be reduced. Obtained negative economic values for ewe body weight were in agreement with the result obtained by Bett et al. (2011). Also, Conington et al. (2004) reported a negative economic value for mature body weight. The highest economic value was obtained for TWWEE. Likewise, Gallivan (1996) pointed out that reproductive traits have high economic value in Canadian sheep. The main influencing factors to be considered in genetic selection strategies are the accuracy of selection, selection intensity, effective population size and mating system (Falconer and Mackay, 1996). Therefore, optimum response to selection may be achieved by maximizing the accuracy of selection, selection intensity and additive genetic variance or effective population size (Falconer and Mackay, 1996). But, all of these factors cannot be maximized, with limited resources. For example, increasing selection intensity decreases effective population size, which consequently resulted in a reduction in response to selection. Likewise, increasing the accuracy of selection by use of selection indices also reduces effective population sizes, which lead to a quicker increase in the rate of inbreeding (Quinton and Smith, 1995). Inbreeding rates are speed up in livestock, and inbreeding causes economic losses through depression in production, growth, health, and fertility (Weigel, 2001). To increase genetic progress and limit the rate of inbreeding, establishment of appropriate selection index based on different herd size and ram ratio were the purpose of this research. The obtained results (Table 4) showed that the mean of inbreeding during 10 generations of selection was reduced with an increase in flock size and ram ratio. The increase in inbreeding negatively influences on genetic variance and decrease the genetic gain, accordingly (Falconer and Mackay, 1996). Therefore, for reducing the mean of inbreeding, large flock size and the adequate number of rams for mating are required. It is obvious that adequate number of rams for mating can reduce effectively the rate of inbreeding over generations with a minimal decrease of selection response. Therefore, the results of the present study will guide Kermani sheep breeders focus on a scheme to decide the starting numbers of basis animals to create a breeding herd.


Table 5 The values of aggregate genotype and economic gain ($) under different flock size and ram ratios of Kermani sheep


I1: included ewe body weight (EBW), annual wool weight (AWW), and total weaning weight of lambs per exposed ewe (TWWEE); I2: included annual wool weight (AWW) and total weaning weight of lambs per exposed ewe (TWWEE); I3: included ewe body weight (EBW) and total weaning weight of lambs per exposed ewe (TWWEE); I4: included total weaning weight of lambs per exposed ewe (TWWEE).


To prevent excessive inbreeding in the flock, small-sized flocks may be combined. The increase of flock size will cause a reduction of inbreeding, and increase of genetic and economic gain, accordingly. Then, decrease of ram up to 2 % and increase of flock size up to 300 may be suggested as the most suitable breeding program for Kermani sheep under extensive production system. The annual response to selection applying a multi-trait selection index will be smaller in each component trait than the selection accomplished for that trait individually. However, it is important to evaluate all traits with economic value, to improve genetic progress and the greatest rate of change in overall economic merit. Selection for reproductive traits should be based on a criterion that is highly similar to the breeding goals; for example, if the selection is based on the number of lambs, the frequency of genes affecting the total weight of the lamb is not sufficiently affected. Snyman et al. (1996) indicated that the selection for the number of lambs born per ewe is only significant for the number of lambs and will not increase the weight of lambs; while the selection for total weaning weight per ewe will also increase individual weight gain of lambs. Therefore, traits that were considered in the full selection index included EBW, AWW, and TWWEE. Response to selection increased with an increase in herd size. However, with an increase in ram ratio, response to selection decreased (Tables 4 and 5). The highest genetic and economic responses were obtained under the first index (I1). Under the second selection index (I2), the least genetic and economic responses to selection were obtained. Among the developed selection indices in the present study, the first one, which included all the three considered traits into account, may be taken as the best one for flock holders of Kermani sheep. Anyhow, if meat quantity to be the first interest of the flock holders of Kermani sheep, the third index (I3) which combines ewe body weight and total weaning weight of lambs per exposed ewe may be an appropriate alternative.



Results of the present study indicated that, for the genetic improvement of Kermani sheep populations, Index 1, with all traits including (ewe body weight, annual wool weight and total weaning weight for each exposed ewe), may be used appropriately. The choice and use of the index depending on the definition of objectives. The results of the study can be used to design local sheep breeding programs in Iran, and in other countries in low-input production systems with low technological progress, such that those found in the tropics.



We are grateful to the farmers that participated in the study, answered all questions and allowed us to work with them and their herds.

Agudelo-Gómez D.A., Agudelo-Trujillo J.H. and Cerón-Muñoz M.F. (2016). Selection index for meat and milk traits of buffaloes in Colombia. Livest. Sci. 191, 6-11.
Bahreini Behzadi M.R., Shahroudi F.E. and Van Vleck L.D. (2007). Estimates of genetic parameters for growth traits in Kermani sheep. J. Anim. Breed. Genet. 124, 296-307.
Bett R.C., Kosgey I.S., Kahi A.K. and Peters K.J. (2011). Definition of breeding objectives and optimum crossbreeding levels for goats in the smallholder production systems. Small Rumin. Res. 96, 16-24.
Campos G.S., Neto B.J., Oaigen R.P., Cardoso F.F., Cobuci J.A., Kern E.L., Campos L.T., Bertoli C.D. and McManus C.M. (2014). Bioeconomic model and selection indices in Aberdeen Angus cattle. J. Anim. Breed. Genet. 131, 305-312.
Conington J., Bishop S.C., Waterhouse A. and Simm G. (2004). A bioeconomic approach to derive economic values for pasture-based sheep genetic improvement programs. J. Anim. Sci. 82, 1290-1304
Falconer D.S. and Mackay T.F.C. (1996). Introduction to Quantitative Genetics. Longmans Green, Harlow, Essex, United Kingdom.
Gallivan C. (1996). Breeding objectives and selection index for genetic improvement of Canadian sheep. Ph D. Thesis. University of Guelph, Canada.
Ghiasi H., Khaldari M. and Taherkhani R. (2016). Comparison of Ordinary and restricted selection indices on genetic gain of body weight in Iranian Moghani sheep. Iranian J. Appl. Anim. Sci. 6, 369-373. 
Harris D.L. (1970). Breeding for efficiency in livestock production: Defining the economic objectives. J. Anim. Sci. 30, 860-865.
Hossein-Zadeh N.G. and Ardalan M. (2008). Comparison of different models for the estimation of genetic parameters of body weight traits in Moghani sheep. Agric. Food Sci. 19, 207-213.
Kargar Borzi N., Ayatollahi Mehrgardi A., Asadi Fozi M. and Vatankhah M. (2017). Determining the appropriate selection index for Rayeni Cashmere goat under pasture-based production system. Anim. Prod. Sci. 58, 1595-1602.
Kosgey I.S., Van Arendonk K.G. and Baker R.L. (2003). Economic values for traits of meat sheep in medium to high production potential areas of the tropics. Small Rumin. Res. 50, 187-202.
Lôbo R.N.B., Pereira I.D.C., Facó O. and McManus C. (2011). Economic values for production traits of Morada Nova meat sheep in a pasture based production system in semi-arid Brazil. Small Rumin. Res. 96, 93-100.
Meyer K. (2006). WOMBAT-A Program for Mixed Model Analyses by Restricted Maximum Likelihood. User Notes. Animal Genetics and Breeding Unit, Armidale.
Mohammadabadi M., Esfandyarpoor E. and Mousapour A. (2017). Using inter simple sequence repeat multi-Loci Markers for studying genetic diversity in Kermani sheep. J. Res. Dev. 5, 154-158.
Mokhtari M.S., Rashidi A. and Mohammadi Y. (2008). Estimation of genetic parameters for postweaning traits of Kermani sheep. Small Rumin. Res. 80, 22-27.
Mokhtari M.S., Rashidi A. and Esmailizadeh A.K. (2010). Estimates of phenotypic and genetic parameters for reproductive traits in Kermani sheep. Small Rumin. Res. 88, 27-31.
Okeno T.O., Magothe T.M., Kahi A.K. and Peters K.J. (2012). Application of risk-rated profit model functions in estimation of economic values for indigenous chicken breeding, Trop. Anim. Health Prod. 44, 1279-1287.
Parish J.A., Rhinehart J.D. and Smith T. (2011). Expected Progeny Differences and Selection Indices for Beef Cattle Selection. Extension Service of Mississippi State University Publication, Mississippi, United States.
Paswan C., Gowane G.R., Prakash V., Prince L.L.L. and Kumar A. (2016). Construction of new selection criteria for Malpura sheep. Indian J. Anim. Sci. 86, 61-64.
Quinton M. and Smith C. (1995). Comparison of evaluation selection system for maximizing genetic response at the same level of inbreeding. J. Anim. Sci. 73, 2208-2212.
Rashidi A., Mokhtari M.S., Safi Jahanshahi A. and Mohammad Abadi M.R. (2008). Genetic parameter estimates of pre weaning growth traits in Kermani sheep. Small Rumin. Res. 74, 165-171.
Seidavi A.R., Bizhannia A.R., Ghanipoor M. and Qotbi A.A.A. (2008). Investigation on improvement possibility of resistance, production and reproduction traits in 3P, 2P and P generations in three Japanese pure lines of silkworm Bombyx mori, using individual selection in 3P generation. Asian J. Anim. Vet. Adv. 3, 443-447.
Snyman M.A., Olivier J.J. and Olivier W.J. (1996). Variance components and genetic parameters for body weight and fleece traits of Merino sheep in an arid environment. South African J. Anim. Sci. 26, 11-14.
Weigel K.A.) 2001(. Controlling inbreeding in modern breeding programs. J. Dairy Sci. 84, 177-184.
Volume 10, Issue 4
December 2020
Pages 631-637
  • Receive Date: 18 December 2018
  • Revise Date: 12 March 2019
  • Accept Date: 16 March 2019
  • First Publish Date: 01 December 2020