Chief EditorK.M.L. Pathak
Print ISSN 0367-6722
Online ISSN 0976-0555
NAAS Rating 6.43
Impact Factor 0.5 (2023)
Morphological Standard of Murrah Buffalo of India based on Multivariate Analysis of the Phenotypic Traits
- Email firstname.lastname@example.org
First Online 02-03-2021|
Methods: Fourteen body biometric traits of 258 she buffaloes were recorded from the core and peripheral breeding tract of the breed. In order to separate the animals into two groups, data were subjected to cluster analysis. Principal component analysis (PCA) and discriminant analysis were used respectively to describe the body conformation and identify the combination of independent traits that distinguished Murrah from grade Murrah.
Result: Our study showed significant differences between the traits of Murrah and grade Murrah (p<0.05). The grade Murrah she buffaloes were smaller than the Murrah counterparts. PCA results showed that five factors explained 67.4% of total variance of studied morphometric traits in Murrah she buffaloes, whereas in Grade Murrah six components explained 72.7% of the total variance. The stepwise discriminant function derived with eight body biometric traits was able to discriminate Murrah and grade Murrah she buffaloes effectively. A total of 93.5% original grouped animals were classified correctly. The present work determines the morphometric standard of the Murrah buffalo through characterization and identification of latent factors of morphometric traits. Subsequent use of this information in marketing, breeding, management and conservation purposes can give significant productivity gains.
Buffaloes in the Indian sub-continent are an integral part of traditional agriculture and hold the greatest promise and potential for food production (Cockrill, 1994). India is the goldmine of farm livestock resources. The world’s 57 percent buffalo population contributing 68 per cent of world buffalo milk production with its 16 registered breeds describe the spectrum of buffalo portfolio the country holds. Among the 16 well-defined breeds, based on their phenotypic characteristics, production performances and geographical distribution (www.nbagr.res.in/), a world famous Murrah breed occupies prominent position being the highest milk producer and improver breed in India as well as abroad. The breed has spread over all parts of the country and is being bred either in pure form or is being used for grading up and substitution of local breeds (George et al., 1988; Sukla et al., 2006). The core home tract of Murrah is Haryana state. The native tract lies between 28o 15' and 30o North latitude and 75o 45' and 70o 80' East longitude (Sadana et al., 2006). The breed is aptly regarded as the ‘Black Gold’ or ‘Holstein-Friesian’ of the buffalo world in view of its global demand, higher milk production capacity, adaptation to different environmental conditions and feed conversion quality (Kumar et al., 2019).
Morphological characterization entails the description and documentation of the physical traits of a breed (Rege, 1992). It is the basis for differentiation of populations and/or breeds and provide support for breeding and conservation. Identification of true to breed Murrah animals is very important as the typical Murrah animals fetch a distinctively higher price in the market. Traditionally, the developmental agencies and farmers characterise Murrah on the basis of several morphometric characteristics (Nivsarkar et al., 2000; Sadana et al., 2006; http://pashudhanharyana.gov.in/). As the Murrah breed is being used for grading up of non-descript buffaloes, sometimes distinguishing Murrah from Murrah grades, becomes cumbersome, especially when their morphometric standards are lacking. This study characterizes the morphometric traits of Murrah buffalo from the core breeding tract and compares it with the surrounding graded Murrah using multivariate analysis of the traits. The findings would lend objectivity in decision making for marketing, genetic and breeding purposes of Murrah buffalo.
MATERIALS AND METHODS
Purposive sampling was done from core breeding tract having typical Murrah animals (Hisar, Bhiwani, Rohtak and Jhajjar districts of Haryana) and peripheral breeding tract having Murrah type or grade Murrah animals (Panipat, Karnal, Kurukshetra, Yamunanagar and Kaithal districts of Haryana). Data on Murrah bulls were recorded from semen stations of Haryana Livestock Development Board at Hisar, Gurugram and Jagadhari. The study was carried out at Indian Council of Agricultural Research-National Bureau of Animal Genetic Resources (ICAR-NBAGR), Karnal (Haryana), India during 2009-2013.
The measured morphometric traits were body length (BL), the distance from the point of the shoulder joint to the point of the pin bone; wither height (WH), the distance from the highest point of withers to the ground; chest girth (CG), the circumference of the chest just behind the elbow joint; paunch girth (PG), the circumference at paunch region just anterior to the hip joint; distance between hip bones (HBD); distance between pin bones (PBD); ear length (EL), distance from the point of attachment of ear to the tip of the ear; tail length (TL), distance from the root of tail droop to the tip of the tail excluding switch; tail whiteness (TW), whiteness of the tail from the distal end; horn length (HoL), external length of the horn; horn diameter (HoD), girth of the horn at base; horn curliness (HoC), chord distance between centre and base of the horn; head length (HL), distance between the tip of the skull to the tip of the nose and head width (HW), the distance between the left temple and the right temple.
Measurements (cm) of 258 she buffaloes and 93 Murrah bulls on 14 morphometric traits were taken with a measuring tape by same person to avoid inter-individual variations. For this data, we were interested in exploring whether the buffaloes from the core breeding tract and peripheral tract formed two distinct groups. Hierarchical cluster analysis using Ward distance metrics was performed to obtain two homogenous groups while discriminant function analysis was performed to ascertain the accuracy of the clustering. SAS and JMP Pro10 software of SAS were used in the analysis. Descriptive statistics for the morphometric traits of two groups were obtained. Means separation was done using Tukey’s HSD (p<0.05). Stepwise discriminant procedure was applied using PROC STEPDISC (SAS 9.3) to determine which morphological traits have more discriminant power than others. The relative importance of the morphometric variables in discriminating the Murrah from grade Murrah was assessed using the level of significance (p<0.05) and partial R2 values ≥0.01. The CANDISC procedure was used for calculating the Mahalanobis distances (Mahalanobis, 1936) of the morphological traits. The DISCRIM procedure was used to assign each individual animal to its group that included the eight morphometric traits. PCA with Varimax Kaiser Normalization method was used to describe the relationship between a large number of measured traits and a small number of unobserved factors. Kaiser criterion extracted the number of principal components (PCs). A variable was considered to be associated to a specific PC if the absolute value of its loading was ≥0.40.
RESULTS AND DISCUSSION
The phenotypic correlations among various morphometric traits of Murrah and Grade Murrah she buffaloes are presented in Table 2, while those of Murrah bulls are given in Table 3. Most of the correlations between the assessed variables in Murrah she buffalo were significant and ranged from 0.16 to 0.68 except tail and horn combinations with other traits. Correlations above 0.50 were obtained for CG, PG, BL and WH (Table 2), related to animal size. Similar pattern of correlation between the assessed traits in Murrah bulls was observed, however, significant correlation values were slightly higher and ranged from 0.21 to 0.73 (Table 3). Some of the traits (ear, tail and horn) exhibited negative correlations with the other traits which were mostly non-significant. Correlation matrix indicates the degree of relationship between the two variables. Our results indicated that in Murrah she buffalo the trait CG was positively associated with PG (0.68), WH (0.50), EL (0.19) and HW (0.21). Murrah bulls showed higher positive correlations of BL with WH (0.73) followed by CG with PG (0.71), BL (0.41) and WH (0.41). Vohra et al., (2015) showed in Gojri buffalo that WH had significantly higher correlation with BL (0.72), HL (0.50) and CG (0.50). Significantly positive associations observed among different traits suggests that indirect measure of a trait can be estimated through the estimations of other correlated trait.
Five PCs explained 67.4% of total variance of studied morphometric traits (Table 4) in Murrah she buffaloes, whereas in Grade Murrah six PCs explained 72.7% of the total variance (Table 5). The communalities found in Murrah she buffaloes ranged from 0.350 to 0.826. The communalities explain how much a particular trait contributes to explain the number of factors being considered (Morrison, 1976). The HoD trait showed lowest commonality, indicating that it contributed little towards explaining the total accumulated variation in the factors. In Murrah bulls, 4 PCs contributed to 62% of total variance of the studied traits, while the communalities varied from 0.352 to 0.793 (Table 6). In Murrah she buffaloes (Table 4), the first principal component contributed 18.6% of total variation and was represented by high loadings for PBD, CG and PG. The second component explained 18% of the total variance with high loading of HBD, BL, WH, EL and HL. These two components explained the general body conformation. The third component explained 12.7% of variance and showed high component loadings for horn characteristics (HoL, HoD and HoC). The fourth component explained 9.9% of the variance and showed high component loadings for tail characteristics (TL and TW). The fifth component explained 8.2% of variance and showed high component loading for HW. The traits that loaded highly onto first two PCs in Murrah loaded highly together onto PC1 in grade Murrah. The pattern of loading of traits in rest of the components in grade Murrah was different from Murrah (Table 4, 5). In Murrah bulls (Table 6), the first component explained 26% of the total variance with high loadings of seven traits (HBD, PBD, CG, PG, BL, WH and HW) explaining general body conformation of the bulls. Shahin et al., (1993) found three PCs explaining 88% of the variability while studying 13 body measurements in Egyptian buffalo bulls. Tolenkhomba et al., (2013) evaluating 18 biometric traits of bulls of local cattle of Manipur, India, reported that 69.77% of total variation was explained by six factors. Vohra et al., (2015), studying Gojri buffalo in India (based on 13 morphometric measurements), found four factors that explained 70.86% of the total variation. Yakubu et al., (2009) and Casanova et al., (2012) reported similar results in cattle. Murrah animals have distinguishable traits of horn, tail, body and forehead conformations and this knowledge is traditionally practiced by the local farmers to judge this breed. The factor loadings confirms this pattern. Our study also elucidates the structure of morphometric traits which were most useful for identification of Murrah animals. These factors could also be exploited in breeding and selection programmes to acquire highly coordinated animal bodies using fewer measurements (Yakubu et al., 2009).
Table 7 shows the results of stepwise discriminant analysis. Based on F-values and Wilk’s lambda, eight traits were significant. These variables in the data set were found to have potential discriminatory power. CG had more discriminant power than the others traits followed by HL as indicated by their higher R2 and F-values. In the canonical discriminant analysis, the canonical variable (CAN1) generated was significant (p<0.0001) (Table 8). Mahalanobis distance between two buffalo groups and per cent (%) of individual classified into source group are given in Table 9. The Mahalanobis distance between Murrah and grade Murrah she buffaloes was 6.92 and was highly significant (p<0.0001). This was substantiated by the classification result. In Murrah, 92% animals were correctly assigned while the corresponding figure in grade Murrah was 97%. A very high percentage of the correct assignment of the animals in their respective classes suggests that the Murrah breed is clearly distinguishable based on its morphological characteristics. Similar results of discriminant analysis using ten morphological traits were reported by Yakubu et al., (2010) in the study of morphometric differentiation of two Nigerian cattle breeds. The authors reported that 96.55% of Sokoto Gudali cattle were classified into their source genetic group, while 85.48% Bunaji counterparts were correctly assigned into their source population. The eight morphological variables (CG, HL, PG, HoC, HBD,TL, HoD, TW) obtained in the present study are more important and informative and could be used to assign the animals into Murrah and grade Murrah populations, thereby reducing the errors of selection in future breeding and selection programmes.
- Casanova, P.M., Parés i, Sniffed, I. and Villalba, D. (2012). Principal component analysis of cephalic morphology to classify some Pyrenean cattle. Animal Genetic Resource Information. 50: 59-64.
- Cockrill, W.R. (1994). Present and future of buffalo production in the world. Proceedings of the Fifth World Buffalo Congress, 27-30 June, Sao Paulo, Brazil.
- Dhillod, S., Kar, D., Patil, C.S., Sahu, S. and Singh, N. (2017). Study of the dairy characters of lactating Murrah buffaloes on the basis of body parts measurements. Veterinary World. 10(1): 17-21.
- George, M., Balaine, D.S., Vij, P.K., Kumar, S. (1988). Conservation and management of buffalo genetic resources of India, In: Buffalo Production and Health, [Nagarcenkar, R. (Ed.)], ICAR, New Delhi, pp. 31-37.
- Hardel, W.K. and Simar, L. (2007). Applied Multivariate Statistical Analysis, Second edn. Springer, Heidelberg.
- Johnson, R.A. and Wichern, D.W. (2007). Applied Multivariate Statistical Analysis, fifth edn. Pearson Education, Inc.
- Kumar, M., Dahiya, S.P., Ratwan, P., Kumar, S and Chitra, A. (2019). Status, constraints and future prospects of Murrah buffaloes in India. Indian Journal of Animal Sciences. 89(12): 1291-1302.
- Legaz, E., Cervantes, I., Pérez-Cabal, M.A., de la Fuente, L.F., Mártinez, R., Goyache, F. and Gutiérrez, J.P. (2011). Multivariate characterisation of morphological traits in Assaf (Assaf. E) sheep. Small Ruminant Research. 100: 122-130.
- Mishra, B.P., Singh, K.P., Chavan, D.B., Sadana, D.K., Kataria, R.S., Kathiravan, P. and Ahlawat, S.P.S. (2009). Characterization of Banni buffalo of Western India. Inimal Genetic Resource Information. 44: 77-86.
- Mahalanobis, P.C. (1936). On the generalized distance in statistics. Proceedings of the National Institute of Science India. 2(1): 49-56.
- Morrison, D.F. (1976). Multivariate statistical methods. 2nd edn. New York, McGraw-Hill Company.
- Nivsarkar, A.E., Vij, P.K. and Tantia, M.S. (2000). Animal Genetic Resources of India: Cattle and Buffalo. Directorate of Information and Publications of Agriculture, ICAR, New Delhi.
- Pundir, R.K., Vij, P.K., Singh, R.V. and Nivsarkar, A.E. (1996). Bhadawari buffaloes in India. Animal Genetic Resource Information. 20: 101-114.
- Rege, J.E. (1992). Animal genetic resources their characterization, conservation and utilization. Research Planning Workshop.
- Rezende, M.P.G. de, Ferraz, P.C, Carneiro, P.L.S. and Malhado, C.H.M. (2017). Phenotypic diversity in buffalo cows of the Jafarabadi, Murrah and Mediterranean breeds. Pesquisa Agropecuaria Brasileira. 52(8): 663-669.
- Statistical Analysis System. (2012). Version 9.3. SAS Institute Inc., Cary, NC, USA.
- Sadana, D.K., Kataria, R.S. and Mishra, B.P. (2006). Buffalo Genetic Resources of India: Murrah. Monograph 25, ICAR-National Bureau of Animal Genetic Resources, P.O. Box No.129, Karnal, Haryana, India.
- Shahin, K.A., Soliman, A.M. and Moukhtar. A.E. (1993). Sources of shared variability for the Egyptian buffalo body shape (conformation). Livestock Production Science. 36: 323-334.
- Sukla, S., Yadav, B.R. and Bhattacharya, T.K. (2006). Characterization of Indian Riverine Buffaloes by Microsatellite Markers. Asian-Australasian Journal of Animal Sciences. 19: 1556-1560.
- Tolenkhomba, T. C., Singh, N. S. and Konsam, D.S. (2013). Principal component analysis of body measurements of bulls of local cattle of Manipur, India. Indian Journal of Animal Science. 83: 281-284.
- Vohra, V., Niranjan, S.K., Mishra, A.K., Jamuna, V., Chopra, A., Sharma, N. and Jeong, D.K. (2015). Phenotypic characterization and multivariate analysis to explain body conformation in lesser known buffalo (Bubalus bubalis) from North India. Asian-Australasian Journal of Animal Sciences. 28: 311-317.
- Yadav, D.K., Jain, A., Kulkarni, V.S., Govindaiah, M.G., Aswathnarayan, T. and Sadana, D.K. (2013). Classification of four ovine breeds of southern peninsular zone of India: Morphometric study using classical discriminant function analysis. Springer Plus. 2: 29. https://doi.org/10.1186/2193-1801-2-29 PMID: 23449966.
- Yadav, D.K., Arora, R. and Jain, A. (2017). Classification and conservation priority of five Deccani sheep ecotypes of Maharashtra, India. PLoS One 12(9): e0184691. https://doi.org/10.1371/journal.pone.0184691.
- Yakubu, A., Ogah, D.M. and Idahor, K.O. (2009). Principal component of the morphostructural indices of White Fulani cattle. Trakia Journal of Sciences. 7: 67-73.
- Yakubu, A., Idahor, K.O., Haruna, H.S., Wheto, M. and Amusan, S. (2010). Multivariate analysis of phenotypic differentiation in bunaji and sokoto gudali cattle. Acta agriculturae Slovenica. 96(2): 76-80. doi:10.2478/v10014-010-0018-9.
- Yakubu, A., Salako, A.E. and Imumorin I.G. (2011). Comparative multivariate analysis of biometric traits of West African Dwarf and Red Sokoto goats. Tropical Animal Health Production. 43: 561-566.
- Yunusa, A.J., Salako, A.E. and Oladejo, O.A. (2013). Morphometric characterization of Nigerian indigenous sheep using multifactorial discriminant analysis. International Journal of Biodiversity Conservation. 5: 661-665.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.