## Services on Demand

## Journal

## Article

## Indicators

## Related links

## Share

## Genetics and Molecular Biology

##
*Print version* ISSN 1415-4757*On-line version* ISSN 1678-4685

### Genet. Mol. Biol. vol.29 no.4 São Paulo 2006

#### http://dx.doi.org/10.1590/S1415-47572006000400012

**ANIMAL GENETICS RESEARCH ARTICLE**

**Genetic trends in the expected progeny difference of the asymptotic weight of Nelore females**

**Analía del Valle Garnero ^{I}; Cintia Righetti Marcondes^{I}; Ricardo José Gunski^{I}; Henrique Nunes de Oliveira^{II}; Raysildo Barbosa Lôbo^{III}**

^{I}Laboratório de Genética, Ciências Biológicas, Campus Porto Nacional, Universidade Federal de Tocantins, Porto Nacional, Tocantins, Brazil

^{II}Departamento de Melhoramento e Nutrição Animal, UNESP - Universidade Estadual Paulista, Botucatu, São Paulo, Brazil

^{III}Departamento de Genética, Faculdade de Medicina de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, São Paulo, Brazil

**ABSTRACT**

There are few studies on weight covering the full life cycle of Zebu cattle, and there is no entire growth description or mean growth pattern for animals belonging to this breed. In order to provide such data, 1,158 Nelore females born between 1985 and 1995 were weighed 14,563 times from birth to full growth maturity, in ten herds spread over seven Brazilian states. The Von Bertalanffy, Brody, logistic and Gompertz non-linear models were used to obtain the asymptotic weights (A) and the maturation rates (K). The (co)variance and breeding value components for A and K were obtained by using the multiple trait derivative free restricted maximum likelihood method under the animal model. Genetic trends were calculated in function of the mean expected progeny differences (EPD) for the trait (A or K) divided by the number of animals according to their year of birth. The genetic trends of the expected progeny difference with reference to the date of birth of the cows were, on average, -6.5g y^{-1} for A and 2.0g y^{-1} for K, close to zero as confirmed by the low (0.0023 to 0.003) coefficient of regression values. The curve parameters are recommended as a selection criterion to reach precocity and avoid adult weight increase in the female herd.

**Key words:** adult weight, bovine(s), genetic trends, maturation rate, simulation with growth curves.

**Introduction**

A selection aimed at increasing the rate of weight gain or weight converted in beef cattle upon slaughter will result in a weight increase in adult animals and, consequently, greater herd maintenance costs. All these aspects together can reduce the economic advantages of weight increase in slaughter cattle (Cartwright, 1970; Dickerson, 1978). On the other hand, some authors recommend selection as one of the means to alter the shape of the growth curve so as to obtain precocious animals that can reach their slaughter weight earlier, without increasing their size (Cartwright, 1970; Fitzhugh, 1976; Carrijo and Duarte, 1999).

