Donkey Populations Research Articles

Morphobiometric characterization of donkey resources in Katsina state

This study was conducted to assess the morphological characteristics of the donkey (Equus asinus) population in Katsina State, Nigeria. One hundred and three (103) donkeys (39 female and 64 male), randomly selected, were grouped into three age categories for this study. Data collected included: body weight and ten (10) body measurements. The statistical analysis of variance was applied for estimation the effects of sex and age on the studied parameters. Variance Inflation Factors (VIF), as a multicollinearity diagnostic tool of the independent variables, were incorporated in the regression models. The principal component analysis with a variance maximizing the orthogonal rotation was used to extract the components. The factor analysis identified four (4) principal components (PCs) of morphometric traits which explained about 70.392% of the total variance, with a PC1 accounted for 37.912%. The results showed that the body weight was correlated with HG, PG, BL, RH, SH, EL, TL and NC (ranging from 0.464 to 0.946; p<0.01) and body weight correlated to NL (r = 0.402; p<0.05) and HL (r = 0.375; P<0.05). The best predictors of body weight (BWT) in female and male revealed that the combination of heart girth, body length and rump height featured prominently in the prediction of BWT. Keywords: body weight, donkeys, morphological, prediction, regression

DWINDLING POPULATION OF INDIAN WORKING EQUINES: A WORRYING TREND

The equine species, which includes donkeys, mules, horses, and ponies, have made major contributions to human civilizations all throughout the world. The Indian livestock sector, after an initial decline in period from 2007 - 2012, staged a remarkable comeback during the period from 2012-2019 with an impressive overall growth rate of 4.8%. However, the total population of horses, ponies, mules and donkeys in India is 0.55 million as per the latest livestock census, a decrease by 51.9% over the previous livestock census (2012). The ‘decreasing utility’ due to increasing vehicles and urbanization has been considered as factors for the drastic dip in population. This paper, based on secondary data, aims to provide a concise summary of the current status, growth rate, distribution patterns, and factors leading to the decline in the population of working equines in India. Tabular presentations and graphical representations have been utilized for easy analysis and comparison of population trends across different years.

Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys.

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

Reproductive Disorders in Donkeys: Current Evidence and Update.

Reproductive disorders in donkeys present a significant challenge to their health and welfare, impacting their roles in agriculture, conservation, and companionship. With the development of large-scale donkey farming in recent years, reproductive disorders have become a limiting factor for the expansion of the donkey population. In general, donkeys suffer from a similar array of diseases like horses, but little is known about the specificities of donkey reproductive disorders. This review synthesizes current knowledge on the pathogenesis, distribution, presentation, diagnosis, treatment, and prognosis of a diverse array of reproductive disorders affecting donkeys. There are similar infectious and non-infectious causes for infertility and pregnancy loss in jennies compared with mares, but a difference in disease susceptibility does exist, which may be attributed to genetic influence, pathogen specificity, the environment, and reproductive management practices. Diagnostic and treatment plans need to be tailored towards the particularities of donkey reproductive disorders to increase donkey populations and to enhance the standard of care for this species. Enhancing reproductive health in donkeys not only ensures their sustainable use but also promotes their welfare and longevity in diverse human-animal interactions.

First identification and isolation of equine herpesvirus type 1 in aborted fetal lung tissues of donkeys.

