Origin and history of mitochondrial DNA lineages in domestic horses

Domestic horses represent a genetic paradox: although they have the greatest number of maternal lineages (mtDNA) of all domestic species, their paternal lineages are extremely homogeneous on the Y-chromosome. In order to address their huge mtDNA variation and the origin and history of maternal lineages in domestic horses, we analyzed 1961 partial d-loop sequences from 207 ancient remains and 1754 modern horses. The sample set ranged from Alaska and North East Siberia to the Iberian Peninsula and from the Late Pleistocene to modern times. We found a panmictic Late Pleistocene horse population ranging from Alaska to the Pyrenees. Later, during the Early Holocene and the Copper Age, more or less separated sub-populations are indicated for the Eurasian steppe region and Iberia. Our data suggest multiple domestications and introgressions of females especially during the Iron Age. Although all Eurasian regions contributed to the genetic pedigree of modern breeds, most haplotypes had their roots in Eastern Europe and Siberia. We found 87 ancient haplotypes (Pleistocene to Mediaeval Times); 56 of these haplotypes were also observed in domestic horses, although thus far only 39 haplotypes have been confirmed to survive in modern breeds. Thus, at least seventeen haplotypes of early domestic horses have become extinct during the last 5,500 years. It is concluded that the large diversity of mtDNA lineages is not a product of animal breeding but, in fact, represents ancestral variability.     

Iberian origins of New World horse breeds

Fossil records, archaeological proofs, and historical documents report that horses persisted continuously in the Iberian Peninsula since the Pleistocene and were taken to the American continent (New World) in the 15th century. To investigate the variation within the mitochondrial DNA (mtDNA) control region of Iberian and New World horse breeds, to analyze their relationships, and to test the historical origin of New World horses, a total of 153 samples, representing 30 Iberian and New World breeds, were analyzed by sequencing mtDNA control region fragments. Fifty-four haplotypes were found and assigned to seven haplogroups. Reduced levels of variation found for the Menorquina, Sorraia, and Sulphur Mustang breeds are consistent with experienced bottlenecks or limited number of founders. For all diversity indices, Iberian breeds showed higher diversity values than South American and North American breeds. Although, the results show that the Iberian and New World breeds stem from multiple origins, we present a set of genetic data revealing a high frequency of Iberian haplotypes in New World breeds, which is consistent with historical documentation.     

Phylogenetic relationships of Argentinean Creole horses and other South American and Spanish breeds inferred from mitochondrial DNA sequences

South American horses constitute a direct remnant of the Iberian horses brought to the New World by the Spanish conquerors. The source of the original horses was Spain, and it is generally assumed that the animals belonged to the Andalusian, Spanish Celtic, Barb or Arabian breeds. In order to establish the relationship between Argentinean and Spanish horses, a portion of the mitochondrial D-loop of 104 animals belonging to nine South American and Spanish breeds was analysed using SSCP and DNA sequencing. The variability found both within and between breeds was very high. There were 61 polymorphic positions, representing 16% of the total sequence obtained. The mean divergence between a pair of sequences was 2.8%. Argentinean Creole horses shared two haplotypes with the Peruvian Paso from Argentina, and the commonest haplotype of the Creole horses is identical to one of the Andalusian horses. Even when there was substantial subdivision between breeds with highly significant Wright's Fixation Index (FST), the parsimony and distance-based phylogenetic analyses failed to show monophyletic groups and there was no clear relationship in the trees between the South American and any of the other horses analysed. Although this result could be interpreted as mixed ancestry of the South American breeds with respect to the Spanish breeds, it is probably indicating the retention of very ancient maternal lineages in the breeds analysed.     

Genetic analysis of the Hispano-Breton heavy horse

Hispano‐Breton (HB) is a horse breed with a recent mixed ancestry. It was developed in the 1930s by crossing local mares with Breton draught horses imported from France. Nowadays it is considered to be in a vulnerable situation due to census decline. To genetically characterize the breed and to set up the basis for a conservation programme, we have employed two types of molecular markers: a 347‐bp D‐loop mitochondrial DNA (mtDNA) fragment and 13 microsatellite loci. A representative sample of 53 HB individuals was analysed together with a sample of 40 Pura Raza Española horses for comparison. Both types of markers revealed a high level of genetic diversity in the HB breed, emphasizing the importance of its conservation. The construction of a phylogenetic network with mtDNA sequences including various Iberian breeds and European heavy horses provided an overall picture of the ubiquitous appearance of HB matrilines with respect to other breeds and revealed the singularity of certain HB maternal lineages. Despite the high allelic richness found in HB horses, microsatellite analysis evidenced a certain degree of inbreeding as a consequence of the type of management generally used for local breeds.     

