Livestock Research for Rural Development 30 (1) 2018 Guide for preparation of papers LRRD Newsletter

Citation of this paper

Effect of inbreeding on traits of economic importance in Kenyan Sahiwal cattle

M Musingi1, T K Muasya, E D Ilatsia2 and A K Kahi

Animal Breeding and Genomics Group, Department of Animal Sciences, Egerton University, P O Box 536, 20115 Egerton, Kenya
muasyakt@yahoo.com
1 Department of Biological Sciences, Egerton University, P O Box 536, 20115 Egerton, Kenya
2 Kenya Agricultural and Livestock Research Organisation, Dairy Research Institute, P O Box 25, 20117 Naivasha, Kenya

Abstract

Pedigree and performance data for the Sahiwal cattle from the National Sahiwal Stud (NSS) in Kenya were used to evaluate the effect of inbreeding on first lactation milk yield (LMY), age at first calving (AFC), calving interval (CI), and lactation length (LL). Effects of inbreeding on the traits were determined by fitting four regression models (linear, quadratic, exponential and Michaelis-Menten) to the errors generated by the animal model. The linear, exponential and Michaelis-Menten models were significant for all the studied traits while the quadratic model was only significant for calving interval. Inbreeding had a positive effect on calving interval, age at first calving, and lactation length, shortening calving both interval and age at first calving and increasing lactation length. The relation between inbreeding and depression of traits was not linear, with greater depression after 15% inbreeding. Genetic evaluation of the Kenyan Sahiwal should account for inbreeding. The results of the current study indicate the need to consider effect of inbreeding on traits of economic importance for the Sahiwal cattle breed when carrying out genetic evaluations. The mating plan for the NSS should be designed so as to control future rates of inbreeding while achieving genetic gain.

Key words: genetic gain,inbreeding depression


Introduction

Increase in inbreeding leads to reduced genetic variability by reducing heterozygosity over many loci (Falconer and Mackay 1996), increased risk of breeding programmes due to variance of genetic gains (Meuwissen, 1991) and increased risk of emergence of lethal recessive homozygous alleles (König and Simianer 2006). An immediate concern to dairy farmers is the reduction of performance of inbred animals, referred to as inbreeding depression due to inbreeding (Falconer and Mackay 1996).

For the Kenya Sahiwal population, Muasya et al (2011) reported a gradual increase in inbred animals from 0 in 1967 to 80% in 2008 with a corresponding increase in mean individual inbreeding coefficient from 0 to about 2.5% in 2008. The proportion of inbred animals increased rapidly from 1% to about 98% in the most recent complete generation. As the proportion of inbred animals increase, reduction of inbreeding through pairing of mates that are less related than the average in the population becomes difficult (Thompson et al 2000). In the Kenyan Sahiwal cattle, inbreeding level increased from 1.2 to 2% as the proportion of inbred individuals increased (Muasya et al 2011). Other studies have reported various estimates of inbreeding for zebu the Sahiwal breed. In Nicaragua an inbreeding level of 13.0% has been reported for Creole cattle breed (Corrales et al 2010). Among the Nellore, Guzerat and Gyr zebu cattle breeds in Brazil, inbreeding levels ranging from 1.8 to 2.8% have been reported (Filho et al 2010; Faria et al2009).

Inbreeding depression decreases cow survival, reproductive performance and milk production and increases rate of disposal or loss of replacement heifers before first calving, age at puberty through reduced growth (du Toit et al 2012). Every 1% increase in inbreeding leads to a 10 kg of milk to 26kg decline in milk production per lactation (Mostert, 2011) and decrease of about 13 days in productive life (CDN, 2008; Smith et al 1998). du Toit et al (2012) reported a similar effect of inbreeding on functional herd life in Jersey cattle. Cows with high inbreeding level have also been reported to have a high risk of being culled (Rokouei et al 2010).

Inbreeding depression depends not only on actual level of inbreeding but also on the rate of inbreeding such that animals with the same level of inbreeding may have different inbreeding depression effects depending on the completeness of their respective pedigrees (González-Recio et al 2007; Gutiérrez et al 2009). The quality of pedigree can be accounted for by estimating the rate of inbreeding (Gonzalez-Recio et al 2007) which indicates the increment in inbreeding regardless of number of known generations in an individual’s pedigree (Gonzalez-Recio et al 2007). The objective of the study was therefore to evaluate the effect of inbreeding on lactation milk yield, lactation length and fertility traits in the Kenyan Sahiwal cattle population.