Equine herpesvirus type 1 (EHV-1) is commonly associated with horse abortion. Currently, there are no reported cases of abortion resulting from EHV-1 infection in donkeys. This was the first survey-based study of Chinese donkeys. The presence of EHV-1 was identified by PCR. This survey was conducted in Chabuchar County, North Xinjiang, China, in 2020. A donkey EHV-1 strain (Chabuchar/2020) was successfully isolated in MDBK cells. Seventy-two of 100 donkey sera were able to neutralize the isolated EHV-1. Moreover, the ORF33 sequence of the donkey-origin EHV-1 Chabuchar/2020 strain showed high levels of similarity in both its nucleotide (99.7‒100%) and amino acid (99.5‒100%) sequences, with those of horse EHV-1 strains. EHV-1 Chabuchar/2020 showed significant consistency and was classified within cluster 1 of horse EHV-1 strains. Further, analysis of the expected ORF30 nucleotide sequence revealed that donkey EHV-1 strains contained guanine at position 2254, resulting in a change to aspartic acid at position 752 of the viral DNA polymerase. Therefore, these strains were classified as horse neuropathogenic strains. Lastly, a phylogenetic tree was constructed using the partial ORF68 nucleotide sequences, showing that the identified donkey EHV-1 strain and the EHV-1 strain found in aborted Yili horses in China comprised a novel independent VIII group. This study showed the first isolation and identification of EHV-1 as an etiological agent of abortions in donkeys. Further analysis of the ORF33, ORF30, and ORF68 sequences indicated that the donkey EHV-1 contained the neuropathogenic genotype of strains in the VIII group. It is thus important to be aware of EHV-1 infection in the donkey population, even though the virus has only been identified in donkey abortions in China.

Equine Infectious Anaemia: The Active Surveillance of an Entire Equid Population Reduces the Occurrence of the Infection

Equine infectious anaemia (EIA) is a life-long viral infection affecting equids, transmitted mechanically by biting flies and iatrogenic means. Despite its global distribution, active surveillance is limited, with passive clinical surveillance or control of specific equine sectors prevailing. In Italy, a national surveillance plan in horse, donkey, and mule populations has been established and includes mandatory passive and active surveillance through annual serological tests. During 2007–2010, the agar gel immunodiffusion (AGID) test served as both screening and confirmatory tests. Since 2011, a three-tier diagnostic pathway was introduced, utilizing the ELISA test for screening, AGID as the confirmatory test, and the immunoblot test for cases where ELISA was positive and AGID was negative. From a total equid population of 406,000 animals, 1,337,899 samples were analysed during 2007–2012, with 2,348 (0.18%) testing positive. EIA seroprevalence significantly decreased across all the species/hybrids during the study period. EIA occurrence was higher in mules (IRR = 48.90) and lower in donkeys (IRR = 0.56) compared to horses. The holding seroprevalence was 1.15%. Spatial analysis revealed clusters of infection in central Italy. These findings demonstrate that active systematic surveillance effectively reduces EIA prevalence in equid populations. Mules and working horses in wooded areas appeared to be at higher risk of infection and act as EIA reservoirs. Surveillance and control should be maintained and strengthened in these species/hybrids and in these areas to effectively control EIA. Passive surveillance alone is insufficient to eradicate the disease, and EIA remains a constant threat for the equine industry if active control is not implemented.

Analysis of the Genetic Relationship and Inbreeding Coefficient of the Hetian Qing Donkey through a Simplified Genome Sequencing Technology.

The Hetian Qing donkey is an excellent local donkey breed in Xinjiang. It is of great significance to accelerate breeding and the speed of breeding and rejuvenation, as well as to understand the genetic basis of the strategies and population. This study collected a total of 4 male donkeys and 28 female donkeys. It then obtained genotype data through Simplified Genomic Sequencing (GBS) technology for data analysis. The results detected a total of 55,399 SNP loci, and the genotype detection rate of individuals was ≥90%. A total of 45,557 SNP loci were identified through quality control, of which 95.5% were polymorphic. The average minimum allele frequency was 0.250. The average observed heterozygosity was 0.347. The average expected heterozygosity was 0.340. The average IBS (state homologous) genetic distance was 0.268. ROH: 49 (homozygous fragments), with 73.47% of the length between 1 and 5 Mb. The average per-strip ROH length was 1.75 Mb. The mean inbreeding coefficient was 0.003. The 32 Hetian green donkeys could be divided into six families. The number of individuals in each family is significant. To sum up, the Hetian Qing donkey population has low heterozygosity, few families, and large differences in the number of individuals in each family, which can easily cause a loss of genetic diversity. In the subsequent process of seed protection, seed selection should be conducted according to the divided pedigree to ensure the long-term protection of the genetic resources of Hetian green donkeys.

