Livestock Research for Rural Development 23 (7) 2011 Notes to Authors LRRD Newsletter

Citation of this paper

Buffalo milk production in Brazil and Colombia: Genotype by environment interaction

Naudin Hurtado-Lugo*, Mario Cerón-Muñoz*, Raul Aspilcueta-Borquis, Roberta Sesana, Lucia Galvão de Albuquerque and Humberto Tonhati

Faculty of Agrarian and Veterinary Sciences (FCAV), University Estadual Paulista, Jaboticabal (UNESP), 14884900, SP, Brazil
* GaMMA, Faculty of Agricultural Sciences, University of Antioquia, Medellín, Colombia
Corresponding author: mceronm@hotmail.es

Abstract

The objective of this study was to determine whether there is a genotype by environment interaction (GxE) for dairy buffaloes in Brazil and Colombia. The (co)variance components were estimated by using a bi-trait repeatability animal model with the REML method. Each trait consisted in the milk yield obtained in both countries. Contemporary group (herd, year and season of parity) and age at parity (linear and quadratic covariate) fixed effects, along with the additive genetic, permanent environment, and the residual random effects were included in the model.

 

Genetic, permanent environmental and residual variance and heritabilities were different for both countries. The genetic correlations for milk yield between Brazil and Colombia were low (between 0.10 and 0.13), indicating a GxE interaction between both countries. Knowing that this interaction influences the genetic progress of buffalo populations in Brazil and Colombia, we recommend choosing sires tested in the country they will be used, along with conducting joint genetic evaluations that consider GxE interaction effects.

Key words: Genetic correlation, Genetic progress, Genotype by environment interaction, Milk yield


Introduction

The  genotype by environment interaction effect has been previously studied, aiming to meet the need of evaluating the genotypes in the environments where they are located (Banos and Smith 1991; Stanton et al 1991; Cienfuegos-Rivas et al 1999; Costa et al 2000; Fikse et al 2003).  This effect reduces the selection response, which negatively impacts milk yield (MY) and thus diminishes the economic return (Stanton et al 1991; Cienfuegos-Rivas et al 1999; Costa et al 2000; Calus et al 2004; Valencia et al 2008).

 

Genetic evaluation of buffaloes in South American countries has been already implemented, especially in Brazil and Colombia. There is a genotype and environmental variability of MY and also exchange of genetic material in both populations (Aspilcueta-Borquis et al 2010a; 2010b; Gutierrez-Valencia et al 2006; Hurtado-Lugo et al 2006).

 

In order to conduct genetic evaluations between Brazil and Colombia, it is necessary to understand its viability and importance. The objective of this study was to investigate the existence of a GxE interaction for milk production in buffaloes for these two countries.

 

 

Materials and Methods

 

Data covering 25877 lactations of 13857 buffalo cows from 12 herds in Sao Paulo and Ceará states (Brazil), and five herds in Cordoba and Magdalena provinces (Colombia) were used. Data was provided by the Brazilian Breeder´s Association of Buffaloes (made available through the Universidade Estatual Paulista, Campus de Jaboticabal-Brazil), and the Colombian Buffalo Breeders Association (available through the Universidad de Antioquia, Colombia) (Table 1). The animals of Brazilian herds were raised on Brachiaria spp and Panicum sp pastures, and supplemented with 1 kg concentrate per 3 kg of milk yielded. Animals were also supplemented with forage, especially sugar cane, and mineral salt, offered ad libitum during the dry season (April to October). The AM/PM milking test controls in were conducted keeping the calves close to the dams in order to stimulate milk production.  Some herds used manual milking system.


Table 1. Distribution of buffalo cattle records for milk yield in Brazil and Colombia.  

 

Number of records

Milk yield records

25,877

Inbred animals

3,530

(with average inbreeding coefficient)

(5.16%)

Sires

806

 

Brazil

Colombia

Herds

12

5

Female buffaloes

5,710

8,147

Lactations

4,888

9,531

Buffaloes with lactations

1,806

3,291

Lactations of daughters of sires common to both countries

1,683

Sires with daughters

272

519

Sires with daughters in both countries

27


The animals of Colombian herds were mostly raised on improved and native grasses (Brachiaria spp,  Panicum sp, among others), and to a lesser extent on forages. Animals were supplemented with hay from native grasses (e.g., Dichanthium annulatum), cropped in the same farms during the dry season (December to March). Milking was manual, once a day during the mornings, keeping the calf close to the cow. To take the milk samples, the calf was kept close to the cow, simulating normal milking conditions.

 

