Population diversity, admixture, and demographic trend of the Sumba Ongole cattle based on genomic data

Article information

Anim Biosci. 2024;37(4):591-599
Publication date (electronic) : 2023 November 2
doi : https://doi.org/10.5713/ab.23.0289
1Faculty of Animal Science, Universitas Gadjah Mada, Yogyakarta, 55281, Indonesia
2Research Center for Animal Husbandry, Research Organization for Agriculture and Food, National Research and Innovation Agency (BRIN), Bogor, 16911, Indonesia
3Research Center for Applied Zoology, Research Organization for Life Sciences and Environment, National Research and Innovation Agency (BRIN), Bogor, 16911, Indonesia
4Indonesian Research Institute for Animal Production, Indonesian Agency for Agricultural Research and Development, Ministry of Agriculture, Bogor, 16002, Indonesia
*Corresponding Author: Dyah Maharani, Tel: +62-274-513363, E-mail: d.maharani@ugm.ac.id
Received 2023 August 7; Revised 2023 September 18; Accepted 2023 October 20.



Sumba Ongole (SO) cattle are valuable breed due to their important role in the development of Indonesian cattle. Despite rapid advances in molecular technology, no genomic studies on SO cattle have been conducted to date. The aim of this study is to provide genomic profile related to the population diversity, admixture, and demographic trends of SO cattle.


Genomic information was gathered from 79 SO cattle using the Illumina Bovine SNP50 v3 Beadchip, and for comparative purposes, additional genotypes from 209 cattle populations worldwide were included. The expected and observed heterozygosity, inbreeding coefficient, pairwise fixation indices between-population, and Nei's genetic distance were examined. Multidimensional scaling, admixture, and treemix analyses were used to investigate the population structure. Based on linkage disequilibrium and effective population size calculations, the demographic trend was observed.


The findings indicated that the genetic diversity of SO cattle was similar to that of other indicine breeds. SO cattle were genetically related to indicines but not to taurines or Bali cattle. The study further confirmed the close relationship between SO, Ongole, and Nellore cattle. Additionally, a small portion of the Ongole mixture were identified dominant in the SO population at the moment. The study also discovered that SO and Bali cattle (Bos javanicus) could have been ancestors in the development of Ongole Grade cattle, which corresponds to the documented history of Ongolization. Our finding indicate that SO cattle have maintained stability and possess unique traits separate from their ancestors.


In conclusion, the genetic diversity of the SO cattle has been conserved as a result of the growing significance of the present demographic trend. Consistent endeavors are necessary to uphold the fitness of the breed.


Sumba Ongole (SO) cattle (Bos indicus) have historically played a significant role in Indonesian cattle development. As the name implies, this cattle breed was developed from Ongole cattle bred on Sumba Island - Indonesia [1]. The high cost of importation and the increasing scarcity of superior Ongole cattle in India have prompted the Ongole cattle breeding program since 1914 in Sumba Island [2]. Subsequently, SO cattle breed has spread throughout Indonesia, especially on the island of Java. This process was dubbed "Ongolization", where these cattle have been used for grading up program [3]. In 1936, a strict regulation was issued stating that SO cattle were being spread through the "Sumba Contract" system. Furthermore, male local cattle had to be castrated, while females were mated with SO bulls [4]. However, many superior SO bulls were later sold outside of the region, causing genetic erosion and quality degradation [1]. Since then, other cattle breeds have been imported to Sumba Island and crossed with SO cattle in order to improve their performance [5]. According to statistics, the beef cattle population in Sumba was expected to be around 43 thousand head by 2022 [6].

A thorough understanding of livestock diversity and population structure is required for maintaining genetic variation and is crucial for developing future breeding strategies, and the use of molecular technology for this purpose is gaining popularity [7]. Previously, the genetic diversity of SO cattle has been investigated through various analyses utilizing low-density microsatellite markers and mitochondrial DNA markers [8,9]. The revolution in molecular approach and the increasing availability of a large number genome-wide single nucleotide polymorphisms (SNPs) of animals have transformed them into the most important markers, and they are now intended to be used for genetic studies [10,11]. Genomic markers have been used to study diversity in Indonesian local cattle in the worldwide population context [12], however SO cattle had not yet been included. By conducting genomic studies on SO cattle, will supplement previous knowledges and hypotheses, as these cattle play a significant role in the formation of Ongole Grade, a valuable cattle breed with the widest distribution in Indonesia [3,13].