Morphometric characterization of udder and teat of Jennies

The present study was conducted in eight districts of Marathwada region, Maharashtra, India. A pre-survey was carried to find the maximum available population of donkeys in the respective area. Two hundred two (202) female donkeys were selected for the study. The data was collected based on morphometry of udder, teat, physiological status, pregnancy, lactation, and milk composition. Udder traits such as udder length, udder depth, udder width, udder circumference and teat traits such as teat length and teat diameter were found highly significant in lactating Jennies when compared with dry and pregnant females. The milk composition parameters were milk fat (MF), milk protein (MP), milk solid not fat (MSNF), total solid (TS), odour (O) and test (T): 1.31±0.08, 1±0.01, 6.93±0.03, 8.31±0.09% and 4593±120 per ml, respectively, for overall lactating jennet population.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to 9 different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified ROH islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights on the evolution of the genomic structure and origin of this important genetic resource.

Mitochondrial DNA diversity and maternal origins of Pakistani donkey.

Domestic donkey plays a key role as a draft animal in rural economy of Pakistan where its population is increasing every year. The complete mtDNA control region of forty randomly sampled donkeys was PCR- amplified and sequenced bi-directionally using specific primers. Distinct mtDNA haplotypes obtained in the current study (KY446001-KY446011) were subjected to haplotype (h) and nucleotide diversity (π) measures using DnaS as well as to phylogenetic, Network, and AMOVA analyses. There were a total 27 polymorphic sites present within 11 unique mtDNA haplotypes from the studied 40 animals from different regions. Neighbor-joining network and median-joining network both illustrated the splitting of all these haplotypes into two well-defined Nubian and Somali lineages, confirming African maternal origin of Pakistani domestic donkey. Diversity parameters h (0.967± 0.037) and π (0.02917± 0.00307) were found to reveal high levels of genetic diversity in Pakistani donkeys. AMOVA demonstrated only 1% of genetic differences between two mtDNA maternal lineages, pointing to lack of population substructure in Pakistani donkeys as is the case with worldwide domestic donkey population. Pakistani donkeys have African maternal origin and high levels of mtDNA diversity. High genetic diversity may be due to non-selective breeding and heteroplasmy. We herein provide the first report on mtDNA diversity of control region in Pakistani domestic donkey.

Feral Donkey Distribution and Ecological Impacts in a Hyper-Arid Region.

The feral donkey (Equus asinus L.) is an invasive species in Saudi Arabia and can cause severe damage to natural and cultural heritage. Over the last 30 years, feral donkeys have become a serious problem, as their abundance and geographic distribution has increased drastically. The impacts of feral donkeys are not well documented, and information about their abundance and distribution is lacking, certainly in Saudi Arabia, which hampers the implementation of effective management plans. Accordingly, we used the minimum population number approach (MPN) to determine the number of feral donkeys in this part of northwest Saudi Arabia. A total of 1135 feral donkeys were encountered in the region. The area around Khaybar harbors ~25% (n = 338) of the feral donkey population, whereas Tayma and AlGhrameel nature reserves were the least-inhabited sites (almost absent). The average population density of feral donkeys was estimated as 1.03 (0.19 SE) donkey/km2. We documented the negative ecological impact of feral donkeys on natural resources, which constituted overgrazing that resulted in habitat fragmentation and competition for resources with native species. We propose urgent actions to control the presence of feral donkeys in the region and suggest humane eradication as the most efficient and applicable to significantly reduce the negative impacts of feral donkeys.

Phenotypic characterization of donkey population in South Omo Zone, Southern Ethiopia

