The evolutionary history of the DMRT3 ‘Gait keeper’ haplotype

A previous study revealed a strong association between the DMRT3:Ser301STOP mutation in horses and alternate gaits as well as performance in harness racing. Several follow‐up studies have confirmed a high frequency of the mutation in gaited horse breeds and an effect on gait quality. The aim of this study was to determine when and where the mutation arose, to identify additional potential causal mutations and to determine the coalescence time for contemporary haplotypes carrying the stop mutation. We utilized sequences from 89 horses representing 26 breeds to identify 102 SNPs encompassing the DMRT3 gene that are in strong linkage disequilibrium with the stop mutation. These 102 SNPs were genotyped in an additional 382 horses representing 72 breeds, and we identified 14 unique haplotypes. The results provided conclusive evidence that DMRT3:Ser301STOP is causal, as no other sequence polymorphisms showed an equally strong association to locomotion traits. The low sequence diversity among mutant chromosomes demonstrated that they must have diverged from a common ancestral sequence within the last 10 000 years. Thus, the mutation occurred either just before domestication or more likely some time after domestication and then spread across the world as a result of selection on locomotion traits.     

Comparison of DMRT3 genotypes among American Saddlebred horses with reference to gait

Horse owners choose whether or not to train American Saddlebred horses (ASHs) to perform the 4‐beat gaits called rack and slow gait. The rack and slow gait are similar to ambling gaits shown to be associated with variation in the DMRT3 gene in other breeds but are trained rather than naturally occurring gaits. A premature stop codon in the DMRT3 gene (DMRT3_Ser301STOP) caused by the Ch23:g.22999655C>A SNP has an effect on the pattern of locomotion in horses and allows for the pacing gait and strong association with performance of ambling gaits in diverse breeds. We used horse show records to identify ASHs trained to perform as 5‐gaited horses and compared them to other Saddlebreds for the prevalence of the A allele of the Ch23:g.22999655C>A SNP of DMRT3. Genomic DNA was typed using a PCR‐RFLP technique for 5‐gaited ASHs (n = 55), 3‐gaited or harness ASHs (n = 64) and ASHs not shown (n = 128). We did not observe differences in the frequencies of the Ch23:g.22999655C>A SNP of DMRT3 when comparing 5‐gaited ASHs to other ASHs (P > 0.22). These results suggest that DMRT3 alleles do not affect the choice of breeders to train and show ASHs as 5‐gaited horses.     

To pace or not to pace: a pilot study of four‐ and five‐gaited Icelandic horses homozygous for the DMRT3 ‘Gait Keeper’ mutation

The Icelandic horse is a breed known mainly for its ability to perform the ambling four‐beat gait ‘tölt’ and the lateral two‐beat gait pace. The natural ability of the breed to perform these alternative gaits is highly desired by breeders. Therefore, the discovery that a nonsense mutation (C>A) in the DMRT3 gene was the main genetic factor for horses' ability to perform gaits in addition to walk, trot and canter was of great interest. Although several studies have demonstrated that homozygosity for the DMRT3 mutation is important for the ability to pace, only about 70% of the homozygous mutant (AA) Icelandic horses are reported to pace. The aim of the study was to genetically compare four‐ and five‐gaited (i.e. horses with and without the ability to pace) AA Icelandic horses by performing a genome‐wide association (GWA) analysis. All horses (n = 55) were genotyped on the 670K Axiom Equine Genotyping Array, and a GWA analysis was performed using the genabel package in r . No SNP demonstrated genome‐wide significance, implying that the ability to pace goes beyond the presence of a single gene variant. Despite its limitations, the current study provides additional information regarding the genetic complexity of pacing ability in horses. However, to fully understand the genetic differences between four‐ and five‐gaited AA horses, additional studies with larger sample materials and consistent phenotyping are needed.     

DMRT3 is associated with gait type in Mangalarga Marchador horses,but does not control gait ability

