Indian Journal of Animal Research

  • Chief EditorK.M.L. Pathak

  • Print ISSN 0367-6722

  • Online ISSN 0976-0555

  • NAAS Rating 6.50

  • SJR 0.263

  • Impact Factor 0.4 (2024)

Frequency :
Monthly (January, February, March, April, May, June, July, August, September, October, November and December)
Indexing Services :
Science Citation Index Expanded, BIOSIS Preview, ISI Citation Index, Biological Abstracts, Scopus, AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus
Indian Journal of Animal Research, volume 57 issue 11 (november 2023) : 1561-1567

Factor and Discriminant Analyses in the Morphostructure of Batur and Wonosobo Sheep Breeds

Alek Ibrahim1,2,*, I Gede Suparta Budisatria1, Endang Baliarti1, Widya Pintaka Bayu Putra2, Endang Tri Margawati2, Bayu Andri Atmoko3, Wayan Tunas Artama4, Rini Widayanti4
1Department of Animal Production, Faculty of Animal Science, Universitas Gadjah Mada, Jl.Fauna No.3, Karangmalang, Yogyakarta-55281, Indonesia.
2Research Center for Applied Zoology, National Research and Innovation Agency, Jl.Raya Bogor Km.46, Cibinong-16911, West Java, Indonesia.
3Research Center for Animal Husbandry, National Research and Innovation Agency, Jl.Raya Jakarta-Bogor, Bogor-16915, Indonesia.
4Department of Biochemistry and Molecular Biology, Faculty of Veterinary Medicine, Universitas Gadjah Mada, Jl.Fauna No.2, Karangmalang, Yogyakarta-55281, Indonesia.
Cite article:- Ibrahim Alek, Budisatria Suparta Gede I, Baliarti Endang, Putra Bayu Pintaka Widya, Margawati Tri Endang, Atmoko Andri Bayu, Artama Tunas Wayan, Widayanti Rini (2023). Factor and Discriminant Analyses in the Morphostructure of Batur and Wonosobo Sheep Breeds . Indian Journal of Animal Research. 57(11): 1561-1567. doi: 10.18805/IJAR.BF-1455.
Background: Batur and Wonosobo sheep are two Indonesian indigenous sheep breeds from Central Java Province, Indonesia. This study aimed to describe and discriminate the morphostructure of Batur sheep and Wonosobo sheep with body measurements and body indices. 

Methods: Fourteen body measurements and nine body indices were measured from 37 Batur and 83 Wonosobo sheep. Two statistical methods of factor and discriminant analyses were computed with SPSS 25.0 program to describe the sheep’s morphostructure and characterization. 

Result: The factor analysis in body measurements was obtained three principal components (PCs) that explain about 76.58% (Wonosobo sheep) and 73.51% (Batur sheep) of the total variance in animals’ morphostructure. Meanwhile, the factor analysis in body indices was obtained on three PCs that explain about 79.16% (Batur sheep) and 95.57% (Wonosobo sheep) of the total variance in animals’ morphostructure. About 95.2% of Wonosobo sheep and 92.7% of Batur sheep can be classified with body measurements (RW, SW, EL, EW, RH and TL). Meanwhile, about 75.7% of Batur sheep and 80.7% of Wonosobo sheep can be classified with body indices (TI and CI). In conclusion, Batur and Wonosobo sheep can be accurately characterized by their body measurements.
Batur sheep and Wonosobo sheep are two Indonesian indigenous sheep breeds from Central Java Province, Indonesia. Batur sheep resulted from a crossbreed between Merino sheep and Thin-Tailed sheep that originated and spread widely in Batur District, Banjarnegara Regency. Wonosobo sheep is a new breed developed from crossing the local sheep and Texel sheep and spread widely in the Wonosobo Regency. Batur and Wonosobo sheep have been cultivated, maintained and raised by smallholder farmers for generations and are suitable for meat production (Ibrahim et al., 2021; Noviani et al., 2013; Hakim et al., 2019).
       
The characterization of livestock is the first approach to sustainable use of its animal genetic resources. Body dimensions have been used to indicate the breed, origin and relationship or shape and size of an individual as they give an idea of body conformation. Breed characterization in goats and sheep can be carried out based on body measurements and body indices (Markovic et al., 2019; Putra and Ilham, 2019). The morphometric characterization of livestock is essential for planning improvement, sustainable utilization, conservation strategies and breeding programs for a breed (FAO, 2012).   
       