This study therefore aimed to investigate the genomic profile of SO cattle. Our current study represents the first comprehensive genome-wide analysis for this population by optimizing the genotypes obtained using the Bovine50SNP beadchip (Illumina Inc., San Diego, CA, USA). To better understand SO population ancestry referring to its historical context, we added genome datasets of pure Ongole cattle and other white-coated cattle from India previously published by Dixit et al [14]. The findings of our study will provide novel insights into the population diversity, admixture, and demographic trend of Indonesian SO cattle, which could be used as fundamental consideration in future breeding strategies to maintain breed fitness.


Ethical statement

All animal procedures were reviewed and approved by the Research Ethics Committee at Faculty of Veterinary Medicine, Universitas Gadjah Mada No.0120/EC-FKH/Eks./2022. Genotypes for the worldwide cattle breeds were obtained from previously published data [14,15], and therefore, no ethical approval was required.

Genomic data sources and quality controls

This study utilized genomic data of SO cattle in Sumba island, consisting of 79 animals with medium density SNPs. Animals were collected from three SO cattle populations on East Sumba island, which had similar conditions and climate. The genomic data was obtained from DNA which were extracted from the blood samples. DNA extraction was conducted at the animal genetics laboratory of the Indonesian Research Institute for Animal Production in Bogor, following the procedure outlined in the Wizard® Promega isolation kit method. DNA evaluations, i.e. concentration quantification and purity, were performed using Picogreen and NanoDrop devices, respectively. DNA samples were deemed to be of sufficient quality if the concentration was at least 20 ng/μL and the absorbance ratio at 260 and 280 nm was more than 1.8. Appropriate samples were then genotyped using the Illumina BovineSNP50 v3 Beadchip. These DNA evaluations and genotyping process were finished by the service of Macrogen Inc. in Republic of Korea. We added genotypes of 209 animals consisted of 11 cattle populations which were selected from previously published genome datasets by Dixit et al [14] and Decker et al [15] to support the analyses in our study, i.e. Bali, Ongole Grade, Ongole, Nellore, Hariana, Tharparkar, Brahman, Limousine, Simmental, Hereford, and Shorthorn. The last two populations, in particular, were only used in the admixture analysis to determine whether a true Brahman cross (BX) mixture exists in the SO cattle population. The detail of datasets information were summarized in Table 1.

Summary of the genomic datasets information

Quality control procedures were performed using PLINK v1.9 software [16] to ensure the reliability and accuracy of the data. SNP variants would be maintained if the Hardy–Weinberg Equilibrium (HWE) value was not <1×10−4. Any samples with a call rate of less than 90% were removed, as low call rates can indicate poor DNA quality or technical issues during genotyping. We also excluded any SNPs with a call rate of less than 90% or a minor allele frequency of less than 0.01, as these could indicate genotyping errors or non-polymorphic markers. Only autosomal SNPs with known genomic position based on UMD3.1 assembly and found in all populations were used in the subsequent analyses.

Genetic diversity, population structure, and admixture analysis

This study began by calculating several diversity parameters, including expected heterozygosity (He), observed heterozygosity (Ho), inbreeding coefficient (FIS), pairwise fixation indices between-population (FST), and Nei’s genetic distance to assess the genetic diversity of the SO and the comparators cattle populations. PLINK v1.9 [16] and R v4.1.0 [17] software were used to calculate these parameters. A genetic relationship matrix (GRM) of SO cattle was estimated using GCTA v1.93.2 [18]. A diagonal GRM value referring to the inbreeding of the animal itself and off-diagonal value referring to relationship between animals in the population.