The study conducted in South Omo Zone, Southern Ethiopia with the aim of investigation within population diversification and characterization morphologically that helps to fill the gap of molecular characterization on this population. The data was collected from 500 adult donkeys from both sexes. Quantitative data was subjected to SAS GLM procedures by fitting districts and sex as the main effects. Qualitative data was subjected to a chi-square test with the district as the main effect. Color graph of donkey was done using Microsoft Excel, 2010. For both qualitative and quantitative data, the significance test was conducted at 5% of the level of error, and Tukey multiple range tests were used to separate the significance levels for the two types of data. CANDISC was used to calculate Mahalanobis distances, DISCRIM was used to cluster observations into predetermined groups, and STEPDISC was used to determine the quantitative characteristics that better differentiate populations. Roan coat color cover highest number compare with other coat color of donkey population. Quantitative traits of donkey has variation (P < 0.05) both in study areas and sex of donkeys. Overtly, except height at wither and height at the back Hammer donkey has mostly better metric value than the rest districts of the study areas. Moreover, CANDISC show variation on Hammer and Dasenech districts of donkey population. Furthermore, the longest (6.32) Mahalanobis distance observed in between Hammer and Dasenech donkey population. The Hammer and Dasenech donkey population is where the study fills in population variation the most. This can be because to management or genetics. Therefore, additional research might be required. Furthermore, morphometric measures show that donkey sex is similar, with the exception of heart girth circumference. This can be the result of poor selection, where superior male donkeys are sold for a higher price. Therefore, sound breeding programs should be used to reverse it.

Spatial and temporal distribution of donkey and mule population in India

The donkeys and mules are members of the Equidae family. Spiti, Halari and Kachchhi are registered breeds of donkeys in India. They are well acclimatized to the different agro-climatic conditions to which they belong. Donkeys and mules constitute about 37% of the total Indian equine population and play multifarious roles in human life, especially for the poorest of the poor section of society. They are mainly reared by underprivileged and weaker communities of the society and used for various purposes like transportation of goods, farm produces, pack animals at brick kilns, transportation of garbage, building materials and by nomadic pastoralists etc. The severe decreasing trend observed in Livestock Census elicited us to understand the present status of donkeys and mules. The total donkey and mule populations decreased by 61.23% and 57.10%, respectively in 2019 over the last census. This can be due to their decreased working utility, as alternate economic sources of mechanization are available. Further, motorized vehicles are replacing donkeys and mules in plains and hilly regions where there is a road network. Presently, donkeys and mules are more in rural area than urban area of the country, indicating their utility still exists among rural masses and rural transport depends on these equids. The density of mules is high in the hilly states of India, indicating their utility in difficult hilly terrains.

Means of Livelihood, Clean Environment to Women Empowerment: The Multi-Faceted Role of Donkeys

Despite the substantial contribution donkeys make to the livelihood of the world’s poorest populations, the existence of donkeys has received little notice worldwide. This article reviews the value of donkeys in a variety of sectors, including agriculture, construction industry, and mining, as well as their role in empowering women and achieving sustainable development goals. However, donkeys and mules are not given enough credit or attention in terms of developing strategies regarding their role in reducing poverty. There is a dearth of information and statistics on their impact across industries, the factors contributing to the donkey population dropping, the socioeconomic status of the dependent communities, and related animal and human welfare issues.

Identification and characterization of single nucleotide polymorphisms in DMRT3 gene in Indian horse (Equus caballus) and donkey (Equus asinus) populations