Materials and methods

Description of the study sites and data collection

Pedigree data for the Sahiwal cattle in Kenya were obtained from the Sahiwal National Stud at Kenya Agricultural and Livestock Research organization (KALRO), Naivasha. A database was created by systematically entering each animal, sire and dam, date of birth and sex into a database. The data were checked for consistency to ensure that all animals were ordered sequentially based on date of birth such that no progeny was born before any of its parents and that there were no cyclic pedigrees (i.e. no animal appeared as both male or female) and no animal appeared as both sire and dam using the Animal Breeder’s Tool Kit (ABTK, Golden et al 1992). Data on lactation milk yield for the first lactation (LMY), lactation length (LL), calving interval (CI) and age at first calving (AFC) were added onto the pedigree database of each animal. Age at first calving was derived as the interval in days from date of birth to first calving. Calving interval was derived as the interval between consecutive calvings for the first three lactations. Means, standard deviations (SD) and coefficients of variation (CV %) for LMY, LL, AFC CI are given in Table 1.

The large CV% for LMY reflects the great variation between individual cows in these traits. Similar results were reported for Egyptian buffalos (Khatab et al 2007).

Table 1. Descriptive statistics of lactation milk yield (LMY), age at
first calving (AFC), calving interval (CI) and lactation length (LL)

Trait

No. of records

Mean

SD

LMY, kg

1841

1287.9

509.6

AFC, months

1767

44.1

7.2

CI, days

1841

474.5

125.5

LL, days

1837

287.1

50.7

Assessment of level of inbreeding

The inbreeding coefficient (Fi) for each individual in the pedigree was calculated as the probability that two alleles are identical by descent according to the method of VanRaden (1992). This method makes it possible to compute the inbreeding coefficients per generation, and assumes that founders are inbred or related, an assumption, which is important when pedigrees are heterogeneous. The inbreeding coefficient of an individual with unknown origin (founder) is equal to half the average genetic relationship between genetic groups of its phantom parents.

Estimation of (co)variances and estimated breeding values

Estimates of variance and estimated breeding values were obtained by performing univariate analyses on first lactation milk yield (LMY), lactation length (LL), age at first calving (AFC) and calving interval (CI) using the following animal model:

Y= Xb + Za + e

where y is a vector of observations for LMY, LL, AFC and CI in each lactation, b and a and e are vector of fixed effects, random animal effects, and random residuals, respectively and X and Z are the incidence matrices relating fixed and random animal effects to observations, respectively. The random effects will be assumed to be normally distributed with a mean=0 and variance as follows: a ~ N (0, A I σ 2a ) and a ~ N (0, I σ 2e ) whereσ 2a and σ 2e are animal and residual variances, respectively, and A and I are the numerator relationship and identity matrices, respectively. The models for LMY, CI and LL included fixed effects of year-season of calving and age at calving. A similar model was used for AFC but with year-season of birth instead of year season of calving. Variance covariance components for the traits under study were estimated using MTDFREML software package (Boldman et al 1995).

Analysis of pedigree, inbreeding coefficient (F) and equivalent complete generations were estimated using ENDOG version 4.5 computer programme (Gutiérrez and Goyache 2005). Equivalent complete generations were estimated as:

where n is the number of generations separating each individual from each known ancestor (Maignel et al 1996). Inbreeding coefficient (F) was computed after Meuwissen and Luo (1992). Rate of inbreeding (∆F) was calculated for each generation as:

where Ft and Ft-1 are the average inbreeding at tth generation. Only animals with ECG of at least included in the analysis of F (2367 animals).

Estimation of inbreeding depression

A vector of errors for each trait were generated from univariate analysis using WOMBAT (Meyer 2007) as follows;

ei = Xiβ + Zia

The effect of inbreeding on the traits studied was determined by fitting the following four regression models:

around the errors. The model of the analyses was ei=Φ(Fi)+ɛ, where Φ(Fi) is the regression function on inbreeding efficient of an animal (Fx) and ɛ represents the deviations of the errors from the predicted errors in the regression function.


Results and discussion