This study also investigated the population structure using two different clustering methods: multi-dimensional scaling (MDS) and admixture analyses, which were performed using PLINK v1.9 [16] and fastSTRUCTURE [19], respectively. To determine the optimal number of subpopulations in the studied populations, the admixture analysis was ran using various estimated ancestral populations (referred to as K values). The cross-validation (CV) analysis was performed to select a minimum error value for K. Both clustering results were plotted in R v4.1.0 [17] to get understandable visualization. TreeMix v1.12 [20] was used to conduct extensive analysis of the relationships between cattle populations. A phylogenetic tree model was generated by positioning Bali cattle as root. The threepop command line was used to perform f3 statistical analyses introduced by Reich et al [21] to identify possible traces of introgression. As possible introgression events of populations B and C into population A, three-population statistical models (A, B, and C) with a significant negative value for both the f3 statistic and the Z-score were chosen.

Estimation of linkage disequilibrium and effective population size

The cattle population's demographic history was reflected in the estimated current and former effective population sizes (Ne). Sved and Hill [22] explained how past Ne was estimated from the linkage disequilibrium (LD) value bases on the following equation: E(r2) = 1/(1+4Nec). LD was annotated as r2 before Ne to measure the correlation of alleles at two loci and c was a recombination fraction in Morgan's unit. To estimate LD and Ne, the default PLINK v1.9 [16] and SNeP v1.0 [23] approaches were used, respectively. The --r2 command was used to get the LD value of SNP pairs, and the --ld- window-r2 command was used to get reports for all pairs. The Ne analysis was performed by using the standard input parameters in GONE approach [24] to infer a population's recent demographic history from SNP data of a small sample of individuals. Following the estimated times in the horizontal ordinate, the historical Ne values were plotted using R v4.1.0 [17].


Population structure and diversity

The genotyping process using the Illumina BovineSNP50 v3 Beadchip yielded a total of 53,218 SNPs on SO cattle. The observed genotyping call rate was 99.28%. Due to the quality control process, 28% of the variants were removed, leaving only 38,310 SNPs to be merged with the other cleaned genotypes of other cattle breeds (Table 1). Combined file consisted of 22,652 variants which have been encountered on 288 animals from 12 cattle breeds. Variability of the samples in SO population was proven by the negative off-diagonal variances in the GRM analysis, i.e. −0.0131. SO cattle showed commensurate level of both He and Ho with the other indicines, where the position is lower than taurines and higher than Bos javanicus (Bali cattle). Ongole Grade cattle had the highest He and Ho among indicines. Discovered in all populations, the Ho was slightly higher than the He. Likewise with the FIS of SO cattle, the value (−0.21) was almost the same as other indicines. Moreover, the FIS of cattle populations observed in this study were all negative.

An initial investigation into genetic diversity is summarized in the pairwise FST and Nei's genetic distance matrix (Table 2). This matrix clarified the distinctions and differences between SO and other cattle breeds. SO cattle were close to indicines but far from taurines and Bali cattle. Among indicines, Nellore cattle were the closest to Ongole from India (FST 0.0025; Nei 0.0061), while SO cattle were the farthest (FST 0.0074; Nei 0.0096). However, SO cattle had the shortest distance to Ongole cattle than to the others.

Pairwise FST (lower diagonal) and Nei’s genetic distances between populations (upper diagonal)

Population cluster and ancestries

Best illustration of the cattle genetic cluster could be explained by the MDS plot (Figure 1). A three-dimensional plot was used to clarify SO cattle position among indicines, taurines, and Bali cattle (Bos javanicus) coordinate. Interestingly, SO cattle populations tend to assemble in one colony but still near to indicines cluster, while Ongole, Nellore, Ongole Grade, Hariana, Tharparkar, and Brahman were close to one another. The admixture diagram (Figure 2) also shows the results of a more detailed cluster analysis. Based on the minimum CV error calculation (Supplementary Figure S1), K = 8 was chosen as the sensible model in this analysis. SO cattle could clearly be distinguished from other indicines, but the majority of their mixtures were similar to those found in Ongole cattle from India. Furthermore, the admixture with K = 3 was still displayed because it can be used to describe the sub species clusters, namely Bos javanicus (black color), Bos indicus (orange color), and Bos taurus (gray color).

Figure 1

Multi-dimensional scaling plot provide an overview of the position of Sumba Ongole cattle againts other indicines (Ongole, Nellore, Ongole Grade, Hariana, Tharparkar, Brahman), taurines (Simmental, Limousine), and Bos javanicus (Bali) clusters.

