This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity.
several rare populations are included in this dataset. The Caspian is one of the oldest breeds in the Middle East and was thought to be extinct until its recent rediscovery in 1965. The Florida Cracker, a now rare breed, was developed in the United States from feral stock of Iberian descent [57]. The sample size of the Florida Cracker limits the conclusions that can be drawn regarding within-breed diversity. However, the Caspian shows high Ne, He, and estimates of f, given its rarity. After rediscovery of the breed, which historically was believed to represent a type of landrace population, [58] describes a three-year survey, which found approximately 50 individuals remaining, noting that many could not be considered “pure.” In addition, [33] were unable to show evidence of a recent bottleneck in the Caspian breed. The diversity observed in what are now considered Caspian horses likely stems from high levels of diversity within those individuals that founded the modern population.
PLoS One. 2013; 8(1): e54997.
Published online 2013 Jan 30. doi: 10.1371/journal.pone.0054997
PMCID: PMC3559798
PMID: 23383025
Jessica L. Petersen, 1 , * James R. Mickelson, 1 E. Gus Cothran, 2 Lisa S. Andersson, 3 Jeanette Axelsson, 3 Ernie Bailey, 4 Danika Bannasch, 5 Matthew M. Binns, 6 Alexandre S. Borges, 7 Pieter Brama, 8 Artur da Câmara Machado, 9 Ottmar Distl, 10 Michela Felicetti, 11 Laura Fox-Clipsham, 12 Kathryn T. Graves, 4 Gérard Guérin, 13 Bianca Haase, 14 Telhisa Hasegawa, 15 Karin Hemmann, 16 Emmeline W. Hill, 17 Tosso Leeb, 18 Gabriella Lindgren, 3 Hannes Lohi, 16 Maria Susana Lopes, 9 Beatrice A. McGivney, 17 Sofia Mikko, 3 Nicholas Orr, 19 M. Cecilia T Penedo, 5 Richard J. Piercy, 20 Marja Raekallio, 16 Stefan Rieder, 21 Knut H. Røed, 22 Maurizio Silvestrelli, 11 June Swinburne, 12 , 23 Teruaki Tozaki, 24 Mark Vaudin, 12 Claire M. Wade, 14 and Molly E. McCue 1
Hans Ellegren, Editor
Abstract
Horses were domesticated from the Eurasian steppes 5,000–6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. FST calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.
Introduction
With a world-wide population greater than 58 million [1], and as many as 500 different breeds, horses are economically important and popular animals for agriculture, transportation, and recreation. The diversity of the modern horse has its roots in the process of domestication which began 5,000–6,000 years ago in the Eurasian Steppe [2]–[4]. Unlike other agricultural species such as sheep [5] and pigs [6], [7], archaeological and genetic evidence suggests that multiple horse domestication events occurred across Eurasia [2], [8]–[12]. During the domestication process, it is believed that gene flow continued between domesticated and wild horses [13] as is likely to also have been the case during domestication of cattle [14], [15]. Concurrent gene flow between domestic and wild horses would be expected to allow newly domestic stock to maintain a larger extent of genetic diversity than if domestication occurred in one or few events with limited individuals.
Prior genetic work aimed at understanding horse domestication has shown that a significant proportion of the diversity observed in modern maternal lineages was present at the time of domestication [2], [8], [16]. The question of mitochondrial DNA (mtDNA) diversity was further addressed by recent sequencing of the entire mtDNA genome. These studies estimate that, minimally, 17 to 46 maternal lineages were used in the founding of the modern horse [2], [17]; however, those data were unable to support prior studies suggesting geographic structure among maternal lineages [9], [18]. Recent nuclear DNA analyses have utilized “non-breed” horses sampled across Eurasia to attempt to understand the population history of the horse. These microsatellite-based studies suggest a weak pattern of isolation by distance with higher levels of diversity in, and population expansion originating from Eastern Asia [13], [19]. High diversity as observed by both mtDNA and microsatellites and the absence of strong geographical patterns is likely a result of continued gene flow during domestication, the high mobility of the horse, and its prevalent use for transportation during and after the time of domestication. Interestingly, while significant diversity is observed in maternal lineages, paternal input into modern horse breeds appears to have been extremely limited as shown by a lack of variation at the Y-chromosome [20], [21].
Diversity in the founding population of the domestic horse has since been exploited to develop a wealth of specialized populations or breeds. While some breeds have been experiencing artificial selection for hundreds of years (e.g. Thoroughbred, Arabian), in general, most modern horse breeds have been developed recently (e.g. Quarter Horse, Paint, Tennessee Walking Horse) and continue to evolve based upon selective pressures for performance and phenotype (Table 1). Horse breeds resulting from these evolutionary processes are generally closed populations consisting of individual animals demonstrating specific phenotypes and/or bloodlines. Each breed is governed by an independent set of regulations dictated by the respective breed association. Not all breeds are closed populations. Some breed registries allow admixture from outside breeds (e.g. Swiss Warmblood, Quarter Horse), and others are defined by phenotype (e.g. Miniature). Finally, some populations that are often referred to as breeds are classified simply by their geographic region of origin and may not be actively maintained by a formal registry (e.g. Mongolian, Tuva) (Table 1). Those breeds that may be free ranging and experience lesser degrees of management may more appropriately be termed “landrace populations.” Therefore, genetic characteristics within horse breeds are expected to differ based upon differences in the definition of the breed, the diversity of founding stock, the time since breed establishment, and the selective pressures invoked by breeders. The extent of gene flow not only varies within breed, but among horse breeds, the direction and level of gene flow is influenced by breed restrictions/requirements, and potentially by geographic distance.
Considering modern breeds, unlike mtDNA, nuclear markers can discern breed membership [12]. However, studies of nuclear genetic diversity of modern breeds to date have most commonly focused on a single population of interest, sets of historically related breeds, or breeds within a specific geographic region [22]–[36]. Additionally, these analyses of nuclear genetic diversity in horse breeds are largely based upon microsatellite loci, which do not often permit consolidation of data across studies. Thus, large, across-breed investigations of nuclear diversity in the modern, domestic horse are lacking.
The Equine Genetic Diversity Consortium (EGDC), an international collaboration of the equine scientific community, was established in an effort to quantify nuclear diversity and the relationships within and among horse populations on a genome-wide scale. The development of this consortium has facilitated the collection of samples from 36 breeds for genotyping on the Illumina 50K SNP Beadchip. The breeds included in this report represent many of the most popular breeds in the world as well as divergent phenotypic classes, different geographic regions of derivation, and varying histories of breed origin (Table 1). The standardized SNP genotyping platform permits the compilation of data across breeds at a level never before achieved. Results of this collaboration now allow for the detailed description of diversity and assessment of the effects of genetic isolation, inbreeding, and selection within breeds, and the description of relationships among breeds. These data will also facilitate future across breed genome-wide association studies as well as investigations into genomic targets of selection.
Results
Samples
Of the 38 populations sampled, two breeds were represented by geographically distinct populations: the Thoroughbred was sampled in the both the United States (US) and the United Kingdom and Ireland (UK/Ire), and the Standardbred was sampled in the US as well as in Norway. Eight Standardbred horses sampled from the US were noted to be pacing horses as opposed to the Norwegian and remaining US individuals that were classified as trotters. In addition, the International Andalusian and Lusitano Horse Association Registry (IALHA) in the US maintains one stud book but designates whether the individual was derived from Spanish (Pura Raza Española) or Portuguese (Lusitano) bloodlines, or a combination of both. Of the Andalusian horses collected in the US, five were noted to have Portuguese bloodlines.
Phenotypic classifications of the horse breeds include those characterized by small stature (Miniature Horse, pony breeds), breeds characterized by large stature and/or large muscle mass in proportion to size (draft breeds), light horse or riding breeds, gaited breeds, rare breeds, breeds founded in the past 80 years, and populations that are relatively unmanaged (“landrace”). The number of samples, sampling location, region of breed origin, and a list of primary breed characteristics are found in Table 1.
After pruning of individuals for genotyping quality and relationships (see methods), and keeping a similar number of individuals per breed, 814 of the 1,060 horses remained in the analysis. Of the horses removed, 12 had known pedigree relationships at or more recent to the grandsire/dam level, 44 individuals were removed at random from overrepresented breeds to equalize sample size across breeds, 4 failed to genotype at a rate greater than 0.90, and 186 were removed due to pi hat values (pairwise estimates of identity by descent) above the allowed threshold. Of those last 186 horses that were removed, 122 were from disease studies where relationships were common due to sampling bias.
Within Breed Diversity
Diversity indices were calculated using 10,536 autosomal SNPs that remained after pruning for minor allele frequency (MAF), genotyping rate, and linkage disequilibrium (LD) across breeds (referred to as the primary SNP set). Diversity indices were also calculated using three other SNP sets, resulting from different levels of LD-based pruning (see methods). Individuals noted as outliers in parsimony and cluster analyses (see below) were excluded from within-breed diversity calculations.
Using the primary SNP set, diversity, as measured by expected heterozygosity (He), ranged from 0.232 in the Clydesdale, to 0.311 in the Tuva (Table 2). Considering the SNP sets pruned less stringently for LD, the diversity within the Thoroughbred increased in relationship to the other breeds, as did that of nine other breeds. Mean and total heterozygosity increased with increased number of loci and less stringent LD pruning (Table 2). Inbreeding coefficients (FIS) calculated on the primary SNP set showed significant excess homozygosity in 17 populations, which was greatest in the Andalusian (0.065). Three of the four lowest FIS values were found in the Thoroughbred samples (UK/Ire, US, and when considered together) (Table 2).
Table 2
Number of samples (N), effective population size (Ne), individual inbreeding estimates (f), inbreeding coefficient (FIS), and expected heterozygosity (He) from four SNP sets pruned based upon varying levels of LD.
aIndividuals from this breed also included in [41];
b20 of these individuals were also reported in [41];
c17 of these individuals were also reported in [41];
d21 of these individuals were also reported in [41];
e19 of these individuals were also reported in [41].
FIS and f were calculated based upon the primary SNP set (10,536 loci). Samples also used in [41] are indicated in the footnotes.
*indicates significance at α<0.05 determined by 10,000 permutations.
Inbreeding coefficients (f) calculated for each individual based upon observed and expected heterozygosity showed several individuals with significant loss of heterozygosity. The highest individual value of f (0.56) was found in an Exmoor pony. Within breeds, average individual estimates of f were greatest in the Clydesdale, Mangalarga Paulista, and Exmoor while the lowest breed means were found in the landrace populations (Table 2).
Effective population size (Ne), as estimated by LD [37] using an autosomal SNP set pruned within each breed for quality, was lowest (143) in the UK/Ire sample of the Thoroughbred (UK/Ire) but also low in the other racing breeds as well as the Clydesdale (Table 2). Highest values of Ne were observed in the Eurasian landrace populations, the Mongolian (743) and Tuva (533), and also in the Icelandic (555), Finnhorse (575), and Miniature (521). Breed-specific decay of LD essentially mirrors the results of the Ne calculation given the relationship between the statistics. A plot of LD across 2 Mb in a subset of the breeds that represent the range of Ne estimates is found in Figure S1.
Parsimony and Principal Component Analyses
With a domestic ass designated as the outgroup, parsimony analysis of 10,066 loci pruned for LD of R2 = 0.2 (see methods) resulted in generally tight clustering and monophyly of samples within breeds, supported by high bootstrap values (Figure 1). Major clades of the tree show grouping of the Iberian breeds (Lusitano and Andalusian), ponies (Icelandic, Shetland, Miniature), Scandinavian breeds (Finnhorse, North Swedish Horse, Norwegian Fjord), heavy draft horses (Clydesdale, Shire, Belgian, Percheron), breeds recently admixed with and/or partly derived from the Thoroughbred (Paint, Quarter Horse, Maremmano, Swiss Warmblood, Hanoverian), modern US breeds (American Saddlebred (hereafter “Saddlebred” and Tennessee Walking Horse), trotting breeds (Standardbred and French Trotter), and Middle Eastern breeds (Akhal Teke and Arabian). Exceptions to monophyly include the Paint and Quarter Horse as well as the Hanoverian and Swiss Warmblood, which are mixed in clades surrounding the Thoroughbred and Maremmano. In addition, the Clydesdale was placed as a clade within the Shire breed and the Shetlands as a clade within the Miniatures. Strong bootstrap support for monophyly is present within a subset each of Lusitanos (83%), and Andalusians (87%); however the remainder of individuals from these breeds were intermixed. No structure was found within the US sample regarding individual Andalusians noted to have Portuguese bloodlines opposed to those with Spanish bloodlines (Figure S2). The Mongolian and most Tuva horses were grouped together while a subset of the Tuvas fell out as a sister clade to the Caspians. Several individuals were not positioned in the clades that represented the majority of the other individuals in the breed (Figure 1). These include three Shires, two Mongolians, a Caspian, and a Norwegian Fjord. In each instance, the outlier status of these individuals was also supported by cluster analysis (see below).
Figure 1
Individual and breed relationships among 814 horses illustrated by parsimony.
Parsimony tree created from 10,066 SNPs and rooted by the domestic ass. Breeds are listed in the legend in order starting from the root and working counterclockwise. Individual outliers with respect to their breeds are noted with arrows. Bootstrap support calculated from 1,000 replicates is shown for major branches when greater than 50%.
Principal component analysis (PCA) also serves to visualize individual relationships within and among breeds. The plot of PC1 vs. PC2 shown in Figure S3 illustrates relationships similar to those shown by parsimony, including the placement of outliers outside of their respective breeds. All Thoroughbred samples, regardless of origin, are separated from the others by PC1 and form a cluster at the top of the figure. Intermediate between the Thoroughbred and central cluster of breeds are the Hanoverian, Swiss Warmblood, Paint, and Quarter Horse. The Shetland, Icelandic, and Miniature split from the remainder of samples in PC2, falling out in the lower left corner, and the British drafts anchor the figure at the lower right. While most breeds cluster tightly, several are dispersed across one or both PCs. The Hanoverian, Swiss Warmblood, Paint, and Quarter Horse, as noted above, are extended along PC1, while the Arabian and Franches-Montagnes show similar spreading, also along PC1. The Tuva, Clydesdale, and Shire individuals also are not as tightly clustered as other populations despite the low within breed diversity of the latter two.
Distance Analysis
An unrooted neighbor joining (NJ) tree of Nei’s distance [38] was constructed using SNP frequencies within breeds from the 10,536 SNP data set (Figure 2). The relative placement of breeds reflects that seen in the parsimony tree with several exceptions. The Paint, Quarter Horse, Swiss Warmblood, Hanoverian, Maremmano, and Thoroughbred, are found in one large branch of the tree, although the Maremmano is placed outside of the clade containing the aforementioned breeds. The position of the Morgan with the Saddlebred and Tennessee Walking Horse also deviates from parsimony analysis but reflects historic records of relationships among these breeds. The Scandinavian breeds remain in one branch of the clade, which also includes the Shetland and Miniature. Unlike the parsimony cladogram, the Caspian falls in a clade with the other Middle Eastern breeds, the Arabian and Akhal Teke. Finally, the Exmoor, a British breed, is placed with another British breed, the New Forest Pony, rather than with the Scandinavian breeds as in the parsimony analysis. Each branch shows support of over 50%, with many clades being supported by over 99% of the 1,000 bootstrap replicates.
Figure 2
Distance based, neighbor joining tree calculated from SNP frequencies in 38 horse populations.
Majority rule, neighbor joining tree created from 10,536 SNP makers using Nei’s genetic distance and allele frequencies within each population. Percent bootstrap support for all branches calculated from 1,000 replicates is shown.
Cluster analysis
Likelihood scores for runs of various K in Structure showed an increase in overall mean ln P(X|K) until K = 35 (Figure S4). A clear “true” value of K is not obvious examining the likelihood scores or using the Evanno method [39] (data not shown); however, variance among runs begins to increase with a diminishing increase in likelihood scores after K = 29, which is near the peak of the curve. The value of the highest proportion (breed average q-value) of assignment of each breed for each value of K, as well as the cluster to which it assigns is shown in Table S1. Additionally, the proportion assignment at K = 29 for each of the breeds is found in Table S2.
The first breeds to have all individuals assign strongly to one cluster are the Thoroughbred and Clydesdale (with Shire) at K = 2, followed by the Shetland at K = 3; these four breeds do not show signs of admixture at any K value analyzed. Evidence of weak geographic grouping is observed at K = 4, which consists of: 1, the Middle Eastern and Iberian breeds (pink); 2, the Thoroughbred and breeds to which it continues to be or was historically crossed (yellow); 3, breeds developed in Scandinavia and Northern Europe (orange); and 4, the British Isles draft breeds (blue) (Figure 3).
Figure 3
Bayesian clustering output for five values of K in 814 horses of 38 populations.
Structure output for five values of K investigated. Each individual is represented by one vertical line with the proportion of assignment to each cluster shown on the y axis and colored by cluster. Other values of K are shown in Figure S1 and a summary of assignment of each breed in Tables S1 and S2.
Middle Eastern and Iberian Breeds
As also observed in the NJ tree, clustering of the Iberian and Middle Eastern breeds with the Mangalarga Paulista, Peruvian Paso, and Puerto Rican Paso Fino (q >0.5) is observed until K = 8, after which point the Mangalarga Paulista assigns with q = 0.93 to another cluster. The remaining breeds cluster together until K = 12, at which time the Middle Eastern breeds (Arabian, Akhal Teke, and Caspian) are assigned to their own cluster, leaving the Iberian breeds clustered with the Peruvian Paso and Puerto Rican Paso Fino. At low values of K (i.e. K <6) the Florida Cracker, Saddlebred, Standardbreds, Morgan, and Tennessee Walking Horse fall into the cluster with the Iberian and Middle Eastern breeds with breed mean q >0.5. At K = 29, each of these breeds is assigned with q >0.72 to an individual cluster with the exception of the Lusitano and Andalusian, which remaining clustered together.
Thoroughbreds and Thoroughbred Crossed Breeds
Relationships described by the NJ tree among the Thoroughbred, Hanoverian, Swiss Warmblood, Paint, Quarter Horse, and Maremmano are also seen in cluster analysis. Clustering of those breeds with the Thoroughbred is observed throughout the values of K examined although at moderate frequencies (Figure 3, Table S1, Figure S5). At K = 29, the Hanoverian and Swiss Warmblood remain assigned to the cluster defined by the Thoroughbred but with assignment probabilities of 0.51 each. The Quarter Horse and Paint also assign to this cluster with q-values of 0.30 and 0.34, respectively. Neither the Quarter Horse, Paint, Hanoverian, or Swiss Warmblood populations assign to any cluster with q >0.62 at K = 29. No evidence of population substructure is observed between the US and UK/Ire Thoroughbreds as also shown by PCA and parsimony analyses (Figure S6).
Scandinavian and Northern European Breeds
As in the NJ and parsimony trees, the Finnhorse, Icelandic, Miniature, North Swedish Horse, Norwegian Fjord, and Shetland are parsed into the same cluster (q-value >0.5) through K = 5. However, unlike the NJ tree, at K = 4, the highest value of assignment places the Belgian and Percheron into this cluster although with q <0.5 (0.42 and 0.38, respectively). The relationship remains until K = 6, at which time the Miniature, Icelandic, and Shetland fall into a different cluster. At K = 10, the Icelandic clusters again with the North Swedish Horse and Norwegian Fjord. The Norwegian and United States Standardbred populations, which at K = 4 assign with q >0.5 to the cluster containing the Scandinavian breeds, separate from the Scandinavian breeds at K = 5. At K = 31, substructure appears in the Standardbred samples, which correlates to those individuals identified as pacers and that fall into an individual clade in the parsimony tree (Figure S7). At K = 29, the Miniature and Shetland continue to be assigned to the same cluster (q-values = 0.55 and 0.95, respectively). The next highest proportions of assignment of the Miniature horse are to the clusters described by the New Forest Pony (q = 0.20) and Icelandic (q = 0.11). No value of K evaluated eliminated signals of admixture from all populations in the dataset at K = 38 (the actual number of populations sampled) or any value of K through 45 (data not shown).
British Isles Draft
The Clydesdale and Shire cluster together, and apart from the other breeds beginning at K = 3. In addition, the Fell Pony, which is placed within the same clade in the NJ and parsimony trees, and proximal to the Clydesdale and Shire in PCA, shows moderate assignment to this cluster (0.29< q <0.41) for several values of K from 4 to 14. At K = 29, the Shire assigns to the same cluster as the Clydesdale with q = 0.69. The individual outliers from the Shire breed also noted in parsimony analysis are evident beginning at K = 3. Excluding these outliers, at K = 29, the proportion of assignment for the Shires to the cluster with the Clydesdale increases to 0.74.
FST
All pairwise FST values calculated between the 37 populations (excluding the Florida Cracker) were significant as tested by 20,000 permutations (Figure 4). The lowest level of differentiation was found between the Paint and Quarter Horse populations (FST = 0.002), while the greatest divergence was observed between the Clydesdale and Mangalarga Paulista (FST = 0.254). The two Thoroughbred populations had an FST value of 0.004, while the two Standardbred populations had 10-fold greater divergence (FST = 0.020) than the minimum observed value in this dataset; this value is similar to that observed between the Lusitano and Andalusian (0.021). An FST value of 0.006 was identified between the Tuva and Mongolian populations. The global FST value was 0.100. AMOVA computed on the set of 37 samples (excluding the outliers identified in Structure and the Florida Cracker) showed that 10.03% of the variance was accounted for among populations (p = <0.001), 0.53% of the variance was among individuals within populations (p = 0.19), and 89.44% of the variation was within individuals (p<0.001).
Figure 4
Pairwise FST values based upon 10,536 SNPs in 37 horse populations.
Pairwise FST values as calculated in Arlequin using 10,536 autosomal SNPs and significance tested using 20,000 permutations. All pairwise values are significantly different from zero. (individual outliers were removed from this analysis).
Discussion
These data are gathered from populations that represent tremendous diversity in phenotype and breed specialization. With breeds sampled across four continents, the resulting relationships observed largely reflect similarities of geographic origin, documented breed histories, and shared phenotypes. In general the highest within breed diversity was observed in breeds that are recently derived, continue to allow introgression of other populations, those that have a large census population size, and landrace populations that experience a lesser degree of controlled breeding. Not surprisingly, low diversity is observed in breeds with small census size, relatively old breeds with closed populations, and those with documented founder effects, whether due to population bottlenecks or selective breeding.
A total of seven individuals were identified by parsimony and cluster analysis as outliers with respect to the breed to which they were assigned. The pedigrees of these individuals were unknown. Because it is possible these horses were unknowingly crossbred or subject to mishandling in the field or laboratory, they were excluded from the within-breed analyses to avoid potential bias in indices of diversity. In addition, the potential impact of SNP ascertainment bias on diversity calculations must be acknowledged. The reference genome is from a Thoroughbred mare [40] and SNP identification was based upon the reference genome and data from seven other horses representing six breeds. Therefore, SNPs are generally derived to identify modern variation within the Thoroughbred as well as between the Thoroughbred and these other breeds. Thus, the SNPs identified may reflect an upward bias in diversity indices in the Thoroughbred and closely related breeds [41]. It seems that ascertainment bias may have particularly influenced the results when considering the data sets that have an increased number of loci resulting from relaxed LD pruning. These results show an increase in the relative diversity of the Thoroughbred, breeds with which the Thoroughbred continues to actively interbreed, and the other SNP discovery breeds, with respect to other breeds in the study. This is opposite of what may be expected given the high levels of genome-wide LD in the Thoroughbred. Without considering SNP ascertainment bias, it is expected that measured diversity would increase in breeds with low LD more quickly than in those with high LD, due to greater independence of markers in the former breeds. These SNPs, derived largely from the Thoroughbred, are apparently detecting a higher proportion of Thoroughbred-specific, rare variants and it appears that as more loci are included, more of these Thoroughbred-based variants are assayed, resulting in the observed increase in variation in the Thoroughbred, Thoroughbred-influenced breeds, and breeds used in SNP ascertainment.
The majority of the analyses were performed using 10,536 SNP markers pruned across breeds for LD of r2>0.2 as well as MAF of 0.05 or above. Even though additional markers could have been used for analysis, many population-level statistics assume independence of loci. The stringent pruning for LD was therefore undertaken to eliminate bias in the test statistics that may result from substantial breed-specific differences in LD [40], [41]. A truncated data set also helped to make calculations, especially cluster analysis, computationally feasible. On the other hand, diversity indices were calculated after pruning the full data set to r2 = 0.2 and r2 = 0.4 (using pairwise correlation), and one replicate setting the threshold to R2<0.1 (using the variance inflation factor), to examine the effect of allowing for varying levels of LD and therefore varying numbers of loci (see methods).
Within-breed Diversity
Even considering SNP ascertainment, low diversity as measured by He was observed in the Thoroughbred as well as the Standardbred, which both experience high selective pressures and are closed populations. Low diversity was also observed in breeds that have undergone a severe population bottleneck, such as the Exmoor and Clydesdale, and breeds that have small census population sizes, such as the Florida Cracker. Although the Thoroughbred is a large population that is widely distributed on a geographic scale, historic records suggest that one sire is responsible for 95% of the paternal lineages in the breed and as few as 30 females make up 94% of maternal lineages [28]. In addition, the population has been largely closed to outside gene flow since the formation of the first stud book in 1791 [42] and individuals within the breed are subject to selective pressure for racing success; therefore low, within-breed diversity is not at all surprising.
Using LD-based calculations, the estimated Ne for the Thoroughbred was similar to that found in a UK sample [43] and among the lowest of the study set despite the large census population size and geographic distribution of this breed. Individual inbreeding values based upon observed vs. expected homozygosity indicate that individual Thoroughbred horses show signs of inbreeding, with a mean loss of heterozygosity of 16.3%. This value is slightly larger than that found in [28] (13.9%). Using the same SNP array, [44] also showed inbreeding in the Thoroughbred, and specifically an increase in inbreeding over time. The only breeds with higher f values were the Exmoor, Clydesdale, Mangalarga Paulista, and Shire. Despite low individual diversity, FIS values do not show significant inbreeding in either of the Thoroughbred populations as a whole, or in the Norwegian Standardbred although FIS is significant in the US Standardbred population (discussed below).
The Clydesdale and Exmoor, in addition to having high individual estimated coefficients of inbreeding, also show the lowest within-breed diversity observed in the dataset. A lack of diversity in the Clydesdale and another British draft breed, the Shire, is likely a result of a severe population bottleneck observed in most draft breeds with the onset of industrialization and after the conclusion of World War II (WWII) as well as selection for size and color [45], [46]. The Exmoor pony, considered to be one of the purest native breeds of Britain, has been naturally selected for survival in harsh winter conditions on the moors in southwest England [45], [47]. Similar to the draft breeds, the Exmoor population decreased significantly after WWII to approximately 50 individuals, undoubtedly influencing the diversity observed in this study. The effect of low population size and selection is also reflected in extremely high individual estimates of f within some individuals. Finally, the Mangalarga Paulista shows low levels of heterozygosity, and as discussed below, the greatest divergence as measured by pairwise FST of all breeds in the study. While these results could be due to geographic distance between this and other breeds, and/or genetic drift, unfortunately these horses were all sampled from only two farms and likely do not represent the entirety of the diversity present in the breed; therefore we cannot rule out sampling error which would inflate the estimated level of divergence between these individuals and the other breeds and result in a decrease in He. However, a lack of diversity in sampling of the breed would not have an effect on estimates of individual inbreeding coefficients, which were among the highest of the entire data set.
Converse to the above examples, high levels of diversity as measured by both He and Ne, accompanied by low estimates of inbreeding (f and FIS), are observed in the Mongolian, Tuva, and New Forest Pony. The Mongolian and Tuva are unique in that they represent landrace populations that are less managed than the popular breeds of Western Europe and North America; they occupy a diverse range of habitat, have been selected for meat and milk in addition to use in transportation, and originate in the region where domestication was likely to have occurred. The population of Mongolian horses is large and individuals are phenotypically diverse [48]. In 1985, approximately two million Mongolian horses of four different types were estimated to live within the country [45]. The Tuva is not as numerous as the Mongolian but is similar in its purpose and also has high within-breed phenotypic diversity. In addition, it is suggested that the Tuva has experienced outcrossing in order to increase its size and stamina [45] as may also be the case in the Mongolian [49]. Similarly, the New Forest Pony was historically a free-ranging population in Great Britain, but was crossbred until the 1930’s. These traits: old populations, large population size, outcrossing, high phenotypic diversity, and lesser artificial selection/management, result in the high levels of genetic diversity observed. This extent of diversity appears to diminish as populations are restricted by selective pressures into formal breeds.
Other population characteristics are likely the cause of the diversity observed in the Finnhorse, Icelandic, and Miniature. In the case of the Icelandic, the high level of diversity was possibly maintained by a large census population size despite isolation for almost a thousand years and several population bottlenecks due to natural disasters [45]. In the case of the Finnhorse, diversity may be due to within-breed substructure into four sections of the studbook established in 1970: the work horse (draft), trotters, riding horse, and pony [50]. Finally, high diversity in breeds such as the Miniature is likely a result of a diverse founding stock [45], [51], [52]; horses of small size from a variety of geographic regions and bloodlines were utilized in founding the breed, which is defined by phenotype.
All of these three factors, large population size, phenotypic diversity within the breed, and a diversity of founding stock, also lead to the relatively high levels of diversity observed in the Paint and Quarter Horse; in addition, these breeds both allow continued outcrossing between themselves and with the Thoroughbred and have experienced a tremendous population expansion since the formal foundation of the breeds within the past 45–75 years. Due to the relative infancy of these populations, it could be argued that the Paint, Quarter Horse, and other, newly-derived breeds, have not yet had time to undergo the evolutionary processes necessary to be genetically distinct populations as is observed in breeds with longer histories and closed studbooks. However, even with high within-breed diversity and large census population sizes (over 1 million and 4 million worldwide for the Paint and Quarter Horse, respectively), Ne for these breeds account for only a fraction of the census size, demonstrating non-random mating and selection. Outcrossing is also continued in the Swiss Warmblood and Hanoverian breeds, which show similar trends in diversity measures as the Quarter Horse and Paint. The relatively low Ne in these breeds, accompanied by moderate He may partially be due to significant crossing with the Thoroughbred, which would contribute long blocks of LD [40], [41], resulting in decreased estimates of Ne.
Of note in breeds such as the Quarter Horse, Lusitano, and Andalusian, is that despite moderate to high relative levels of He, and low to moderate estimates of f, FIS values in each breed are significantly positive. Significant FIS was also previously observed in the Iberian breeds using microsatellite markers [53]. While selection and inbreeding may be responsible for significant values of FIS in some of these breeds, another instance in which FIS may be significantly positive is in the presence of subpopulation structure within the sample. Evidence of this in the Lusitano and Andalusian is present in parsimony analysis where individuals of the two breeds fall into one clade, but within that clade are two highly supported branches represented by a subset of each breed. In addition, when forcing high values of K in Structure, such as observed at K = 35, Andalusian and Lusitano individuals fall into one of two clusters with q-value >0.5 in a nonbreed-specific manner (data not shown). These results support [54], which showed potential subpopulation structure in the Lusitano via microsatellite analysis. In the Quarter Horse, subpopulation structure is evident through the evaluation of bloodlines and the selection of popular sires for diverse performance classes. This population substructure in the Quarter Horse has also been demonstrated by marked differences in allele frequencies among performance types (cutting, western pleasure, halter, racing, etc.) [55]. A similar instance is found in the US population of the Standardbred, which also has significant excess homozygosity (FIS). Unlike Standardbreds in Europe, which are raced at a trot, those in the US are divergently selected for racing at either the pace or the trot, creating structure within the breed [56].
Finally, several rare populations are included in this dataset. The Caspian is one of the oldest breeds in the Middle East and was thought to be extinct until its recent rediscovery in 1965. The Florida Cracker, a now rare breed, was developed in the United States from feral stock of Iberian descent [57]. The sample size of the Florida Cracker limits the conclusions that can be drawn regarding within-breed diversity. However, the Caspian shows high Ne, He, and estimates of f, given its rarity. After rediscovery of the breed, which historically was believed to represent a type of landrace population, [58] describes a three-year survey, which found approximately 50 individuals remaining, noting that many could not be considered “pure.” In addition, [33] were unable to show evidence of a recent bottleneck in the Caspian breed. The diversity observed in what are now considered Caspian horses likely stems from high levels of diversity within those individuals that founded the modern population.
Among-breed Diversity
The expectation of homogeneity within breeds due to closed populations and selection is supported by the results of AMOVA, which show significant variation present among populations, but a non-significant proportion of variance within. However, the variation among samples lends information about current and historic relationships. Observed trends include patterning based upon geographic origin and/or phenotypic similarities, and relatively low divergence observed in comparisons that include breeds with high within-breed diversity. In Structure analysis, K = 29 was chosen as the most likely value of K; however no single value stood out as the “best” number of clusters to describe these data. Regardless, patterns observed in clustering were also supported by pairwise FST values, parsimony, PCA, and NJ dendograms.
High Diversity and Low Divergence – Landrace Breeds
The Mongolian and Tuva populations are believed to have been influential in the spread of horses across Asia and Europe [45], [59]; these landrace populations, harboring high levels of within-breed diversity, were found to be similar to one another, with a pairwise FST value of 0.006. In addition, with the exception of six Tuva individuals that fell into a clade with the Caspians, both parsimony and NJ analyses place the Tuva and Mongolian into the same clade of each tree. Examining all comparisons, low FST values were observed between the Mongolian and Tuva compared to the other breeds in this study, supporting the potential role of Eurasian horses of similar type in founding modern stocks. This also aligns with high microsatellite diversity observed in Eastern Eurasian “non-breed” (landrace) populations in [19]. On the other hand, breeds with high diversity in general show lower levels of divergence as measured by FST, while those with low diversity show higher values of FST. Low divergence in breeds with high diversity is expected as variation within a breed may indicate outcrossing with other populations, and high variation also makes these breeds more likely to share variation with others by chance. In contrast, if a breed has little within-breed variation, it is less likely to share genetic variation with another breed by chance, especially with another breed that is relatively homogeneous itself. As demonstrated in human literature, source populations are expected to contain greater diversity than those populations which they found [60], [61]; this is also suggested in the horse by [11], which showed greater mtDNA diversity in Iberian breeds than the recently founded American breeds. If the argument is made that the low FST values of Tuvas and Mongolians supports their role in founding modern breeds, the same argument could be made for the Quarter Horse, which also shows low levels of pairwise divergence; however that argument would be unreasonable as the Quarter Horse was developed in only the past century. Therefore, the relative values of FST are informative, but these FST values and data, which represent modern breeds generally derived from limited founding stock, and subjected to intense artificial selective pressures, cannot be used independently to elucidate the evolution of the modern horse. Also, as is the case for other analyses, while the relationships observed can shed light on the history of breeds, they cannot distinguish between recent admixture and shared ancestry.
Thoroughbred-influenced Breeds
The Thoroughbred is believed to have founding sires of Arabian, Turk, and Barb ancestry [42], and [28] found that two sires, noted as being Arabian (Godolphin Arabian and Darley Arabian) together contributed to over 20% of the modern population. However, it is likely that the “Arabian” foundation stallions were not Arabians as the breed is known today. It is noted in [62] that the Godolphin Arabian was a Turkoman stallion with partial Arabian blood, while in other work it is suggested he was a Barb [45]. Regardless of the true ancestry of these stallions, restrictions placed upon the export of purebred Arabians during the 16th and 17th centuries, as well as the general use of the term “Arab” for horses of Middle Eastern descent, it is likely other “Arabian” horses with influence on the Thoroughbred breed also had Turkoman, and Barb bloodlines [45], [62]. The pairwise FST values between the Thoroughbred and the Arabian do not suggest any less divergence than observed between the Thoroughbred and a majority of the other breeds. In addition, at K = 29 the Arabian assigns to the Thoroughbred cluster at only 2.3%. If the Arabian did have significant influence on the Thoroughbred breed, there are several possible explanations for why the supposed Arabian influence is not more apparent. The first is related to SNP ascertainment and the bias of SNPs toward modern variation in Thoroughbred. It is possible that the genes derived from the Arabian are at or near fixation in the Thoroughbred, which would reduce the chance that these SNPs, and the variation described within them are present in the dataset. Another possibility is that the current Arabian sample, taken from the United States, may not reflect the Arabian lineage(s) that were influential in the founding of the Thoroughbred. Finally, as noted above and also suggested elsewhere [63], it may simply be that Arabian bloodlines were not as instrumental in the Thoroughbred breed as once thought or that the initial Arabian influence (and genes) have been selected against or lost to drift during the development of the modern Thoroughbred racehorse.
Within the Thoroughbred itself, divergence between the US and European samples had a significant FST of 0.004, similar to that observed between the Hanoverian and Swiss Warmblood (0.008) and Mongolian and Tuva (0.006), but larger than the minimally observed value seen between the Paint and the Quarter Horse (0.002). Although artificial insemination is prohibited in the Thoroughbred and would be anticipated to limit gene flow to some extent, the founder effect in the original European Thoroughbred by few high-impact sires and dams, accompanied by shared selective pressures, relatively recent importation of the breed to the United States, and ongoing shipment of horses between continents are likely contributing to the lack of geographic population structure identified by parsimony and cluster analyses.
While within-breed diversity of the Thoroughbred was relatively low, and a notable relationship with the Arabian was not observed, among-breed analysis shows a clear influence of the Thoroughbred on many other breeds. Placed with the Thoroughbred in parsimony analysis are the Hanoverian, Maremmano, and Swiss Warmblood. The Maremmano, an Italian breed, shows a q-value of assignment to the Thoroughbred clade of 0.26 at K = 29. This is not surprising given reports that Thoroughbreds contributed over 13% of the maternal lineages to the stallion lines within the stud book [64]. Low differentiation of the Maremmano compared to the Hanoverian was also reported in [29], which is logical given similar influence of the Thoroughbred on the Hanoverian. This and the allowed crossing of Thoroughbreds into the Swiss Warmblood population is reflected in minimal measures of divergence between these samples. The continued influence of the Thoroughbred on the Paint and Quarter Horse is also reflected both by low FST values as well as greater than 30% assignment of the Paint and Quarter Horse to the Thoroughbred cluster. As each of these breeds experience continued gene flow from the Thoroughbred, and had Thoroughbred founding stock, these results are not unexpected. While outcrossing of these breeds is allowed (with restrictions), it is likely that even in cases of breeds with closed stud books, some outcrossing, intentional, unintentional, and/or undocumented, has occurred; it has been demonstrated that historical pedigrees, while helpful, are not always accurate [65]. Issues of outcrossing and individual identification can now be more easily addressed using genetic testing and have the potential to assist managers in decisions regarding breeding and registration.
Finally, the Standardbred samples, which also represent two continents, had a significant pairwise FST value of 0.020, five-fold greater than that observed between the two Thoroughbred samples. This comparison may reflect geographic structure, the influence of French Trotter bloodlines in the European sample, and selection for pacing horses in the US population. Within the horses included in this study, eight were noted to be pacers. These pacers all fall within one clade of the parsimony tree, supported by a bootstrap value of 98%. The limited sample size does not allow a thorough comparison of the pacing vs. trotting Standardbreds, however significant genetic differentiation between horses of the two racing groups has been reported [56].
Middle Eastern and Iberian Breeds
The Middle Eastern breeds, Arabian, Akhal Teke, and Caspian, were placed into a single clade of the NJ tree, with the Arabian and Akhal Teke in their own, highly supported clade. This relationship was supported by low values of K in Structure, which placed the Iberian and Middle Eastern breeds into the same cluster until K = 12. However, parsimony analysis did not support this relationship between the Middle Eastern and Iberian breeds as the Arabian and Akhal Teke individuals were placed into one clade, apart from the Iberian samples and from the Caspians.
In Europe, the Iberian breeds (Lusitano and Andalusian) have only recently been distinguished from one another depending upon the region in which they are bred and divergent selective pressures. In the US, horses of each breed are occasionally interbred [51]. The close relationship between the Andalusian and Lusitano samples in this study is reflected in the minimum FST value observed in either breed, 0.021. The parsimony tree and PCA also shows that individuals cannot necessarily be distinguished from one another regardless of whether they were sampled in the US or Portugal. Of the two clades that appear to suggest population substructure, one includes only Portuguese Lusitanos while the other includes only US samples, although two are of Portuguese ancestry.
Iberian and Gaited Breeds
A horse that is considered to be “gaited” naturally moves in a means other than the traditional walk, trot, canter, and gallop. Alternative gaits in horses are distinguished from traditional gaits by their unique footfall pattern and/or rhythm. The genetic basis of gait has recently been investigated and suggests that all modern gaited breeds share a common ancestor as supported by a shared, extended haplotype spanning a variant significantly associated with the ability to pace [66]; [67]. There is a great deal of historical evidence that the shared ancestry of gaited breeds traces back to Iberian bloodlines, in particular to the Spanish Jenette [45], [51]. Influence of the Iberian breeds on modern gaited breeds is seen in early clustering in Structure analysis as well as in the NJ tree where the Puerto Rican Paso Fino and Peruvian Paso are placed on the same branch as the Andalusian and Lusitanos. In addition to Iberian lines, the Narragansett Pacer is often named as instrumental in the founding of American breeds that may gait including the Saddlebred, Standardbred, and Tennessee Walking Horse [45]. Within those breeds, the Tennessee Walking Horse was documented to be greatly influenced by the Saddlebred, Standardbred, and Morgan [68], [69]; and the Standardbred itself had influence from the Thoroughbred and Morgan (among others) [70]. While we do not have samples of the now extinct Narragansett Pacer, our data set does support a close relationship between the Tennessee Walking Horse and Saddlebred, as observed in the NJ and parsimony analysis as well as with the Morgan and Standardbred. At low values of K in cluster analysis, the Tennessee Walking Horse and Saddlebred cluster strongly. In NJ analysis, the Florida Cracker, which has many individuals that demonstrate the ability to gait, is found intermediate to the Iberian and the modern US gaited breeds. Interestingly, the Icelandic, a four- or five-gaited breed, does not show any significant affinity to the other gaited breeds although they share the recently identified major locus that appears to be essential to the ability to perform alternate gaits [66]; [67]. It thus seems that the genetic variant associated with the gait phenotype arose well before the separation of breeds. Instead of clustering with the other gaited breeds, the Icelandic clusters with the Shetland through K = 16 and also is within a highly supported branch of the NJ tree with the Shetland and Miniature. Finally, the influence of the Shetland on the development of the Miniature is observed at all values of K as well as in the parsimony tree where both breeds occupy the same clade.
Drafts
The Shire and Clydesdale populations share assignment to the same cluster throughout Structure runs. The similarity between these breeds is also seen in a lack of monophyly in the parsimony tree, the sharing of a branch of the NJ tree, positioning in PCA, and a pairwise FST value of 0.037. The Fell Pony, a British breed, falls out as sister taxa to the British draft horses, the Shire and Clydesdales. However FST values show that divergence between the Fell Pony and either the Clydesdale or Shire is not significantly less than seen with most other populations. The other branch of the “draft” clade of the NJ tree contains the breeds from the European mainland, the Belgian, Percheron, and Franches-Montagnes; each of these breeds shows monophyly in parsimony analysis. In addition, similarities among draft and light draft breeds are reflected in cluster analyses at K = 6, which show the populations from the European mainland and Scandinavian Peninsula (Belgian, Finnhorse, Franches-Montagnes, North Swedish Horse, Norwegian Fjord, and Percheron) assign to one cluster with q >0.5. The grouping of the Scandinavian breeds is similar to that previously reported [31]. Geographic relationships are also suggested by the two British breeds, the New Forest Pony and Exmoor that fall just basal to the draft clade in the NJ tree.
Summary
This data set resulting from a large international collaboration represents the first study in the horse to provide an extensive overview of nuclear genetic diversity within, and relationships among a diverse sample of breeds and landrace populations. These data are now available for use in subsequent studies of population-level relationships and provide a baseline for monitoring changes in breed diversity. With high mtDNA diversity but limited paternal input during domestication, this increased understanding of nuclear diversity within the horse will allow for the identification of genomic regions of importance to breed derivation and will be instrumental in guiding across-breed gene discovery projects.
Methods
Ethics Statement
DNA sampling was limited to the collection of blood by jugular venipuncture performed by a licensed veterinarian or from hairs pulled from the mane or tail by the horse owner or researcher. All animal work was conducted in accordance with and approval from the international and national governing bodies at the institutions in which samples were collected (the University of Minnesota Institutional Animal Care and Use Committee (IACUC); the University of Kentucky IACUC; the University College Dublin, Animal Research Ethics Committee; Swiss Law on Animal Protection and Welfare; the Ethical Board of the University of Helsinki; the Animal Health Trust Clinical Research Ethics Committee; Norwegian Animal Research Authority; UK Home Office License; and the Lower Saxon state veterinary office).
Samples and Genotyping
Tissue samples and previously collected genotypes from 1,060 horses were obtained from members of the EGDC or were obtained by our laboratory. 814 samples representing the 38 populations included in this study were selected from the EDGC sample collection with the goal of obtaining as random of a sample as possible and to minimize close relationships among individuals. In some cases, genotypes were available from breeds collected for genome-wide association studies (GWAS). In all cases, when pedigree information was available, no relationships were allowed at or more recent to the grandsire/dam level. If no pedigrees were available, once genotyping was performed, individuals were removed from the analyses to reduce genome sharing as measured by autosomal estimates of identity by descent (pi hat) values in PLINK [71] greater than 0.3 (after pruning for MAF>0.05). In samples that were obtained as a result of GWAS, “control” individuals were preferentially chosen for inclusion in these analyses. When necessary, DNA isolation from hair roots took place using a modification of the Puregene (Qiagen) protocol for DNA purification from tissue. Modifications include the addition of 750 µl of isopropanol rather than 300, increasing the precipitation spin time to 15 m at 4°C, and washing the pellet twice. Approximately 1 µg of DNA was used for SNP genotyping using the Illumina SNP50 Beadchip according to the manufacturer’s protocol. All genotype calls were extracted from the raw intensity data using GenomeStudio (Illumina) with the minimum gencall score threshold of 0.15. The raw intensity scores were available for all populations with the exception of the Lusitano and Maremmano.
Data Pruning
SNP discovery was conducted using horses from seven breeds (Akhal Teke, Andalusian, Arabian, Icelandic, Quarter Horse, Standardbred, Thoroughbred) as well as the reference genome of a Thoroughbred mare [41]. To eliminate ascertainment bias as much as possible, horses from the discovery breeds were removed from the dataset, which was then pruned to exclude SNPs with MAF less than 0.05. All horses were then replaced and those SNPs removed from all analyses. In this new, complete data set, SNP markers that failed to genotype in at least 99% of the individuals and SNPs that had a MAF of 0.05 or less across all samples were removed as well as SNPs on ECAX. SNPs that were in LD across breeds were also removed; files used for basic diversity indices were pruned for r2<0.1 in PLINK [71] considering 100 SNP windows and moving 25 SNPs per set (–indep-pairwise 100 25 1.11). Allowing for additional LD, data sets were also created for r2<0.2 and 0.4. An additional data set, used for Structure analysis was pruned for R2<0.1 in Plink (–indep; R = multiple correlation coefficient), which is similar to the above method but instead of analyzing pairwise relationships of SNPs as in the former method, uses a multiple regression approach upon the SNPs in the analysis window. Files were converted for usage between analyses programs using PLINK, perl script, CONVERT [72] and/or PGDSpider 2.0.1.4 [73].
Within-breed Diversity
Expected heterozygosity (He) and AMOVA were calculated in Arlequin3.5 [74] on all four data sets. AMOVA was conducted on the primary data set with breeds designated as populations and excluding the Florida Cracker due to small sample size. Analyses were also conducted grouping the two Thoroughbred and Standardbred samples together by breed. FIS was calculated and significance tested on the primary data set, with 10,000 permutations of the data in Genetix [75]. Individual inbreeding coefficients (f) were calculated in PLINK based upon loss of heterozygosity (–het).
Among Breed Relationships
Pairwise FST values were calculated on the primary, 10,536 SNP dataset in Arlequin3.5 [74] using Reynolds’ distance [76] with significance tested using 20,000 permutations.
A neighbor joining (NJ) cladogram was built using breed allele frequencies calculated from the primary SNP set using the packages seqboot, gendist, neighbor, and consense and Nei’s genetic distance [38] in PHYLIP ver3.69 [77]. Bootstrap support from 1,000 iterations of the data was used to assess support for the resulting majority rule consensus cladogram.
A parsimony cladogram was constructed using 10,066 SNP markers pruned from the original data set using the MAF and genotyping rate criteria as above and allowing for R2<0.2 (–indep-pairwise). The domestic ass was included as an outgroup for traditional and new-technology searches in TNT [78].
Principal component analysis (PCA) was conducted in snpStats in R (http://cran.r-project.org) on the full SNP set consisting of all 814 individuals and 38,755 autosomal SNPs (pruned only for MAF and genotyping rate).
Cluster Analysis
Clustering of breeds into genetic groups was examined using the program Structure 2.3.3 [79], [80] assuming K = 1 to 45. The Structure algorithm included the admixture model and correlated allele frequencies. Three iterations of each K value were conducted with 35,000 MCMC repetitions (15,000 burn-in). The convergence of Structure runs was evaluated by equilibrium of alpha and likelihood scores. The value of K most suitable to explain the diversity in these data was predicted by the highest mean estimated ln P(X|K) while minimizing variance and also making biological sense [80], [81]. The replicates from each run of K were input into CLUMPP [82] and the average cluster membership calculated using the LargeK Greedy algorithm. Output from CLUMPP was visualized in Distruct [83].
Effective Population Size
To estimate effective population size (Ne), the full set of 54,602 SNP markers was pruned within each population to remove those with MAF <0.01 and genotyping rate of <0.05. Pairwise r2 values between remaining SNPs were calculated in Haploview [84], for each population considering intermarker distances from 0 to 4 Mb in 50 kb increments. Values of Ne were calculated using the method of [37], which includes a correction for small sample size and the assumption that 1 Mb = 1 cM.
Data Access
All SNP genotype data are available at the NAGPR Community Data Repository (animalgenome.org) for the purpose of reconstructing the analyses. The only exception is the data collected from the Tennessee Walking Horse, which, under agreement from the granting agency (to the University of Minnesota from the Foundation for the Advancement of the Tennessee Walking Show Horse (FAST) and the Tennessee Walking Horse Foundation (TWHF)), is only available under a Material Transfer Agreement (MTA) between interested individuals and the University of Minnesota.
Acknowledgments
The computational support provided by Rob Schaefer and Aaron Rendahl is gratefully acknowledged. Finally, the EGDC thanks the Havemeyer Foundation for their continued support of equine genomics research.
Funding Statement
This work was supported by National Research Initiative Competitive Grants 2008-35205-18766, 2009-55205-05254, and 2012-67015-19432 from the United States Department of Agriculture-National Institute of Food and Agriculture (USDA-NIFA); Foundation for the Advancement of the Tennessee Walking Show Horse and Tennessee Walking Horse Foundation; National Institutes of Health-National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIH-NAIMS) grant 1K08AR055713-01A2 (MEM salary support), and 2T32AR007612 (JLP salary support), American Quarter Horse Foundation grant “Selective Breeding Practices in the American Quarter Horse: Impact on Health and Disease 2011–2012”; Morris Animal Foundation Grant D07EQ-500; The Swedish Research Council FORMAS (Contract 221-2009-1631 and 2008-617); The Swedish-Norwegian Foundation for Equine Research (Contract H0847211 and H0947256); The Carl Tryggers Stiftelse (Contract CTS 08∶29); IBB-CBA-UAç was supported by FCT and DRCT and MSL by FRCT/2011/317/005; Science Foundation Ireland Award (04/Y11/B539) to EWH; Volkswagen Stiftung und Niedersächsisches Ministerium für Wissenschaft und Kultur, Germany (VWZN2012). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
1. FAOSTAT (2010). Available: http://www.faostat.fao.org.Accessed 2012 Jun 5.
2. Lippold S, Matzke NJ, Reissmann M, Hofreiter M (2011) Whole mitochondrial genome sequencing of domestic horses reveals incorporation of extensive wild horse diversity during domestication. BMC Evol Biol 11: 328. [PMC free article] [PubMed]
3. Ludwig A, Pruvost M, Reissmann M, Benecke N, Brockmann GA, et al. (2009) Coat color variation at the beginning of horse domestication. Science 324: 485. [PMC free article] [PubMed]
4. Outram AK, Stear NA, Bendrey R, Olsen S, Kasparov A, et al. (2009) The earliest horse harnessing and milking. Science 323: 1332–1335. [PubMed]
5. Pedrosa S, Uzun M, Arranz JJ, Gutierrez-Gill B, Primitivo FS, et al. (2005) Evidence of three maternal lineages in near eastern sheep supporting multiple domestication events. P R Soc B 272: 2211–2217. [PMC free article] [PubMed]
6. Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, et al. (2005) Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science 307: 1618–1621. [PubMed]
7. Wu GS, Yao YG, Qu KX, Ding ZL, Li H, et al. (2007) Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia. Genome Biol 8: R245. [PMC free article] [PubMed]
8. Cieslak M, Pruvost M, Benecke N, Hofreiter M, Morales A, et al. (2010) Origin and history of mitochondrial DNA lineages in domestic horses. PLOS One 5: e15311. [PMC free article] [PubMed]
9. Jansen T, Forster P, Levine MA, Oelke H, Hurles M, et al. (2002) Mitochondrial DNA and the origins of the domestic horse. PNAS 99: 10905–10910. [PMC free article] [PubMed]
10. Lei CZ, Su R, Bower MA, Edwards CJ, Wang XB, et al. (2009) Multiple maternal origins of native modern and ancient horse populations in China. Anim Genet 40: 933–944. [PubMed]
11. Lira J, Linderholm A, Olaria C, Brandstrom Durling M, Gilbert MT, et al. (2010) Ancient DNA reveals traces of Iberian Neolithic and Bronze Age lineages in modern Iberian horses. Mol Ecol 19: 64–78. [PubMed]
12. Vila C, Leonard JA, Gotherstrom A, Marklund S, Sandberg K, et al. (2001) Widespread origins of domestic horse lineages. Science 291: 474–477. [PubMed]
13. Warmuth V, Eriksson A, Bower MA, Barker G, Barrett E, et al. (2012) Reconstructing the origin and spread of horse domestication in the Eurasian steppe. PNAS 109: 8202–8206. [PMC free article] [PubMed]
14. Beja-Pereira A, Caramelli D, Lalueza-Fox C, Vernesi C, Ferrand N, et al. (2006) The origin of European cattle: evidence from modern and ancient DNA. PNAS 103: 8113–8118. [PMC free article] [PubMed]
15. Gotherstrom A, Anderung C, Hellborg L, Elburg R, Smith C, et al. (2005) Cattle domestication in the Near East was followed by hybridization with aurochs bulls in Europe. P R Soc B 272: 2345–2350. [PMC free article] [PubMed]
16. Keyser-Tracqui C, Blandin-Frappin P, Francfort HP, Ricaut FX, Lepetz S, et al. (2005) Mitochondrial DNA analysis of horses recovered from a frozen tomb (Berel site, Kazakhstan, 3rd Century BC). Anim Genet 36: 203–209. [PubMed]
17. Achilli A, Olivieri A, Soares P, Lancioni H, Kashani BH, et al. (2012) Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication. PNAS 109: 2449–2454. [PMC free article] [PubMed]
18. McGahern A, Bower MA, Edwards CJ, Brophy PO, Sulimova G, et al. (2006) Evidence for biogeographic patterning of mitochondrial DNA sequences in Eastern horse populations. Anim Genet 37: 494–497. [PubMed]
19. Warmuth V, Manica A, Eriksson A, Barker G, Bower M (2012) Autosomal genetic diversity in non-breed horses from eastern Eurasia provides insights into historical population movements. Anim Genet. In press. doi:10.1111/j.1365–2052.2012.02371.x. [PubMed]
20. Lindgren G, Backstrom N, Swinburne J, Hellborg L, Einarsson A, et al. (2004) Limited number of patrilines in horse domestication. Nat Genet 36: 335–336. [PubMed]
21. Ling Y, Ma Y, Guan W, Cheng Y, Wang Y, et al. (2010) Identification of y chromosome genetic variations in Chinese indigenous horse breeds. J Hered 101: 639–643. [PubMed]
22. Aberle KS, Hamann H, Drogemuller C, Distl O (2004) Genetic diversity in German draught horse breeds compared with a group of primitive, riding and wild horses by means of microsatellite DNA markers. Anim Genet 35: 270–277. [PubMed]
23. Achmann R, Curik I, Dovc P, Kavar T, Bodo I, et al. (2004) Microsatellite diversity, population subdivision and gene flow in the Lipizzan horse. Anim Genet 35: 285–292. [PubMed]
24. Bjornstad G, Gunby E, Roed KH (2000) Genetic structure of Norwegian horse breeds. J Anim Breed Genet 117: 307–317.
25. Bömcke E, Gengler N, Cothran EG (2010) Genetic variability in the Skyros pony and its relationship with other Greek and foreign horse breeds. Genet and Mol Biol 34: 68–76. [PMC free article] [PubMed]
26. Canon J, Checa ML, Carleos C, Vega-Pla JL, Vallejo M, et al. (2000) The genetic structure of Spanish Celtic horse breeds inferred from microsatellite data. Anim Genet 31: 39–48. [PubMed]
27. Cothran EG, Canelon JL, Luis C, Conant E, Juras R (2011) Genetic analysis of the Venezuelan Criollo horse. Genet Mol Res 10: 2394–2403. [PubMed]
28. Cunningham EP, Dooley JJ, Splan RK, Bradley DG (2001) Microsatellite diversity, pedigree relatedness and the contributions of founder lineages to thoroughbred horses. Anim Genet 32: 360. [PubMed]
29. Felicetti M, Lopes MS, Verini-Supplizi A, Machado Ada C, Silvestrelli M, et al. (2010) Genetic diversity in the Maremmano horse and its relationship with other European horse breeds. Anim Genet 41 Suppl 2 53–55. [PubMed]
30. Glowatzki-Mullis ML, Muntwyler J, Pfister W, Marti E, Rieder S, et al. (2005) Genetic diversity among horse populations with a special focus on the Franches-Montagnes breed. Anim Genet 37: 33–39. [PubMed]
31. Leroy G, Callede L, Verrier E, Meriaux JC, Ricard A, et al. (2009) Genetic diversity of a large set of horse breeds raised in France assessed by microsatellite polymorphism. Genet Sel Evol 41: 5. [PMC free article] [PubMed]
32. Marletta D, Tupac-Yupanqui I, Bordonaro S, Garcia D, Guastella AM, et al. (2006) Analysis of genetic diversity and the determination of relationships among western Mediterranean horse breeds using microsatellite markers. J Anim Breed Genet 123: 315–325. [PubMed]
33. Shasavarani H, Rahimi-Mianji G (2010) Analysis of genetic diversity and estimation of inbreeding coefficient within Caspian horse population using microsatellite markers. Afr J Biotechnol 9: 293–299.
34. Solis A, Jugo BM, Meriaux J-C, Iriondo M, Mazon LI, et al. (2005) Genetic diversity within and among four South European native horse breeds based on microsatellite DNA analysis: Implications for conservation. J Hered 96: 670–678. [PubMed]
35. Tozaki T, Takezaki N, Hasegawa T, Ishida N, Kurosawa M, et al. (2003) Microsatellite variation in Japanese and Asian horses and their phylogenetic relationships using a European horse outgroup. J Hered 94: 374–380. [PubMed]
36. Luis C, Juras R, Oom MM, Cothran EG (2007) Genetic diversity and relationships of Protuguese and other horse breeds based on protein and microsatellite loci variation. Anim Genet 38: 20–27. [PubMed]
37. Weir BS, Hill WG (1980) Effect of mating structure on variation in linkage disequilibrium. Genetics 95: 477–488. [PMC free article] [PubMed]
38. Nei M (1972) Genetic distance between populations. Amer Nat 106: 283–292.
39. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620. [PubMed]
40. Wade CM, Giulotto E, Sigurdsson S, Zoli M, Gnerre S, et al. (2009) Genome sequence, comparative analysis, and population genetics of the domestic horse. Science 326: 865–867. [PMC free article] [PubMed]
41. McCue ME, Bannasch DL, Petersen JL, Gurr J, Bailey E, et al. (2012) A high density SNP array for the domestic horse and extant perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies. PLOS Genet 8: e1002451. [PMC free article] [PubMed]
42. Weatherby and Sons (1791) An Introduction to a General Stud Book. London, UK.
43. Corbin LJ, Blott SC, Swinburne JE, Vaudin M, Bishop SC, et al. (2010) Linkage disequilibrium and historical effective population size in the Thoroughbred horse. Anim Genet 41 Suppl 2 8–15. [PubMed]
44. Binns MM, Boehler DA, Bailey E, Lear TL, Cardwell JM, et al. (2012) Inbreeding in the Thoroughbred horse. Anim Genet 43: 340–342. [PubMed]
45. Hendricks BL (2007) International Encyclopedia of Horse Breeds. Norman, OK: University of Oklahoma Pres. 486 p.
46. Weatherley L (1978) Great horses of Britain. Hindhead: Spur Publications. viii, 269 p.
47. Gates S (1979) A study of the home ranges of free-ranging Exmoor ponies. Mammal Rev 9: 3–18.
48. Peilieu C (1984) Livestock Breeds of China. FAO: Animal Production and Health paper 46: 1–217.
49. Hund A (2008) The Stallion’s Mane: The Next Generation of Horses in Mongolia. Available: http://digitalcollections.sit.edu/cgi/viewcontent.cgi?article=1543&context=isp_collection. Accessed 2012 Dec 31.
50. Ticklen M, editor (2006) Get to know the Finnhorse. Available: http://www.suomenhevonen.info/hippos/sh2007/pdf/SHjulkaisu_englanti_nettiin.pdf. Accessed 2012 Dec 31.
51. Bowling AT, Ruvinsky A (2000) The genetics of the horse. Wallingford: CABI Pub. viii, 527 p.
52. Lynghaug F (2009) The Official Horse Breeds Standards Guide: The Complete Guide to the Standards of All North American Equine Breed Associations. Stillwater, MN: Voyageur Press.
53. Conant EK, Juras R, Cothran EG (2012) A microsatellite analysis of five Colonial Spanish horse populations of the southeastern United States. Anim Genet 43: 53–62. [PubMed]
54. Lopes MS (2011) Molecular tools for the characterisation of the Lusitano horse. Angra do Heroísmo: University of Azores. 219 p.
55. Tryon RC, Penedo MC, McCue ME, Valberg SJ, Mickelson JR, et al. (2009) Evaluation of allele frequencies of inherited disease genes in subgroups of American Quarter Horses. J Am Vet Med Assoc 234: 120–125. [PubMed]
56. Cothran EG, MacCluer JW, Weitkamp LR, Bailey E (1987) Genetic differentiation associated with gait within American standardbred horses. Anim Genet 18: 285–296. [PubMed]
57. Florida Cracker Horse Association. Available: http://www.floridacrackerhorses.com/index.htm. Accessed 2012 Jun 5.
58. Firouz L (1969) Conservation of a domestic breed. Biol Conserv 1: 1–2.
59. Bjornstad G, Nilsen NO, Roed KH (2003) Genetic relationship between Mongolian and Norwegian horses? Anim Genet 34: 55–58. [PubMed]
60. Conrad DF, Jakobsson M, Coop G, Wen X, Wall JD, et al. (2006) A worldwide survey of haplotype variation and linkage disequilibrium in the human genome. Nat Genet 38: 1251–1260. [PubMed]
61. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, et al. (2008) Worldwide human relationships inferred from genome-wide patterns of variation. Science 319: 1100–1104. [PubMed]
62. Mackay-Smith A (2000) Speed and the thoroughbred : the complete history. Lanham, MD: Derrydale Press. xxvii, 193 p.
63. Bower MA, Campana MG, Whitten M, Edwards CJ, Jones H, et al. (2011) The cosmopolitan maternal heritage of the Thoroughbred racehorse breed shows a significant contribution from British and Irish native mares. Biol Letters 7: 316–320. [PMC free article] [PubMed]
64. Silvestrelli M (1991) The Maremmano horse. Animal Genetic Resources Information (FAO/UNEP) 8: 74–83.
65. Bower MA, Campana MG, Nisbet RER, Weller R, Whitten M, et al. (2012) Truth in the bones: Resolving the identify of the founding elite Thoroughbred racehorses. Archaeometry 54: 916–925.
66. Andersson LS, Larhammar M, Memic F, Wootz H, Schwochow D, et al. (2012) Mutations in DMRT3 affect locomotion in horses and spinal circuit function in mice. Nature 488: 642–646. [PMC free article] [PubMed]
67. Petersen JL, Mickelson JR, Rendahl AK, Valberg SJ, Andersson LS, et al. . (2012) Genome-wide analysis reveals selection for important traits in domestic horse breeds. PLOS Genet. In press. [PMC free article] [PubMed]
68. Fletcher JL (1946) A study of the first fifty years of Tennessee walking horse breeding. J Hered 37: 369–373. [PubMed]
69. Tennessee Walking Horse Breeders and Exhibitors Association. Available: http://www.twhbea.com/breed/history.php. Accessed 2012 Jun 5.
70. MacCluer JW, Boyce AJ, Dyke B, Weitkamp LR, Pfennig DW, et al. (1983) Inbreeding and pedigree structure in Standardbred horses. J Hered 74: 394–399.
71. Purcell S, Neale B, Tood-Brown K, Thomas L, Ferreira MAR, et al. (2007) PLINK: a toolset for whole-genome association and population-based linkage analysis. Amer J Hum Genet 81: 559–575. [PMC free article] [PubMed]
72. Glaubitz JC (2004) CONVERT: A user-friendly program to reformat diploid genotypic data for commonly used population genetic software packages. Mol Ecol Notes 4: 309–310.
73. Lischer H, Excoffier L (2012) An automated data conversion tool for connecting population genetics and genomics programs. Bioinformatics 28: 298–299. [PubMed]
74. Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform 1: 47–50. [PMC free article] [PubMed]
75. Belkhir KP, Borsa P, Chikhi L, Raufaste N, Bonhomme F (2001) GENTETIX, logiciel sous Windows TM pour la genetique des population. Universite de Montpellier II, Montpellier, France.
76. Reynolds J, Weir BS, Cockerham CC (1983) Estimation of the coancestry coefficient: basis for a short-term genetic distance. Genetics 105: 767–779. [PMC free article] [PubMed]
77. Felsenstein J (1989) PHYLIP-Phylogeny Inference Package (version 3.2) Cladistics. 5: 164–166.
78. Goloboff PA, Farris JS, Nixon K (2003) TNT: tree analysis using new technology. Syst Biol 54: 176–178.
79. Falush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164: 1567–1587. [PMC free article] [PubMed]
80. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959. [PMC free article] [PubMed]
81. Pritchard JK, Wen X, Falush D (2010) Documentation for Structure software: version 2.3.38.
82. Jakobsson M, Rosenberg NA (2007) CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23: 1801–1806. [PubMed]
83. Rosenberg NA (2004) DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes 4: 137–138.
84. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265. [PubMed]