Milk yield was adjusted to 240 (MY240) and 270 (MY270) days in milk (DIM) using the ICAR method (ICAR, 2002), with adjustment factors suggested by Gutierrez-Valencia et al (2006) and Tonhati et al (2004) for Colombia and Brazil, respectively.

 

The bi-trait analysis was conducted using an animal model with repeated measurements, including the derivative-free restricted maximum likelihood method. The (co)variance components were estimated using the MTDFREML software (Boldman et al 1995). The fixed effects of contemporary group (herd and year of calving), the calving season (Colombia: January-April, May-July, August-October and November-December; and Brazil: April- September, October-March) and age of calving (linear and quadratic effects), as well as the random effects of animal, permanent  environment, and residual, were included. The bi-trait model for 240MY between the 2 countries can be represented in a matrix notation as:

 

 


Results and Discussion

The 240MY average was 871±323 and 1638±652 kg for Colombia and Brazil, respectively. Regarding 270MY, the average was 917±356 and 1722±703 kg, for Colombia and Brazil, respectively (Table 2). Other studies in dairy buffaloes in the same countries reported that 240MY and 270MY ranged from 1025 to 1712 kg (Tonhati et al 2004; Hurtado-Lugo et al 2006; Gutierrez-Valencia et al 2006). The average milk yield values reported in this study are in agreement with those for dairy buffalo herds in Asia, India and Europe, which ranged from 796 to 2545 kg (Gogoi et al 1985; Singh and Yadav 1987; Tien and Tripathi 1991; Mathur and Mathur 1992; Rosati and Van Vleck 2002). Variance components, heritabilities, and repeatabilities were within the values reported by others researchers (Aspilcueta-Borquis et al 2010a; 2010b; Rosati and Van Vleck 2002; Tonhati et al 2004; Hurtado-Lugo et al 2006).


Table 2. Estimates of (co)variance, heritabilities, repeatabilities and genetic correlations for milk yield at 240 and 270 DIM for dairy buffaloes in Colombia and Brazil.

 

Milk yield at 240 DIM

 

Milk yield at 270 DIM

 

Colombia

Brazil

 

Colombia

Brazil

Milk yield average (kg)

871

1,638

 

917

1,722

Standard deviation

323

652

 

356

703

Coefficient of variation (%)

37

40

 

39

41

Additive genetic variance (kg2)

13,064

43,731

 

15,706

50,650

Environment permanent variance (kg2)

11,818

53,366

 

13,180

64,287

Residual variance (Kg2)

38,839

117,144

 

44,663

143,107

Heritability

0.21

0.20

 

0.21

0.20

Repeatability

0.39

0.45

 

0.39

0.45

Additive genetic covariance

2,332

 

3,651

Genetic correlation

0.10

 

0.13


Means and genetic, environment and residual variance of milk yield were higher in Brazil, but heritabilities and repeatabilities were similar (Table 2). Several researchers have indicated that variance heterogeneity is due to differences between production systems, environmental conditions specific to each region, herd size and management, and number of daughters per sire (Calus et al 2004; Konig et al 2005; Valencia et al 2008; Stanton et al 1991). In this study, the low number of sires common to both countries may be increasing the heterogeneity of variances. Consequently, this factor should be considered for future joint genetic evaluations.

 

Genetic correlations between both countries varied from 0.10 to 0.13 (Table 2), indicating that a genotype-environment interaction exists, because there were differences in genetic values and changes in rank of sires (Stanton et al 1991; Hill et al 1983; Calus and Veerkamp 2003), as shown in Figure 1. Because of this, changes were observed in the genetic values ​​and the animals tested were reclassified, implying that the sires evaluated as superior in Brazil may not improve the productive performance of their descendants in Colombia (Table 3).


Figure 1. Breeding values for milk yield (kg) at 240 (left) and 270 (right) days for several buffalo sires in Colombia and Brazil


 Table 3. Genetic values (BV) for milk yield at 240 (PL240) and 270 (PL240) days of lactation, re-classification (C) and reliabilities (R) for buffaloes in Colombia and Brasil, obtained trough bi-carácter analysis. 

PL240

Buffaloes siries

