Genome-wide single-nucleotide polymorphism data and mitochondrial hypervariable region 1 nucleotide sequence reveal the origin of the Akhal-Teke horse

Article information

Anim Biosci. 2023;36(10):1499-1507

Publication date (electronic) : 2023 May 4

doi : https://doi.org/10.5713/ab.23.0044

Zhoucairang Kang¹, Jinping Shi¹, Ting Liu¹, Yong Zhang², Quanwei Zhang³, Zhe Liu¹, Jianfu Wang¹, Shuru Cheng^1,

¹College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China

²College of Veterinary Medicine, Gansu Agricultural University, Lanzhou 730070, China

³College of Life Science and Biotechnology, Gansu Agricultural University, Lanzhou 730070, China

^*Corresponding Author: Shuru Cheng, Tel: +86-0931-7632509, E-mail: chengsr@gsau.edu.cn

Received 2023 February 8; Revised 2023 March 13; Accepted 2023 April 10.

Abstract

INTRODUCTION

The Akhal-Teke horse is estimated to be at least 3,000 years old and is one of the oldest horse breeds in the world. Native to the southern region of Turkmenistan, the Akhal-Teke horse is adapted to harsh climatic conditions [1] and is famous for its beautiful appearance, familiarity with people, and shining coat color. The Akhal-Teke horse is one of the finest horse breeds in the world. Today, the Akhal-Teke horse population is slowly declining and is endangered. Currently, the Akhal-Teke is found mainly in Turkmenistan [2].

The Akhal-Teke is the best of a number of excellent horse breeds. The origin of this horse breed is highly controversial. Several hypotheses regarding the origin of the Akhal-Teke horse have been proposed based on archaeological and documentary evidence. The most plausible hypothesis is that the Akhal-Teke horse originated in Central Asia, in the area between the Caspian coast and the Fergana region [3]. It is thought that the earliest domestication of horses occurred in the Eurasian steppe around 3500 BCE, from where they spread across the Eurasian continent [4]. A study of DNA microsatellite genotyping of 2,024 Akhal-Teke horses from around the world showed that this horse breed has a high level of genetic diversity with abundant rare alleles at microsatellite DNA loci. Further analysis of the rare alleles showed that the Akhal-Teke horse is very similar to Tuvinia and Khakassia horses [5]. These results support the hypothesis of a Central Asian origin; studies investigating population changes associated with domestication from 273 ancient horse genomes also suggest that the western Eurasian steppe region is the home of the Akhal-Teke horse [6]. In contrast, there is a great deal of evidence indicating that the Akhal-Teke horse originated in the Middle East. The horse was widely documented in Middle Eastern epics and was referred to as the “supernatural horse,” “celestial horse,” or “bleeding stallion” [7]. Various molecular markers have been used in previous studies, which suggested that the Akhal-Teke horse belongs to the Middle Eastern horse breed clade [8–11]. Other archaeological studies at various sites suggest that the Akhal-Teke horse originated in Southwest Asia and migrated from Central Asia to the Middle East, or originated as a cross between the Tarpan and a domesticated native horse from Central Asia [12]. Thus, although several hypotheses have been proposed for the origin of the Akhal-Teke horse, many aspects surrounding its origin are unknown. In this study, we combined genotypic data from genome-wide high-density single-nucleotide polymorphism (SNP) arrays with equine mitochondrial hypervariable regions 1 (HVR-1) sequence data to elucidate the origin of the Akhal-Teke. The results provide new insights into the origin of the Akhal-Teke horse.

MATERIALS AND METHODS

The experiment was approved by the Institutional Animal Care and Use Committee of the Gansu Agricultural University under permit number No.GAU-LC-2018-12.

Datasets

The genome-wide SNP genotype data for Kazakh and Przewalskii horses were obtained from our previous study [13] that used the Illumina horse SNP 70 bead array (Illumina, San Diego, CA, USA) with the Illumina Equine SNP 50 bead array data (animalgenome.org) (Illumina, USA) from 20 Western horse breeds genotyped by Petersen et al [14,15]. The SNP70 and SNP50 arrays shared a total of 43,782 common markers and were used in subsequent analyses (Table 1). Quality control of the raw data was performed with PLINK v1.07 [16] and SNPs with call rates lower than 0.90, p<10⁻⁵ or MAF<0.01 were excluded. Finally, 36,734 sites remained after quality control.

Table 1

Information on the horse breeds included in this study

RESULTS

Genetic relationships between the Akhal-Teke horse and other horse breeds of the world