Figure 2

Admixture diagram depicts the estimated ancestry models composed on every individual in the populations. Mixture with K = 8 was chosen as the sensible model due to the minimum cross-validation error (0.523; Supplementary S1). The dominant mixture in Sumba Ongole cattle resembled a small portion of the Ongole cattle ancestry pattern.

The phylogenetic tree was developed specially for indicines and Bali cattle positioned as a root (Figure 3). Nellore and Ongole cattle from India were directly linked in the tree and SO cattle were closely related to those branch, while Ongole Grade cattle were in line after SO branch. The other indicines, i.e. Brahman, Tharparkar, and Hariana were in the different direction. Moreover, only one migration edge was observed, originating from the direction of Bali to Ongole Grade cattle. The f3 statistic was generated and the most significant results marked by negative value of both f3 statistics and Z-score were summarized in Table 3. This analysis only detect the possible introgression of Bali cattle (Bos javanicus) and indicines into Ongole Grade cattle. Here, SO cattle were the most likely candidate representing indicine ancestry (f3 statistic −0.0055; Z-score −11.4887). We looked more closely at the f3 statistics in which the Population A was SO cattle. This step was completed in order to identify the cattle breeds that have an impact to the SO population. Despite the fact that the values of both f3 statistic and Z-score were not negative, we select the combination of three-population statistical models with the lowest value. As a result, we identified that Ongole or Nellore (indicines originated from India) and Brahman cattle were the populations that suitable in the columns of Population B and C, respectively.

Figure 3

Phylogenetic tree of Sumba Ongole cattle among indicines with Bali cattle (Bos javanicus) as root. Sumba Ongole cattle were definitely associated with Ongole and Nellore cattle. Only one migration edge was observed, originating from Bali, and heading to Ongole Grade cattle.

The most significant f3 statistics and the top results related to Sumba Ongole

Demographic trends

An estimation of LD values based on r2 equation were done for all of SNP pairs and depicted following interallelic distances (Supplementary Figure S2). The observed LD values were high over a short allele distance and decayed steadily as the allele distance increased, while the average LD of adjacent markers (~100 kb) was 0.16±0.25. The historical Ne was calculated based on the LD value across the genome and was used to represent demographic changes in SO cattle population (Figure 4; Supplementary Table S1). Results from the SNeP software started the Ne estimation from 13th generation ago, while in GONE started from more recent time, i.e. 1st generation ago. Given in the similar generation time, both analyses were resulted in different values of Ne, e.g. Ne values in 13th generation ago were 922.48 and 342 for GONE and SNeP, respectively. However, both GONE and SNeP were detect similar pattern of dramatic Ne decline at almost in the simultaneously generation time estimates, i.e. GONE at 34 to 26 generations ago and SNeP at 54 to 13 generations ago. Moreover, according to the GONE estimation results, the values of Ne from the 25th generation ago to the current generation were stable with a slight increase.

Figure 4

The historical effective population size (Ne) estimation resulted from SNeP (A) and GONE (B) softwares. Both A and B reflect demographic changes over past generations of Sumba Ongole cattle. The similar patterns of dramatic Ne decline were identified, i.e., at 54 to 13 (A) and 34 to 26 (B) generations ago (Supplementary Table S1).


Historical events in the development of cattle in Indonesia have resulted in the formation of unique local cattle populations, one of which is the SO cattle. Despite historical records indicating that the SO cattle have an ancestral connection to the Ongole breed from India, our investigations have revealed genetic characteristics in the SO population that differentiate it from its parent population. Thus, the off-diagonal GRM value was negative which means that each individual in the population was not genetically related [18]. Moreover, the amount of variant that remained after the quality control process is very common and feasible in other cattle breeds (Table 1). However, we used only 22,652 SNPs found in all cattle breed populations studied, ensuring that all analysis parameters can be applied equally to each population.