Summary of number of animals, mean inbreeding levels by groups of birth year and number of animals in each inbreeding class are shown in Table 2. In the 1965-1974 birth year group 79% of the animals were not inbred, while in the last year group only 11% were non-inbred. Mean inbreeding level across the year of birth groups increased from 0.9 to 1.6%. A similar trend was reported in Jersey cattle in South Africa (du Toit et al 2012).

The mean inbreeding level and annual rate of inbreeding for the population studied were 0.7 and 0.33%, respectively. Higher average inbreeding levels of 15% in Guzerat cattle (Panetto et al 2010), 1.75 to 2.28% in Brazilian zebu cattle (Faria et al 2009) and 1.73% in the Nellore breed (Brito et al 2013) have been reported. Rates of increase in inbreeding lower than that reported in the current study ranging from 0.04 to 0.08% were reported in dairy cattle breeds in Canada (CDN 2010).

Table 2. Number of animals, mean inbreeding by birth year group and
distribution of animals in each inbreeding class for cows

Birth year

No.
animals

Inbred
(%)

Mean inbreeding
level (%)

1965-1974

732

11.1

0.9

1975-1984

733

12.1

1.0

1985-1994

480

50.8

1.5

1994-2008

422

78.8

1.6

The general increase in inbreeding among the cows in the current study could be linked to intense use of a small number of superior sires and more complete pedigree information in more recent years (Muasya et al 2011).

Table 3. Least square means (standard errors) of the effect of inbreeding coefficients on lactation milk yield (LMY), age at first calving (AFC), calving interval (CI) and lactation length (LL) in Sahiwal cows.

Inbreeding coefficient

No. of animals

LMY, kg

LL, days

AFC, days

CI, days

Zero

2015

1227.9 (15.3)a

287.4(1.4) a

1387.8(5.3)a

474.2 a

Non-zero

352

970.4 (128.9)b

283.8(11.5) a

1388.3(43.6) a

538.4b

For LMY, non-inbred animals produced more milk yield and longer CI (P<0.05) compared to inbreds. For the other traits, non-inbred animals performed better than inbred animals, though the means were not different (P>0.05). The depression in performance or inbreeding depression is attributed to inbreeding depression caused by increased homozygosity of individuals. Increased homozygosity lowers fitness through increase in the frequency of homozygous recessive detrimental mutations, and increased homozygosity for alleles at loci with heterozygote advantage (Griffiths et al 1999; Charlesworth and Willis 2009). There is emerging evidence that lowly heritable traits such as fertility (e.g. CI, AFC) exhibit more depression due to inbreeding because of their low genetic variation (Kristensen et al 2005), though results from the current study showed significant depression in LMY only (Table 3).

Effect of inbreeding on the traits studied is shown in Table 4. The linear, exponential and Michaelis-Menten models were significant (P<0.05) for all traits studied (Table 4). The quadratic model was significant (P<0.05) for CI only. Carrillo and Siewerdt (2010) and Mahlado et al (2013) reported a similar pattern, where the linear, exponential and Michaelis-Menten models were significant for all traits studied, whereas the quadratic model was not significant for any traits studied.

Table 4. Estimates of the inbreeding regression coefficients for first lactation milk yield (LMY), lactation length (LL), calving interval (CI) and age at first calving (AFC) of Kenyan Sahiwal cattle

Models

Linear

Quadratic

Exponential

Michaelis-Menten

LMY (kg)

-16.7*

-0.42ns

0.74**

-0.198***

LL (days)

15.17*

0.33ns

1.02*

-0.48***

CI (days)

-13.47*

-0.09*

1.45**

0.23***

AFC (mo)

-0.18*

-0.05ns

0.95*

0.06***

ns=not significant, *=P<0.05; **=P<0.01; ***P<0.001

The residual variances for the three models were similar (Table 5) implying equivalence in goodness of fit. Similar results were reported by Carrillo and Siewerdt (2010) and Mahlado et al (2013), meaning that any of the three models would be adequate in estimating inbreeding depression. Of these models, the linear model, which is to fit and the fact that its parameters have direct interpretation, would suffice. However, its limitations are its lack of flexibility and its inability to cope with departure from linearity at higher inbreeding levels (Figure 1).

Table 5. Error means square from linear, exponential and Michaelis-Menten for first lactation milk yield (LMY), lactation length (LL), age at first calving (AFC) and calving interval (CI) of Kenyan Sahiwal cattle