Colombia

Brazil

BV

R

C

BV

R

C

14004

182

0,89

40

391,32

0,91

1

14610

59,34

0,44

3342

341,11

0,7

4

14180

191,73

0,54

25

282,91

0,8

36

14262

115,56

0,87

497

193,81

0,83

261

4509

181,25

0,78

41

152,05

0,89

487

4494

105,75

0,67

760

104,81

0,58

935

20942

-4,96

0,08

9204

101,97

0,39

995

23236

115,88

0,78

494

86,24

0,77

1455

12883

143,72

0,76

168

67,59

0,95

1771

11027

136,36

0,48

220

33,69

0,49

2976

21297

143,49

0,45

173

17,79

0,69

3650

23238

75,25

0,36

2150

8,87

0,43

4153

26309

126,05

0,69

329

-1,52

0,7

7124

23239

106,72

0,4

733

-15,28

0,44

8281

23226

144,79

0,69

159

-45,64

0,33

9071

14439

230,17

0,75

6

-47,56

0,66

9105

PL270

Buffaloes siries

Colombia

Brasil

BV

R

C

BV

R

C

14004

203,26

0,89

31

420,93

0,91

1

14610

62,97

0,45

3461

362,21

0,7

5

14180

190,98

0,55

51

308,18

0,81

34

14262

118,07

0,88

650

180,86

0,84

374

4509

172,61

0,79

86

150,97

0,9

551

4494

116,97

0,68

671

109

0,59

999

20942

2,81

0,04

7738

133,23

0,39

697

23236

118,85

0,79

632

90,19

0,78

1491

12883

157,55

0,77

148

64,98

0,95

2004

11027

151,18

0,48

184

24,76

0,5

3477

21297

164,24

0,46

109

35,68

0,7

2963

23238

98,52

0,36

1343

7,38

0,44

4301

26309

136,35

0,69

317

-20,87

0,7

8335

23239

124,96

0,4

501

-10,94

0,45

8062

23226

185,14

0,7

58

-55,55

0,33

9119

14439

209,87

0,77

26

-58,11

0,67

9165

             

Conclusions


Implications


Acknowledgments

The authors wish to acknowledge the logistic support from the Colombian and Brazilian buffalo breeders, and financial support of the Ministerio de Agricultura y Desarrollo Rural de Colombia, Asociación Colombiana de Criadores de Búfalos, Federación Colombiana de Ganaderos (Fedegan), Ministerio Brasileiro da Agricultura, Pesca e Abastecimento, Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and to Fundação de Apoio à Pesquisa do Estado de São Paulo (FAPESP–Process Nº 2006/59270-3).


References

Aspilcueta-Borquis R R, Araujo Neto F, Baldi F, Bignardi A B, Albuquerque L G y  Tonhati H 2010a Genetic parameters for buffalo milk yield and milk quality traits using Bayesian inference.  Journal of Dairy Science 93:2195-2201.http://www.journalofdairyscience.org/article/S0022-0302(10)00210-9/fulltext.Accessed Apr. 7, 2010.

 

Aspilcueta-Borquis RR, Sesana R, Berrocal M, Seno O, Bignardi AB, El Faro L, Albuquerque L, Camargo G and Tonhati H 2010b Genetic parameters for milk, fat and protein yields in Murrah buffaloes. Genetics Molecular Biology 33(1):71-77. http://www.scielo.br/pdf/gmb/v33n1/2008-338.pdf Accessed March. 10, 2010.

 

Banos G and Smith C 1991 Selecting bulls across countries to maximize genetic improvement in dairy cattle. Journal of Animal Breeding and Genetics 108:174-181.

 

Boldman K, Kriese L D and Van Vleck C 1995 A manual for use of MTDFREML. Lincoln: Department of Agriculture / Agricultural Research Service, 1995. 120, p.

 

Calus M, Bijma P and Veerkamp R 2004 Effects of data structure on the estimation of covariance functions to describe genotype by environment interactions in a reaction norm model. Genetics Selection Evolution 36: 489-507. http://www.gse-journal.org/index.php?option=com_toc&url=/articles/gse/abs/2004/05/contents/contents.html. Accessed July. 10, 2007

 

Calus M and Veerkamp R 2003 Estimation of environmental sensitivity of genetic merit for milk production traits using a random regression model. Journal of Dairy Science 86: 3756-3764. http://www.journalofdairyscience.org/article/S0022-0302(03)73982-4/abstract Accessed June. 11, 2007.

 