In terms of heterozygosity, SO cattle have a comparable value with other indicines (Table 1). It has been consistently reported that the heterozygosity value of the Bos indicus sub species is in the midst of Bos taurus and Bos javanicus [12]. However, some argued that Bos indicus should have a high level of heterozygozity. This is most likely explained by the ascertainment bias, which causes the heterozygozity level of indicine to be underestimated [25]. In contrast to the previous statement, it should be noted that there are a limited number of appropriate bulls available in the area where indicine cattle frequently reared [26], which may account for the genetic erosion and low heterozygozity. Only favorite bulls are used in mating, and this has an impact on heterozygozity which decreases over time. Moreover, the Ho was reported higher than He indicating that there are polymorphisms in the population. It also suggests that the population has an excess of heterozygotes. The ratio of heterozygous alleles to homozygous alleles can be a factor of deviation in the HWE, causing the appearance of a negative value in FIS [27], as found in this study. Low FIS value suggested a beneficial effect of regulated genetic inbreeding in the examined populations.

The measurements of genetic divergence between SO and the comparator populations in this study were proven in our consecutive analyses. The FST and Nei's genetic distance analyses (Table 2) indicate that SO cattle are genetically closer to indicines than taurines and Bos javanicus, but among indicines, this cattle breed were not most similar (in term of allele frequencies) to Ongole cattle native to India. We noticed that Nellore cattle have approximately closest genetic relationship to Ongole cattle. In this study, Ongole Grade cattle were the population shown to have a close distance to both Ongole and Nellore. While in our previous hypothesis, Ongole Grade population in Kebumen regency were proposed as the maintained pure Ongole cattle breed in Indonesia [12]. Three-dimensional MDS plot clearly depicts the description of SO cattle position between the coordinates of indicines, taurines, and Bos javanicus. Even though the SO cattle were classified as zebu, they do not have a coordinate point that is close to the other indicines (Figure 1), whereas SO cattle were clearly on the same branch as Ongole and Nellore cattle in the phylogenetic tree (Figure 3).

The results of the ancestry estimation can explain such SO conditions. Only a small portion of the Ongole cattle ancestry pattern is now the dominant mixture in SO cattle (Figure 2). This is possible because the Ongole breeding, and purification program were carried out on the island of Sumba over a long period of time [28]. The dominant characters that SO cattle have today may be the result of adaptation process and selection patterns that have occurred within their environment. In the admixture analysis, we were unable to identify the BX mixture in the SO population. Therefore, the study was failed to verify the information provided by Sumadi and Siliwolu [5] that a cross between SO and BX cattle had occurred on Sumba island. This could be because BX cattle were only introduced outside of the SO breeding area. Further study must be carried out directly including BX cattle genome data, keeping in mind that we currently only use Brahman, Hereford, and Shorthorn cattle genotypes as a mixture representation of BX cattle.

Furthermore, we were able to provide genetical evidence that SO cattle were used to form Ongole Grade cattle. The results of the f3 statistical analysis significantly estimated that SO cattle were an introgressive population to Ongole Grade cattle (Table 3). This finding supports historical confirmation that SO cattle were used in Indonesia's Ongolization program [3]. In addition, the study has also detected evidence confirming the occurrence of introgression from Ongole cattle to SO cattle, but its value has since declined. These findings align with the results obtained from the admixture analysis. As a result, the study verifies that, despite the SO cattle are descendants of India's Ongole breed, they have exhibited stability and possess unique characteristics that differentiate them from their ancestors. The demographic trend analysis provides interesting points. Patterns of drastic Ne decline in SO population detected by both SNeP and GONE approaches (Figure 4) were common and found in most cattle populations. Breed formation, importation of animals, the isolation of population, merging, and eventually breed separation can all be attributed to the extraordinary changes in the demographic pattern [29]. If we look attentively, we can observe that the initial decline of Ne occurred during the formation of SO. While the stable curve condition with slight increase obtained from GONE analysis since the 25th generation ago up to the present may reflect SO cattle development efforts on the island of Sumba, to maintain genetic variation.


The genetic analyses show that the SO cattle population in Indonesia has distinct characteristics that set it apart from its Indian ancestor, the Ongole cattle. This study provides genetic evidence for historical records showing that Ongole Grade cattle originated through introgression from the SO population, and the SO cattle have undergone adaptation and selection processes over time. Despite the decline in effective population size, efforts have been made to maintain its genetic variation.