Molecular phylogeny of Indian horse breeds with special reference to Manipuri pony based on mitochondrial D-loop

Manipuri pony is the geographically distant breed of horse from the five recognized horse breeds found in the Indian subcontinent. The phylogenetic relationship of Manipuri pony with the other breeds is unknown. The diversity in the mitochondrial (mt) DNA D-loop region is employed as an important tool to understand the origin and genetic diversification of domestic horses and to examine genetic relationships among breeds around the world. This study was carried out to understand the maternal lineages of Manipuri pony using the 247 bp region of the mtDNA D-loop. The dataset comprised of eleven numbers of self developed sequences of Manipuri pony, 59 and 35 number of retrieved sequences of Indian horse breeds and other worldwide breeds respectively. A total of 35 haplotypes was identified with a high level of genetic diversity in the Indian breeds. A total of seven major mtDNA haplogroups (A–G) was identified in the Indian horse breeds that indicated the abundance of mtDNA diversity and multiple origins of maternal lineages in them. The majority of the studied sequences of Indian breeds (33.3 %) were grouped into haplogroup D and least (3.9 %) in haplogroup E. The Manipuri breed showed the least F_ST distance (0.03866) with the most diverged Indian breeds and with Thoroughbred horse among the worldwide. This study indicated a close association between Manipuri pony and Thoroughbred.     

Maternal and paternal genetic variation in Estonian local horse breeds in the context of geographically adjacent and distant Eurasian breeds

The maternal and paternal genetic variation of horse breeds from the Baltic Sea region, including three local Estonian breeds, was assessed and compared with that of Altai and Yakutian horses. In the mtDNA D‐loop region, 72 haplotypes assigned to 20 haplogroups in the nine breeds were detected. In Estonian local breeds, 38 mtDNA haplotypes were found, and five of them were shared by the three breeds. More than 60% of all identified haplotypes were rare. Compared with the Estonian Native and Estonian Heavy Draught breeds, a higher haplotypic diversity was found in the Tori breed (h = 0.969). Moreover, four haplotypes shared among Finnish and Estonian local horse breeds indicated ancient ancestry, and of these, H30 (haplogroup D3) showed global sharing and genetic links between modern Baltic Sea region and Siberian horses, specifically. The studied breed set showed high variability in maternal inheritance and mixed patterns of the international and native breeds of the Siberian and Baltic regions. No variation was found in paternally inherited markers among horse breeds in the Baltic Sea region.     

Characterization of cytochrome b diversity in Chinese domestic horses

Previous mitochondrial DNA (mtDNA) D‐loop and microsatellite studies have shown that Chinese horses have multiple maternal origins and high genetic diversity. To better characterize maternal genetic origins and diversity of Chinese domestic horses, we conducted a comprehensive analysis of 407 complete 1140 bp sequences of the horse mitochondrially encoded cytochrome b (CYTB) gene, including 323 horses from 13 Chinese indigenous breeds and 84 reference sequences from GenBank. A total of 114 haplotypes were identified, of which 73 appeared among the 13 Chinese horse breeds. The high mitochondrially encoded cytochrome b haplotypic diversity suggests multiple maternal origins in Chinese horses.     

Genetic analysis of seven Italian horse breeds based on mitochondrial DNA D‐loop variation

To understand the origin and genetic diversity of Italian horses, mitochondrial DNA D‐loop sequences were generated for 163 horses from seven breeds. Sequence analysis of a 480‐bp segment revealed a total of 84 haplotypes with 57 polymorphic sites, indicating multiple maternal origins and high genetic diversity. Comparison of the haplotypes with the equine mtDNA haplotype/haplogroup nomenclature showed a haplogroup distribution in the Italian breeds more similar to that found in the Middle East breeds than in the European breeds, probably due to the economic and cultural relationship with the Middle East in the past centuries.     