Female herds have higher beef cattle production costs due to a lower reproduction rate per animal sold. The two most important components determining the economic efficiency of beef cows are milk production and mature weight (Dickerson, 1978; McMorris and Wilton, 1986; Montaño-Bermudez *et al*., 1990; Arango and Van Vleck, 2002). Therefore, studies on growth and maturity traits deriving from studies on growth curves should be used as additional information in genetic breeding programs (Fitzhugh, 1976; Kaps *et al*., 2000; Arango and Van Vleck, 2002). Most sire summaries from different breeding programs provide the expected progeny differences (EPD) for weights at different ages. However, some summaries only provide the adult weight EPD, e.g. the 2002 US Angus sire summary (http:// www.angussiresearch.com/index.html) and the Nelore sire summary under the Nelore Genetic Breeding Program (http://www.ancp.org.br/ sumarios/index.asp).

In Brazil, there are only a few studies on the full life cycle (from birth to maturity) weight of Zebu cattle, one such study being that of Oliveira *et al*. (2000) on Guzerá females and another being the study by Garnero *et al*. (2005) on Nelore females. Most of the studies have used data from cattle whose weighing was interrupted before maturity (Carrijo and Duarte, 1999; Santoro *et al*., 2005) and, therefore, represent neither a general description of the growth rate nor an average growth pattern for Zebus.

In most growth curve studies, different weights are used to plot the curve, and their parameters are estimated with the aim of using them as selection criterion. Nevertheless, in their study of Guzerá females Oliveira *et al*. (2000) concluded that it would be difficult to alter the growth curve through selection practiced under the parameters of the curve, since the maturation rate and the asymptotic weight showed high and negative correlations (genetic, environmental and phenotypic). Carrijo and Duarte (1999) observed similar results when analyzing Nelore and Chianina specimens. Nevertheless, according to the study of various taurine breeds by Fitzhugh (1976) the alteration of the growth curve depends on the degree of genetic flexibility of the shape of the curve, suggesting the use of an index which includes size and maturity traits.

For taurine breeds, using growth curves to describe the female growth process is a well-known and established practice (Brown *et al*. 1976; Fitzhugh, 1976; Doren and Baker, 1989). The Richards and Brody models have generally proved to be the most adequate, although for the Frisian x Zebu genetic groups the Von Bertalanffy and Richards models have proved to be more suitable (Perotto *et al*., 1992; Val *et al*., 2004). Oliveira *et al*. (2000) and Garnero *et al*. (2005) concluded that the Brody, Von Bertalanffy, logistic and the Gompertz models are adequate to describe the growth of Zebu females, with the Von Bertalanffy model being the most suited for Guzerá females and the Brody model for Nelore females.

Using different models, the Nelores have shown great variation in their growth curve parameters and variance (Carrijo and Duarte, 1999; Garnero *et al*., 2005). This evidence confirms the broad variability in genetic resources observed in Zebu cattle, which could lead to different demands affecting production systems, environmental.

Since the main goal of any animal breeding strategy is to change the population average in order to match the desired character, it is necessary to adopt a rational selection program that could result in permanent and stable changes in animal genotypes to accomplish this task. Genetic trends are statistical resources used by researchers to quantify the effects of selection in trait improvement (Falconer and Mackay, 1997). Genetic trends, therefore, can quantify variations in productive traits per unit of time. This is obtained from the alterations in the mean breeding values for sires and dams. Rosa *et al.* (2001) studied Nelore females born between 1982 and 1993 and observed a 0.42 kg y^{-1} increase in adult weight. However, the author pointed out that the genetic increase observed in adult weight had been smaller than the genetic progress reached for products at 365 and 550 days of age, 0.53 kg y^{-1} and 0.654 kg y^{-1} respectively. These results made it very clear that the selection criteria applied to the herds could be considered adequate.

The aim of this study was to predict the breeding values for asymptotic weight and maturity rate traits for Nelore females and to estimate the genetic trends of these curve parameters.

**Material and Methods**

The data set consisted of 14,563 weights and ages describing 1,158 Nelore females born between 1985 and 1995 in ten herds that were included in the Nelore Genetic Breeding Program (PMGRN). These herds were spread throughout seven Brazilian states (Bahia (11), Goiás (73), Mato Grosso do Sul (476), Minas Gerais (320), Pará (3), Paraná (2) and São Paulo (273)).

The animals included in the PMGRN were identified in order to ascertain their pedigree information, and submitted to measurements every three months. Young Nelore females were measured from birth to 21 months, and cows were also weighed from birth to the weaning of their offspring, as well as in April and October of each year. Herds adopt a breeding season, which lasts from 60 to 90 days between October and March, according to the region and year. The use of artificial insemination was frequent. Calves were generally weaned from seven to eight months of age, and all remaining health care and vaccination procedures fully complied with the regulations of the veterinary sanitary control board of each state.

To describe the growth curve, four non-linear mathematical models were used, under SAS's NLIN procedure (SAS, 1996): Brody, *Y _{t}* =

*A*(1 -

*Be*) + e (Brody, 1945); Von Bertalanffy,

^{-Kt}*Y*=

_{t}*A*(1 -

*Be*)

^{-Kt}^{3}+ e (Bertalanffy, 1957); logistic,

*Y*=

_{t}*A*(1 +

*e*)

^{-Kt}*+ e (Nelder, 1961) and Gompertz,*

^{-M}*Y*+ e (Laird, 1966).

_{t}= Ae^{–Be–Kt}In these expressions, Y represents the weight of the animal at a definite age (*t*); A is the asymptotic value for size as *t* ® ¥, generally interpreted as average size at maturity independently of short-term fluctuations in size due to extraneous environmental effects of climate and food supply; B is a scaling parameter (integration constant) established by the initial values of *y _{o}* and

*t*(this parameter is adjusted for the situation where

*y*¹ 0 and/or

_{o}*t*¹ 0 - for example, when only postnatal size records are available and

_{o}*t*is taken as birth); K is a function of the ratio between maximum growth rate and mature size, commonly referred to as a maturing index (the speed in which an animal approaches its adult size);

_{o}*M*is a variable in the general indeterminate form, but is the defining constant for special cases (

*M*= 1 for Brody,

*M*= 3 for Von Bertalanffy,

*M*® ¥ for Gompertz,

*M*= variable for logistic);

*e*is the logarithm in a natural base; e represents the random error associated with each observation (Fitzhugh, 1976).

Simple criteria were used to compare the fit of the data to the models. Then, the number of iterations, the mean square error (MSE) and the coefficient of determination were calculated.

To acquire the (co)variance and breeding values for asymptotic weight (*A*) and maturity rate (K), the multiple trait derivative free restricted maximum likelihood method (MTDFREML) was used, fitting an animal model and using the *Multiple Trait Derivative Free Restricted Maximum Likelihood* program (Boldman *et al*., 1995) in single-trait and bi-trait analyses. This model included the effects of the contemporary group (CG) as fixed, and the additive direct genetic effect as random.

The herd effect was not significant (p < 0.01). Therefore, the CGs were formed by animals born in the same state, in the same quarter (January-March, April-June, July-September, or October-December), in the same year and with the same feed (pasture, semi-intensive or intensive).

Animals that were not weighed at least once from 24 to 36 months of age, once between 36 and 48 months of age, and once after age 54 months were eliminated. All CGs that had only one observation were also eliminated, with a remaining universe of 153 group containing 1,087 animals for analyses. The relationship matrix had 1,900 animals.

In matrix terms, the animal model may be described as:

where *y* is an observation vector (*A* or K parameters); *X* is a fixed effect incidence matrix (CG); b is a fixed effect vector; Z is an incidence matrix of each animal additive direct genetic effect; a is a vector of random direct genetic effects; and e is a vector of random residual effects.

The heritabilities obtained by Garnero *et al*. (2005) were considered to predict breeding values (BV) and expected progeny differences (EPD). The estimates were 0.39 (Von Bertalanffy and Gompertz), 0.42 (Brody) and 0.40 (logistic) for A; and 0.42; .44; 0.41 and 0.39 for K (Von Bertalanffy, Brody, logistic and Gompertz, respectively).

Asymptotic weight and maturation rate genetic trends were obtained by regressing breeding values on the year of birth of the cow. The data used referred to females born from 1988 onwards (N = 981), the year in which PMGRN established measurement rules.

The animals with the lowest, highest, and next to average EPD (for these same traits) were identified, as well as those with values A, B and K. Growth curve and genetic trend graph simulations were created on Microsoft Excel for Windows (2000).

**Results and Discussion**

In general, all the models were adjusted to the weights observed. The Gompertz model was only applied when an animal did not reach the convergence criteria. The number of iterations was 10 for Brody, 13 for logistic, 14 for Von Bertalanffy and 17 for Gompertz. The mean square error was similar in all models: Brody (2572.7), Von Bertalannfy (2633.8), Gompertz (2696.1) and logistic (2763.3). The coefficient of determination was also similar in all models (0.977 or 0.978).

The mean breeding values, as well as the maximum and minimum values, showed high variability, both for the asymptotic weight and the maturation rate (Table 1). Standard deviations of around 27 kg may be observed in the Brody model for asymptotic weight (*A*), with differences from -85 kg to 230 kg from the predicted breeding values. For the other models - Von Bertalanffy, Gompertz and logistic - the magnitude of the values for the same trait was lower, the standard deviation varied between 16 and 18 kg, and the predicted breeding values ranged from -57 to 105 kg. The values predicted for the maturation rate (K) were similar in all models, with standard deviations of approximately 0.005 kg and predicted breeding values between -0.05 and 0.08 kg. Breeding values that are too low or too high can indicate a lack of direct selection on the traits.

The genetic trends of the breeding value (BV) for the asymptotic weight (over the cow's birth date) were negative, except in the Brody model, where it was close to zero (Figure 1). These results were confirmed by the low value of (Table 2). The asymptotic weight drop an average -0.12 kg/year, with the yearly average reduction of -0.15 kg when estimated according to the Von Bertalanffy model; 0.05 kg based on the Brody model; -0.18 kg based on the Gompertz model; and -0.20 kg based on the logistic model. The results were the opposite of those estimated by Rosa *et al*. (2001). Although these authors established the adult weight based on the first time the cow was weighed, they reached their limits from 4 to 12 years of age, not on the growth curve.

On the other hand, the genetic trends for the maturation rate BV (over the year of the cows' birth) were positive for all the models studied. The maturation rate increased an average of 0.04 g d^{-1}, with the average yearly growth being 0.04 g d^{-1} when estimated according to the Von Bertalanffy model, 0.01 g d^{-} for the Brody model, 0.05 g d^{-1} for Gompertz model and 0.06g d^{-1} for the logistic model. Therefore, we can see that the cows' adult weight did not increase (reduction of only 120 g y^{-1}) and the maturation rate is not being used in the selection, as a 15-gram increase per year is a very low value (Figure 2).

To illustrate the application of the growth curves in simulations with high, average and low EPD for both traits, results were obtained using the Brody model, and compared to the average growth curve (Figures 3 and 4).

The use of these traits as selection criteria is viable, not only for heritabilities close to 0.40 (Garnero *et al*., 2005) and for the high variability observed in breeding values, but also for the genetic progress that could be reached. For example, in a selection against the asymptotic weight, the adult weight would be reduced by approximately 50 kg (Figure 3); as a result, there would be a reduction in maintenance costs. In addition, in the selection in favor of the maturation rate, animals with faster growth would be identified earlier both in terms of growth and sexual aspects (approximately 10 months), at the same time, and would turn into lighter adults (Figure 4). Regardless of how great the correlation between weights at different ages (Marcondes *et al*, 2000), it is possible to make a selection for maturation rate, due to its great variability and maintenance of low adult weight.

The predicted breeding values for asymptotic weight and maturity rate in the analyzed population indicated - with a high degree of variation - the existence of genetic resources that can be adjusted for different productive systems. The negative or close-to-zero genetic trends referring to asymptotic weight and maturation rate may be indicating a lack of direct selection on these traits in Nelore beef cows.

**Acknowledgements**

The authors would like to thank the PMGRN breeders and the following Brazilian research funding entities: FINEP/BID, FAPESP, ANCP, CNPq and CAPES.

**References**

Arango JA and Van Vleck LD (2002) Size of beef cows: Early ideas, new developments. Gen Mol Res 1:51-63. [ Links ]

Bertalanffy L von (1957) Quantitative laws in metabolism and growth. The Quarterly Review of Biology 32:270-230. [ Links ]

Boldman KG, Kriese LA, Van Vleck LD, Van Tasssel CP and Kachman SD (1995) A Manual for Use of MTDFREML: A Set of Programs to Obtain Estimates of Variance and Covariance [DRAFT]. Agricultural Research Service, Lincoln, 120 pp. [ Links ]

Brody S (1945) Bioenergetics and Growth. Reinhold Publishing Corp., New York, 1023 pp. [ Links ]

Brown JE, Fitzhugh Jr HA and Cartwright TC. (1976) A comparison of nonlinear models for describing weight-age relationships in cattle. J Anim Sci 42:810-818. [ Links ]

Carrijo SM and Duarte FAM (1999) Description and comparison of growth parameters in Chianina and Nelore cattle breeds. Gen Mol Biol 22:187-196. [ Links ]

Cartwright TC (1970) Selection criteria for beef cattle for the future. J Anim Sci 30:706-711. [ Links ]

Dickerson GE (1978) Animal size and efficiency: Basic concepts. Anim Prod 27:367-379. [ Links ]

Doren PE and Baker JF (1989) Estimating parameters growth curves of Bulls. J Anim Sci 67:1432-1445. [ Links ]

Falconer DS and Mackay TFC (1997) Introduction to Quantitative Genetics. 4^{th} edition. Longman, Edinburgh, 464 pp. [ Links ]

Fitzhugh HA Jr (1976) Analysis of growth curves and strategies for altering their shape. J Anim Sci 42:1036-1051. [ Links ]

Garnero A del V, Marcondes CR, Bezerra LAF, Oliveira HN and Lôbo RB (2005) Parâmetros genéticos da taxa de maturação e do peso assintótico de fêmeas da raça Nelore. Arq Bras Med Vet Zootec 57:652-662. [ Links ]

Kaps M, Herring WO and Lamberson WR (2000) Genetic and environmental parameters for traits derived from the Brody growth curve and their relationships with weaning weight in Angus cattle. J Anim Sci 78:1436-1442. [ Links ]

Laird AK (1966) Dynamics of relative growth. Growth, Bar Habor 29:249-263. [ Links ]

Marcondes CR, Bergmann JAG, Eler JP, Pereira JCC, Ferraz JBS and Penna VM (2000) Análise de alguns critérios de seleção para características de crescimento na raça Nelore. Arq Bras Med Vet Zootec 52:83-89. [ Links ]

McMorris MR and Wilton JW (1986) Breeding system, cow weight and milk yield effects on various biological variables in beef production. J Anim Sci 63:1361-1372. [ Links ]

Montaño-Bermudez M, Nielsen MK and Deutscher GH (1990) Energy requirements for maintenance of crossbred beef cattle with different genetic potential for milk. J Anim Sci 68:2279-2288. [ Links ]

Nelder JA (1961) The fitting of a generalization of the logistic curve. Biometrics 17:89-110. [ Links ]

Oliveira HN, Lôbo RB and Pereira CS (2000) Comparação de modelos não lineares para descrever o crescimento de fêmeas da raça Guzerá. Pesq Agropec Bras 35:1843-1851. [ Links ]

Perotto D, Cue RI and Lee AJ (1992) Comparison of nonlinear functions for describing the growth curve of three genotypes of dairy cattle. Can J Anim Sci 72:773-782. [ Links ]

Rosa AN, Lôbo RB, Oliveira HN, Bezerra LAF and Borjas AR (2001) Peso adulto de matrizes em rebanhos de seleção da raça Nelore no Brasil. Rev Bras Zootec 30:1027-1036. [ Links ]

Santoro KR, Barbosa SBP, Brasil LHA and Santos ES (2005) Estimativas de parâmetros de curvas de crescimento de bovinos Zebu, criados no estado de Pernambuco. R Bras Zootec 34:2262-2279. [ Links ]

SAS (1996) Statistical Analysis System, 6.12. Versão for Windows®, SAS Institute, Cary. [ Links ]

Val JE, Freitas MAR, Oliveira HN, Cardoso VL, Machado PF and Paneto JCC (2004) Indicares de desempenho em rebanho da raça Holandesa: Curvas de crescimento e altura, características reprodutivas, produtivas e parâmetros genéticos. Arq Bras Med Vet Zootec 56:86-93. [ Links ]

**Send correspondence to**

Analía del Valle Garnero

Laboratório de Genética

Ciências Biológicas

Campus Porto Nacional, UFT

Rua 03, Q17, Jardim dos Ipês

77500-000, Porto Nacional, Tocantins, Brazil

E-mail: analiagarnero@yahoo.com.br

Received: April 13, 2005; Accepted: May 24, 2006.

Associate Editor: Pedro Franklin Barbosa