Summary Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN “speed gene,” only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

There are three native Portuguese horse breeds: Lusitano, Sorraia and Garrano. This study compares diversity patterns of 17 protein and 12 microsatellite markers in these three as well as 30 other breeds to infer relationships among the breeds and to compare levels of polymorphism of these breeds for use in conservation efforts. The Garrano and the Lusitano showed a high level of genetic diversity, similar to that observed for most of the other analysed breeds, while the Sorraia and Friesian breeds showed low levels of variation for both genetic marker types. The combined protein and microsatellite data produced a tree that fit historical records well and with greater confidence levels than those for either data set alone. The combined genetic diversity and relationship information provides important baseline data for future breed conservation efforts, especially for a critically endangered breed such as the Sorraia.

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.

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.

The Sorraia horse is a closed breed with reduced effective population size and considered in critical maintained risk status. The breed exists in 2 main breeding populations, one in Portugal and one in Germany, with a smaller population size. A set of 22 microsatellite loci was used to examine genetic diversity and structure of the Sorraia horse breed and to compare individual inbreeding coefficient F, estimated from pedigree data, with individual heterozygosity and mean d(2). The Sorraia horse shows lower levels of microsatellite diversity when compared with other horse breeds. Due to management strategies, there are clear differences in the genetic structure of the 2 main Sorraia horse populations. Individual heterozygosity was shown to be a good estimator, used together with or as an alternative to inbreeding coefficient, in predicting fitness and evaluating the inbreeding level of the Sorraia horse. The information gathered in this study, combined with information available from previous studies, offers an important and wide information base for the future development of an effective breeding management of the Sorraia horse in order to preserve this endangered breed.

The mitochondrial control region sequence was compared between the two surviving matrilineal lineages of the Sorraia horse. Two haplotypes were found, showing three differences between them. The reduced number of surviving maternal lineages emphasizes the importance of establishing a conservation plan for this endangered breed.

The present study intends to survey the genetic variability of an endangered semiferal Portuguese native pony breed, the Garrano. Thirteen microsatellite markers were examined in 277 animals born in 1998, belonging to eight subpopulations corresponding to eight northern Portuguese geographic regions. Mean heterozygosity (H(o)) in the Garrano breed was 0.732, ranging from 0.531 to 0.857 across subpopulations. Allelic frequencies and diversity differed significantly between regions, suggesting the existence of genetic differentiation within the breed confirmed by the population differentiation estimator F(ST). Allele sharing genetic distance (Dps) was used to determine the relationship between the analyzed subpopulations, some of which are diverging significantly from the others. Relationships among six Iberian horse breeds (including the Garrano) were assessed through the restricted maximum likelihood method, which clusters the Garrano with another Celtic pony, the Spanish Asturcon, traditionally bred in geographical continuity with the Garrano. Results reveal a statistically significant deficit of heterozygotes within the Garrano breed (F(IT) = 0.031, p < 0.05). Regular monitoring should therefore be implemented so that the effect of genetic drift within subpopulations, enhanced by inbreeding, may be successfully minimized.

The use of DNA technology for parentage testing is increasing every day. Most laboratories have improved their exclusion probabilities (PE) by the addition of DNA microsatellite loci to standard blood-typing results. The efficacy of each locus depends on the number of alleles detected and their frequencies in the breeds tested. Here we analyzed the usefulness of six microsatellites for routine parentage testing in three Portuguese autochthonous horse breeds: Lusitano, Sorraia and Garrano. The DNA loci analyzed -ASB2, HMS3, HMS7, HTG4, HTG10, and VHL20 -were chosen based upon the polymorphism detected in other breeds with Iberian horse influence. The estimated probability of exclusion of wrongly named parents (PE) was high, with values ranging from 88.5% to 99.6%.