Multiple genotypes of mitochondrial DNA within a horse population from a small region in Yunnan province of China

mtDNA genotypes of six domestic horses (three adult short horses whose heights are under 1 m and three common domestic horses) from a small region of 15 km² in Malipo county of Yunnan province of China were investigated by the technique of restriction fragment length polymorphism (RFLP) with 16 restriction endonucleases which recognize 6-bp sequences. An average of 56 fragments for an individual was obtained. Unlike other domestic animals, this population of horses exhibits high mtDNA genetic diversity. Each of the six horses has a specific mtDNA genotype showing a pattern of multiple maternal origins, as suggested by fossil and literature records. We think the population of horses is an amazing seed-resource pool of horses and hence deserves to be paid more attention from the view of conservation genetics. However, it is also remarkable that we did not find any typical mtDNA genetic markers which would discriminate between short horses and common domestic horses.     

Mitochondrial D-loop sequence variation among Italian horse breeds

The genetic variability of the mitochondrial D-loop DNA sequence in seven horse breeds bred in Italy (Giara, Haflinger, Italian trotter, Lipizzan, Maremmano, Thoroughbred and Sarcidano) was analysed. Five unrelated horses were chosen in each breed and twenty-two haplotypes were identified. The sequences obtained were aligned and compared with a reference sequence and with 27 mtDNA D-loop sequences selected in the GenBank database, representing Spanish, Portuguese, North African, wild horses and an Equus asinus sequence as the outgroup. Kimura two-parameter distances were calculated and a cluster analysis using the Neighbour-joining method was performed to obtain phylogenetic trees among breeds bred in Italy and among Italian and foreign breeds. The cluster analysis indicates that all the breeds but Giara are divided in the two trees, and no clear relationships were revealed between Italian populations and the other breeds. These results could be interpreted as showing the mixed origin of breeds bred in Italy and probably indicate the presence of many ancient maternal lineages with high diversity in mtDNA sequences.     

Historical demographic profiles and genetic variation of the East African Butana and Kenana indigenous dairy zebu cattle

Butana and Kenana breeds from Sudan are part of the East African zebu Bos indicus type of cattle. Unlike other indigenous zebu cattle in Africa, they are unique due to their reputation for high milk production and are regarded as dairy cattle, the only ones of their kind on the African continent. In this study, we sequenced the complete mitochondrial DNA (mtDNA) D‐loop of 70 animals to understand the maternal genetic variation, demographic profiles and history of the two breeds in relation to the history of cattle pastoralism on the African continent. Only taurine mtDNA sequences were identified. We found very high mtDNA diversity but low level of maternal genetic structure within and between the two breeds. Bayesian coalescent‐based analysis revealed different historical and demographic profiles for the two breeds, with an earlier population expansion in the Butana vis a vis the Kenana. The maternal ancestral populations of the two breeds may have diverged prior to their introduction into the African continent, with first the arrival of the ancestral Butana population. We also reveal distinct demographic history between the two breeds with the Butana showing a decline in its effective population size (N_e) in the recent past ~590 years. Our results provide new insights on the early history of cattle pastoralism in Sudan indicative of a large ancient effective population size.     

Genetic characterization of southwestern European bovine breeds: a historical and biogeographical reassessment with a set of 16 microsatellites

The origin of Iberian cattle has been suggested by some authors to be the product of European and north African cattle entrances during the last few thousands of years. However, these hypotheses were mainly based on morphological similarities. This study analyzed 889 unrelated individuals from 15 representative Iberian breeds and 3 French breeds for 16 microsatellite loci. Statistical tests were used to calculate interpopulation genetic distances (D(A)) and principal components analysis (PCA). To visualize the geographical distribution of the genetic differentiation between Iberian cattle breeds, data from the PCA analysis were used to construct synthetic maps. Genetic similarity among neighboring Iberian breeds is mainly caused by gene flow. However, recent demographic fluctuations and reproductive isolation in Alistana, Mirandesa, and Tudanca has increased genetic drift, which may be the main cause for the relatively high differentiation of these populations. The synthetic maps constructed with the first and second PCs revealed (1) a large differentiation between Northern Iberian breeds rather than between more geographically distant breeds, and (2) a clear east-west gradient that may be related with the model of demic diffusion of agriculture. Finally, we detected no strong evidence for an African genetic influence in the Iberian cattle breeds analyzed in this study.     