The results of the nearest neighbor-joining tree showed that Akhal-Teke horses and Arabian horses form a cluster together, indicating that they are closely related. However, no close kinship was found between the Akhal-Teke horse and Caspian or Tuwa horses (Figure 1A), and principal component analysis results further confirmed the close genetic link between the Akhal-Teke horse and Middle Eastern horse populations (Figure 1B). Estimation of the genetic composition of the Akhal-Teke horse and other horses showed a minimum cross-validation statistic when K was set to 15. The results showed that the Akhal-Teke horse had a similar genetic background to Middle Eastern horse breeds, especially Arabian horses (Figure 1D). The Treemix results reflect several well-known gene flows (Figure 2A). For example, Thoroughbreds were observed to have had a significant genetic influence on Quarter horses, and Caspian horses were influenced by gene flow from Arabian horses. However, no gene exchange between Akhal-Teke and other horse breeds was detected and LD decay analysis showed that Akhal-Teke horses had LD decay levels similar to those of Tuva, Caspian, and Arabian horses. The results of this analysis indicated a moderate level of LD decay in the Akhal-Teke horses(Figure 2B).

Figure 1

Genetic relationships between the Akhal-Teke horse and other horse breeds of the world. (A) Relationship between the Akhal-Teke horse and other horse populations shown by a neighbor-joining phylogenetic tree. Study populations are shown in different colors. (B) Principal component analysis of the Akhal-Teke horse and other horse populations. The different study populations are indicated by different colors. (C) Cross-validation under different K-values. (D) ADMIXTURE analysis of the Akhal-Teke horse and world horse populations. (K) Admixture plots of breeds analyzed with different numbers of assumed ancestors. Each column represents an individual. Columns represent individuals, groups of columns represent horse populations.

Figure 2

Population gene flow and the degree of domestication. (A) Migration showing four migration events in the Akhal-Teke horse and worldwide horse populations, shown by the TreeMix program; (B) linkage disequilibrium (LD) decay analysis of 22 horse populations from around the world, including Akhal-Teke horses. The different horse populations are represented by different colors.

IBD analysis of the Akhal-Teke horse

The results of the IBD shared-fragment analysis indicated the presence of relatively large shared fragments between the Akhal-Teke horse and Arabian horses and their derivatives (Thoroughbreds and Quarter horses). On the other hand, the length of shared fragments between the Tuva and Caspian horses and the Akhal-Teke horses was small. Nordic horses, especially the Exmoor breed, shared the fewest fragments with the Akhal-Teke horse (Figure 3).

Figure 3

Estimated lengths of identity by descent (IBD) segments common to the Akhal-Teke and other horse breeds from around the world.

Estimation of possible ancestry

While the above analysis indicated a close genetic relationship between the Akhal-Teke horse and the Arabian horse and several other breeds, its ancestry was still undetermined. Therefore, the ADMIXTOOLS software package was applied to analyze the best ancestral composition. Exmoor and Thoroughbred horses were used as reference populations and the Akhal-Teke horse was used as the target population (X). The higher the output f4 ratio, the greater the likelihood of the reference population being ancestral to the target population. The results showed that Middle Eastern horses had a larger f4 ratio than the Central Asian horse population (Table 2).

Table 2

Analysis of the F4 ratio of Akhal-Teke horses

Evaluation of ancestry ratios of Akhal-Teke horses

The horse groups most likely to be ancestral to the Akhal-Teke horse were selected from the Middle Eastern and Central Asian horse populations. Alder v103 software was used to estimate the ancestral ratios of Arabian horses and Caspian horses in the Akhal-Teke horse. The Caspian and Arabian horses were used as the two reference groups, and the Akhal-Teke horse was used as the target group. The results showed that Arabian horses accounted for 26% of the ancestry of Akhal-Teke horses and Caspian horses accounted for 18% of the ancestry of Akhal-Teke horses (Figure 4).

Figure 4

Evaluation of ancestry ratios of the Akhal-Teke horse. Arabian and Caspian horses were used as the reference populations, and the Akhal-Teke horse as the target population. Ancestry ratios are indicated by a percent sign.

Mitochondrial HVR-1 regions

The 544 hypervariable regions of the 24 breeds were classified into seven major groups according to the geographical location of the population. These groups were the Middle Eastern, Akhal-Teke, Central Asian, European, Tarpan, Przewalski, and Thoroughbred horse groups. A total of 175 haplotypes were detected (Figure 5A). The results of the central join network showed that Akhal-Teke horses shared several haplotypes with horses from the Middle East and Central Asia. However, there were no common haplotypes between Akhal-Teke horses and Tarpan horses. The 175 haplotypes were also divided into five haplotype groups, A–E, with Tarpan horses found only in the E haplotype group and other horses marked for inclusion. The results of genetic differentiation between groups (F_ST) analysis showed that the F_ST values between the Akhal-Teke and Tarpan, Przewalski, and Andalusian horses were relatively high compared to the other horse populations (Figure 5B). Mismatch analysis of the Akhal-Teke horse population showed multiple peaks (Figure 5C), indicating that there was no expansion in the history of this breed. Ultimately, the hypotheses that the Akhal-Teke horse is a cross between a local Central Asian breed and the Tarpan horse or that the Akhal-Teke horse migrated from Central Asia to the Middle East were not supported.