The Mangalarga Marchador (MM) is a Brazilian horse breed known for a uniquely smooth gait. A recent publication described a mutation in the DMRT3 gene that the authors claim controls the ability to perform lateral patterned gaits (Andersson et al. 2012). We tested 81 MM samples for the DMRT3 mutation using extracted DNA from hair bulbs using a novel RFLP. Horses were phenotypically categorized by their gait type (batida or picada), as recorded by the Brazilian Mangalarga Marchador Breeders Association (ABCCMM). Statistical analysis using the plink toolset (Purcell, 2007) revealed significant association between gait type and the DMRT3 mutation (P = 2.3e‐22). Deviation from Hardy–Weinberg equilibrium suggests that selective pressure for gait type is altering allele frequencies in this breed (P = 1.00e‐5). These results indicate that this polymorphism may be useful for genotype‐assisted selection for gait type within this breed. As both batida and picada MM horses can perform lateral gaits, the DMRT3 mutation is not the only locus responsible for the lateral gait pattern.     

Single nucleotide polymorphisms of myostatin gene in Chinese domestic horses

The myostatin gene (MSTN) is a genetic determinant of skeletal muscle growth. Single nucleotide polymorphisms (SNP) in MSTN are of importance due to their strong associations with horse racing performances. In this study, we screened the SNPs in MSTN gene in 514 horses from 15 Chinese horse breeds. Six SNPs (g.26 T > C, g.156 T > C, g.587A > G, g.598C > T, g.1485C > T, g.2115A > G) in MSTN gene were detected by sequencing and genotyped using PCR-RFLP method. The g.587A > G and g.598C > T residing in the 5′UTR region were novel SNPs identified by this study. The g.2115A > G which have previously been associated with racing performances were present in Chinese horse breeds, providing valuable genetic information for evaluating the potential racing performances in Chinese domestic breeds. The six SNPs together defined thirteen haplotypes, demonstrating abundant haplotype diversities in Chinese horses. Most of the haplotypes were shared among different breeds with no haplotype restricted to a specific region or a single horse breed. AMOVA analysis indicated that most of the genetic variance was attributable to differences among individuals without any significant contribution by the four geographical groups. This study will provide fundamental and instrumental genetic information for evaluating the potential racing performances of Chinese horse breeds.     

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.     

Glycogen branching enzyme (<Emphasis Type="Italic">GBE1</Emphasis>) mutation causing equine glycogen storage disease IV

Comparative biochemical and histopathological evidence suggests that a deficiency in the glycogen branching enzyme, encoded by the GBE1 gene, is responsible for a recently identified recessive fatal fetal and neonatal glycogen storage disease (GSD) in American Quarter Horses termed GSD IV. We have now derived the complete GBE1 cDNA sequences for control horses and affected foals, and identified a C to A substitution at base 102 that results in a tyrosine (Y) to stop (X) mutation in codon 34 of exon 1. All 11 affected foals were homozygous for the X34 allele, their 11 available dams and sires were heterozygous, and all 16 control horses were homozygous for the Y34 allele. The previous findings of poorly branched glycogen, abnormal polysaccharide accumulation, lack of measurable GBE1 enzyme activity and immunodetectable GBE1 protein, coupled with the present observation of abundant GBE1 mRNA in affected foals, are all consistent with the nonsense mutation in the 699 amino acid GBE1 protein. The affected foal pedigrees have a common ancestor and contain prolific stallions that are likely carriers of the recessive X34 allele. Defining the molecular basis of equine GSD IV will allow for accurate DNA testing and the ability to prevent occurrence of this devastating disease affecting American Quarter Horses and related breeds.The nucleotide sequence data reported in this article have been submitted to GenBank and have been assigned the accession numbers AY505107–AY505110.     

Whole-Genome SNP Association in the Horse: Identification of a Deletion in Myosin Va Responsible for Lavender Foal Syndrome