The Batur and Wonosobo sheep were still considered non-descript sheep in the country and many things of morphostructures have not been studied. While some studies reported that the Batur and Wonosobo have similar ancestors, namely Thin-Tailed sheep (Hakim et al., 2019; Ibrahim et al., 2020b; Prayitno et al., 2008). Genetically, the Batur, Wonosobo and Thin-Tailed sheep have a close genetic relationship (Ibrahim et al., 2020b). Certainly, with similar ancestors and relatively close genetic distance, it is expected that the sheep provide performance and productivity so close different. Therefore, the present study was undertaken to describe and discriminate the morphostructure of Batur and Wonosobo sheep with body measurements and body indices.
 
 
Research site and animals
 
The samples were collected from 2019 to 2020 in Wonosobo and Banjarnegara Regencies, Central Java Province, Indonesia. However, this study was conducted in July-October 2021 in the Laboratory of Meat, Drought and Companion Animals, Universitas Gadjah Mada. The body measurements of Wonosobo sheep (22 rams and 61 ewes) and Batur sheep (13 rams and 24 ewes) were collected from smallholder farmers in Wonosobo and Banjarnegara Regencies, respectively. All animals in this study have two permanent pairs of incisors conditions with an average of 1.5 years of age (non-pregnancy stage). The phenotypic characteristic of Wonosobo sheep and Batur sheep and their breeding site were illustrated in Fig 1.
 

Fig 1: Wonosobo and Batur sheep physical characteristics and breeding site.


 
Animal’s measurements
 
The body measurements were taken from animals in a standing position with a raised head. Body measurements of animals were performed using a measuring stick and flexible measuring tape and taken based on FAO (2012) Guidelines. Fourteen body measurements of body length (BL), withers height (WH), rump height (RH), rump width (RW), chest circumference (CC), chest depth (CD), shoulder-width (SW), head length (HL), head width (HW), ear length (EL), ear width (EW), tail length (TL), tail width (TW) and tail circumference (TC) were measured in each animal as described in Fig 2. Moreover, nine body indices of length index (LI), thoracic index (TI), depth index (DI), height index (HI), thoracic development (TD), conformation index (CI), body index (BI), proportionality (Pr) and area index (AI) were computed in this study with mathematical formula (Boujenane, 2015; Khargharia et al., 2015) as follows:
 

Fig 2: The body measurements scheme of BL (a), WH (b), RH (c), RW (d), CC (e), CD (f), SW (g), HL (h), HW (i), EL (j), EW (k), TL (l), TW (m) and TC (n) in a sheep.


 
Data analysis
 
The statistics parameter of mean, standard deviation and Pearson’s coefficient of correlation (r) for body measurements and body indices were calculated with SPSS 25.0 computer program (IBM, USA). Therefore, factor analysis [Principal component analysis (PCA)] and discriminant analysis [Canonical discriminant analysis (CDA)] were computed using a similar computer program.
 
Animal morphostructure and Pearson’s coefficient of correlation
 
Nine body measurements of RW, CC, HL, HE, EL, EW, TL, TW and TC in Batur sheep were higher than in Wonosobo sheep (P<0.05) as presented in Table 1. Two body indices of HI and BI in Wonosobo sheep were higher than in Batur sheep (P<0.05), as presented in Table 2. Therefore, body indices of CI in Wonosobo sheep were lower than in Batur sheep. The Pearson’s coefficient of correlation (r) value between WH and RH in both sheep breeds was highest than the other correlations (0.97 in Wonosobo sheep and 0.93 in Batur sheep), as presented in Table 3. In addition, the highest r value among body indices showed in Pr-AI (-0.99) for both sheep breeds, as presented in Table 4.

Table 1: Means (±SD) of body measurements in Wonosobo and Batur sheep.



Table 2: Means (±SD) of body indices in Wonosobo and Batur sheep.



Table 3: Pearson’s coefficient of correlation (r) among body measurements in Wonosobo (above diagonal) and Batur (under diagonal) sheep.



Table 4: Pearson’s coefficient of correlation (r) among body indices in Wonosobo sheep (above diagonal) and Batur sheep (under diagonal).


       
The body measurements in the present study were higher than the previous study reported by Haryanti et al. (2015) and Haren et al. (2018) in Wonosobo and Batur sheep, respectively. According to the BI value, Wonosobo and Batur sheep included brevigline (BI<85) animals type similar to Ripollesa sheep (83.90±6.10) (Esquivelzeta et al., 2011). Handiwirawan et al. (2011) obtained LI value of 109.00±2.00 in Barbados Black Belly Cross (BC) sheep and lower than Wonosobo and Batur sheep. Therefore, the DI value in BC sheep was 52.0±0.9 and close to Batur sheep. The r-value of WH-RH in Uda (0.90) sheep (Yakubu and Akinyemi, 2010) was close to the present study.
 
Factor analysis
 