European domestic horses originated in two holocene refugia

The role of European wild horses in horse domestication is poorly understood. While the fossil record for wild horses in Europe prior to horse domestication is scarce, there have been suggestions that wild populations from various European regions might have contributed to the gene pool of domestic horses. To distinguish between regions where domestic populations are mainly descended from local wild stock and those where horses were largely imported, we investigated patterns of genetic diversity in 24 European horse breeds typed at 12 microsatellite loci. The distribution of high levels of genetic diversity in Europe coincides with the distribution of predominantly open landscapes prior to domestication, as suggested by simulation-based vegetation reconstructions, with breeds from Iberia and the Caspian Sea region having significantly higher genetic diversity than breeds from central Europe and the UK, which were largely forested at the time the first domestic horses appear there. Our results suggest that not only the Eastern steppes, but also the Iberian Peninsula provided refugia for wild horses in the Holocene, and that the genetic contribution of these wild populations to local domestic stock may have been considerable. In contrast, the consistently low levels of diversity in central Europe and the UK suggest that domestic horses in these regions largely derive from horses that were imported from the Eastern refugium, the Iberian refugium, or both.     

Phylogenetic relationships of Cheju horses to other horse breeds as determined by mtDNA D-loop sequence polymorphism

Historical records suggest that horses inhabiting the island of Cheju in Korea are descendants of Mongolian horses introduced in 1276. Other studies, however, suggest that horses may have been present on the island prior to the Mongolian introduction. To determine the origin of the Cheju horses we used a phylogenetic analysis of sequences of the mitochondrial DNA (mtDNA) D-loop region, including tRNA Pro and parts of tRNA thr and tRNA Phe sequences (1102-bp excluding the tandem repeat region). Maximum parsimony and neighbor-joining trees were constructed using sequences determined for seven Cheju, four Mongolian, one Przewalskii and two Chinese Yunnan horses, and published sequences for one Swedish and three Thoroughbred horses. Donkey mtDNA was used as an outgroup. We found that the mtDNA D-loop sequence varies considerably within Mongolian, Cheju and Thoroughbred horse breeds, and that Cheju horses clustered with Mongolian horses as well as with horses from other distantly related breeds. On the basis of these findings we propose that horses on Cheju Island are of mixed origin in their maternal lineage, and that horses may have existed and been traded on the island before the Mongolian introduction.     

History of Lipizzan horse maternal lines as revealed by mtDNA analysis

Sequencing of the mtDNA control region (385 or 695 bp) of 212 Lipizzans from eight studs revealed 37 haplotypes. Distribution of haplotypes among studs was biased, including many private haplotypes but only one haplotype was present in all the studs. According to historical data, numerous Lipizzan maternal lines originating from founder mares of different breeds have been established during the breed''s history, so the broad genetic base of the Lipizzan maternal lines was expected. A comparison of Lipizzan sequences with 136 sequences of domestic- and wild-horses from GenBank showed a clustering of Lipizzan haplotypes in the majority of haplotype subgroups present in other domestic horses. We assume that haplotypes identical to haplotypes of early domesticated horses can be found in several Lipizzan maternal lines as well as in other breeds. Therefore, domestic horses could arise either from a single large population or from several populations provided there were strong migrations during the early phase after domestication. A comparison of Lipizzan haplotypes with 56 maternal lines (according to the pedigrees) showed a disagreement of biological parentage with pedigree data for at least 11% of the Lipizzans. A distribution of haplotype-frequencies was unequal (0.2%–26%), mainly due to pedigree errors and haplotype sharing among founder mares.     

Mitochondrial DNA D‐loop sequence variation in maternal lineages of Iranian native horses