Figure 5

Results of the mitochondrial HVR-1 nucleotide sequence analysis. (A) Median-joining network between groups constructed in the mitochondrial HVR-1 region. Study populations are represented by different colors. All haplotypes were divided into five haplotype subgroups, with each subgroup circled and represented by a different letter. (B) Indicators of genetic differentiation between horse groups (F_ST). Color intensity indicates the degree of differentiation between groups. (C) Mismatch analysis diagram of the Akhal-Teke horse. HVR-1, hypervariable region 1.

DISCUSSION

The Akhal-Teke horse is one of the oldest and most famous horse breeds. It is now mainly distributed in Turkmenistan [2]. Most of the horses have golden chestnut hair, sometimes black or brown. The shoulder height is 150 to 163 cm. Akhal-Teke horses belong to the riding-horse category due to their sporting characteristics [23]. However, despite being one of the few purebred horse breeds in the world, the origin of the Akhal-Teke horse is largely unknown. The most likely hypothesis is that the Akhal-Teke horse originated in Central Asia, adjacent to the Eurasian steppe, where the earliest domestication of horses took place [24]. Conversely, the hypothesis that the Akhal-Teke originated in the Middle East is less plausible, given the paucity of archaeological evidence associated with Akhal-Teke horses found in the Middle East [12]. As the Turkmen horse [25] is the most important ancestor of the modern Akhal-Teke horse [26], another proposal is that the Akhal-Teke horse originated in Central Asia and was introduced to the Middle East. To further test this hypothesis, this study applied genome-wide SNP data from 22 breeds from around the world to analyze the genetic relationship between the Akhal-Teke horse and other breeds. The results suggested that the Akhal-Teke horse is not closely related to the Caspian or Tuva horse and that there was little genetic exchange between the two, which does not support the longstanding hypothesis that the Akhal-Teke horse originated in Central Asia [27]. However, closer genetic links between the Akhal-Teke horse and horses from the Middle East were found. Genetic exchanges between the Akhal-Teke horse and other horse breeds could not be detected by the Treemix analysis, which may be due to the relatively limited geographic distribution of the Akhal-Teke horse. Shared IBD fragment analysis of the Akhal-Teke horse provided further evidence that Middle Eastern horses are more likely to be the ancestors of the Akhal-Teke horse than other horses with the admixture results indicating that both Caspian and Arabian horses have similar genetic backgrounds to the Akhal-Teke horse. Furthermore, the ancestry estimates showed that Arabian horses contributed a greater proportion to the ancestry of the Akhal-Teke horse. The geographical distributions of the Arabian and Akhal-Teke horses almost overlap, which would have facilitated genetic exchange between the two breeds. These results are consistent with fossil finds reported in the literature that the Libyan peninsula, adjacent to the Middle East, played an essential role in the process of horse domestication as a major center of domestication [24]. The fact that the modern Akhal-Teke breed shares haplogroups with horses from Turkmenistan and Kazakhstan is more clearly delineated by the mitochondrial data, which also points to the likelihood that the modern Middle East is the likely origin of the Akhal-Teke horse [2,10].

Due to the lack of whole-genome SNP array data for the Tarpan horse, data on the mitochondrial HVR-1 region of the Tarpan horse [16,28–30] was obtained from NCBI. The genetic relationship between the Akhal-Teke and Tarpan horses was examined. The haplotypes were divided into five groups. All Tarpan horses were found in the E haplotype group, and no haplotypes were shared between Tarpan horses and other horses [29]. Therefore, our mtDNA results do not support the hypothesis that the Akhal-Teke horse is a hybrid between the Tarpan horse and a local Central Asian horse. This result is also consistent with the results from the genome-wide SNP data analysis in this study. Furthermore, the mtDNA data also showed that there was no historical group expansion of the Akhal-Teke horse, which did not support the hypothesis that the Akhal-Teke migrated from Central Asia to the Middle East. Thus, the results from both the genomic and mtDNA data suggest that the Akhal-Teke most likely originated in the Middle East.

CONCLUSION

In this study, data from horse whole-genome SNP array chips and mitochondrial hypervariable regions indicated that the Akhal-Teke horse may have originated in the Middle East rather than in Central Asia. The study also found no support for the hypothesis that the breed was derived from a historical migration from Central Asia to the Middle East, or that it originated from a hybrid between the Tarpan horse and a domesticated local horse breed from Central Asia.

DATA AVAILABILITY STATEMENT

The access numbers of 544 mitochondrial HVR-1 sequences of 24 breads from NCBI are summarized in the excel table. Mitochondrial HVR-1 Nucleotide Sequence.xlsx

Article information Continued