- Research article
Maternal phylogenetic relationships and genetic variation among Arabian horse populations using whole mitochondrial DNA D-loop sequencing
BMC Geneticsvolume 14, Article number: 83 (2013)
Maternal inheritance is an essential point in Arabian horse population genetics and strains classification. The mitochondrial DNA (mtDNA) sequencing is a highly informative tool to investigate maternal lineages. We sequenced the whole mtDNA D-loop of 251 Arabian horses to study the genetic diversity and phylogenetic relationships of Arabian populations and to examine the traditional strain classification system that depends on maternal family lines using native Arabian horses from the Middle East.
The variability in the upstream region of the D-loop revealed additional differences among the haplotypes that had identical sequences in the hypervariable region 1 (HVR1). While the American-Arabians showed relatively low diversity, the Syrian population was the most variable and contained a very rare and old haplogroup. The Middle Eastern horses had major genetic contributions to the Western horses and there was no clear pattern of differentiation among all tested populations. Our results also showed that several individuals from different strains shared a single haplotype, and individuals from a single strain were represented in clearly separated haplogroups.
The whole mtDNA D-loop sequence was more powerful for analysis of the maternal genetic diversity in the Arabian horses than using just the HVR1. Native populations from the Middle East, such as Syrians, could be suggested as a hot spot of genetic diversity and may help in understanding the evolution history of the Arabian horse breed. Most importantly, there was no evidence that the Arabian horse breed has clear subdivisions depending on the traditional maternal based strain classification system.
The fact that the Arabian horse breed is one of the oldest pure breeds and among the most widespread on earth is not disputed . It is the most influential breed throughout the world  and has been used to genetically improve several other breeds like the Thoroughbred  and the Lipizzan . Furthermore, Western Arabian populations were created in Europe and the USA using stock originally imported from Middle Eastern Arabian populations from sources such as Syria and the Arabian peninsula no longer than 200 years ago [2, 5].
The advent of mitochondrial DNA (mtDNA) sequencing in population genetics during the 1970s produced a revolutionary change regarding historical, biogeographic and phylogenetic perspectives on intra- and inter-species genetic structure . Single base pair substitutions have a frequency of up to 10 times higher than nuclear DNA . In addition, strict maternal inheritance of mtDNA  and lack of recombination make sequence polymorphism identification a unique application for domestic animals genetic studies not provided by nuclear genes . Therefore, mtDNA has been widely used as a highly informative tool to infer intra- and inter-species phylogenetic relationships , and it has successfully been applied to characterize intra-breed variation and the origin of many horse breeds [2, 5, 10–17]. Furthermore, it can be used to track breed migration and distribution by comparing maternal lines among different populations [18, 19].
Unlike the protein-coding gene region, the displacement loop hyper-variable region of mtDNA (D-loop), represented by about 1200 bp, is of particular interest because of its high level of sequence variation . D-loop also is known to be an under-represented mtDNA in Nuclear Mitochondrial sequences (NUMTs) in mammalian species so that it seldom produces NUMTs [21, 22]. The D-loop region in horses contains two highly variable segments (HVR1 and HVR2), four conserved blocks (CSB), and variable repeats of 8 bp motif [23, 24].
The traditional pattern of breeding Arabian horses affords special opportunities to evaluate variation in matrilineal markers, such as mitochondrial DNA. From a glance of historical records, the Arabian horse breed, in the desert, consists of five strains (RASANs) based upon dam line: Kahlila, Saklawia, Abiah, Shweemat, and Muanakii (some breeders and historians refer to an additional three RASANs which are Hamadania, Dahmaa and Hadbaa). The traditional breeders in the Middle East desert (Bedouins) have preserved the purity of the Arabian by avoiding any cross-breeding between the Arabians and non-Arabians and maintaining strictly separated RASANs. Consequently, all individuals within a RASAN are expected to share the same maternal family line, and they should have similar mtDNA haplotype.
While many studies have been done in horses using mtDNA, only a few included Arabians [2, 5, 11, 26]. Also, the Arabians used were mainly collected from Western populations. Most of the previous studies related to Arabian population genetics used only about 400 bp out of 1200 bp of the mtDNA D-loop. In the present study, we sequenced the whole mtDNA D-loop of Arabian horses collected from the Middle East as well as from Western populations. Our study was designed to investigate the maternal diversity and phylogenetic relationships of Arabian populations and to examine the traditional classification system of the Arabian breed (RASANs system) that depends upon maternal family lines.
Table 1 shows the number of haplotypes (NHap), haplotype diversity (HapD), average number of nucleotide differences (k), the number of polymorphic sites (NPS) and nucleotide diversity (π) for each population. A total of 74 haplotypes from 60 polymorphic sites were found in 271 horses from 11 populations by using the HVR1. NHap increased to 97 using the whole D-loop sequences. Although π decreased from 0.022 to 0.015, NPS increased from 60 to 99 and k increased from 9.7 to 14.5 comparing the HVR1 to the whole D-loop, respectively. The highest HapD values among all tested Arabian populations were in the Syrian, Shagya Arabian and Iranian Arabian populations 0.97, 0.97, 0.96, respectively. The non-Arabian populations also showed high values of HapD, 1.0 in Mongolian and 0.93 in Caspian. All American-Arabian populations -USA-Saudi, USA-Egyptian, USA-Egyptian & Saudi mix and Davenport- showed relatively low HapD ranging between 0.74 and 0.83.
The tested samples were then grouped into strains according to pedigree records and regardless of their populations. We could assign 191 out of 271 samples into seven strains (RASANs). Table 2 shows NHap, HapD, k and π for each strain. A total of 44 haplotypes from 52 polymorphic sites were found in these 191 horses of the seven strains using the HVR1 part of the D-loop. The NHap increased to 55 using the whole D-loop sequences. Only Shweemat strain had all individuals with a single haplotype. Hadbaa and Dahmaa also had low NHap, 3 and 2, respectively. Kahlila was the most variable strain showing 26 haplotypes. The total NHap calculated from all individuals together (NHap = 55) was less than the sum of NHap calculated from each strain separately due to some shared haplotypes among strains.
Figure 1 shows the consensus Neighbor-joining tree of the 97 haplotypes found from all tested populations. No single population was found only in one cluster and different populations shared haplotypes. Fifteen haplotypes: 24, 15, 29, 30, 44, 26, 11, 22, 33, 14, 23, 55, 16, 4 and 28 were found in at least two Arabian populations; for example, haplotype 24 was found in two populations: USA-Saudi and USA-Egyptian; haplotype 14 appeared in five populations: USA-Saudi, Syrian, USA-Egyptian, USA-Egyptian & Saudi mix and Polish Arabian. In addition, haplotype 4 was found in Arabian and non-Arabian populations: Syrian and Iranian Arabian with Caspian. The dendrogram gave seven main clades plus the out-group. Syrian population was the most variable among all populations with individuals found in all clades.
Figure 2 shows the consensus Neighbor-joining tree of the 55 haplotypes found in the individuals who were assigned to their strains. None of the tested strains, except Shweemat, was represented by a single haplotype or phylogenetically close haplotypes. Each one of the thirteen haplotypes: 16, 23, 22, 12, 27, 14, 29, 15, 18, 75, 11, 74 and 17 was found in at least two strains. For example, haplotype 16 was present in two strains -Hamadania and Dahmaa- and haplotype 23 in three strains -Abiah, Kahlila and Hamadania-. The most frequent mixing was noticed between Saklawia and Kahlila strains. The Kahlila strain was the most variable among all strains and its individuals were distributed among all clades. The bootstrap values were not high in the trees depicted in Figure 1 and Figure 2. However, a similar pattern of phylogeny resulted by running Bayesian approaches (data not shown).
After taking into account mutational hot spots for the median-joining network (MJ network), the number of haplotypes dropped from 97 to 86. Figures 3 and 4 show the MJ network based on 951 bp of the D-loop representing 271 samples by 86 haplotypes. While in Figure 3 each haplotype is shown by the proportion of the different populations included in this haplotype, in Figure 4 each haplotype is shown by the proportion of different strains. The MJ network showed 14 haplogroups -A, B, C, D, E, G, I, J, L, M, N, P, Q and R- as was defined by Achilli et al. .
As shown in Figure 3, each of the 13 haplogroups -A, B, C, D, E, G, I, L, M, N, P, Q and R- contained identical or very close haplotypes from at least two populations. The highest number of populations was found in the haplogroup L. The Arabian populations were represented in all haplogroups except J. The non-Arabian samples were placed in the haplogroups -A, E, I, L, M, N, Q and R- and -A, B, C, J, and P- for the Caspian and the Mongolian populations, respectively. The Syrian population was the most variable with individuals distributed across all haplogroups except J and R. Individuals from the Syrian population had identical or very close haplotypes to individuals from all other Arabian and non-Arabian populations. The Davenport was the least variable Arabian population with only three haplogroups -I, L and P-.
Figure 4 shows that individuals from different strains shared a single haplotype. Identical matching between two or more individuals from different strains was seen in 13 cases. Also, matching was found between known strains and other Arabian groups -Polish Arabian, Shagya Arabian and Iranian Arabian- and non-Arabian populations -Caspian and Mongolian-. In addition, individuals from a single strain were found in distinctive haplogroups (for example: strain Hamadania in haplogroups P, C and R). Kahlila strain was the most variable with individuals distributed across all haplogroups except J and R. All of the unknown-strain samples were identical or very close to samples of known strains. Although Shweemat strain was the least variable with only haplogroup L, it was very close to samples from some other strains.
Figure 5 shows the PCoA plot of the two first axes which explain 26.52% and 18.77% of the variability, respectively, and grouped the 98 haplotypes into five clusters. Cluster I included a combination of three haplogroups -M, N and R-. Cluster II consisted of two haplogroups (P and Q). Cluster III included seven haplogroups -A, B, C, E, G, I, and J-. Cluster IV had only haplogroup D, and Cluster V included only haplogroup L. The clustering by PCoA did not show any differentiation among haplotypes that came from different populations (or different strains) but it showed that each cluster contained a mixture of individuals that represented different populations (or strains).
AMOVA showed that the proportion of the variation among populations was 8.25% and the frequency of the variation within populations was 91.75%. The fixation index was equal to 0.083.
The pairwise FST values are shown in Table 3. Out of 55 pairwise FST values 28 comparisons had FST values between 0 and 0.05 showing little genetic differentiation while 21 comparisons had FST values between 0.05 and 0.15 showing moderate genetic differentiation, six comparisons had FST values between 0.15 and 0.25 showing great genetic differentiation. Negative FST values were recorded in some comparisons and these equate to zero FST values. While most of the lowest FST values were seen between Syrian and eight other populations -Caspian, USA-Egyptian, USA-Egyptian & Saudi mix, Saudi, Shagya Arabian, Polish Arabian, Mongolian and Iranian Arabian-, the highest FST values were between the Davenport and five other populations -USA-Egyptian, Saudi, USA-Saudi, Polish Arabian and Mongolian-. None of the comparisons showed values corresponding to very great genetic differentiation.
This study presents the first description of maternal genetic diversity based upon the whole mtDNA D-loop of native Arabian horses sampled from Syria, Iran and Saudi Arabia, as well as of Western Arabian populations. One of the unique aspects of this study is the inclusion of the traditional classification system (RASANs or strains system) of native Arabians.
HVR1 and the whole mtDNA comparison
Most previous maternal diversity studies of horses are based upon sequencing of the HVR1 [5, 12, 15–17, 26–30]. Our results of the comparison between HVR1 and the entire mtDNA D-loop showed that the variability in the upstream region of the D-loop revealed differences among 22 additional haplotypes which had identical sequences in the HVR1. This agrees with Kavar et al.  where such a pattern of variability had been found in the D-loop of the Lipizzan horse breed. Higher haplotype diversity (HapD = 0.98) and average number of nucleotide differences (k = 14.5) were found by using the whole mtDNA D-loop compared with the HVR1 (HapD = 0.98 and k = 9.5) (Table 1). Thus, using the whole mtDNA D-loop is more robust and powerful than using the HVR1 alone for analysis of genetic diversity of the mtDNA in horses. Similar results have been reported in goats .
Population genetic diversity
Maternal genetic diversity of the Arabian populations described in this study was similar to that reported in some other breeds [17, 33]. Although Syrian, Shagya Arabian and Iranian Arabian populations had equally high HapD values (Table 1), the Syrian population was the most variable based on the consensus Neighbor-joining tree (Figure 1) where the Syrian individuals were found in eight clades compared to the Iranian Arabian and Shagya Arabian individuals found only in five and three clades, respectively. This result also was supported by the MJ-network (Figure 3) where the Syrian population was represented in 12 haplogroups compared to the Iranian Arabian and Shagya Arabian with six and five haplogroups, respectively. According to Achilli et al.  there is a total of 18 major haplogroups of horses throughout Asia, Middle East, Europe and America; our results showed that the Syrian population covers 12 of those 18 haplogroups showing a wide maternal genetic diversity. In our opinion, the huge diversity of the Syrian population is not a consequence of recent animal breeding or outcrossing but instead a feature that was already present in this very old population, and this is supported by results of Cieslak et al. . In addition, the huge diversity in Arabian populations is consistent with the multiple origins in the maternal lineages of domestic horse breeds reported by other studies [15, 16, 28, 34]. Some of the Syrian individuals were represented in haplogroup D (Figure 3), haplogroup E according to Jansen et al. , that was reported as a very rare and old haplogroup which may date back as far as Bronze age [15, 29, 35].
The American-Arabian populations showed relatively low HapD values and were represented in a limited number of haplogroups. Davenport population was the least variable with only three haplogroups -I, L and P-. The low maternal diversity found in the American-Arabian populations is probably due to the founder effect. This result is supported by our previous work done by using microsatellite markers where American-Arabian populations showed less genetic variability compared with Middle Eastern populations . Also, Polish Arabians did not show a very high genetic diversity with only 6 haplotypes distributed in four haplogroups. This result did not match with Glazewska et al.  where 14 distinct haplotypes were reported. This could be due to sample size or because the horses we used came from close maternal lines.
Population relationships and genetic structure
The low bootstrap values of the Neighbor-joining trees in Figure 1 and Figure 2 are primarily due to the overall high degree of relationship among horses . Low bootstrap values have been reported in many mtDNA studies in horses [10, 14, 28, 30, 38]. Although bootstrap values were low, the populations consistently fell into the same groupings in the trees. The consensus Neighbor-joining tree (Figure 1) and the MJ-network (Figure 3) show that individuals from different populations share identical haplotypes. This indicates possible gene flow among those populations or common ancestry. Identical maternal lines were found between the Syrian and Polish Arabian populations revealing that Syrian mares were probably part of Polish Arabian founders, or some horses were recently introduced to this population. The identical maternal lines that were found between the American Arabian populations -USA-Saudi, USA-Egyptian & Saudi mix and USA-Egyptian- and populations from the Middle East -Syrian, Saudi and Iranian Arabian- confirms that the current registered Arabian horses in America have been primarily founded by mares exported from the Middle East . While Shagya Arabian population is thought to be descended from a Syrian stallion , our results show some shared maternal lines between Shagya Arabians and Syrians suggesting a maternal contribution of Syrian horses in Shagya Arabian population, or possibly recent gene flow between these two populations. Furthermore, the phylogenetic analysis revealed that different populations, including Arabian and non-Arabian, often had very close haplotypes, and none of these populations formed a distinct clade. These results together reveal the mixed origin and/or a likely common ancestor of these populations. The genetic clustering analysis using both phylogenic (Figures 1 and 3) and PCoA (Figure 5) did not show any clear pattern of differentiation among all populations. Haplotypes within a population were found in separate haplogroups. Similar results have been reported in other studies of horse mtDNA [10, 12, 16]. FST analysis supports this unclear pattern of differentiation showing high rates of mtDNA sharing between populations, with in some cases negative FST values. FST values sometimes are produced by software which uses algorithms that include sampling error corrections, such as Arlequin, when the true FST values are close to zero , and usually appear when there are great differences between two random individuals from the same population rather than between two random individuals from different populations . These negative values represent program idiosyncrasies and are effectively zero . AMOVA results also support high within group variation with 91.75% of variability as within population variation.
Strain relationships and classification system
In the Middle East, strain breeding is still an important factor in the Arabian horse breed . According to Bedouin breeding traditions, Arabian horses were subdivided into strains depending on the maternal lineage. The phylogenetic and principle coordinate analyses in our study using 191 samples, of known strains, showed no evidence that the Arabian breed has clear divisions based upon traditional strain classification. There are four points that support this finding. First, 13 cases revealed that individuals from different strains shared a single haplotype. For example, haplotype 23 was found in individuals that came from three different strains -Abiah, Kahlila and Hamadania-; haplotype 29 was in individuals from three strains -Abiah, Kahlila and Saklawia- (Figure 2). Second, individuals from different strains were found in a single haplogroup. For example, haplogroup P was seen in five strains -Kahlila, Saklawia, Abiah, Hadbaa, and Hamadania- (Figure 4). Third, each of the strains -Kahlila, Saklawia, Abiah, Dahmaa, Hadbaa and Hamadania- was represented in clearly separated haplogroups. For example, Kahlila was found in 12 haplogroups (Figure 4). Finally, PCoA did not show any pattern of clustering that fits strains subdivision (Figure 5). Our results agree with the conclusion reported by Bowling et al.  about American Arabian horses.
It is possible to have some minor mistakes in the pedigree records of any breed , but with our results we can confirm that these mistakes, if they existed in the records that we used, cannot be the reason behind having the huge admixture among tested strains. We do not suspect admixture into the Arabian horse breed, but it is clear that the pedigree records of the Arabian breed were not built using robust genetic tools that can recognize distinct maternal lines in the establishment of the pedigree.
Another important factor in the Bedouin breeding traditions is the sub-strain subdivisions (MARBATT) that subdivides each Arabian strain into related groups depending on the tribe’s or owner’s name . Although we did not test the sub-strain subdivisions of Arabians in our study because of a lack of information, we can say that the sub-strain system might be able to partially explain the third point mentioned above (related to the differences among individuals from same strain), but it does not answer the other questions.
The maternal phylogenetic analysis of native Arabian horses in our study revealed (1) that the whole mtDNA D-loop sequence was more powerful to analyze the genetic diversity in Arabian horses than using just the HVR1, (2) that the maternal genetic diversity was wide in the Arabian horse populations especially in the Syrian population, (3) that there was no clear pattern of differentiation among all tested populations, (4) that the Syrian mares probably had maternal contributions to the Polish Arabian and Shagya Arabian populations, (5) that the current registered Arabian horses in America have been primarily founded by mares exported from the Middle East, and most importantly (6) that there was no evidence that the Arabian breed has clear subdivisions depending on the traditional strain classification system.
Sampling and DNA extraction
Hair samples were collected from 271 horses representing Middle Eastern Arabian, Western Arabian and non-Arabian populations. All tested horses were unrelated from the mothers’ side at least for 3 generations based upon their pedigree records. Tables 1 and 2 show the number of animals used in relation to their populations and strains, respectively. Additional file 1: Table S1 included a summary about populations sampling background and sampling permissions. Total genomic DNA was extracted from hair follicles using the PUREGENE® DNA purification kit following the manufacturer’s instructions.
Whole D-loop sequencing and data analysis
We designed two pairs of primers based upon the horse mtDNA sequence reference X79547 . We also considered the outcomes reported by Nergadze et al.  to minimize the possible amplification of NUMTs that may overlap with D-loop. The designed primers were used to amplify the upstream part between sites 15440 and 16108 (Forward: 5′-AGCTCCACCATCAACACCCAAA-3′. Reverse 5′-CCATG GACTGAATAACACCTTATGGTTG-3′) and the downstream part between sites 16377 and 16642 (Forward 5′-ACCTACCCGCGCAGTAAGCAA-3′. Reverse 5′-AC GGGGGAAGAAGGGTTGACA-3′). Polymerase chain reactions were done for each part separately using the protocol described by Cothran et al. . A total of four sequencing reactions for each sample, including both strands in each part, were carried out using the BigDye® Terminator v1.1 Cycle Sequencing Kit (Applied Biosystems, USA). Sequencing products were purified with BigDye® XTerminator™ Purification Kit (Applied Biosystems, USA). DNA sequences were determined using the ABI 3130 xl Genetic Analyzer (Applied Biosystems, USA). Editing and aligning all sequences were carried out by MEGA 4  using the horse mtDNA sequence X79547 as a reference. Haplotype sequences included in this study were entered into the National Center for Biotechnology Information (NCBI) GenBank database available at http://www.ncbi.nlm.nih.gov/ with the accession numbers [NCBI: KC840701-KC840797]. The statistical quantities for the DNA sequences, including number of haplotypes and haplotype diversity and nucleotide diversity, were carried out using DnaSP 5.10.1 . The statistical analysis was done for each population, as well as for each strain, using two sources of data HVR1, (450 sites) and whole D-loop sequences (951 sites).
Phylogenetic analysis of the haplotypes using a whole D-loop sequence was carried out with the PHYLIP software package  based on the Kimura 2-parameter model to calculate genetic distances on the assumption of an equal substitution rate per site . A consensus tree was also constructed with PHYLIP using the Neighbor-joining method  with 1000 bootstrap repetitions. The donkey (Equus asinus) mtDNA sequence [NCBI: nc_001788]  was used as an out-group [10, 11].
Another approach for phylog enetic analysis was carried out by drawing the median-joining network (MJ network)  in accordance with the haplotype sequences of the whole D-loop using the NETWORK 4.6.1 software (available at http://fluxus-engineering.com). Default settings were applied (r = 2, ϵ = 0) , and preliminary trials were done in order to determine the mutational hotspots. Four mutational hot spots were excluded and an additional four were down-weighted into 0.5 [15, 16]. Each haplotype in the MJ network was shown by color codes representing the proportions of different strains (or populations) depending on the individual frequencies in each haplotype. Furthermore, the haplotype sequences were compared to the NCBI database using the BLAST search as implemented in MEGA 4, and haplogroups were named as defined by Achilli et al. .
To represent the genetic structure and differentiation of tested populations, principal coordinate analysis (PCoA), analysis of molecular variance (AMOVA) and pair-wise FST were applied. PCoA of the dissimilarity matrix according to Kimura  based on 951 bp of the 98 haplotypes sequences was carried out using DARwin 5.0 [51, 52]. AMOVA and pair-wise FST were done using the Kimura 2-parameter model with 1000 permutations and were carried out with Arlequin 3 . For the interpretation of pair-wise FST results, we followed the suggestion that refers that a value between 0–0.05 indicates little genetic differentiation; a value between 0.05 and 0.15, moderate differentiation; a value between 0.15 and 0.25, great differentiation; and values above 0.25, very great genetic differentiation [54, 55].
Hendricks BL: International Encyclopedia of Horse Breeds. 1995, Norman, OK, USA: University of Oklahoma Press
Głażewska I: Speculations on the origin of the Arabian horse breed. Livest Sci. 2010, 129 (1–3): 49-55.
Bowling A, Ruvinsky A: Genetic aspects of domestication. The Genetic of the Horse. 2000, UK: CAB International
Zechner P, Solkner J, Bodo I, Druml T, Baumung R, Achmann R, Marti E, Habe F, Brem G: Analysis of diversity and population structure in the Lipizzan horse breed based on pedigree information. Livest Prod Sci. 2002, 77 (2–3): 137-146.
Bowling AT, Del Valle A, Bowling M: A pedigree-based study of mitochondrial D-loop DNA sequence variation among Arabian horses. Anim Genet. 2000, 31 (1): 1-7. 10.1046/j.1365-2052.2000.00558.x.
Avise JC: Molecular Markers, Natural History and Evolution. 1994, New York: Chapman & Hall
Brown WM, George M, Wilson AC: Rapid Evolution of Animal Mitochondrial-DNA. Proc Natl Acad Sci USA. 1979, 76 (4): 1967-1971. 10.1073/pnas.76.4.1967.
Hutchison CA, Newbold JE, Potter SS, Edgell MH: Maternal inheritance of mammalian mitochondrial DNA. Nature. 1974, 251 (5475): 536-538. 10.1038/251536a0.
Mirol PM, Garcia PP, Vega-Pla JL, Dulout FN: Phylogenetic relationships of Argentinean Creole horses and other South American and Spanish breeds inferred from mitochondrial DNA sequences. Anim Genet. 2002, 33 (5): 356-363. 10.1046/j.1365-2052.2002.00884.x.
Vila C, Leonard JA, Gotherstrom A, Marklund S, Sandberg K, Liden K, Wayne RK, Ellegren H: Widespread origins of domestic horse lineages. Science. 2001, 291 (5503): 474-477. 10.1126/science.291.5503.474.
Achilli A, Olivieri A, Soares P, Lancioni H, Kashani BH, Perego UA, Nergadze SG, Carossa V, Santagostino M, Capomaccio S, et al: Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication. Proc Natl Acad Sci USA. 2012, 109 (7): 2449-2454. 10.1073/pnas.1111637109.
Cothran EG, Juras R, Macijauskiene V: Mitochondrial DNA D-loop sequence variation among 5 maternal lines of the Zemaitukai horse breed. Genet Mol Biol. 2005, 28 (4): 677-681.
Kavar T, Dovc P: Domestication of the horse: Genetic relationships between domestic and wild horses. Livest Sci. 2008, 116 (1–3): 1-14.
Lippold S, Matzke NJ, Reissmann M, Hofreiter M: Whole mitochondrial genome sequencing of domestic horses reveals incorporation of extensive wild horse diversity during domestication. Bmc Evol Biol. 2011, 11: 328-10.1186/1471-2148-11-328.
Cieslak M, Pruvost M, Benecke N, Hofreiter M, Morales A, Reissmann M, Ludwig A: Origin and History of Mitochondrial DNA Lineages in Domestic Horses. Plos One. 2010, 5: 12-
Jansen T, Forster P, Levine MA, Oelke H, Hurles M, Renfrew C, Weber J, Olek K: Mitochondrial DNA and the origins of the domestic horse. Proc Natl Acad Sci USA. 2002, 99 (16): 10905-10910. 10.1073/pnas.152330099.
Guastella AM, Zuccaro A, Criscione A, Marletta D, Bordonaro S: Genetic Analysis of Sicilian Autochthonous Horse Breeds Using Nuclear and Mitochondrial DNA Markers. J Heredity. 2011, 102 (6): 753-758. 10.1093/jhered/esr091.
Kivisild T, Reidla M, Metspalu E, Rosa A, Brehm A, Pennarun E, Parik J, Geberhiwot T, Usanga E, Villems R: Ethiopian mitochondrial DNA heritage: tracking gene flow across and around the gate of tears. Am J Hum Genet. 2004, 75 (5): 752-770. 10.1086/425161.
Matisoo-Smith E, Robins J: Mitochondrial DNA evidence for the spread of Pacific rats through Oceania. Biol Invasions. 2009, 11 (7): 1521-1527. 10.1007/s10530-008-9404-1.
Aquadro CF, Greenberg BD: Human mitochondria DNA variation and evolution: analysis of nucleotide sequences from seven individuals. Genetics. 1983, 103: 287-312.
Tsuji J, Frith MC, Tomii K, Horton P: Mammalian NUMT insertion is non-random. Nucleic Acids Res. 2012, 40 (18): 9073-9088. 10.1093/nar/gks424.
Mourier T, Hansen AJ, Willerslev E, Arctander P: The Human Genome Project reveals a continuous transfer of large mitochondrial fragments to the nucleus. Mol Biol Evol. 2001, 18 (9): 1833-1837. 10.1093/oxfordjournals.molbev.a003971.
Xu XF, Arnason U: The Complete Mitochondrial-DNA Sequence of the Horse, Equus-Caballus - Extensive Heteroplasmy of the Control Region. Gene. 1994, 148 (2): 357-362. 10.1016/0378-1119(94)90713-7.
Ishida N, Hasegawa T, Takeda K, Sakagami M, Onishi A, Inumaru S, Komatsu M, Mukoyama H: Polymorphic Sequence in the D-Loop Region of Equine Mitochondrial-DNA. Anim Genet. 1994, 25 (4): 215-221.
Pruski W: Dwa wieki polskiej hodowli koni arabskich (1778–1978) i jej sukcesy na świecie. 1983, Warszawa: PWRiL
Glazewska I, Wysocka A, Gralak B, Sell J: A new view on dam lines in Polish Arabian horses based on mtDNA analysis. Genet Sel Evol. 2007, 39 (5): 609-619. 10.1186/1297-9686-39-5-609.
Ivankovic A, Ramljak J, Konjacic M, Kelava N, Dovc P, Mijic P: Mitochondrial D-loop sequence variation among autochthonous horse breeds in Croatia. Czech J Anim Sci. 2009, 54 (3): 101-111.
Georgescu SE, Manea MA, Dudu A, Costache M: Phylogenetic relationships of the Hucul horse from Romania inferred from mitochondrial D-loop variation. Genet Mol Res. 2011, 10 (4): 4104-4113. 10.4238/2011.October.31.7.
Prystupa JM, Hind P, Cothran EG, Plante Y: Maternal Lineages in Native Canadian Equine Populations and Their Relationship to the Nordic and Mountain and Moorland Pony Breeds. J Heredity. 2012, 103 (3): 380-390. 10.1093/jhered/ess003.
Cozzi MC, Strillacci MG, Valiati P, Bighignoli B, Cancedda M, Zanotti M: Mitochondrial D-loop sequence variation among Italian horse breeds. Genet Sel Evol. 2004, 36 (6): 663-672. 10.1186/1297-9686-36-6-663.
Kavar T, Habe F, Brem G, Dovc P: Mitochondrial D-loop sequence variation among the 16 maternal lines of the Lipizzan horse breed. Anim Genet. 1999, 30 (6): 423-430. 10.1046/j.1365-2052.1999.00557.x.
Kang J, Li X, Zhou R, Li L, Zheng G, Zhao H: Genetic diversity and differentiation of four goat lineages based on analysis of complete mtDNA d-loop. Front Agric China. 2011, 5 (1): 87-93. 10.1007/s11703-010-1061-3.
Cai DW, Tang ZW, Han L, Speller CF, Yang DYY, Ma XL, Cao JE, Zhu H, Zhou H: Ancient DNA provides new insights into the origin of the Chinese domestic horse. J Archaeol Sci. 2009, 36 (3): 835-842. 10.1016/j.jas.2008.11.006.
Aberle KS, Hamann H, Drogemuller C, Distl O: Phylogenetic relationships of German heavy draught horse breeds inferred from mitochondrial DNA D-loop variation. J Anim Breed Genet. 2007, 124 (2): 94-100. 10.1111/j.1439-0388.2007.00636.x.
Kakoi H, Tozaki T, Gawahara H: Molecular analysis using mitochondrial DNA and microsatellites to infer the formation process of Japanese native horse populations. Biochem Genet. 2007, 45 (3–4): 375-395.
Khanshour AM, Juras R, Conant E, Cothran EG: Microsatellite analysis of genetic diversity and population structure of Arabian horse populations. J Hered. 2013, doi:10.1093/jhered/est003
Cothran EG, Luis C: Genetic distance as a tool in the conservation of rare horse breeds. Eaap Public. 2005, 116: 55-71.
Kim KI, Yang YH, Lee SS, Park C, Ma R, Bouzat JL, Lewin HA: Phylogenetic relationships of Cheju horses to other horse breeds as determined by mtDNA D-loop sequence polymorphism. Anim Genet. 1999, 30 (2): 102-108. 10.1046/j.1365-2052.1999.00419.x.
Musick JA: Genetic: Stock identification. Management Techniques for Elasmobranch Fisheries. Volume 474. Edited by: Musick JA, Bonfil R. 2005, Rome: FAO, 62-63.
Arnason E, Palsson S: Mitochondrial cytochrome b DNA sequence variation of Atlantic cod Gadus morhua, from Norway. Mol Ecol. 1996, 5 (6): 715-724. 10.1111/j.1365-294X.1996.tb00368.x.
Humphries EM, Winker K: Discord reigns among nuclear, mitochondrial and phenotypic estimates of divergence in nine lineages of trans-Beringian birds. Mol Ecol. 2011, 20 (3): 573-583. 10.1111/j.1365-294X.2010.04965.x.
Hill EW, Bradley DG, Al-Barody M, Ertugrul O, Splan RK, Zakharov I, Cunningham EP: History and integrity of thoroughbred dam lines revealed in equine mtDNA variation. Anim Genet. 2002, 33 (4): 287-294. 10.1046/j.1365-2052.2002.00870.x.
Nergadze SG, Lupotto M, Pellanda P, Santagostino M, Vitelli V, Giulotto E: Mitochondrial DNA insertions in the nuclear horse genome. Anim Genet. 2010, 41: 176-185.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.
HYLIP (phylogeny inference package). http://evolution.genetics.washington.edu/phylip.html,
Kimura M: A simple methodfor estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Saitou N, Nei M: The Neighbor-Joining Method a New Method for Reconstructing Phylogenetic Trees. Mol Biol Evol. 1987, 4 (4): 406-425.
Xu XF, Janke A, Arnason U: The complete mitochondrial DNA sequence of the greater Indian rhinoceros, Rhinoceros unicornis, and the phylogenetic relationship among Carnivora, Perissodactyla, and Artiodactyla (plus Cetacea). Mol Biol Evol. 1996, 13 (9): 1167-1173. 10.1093/oxfordjournals.molbev.a025681.
Bandelt HJ, Forster P, Sykes BC, Richards MB: Mitochondrial Portraits of Human-Populations Using Median Networks. Genetics. 1995, 141 (2): 743-753.
Perrier X, Flori A, Bonnot F: Data analysis methods. Genetic diversity of cultivated tropical plants. Edited by: Hamon P, Seguin M, Perrier X, Glaszmann JC. 2003, Montpellier: Enfield, Science Publishers
DARwin software. http://darwin.cirad.fr/darwin,
Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform. 2005, 1: 47-50.
Wright S: Evolution and the genetics of population, variability within and among natural populations. 1978, Chicago: The University of Chicago
Hartl D, Clark A: Principles of Population Genetics. 1997, Sunderland, MA: Sinauer Associates Inc.
We would like to thank all the Syrian horse breeders who provided us their horse samples, and Esa S. Allil, Mustafa Haroon and Sultan Alahmad for helping in the Syrian samples collection, Edie Booth for providing Arabian samples with their pedigree information, Dr. G. Cholewinski who provided the Polish Arabian samples and F. Maloufi who provided the Persian Arab samples. A.M. Khanshour was supported by a PhD scholarship from the Islamic Development Bank (IDB MERIT SCHOLARSHIP PROGRAMME). Many thanks to our lab team Dr. Rytis Juras, Dr. Eleanore Conant and Smriti Shankar for their positive support.
We declare that we have no competing interests related to this study.
AMK and EGC designed the experiment. AMK collected the Syrian Samples, performed the lab works and data analysis and wrote the manuscript. EGC provided the other tested samples, revised the manuscript and supervised the project. AMK and EGC read and approved the final manuscript.