Linear

Quadratic

Exponential

Michaelis-Menten

First lactation Milk yield

120575

120616

120625

120313

Lactation length

594.9

597.5

595.5

579.0

Calving interval

7459.8

7438.0

7434.0

7489.0

The increase in inbreeding led to a decrease in first lactation milk yield, calving interval and age at first calving (Figure 1) implying a depression in milk yield but an improvement in CI and AFC.

Figure 1. Effects of inbreeding depression on LMY, CI, LL and AFC for Kenyan Sahiwal cattle

Similar findings have been reported in previous studies. Mahlado et al (2013) found that the exponential and the Michaelis-Menten models were significant for all traits studied. However, unlike in the current study, the exponential model was not significant for any of the traits studied. In the current study, for inbreeding levels up to 0.15, inbreeding was approximately linear. Beyond that level of inbreeding, there seemed to be a higher inbreeding depression in all traits (Figure 1). Similar results have been reported by Mahlado et al (2013) and Santana Junior et al (2011). In the study by Santana Junior et al, animals with an average F=25% produced less milk by 107 kg compared to non inbreds. In the current study, inbred animals (F=15%) had a depression in milk production of 258 kg in the first lactation. In Mahlado et al (2013) study, at an average F=0.25, animals produced 50.4 kg less. Others authors that have reported a decline in milk yield due to increase in inbreeding include Panetto et al (2010) who found a reduction of 15.25 kg in milk yield in Guzerat cattle when half sibs were mated.

In populations with inbreeding levels of 12.5%, high inbreeding depression for milk yield ranging from 345 to 480 kg in Holstein cows translating to 27.6 kg to 38.4 kg per 1% increase in inbreeding have been reported (Smith et al 1998; Thompson et al 2000). Other studies have found a non-linear relationship between inbreeding and milk yield resulting in a depression of 380 kg in 305 day milk yield among highly inbred cows (F=12.5%) (Gulisija et al 2007). For cows with a similar inbreeding level, Maiwashe et al (2006) reported decline of 170kg in 305 milk yield. Panetto et al (2010) reported a positive effect of inbreeding on daily milk yield, contrary to that reported in the current study. Among Buffaloes in Brazil, Santana Junior et al (2011) reported a decrease of 4.3 kg while Malhado et al (2013) found a value of 14.8 kg for 1% increase in inbreeding; whereas Mirhabibi et al (2007) reported a decrease of 18.2 kg.

Lactation length deteriorated by 0.3 days per 1% increase in inbreeding, though the effects were not significant (P<0.05). To the contrary, inbred animals had slightly better fertility compared to non-inbreds (Table 5). Similar to the current study, Malhado et al (2013) found a favourable relationship between inbreeding and AFC among Brazilian Buffaloes. In this study, a 1% increase in inbreeding led to a shortening of AFC by 0.76 days. Thompson et al (2000) reported that low to moderate levels of inbreeding (F<0.07) reduced AFC by 3 to 5 days compared to non-inbred animals. Contrary to the positive effects of inbreeding on CI and AFC reported in the current study, Panetto et al (2010) found that these traits increased by up to 14 day and 11 days for every 1% increase in inbreeding in a Guzerat dairy herd. Other studies reported increase in CI and AFC. Hudson and Van Vleck (1984) reported an expected increase in CI by 2 day whereas AFC was expected to increase by 0.4 days for every 1% increase in inbreeding. Among Irish Holstein-Friesian cows, McParland et al (2007) reported that mating between half-sibs (F=12.5%) would result in an increase in AFC and CI of 2.5 days and 8.8 days, or 0.2 and 0.7 days per 1% increase in inbreeding, respectively. A smaller effect of inbreeding depression of 0.263 days on CI was reported for Alentejana cattle (Carolino and Gama 2008). Other studies that have reported a decline in fertility due to inbreeding include González-Recio et al (2007), Malhado et al (2013). In Iranian Holstein cattle, Rokouei et al (2010) reported a non-significant effect of inbreeding on CI in first and second lactations.


Conclusion


Acknowledgements

The authors are grateful to the Kenya Agriculture and Livestock research organization (KALRO) for provision of data and Egerton University, Kenya for financial support and provision of computing facilities.


References