Lavender Foal Syndrome (LFS) is a lethal inherited disease of horses with a suspected autosomal recessive mode of inheritance. LFS has been primarily diagnosed in a subgroup of the Arabian breed, the Egyptian Arabian horse. The condition is characterized by multiple neurological abnormalities and a dilute coat color. Candidate genes based on comparative phenotypes in mice and humans include the ras-associated protein RAB27a (RAB27A) and myosin Va (MYO5A). Here we report mapping of the locus responsible for LFS using a small set of 36 horses segregating for LFS. These horses were genotyped using a newly available single nucleotide polymorphism (SNP) chip containing 56,402 discriminatory elements. The whole genome scan identified an associated region containing these two functional candidate genes. Exon sequencing of the MYO5A gene from an affected foal revealed a single base deletion in exon 30 that changes the reading frame and introduces a premature stop codon. A PCR–based Restriction Fragment Length Polymorphism (PCR–RFLP) assay was designed and used to investigate the frequency of the mutant gene. All affected horses tested were homozygous for this mutation. Heterozygous carriers were detected in high frequency in families segregating for this trait, and the frequency of carriers in unrelated Egyptian Arabians was 10.3%. The mapping and discovery of the LFS mutation represents the first successful use of whole-genome SNP scanning in the horse for any trait. The RFLP assay can be used to assist breeders in avoiding carrier-to-carrier matings and thus in preventing the birth of affected foals.     

Characteristics of Adopted Thoroughbred Racehorses in Second Careers

The unwanted horse issue continues to be a major concern in the U.S. equine industry. Nonprofit organizations dedicated to rescuing, retraining, and rehoming unwanted horses are critical in minimizing this problem. This study utilized data collected nationwide from organizations that provide these services for thoroughbreds retired from racing to identify individual horse characteristics that influenced length of stay at the adoption facility as well as characteristics that increased the probability that an adopted horse would be returned to the facility. The results suggested that horses with fewer activity limitations were rehomed more quickly (p < .01), as were gray horses (relative to bays, p < .03) and stallions (relative to geldings, p < .04). Older horses took longer to rehome (p < .05). Interestingly, the results also suggested that gray horses were more likely to be returned to the facility postadoption (p < .02). Results from this study could benefit thoroughbreds retired from racing, nonprofit organizations, end consumers, and the thoroughbred racing industry.     

Epidemiology of Exertional Rhabdomyolysis Susceptibility in Standardbred Horses Reveals Associated Risk Factors and Underlying Enhanced Performance

Background

Exertional rhabdomyolysis syndrome is recognised in many athletic horse breeds and in recent years specific forms of the syndrome have been identified. However, although Standardbred horses are used worldwide for racing, there is a paucity of information about the epidemiological and performance-related aspects of the syndrome in this breed. The objectives of this study therefore were to determine the incidence, risk factors and performance effects of exertional rhabdomyolysis syndrome in Standardbred trotters and to compare the epidemiology and genetics of the syndrome with that in other breeds.

Methodology/Principal Findings

A questionnaire-based case-control study (with analysis of online race records) was conducted following identification of horses that were determined susceptible to exertional rhabdomyolysis (based on serum biochemistry) from a total of 683 horses in 22 yards. Thirty six exertional rhabdomyolysis-susceptible horses were subsequently genotyped for the skeletal muscle glycogen synthase (GYS1) mutation responsible for type 1 polysaccharide storage myopathy. A total of 44 susceptible horses was reported, resulting in an annual incidence of 6.4 (95% CI 4.6–8.2%) per 100 horses. Female horses were at significantly greater risk than males (odds ratio 7.1; 95% CI 2.1–23.4; p = 0.001) and nervous horses were at a greater risk than horses with calm or average temperaments (odds ratio 7.9; 95% CI 2.3–27.0; p = 0.001). Rhabdomyolysis-susceptible cases performed better from standstill starts (p = 0.04) than controls and had a higher percentage of wins (p = 0.006). All exertional rhabdomyolysis-susceptible horses tested were negative for the R309H GYS1 mutation.