PS, HH, and DM were contributed to the study conception and design. Material preparation and blood collection were conducted by PS and HH. The DNA isolation was finished by BDPS, AARH, and SA. Genomic analyses were performed by PS. The manuscript was written by PS. Initial reviews of the manuscript were provided by DM, TSMW, and SB. All authors read and approved the final manuscript.


We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.


All process passed during data analysis, manuscript drafting, and publication were supported by the RTA program Universitas Gadjah Mada with Grant No. 5075/UN1.P.II/Dit-Lit/PT.01.01/2023.


We thank the animal owners and scientists who willingly shared their genetic information with the general public. Sample collection and genotyping of SO cattle were supported by the Indonesian Agency for Agricultural Research and Development through KP4S scheme.


The datasets analysed during the current study are available from the corresponding author on reasonable request.


Supplementary file is available from: https://doi.org/10.5713/ab.23.0289

Supplementary Figure S1. Plot of the cross-validation error estimation on the admixture analysis.


Supplementary Figure S2. Linkage disequilibrium decay of Sumba Ongole cattle.


Supplementary Table S1. Linkage disequilibrium (r2) and the historical effective population size (Ne) of Sumba Ongole cattle resulted in SNeP and GONE softwares



1. Hardjosubroto W. Alternative policies for the sustainable management of local beef cattle genetic resources within the national livestock breeding system. Wartazoa 2004;14:93–7.
2. Barwegen M. Golden Horns: The History of Livestock Farming on Java, 1850–2000 Ann Arbor, MI, USA: ProQuest LLC; 2020.
3. Astuti M. Potential and genetic resources diversity of Ongole Grade cattle. Wartazoa 2004;14:98–106.
4. Diwyanto K. Utilization of local resources and technological innovation to support the development of beef cattle in Indonesia. Pengembangan Inovasi Pertanian 2008;1:173–88.
5. Sumadi Siliwolu. Research on the genetic quality of Ongole and Brahman cattle in East Sumba Regency, East Nusa Tenggara. In : Proceedings of National Beef Cattle Workshop. Yogyakarta: Indonesian Center for Animal Research and Development; 2004. p. 31–41.
6. Statistics of Sumba Timur. Sumba Timur in Figures 2023 Waingapu: Indonesia: BPS-Statistics of Sumba Timur; 2023.
7. Lenstra JA, Groeneveld LF, Eding H, et al. Molecular tools and analytical approaches for the characterization of farm animal genetic diversity. Anim Genet 2012;43:483–502. https://doi.org/10.1111/j.1365-2052.2011.02309.x.
8. Agung PP, Anwar S, Wulandari AS, Sudiro A, Said S, Tappa B. The potency of Sumba Ongole (SO) cattle: A study of genetic characterization and carcass productivity. J Indones Trop Anim Agric 2015;40:71–8. https://doi.org/10.14710/jitaa.40.2.71-78.
9. Jakaria J, Musyaddad T, Rahayu S, Muladno M, Sumantri C. Diversity of D-loop mitochondrial DNA (mtDNA) sequence in Bali and Sumba Ongole cattle breeds. J Indones Trop Anim Agric 2019;44:335–45. https://doi.org/10.14710/jitaa.44.4.335-345.
10. Vignal A, Milan D, SanCristobal M, Eggen A. A review on SNP and other types of molecular markers and their use in animal genetics. Genet Sel Evol 2002;34:275. https://doi.org/10.1186/1297-9686-34-3-275.
11. Wilkinson S, Wiener P. Population genomics of animal domestication and breed development. In : Rajora OP, ed. Population genomics: concepts, approaches and applications Cham, Switzerland: Springer Nature; 2019. p. 709–53. https://doi.org/10.1007/13836_2017_8.
12. Sudrajad P, Subiharta S, Adinata Y, et al. An insight into the evolutionary history of Indonesian cattle assessed by whole genome data analysis. PLoS ONE 2020;15:e0241038. https://doi.org/10.1371/journal.pone.0241038.
13. Adinata Y, Noor RR, Priyanto R, Cyrilla L, Sudrajad P. Comparison of growth curve models for Ongole Grade cattle. Trop Anim Health Prod 2022;54:252. https://doi.org/10.1007/s11250-022-03254-z.
14. Dixit SP, Singh S, Ganguly I, et al. Genome-wide Runs of Homozygosity revealed selection signatures in Bos indicus. Front Genet 2020;11:92. https://doi.org/10.3389/fgene.2020.00092.
15. Decker JE, McKay SD, Rolf MM, et al. Worldwide patterns of ancestry, divergence, and admixture in domesticated cattle. PLoS Genet 2014;10:e1004254. https://doi.org/10.1371/journal.pgen.1004254.
16. Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 2015;4:s13742–015–0047–8. https://doi.org/10.1186/s13742-015-0047-8.
17. R Core Team. R: A language and environment for statistical computing R Foundation for Statistical Computing; Vienna, Austria: 2021. [cited 2023 Mar 10]. Available from: https://www.R-project.org.
18. Yang J, Lee SH, Goddard ME, Visscher PM. Genome-wide complex trait analysis (GCTA): methods, data analyses, and interpretations. In : Gondro C, van der Werf J, Hayes B, eds. Genome-wide association studies and genomic prediction London, UK: Springer Science+Business Media; 2013. p. 215–36.
19. Raj A, Stephens M, Pritchard JK. fastSTRUCTURE: Variational inference of population structure in large SNP data sets. Genetics 2014;197:573–89. https://doi.org/10.1534/genetics.114.164350.
20. Pickrell JK, Pritchard JK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet 2012;8:e1002967. https://doi.org/10.1038/npre.2012.6956.1.
21. Reich D, Thangaraj K, Patterson N, Price AL, Singh L. Reconstructing Indian population history. Nature 2009;461:489–94. https://doi.org/10.1038/nature08365.
22. Sved JA, Hill WG. One hundred years of linkage disequilibrium. Genet 2018;209:629–36. https://doi.org/10.1534/genetics.118.300642.
23. Barbato M, Orozco-terWengel P, Tapio M, Bruford MW. SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front Genet 2015;6:109. https://doi.org/10.3389/fgene.2015.00109.
24. Santiago E, Novo I, Pardiñas AF, Saura M, Wang J, Caballero A. Recent demographic history inferred by high-resolution analysis of linkage disequilibrium. Mol Biol Evol 2020;37:3642–53. https://doi.org/10.1093/molbev/msaa169.
25. Utsunomiya YT, Milanesi M, Fortes MRS, et al. Genomic clues of the evolutionary history of Bos indicus cattle. Anim Genet 2019;50:557–68. https://doi.org/10.1111/age.12836.
26. Hegde NG. Livestock development for sustainable livelihood of small farmers. Asian J Res Anim Vet Sci 2019;3:1–17.
27. Movahedin MR, Amirinia C, Noshary A, Mirhadi SA. Detection of genetic variation in sample of Iranian proofed Holstein cattle by using microsatellite marker. Afr J Biotechnol 2010;9:9042–5.
28. Watuwaya BK, Syamsu JA, Budiman , Useng D. Analysis of the potential development of beef cattle in East Sumba Regency, East Nusa Tenggara Province, Indonesia. IOP Conf Ser: Earth Environ Sci 2020;492:012153. https://doi.org/10.1088/1755-1315/492/1/012153.
29. Adepoju D. Estimating the effective population size of Swedish native cattle [thesis] Uppsala, Sweden: Swedish University of Agricultural Sciences; 2022.