Cienfuegos-Rivas E, Oltenacu P, Blake R, Schwager S, Castillo-Juarez H and Ruíz H 1999 Interaction between milk yield of Holstein cows in México and the United States. Journal of Dairy Science 82: 2218-2223. http://www.journalofdairyscience.org/article/S0022-0302(99)75468-8/abstract. Accessed Feb 21, 2007..

 

Costa C, Blake R, Pollak E, Oltenacu A, Quaas L and Searle R 2000 Genetic analysis of Holstein cattle populations in Brazil and the United States. Journal of Dairy Science 83: 2963-2974. http://www.journalofdairyscience.org/article/S0022-0302(00)75196-4/abstract. Accessed Agus 03, 2007.

 

Fikse W, Rekaya R and Weigel A 2003 Genotype x Environment interaction for milk production in Guernsey Cattle. Journal of Dairy Science 86:1821-1827. http://www.journalofdairyscience.org/article/S0022-0302(03)73768-0/abstract Accessed Jan. 11, 2007.

 

Gogoi P, Johar K and Singh A 1985 Genetic analysis of milk yield in Murrah buffaloes. Indian Veterinary  Journal 62:970-975.

 

Gutierrez-Valencia A, Cerón-Muñoz M, Hurtado-Lugo N 2006 Estimativas de factores de  corrección para duración de la lactancia, edad y época de parto en búfalas de la Costa Atlántica Colombiana. Livestock Research Rural Development 18: (4). http://www.lrrd.org/lrrd18/4/guti18050.htm  Accessed Feb. 15, 2009.

 

Hill W, Edwards M, Thompson R 1983 Heritability of milk yield and composition at different levels and variability of production. Animal Production 36: 59-68.

 

Hurtado-Lugo N, Cerón-Muñoz M, Gutierrez-Valencia A 2006 Estimación de parámetros genéticos para la producción de leche en el día del control en búfalos de la Costa Atlántica Colombiana. Livestock Research Rural Development. 18 (3). http://www.lrrd.org/lrrd18/3/hurt18039.htm  Accessed Feb. 19, 2009.

 

ICAR 2002 Guidelines for buffalo milk recording for low to medium and medium to high input production systems. http://www.icar.org/pages/recording_guidelines.htm  Accessed Jan. 20, 2009.

 

Konig S, Dietl G, Raeder I, Swalve H 2005 Genetic relationships for dairy performance between large-scale and small-scale farm conditions. Journal of Dairy Science 88: 4087-4096. http://www.journalofdairyscience.org/article/S0022-0302(05)73093-9/abstract. Accessed May. 19, 2007.

 

Mathur AK and Mathur BS 1992 Murrah buffaloes in their breeding tracts. Indian J Animal Science 62: 961-967.

 

Rosati A and Van Vleck L 2002 Estimation of genetic parameters for milk, fat, protein and mozzarella cheese production in the Italian river buffalo population. Livestock Production Science 74: 82-85.

 

Singh C and Yadav C 1987 Effect of genetics and non-genetics factors on life time production traits in Indian buffaloes. Indian J Anim Sci 57:555-558.

 

Stanton T, Blake R, Quaas D, Van Vleck L and Carabaño M 1991 Genotype by environment interaction for Holstein milk yield in Colombia, Mexico, and Puerto Rico. Journal of Dairy Science 74:1700-1714. http://www.journalofdairyscience.org/article/S0022-0302(91)78333-1/abstract. Accessed May. 19, 2007.
 

Tien NQ and Tripathi VN 1991 Effect of age and weight at first calving on first lactation production in Murrah buffaloes.  Buffalo Bulletin 10:3-7.

 

Tonhati H, Muñoz F, Duarte J, Reichert R, Oliveira R, Lima A 2004 Estimates of correction factors for lactation length and genetic parameters for milk yield in buffaloes. Arquivo Brasileiro de Medicina Veterinária e Zootecnia 251-257.

 

Valencia M, Montaldo H, Ruíz F 2008 Interaction between genotype and geographic region for milk production in Mexican Holstein cattle. Archivo de Zootecnia 57:457-463.



Received 14 March 2011; Accepted 4 May 2011; Published 1 July 2011

Go to top