To understand the origin and genetic diversity of Iranian native horses, mitochondrial DNA (mtDNA) D‐loop sequences were generated for 95 horses from five breeds sampled in eight geographical locations in Iran. Sequence analysis of a 247‐bp segment revealed a total of 27 haplotypes with 38 polymorphic sites. Twelve of 19 mtDNA haplogroups were identified in the samples. The most common haplotypes were found within haplogroup X2. Within‐population haplotype and nucleotide diversities of the five breeds ranged from 0.838 ± 0.056 to 0.974 ± 0.022 and 0.011 ± 0.002 to 0.021 ± 0.001 respectively, indicating a relatively high genetic diversity in Iranian horses. The identification of several ancient sequences common between the breeds suggests that the lineage of the majority of Iranian horse breeds is old and obviously originated from a vast number of mares. We found in all native Iranian horse breeds lineages of the haplogroups D and K, which is concordant with the previous findings of Asian origins of these haplogroups. The presence of haplotypes E and K in our study also is consistent with a geographical west–east direction of increasing frequency of these haplotypes and a genetic fusion in Iranian horse breeds.     

Genetic diversity of Cheju horses (Equus caballus) determined by using mitochondrial DNA D-loop polymorphism

We used sequence polymorphism of the mitochondrial DNA D-loop (968 bp excluding the tandem repeat region) to determine genetic diversity of horses inhabiting Cheju (a southern island of Korea). Seventeen haplotypes with frequencies from 1.5 to 21.5% were found among 65 Cheju horse samples. Genetic diversity (h) of the 17 haplotypes was calculated to be 0.91, indicating that the extant Cheju horse population consists of diverse genetic groups in their maternal lineage. Phylogenetic analysis showed that 17 types of Cheju (D-loop sequences determined), 5 Mongolian, 6 Arabian, 3 Belgian, 2 Tsushima, 2 Yunnan, 1 Przewalskii, and 3 Thoroughbred horses (published sequences for the latter seven breeds) showed that Cheju horses were distributed into many different clusters in the tree. Four Mongolian horses clustered with separate Cheju horse groups, showing that some Cheju horses are clearly of Mongolian origin. The analysis of partial sequences (284 bp) of the D-loop of 109 horses showed that Thoroughbred, Mongolian, Lipizzan, and Arabian breeds are as diverse as Cheju horses. Our data together with others' suggest that most horse breeds tested with reasonably sufficient numbers of samples are diverse in their maternal lineages and also are not uniquely different from each other.     

Extensive mtDNA diversity in horses revealed by PCR–SSCP analysis

The hypervariable D-loop region of mitochondrial DNA (mtDNA) was amplified with the polymerase chain reaction using total horse DNA samples. Analysis of single strand conformation polymorphism (SSCP) of denatured amplification products was carried out by native polyacrylamide (8%) gel electrophoresis followed by silver staining. As many as 15 distinct SSCP variants were revealed when screening a total of 78 maternally unrelated horses representing five different breeds. All breeds showed a high degree of polymorphism and the estimated probability (PI_mt) that two maternally unrelated individuals have, by chance, identical SSCP variants varied between 0.14 and 0.30. We detected no heteroplasmy or deviations from strict and stable maternal inheritance when examining four maternal lineages, each represented by six to eight horses, separated by up to five generations from a common ancestral mare. The study establishes a simple screening method for detecting equine mtDNA types, which can be applied for tracing maternal genealogies and for association studies.     

Mitochondrial diversity and the origin of Iberian sheep

Mitochondrial DNA diversity was analysed in 19 Iberian and six foreign sheep breeds. Three mtDNA lineages (B, A and C) were found in the Iberian sheep, with type B clearly predominating over the others. The results were analysed for each of the morphologically determined breed groups in Iberian sheep: Merino, Entrefino, Churro and Iberian trunks. MtDNA lineage C was found only in the Iberian trunk composed of Montesina and Ojalada. These two populations had high mtDNA variability, and in the Iberian sheep only Merino Branco had more variation. The other three Merino types studied showed moderate variability, including the most authentic Merino, the Spanish Merino. These three Merinos clustered closely in a multidimensional scaling representation of distances, while the fourth breed (Merino Branco) showed a clear separation. As for the other two trunks, breeds from the Churro group showed greater maternal uniformity while results for populations included in the so-called Entrefino trunk seemed to have a more heterogeneous maternal origin. The results obtained are discussed with available data from nuclear markers and with morphological classifications, and all this information is analysed in relation to the origin of the different Iberian sheep breeds.