ABSTRACT
The present study was carried out to evaluate the phenotypic diversity between fighting cocks; roosters of the naturalized ecotypes Graúna Dourada, Brejeira, Teresina, and Nordestina; roosters of the Sertaneja genetic (synthetic) group; and family-farmed roosters via descriptive and multivariate analyses. The qualitative traits evaluated by descriptive analysis indicated a prevalence of single comb in the Graúna Dourada, Teresina, Nordestina, and Sertaneja ecotypes, whereas the Brejeira ecotype, fighting cocks, and family-farmed roosters have a pea-type comb, which can be viewed as a sign of introgression of the fighting cock in flocks of the region. Dark plumage prevailed across the groups. The quantitative traits subjected to multivariate analysis by clustering criteria placed the four ecotypes and the synthetic type in a single group, the family-farmed animals were allocated to another group, and, lastly, the fighting cocks were placed in a third group. These results were similar to those generated by principal component analysis, which indicated body length, wingspan, and body weight as the traits of greatest variability to discriminate the different groups. It is important to stress the risk of genetic erosion to which naturalized chicken populations are exposed as a consequence of the introgression of fighting cocks in naturalized-chicken rearing facilities.
Keywords:
cluster analysis; morphometry; multivariate analysis; principal component analysis
Introduction
In Brazil, chicken considered “naturalized” are descendants of birds that have been introduced in the country since its discovery and colonization. Over time, these animals have undergone natural selection and adapted to the environment, exhibiting original and acquired traits.
The importance of regional birds has grown in the last few years. This is especially true for those from the northeast region of Brazil, which are also important for the selective breeding of flocks that constitute family poultry farming (Clementino et al., 2010Clementino, C. S.; Barbosa, F. J. V.; Carvalho, A. M. F.; Costa-Filho, R. A. R.; Silva, G. R.; Campelo, J. E. G.; Britto, F. B. and Diniz, F. M. 2010. Microsatellite DNA loci for population studies in Brazilian chicken ecotypes. International Journal of Poultry Science 9:1100-1106. https://doi.org/10.3923/ijps.2010.1100.1106
https://doi.org/10.3923/ijps.2010.1100.1...
). Therefore, the interest in maintaining groups with a defined genetic identity has increased, and conservation programs have thus been implemented. Before this process, however, producers must characterize the variability and potential to guide genetic breeding programs. This step is undertaken in accordance with methodologies proposed by Food and Agriculture Organization (FAO, 1981FAO. 1981. Descriptores de especies avícolas. p.13-15. In: Banco de datos de recursos genéticos animales. Roma, Italia.) and Francesch et al. (2011)Francesch, A.; Villalba, I. and Cartañá, M. 2011. Methodology for morphological characterization of chicken and its application to compare Penedesenca and Empordanesa breeds. Animal Genetic Resources 48:79-84. https://doi.org/10.1017/S2078633610000950
https://doi.org/10.1017/S207863361000095...
and works developed by Méndez et al. (2011)Méndez, Y.; Pons, A. and Francesch, A. 2011. Comparison of zoometric measurements in Balearic hens. Archivos de Zootecnia 60:445-448. https://doi.org/10.4321/S0004-05922011000300032
https://doi.org/10.4321/S0004-0592201100...
, based on morphometric measurements.
Information about the extent of genetic diversity in poultry is essential for the development of strategies to preserve and exploit these resources (Zaragoza et al., 2013Zaragoza, M. L.; Rodríguez H., J. V.; Hernández Z., J. S.; Perezgrovas, G. R.; Martínez C., B. and Méndez E., J. A. 2013. Characterization of hens batsi alak in the Highlands of Southeast Mexico. Archivos de Zootecnia 62:321-332. https://doi.org/10.4321/S0004-05922013000300001
https://doi.org/10.4321/S0004-0592201300...
).
No less important is the risk of genetic erosion to which local populations are exposed, or the extinction of isolated groups of those birds (Riztyan et al., 2011Riztyan; Katano, T.; Shimogiri, T.; Kawabe, K. and Okamoto, S. 2011. Genetic diversity and population structure of Indonesian native chickens based on single nucleotide polymorphism markers. Poultry Science 90:2471-2478. https://doi.org/10.3382/ps.2011-01450
https://doi.org/10.3382/ps.2011-01450...
). A situation that may be considered relevant in poultry facilities is the possibility of mating females belonging to family farms to fighting cocks. Both are genetic groups of regional importance, and when the useful lifespan of a fighting cock is over, the animal is incorporated into the breeding grounds of family farms. Because of their aggressiveness, fighting cocks overpower the males of naturalized groups. As a consequence, the genetic constitution of these populations is affected.
Because males remain in breeding facility for a longer period and leave more descendants, they have a much more impacting contribution to the genetic composition of the flock than females. Therefore, the use of information of males in the evaluation of the variability of a population is of paramount importance. The present study thus proposes to examine the phenotypic diversity of roosters of different genetic groups reared in family farming systems.
Material and Methods
Research on animals was conducted according to the institutional committee on animal use (case no. 404/17).
The study was conducted with chickens reared in the mid-North region of Brazil. Data were collected only from adult males. Morphological (plumage color, comb type, presence of crest, and ear tucks) and morphometric (lengths, circumferences, and heights of body parts) body traits were measured in four ecotypes and one synthetic group of naturalized birds and also in fighting cocks and in roosters from family farming facilities in the municipalities of Teresina (5°05'21.0" S and 42°48'07.0" W) and Campo Maior (4°49'40.0" S and 42°10'07.0" W), Piauí (PI), Brazil.
Data collected from the Graúna Dourada, Brejeira, Teresina, and Nordestina ecotypes and the Sertaneja synthetic group were measured in 34, 35, 28, 32, and 10 animals, respectively, located in Teresina.
The family farms were randomly sampled in the rural area of the municipalities of Teresina and Campo Maior, where measurements were taken in 30 and 10 animals, respectively. We also collected data of 106 fighting cocks rescued from fighting in Teresina.
The following qualitative traits were obtained in a visual assessment: plumage color, comb type, eye color, beak color, feet color, presence of crest, presence of naked neck, and presence of feathered feet.
Quantitative traits were measured following the methodology suggested by Francesch et al. (2011)Francesch, A.; Villalba, I. and Cartañá, M. 2011. Methodology for morphological characterization of chicken and its application to compare Penedesenca and Empordanesa breeds. Animal Genetic Resources 48:79-84. https://doi.org/10.1017/S2078633610000950
https://doi.org/10.1017/S207863361000095...
(Table 1). The data were tabulated on electronic spreadsheets for later editing of the work files and statistical analyses.
Qualitative traits were described by frequency analysis. Analysis of variance was applied for the quantitative traits, according to the statistical model below:
in which Yijk is the expected response of the dependent variable, μ is the population mean, αi is the effect of genetic group i, βj is the effect inherent to handler j, (αβ)ij is the interaction effect between the two factors, and eijk is the random error. The means of the traits for the model effects were compared with each other using Scott Knott's test (P<0.05).
To evaluate the phenotypic diversity of the genetic groups, the means of each body morphometric trait and live weight of each genetic group were estimated. Means were subjected to principal component analysis to estimate the eigenvalues and eigenvectors and then identify the traits that allow for the differentiation between genetic groups.
Next, the means of the traits were standardized by subtracting their average and dividing by the standard deviation. Results were used to calculate the Euclidean distance between genetic groups and were subjected to the Average Linkage (UPGMA) and Ward's (1963) clustering methods and Tocher's optimization method, using R software version 3.4.2 (R Core Team, 2017R Core Team. 2017. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.).
In the cluster analyses using the first two criteria, we calculated the AU (approximately unbiased) values by multiscale bootstrap resampling and BP (bootstrap probability), which is less precise than the AU value. Clusters with high AU values (e.g., 95%) can be considered as strongly supported by the data (Suzuki and Shimodaira, 2006Suzuki, R. and Shimodaira, H. 2006. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22:1540-1542. https://doi.org/10.1093/bioinformatics/btl117
https://doi.org/10.1093/bioinformatics/b...
).
The same software was adopted to calculate the cophenetic correlation coefficient to evaluate the consistency of clustering and of the frequency homogeneity analysis of qualitative traits and to calculate the estimates of frequencies, analysis of variance, and Scott Knott's test of means.
Results
There was no significant interaction between the effects of genetic groups and the handlers who took the measurements. In the evaluation of feathering type, no “naked neck” was observed, and thus full feathering prevailed (Figure 1).
Photographs of Graúna Dourada (A), Nordestina (B), Brejeira (C), and Teresina (D) ecotype roosters; Sertaneja genetic group rooster (E); family-farmed rooster (F); and fighting cock (G).
Dark-colored plumage prevailed among the roosters (Table 2). The Graúna Dourada ecotype showed 100% black plumage, whereas the Teresina ecotype had 100% red plumage. Red plumage prevailed in the Brejeira ecotype (42.86%), while the Nordestina group was distinct from the others for having a higher percentage of animals with light plumage, in which the white color was the most common (50%). The fighting cocks showed a higher percentage of black (35.85%) and red (34.91%) plumage.
Qualitative phenotypic traits of roosters of genetic groups reared in mid-Northern Brazil, 2018
There was a predominance of birds with light eyes and skin, except for the Sertaneja group, the roosters reared in Campo Maior, and the Grauna Dourada ecotype, in which 70, 60, and 52.94% of the animals were dark-eyed, respectively. It is noteworthy that fighting cocks had a reddish drumstick. In terms of beak color, the Graúna Dourada ecotype roosters showed dark (brown and black) beaks, whereas the other groups had a larger proportion of yellow beak and the Nordestina ecotype had a light-colored beak. Another noteworthy finding was the prevalence of yellow shank with reddish tones.
Lacking crest, presence of ear tuft, naked neck, and feathered shank prevailed in all genetic groups evaluated. The frequency of yellow-pinkish shank was homogenous across the genetic groups.
The single- and pea-type combs were observed in the studied genetic groups. Single combs were found in 100% of the Graúna Dourada, Nordestina, and Teresina ecotypes and in the synthetic group Sertaneja. Pea-type combs occurred in 100% of the Brejera animals and fighting cocks. The same type of comb occurred with a greater frequency in the family-farmed animals evaluated in Teresina (96.67%) and Campo Maior (90%).
The high coefficient of variation for comb height (60.13%; Table 3) clearly demonstrates the existence of differences between the comb types present in the roosters of the evaluated groups. The mean value of this trait for those animals was lower and similar to that of the roosters from Teresina.
The highest live weight was obtained by fighting cocks, whose mean values were similar to those of the family-farmed animals. The average live weights of the Graúna Dourada, Nordestina, Brejeira, and Teresina ecotype roosters were 2.309, 2.297, 2.212, and 2.404 kg, respectively.
Live weight showed large variations, differing between the evaluated genetic groups. Principal component analysis of the traits (Table 3) showed that the first component explained 97.91% of the existing variability between genetic groups, and body weight was the trait that most contributed in this component.
Given the elevated percentage variability attributed to the live weight of the animals in the first component and the inexpressive participation of the other traits, in addition to the fact that this trait was influenced by the environment, a second analysis was undertaken excluding body weight. The analysis revealed that the first component explained 80.16% of the variability present in the groups (Figure 2). Together, the first two components explained 91.47% of the variability. The traits with the highest percentage variability in the first and in the second components were wingspan and body length.
Scatterplot derived from principal component analysis of fighting cocks and naturalized roosters of the mid-North region of Brazil, 2018.
Figure 3 depicts the dendrograms resulting from the Ward's and Average cluster analyses, which were obtained with the same traits as those of the first principal component analysis.
Dendrograms of phenotypic similarity between fighting cocks, naturalized ecotypes, and roosters from municipalities of mid-Northern Brazil according to Ward's method (Ward, 1963Ward, J. H. 1963. Hierarchical grouping to optimize an objective function. Journal of the American Statistical Association 58:236-244. https://doi.org/10.1080/01621459.1963.10500845
https://doi.org/10.1080/01621459.1963.10... ; A) and the Average method (B) and respective approximately unbiased (au) and bootstrap probability (bp) values.
By observing the P-values of AU, we can visualize the animals in two big groups. The first group is composed of the four ecotypes of naturalized chickens and the Sertaneja genetic group, whereas the second contains fighting cocks and family-farmed roosters from Teresina and Campo Maior. A similar finding was obtained when we removed the body weight trait. The qualitative traits showed little variation, which was why they were not included in those analyses.
In this study, cophenetic correlation coefficients of 0.98, 0.96, and 0.95 were estimated according to the Average, Ward's, and Tocher's optimization methods, respectively.
Discussion
According to Valdés Corrales et al. (2010)Valdés Corrales, R. J.; Pimentel, O.; Martínez, K. and Ferro, E. M. 2010. Phenotypic characterization of creole poultry genetic basis from San Andrés, Pinar del Río, Cuba. Archivos de Zootecnia 59:597-600., single comb and normal feathering are characteristics associated with recessive genes and also with greater ability to adapt to higher environmental temperatures, revealing an adaptive advantage of these animals to climatic adversities.
The prevalence of these colors is similar to what was found in a study of characterization of native birds in Cuba, in which 77.1% of the studied animals had dark plumage. This plumage tone was associated with the advantage of better camouflage in the presence of predators and with the inexistence of crosses between these populations and commercial strains (Pérez et al., 2004Pérez, A.; Polanco, G. and Pérez, Y. 2004. Morphological characteristics of local chicken ecotypes in Villa Clara Province in Central Cuba. Livestock Research for Rural Development 16(10).).
The largest number of colors was found in the family-farmed animals, which is similar to the results observed in birds from districts in the mountains in the south of Tanzania, where animals with multi-colored plumage are most frequently found (50.8%). The other plumages are present in low numbers, such as those composed of the colors black (18.6%), maroon (9.81%), white (8.37%), gray (6.35%), and red (2.16%) (Guni and Katule, 2013Guni, F. S. and Katule, A. M. 2013. Characterization of local chickens in selected districts of the Southern Highlands of Tanzania: I. Qualitative characters. Livestock Research for Rural Development 25(9).).
The lack of feathering on the feet may indicate a lack of influence of Asian ancestors. In this regard, Rodrigues et al. (2006)Rodrigues, F. P.; Queiroz, S. A. and Duarte, J. M. B. 2006. Genetic relatedness among wild, domestic and Brazilian fighting roosters. Revista Brasileira de Ciência Avícola 8:83-87. https://doi.org/10.1590/S1516-635X2006000200003
https://doi.org/10.1590/S1516-635X200600...
evaluated fighting animals in Bahia state and confirmed a common ancestry between domestic Brazilian chickens and fighting cocks, which was a product of the cross between Asian and Indian strains.
The phenotypic uniformity for the traits comb type and presence of ear tuft in the naturalized ecotypes may be regarded as an indicator of maintenance of phenotypic identity in those groups, which were evaluated only in a conservation center.
If we consider the type of heritability of the comb type trait, it will not be possible to monitor the occurrence of genetic erosion caused by fighting cocks based on visual analysis. This is due to the random or non-random matings in the Brejeira ecotype, which is the group most similar to the fighting cocks when seen from the side.
The hypothesis of presence of fighting cocks in the composition of animals from the family farms in Teresina and Campo Maior is consistent, given the observed high percentage of animals with a pea-type comb. It should be emphasized that the comb type of the animals from the above-mentioned municipalities might have originated from crosses with strains from other regions.
Another noteworthy fact in this case is that the comb type of the Teresina ecotype does not prevail in those breeding farms, since the group was thus named due to its origins in breeding farms in Teresina, selected from nine birds that showed phenotypic uniformity. This ecotype has been conserved with a view to preserving some originally found traits.
It can thus be assumed that the frequency of occurrence of single comb in birds that had been reared in the region is high. However, because this was not observed, it is possible that matings with fighting cocks had a marked contribution to the frequency of pea-type comb alleles in the two sampled locations as well as alleles of other body traits. In this way, there is a greater possibility of the animals from Teresina and Campo Maior exhibiting similar combs and body characteristics to those of fighting animals or that these animals are potentially strong candidates to cause genetic erosion in extensive breeding flocks of the region.
The presence of cockfighting in the region is expressive, even though it is a non-officially recognized activity. Fighting cocks were brought to Brazil around the year 1530 by Spaniards and quickly spread across the country (Escobar et al., 2014Escobar, M. L.; Aguiar, J. O. and Zagui, P. A. 2014. Galos em combate na Paraíba: O descumprimento da legislação ambiental. Revista Direitos Humanos e Democracia 2:143-165. https://doi.org/10.21527/2317-5389.2014.4.143-165
https://doi.org/10.21527/2317-5389.2014....
). The aggravating factor from the perspective of regional genetic resource conservation is how they are “discarded” after their useful lifespan as fighting cocks. Their destination is usually family farms, which is the prevalent chicken farming system in the region. Because these roosters are aggressive, they easily impose themselves as the main breeders of the breeding ground, consequently having great potential to cause genetic erosion in the populations of the region.
From the economic standpoint, this is an unfavorable trait in birds, because birds with a small comb tend to have an increase in body temperature caused by a less efficient heat dissipation system, which can interfere negatively with the fertility of breeder chickens (Rosa et al., 2010Rosa, A. P.; Pilecco, M.; Lucca, W. and Uttpatel, R. 2010. Níveis de cálcio no desempenho de machos reprodutores de corte com e sem crista. Ciência Rural 40:2174-2180. https://doi.org/10.1590/S0103-84782010005000155
https://doi.org/10.1590/S0103-8478201000...
).
In comparison with other native chickens, live weight values indicate that those ecotypes are classified as medium-sized. This trait has a genetic component but is also greatly influenced by the environment. This result was similar to what occurred with Peloco chickens, a genetic group naturalized in the state of Bahia that showed high phenotypic variability for carcass traits (Almeida et al., 2013Almeida, E. C. J.; Carneiro, P. L. S.; Wenceslau, A. A.; Farias Filho, R. V. and Malhado, C. H. M. 2013. Características de carcaça de galinha naturalizada Peloco comparada a linhagens de frango caipira. Pesquisa Agropecuária Brasileira 48:1517-1523. https://doi.org/10.1590/S0100-204X2013001100013
https://doi.org/10.1590/S0100-204X201300...
; Almeida et al., 2014Almeida, E. C. J.; Carneiro, P. L. S.; Malhado, C. H. M.; Pereira, A. H. R.; Farias, R. V.; Campos, B. M. and Silva, O. L. 2014. Relative contribution of carcass biometric characteristics phenotypic in diversity among native chickens Peloco and genotypes comercial. Actas Iberoamericanas de Conservación Animal 4:9-11.).
The results obtained with principal component analysis showed the animals distributed into three groups, in the analyses with and without live weight. Fighting cocks were set apart from the others. The second group was formed by family-farmed animals from Teresina and Campo Maior, and the third group consisted of the naturalized ecotypes and the synthetic group Sertaneja. The same result was obtained by cluster analysis using Tocher's optimization method, which showed consistency in the dispersion across the groups visualized in principal component analysis. It can be observed that the fighting cocks had the greatest body length and the roosters from Teresina and Campo Maior showed the largest wingspan.
Cophenetic correlation values greater than 0.70 indicate lower distortions in the formed clusters (Sokal and Rohlf, 1962Sokal, R. R. and Rohlf, F. J. 1962. The comparison of dendrograms by objective methods. Taxon 11:33-40. https://doi.org/10.2307/1217208
https://doi.org/10.2307/1217208...
). However, these estimated coefficients close to 1 might have been influenced by the utilized data, since we used the means of groups instead of the individual values of the animals belonging to each genetic group.
The existing similarity between the animals from Teresina and Campo Maior and of the fighting cocks reinforces the affirmation of the presence of those animals in family farms. As observed in loco, because these are aggressive animals, they tend to overpower the present roosters as breeders.
Conclusions
There is phenotypic variability between the roosters of the genetic groups evaluated by quantitative and qualitative traits. The greater phenotypic similarity between fighting cocks and the animals belonging to family farms in comparison with the evaluated ecotypes may indicate a risk of genetic erosion to which the populations of naturalized chickens are exposed due to the destination of fighting cocks in the region.
Acknowledgments
The authors thank the Núcleo de Conservação de Galinhas Naturalizadas do Meio-Norte do Brasil (NUGAN-MN), financed by Banco do Nordeste, for the provided data.
References
- Almeida, E. C. J.; Carneiro, P. L. S.; Malhado, C. H. M.; Pereira, A. H. R.; Farias, R. V.; Campos, B. M. and Silva, O. L. 2014. Relative contribution of carcass biometric characteristics phenotypic in diversity among native chickens Peloco and genotypes comercial. Actas Iberoamericanas de Conservación Animal 4:9-11.
- Almeida, E. C. J.; Carneiro, P. L. S.; Wenceslau, A. A.; Farias Filho, R. V. and Malhado, C. H. M. 2013. Características de carcaça de galinha naturalizada Peloco comparada a linhagens de frango caipira. Pesquisa Agropecuária Brasileira 48:1517-1523. https://doi.org/10.1590/S0100-204X2013001100013
» https://doi.org/10.1590/S0100-204X2013001100013 - Clementino, C. S.; Barbosa, F. J. V.; Carvalho, A. M. F.; Costa-Filho, R. A. R.; Silva, G. R.; Campelo, J. E. G.; Britto, F. B. and Diniz, F. M. 2010. Microsatellite DNA loci for population studies in Brazilian chicken ecotypes. International Journal of Poultry Science 9:1100-1106. https://doi.org/10.3923/ijps.2010.1100.1106
» https://doi.org/10.3923/ijps.2010.1100.1106 - Escobar, M. L.; Aguiar, J. O. and Zagui, P. A. 2014. Galos em combate na Paraíba: O descumprimento da legislação ambiental. Revista Direitos Humanos e Democracia 2:143-165. https://doi.org/10.21527/2317-5389.2014.4.143-165
» https://doi.org/10.21527/2317-5389.2014.4.143-165 - FAO. 1981. Descriptores de especies avícolas. p.13-15. In: Banco de datos de recursos genéticos animales. Roma, Italia.
- Francesch, A.; Villalba, I. and Cartañá, M. 2011. Methodology for morphological characterization of chicken and its application to compare Penedesenca and Empordanesa breeds. Animal Genetic Resources 48:79-84. https://doi.org/10.1017/S2078633610000950
» https://doi.org/10.1017/S2078633610000950 - Guni, F. S. and Katule, A. M. 2013. Characterization of local chickens in selected districts of the Southern Highlands of Tanzania: I. Qualitative characters. Livestock Research for Rural Development 25(9).
- Méndez, Y.; Pons, A. and Francesch, A. 2011. Comparison of zoometric measurements in Balearic hens. Archivos de Zootecnia 60:445-448. https://doi.org/10.4321/S0004-05922011000300032
» https://doi.org/10.4321/S0004-05922011000300032 - Pérez, A.; Polanco, G. and Pérez, Y. 2004. Morphological characteristics of local chicken ecotypes in Villa Clara Province in Central Cuba. Livestock Research for Rural Development 16(10).
- R Core Team. 2017. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
- Rodrigues, F. P.; Queiroz, S. A. and Duarte, J. M. B. 2006. Genetic relatedness among wild, domestic and Brazilian fighting roosters. Revista Brasileira de Ciência Avícola 8:83-87. https://doi.org/10.1590/S1516-635X2006000200003
» https://doi.org/10.1590/S1516-635X2006000200003 - Rosa, A. P.; Pilecco, M.; Lucca, W. and Uttpatel, R. 2010. Níveis de cálcio no desempenho de machos reprodutores de corte com e sem crista. Ciência Rural 40:2174-2180. https://doi.org/10.1590/S0103-84782010005000155
» https://doi.org/10.1590/S0103-84782010005000155 - Riztyan; Katano, T.; Shimogiri, T.; Kawabe, K. and Okamoto, S. 2011. Genetic diversity and population structure of Indonesian native chickens based on single nucleotide polymorphism markers. Poultry Science 90:2471-2478. https://doi.org/10.3382/ps.2011-01450
» https://doi.org/10.3382/ps.2011-01450 - Sokal, R. R. and Rohlf, F. J. 1962. The comparison of dendrograms by objective methods. Taxon 11:33-40. https://doi.org/10.2307/1217208
» https://doi.org/10.2307/1217208 - Suzuki, R. and Shimodaira, H. 2006. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22:1540-1542. https://doi.org/10.1093/bioinformatics/btl117
» https://doi.org/10.1093/bioinformatics/btl117 - Valdés Corrales, R. J.; Pimentel, O.; Martínez, K. and Ferro, E. M. 2010. Phenotypic characterization of creole poultry genetic basis from San Andrés, Pinar del Río, Cuba. Archivos de Zootecnia 59:597-600.
- Ward, J. H. 1963. Hierarchical grouping to optimize an objective function. Journal of the American Statistical Association 58:236-244. https://doi.org/10.1080/01621459.1963.10500845
» https://doi.org/10.1080/01621459.1963.10500845 - Zaragoza, M. L.; Rodríguez H., J. V.; Hernández Z., J. S.; Perezgrovas, G. R.; Martínez C., B. and Méndez E., J. A. 2013. Characterization of hens batsi alak in the Highlands of Southeast Mexico. Archivos de Zootecnia 62:321-332. https://doi.org/10.4321/S0004-05922013000300001
» https://doi.org/10.4321/S0004-05922013000300001
Publication Dates
-
Publication in this collection
27 June 2019 -
Date of issue
2019
History
-
Received
01 Nov 2018 -
Accepted
08 Apr 2019