Conclusions/Significance

Exertional rhabdomyolysis syndrome in Standardbred horses has a similar incidence and risk factors to the syndrome in Thoroughbred horses. If the disorder has a genetic basis in Standardbreds, improved performance in susceptible animals may be responsible for maintenance of the disorder in the population.     

Genome‐wide association mapping and examination of possible maternal effect for the pace trait of horses

Previously, a single nucleotide polymorphism (SNP) related to gait type was identified at position 22 999 655 of chromosome 23 in the coding region of DMRT3 (DMRT3:Ser301Ter) by showing that a cytosine (C) to adenine (A) mutation of this SNP induced pace in the Icelandic horse. We investigated the effect of DMRT3:Ser301Ter on the gait of Hokkaido Native Horses, a Japanese native breed, and examined genetic factors other than DMRT3 by exploring genome‐wide SNPs related to gait determination. All animals exhibiting pace were AA for DMRT3:Ser301Ter, confirming the association of DMRT3:Ser301Ter with gait determination; however, 14.3% of the animals exhibiting trot also had AA for DMRT3:Ser301Ter, suggesting the presence of another factor(s) cooperatively working with DMRT3:Ser301Ter for gait determination. SNPs on chromosomes 13 and 23 were detected by genome‐wide association analysis (false discovery rate <0.05), although SNPs on chromosome 23 were all located in the vicinity of DMRT3:Ser301Ter, confirming the association with DMRT3. A genome‐wide association study targeting only animals with AA for DMRT3:Ser301Ter to examine genetic factors cooperatively working with DMRT3:Ser301Ter for gait determination suggested associations of 23 SNPs on six chromosomes. In a series of analyses of the effect of a maternal factor (dam's gait) on gait determination, the effect was suggested in comparison of the frequencies of exhibiting pace in gait checks in only two animal groups having dams with different DMRT3:Ser301Ter genotypes (P < 0.05), suggesting that the gait of the dam does not have a major effect on whether progeny homozygous for the DMRT3:Ser301Ter mutation will preferentially pace or trot.     

Haplotype diversity in the equine myostatin gene with focus on variants associated with race distance propensity and muscle fiber type proportions

Two variants in the equine myostatin gene (MSTN), including a T/C SNP in the first intron and a 227‐bp SINE insertion in the promoter, are associated with muscle fiber type proportions in the Quarter Horse (QH) and with the prediction of race distance propensity in the Thoroughbred (TB). Genotypes from these loci, along with 18 additional variants surrounding MSTN, were examined in 301 horses of 14 breeds to evaluate haplotype relationships and diversity. The C allele of intron 1 was found in 12 of 14 breeds at a frequency of 0.27; the SINE was observed in five breeds, but common in only the TB and QH (0.73 and 0.48 respectively). Haplotype data suggest the SINE insertion is contemporary to and arose upon a haplotype containing the intron 1 C allele. Gluteal muscle biopsies of TBs showed a significant association of the intron 1 C allele and SINE with a higher proportion of Type 2B and lower proportion of Type 1 fibers. However, in the Belgian horse, in which the SINE is not present, the intron 1 SNP was not associated with fiber type proportions, and evaluation of fiber type proportions across the Belgian, TB and QH breeds shows the significant effect of breed on fiber type proportions is negated when evaluating horses without the SINE variant. These data suggest the SINE, rather than the intron 1 SNP, is driving the observed muscle fiber type characteristics and is the variant targeted by selection for short‐distance racing.     

Does heterozygosity at the DMRT3 gene make French trotters better racers?

Background

Recently, a mutation was discovered in the DMRT3 gene that controls pacing in horses. The mutant allele A is fixed in the American Standardbred trotter breed, while in the French trotter breed, the frequency of the wild-type allele C is still 24%. This study aimed at measuring the effect of DMRT3 genotypes on the performance of French trotters and explaining why the polymorphism still occurs in this breed. Using a mixed animal model, genetic parameters and environmental effects on performance traits were estimated from data on 173 176 French trotter races. The effect of the DMRT3 gene was then estimated by the effect of genotype at the highly linked SNP BIEC2-620109 (C-C, A-T) for 630 horses. A selection scheme that included qualification and racing performances was modeled to (1) verify if the observed superiority of heterozygous CT horses at this SNP could be explained only by selection and (2) understand why allele C has not disappeared in French trotters.

Results

Heritability of racing performance traits was high for qualification test (0.56), moderate for annual earnings per finished race (0.26 to 0.31) and low for proportion of disqualified races (0.06 to 0.09). Genotype CC was always unfavorable compared to genotype TT for qualification: the probability to be qualified was 20% for CC vs. 48% for TT and earnings were -0.96 σ_y lower for CC than for TT. Genotype CT was also unfavorable for qualification (40%) and earnings at 3 years (-0.21 σ_y), but favorable for earnings at ages greater than 5 years: +0.41 σ_y (P = 7.10⁻⁴). Selection on qualification could not explain more than 19% of the difference between genotypes CC and CT in earnings at ages greater than 5 years. Only a scenario for which genotype CT has a favorable effect on the performance of horses older than 5 years could explain that the polymorphism at the DMRT3 gene still exists in the French trotter breed.

Conclusions

The use of mature horses in the French racing circuit can explain that the CA genotype is still present in the French trotter horses.     

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.     

基于基因组SNPs对东亚家马不同群体遗传多样性的评估

家马地方品种的遗传多样性是动态演变的, 并与其育种模式息息相关; 在过去约300年中, 育种者培育了具突出表型特征的标准化培育品种, 该类型品种对家马地方品种的繁育产生了巨大影响, 这是导致家马地方品种遗传多样性下降及遗传分化的主要因素之一。本研究采集了5个东亚家马地方品种(蒙古马、哈萨克马、河曲马、藏马和西南马)、2个西亚家马品种(阿拉伯马和阿哈尔捷金马)以及2个欧洲家马品种(设特兰矮马和克莱斯黛尔马)的70个样品, 并整合之前所发表的100匹内蒙古地区家马不同群体单核苷酸多态性(SNPs)数据集, 通过全基因组重测序和生物信息学方法分析了东亚家马不同群体的遗传多样性。研究发现东亚家马地方品种具有较丰富的遗传多样性, 相比欧洲和西亚品种产生了显著的遗传分化, 其中蒙古高原群体的遗传多样性最为丰富; 受杂交改良影响, 内蒙古不同群体间产生了一定程度的遗传分化; 河曲马和藏马的遗传背景最为单一, 受引种杂交繁育影响较小。本研究评估了东亚家马地方品种的遗传多样性分布格局和演化特征, 可为建立家马地方种质核心保护群体以及培育新品种提供遗传学理论支持。     

A missense mutation in the endothelin-B receptor gene is associated with Lethal White Foal Syndrome: an equine version of Hirschsprung Disease

Lethal White Foal Syndrome is a disease associated with horse breeds that register white coat spotting patterns. Breedings between particular spotted horses, generally described as frame overo, produce some foals that, in contrast to their parents, are all white or nearly all white and die shortly after birth of severe intestinal blockage. These foals have aganglionosis characterized by a lack of submucosal and myenteric ganglia from the distal small intestine to the large intestine, similar to human Hirschsprung Disease. Some sporadic and familial cases of Hirschsprung Disease are due to mutations in the endothelin B receptor gene (EDNRB). In this study, we investigate the role of EDNRB in Lethal White Foal Syndrome. A cDNA for the wild-type horse endothelin-B receptor gene was cloned and sequenced. In three unrelated lethal white foals, the EDNRB gene contained a 2-bp nucleotide change leading to a missense mutation (I118K) in the first transmembrane domain of the receptor, a highly conserved region of this protein among different species. Seven additional unrelated lethal white foal samples were found to be homozygous for this mutation. No other homozygotes were identified in 138 samples analyzed, suggesting that homozygosity was restricted to lethal white foals. All (40/40) horses with the frame overo pattern (a distinct coat color pattern that is a subset of overo horses) that were tested were heterozygous for this allele, defining a heterozygous coat color phenotype for this mutation. Horses with tobiano markings included some carriers, indicating that tobiano is epistatic to frame overo. In addition, horses were identified that were carriers but had no recognized overo coat pattern phenotype, demonstrating the variable penetrance of the mutation. The test for this mutant allele can be utilized in all breeds where heterozygous animals may be unknowingly bred to each other including the Paint Horse, Pinto horse, Quarter Horse, Miniature Horse, and Thoroughbred. Received: 25 November 1997 / Accepted: 3 February 1998     

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.     

Allelic heterogeneity at the equine KIT locus in dominant white (W) horses

White coat color has been a highly valued trait in horses for at least 2,000 years. Dominant white (W) is one of several known depigmentation phenotypes in horses. It shows considerable phenotypic variation, ranging from ~50% depigmented areas up to a completely white coat. In the horse, the four depigmentation phenotypes roan, sabino, tobiano, and dominant white were independently mapped to a chromosomal region on ECA 3 harboring the KIT gene. KIT plays an important role in melanoblast survival during embryonic development. We determined the sequence and genomic organization of the ~82 kb equine KIT gene. A mutation analysis of all 21 KIT exons in white Franches-Montagnes Horses revealed a nonsense mutation in exon 15 (c.2151C>G, p.Y717X). We analyzed the KIT exons in horses characterized as dominant white from other populations and found three additional candidate causative mutations. Three almost completely white Arabians carried a different nonsense mutation in exon 4 (c.706A>T, p.K236X). Six Camarillo White Horses had a missense mutation in exon 12 (c.1805C>T, p.A602V), and five white Thoroughbreds had yet another missense mutation in exon 13 (c.1960G>A, p.G654R). Our results indicate that the dominant white color in Franches-Montagnes Horses is caused by a nonsense mutation in the KIT gene and that multiple independent mutations within this gene appear to be responsible for dominant white in several other modern horse populations.     

Identification of Genetic Variation on the Horse Y Chromosome and the Tracing of Male Founder Lineages in Modern Breeds

The paternally inherited Y chromosome displays the population genetic history of males. While modern domestic horses (Equus caballus) exhibit abundant diversity within maternally inherited mitochondrial DNA, no significant Y-chromosomal sequence diversity has been detected. We used high throughput sequencing technology to identify the first polymorphic Y-chromosomal markers useful for tracing paternal lines. The nucleotide variability of the modern horse Y chromosome is extremely low, resulting in six haplotypes (HT), all clearly distinct from the Przewalski horse (E. przewalskii). The most widespread HT1 is ancestral and the other five haplotypes apparently arose on the background of HT1 by mutation or gene conversion after domestication. Two haplotypes (HT2 and HT3) are widely distributed at high frequencies among modern European horse breeds. Using pedigree information, we trace the distribution of Y-haplotype diversity to particular founders. The mutation leading to HT3 occurred in the germline of the famous English Thoroughbred stallion “Eclipse” or his son or grandson and its prevalence demonstrates the influence of this popular paternal line on modern sport horse breeds. The pervasive introgression of Thoroughbred stallions during the last 200 years to refine autochthonous breeds has strongly affected the distribution of Y-chromosomal variation in modern horse breeds and has led to the replacement of autochthonous Y chromosomes. Only a few northern European breeds bear unique variants at high frequencies or fixed within but not shared among breeds. Our Y-chromosomal data complement the well established mtDNA lineages and document the male side of the genetic history of modern horse breeds and breeding practices.