Boldman K G, Kriese L A, Van Vleck L D, Van Tassel C P and Kachman 1995 A Manual for Use of MTDFREML: A Set of Programs to Obtain Estimates of Variances and Covariances [Draft]. ARS, USDA, Washington, DC.

Brito F V, Sargolzaei M, Braccini Neto J, Cobuci J A, Pimentel C M, Barcellos J and Schenkel F S 2013 In-depth pedigree analysis in a large Brazilian Nellore herd.Genetics and Molecular Research 12 (4): 5758-5765: https://www.geneticsmr.com/

Canadian Dairy Network (CDN) 2008 Subject: Quantifying inbreeding depression. https://www.cdn.ca. Accessed March, 2017.

Canadian Dairy Network (CDN) 2010 Subject: Inbreeding Update. https://www.cdn.ca. Accessed: Aug. 2017

Carolino N and Gama L T 2008 Inbreeding depression on beef cattle traits: estimate, linearity of effects and heterogeneity among sire-families. Genetics Selection and Evolution 40: 511– 527: https://gsejournal.biomedcentral.com/

Carrillo J A and Siewerdt F 2010 Consequences of long-term inbreeding accumulation on preweaning traits in a closed nucleus Angus herd. Journal of Animal Science 88: 87–95: www.animalsciencepublications.org/publications/jas

Charleswoth D and Willis J H 2009 The genetics of inbreeding depression. Nature Review Genetics 10: 783–796: www.nature.com/nrg/

Corrales R, Näsholm A, Malmfors B and Philipsson J 2010 Population structure of Reyna Creole cattle in Nicaragua. Tropical Animal Health and Production 42: 1427–1434: https://link.springer.com/

du Toit J, van Wyk J B and Maiwashe A 2012 Assessment of inbreeding depression for functional herd life in the South African Jersey breed based on level and rate of inbreeding. South African Journal of Animal Science 42 (No. 1): 55-62: http://www.sasas.co.za/journals

Falconer D S and Mackay T F C 1996 Introduction to quantitative genetics. 4th edition. Longman Scientific and Technical, Harlow, UK.

Faria F J C, Fihlo A E V, Madalena F E and Josahkian L A 2009 Pedigree analysis in the Brazilian Zebu breeds. Journal of Animal Breeding and Genetics 126: 148 – 153: http://onlinelibrary.wiley.com/advanced/search/results

Filho J C R, Lopes P S, Verneque R S, Torres R A, Teodoro R L and Carneiro P L S 2010 Population structure of Brazilian Gyr dairy cattle. Revista Brasileira de Zootecnia 39(12): 2640 – 2645: http://www.scielo.br/

Golden B L, Snelling W M and Mallinckrodt C H 1992 Animal breeder’s tool kit user’s guide and reference manual. Colorado State University Agricultural Experiment Station Technical Bulletin LTB 92-2: www.ccsi.ca/abtk/

González-Recio O, López de Maturana E and Gutiérrez J P 2007 Inbreeding depression on female fertility and calving ease in Spanish dairy cattle. Journal Dairy Science 90: 5744–5752: http://jds.fass.org/

Griffiths A J F, Jeffrey H M, David T S, Richard C L and William M G 1999 An Introduction to Genetic Analysis. Freeman WH, New York. pp. 726–727.

Gulisija D, Gianola D and Weigel K A 2007 Nonparametric analysis of inbreeding on production in Jersey cows. Journal of Dairy Science 90: 493–500: http://jds.fass.org/

Gutiérrez J P and Goyache F 2005 A note on ENDOG: a computer program for analysing pedigree information. Journal of Animal Breeding and Genetics 122: 172- 176: http://onlinelibrary.wiley.com/

Gutiérrez J P, Cervantes I and Goyache F 2009 Improving the estimation of realized effective population sizes in farm animals. Journal of Animal Breeding and Genetics 126: 327–332: http://onlinelibrary.wiley.com/

Hudson G F S and Van Vleck L D 1984 Effects of inbreeding on milk and fat production, stayability, and calving interval of registered Ayrshire cattle in Northeastern United States. Journal of Dairy Science 67: 171–179: http://jds.fass.org/

Khattab A S and Kawthar A M 2007 Inbreeding and it is effect on some productive and reproductive traits in a herd of Egyptian buffaloes. Italian Journal of Animal Science 6: 275–278: http://www.tandfonline.com/