Equines’ ability in racing and riding as well as gaitedness have influenced the human civilization. Aim of this study was to identify and characterize the novel polymorphisms or SNPs in DMRT3 gene in Indian horse and donkey breeds. In this study, the DMRT3 gene was sequenced and characterized in 72 Indian horses’ and 33 Indian donkeys’ samples. One SNP (A > C) at 878 was found in studied horses while identical SNPs (A > C) at two different nucleotide positions i.e., 878 and 942 in DMRT3 gene (chromosome 23) were observed in studied Indian donkey breeds. Horses and donkeys both have a non-synonymous mutation (A > C) at nucleotide 878 (codon 61) that converts a Stop codon (TAG > TCG) to coding codon Serine, whereas donkeys have a synonymous mutation at nucleotide 942 (codon 82) that converts Serine (TCA > TCC) into Serine. A phylogenetic tree indicated that the DMRT3 gene was equally distributed among the equine breeds. Most of the donkey breeds have been shown high levels of genetic diversity while horse breeds and Halari donkey showed the least genetic diversity. Mutation in DMRT3 has a major impact on gaitedness in horses and is presented at a high frequency in gaited breeds and in horses breed for harness racing.

Genome-wide association study for numbers of vertebrae in Dezhou donkey population reveals new candidate genes

Numbers of vertebrae is an important economic trait associated with body size and meat productivity in animals. However, the genetic basis of vertebrae number in donkey remains to be well understood. The aim of this study was to identify candidate genes affecting the number of thoracic (TVn) and the number of lumbar vertebrae (LVn) in Dezhou donkey. A genome-wide association study was conducted using whole genome sequence data imputed from low-coverage genome sequencing. For TVn, we identified 38 genome-wide significant and 64 suggestive SNPs, which relate to 7 genes (NLGN1, DCC, SLC26A7, TOX, WNT7A, LOC123286078, and LOC123280142). For LVn, we identified 9 genome-wide significant and 38 suggestive SNPs, which relate to 8 genes (GABBR2, FBXO4, LOC123277146, LOC123277359, BMP7, B3GAT1, EML2, and LRP5). The genes involve in the Wnt and TGF-β signaling pathways and may play an important role in embryonic development or bone formation and could be good candidate genes for TVn and LVn.

A survey report on the donkey original breeding farms in China: Current aspects and future prospective.

The number of the large-scale donkey breeding farms in China has increased dramatically. However, information regarding the situation of a Chinese donkey population under large-scale donkey breeding farms is limited. This survey report was conducted using questionnaires online to investigate the current situation of the donkey original breeding farms in China, in terms of donkey stock, local breeds, reproductive parameters, growth and lactation performance, and future perspectives. China has developed the donkey reserve system based on national, provincial and non-governmental (self-own) donkey original breeding farms. In the present study, a total of 38 donkey original breeding farms concentrated in Northern of China were studied, and 52% of them keep their donkeys with a stocking density of 100-500 donkeys. China is rich in various local donkey breeds, and 16 local donkey breeds including large-sized, medium-sized and small-sized breeds were collected in our survey. Dezhou donkey are prevalent with a percentage of more than 57% of the total donkeys, while the Cullen donkeys belong to small-sized donkey breeds are scare. The reproductive efficiency and donkey productivity were different across donkey farms, indicating potential differences in management and breeding practices between different donkey original breeding farms. The artificial insemination has been performed in these donkey farms with an average proportion of 73%. Regarding the donkey productivity, the national and provincial donkey original breeding farms showed a higher birthweight and fat content in donkey milk than self-own farms. Furthermore, our results indicate that the donkey breeds with different body size also have important influence on the reproduction parameters and donkey productivity, with the large-sized donkeys had better performance compared to the small-sized donkeys. In summary, our survey provided valuable baseline information on the situation of donkey population dynamics in the donkey original breeding farms. However, further study is required in the future to investigate the factors such as donkey health care, management and nutrition during breeding, fattening and lactation that influence donkey productivity under large-scale farm systems.

The value of donkeys to livelihood provision in northern Ghana.