The factor analysis (PCA) of body measurements reveals three principal components (PCs) that explain 76.58% and 73.51% of the total variance in Wonosobo and Batur sheep morphostructures, respectively (Table 5). The PCA of body indices in sheep study reveals three PCs that explain 89.57% (Wonosobo) and 79.16% (Batur) of the total variance in sheep morphostructure (Table 6). Previous studies were revealed PCA of body measurements with a total variance of sheep morphostructure of about 59.00% (2PCs) in Ripollesa (Esquivelzeta et al., 2011); 61.53% (4PCs) in Rampur-Bushair (Sankhyan et al., 2018) and 96.65% (3PCs) in Pramenka (Markovic et al., 2019) sheep. Unfortunately, a study of the PCA of sheep’s body indices is limited. Previous studies reported that the PCA of body indices accounted for 86.84% (4 PCs) in Katjang does and 89.38% (4 PCs) in Pasundan cows (Putra et al., 2020; Putra and Ilham, 2019), showing higher than in sheep study.

Table 5: Rotated component matrix, eigenvalues, total variance, cumulative, communalities, Kaiser-Meyer-Olkin (KMO) measure adequacy, ad Bartlett’s test of sphericity in the body measurements of Wonosobo and Batur sheep.



Table 6: Rotated component matrix, eigenvalues, total variance, cumulative, communalities, Kaiser-Meyer-Olkin (KMO) measure adequacy and Bartlett’s test of sphericity in the body indices of Wonosobo and Batur sheep.


 
Discriminant analysis
 
The discriminant analysis (CDA) reveals six body measurements (RW, SW, EL, EW, RH and TL) and two body indices (TI and CI) that were selected as the discriminant variables (Table 7). According to the CDA, about 95.2% (Wonosobo) and 97.3% (Batur) of sheep can be classified into their original breeds group based on body measurements (Table 8). Meanwhile, about 80.7% (Wonosobo) and 75.7% (Batur) of sheep’s study can be classified into their original breeds based on body indices. The discriminant plots of body measurements and body indices were illustrated in Fig 3 and 4, respectively.

Table 7: Factor selected by stepwise discriminant analysis to characterize Wonosobo and Batur sheep.



Table 8: Percentage (%) of individual classification per breed based on discriminant analysis.



Fig 3: The discriminant plots of body measurements in Wonosobo and Batur sheep.



Fig 4: The discriminant plots of body indices in Wonosobo and Batur sheep.


       
The CDA showed that body measurements were more accurate for discriminating Wonosobo and Batur sheep than body indices. Ofori et al. (2021) obtained the rc value of 0.62 (high) in the first canonical factor of body measurement in West African Dwarf goats and lower than in the present study. Thus, (Asamoah-Boaheng and Sam, 2016) obtained the rc value of 0.89 (very high) in the first canonical factor of body measurements in sheep and similar to the present study. More than 90% of the morphostructure of sheep study were successfully classified into their original breeds group with six selected discriminant variables of body measurements. A previous study reported that the body measurements could classify many sheep breeds into their original breeds group (Asamoah-Boaheng and Sam, 2016). Moreover, 58% of Uda rams and 70.8% of Uda ewes can be classified into their original sex group based on two discriminant variables of HW and HL measurements (Yakubu and Akinyemi, 2010). The difference in body measurements, genetics (breed), location (geographical area), agro-climatic conditions and management system of animals can be caused by different results compared to previous studies.
 
 
In conclusion, about 95.2% of Wonosobo sheep and 92.7% of Batur sheep can be classified with body measurements (RW, SW, EL, EW, RH and TL). Meanwhile, about 80.7% of Wonosobo sheep and 75.7% of Batur sheep can be classified with body indices (TI and CI). Wonosobo sheep and Batur sheep can be accurately characterized by their body measurements.
 
The authors are thankful to the Direktorat Penelitian dan Tim Peningkatan Reputasi UGM menuju World Class University-Kantor Jaminan Mutu UGM, for funding this study with the Post-Doctoral Program (grant no. 6144/UN1.P.III/DIT-LIT/PT/2021). The authors also thank the Faculty of Animal Science Universitas Gadjah Mada, the National Research and Innovation Agency, the Department of Agriculture, Fisheries and Animal Husbandry of Banjarnegara Regency and the Department of Food, Agriculture and Fisheries of Wonosobo Regency for the support and thank all farmers.
 