Article information Continued

Figure 1

Multi-dimensional scaling plot provide an overview of the position of Sumba Ongole cattle againts other indicines (Ongole, Nellore, Ongole Grade, Hariana, Tharparkar, Brahman), taurines (Simmental, Limousine), and Bos javanicus (Bali) clusters.

Figure 2

Admixture diagram depicts the estimated ancestry models composed on every individual in the populations. Mixture with K = 8 was chosen as the sensible model due to the minimum cross-validation error (0.523; Supplementary S1). The dominant mixture in Sumba Ongole cattle resembled a small portion of the Ongole cattle ancestry pattern.

Figure 3

Phylogenetic tree of Sumba Ongole cattle among indicines with Bali cattle (Bos javanicus) as root. Sumba Ongole cattle were definitely associated with Ongole and Nellore cattle. Only one migration edge was observed, originating from Bali, and heading to Ongole Grade cattle.

Figure 4

The historical effective population size (Ne) estimation resulted from SNeP (A) and GONE (B) softwares. Both A and B reflect demographic changes over past generations of Sumba Ongole cattle. The similar patterns of dramatic Ne decline were identified, i.e., at 54 to 13 (A) and 34 to 26 (B) generations ago (Supplementary Table S1).

Table 1

Summary of the genomic datasets information

