Microsoft Word - 22-Agra_26496_Vet 1475 Original Article Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 SEXUAL DIMORPHISM OF LABRADOR RETRIEVER DOGS BY MORPHOMETRY DIMORFISMO SEXUAL DE CÃES DA RAÇA RETRIEVER DO LABRADOR POR MEIO DA MORFOMETRIA Murilo Antonio Oliveira THULLER1; Marcelo JANGARELLI2; Douglas Mena do COUTO3; Alexandre Herculano Borges ARAÚJO4 1. Mestrando em Produção Animal, Departamento de Zootecnia da Universidade Federal Rural do Rio de Janeiro, Seropédica, RJ, Brasil. murilothuller@zootecnista.com.br; 2. Professor, Doutor, Departamento de Matemática, Instituto de Ciências Exatas da Universidade Federal Rural do Rio de Janeiro, Seropédica, RJ, Brasil; 3. Mestrando em Produção Animal, Departamento de Zootecnia da Universidade Federal Rural do Rio de Janeiro, Seropédica, RJ, Brasil; 4. Professor, Doutor, Departamento de Nutrição e Alimentação Animal, Instituto de Zootecnia da Universidade Federal Rural do Rio de Janeiro, Seropédica, RJ, Brasil. ABSTRACT: The domestic dog (Canis familiaris) is the species of greatest morphological diversity among mammals. Seventy-four Labrador Retriever dogs- 27 males and 47 females – were used in this experiment. Thirty quantitative biometric characteristics, related to morphology were measured. The objective of this study was to evaluate the morphometric traits of the Labrador Retriever breed to establish descriptive biometric attributes that may show sexual dimorphism through principal component analysis (PCA) and discriminant analysis (DA). The PCA was processed using all the variables and performing a pre-selection of the most correlated variables. The DA was performed for the 30 variables and also for the five most correlated variables with the first component (CP1), in order to classify new individuals. The PCA was able to identify sexual dimorphism in size, with both the 30 original variables as with the pre- selected variables, the latter optimized the reduction to two principal components. The DA was able to discriminate the two populations, both for 30 variables as for the five variables most correlated with the CP1. The functions with five variables can be used to classify other purebred dogs for sex, with an error of about 6.75%. PALAVRAS-CHAVE: Principal component analysis. Canis familiaris. Morphology. INTRODUCTION Sexual dimorphism can be considered a key evolutionary trait that can lead to important biological discoveries (AJAYI et al. 2012). Sexual dimorphism evolved in mammals to ensure greater reproductive success for individuals, especially males. Secondary sexual characteristics are attributes that facilitate mate choice, even if this feature is of little benefit or no apparent addition to its survival. These distinctions tend to be more pronounced in polygamous, with day habits and living in an open habitat (MCPHERSON; CHENOWETH, 2012) species. These features are considered as physical or behavioral attributes genetically transmitted to offspring. According to Polák e Frynta (2010), the artificial breeding during domestication implies pressures of different types of selection that occurs under natural conditions, and therefore it has different consequences on the body size. Even more, when breeders do not select by the size of the body, but by other morphological and behavioral characteristics. Therefore, in domestic breeds the favoring of larger males in sexual selection occurring in nature suffers a relaxation as a result of the domestication process. Lark, Chase and Sutter (2006) reported that the population of dogs is managed by man for longer than any other domesticated animal, providing enough time and opportunity to select for new phenotypic variations. The geographical isolation and selection for various tasks such as herding, guarding, hunting, and company have created specialized subtypes within the species. The species that exhibits the greatest morphological diversity among mammals is the domestic dog (Canis familiaris), and may exhibit sexual dimorphism in size of individuals (OSTRANDER; LINBLAD-TOH, 2006). In this species the variation in the skull and skeleton is superior to all other species of the Canidae family (WAYNE, 1986). Shultz et al. (2009) asserted that current studies of sexual dimorphism focus mainly on body size. Although some studies include other parts of the body (FAIRBAIRN, 2005), they use linear measurements and necessarily the size of the body is incorporated in these data. According to Dias e Barros (2009), morphometry studies only measurements of morphological characters of living beings like heights, widths and lengths, among others. Silva et. al. (2007) assigned that phenotypic characterization of a particular racial group may be aided by morphometric measurements. Received: 08/05/14 Accepted: 20/01/15 1476 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 Labrador Retriever is a British breed of dogs widespread in Brazil. It has several functions to men, such as assistance dog, rescue, guide, sniffer police dog and, primarily, a companion dog. According to Palika (2008) in the United States, the breed has achieved great popularity, which led to an increase of rearing corroborating to the lack of racial patterns and hence great variation in the breed. The objective of this study was to evaluate morphometric characteristics of the breed Labrador Retriever to establish descriptive biometric attributes that can differentiate these individuals sexually by means of multivariate statistical techniques. MATERIAL AND METHODS Seventy-four Labrador Retriever dogs - 27 males and 47 females - , obtained from commercial kennels located in the states of Rio de Janeiro, São Paulo and Minas Gerais were used in this experiment. The dogs used met the following requirements: be an adult over two years old, breeder or matrix animal, own record in the Brazilian Confederation of Cinofilia (CBKC) or the American Kennel Club (AKC), and no noticeable signs of pregnancy or lactation. There were no dogs in the same kennel with the same paternal and maternal ancestry, ie, no full siblings evaluated. The objective of this procedure was to maximize genetic variability besides a more reliable morphometric characterization of the Labrador breed in Brazil. Thirty quantitative biometric characteristics, continuously variable, relating to the morphology of the head, trunk and anterior and posterior limbs (Table 1) were measured. The biometric characteristics were measured using anthropometric ruler (± 0.1 cm), measuring tape (± 0.1cm) tape to measure circumferences (±0.1 cm) and a caliper (± 0.1 cm). Measurements of heights, member circumferences and body lengths were made with the standing animal, minimizing errors due to variations in the positions of animals. The animals were placed with the body weight distributed evenly among the four members whose axes (forearm to prior members and shuttle to subsequent) must remain upright. All measurements were taken in centimeters, because it is the unit of measure most used in Brazil. Table 1. Biometric characters measured Measured morphometric characteristics Total lenght (CTT) Height at withers(AGT) Width of the base of the tail (LBC) Length of skull ( CCB) Height of substernal emptiness (AVS) Width of the tip of the tail (LPC) Width of the skull (LCR) Height of the chest (APT) Circumference of the tail (PCD) Length of the muzzle (CCH) Height to the middle of torso (AMD) Length of the tail (CCD) Circumference of the muzzle (PCH) Width of the chest (LPT) Height of the elbow (ACV) Width of the bridge of the nose in base(LCB) Circumference of the chest (PTC) Circumference of forearm (PCN) Width of the bridge of the nose tip (LCP) Height of the rump (AGP) Circumference of metacarpus (PMT) Width of the ear (LOR) Width of the rump (LGP) Height of knee (AJH) Length of the ear (COR) Length of the rump (CGP) Circumference of knee (PJH) Circumference of the neck (PPC) Height of the insertion of the tail (AIC) Length of the trunk (CTR) To analyze the data obtained from the measurements and the existence of sexual dimorphism related to morphometry, the Principal Component Analysis (PCA) and discriminant analysis were used. The PCA is a multivariate statistical technique useful in many areas, to allow simplification in the data set. It summarizes the information originally contained in a group of variables v (characters) in a few components, which have the property of retaining the maximum variation originally available and independent. The technique of principal components is based only on the individual information of each sampling element 1477 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 (animal) without the need for data repetition. Thus, the procedure requires only one set of data values for each animal in relation to the biometric characters studied. Each principal component is a linear combination of all the original variables. The relative importance by which the main components are estimated is assessed by the total variance that they explain, being arranged so that the first component explains most of the total variation and the last component the lowest part (CRUZ; REGAZZI, 2008). According to Ferreira (2008), discriminant Analysis is a multivariate statistical technique used to discriminate and classify objects and/or individuals. This technique examines the separation of objects from the same population into two or more subpopulations and consists of two procedures. One of the procedures is discrimination or separation by means of functions of observed variables (discriminant functions) that are responsible or that explain the differences between subpopulations. The other procedure, allocation or classification, is the use of discriminant functions to classify new objects or individuals in one of the subpopulations functions. Statistica 6.0 (Statsoft Inc.) statistical software, described by Souza e Vicini (2005), for carrying out the multivariate analysis was used. RESULTS AND DISCUSSION Sexual dimorphism in mammals can occur in several ways, including anatomical and physiological characteristics that can serve for different purposes in the genesis, development and raising offspring. Another divergence area includes body size (MCPHERSON; CHENOWETH, 2012). To compare the body size between sexes in Labrador Retriever, the means and standard deviations were estimated from 30 morphometric variables (Table 2). Table 2. Means and standard deviations, in centimeters, of 30 measured variables. Variables Males Females Variables Males Females Total lenght (CTT) 119, 9 ± 4,88 111,20 ± 5,42 Circumference of the chest (PTC) 83,02 ± 7,06 79,58 ± 5,06 Length of skull (CCB) 24,15 ± 1,40 22,40 ± 0,95 Height of the rump (AGP) 54,34 ± 3,27 51,40 ± 3,11 Width of the skull (LCR) 13,12 ± 0,92 11,90 ± 0,74 Width of the rump (LGP) 19,78 ± 2,07 19,45 ± 2,38 Length of the muzzle (CCH) 9,40 ± 0,74 8,91 ± 0,58 Length of the rump (CGP) 11,53 ± 2,34 10,99 ± 1,16 Circumference of the muzzle (PCH) 29,19 ± 1,55 25,86 ± 1,37 Height of the insertion of the tail (AIC) 50,50 ± 3,76 46,89 ± 4,86 Width of the bridge of the nose in base (LCB) 5,08 ± 0,50 4,53 ± 0,45 Width of the base of the tail (LBC) 4,71 ± 0,62 4,16 ± 0,50 Width of the bridge of the nose tip (LCP) 3,99 ± 0,49 3,50 ± 0,28 Width of the tip of the tail (LPC) 1,12 ± 0,19 1,01 ± 0,17 Width of the ear (LOR) 10,54 ± 2,06 10,36 ± 1.13 Circumference of the tail (PCD) 18,43 ± 1,90 17,41 ± 1,54 Length of the ear (COR) 13,71 ± 1,18 12,87 ± 1,36 Length of the tail (CCD) 31,89 ± 2,56 29,31 ± 3,88 Circumference of the neck (PPC) 50,14 ± 3,45 45,89 ± 3,36 Height of the elbow (ACV) 28,30 ± 3,65 25,94 ± 2,81 Height at withers (AGT) 54,94 ± 3,54 51,26 ± 3,03 Circumference of forearm (PCN) 16,81 ± 1,37 15,51 ± 0,81 Height of substernal emptiness (AVS) 26,06 ± 1,68 22,96 ± 2,50 Circumference of metacarpus (PMT) 14,96 ± 0,82 13,99 ± 0,73 Height of the chest (APT) 35,77 ± 4,46 32,49 ± 4,30 Height of knee (AJH) 31,11 ± 1,61 28,36 ± 2,35 Height to the middle of torso (AMD) 53,38 ± 2,74 49,64 ± 3,10 Circumference of knee (PJH) 30,39 ± 2,46 28,03 ± 2,50 Width of the chest (LPT) 24,28 ± 6,17 21,92 ± 1,56 Length of the trunk (CTR) 65,22 ± 2,38 61,70 ± 2,75 1478 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 To demonstrate that quantitative measures of character inherent in skeletal morphology are able to show sexual dimorphism of Labrador Retriever dog, analysis of main components was performed in two ways. Initially we used the 30 measured variables and subsequently the analysis was performed with 23 pre-selected variables. Initially, ACP was processed using all 30 variables. Only the first three principal components were used to demonstrate sexual dimorphism, as they explained about 50% of the total variation (Table 3), summarizing much of the original information. Jackson (1993) reported that there is no optimum value for percentage of explanation by the first two or three axes. Generally, the larger the number of variables, the lower the percentage of explanation for these first components (MELO; HEPP, 2008). Table 3. Principal Components (CPi), eigenvalues (λi) and percentage of variance explained by components (% VCP) with 30 variables. CPi Λi %VCP %VCP (accumulated) 1(1) 10.78650 35.95500 35.9550 2(1) 2.28035 7.60117 43.5562 3(1) 1.98219 6.60730 50.1635 4 1.60749 5.35830 55.5218 5 1.54031 5.13436 60.6561 6 1.31823 4.39408 65.0502 7 1.16917 3.89724 68.9474 8 1.04406 3.48020 72.4276 9 1.02450 3.41502 75.8427 10 0.76762 2.55872 78.4014 11 0.71454 2.38178 80.7832 12 0.62375 2.07917 82.8623 13 0.58365 1.94550 84.8078 14 0.53411 1.78037 86.5882 15 0.46108 1.53693 88.1251 16 0.41593 1.38645 89.5116 17 0.39798 1.32661 90.8382 18 0.38255 1.27517 92.1134 19 0.35170 1.17235 93.2857 20 0.30007 1.00022 94.2859 21 0.28954 0.96515 95.2511 22 0.25117 0.83724 96.0883 23 0.23676 0.78921 96.8775 24 0.22374 0.74581 97.6233 25 0.20245 0.67485 98.2982 26 0.13687 0.45624 98.7544 27 0.12562 0.41875 99.1732 28 0.11450 0.38167 99.5548 29 0.07158 0.23860 99.7934 30 0.06197 0.20656 100.0000 (1) Three components needed to explain 50% of the total variation. Each principal component is a linear combination of all variables, and some of these variables are more representative within this component and others less (FERREIRA, 2008). According to Abreu et al. (1999), the variables strongly correlated with a particular main component have more importance for this component. The variables and their correlations with the first three principal components are presented (Table 4). The CP1 is the most representative component, accounting for approximately 36% of the total variance. The correlations of this component with the original variables show an inversely proportional relationship in view of the negative correlations nature. The CP1 is directly connected to the body size of the animals, because it is highly correlated with variables related to body size directly correlated with each other, such as neck perimeter, height to the middle of the torso, 1479 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 perimeter of the muzzle, height at withers and length of the trunk, these with correlation with approximately 75% or more with CP1 (Table 4). Table 4. Correlations between the original variables and the first three principal components. Variable CP1 CP2 CP3 Total lenght -0.71080 0.10929 0.18247 Length of skull -0.68788 0.18703 0.07873 Width of the skull -0.67838 0.12622 0.20279 Length of the muzzle -0.47981 0.12385 0.11567 Circumference of the muzzle(1) -0.81952 0.03103 0.02748 Width of the bridge of the nose in base -0.63965 0.32182 0.11598 Width of the bridge of the nose tip -0.68830 0.08562 0.32123 Width of the ear 0.07812 -0.28209 0.59751 Length of the ear -0.42056 -0.34168 0.11195 Circumference of the neck (1) -0.80831 -0.25828 0.07177 Height at withers (1) -0.76955 0.05845 -0.21046 Height of substernal emptiness -0.53980 0.42866 0.00605 Height of the chest -0.38969 0.20541 -0.38852 Height to the middle of torso (1) -0.80875 0.32297 -0.23832 Width of the chest -0.46460 -0.44388 0.14272 Circumference of the chest -0.65426 -0.28501 0.10512 Height of the rump -0.62097 0.25162 -0.23407 Width of the rump -0.16176 0.63172 0.02805 Length of the rump -0.34409 0.32327 0.12792 Height of the insertion of the tail -0.68002 0.01173 -0.04055 Width of the base of the tail -0.63101 -0.24733 0.33534 Width of the tip of the tail -0.29709 -0.39371 -0.41870 Circumference of the tail -0.60767 -0.50226 -0.00073 Length of the tail -0.30833 0.18206 0.72792 Height of the elbow -0.35783 -0.09189 -0.50267 Circumference of forearm -0.66771 -0.17618 -0.05659 Circumference of metacarpus -0.69061 -0.06699 -0.34981 Height of knee -0.67155 0.08935 0.02197 Circumference of knee -0.61003 -0.29146 -0.15047 Length of the trunk (1) -0.74765 -0.10305 -0.04809 (1) Five more correlated variables with the CP1. The principal component analysis was able to characterize sexual dimorphism. The CP1 combined with either one of the other two components was effective in trying to make the distinction of male and female Labrador Retriever (Figure 1). The use of three principal components to represent the sexual dimorphism is not required in order that such a distinction is displayed using only the first two components, which together explain approximately 43% of the total variation. This fact confirms the capacity of key components in reducing the dimension of the space of the original 30 variables to only two components. The discrepancy between the positions of females and males in relation to CP1 (y-axis) allows us to suggest a difference in size between males and females in Labrador Retriever (Figure 1). All variables are negatively correlated with the CP1, which establishes higher scores for females, because of the lower estimates for the measured characteristics. When placed on the Cartesian plane the scores of females provide superiority in positioning compared to the y axis (CP1), in analogy to the males, who have higher morphometric estimates and, consequently, lower scores on the CP1. 1480 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 Figure 1. Sexual dimorphism by the scores of principal components (CPi) 1 and 2 (a) and 1 and 3 (b), using 30 variables. Males (triangles) and females (circles). Sexual dimorphism in size may be related to male-male competition, ensuring the reproductive success of these animals, resulting in the selection of larger males (LANDE, 1980). According to Sutter et al. (2008), sexual dimorphism in size is present in most domestic dogs, being evidenced by the height at the withers in some breeds, in which males are generally larger than females. The pre-selection of variables aimed to better effect in the data reduction, and was performed by dispersing 30 variables of the first three principal components, which together represent 50% of the total variation. Wold et al. (1987) assigned that the principal component analysis estimates the correlation structure of the variables, in which the importance of a variable is determined by the size of their residual variance and this can be used for variable selection. Analyzing the graphs of correlation of the variables with the components (Figure 2), we selected the variable with distance of at least 0.6 of the origins of the axes in both graphs, since the variables with the highest correlation with a particular component are those that further get from the origin of the axis of the respective component. The pre-selection of the most correlated variables was performed considering that the PCA assumes that the variables are correlated (VICINI; SOUZA, 2005). Figure 2. Correlations of the 30 variables with the first three principal components. 01 = total length, 02 = length of skull, 03 = width of the skull, 04 = length of the muzzle, 05 = circumference of the muzzle, 06 = width of the bridge of the nose in base, 07 = width of the bridge of the nose tip, 08 = width of the ear, 09 = length of the ear, 10 = circumference of the neck, 11 = height at withers, 12 = height of substernal emptiness, 13 = Height of the chest, 14 = height to the middle of torso, 15 = width of the chest, 16 = circumference of the chest, 17 = height of the rump, 18 = width of the rump, 19 = length of the rump, 20 = height of the insertion of the tail, 21 = width of the base of the tail, 22 = width of the tip of the tail, 23 = circumference of the tail, 24 = length of the tail , 25 = height of the elbow, 26 = circumference of forearm, 27 = circumference of metacarpus, 28 = height of knee, 29 = circumference of knee, 30 = length of the trunk. 1481 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 Using the 23 pre-selected variables through the graphic of dispersion of the variables in main components, PCA noted eigenvalues, the percentage of variance explained by each component and the proportions of the accumulated variance (Table 5). Pre-selection was able to reduce the need for three main components to explain 50% of the total variation, as observed in the first procedure (30 variables) for only two main components. Table 5. Principal Components (CPi), eigenvalues (λi) and percentage of variance explained by components (% PCV) with 23 variables. CPi Λi %VCP %VCP (accumulated) 1(1) 9.864976 42.89120 42.8912 2(1) 1.845659 8.02460 50.9158 3 1.757124 7.63967 58.5555 4 1.308164 5.68767 64.2431 5 1.087780 4.72948 68.9726 6 0.887864 3.86028 72.8329 7 0.772471 3.35857 76.1915 8 0.755221 3.28357 79.4750 9 0.633274 2.75336 82.2284 10 0.558615 2.42876 84.6572 11 0.487826 2.12098 86.7781 12 0.431710 1.87700 88.6551 13 0.415605 1.80698 90.4621 14 0.350779 1.52513 91.9873 15 0.326584 1.41993 93.4072 16 0.294963 1.28245 94.6896 17 0.249527 1.08490 95.7745 18 0.235035 1.02189 96.7964 19 0.196062 0.85244 97.6489 20 0.188818 0.82095 98.4698 21 0.157825 0.68619 99.1560 22 0.116452 0.50631 99.6623 23 0.077666 0.33768 100.0000 (1) Two components needed to explain 50% of the variation The correlations of the variables with the components were obtained for the first two components, which together represented approximately 50% of the total variation (Table 6). For the second procedure (pre-selection of variables), the CP1 was the most representative, accounting for approximately 42% of the total variance. The correlations of this component with the original variables also show an inversely proportional relationship in view of the negative nature of correlations. The principal component analysis with pre- selected data was also able to characterize sexual dimorphism. The CP1 combined with CP2 explains just over 50% of the total variation and was efficient in trying to make the distinction between males and females in Labrador Retriever (Figure 3). This fact confirms the capacity of key components in reducing the dimension of the original space of 23 variables to only two components. The five variables most correlated with CP1 in both procedures (30 or 23 variables) are equivalent, which validates the relevance of the perimeter of the muzzle, neck perimeter, height to the middle of the torso, height at withers and trunk length in sexual dimorphism. González et al. (2011) asserted that morphometric analysis can be used to evaluate the homogeneity/heterogeneity of the indicated features by a particular group of animals and to establish a gradient of distances between them. Therefore, 30 variables measured for the discriminant analysis on the basis of sexual dimorphism of the Labrador Retriever dogs (two populations: males and females) were used. 1482 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 Table 6. Correlations between original variables and the two principal components. Variable CP1 CP2 Total lenght -0.720206 0.116921 Length of skull -0.682759 0.033423 Width of the skull -0.689467 0.071850 Length of the muzzle -0.813714 0.043017 Circumference of the muzzle(1) -0.640794 0.211766 Width of the bridge of the nose in base -0.696813 0.293845 Width of the bridge of the nose tip -0.797732 -0.020469 Width of the ear -0.767229 0.009110 Length of the ear -0.551771 0.140062 Circumference of the neck (1) -0.811652 -0.018611 Height at withers (1) -0.634086 0.031707 Height of substernal emptiness -0.638894 -0.063417 Height of the chest -0.175374 0.369789 Height to the middle of torso (1) -0.702526 -0.064674 Width of the chest -0.631211 0.109675 Circumference of the chest -0.601127 -0.144694 Height of the rump -0.321296 0.835299 Width of the rump -0.353307 -0.748526 Length of the rump -0.671565 -0.140012 Height of the insertion of the tail -0.697979 -0.424378 Width of the base of the tail -0.684659 -0.019991 Width of the tip of the tail -0.601451 -0.049785 Circumference of the tail(1) -0.753093 -0.180807 (1) Five more correlated variables with the CP1. Figure 3. Sexual dimorphism by the scores of principal components (CPi) 1 and 2 (a) and 1 and 3 (b), using the pre-selected variables (23 variables). Males (triangles) and females (circles). Cruz and Regazzi (2008) pointed that the classification is assigned replacing the measures observed in each variable for each of the individuals in the generated functions. The function that presents the highest value after replacing the measurements of each animal will be the population to which the individual belongs. The allocation or classification of new individuals may occur using the discriminant functions for each population (Male and Female) (Table 7). 1483 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 Table 7. Coefficients of the discriminant functions using 30 variables. Variable Female dogs (p=0.50) Male dogs (p=0.50) Total lenght 3.86 4.29 Length of skull 3.90 4.11 Width of the skull 3.45 3.51 Length of the muzzle 15.05 15.10 Circumference of the muzzle 14.16 16.28 Width of the bridge of the nose in base -31.43 -32.03 Width of the bridge of the nose tip 19.24 22.70 Width of the ear -4.58 -4.99 Length of the ear 14.83 15.53 Circumference of the neck -14.04 -14.60 Height at withers 4.70 4.80 Height of substernal emptiness 12.06 13.22 Height of the chest -4.25 -4.53 Height to the middle of torso -1.89 -2.11 Width of the chest -0.28 -0.28 Circumference of the chest 1.21 1.08 Height of the rump 2.32 2.18 Width of the rump 2.04 2.18 Length of the rump -2.26 -2.53 Height of the insertion of the tail -7.79 -8.41 Width of the base of the tail -35.53 -36.95 Width of the tip of the tail 153.04 169.23 Circumference of the tail 0.60 -0.12 Total lenght 4.36 4.59 Length of skull 6.47 6.99 Width of the skull 4.19 4.38 Length of the muzzle 35.05 36.88 Circumference of the muzzle -3.82 -3.90 Width of the bridge of the nose in base 0.33 0.60 Width of the bridge of the nose tip 13.16 14.05 Constant -1103.75 -1282.31 The Hotelling test ( =24.69) was significant (p<0.01), which indicates the difference between the mean vectors of both sexes. According to Mingoti (2005), it is necessary to assess the quality of the discriminant function built, being able to assess their significance and the probability of misclassification. The Hotelling test is a multivariate significance test that makes a comparison between the vectors of means among normal multivariate and independent populations. Ferreira (2008) stated that there are several methods for estimating the costs of misclassification, one being re-substitution method. Mingoti (2005) reported that the method of re- substitution or consistency analysis uses the information generated based on a contingency table, thus it was generated a contingency table with only two individuals misclassified (Table 8). Table 8. Contingency table with 30 variables. Populations % correct Female dogs (p=0.50) Male dogs (p=0.50) Total Female 95.7447 45 2(1) 47 Male 100 0(1) 27 27 Total 97.2973 45 29 74 (1) Animals misclassified. 1484 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 The same probability (p = 0.50) was assigned for the two populations (males and females), when processed discriminant analysis, since there is no difference in the importance of populations. The Chi-Square test ( =0.016) was not significant (p>0.05), therefore, pointed out the similarity between the observed frequencies (animals classified by discriminant functions for each sex), and the expected frequencies (animals measured for each sex). The consistency analysis technique was used to estimate the probability of misclassification of sexes. The estimate is given by the apparent error rate, which was 2.7% for discriminant analysis with 30 variables. The apparent error rate is underestimated, as these animals used to generate the classification function were used to estimate the error rate (FERREIRA, 2008). However, with a rate of 2.7% discriminant analysis was able to demonstrate sexual dimorphism for dogs from the Labrador Retriever breed using the 30 measured variables. The discriminant analysis was also processed using only five variables that were most correlated with CP1 in both previous procedures (PCA with 30 and 23 variables), since this component is the one with higher capacity of explanation of the total variation. This study was carried out in order to, not only get functions that discriminate males and females, but also enable the classification of other breed animals for sex, using a smaller number of characteristics (Table 9). Table 9. Coefficients of the discriminant functions using five correlated variables. Variable Female dogs (p=0.50) Male dogs (p=0.50) Circumference of the muzzle 8.969 10.505 Length of the muzzle -1.840 -1.881 Height to the middle of torso 1.103 1.161 Height at withers 1.226 1.158 Length of the trunk 8.171 8.531 Constant -385.319 -447.809 With the replacement of the measures in the two discriminant functions (males and females), the function that has the highest result (score) is that the individual should be allocated (PINTO et al., 2008). The Hotelling test ( =201.48) was significant (p<0.01), indicating the difference between the vectors of mean of both sexes. Thus, we generated a new contingency table (Table 10). Table 10. Contingency table with five correlated variables. Populations % correct Female dogs (p=0.50) Male dogs (p=0.50) Total Female 91.48936 43 4(1) 47 Male 96.29630 1(1) 26 27 Total 93.24324 44 30 74 (1) Animals misclassified. The Chi-Square test ( =0.022) was not significant (p>0.05), however, noted the similarity between the observed frequencies (animals classified by discriminant functions for each sex) and expected frequencies (animals measured for each sex). With a smaller number of variables, larger number of wrong classifications was observed. In the analysis of the consistency, apparent error rate (AER) equals to 6.75% was obtained. The pre-selection of variables by principal components analysis also allowed discrimination of animals in relation to sex, and still generally reduce the number of variables used for classification with an apparent error rate of approximately 6.75%. According to Purzyc, Kobrynczuk and Bojarski (2010), the correct classification of 70% of the animals by the discriminant analysis was satisfactory. The accuracy rate obtained in this study was approximately 93.25%, therefore, it is considered satisfactory. CONCLUSIONS When processed with 30 variables, principal components analysis was able to reduce the 1485 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 dimension to three components accounting for approximately 50% of the total variation, demonstrating aptitude in identifying the sexual dimorphism of the Labrador Retriever dogs. When processed with 23 variables (pre-selection), the PCA reduced to two principal components needed to explain approximately 50% of the total variation, corroborating its efficiency in indicating sexual dimorphism of dogs. Therefore, the pre-selection of variables was able to optimize the purpose of principal components analysis, given its superiority in reducing the dimensionality. In both studies the variables muzzle perimeter, neck perimeter, height to the middle of the torso, height at withers and trunk length showed greater importance in morphometric differentiation of the sexes. The discriminant analysis of Anderson was able to differentiate the two populations (males and females) for both the 30 original variables as for the five most correlated variables, selected by principal component analysis. The pre-selection of data performed with PCA was able to reduce the number of variables for use in discriminant analysis, without causing changes impact on the proportion of bad rating. However, the use of a greater number of variables made possible to reduce the probability of misclassification. Both the discriminant function obtained by using the 30 variables, as the one obtained when using the five most correlated variables, can be used in other papers/studies in an attempt to classify other dogs breed by sex. When contemplated five correlated variables, the error rate is approximately 6.75%. RESUMO: O cão doméstico (Canis familiaris) é a espécie de maior diversidade morfológica entre os mamíferos. Foram utilizados 74 animais da raça Retriever do Labrador, 27 machos e 47 fêmeas. Foram mensuradas 30 características biométricas quantitativas, relativas à morfologia. Objetivou-se avaliar as características morfométricas da raça Retriever do Labrador para estabelecer atributos biométricos descritivos que possam evidenciar o dimorfismo sexual por meio da análise de componentes principais (ACP) e da análise discriminante (AD). A ACP foi processada utilizando todas as variáveis e realizando uma pré-seleção das variáveis mais correlacionadas. A AD foi realizada para as 30 variáveis e também para as cinco variáveis mais correlacionadas com o primeiro componente (CP1), com intuito de classificar novos indivíduos. A ACP foi capaz de identificar o dimorfismo sexual de tamanho, tanto com as 30 variáveis originais quanto com as variáveis pré-selecionadas, esta última otimizou a redução para dois componentes principais. A AD foi capaz de discriminar as duas populações, tanto para 30 variáveis quanto para as cinco variáveis mais correlacionadas com o CP1. As funções com cinco variáveis podem ser utilizadas para classificar outros cães da raça quanto ao sexo, com um erro de aproximadamente 6,75%. PALAVRAS-CHAVE: Análise de componentes principais. Canis familiaris. Morfologia. REFERENCES ABREU, V. M. N.; SILVA, M. A.; CRUZ, C. D.; FIGUEIREDO, E. A. P.; ABREU, P. G. Capacidade de Combinação de Ovos de Linhagens de Matrizes de Frango de Corte, Usando Análise de Componentes Principais. Revista Brasileira de Zootecnia, Viçosa, v. 28, n. 5, p. 955-959, 1999. AJAYI, O. O.; YAKUBU, A.; JAYEOLA, O. O.; IMUMORIM, I. G.; TAKEET, M. I.; OZOJE, M. O.; IKEOBI, C. O. N.; PETERS, S. O. Multivariate analysis of sexual size dimorphism in local turkeys (Meleagris gallopavo) in Nigeria. Tropical Animal Health and Production, Boston, v. 44, n. 5, p. 1089-1095, 2012. http://dx.doi.org/10.1007/s11250-011-0044-6 CRUZ, C. D.; REGAZZI, A. J. Modelos biométricos aplicados ao melhoramento genético. 2. ed. Viçosa: UFV, 2008. 390p. DIAS, L. A. S.; BARROS, W. S. Biometria Experimental. 1. ed. Viçosa: UFV, 2009. 408p. FAIRBAIRN, D. J. Allometry for sexual size dimorphism: testing two hypotheses for Rensch’s Rule in the water strider Aquarius remigis. American Naturalist, Chicago, v. 166, p. 69–84, 2005. http://dx.doi.org/10.1086/444600 1486 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 FERREIRA, D. F. Estatística multivariada. 1.ed. Lavras: UFLA, 2008. 662p. GONZÁLEZ, A.; LUQUE, M.; RODERO, E.; GONZÁLEZ, C.; AGUILERA, R.; JIMÉNEZ, J.; SEPÚLVEDA, N.; BRAVO, S.; HERRERA, M. Use of Morphometric Variables for Differentiating Spanish Hound Breeds. International Journal of Morphology, Temuco, v. 29, p. 1248-1255, 2011. http://dx.doi.org/10.4067/S0717-95022011000400030 JACKSON, D. A. Stopping rules in principal component analysis: A comparison of heuristical and statistical approaches. Ecology, Washington, v. 74, p. 2204-2214, 1993. http://dx.doi.org/10.2307/1939574 LANDE, R. Sexual dimorphism, sexual selection, and adaptation in polygenic characters. Evolution, Lancaster, v. 34, p. 292–305, Mar. 1980. http://dx.doi.org/10.2307/2407393 LARK, K. G.; K. CHASE; N. B. SUTTER. Genetic architecture of the dog: sexual size dimorphism and functional morphology. Trends Genet, Cambridge, v. 22, p. 537–544, 2006. http://dx.doi.org/10.1016/j.tig.2006.08.009 MCPHERSON, F. J.; CHENOWETH P.J. Mammalian sexual dimorphism. Animal Reproduction Science, Amsterdam, v. 131, n. 3, p. 109-122, Abr. 2012. http://dx.doi.org/10.1016/j.anireprosci.2012.02.007 MELO, A. S.; HEPP L. U. Ferramentas Estatísticas para análise de dados provenientes de biomonitoramento. Oecologia Brasiliensis, Rio de Janeiro, v. 12, n. 3, p. 463-486, 2008. MINGOTI, S. A Análise de Dados Através de Métodos de Estatística Multivariada: Uma Abordagem Aplicada. 1 ed. Belo Horizonte: UFMG, 2005. 297p. OSTRANDER U. G. E. A.; LINDBLAD-TOH K. The dog and its genome. 1. ed. Cold Spring Harbor: Cold Spring Harbor Laboratory Press, 2006. 584p. PALIKA, L. Labrador Retriever. 2.ed. Hoboken: Wiley pub, 2008. 139p. PINTO, L. F. B.; ALMEIDA, F. Q.; QUIRINO, C. R.; AZEVEDO, P. C. N.; CABRAL, G. C.; SANTOS, E. M.; CORASSA, A. Evaluation of the sexual dimorphism in Mangalarga Marchador horses using discriminant analysis. Livestock Sicence, Amsterdam, v. 119, p. 161-166, Mar. 2008. http://dx.doi.org/10.1016/j.livsci.2008.03.014 POLÁK, J.; FRYNTA, D. Patterns of sexual size dimorphism in cattle breeds support Rensch’s rule. Evolutionary-Ecology Research, London, v.24, p.1255-1266. 2010. http://dx.doi.org/10.1007/s10682-010- 9354-9 PURZYC, H.; KOBRYNCZUK, F.; BOJARSKI, J. Sexual dimorphism in hucul horses using discriminant analysis. Animal, Cambridge, v. 5, n. 4, p. 506-511, Nov. 2010. http://dx.doi.org/10.1017/S1751731110002223 SOUZA, A. M.; VICINI, L. Análise Multivariada da Teoria à Prática. Santa Maria: UFSM, CCNE, 2005, 215p. SHULTZ, H.; POLLY, P. D.; KRIEGER, J. D.; GURALNICK, R. P. Diferential sexual dimorphism: size and shape in the cranium and pelvis of grey foxes (Urocyon). Biological Journal of the Linnean Society, London, v. 96, p. 339-353, 2009. http://dx.doi.org/10.1111/j.1095-8312.2008.01132.x SILVA, N. V.; FRAGA, A. B.; FILHO, J. T. A.; NETO, C. C. C.; SILVA, F. L.; COSTA, P. P. S.; JÚNIOR, W. B. L. Caracterização Morfométrica de Ovinos Deslanados Cabugi e Morada Nova. Revista Científica de Produção Animal, Fortaleza, v. 9, n. 1, 2007. 1487 Sexual dimorphism of Labrador… THULLER, M. A. O. et al. Biosci. J., Uberlândia, v. 31, n. 5, p. 1475-1487, Sept./Oct. 2015 SUTTER, B.N.; MOSHER, S.N.; GRAY, M.M.; OSTRANDER, A. E. Morphometrics within dog breeds are highly reproducible and dispute Rensch’s rule. Mammalian Genome Journal, New York, v. 19, p.713-723, Nov. 2008. http://dx.doi.org/10.1007/s00335-008-9153-6 WAYNE, R. K. Cranial morphology of domestic and wild canids: the influence of development on morphological change. Evolution, Lancaster, v. 40, p. 243–261, 1986. http://dx.doi.org/10.2307/2408805 WOLD, S.; ESBENSEN, K.; GELADI, P. Principal Component Analysis. Chemometrics and Intelligent Laboratory Systems Journal, Amsterdam, v. 2, p. 37–52, 1987. http://dx.doi.org/10.1016/0169- 7439(87)80084-9