None.
 

  1. Asamoah-Boaheng, M. and Sam, E.K. (2016). Morphological characterization of breeds of sheep: A discriminant analysis approach. Springer Plus. 5: 1-12. 

  2. Boujenane, I. (2015). Multivariate characterisation of Oulmes Zaer and Tidili cattle using the morphological traits. Iranian Journal of Applied Animal Science. 5: 293-299.

  3. Esquivelzeta, C., Fina, M., Bach, R., Madruga, C., Caja, G., Casellas, J., Piedrafita, J. (2011). Morphological analysis and subpopulation characterization of Ripollesa sheep breed. Animal Genetic Resources. 49: 9-17. 

  4. FAO (2012). Phenotypic Characterization of Animal Genetic Resources. FAO Animal Production and Health Guidelines (Vol. 11). Food and Agriculture Organization of the United Nations, Rome.

  5. Hakim, F.R., Arifin, M., Rianto, E. (2019). Growth Pattern and Productivity of Female Wonosobo Sheep in Wonosobo District, Central Java Province, Indonesia. IOP Conference Series: Earth and Environmental Science. 247: 012044. 

  6. Handiwirawan, E., Noor, R.R., Sumantri, C., Subandriyo. (2011). The differentiation of sheep breed based on the body measurements. Journal of the Indonesian Tropical Animal Agriculture. 36: 1-8.

  7. Haren, H.I.H., Saleh, D.M., Sumaryadi, M.Y., Purwantini, D., Prayitno. (2018). Morphometric Traits and Body Weight at Different Age of Batur Sheep in Banjarnegara. Prosiding Seminar Teknologi dan Agribisnis Peternakan VI. pp. 186-19.

  8. Haryanti, Y., Kurnianto, E., Lestari, C.M.S. (2015). Estimation of body weight using body measurements on Wonosobo sheep. Jurnal Sain Peternakan Indonesia. 10: 1-6. 

  9. Ibrahim, A., Artama, W.T., Budisatria, I.G.S., Yuniawan, R., Atmoko, B.A., Widayanti, R. (2021). Regression model analysis for prediction of body weight from body measurements in female Batur sheep of Banjarnegara District, Indonesia. Biodiversitas. 22: 2723-2730. 

  10. Ibrahim, A., Budisatria, I.G.S., Widayanti, R., Artama, W.T. (2020a). The genetic profiles and maternal origin of local sheep breeds on Java Island (Indonesia) based on complete mitochondrial DNA D-loop sequences. Veterinary World. 13: 2625-2634. 

  11. Ibrahim, A., Budisatria, I.G.S., Widayanti, R., Atmoko, B.A., Yuniawan, R., Artama, W.T. (2020b). On-farm body measurements and evaluation of Batur Sheep on different age and sex in Banjarnegara Regency, Indonesia. Advances in Animal and Veterinary Sciences. 8: 1028-1033. 

  12. Khargharia, G., Kadirvel, G., Kumar, S., Doley, S., Bharti, P.K., Das, M. (2015). Principal component analysis of morphological traits of Assam Hill goat in Eastern Himalayan India. The Journal of Animal and Plant Sciences. 25: 1251-1258.

  13. Markovic, B., Dovè, P., Markoviæ, M., Radonjiæ, D., Adakaliæ, M., Simèiè, M. (2019). Differentiation of some Pramenka sheep breeds based on morphometric characteristics. Archives Animal Breeding. 62: 393-402. 

  14. Noviani, F., Sutopo, Kurnianto, E. (2013). Hubungan genetik antara domba Wonosobo (Dombos), domba Ekor Tipis (DET) dan domba Batur (Dombat) melalui analisis polimorfisme protein darah. Sains Malaysiana. 11: 1-9.

  15. Ofori, S.A., Hagan, J.K., Kyei, F. (2021). Morphometric characterization and differentiation of West African Dwarf goat populations in Ghana. Tropical Animal Health and Production. 53: 1-14. 

  16. Prayitno, Hartatik, T., Pratiwi, R., Artama, W.T. (2008). Genetic relatedness between Batur, Merino and local sheep based on random amplyfied polymorphism DNA marker. Animal Production. 13: 30-38.

  17. Putra, W.P.B. and Ilham, F. (2019). Principal component analysis of body measurements and body indices and their correlation with body weight in Katjang does of Indonesia. Journal of Dairy, Veterinary and Animal Research. 8: 124- 134. 

  18. Putra, W.P.B., Said, S., Arifin, J. (2020). Principal component analysis (PCA) of body measurements and body indices in the Pasundan cows. Black Sea Journal of Agriculture. 3: 49-55.

  19. Sankhyan, V., Thakur, Y.P., Katoch, S., Dogra, P.K., Thakur, R. (2018). Morphological structuring using principal component analysis of Rampur-Bushair sheep under transhumance production in western Himalayan region, India. Indian Journal of Animal Research. 52: 917-922. 

  20. Yakubu, A. and Akinyemi, M.O. (2010). An evaluation of sexual size dimorphism in Uda sheep using multifactorial discriminant analysis. Acta Agriculturae Scand Section A. 60: 74-78.

Editorial Board

View all (0)