Cattle breed Sub species Samples Initial SNPs Cleaned SNPs He Ho FIS
Sumba Ongole B. indicus 79 53,218 38,310 0.31 0.39 −0.21
Ongole Grade1) B. indicus 20 42,885 30,210 0.33 0.43 −0.23
Ongole2) B. indicus 17 470,703 20,016 0.32 0.41 −0.20
Nellore1) B. indicus 20 42,885 30,210 0.31 0.40 −0.22
Hariana2) B. indicus 18 470,703 22,454 0.31 0.40 −0.21
Tharparkar2) B. indicus 17 470,703 22,454 0.31 0.40 −0.22
Brahman1) B. indicus 20 42,885 30,210 0.33 0.42 −0.22
Bali1) B. javanicus 20 42,885 30,210 0.14 0.21 −0.21
Limousine1) B. taurus 20 42,885 30,210 0.41 0.56 −0.32
Simmental1) B. taurus 20 42,885 30,210 0.41 0.56 −0.33
Hereford1) B. taurus 20 42,885 39,275 0.42 0.59 −0.35
Shorthorn1) B. taurus 17 42,885 39,275 0.39 0.56 −0.34

SNP, single nucleotide polymorphisms; He, expected heterozygosity; Ho, observed heterozygosity; FIS, inbreeding coefficients.


Decker et al [15].


Dixit et al [14].

Table 2

Pairwise FST (lower diagonal) and Nei’s genetic distances between populations (upper diagonal)

Cattle breed Sumba Ongole Ongole Grade Ongole Nellore Hariana Tharparkar Brahman Bali Limousine Simmental
Sumba Ongole 0 0.0102 0.0096 0.0107 0.0102 0.0109 0.0106 0.0492 0.0313 0.0331
Ongole Grade 0.0083 0 0.0076 0.0073 0.0077 0.0087 0.0076 0.0394 0.0273 0.0291
Ongole 0.0074 0.0040 0 0.0061 0.0083 0.0092 0.0088 0.0497 0.0312 0.0331
Nellore 0.0089 0.0041 0.0025 0 0.0084 0.0093 0.0085 0.0446 0.0304 0.0322
Hariana 0.0082 0.0043 0.0047 0.0051 0 0.0072 0.0074 0.0446 0.0317 0.0331
Tharparkar 0.0088 0.0052 0.0055 0.0059 0.0038 0 0.0083 0.0462 0.0322 0.0340
Brahman 0.0088 0.0042 0.0050 0.0052 0.0039 0.0046 0 0.0439 0.0284 0.0301
Bali 0.0495 0.0399 0.0512 0.0452 0.0453 0.0469 0.0443 0 0.0587 0.0601
Limousine 0.0303 0.0242 0.0278 0.0275 0.0284 0.0287 0.0253 0.0597 0 0.0145
Simmental 0.0323 0.0262 0.0299 0.0295 0.0301 0.0308 0.0272 0.0613 0.0099 0

Table 3

The most significant f3 statistics and the top results related to Sumba Ongole

Population A Population B Population C f3 statistics Standard error Z-score
A. Significant f3 statistics
 Ongole Grade Bali Sumba Ongole −0.0055 0.0005 −11.4887
 Ongole Grade Bali Ongole −0.0044 0.0008 −5.6271
 Ongole Grade Bali Tharparkar −0.0030 0.0008 −3.7589
 Ongole Grade Nellore Bali −0.0028 0.0007 −3.7935
 Ongole Grade Bali Brahman −0.0019 0.0007 −2.5775
 Ongole Grade Bali Hariana −0.0008 0.0008 −1.0848
B. Top observed results related to Sumba Ongole cattle
 Sumba Ongole Ongole Brahman 0.0104 0.0011 9.7436
 Sumba Ongole Nellore Brahman 0.0105 0.0011 9.6011