Estimation of genetic and phenotypic trends for dairy traits of Gyr and Guzera breeds/Estimativa de tendencias geneticas e fenotipicas para caracteristicas leiteiras de animais das racas Gir e Guzera.
As production efficiency becomes fundamentally important to the production system, it is important to assess the genetic progress achieved throughout the selection process. Knowing the genetic progress of a population allows the breeder to analyze the results obtained with the breeding program (Malhado, Carneiro, Pereira, & Martins Filho, 2008), to know the distance between the selection objectives and the gains obtained over the years and to make the necessary adjustments.
Milk production in general is the most economically important trait in dairy cattle breeding programs. In addition, this trait correlates with other important productive and functional traits in the production system, so the direct selection applied to milk production can produce indirect selection for other correlated traits, such as milk fat and protein percentages (Boligon et al., 2005).
When choosing the traits to be included in a selection program, not only the heritability of the trait and its correlations must be taken into account, but also its economic importance in relation to the general economic performance. Thus, the results depend on the environmental, genetic and value relationships between the traits involved. These relationships, if neglected, may result in losses in productivity and response to selection in each trait, and in the remuneration of these by industry (Lobo, Madalena, & Penna, 2000; Madalena, 2000)
The Gyr breed was a pioneer among the zebu breeds in implanting a genetic breeding program for milk, beginning precisely in 1985. The Guzera breed, in turn, had the genetic breeding program installed 09 years after the Gyr breed. The selection for milk production of these breeds, since then, has shown results of positive genetic trends (Santos, Leite, & Lobo, 1990; Balieiro et al., 2000; Peixoto, Verneque, Teodoro, Penna, & Martinez, 2006). Therefore, it is clear that the genetic improvement tools, if correctly used, are effective for changing the productive performance of herds (Santana Junior et al., 2016).
The objective was to verify the genetic and phenotypic trends of Gyr and Guzera breeds for dairy traits in the state of Rio Grande do Norte.
Material and methods
For the present study, data on the milk, fat and protein (Kg) production traits of the Gyr and Guzera herds, belonging to the Agricultural Research Company of Rio Grande do Norte (EMPARN), and the genetic merit (DEP in Kg) estimated in the genetic evaluations of Gyr and Guzera animals conducted by Embrapa Dairy Cattle from information on the growth performance of pure and crossbred animals in the national breeding programs database. The EMPARN herds are reared at the Rockfeller Experimental Station, with an area of 430 ha and located in the municipality of Sao Goncalo do Amarante, State of Rio Grande do Norte. The property is located on the coast of the State of Rio Grande do Norte, with average annual rainfall, temperature and relative humidity of 1500 mm, 26[degrees]C and 78%, respectively.
At the experimental station, the animals were kept on grazing, with a forage diet variation according to the time of year. In general, in the rainy season, the forage supply was based on Brachiaria sp. pasture and in the dry season, on the remaining pasture of the rainy season and sorghum or corn silage. In both periods, the cows were supplemented with concentrate composed of soybean meal and soybean hull, corn grain, cotton cake, urea, corn meal, salt and mineral supplement.
Milk samples were taken directly from the cooling tank, always after the homogenization of the milk. Each sample was collected in duplicate; one of the samples, containing Bronopol[R], was sent to the laboratory accredited by the Brazilian Milk Quality Network (RBQL), in which somatic cell counts and milk composition analysis were performed. Fat, protein, lactose, total solids analyses were performed by the Fourier transform spectroscopy method (Bentley Nexgen[R], Bentley Instruments, USA).
Data on milk production and constituents used in the genetic evaluations were obtained in monthly dairy controls, from which the cumulative and corrected yields were calculated for 305 days of lactation. For Gyr, the analysis period was from 1985 to 2010 and for Guzera from 1994 to 2009. The DEPs for milk, fat and protein production were estimated annually up to 2012.
It is important to point out that the National Breeding Program of Dairy Gyr (PNMGL) began in 1985 with the progeny test; however, only in 1994 the first genetic evaluation was carried out, based on the dairy control data of the heifers of the first batch of the test bulls of ABCGIL and the ABCZ database. In 1999, analyses for protein, lactose and total solids of the milk were started for progeny test animals, although some of this information was already available at the ABCZ database. The phenotypic and genetic trends for milk constituents are, therefore, presented from the year 1994.
The statistical model used in the evaluation of animals for milk traits included the fixed effects of herd, year of calving, calving season, genetic group and the age of the dam at calving; as random effects, in addition to the residual effect, we considered the effect of animal (cow, mother and father) and the permanent effect of environment.
Yijklm= [mu] + ADi + EPi + GCj+ Ek + CGl +
+ b1(Iijklm-I) + b2(Iijklm-I)2 + [epsilon]ijklm
Yijklm= observation of the nth cow;
[mu] = overall mean;
ADi = direct additive random effect;
EPi = permanent environmental random effect;
GCj = fixed effect of contemporary group (herd-year of calving);
Ek = fixed effect of the calving season;
CGl = effect of the genetic group;
Iijkl = effect of cow's age at calving as a covariate;
b1 = linear regression coefficient for cow's age;
b2 = coefficient of quadratic regression for cow's age;
[epsilon]ijkl = random residue associated with each observation.
The DEPs for milk, fat and protein production are estimated annually from equations of the mixed models, with inclusion of the kinship matrix, using the derivative-free restricted maximum likelihood method, allowing to obtain BLUP estimates (Peixoto et al., 2013; Verneque et al., 2013). The statistical model used in the genetic evaluation of animals for milk traits included the fixed effects of herd, year of calving, calving season, genetic composition and age of the dam at calving; and, as random effects, in addition to the residue, the effect of animal (cow, mother and father) and the permanent effect of environment. For the traits of fat and protein production, a bitrait analysis was applied to milk production. For this, the MTDFReML software was used (Boldman, Kriese, & Van Vleck, 1995).
Results and discussion
During the execution period of breeding programs of the breeds, there was a significant difference in the means of fat and lactose between Gyr (1985 to 2010) and Guzera (1994 to 2009) breeds (Table 1).
These results were presented to characterize the milk level of the breeds of the EMPARN population and to illustrate the advances. The results are below those reported in the literature (Rangel, Guedes, Albuquerque, Novaes, & Lima Junior, 2009; Ribeiro, Tinoco, Lima, Guilhermino, & Rangel, 2009; Galvao Junior et al., 2010; Rompa & Cavalari, 2012) for the breeds in question, which can be attributed mainly to environmental differences.
Figures 1, 2 and 3 illustrate the trends for milk production at 305 days of lactation (PL305), fat production (PG305) and protein production (PP305) accumulated in 305 days of lactation of the Gyr herd of EMPARN participating in the National Breeding Program of Dairy Gyr (PNMGL).
The annual genetic and phenotypic gains for PL305 in the Gyr herd were 26.8 and 88.7 kg [year.sup.-1], respectively. Although positive and expressive, the genetic gain (1.8%) of the Gyr herd still falls short of the possible annual genetic gains with selection for milk production traits of up to 3% of the population mean (Smith, 1985). This result suggests that the selection pressure on the herd cows and the bulls that will be progenitors of the next generation can be increased in order to allow the animals of greater genetic merit for dairy traits to have a greater chance of remaining in the herd and to maximize genetic changes. Phenotypic gain, sum of the effects of genetic and environmental factors, reveals that in addition to the genetic changes, favorable changes in milk production occurred in the environmental conditions. Balieiro et al. (2000), obtained annual genetic and phenotypic gains of 10.46 and 66.93 kg [year.sup.-1], respectively, for the population of Gyr cows selected for milk in 19 herds, which were lower than those achieved by the EMPARN herd for the milk production trait. The highest genetic gain for PL305 was verified in the period from 2001 to 2007.
In the trend curve for fat (Figure 2), there was a positive trend for both the mean DEP and the phenotypic mean for fat production, that is, both presented growth, despite annual fluctuations. The phenotypic gain for fat production was 5.85 kg [year.sup.-1] and the genetic gain was 1.08 kg [year.sup.-1]. These values of selection for production, by selection pressure, promote changes in the production of milk components. The genetic improvements in the direct selection for milk production advance much more than the improvements in the production of the constituents and, therefore, when we deal with the production of constituents as a percentage, we have a reduction. Boligon et al. (2005) report a negative genetic correlation between milk yield and fat percentage, that is, the increase in milk production is not accompanied by an increase in milk fat content.
In the same way, the production of protein had an expressive increase, with a positive genetic trend. Annually, the phenotypic and genetic gains were respectively 4.10 and 0.45 kg. The phenotypic gains for protein production were greater than the genetic gains due to favorable environmental conditions in which the animals were bred.
Despite the annual fluctuations, our findings showed that, during the last decades, there was a positive trend for genetic progress for PL305, PG305 and PP305 in the EMPARN Gyr breeding herd (Figures 1, 2 and 3). This scenario is similar to that occurring throughout Brazil, where the Gyr breed has undergone vigorous selection for productive traits (Santana Junior et al., 2014).
Figures 4, 5 and 6 show, respectively, trends for milk production (PL305), fat production (PG305) and protein production (PP305) accumulated in 305 days of lactation of the EMPARN Guzera herd participant in the National Breeding Program of Guzera (PNMGul), since its inception in 1994. It should be noted that the first bulls summary of the progeny test was published in 2000. Therefore, this year is used as the highlight for the beginning of the trend curve presented in the figures below.
A positive genetic and phenotypic gain was observed for milk production from the year 2000, which has been maintained over the years. The annual phenotypic gain was 37.68 kg and the genetic gain was 8.44 kg per year. Although lower than in the Gyr breed, the genetic and phenotypic gains are positive for the Guzera breed. These results also suggest that the selection pressure on the progenitors can be increased.
In the trend curve for fat (Figure 5), a positive trend line was verified for the mean of DEP, but for the phenotypic mean it was negative. The genetic gain was 0.36 kg year and the phenotype gain was 1.79 kg year.
The protein production (Figure 6) also showed an expressive increase for the genetic trend, since the phenotypic trend remained low and constant, the phenotypic and genetic gains were respectively -0.94 kg [year.sup.-1] and 0.25 kg.
The positive response registered in the genetic gain for the production of milk constituents resulted from the high genetic correlation with milk production, the focus of the selection. The same was not verified for the negative phenotypic trend (fat: -2.3 kg [year.sup.-1] and protein: -1.1 kg [year.sup.-1]). For Peixoto et al. (2010), the intensive use of some breeding bulls has led to an increase in inbreeding coefficient and genetic drift in the population, a fact that may be associated with a negative phenotypic trend.
Nevertheless, the progress of the herd occurred under harsh conditions, representative of the Northeast region and, therefore, the EMPARN herd constitutes an important genetic resource for the improvement of commercial dairy herds in the state of Rio Grande do Norte.
The selection applied in the last decades to the herds of the Gyr and Guzera breeds in the EMPARN, RN, herds under adverse environmental conditions and with low input use, has resulted in satisfactory genetic and phenotypic advances for dairy traits, milk, fat and protein production.
Balieiro, E. S., Pereira, J. C. C., Valente, J., Verneque, R. S., Balieiro, J. C. C., & Ferreira, W. J. (2000). Estimativas de parametros geneticos e de tendencias fenotipica, genetica e de ambiente de algumas caracteristicas produtivas da raca Gir leiteiro. Arquivo Brasileiro de Medicina Veterinaria e Zootecnia, 52(3), 266-275.
Boldman, K. G., Kriese, L. A., & Van Vleck, L. D. (1995). A manual for use of MTDFREML: a set of programs to obtain estimates of variances and covariances (DRAFT). Lincoln, NE: Department of Agriculture/Agriculture Research Service.
Boligon, A. A., Rorato, P. R. N., Ferreira, G. B. B., Weber, T., Kippert, C. J., & Andreazza, J. (2005). Herdabilidade e tendencia genetica para as producoes de leite e de gordura em rebanhos da raca holandesa no estado do Rio Grande do Sul. Revista Brasileira de Zootecnia, 34(5), 1512-1518.
Galvao Junior, J. G. B., Rangel, A. H. N., Medeiros, H. R., Silva, J. B. A., Aguiar, E. M., ... Lima Junior, D. M. (2010). Efeito da producao diaria e da ordem de parto na composicao fisico-quimica do leite de vacas de racas zebuinas. Acta Veterinaria Brasilica, 4(1), 25-30.
Lobo, R. N. B., Madalena, F. H., & Penna, V. M. (2000). Avaliacao de esquemas de selecao alternativos para bovinos zebus de dupla aptidao. Revista Brasileira de Zootecnia, 29(5), 1361-1370.
Madalena, F. H. (2000). Valores economicos para a selecao de gordura e proteina do leite. Revista Brasileira de Zootecnia, 29(3), 678-684.
Malhado, C. H. M., Carneiro, P. L. S., Pereira, D. G., & Martins Filho, R. (2008). Progresso genetico e estrutura populacional do rebanho Nelore no Estado da Bahia. Pesquisa Agropecuaria Brasileira, 43(9), 1163-1169.
Peixoto, M. G. C. D., Poggian, C. F., Verneque, R. S., Egito, A. A., Carvalho, M. R. S., Penna, V. M., ... Machado, M. A. (2010). Genetic basis and inbreeding in the Brazilian Guzerat (Bos indicus) subpopulation selected for milk production. Livestock Science, 131 (2-3), 168-174.
Peixoto, M. G. C. D., Santos, G. G., Bruneli, F. A. T., Penna, V. M., Verneque, R. S., Panetto, J. C. C., ... Carvalho, M. R. S. (2013). Programa Nacional de Melhoramento do Guzera para Leite: resultados do Teste de Progenie, do Programa de Melhoramento Genetico de Zebuinos da ABCZ e do Nucleo MOET. Juiz de Fora, MG: Embrapa Gado de Leite.
Peixoto, M. G. C. D., Verneque, R.S., Teodoro, R. L., Penna, V. M., & Martinez, M. L. (2006). Genetic trend for milk yield in Guzerat herds participating in progeny testing and MOET nucleus schemes. Genetics and Molecular Research, 5(3), 454-465.
Rangel, A. H. N., Guedes, P. L. C., Albuquerque, R. P. F., Novaes, L. P., & Lima Junior, D. M. (2009). Desempenho produtivo leiteiro de vacas Guzera. Revista Verde de Agroecologia e Desenvolvimento Sustentavel, 4(1), 85-89.
Ribeiro, A. B., Tinoco, A. F. F., Lima, G. F. C., Guilhermino, M. M., & Rangel, A. H. N. (2009). Producao e composicao do leite de vacas Gir e Guzera nas diferentes ordens de parto. Revista Caatinga, 22(3), 46-51.
Rompa, P., & Cavalari, C. H. M. (2012). Volume e composicao do leite das racas Zebuinas. Cadernos de Pos-Graduacao da FAZU, 2(1), 1-7.
Santana Junior, M. L., Pereira, R. J., Bignardi, A. B., Ayres, D. R., Menezes, G. R. O., Silva, L. O. C., . Albuquerque, L. G. (2016). Structure and genetic diversity of Brazilian Zebu cattle breeds assessed by pedigree analysis. Livestock Science, 187(1), 6-15.
Santana Junior, M. L., Pereira, R. J., Bignardi, A. B., El Faro, L., Tonhati, H., & Albuquerque, L. G. (2014). History, structure, and genetic diversity of Brazilian Gir cattle. Livestock Science, 163(1), 26-33.
Santos, E. S., Leite, P. R. M., & Lobo, R. B. (1990). Progresso genetico do rebanho Gir Leiteiro de Umbuzeiro, PB. Pesquisa Agropecuaria Brasileira, 25(1), 1339-1348.
Smith, C. (1985). Rates of genetic change in farm livestock. Research Development Agricultural, 1(2), 79-85. Verneque, R. S., Panetto, J. C. C., Teixeira, R. B., Peixoto, M. G. C. D., Bruneli, F. A. T., Santos, G. G., ... Fernandes, A. R. (2013). Programa Nacional de Melhoramento do Gir
Leiteiro--sumario brasileiro de touros--resultado do teste de progenie. Juiz de Fora, MG: Embrapa Gado de Leite.
Received on April 12, 2017.
Accepted on June 9,2017.
Adriano Henrique do Nascimento Rangel (1), Thalita Polyana Monteiro Araujo (2), Guilherme Ferreira da Costa Lima (3), Luiz Henrique Fernandes Borba (1), Maria Gabriela Campolina Diniz Peixoto (4) and Dorgival Morais de Lima Junior (5) *
(1) Unidade Academica Especilaizada em Ciencias Agrarias, Universidade Federal do Rio Grande do Norte, Macaiba, Rio Grande do Norte, Brazil. (2) Departamento de Zootecnia, Universidade Federal Rural de Pernambuco, Recife, Pernambuco, Brazil. (3) Empresa de Pesquisa Agropecuaria do Rio Grande do Norte, Parnamirim, Rio Grande do Norte, Brazil. (4) Centro Nacional de Pesquisa em Gado de Leite, Empresa Brasileira de Pesquisa Agropecuaria, Juiz de Fora, Minas Gerais, Brazil. (5) Universidade Federal de Alagoas, Av. Manoel Severino Barbosa, s/n, 57309-005, Arapiraca, Alagoas, Brazil. * Author for correspondence: firstname.lastname@example.org and reproduction in any medium, provided the original work is properly cited.
Caption: Figure 1. Phenotypic and genetic trends of 305-d milk production for the Gyr herd.
Caption: Figure 2. Phenotypic (Pfat) and genetic (DEFP) trends for fat production of the Gyr herd.
Caption: Figure 3. Phenotypic (Pro) and genetic (DEPP) trends for protein production of the Gyr herd.
Caption: Figure 4. Phenotypic and genetic trends of 305-d milk production for the Guzera herd.
Caption: Figure 5. Phenotypic trend (Pfat) and genetic (DEFP) for fat production of the Guzera herd.
Caption: Figure 6. Phenotypic trend (Pro) and genetic (DEPP) for protein production of the Guzera herd.
Table 1. Means and standard deviations in milk composition for Gyr and Guzera breeds belonging to the herds of EMPARN. Breed Milk composition (%) Fat Protein Gyr 3.56 [+ or -] 0.79 (a) 2.91 [+ or -] 0.34 (a) Guzera 3.92 [+ or -] 0.92 (b) 2.95 [+ or -] 0.35 (a) Breed Lactose Dry extract Gyr 4.14 [+ or -] 0.25 (a) 11.77 [+ or -] 1.00 (a) Guzera 3.75 [+ or -] 0.30 (b) 11.07 [+ or -] 0.99 (a) Significance level p < 0.05 by Tuke'ys test.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||ANIMAL BREEDING AND REPRODUCTION|
|Author:||Rangel, Adriano Henrique do Nascimento; Araujo, Thalita Polyana Monteiro; Lima, Guilherme Ferreira d|
|Publication:||Acta Scientiarum. Animal Sciences (UEM)|
|Date:||Jan 1, 2018|
|Previous Article:||Supplementation with microencapsulated lemongrass essential oil improves protein deposition and carcass yield in silver catfish (Rhamdia quelen)/A...|
|Next Article:||Effect of simulated acid rain on the survival, mortality, behaviour and morphology of African mud catfish Clarias gariepinus (Burchell, 1822)/Efeito...|