König S and Simianer H 2006 Approaches to the management of inbreeding and relationship in the German Holstein dairy cattle population. Livestock Science 103: 40 – 53.

Kristensen T N, Sørensen A C, Sorensen D, Pedersen K S, Sørensen J G and Loeschcke V 2005 A test of quantitative genetic theory using Drosophila- effects of inbreeding and rate of inbreeding on heritabilities and variance components. Journal of Evolutionary Biology 18: 763–770: www.ncbi.nlm.nih.gov/

Maignel L, Boichard D and Verrier E 1996 Genetic variability of French dairy breeds estimated from pedigree information. Interbull Bulletin 14: 49-54: https://journal.interbull.org

Maiwashe A, Nephawe K A, Van der Westhuizen R R, Mostert B E and Theron H E 2006 Rate of inbreeding and effective population size in four major South African dairy cattle breeds. South African Journal of Animal Science 36: 50-57: http://www.sasas.co.za/

Malhado C H M, Malhado A C M, Carneiro P L S, Ramos A A, Carrillo J A and Pala A 2013 Inbreeding depression on production and reproduction traits of buffaloes from Brazil. Animal Science Journal 84: 289–295: http://onlinelibrary.wiley.com/

Mcparland S, Kearney J F, Rath M and Berry D P 2007 Inbreeding effects on milk production, calving performance, fertility, and conformation in Irish Holstein-Friesians. Journal of. Dairy Science 90: 4411–4419: http://jds.fass.org/

Meuwissen T H E 1991 Expectation and variance of genetic gain in open and closed nucleus and progeny testing schemes. Animal Production 53: 133-141.

Meuwissen T I and Luo Z 1992 Computing inbreeding coefficients in large populations Genetics, Selection and Evolution 24: 305-313: https://gsejournal.biomedcentral.com/

Meyer K 2007 WOMBAT- A tool for mixed model analyses in quantitative genetics by REML. Journal of Zhejiang University Science B 8: 815-821: http://dideridoo.une.edu.au

Mirhabibi S, Manafiazar G H, Qaravisi S and Mahmoodi B 2007 Inbreeding and its effect on some productive traits in buffaloes of South Iran. Italian Journal of Animal Science 6: 372–374: http://www.tandfonline.com/

Mostert B E 2011 Inbreeding facts of the South African Jersey breed. South African Jersey Journal 60 (1): 40-45.

Muasya T K, Kariuki J N and Muia J M K 2011 Population structure of the Sahiwal breed in Kenya. Livestock Research for Rural Development. Volume 23, Article #186. Retrieved August 11, 2015, from http://www.lrrd.org/lrrd23/9/muas23186.htm.

Panetto J C C, Gutiérrez J P, Ferraz J B S, Cunha D G and Golden B L 2010 Assessment of inbreeding depression in a Guzerat dairy herd: Effects of individual increase in inbreeding coefficients on production and reproduction. Journal of Dairy Science 93: 4902–4912: http://jds.fass.org/

Rokouei M, Torshizi R V, Shahrbabak M M, Sargolzaei M and Sørensen A C 2010 Monitoring inbreeding trends and inbreeding depression for economically important traits of Holstein cattle in Iran. Journal of Dairy Science 93: 3294–3302: http://jds.fass.org/

Santana Júnior M L, Aspiculeta-Bosquis R R, Bignardi A B, Alburquerque L G and Tonhati H 2011 Population structure and effects of inbreeding on Milk yield and quality of Murrah bufalloes. Journal of Dairy Science 94: 5204–52011: http://jds.fass.org/

Smith L A, Cassel B G and Pearson R E 1998 The effects of inbreeding on the lifetime performance of dairy cattle. Journal of. Dairy Science 81: 2729–2737: http://jds.fass.org/

Thompson J R, Everett R W and Hammerschmidt N L 2000 Effects of inbreeding on productive and survival in Holsteins. Journal of. Dairy Science, Volume, 83, 1856–1864: http://jds.fass.org/

VanRaden P M 1992 Accounting for inbreeding and crossbreeding in genetic evaluation of large populations. Journal of Dairy Science 75; 3136 – 3144: http://jds.fass.org/


Received 14 October 2017; Accepted 8 December 2017; Published 1 January 2018

Go to top