Increased demand for the supply of donkey hides for use in the Traditional Chinese Medicine e'jiao, is leading to a re-appraisal of donkeys' contributions to livelihoods across the world. This research aimed to understand the utilitarian value donkeys provide to poor small holder farmers, especially women, in their efforts to make a living in two rural communities in northern Ghana. Uniquely, children and donkey butchers were interviewed for the first time about their donkeys. A qualitative thematic analysis was undertaken of data disaggregated by sex, age and donkey-ownership. The majority of protocols were repeated during a second visit, ensuring comparative data between one wet, and one dry season. Donkeys are more important in people's lives than had previously been recognised and are highly valued by their owners for their help in reducing drudgery and the multi-functional services they offer. Hiring out donkeys to generate income is a secondary role for people who own donkeys, especially women. However, for financial and cultural reasons the way donkeys are kept results in the loss of a certain percentage of the animals to the donkey meat market, as well as the global hides trade. Increasing demand for donkey meat, coupled with increasing demand for donkeys for farming, is leading to donkey price inflation and theft of donkeys. This is putting pressure on the donkey population of neighbouring Burkina Faso and pricing resource-poor non-donkey owners out of the market. E'jiao has put the spotlight on the value of dead donkeys for the first time, especially to governments and middlemen. This study shows that the value of live donkeys to poor farming households is substantial. It attempts to understand and document this value thoroughly, should the majority of donkeys in West Africa be rounded up and slaughtered for the value of their meat and skin instead.

Donkey production systems and breeding practices in selected districts of South Omo Zone, southern Ethiopia

The study was conducted to understand the production system and breeding practices of Omo donkeys and provide baseline information for future production strategies. The data were collected through individual interviews, focus group discussions, and personal observation. A semi-structured questionnaire using 2012 Food and Agriculture Organisation guidelines was used to obtain respondent views. On the basis of donkey population, production potential, and road accessibility, 76 households from two districts were randomly selected and interviewed. Donkey growth performance and body size were given a higher priority in breeding selection. Most respondents selected their breeding donkeys based on growth performance, body size, and disease resistance. The criteria for culling were old age, body size, poor fertility, temperament, and colour. Variables that were given a higher priority in trait preference were disease resistance and pulling ability. Water shortage and disease were significant constraints to donkey production. Therefore, to address these constraints, the concerned bodies should invest in infrastructure to increase access to water and other services so as to support the livelihoods of donkey-dependent communities. Moreover, funding for animal healthcare and nutrition initiatives should be provided to ensure the health and well-being of animals.

ESMC: a statistical model to infer admixture events from individual genomics data

BackgroundInferring historical population admixture events yield essential insights in understanding a species demographic history. Methods are available to infer admixture events in demographic history with extant genetic data from multiple sources. Due to the deficiency in ancient population genetic data, there lacks a method for admixture inference from a single source. Pairwise Sequentially Markovian Coalescent (PSMC) estimates the historical effective population size from lineage genomes of a single individual, based on the distribution of the most recent common ancestor between the diploid’s alleles. However, PSMC does not infer the admixture event.ResultsHere, we proposed eSMC, an extended PSMC model for admixture inference from a single source. We evaluated our model’s performance on both in silico data and real data. We simulated population admixture events at an admixture time range from 5 kya to 100 kya (5 years/generation) with population admix ratio at 1:1, 2:1, 3:1, and 4:1, respectively. The root means the square error is pm 7.61 kya for all experiments. Then we implemented our method to infer the historical admixture events in human, donkey and goat populations. The estimated admixture time for both Han and Tibetan individuals range from 60 kya to 80 kya (25 years/generation), while the estimated admixture time for the domesticated donkeys and the goats ranged from 40 kya to 60 kya (8 years/generation) and 40 kya to 100 kya (6 years/generation), respectively. The estimated admixture times were concordance to the time that domestication occurred in human history.ConclusionOur eSMC effectively infers the time of the most recent admixture event in history from a single individual’s genomics data. The source code of eSMC is hosted at https://github.com/zachary-zzc/eSMC.