Background: Equine Multiple Congenital Ocular Anomalies (MCOA) syndrome consists of a diverse set of abnormalities predominantly localized to the frontal part of the eye. The disease is in agreement with a codominant mode of inheritance in our horse material. Animals presumed to be heterozygous for the mutant allele have cysts originating from the temporal ciliary body, peripheral retina and/or iris. In contrast, animals predicted to be homozygous for the disease-causing allele possess a wide range of multiple abnormalities, including iridociliary and/or peripheral retinal cysts, iridocorneal angle abnormalities, cornea globosa, iris hypoplasia and congenital cataracts. MCOA is most common in the Rocky Mountain horse breed where it occurs at a high frequency among Silver colored horses. The Silver coat color is associated with mutations inPMEL17that resides on ECA6q23. To map theMCOAlocus we analyzed 11 genetic markers on ECA6q and herein describe a chromosome interval for theMCOAlocus.Results: We performed linkage analysis within 17 paternal half-sib families of the Rocky Mountain horse breed. More than half of the 131 offspring had the Cyst phenotype and about one third had MCOA. Segregation data were obtained by genotyping 10 microsatellite markers most of which are positioned on ECA6q22-23, as well as the missense mutation for the Silver phenotype inPMEL17. Significant linkage was found between theMCOAlocus and eight of the genetic markers, where markerUPP5(Theta = 0, z = 12.3),PMEL17ex11(Theta = 0, z = 19.0) andUPP6(Theta = 0, z = 17.5) showed complete linkage with theMCOAlocus. DNA sequencing ofPMEL17in affected and healthy control individuals did not reveal any additional mutations than the two mutations associated with the Silver coat color.Conclusion: TheMCOAlocus can with high confidence be positioned within a 4.9 megabase (Mb) interval on ECA6q. The genotype data onUPP5,PMEL17ex11andUPP6strongly support the hypothesis that horses with the Cyst phenotype are heterozygous for the mutant allele and that horses with the MCOA phenotype are homozygous for the mutant allele.

Curly horses present a variety of curl phenotypes that are associated with various degrees of curliness of coat, mane, tail and ear hairs. Their origin is still a matter of debate and several genetic hypotheses have been formulated to explain the diversity in phenotype, including the combination of autosomal dominant and recessive alleles. Our purpose was to map the autosomal dominant curly hair locus and identify the causal variant using genome-wide association study (GWAS) and whole-genome sequencing approaches. A GWAS was performed using a Bayesian sparse linear mixed model, based on 51 curly and 19 straight-haired French and North American horses from 13 paternal families genotyped on the Illumina EquineSNP50 BeadChip. A single strong signal was observed on equine chromosome 11, in a region that encompasses the type I keratin gene cluster. This region was refined by haplotype analysis to a segment including 36 genes, among which are 10 keratin genes (KRT-10, -12, -20, -23, -24, -25, -26, -27, -28, -222). To comprehensively identify candidate causal variants within all these genes, whole-genome sequences were obtained for one heterozygous curly stallion and its straight-haired son. Among the four non-synonymous candidate variants identified and validated in the curly region, only variant g.21891160G>A in the KRT25 gene (KRT25:p.R89H) was in perfect agreement with haplotype status in the whole pedigree. Genetic association was then confirmed by genotyping a larger population consisting of 353 horses. However, five discordant curly horses were observed, which carried neither the variant nor the main haplotype associated with curliness. Sequencing of KRT25 for two discordant horses did not identify any other deleterious variant, which suggests locus rather than allelic heterogeneity for the curly phenotype. We identified the KRT25:p.R89H variant as responsible for the dominant curly trait, but a second dominant locus may also be involved in the shape of hairs within North American Curly horses.

As part of the requirements of the Convention on Biological Diversity, Canada has been investigating the genetic diversity of its native equine and pony populations. Along with examining four indigenous Canadian equine populations (Canadian horse, Lac La Croix pony, Newfoundland pony and Sable Island population), another 10 Mountain and Moorland, three Nordic, four horse and two feral equine populations (thought to have influenced some pony breeds) were also investigated. In total, 821 individuals were genotyped at 38 microsatellite loci. Results of the analysis of molecular variance indicated that 13.3% of genetic diversity was explained by breed differences, whereas 84.6% and 2.1% of diversity came from within and among individuals, respectively. The average effective number of alleles and allelic richness was the lowest in the Eriskay (2.51 and 3.98) and Lac La Croix (2.83 and 4.01) populations, whereas it was highest in the New Forest (4.31 and 6.01) and Welsh (4.33 and 5.87) breeds, followed closely by the Newfoundland-CDN (4.23 and 5.86) population. Expected heterozygosities varied from 0.61 in the Lac La Croix to 0.74 in the Welsh and in Newfoundland. Observed heterozygosities ranged from 0.57 in the Exmoor and 0.58 in the Sable Island herd to 0.77 in the Kerry Bog and 0.76 in the New Forest breeds. Structure and admixture analyses revealed that the most likely number of clusters was 21, although some substructure was also observed when K = 16, compared with the 24 predefined populations. Information gathered from this study should be combined with other available phenotypic and pedigree data to develop, or amend, a suitable conservation strategy for all populations examined.

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

The Estonian Native Horse (ENH) is a medium-size pony found mainly in the western islands of Estonia and is well-adapted to the harsh northern climate and poor pastures. The ancestry of the ENH is debated, including alleged claims about direct descendance from the extinct Tarpan. Here we conducted a detailed analysis of the genetic makeup and relationships of the ENH based on the genotypes of 15 autosomal short tandem repeats (STRs), 18 Y chromosomal single nucleotide polymorphisms (SNPs), mitochondrial D-loop sequence and lateral gait allele in DMRT3. The study encompassed 2890 horses of 61 breeds, including 33 ENHs. We show that the expected and observed genetic diversities of the ENH are among the highest within 52 global breeds, and the highest among 8 related Northern European ponies. The genetically closest breeds to the ENH are the Finn Horse, and the geographically more distant primitive Hucul and Konik. ENH matrilines are diverse and relate to draught and Pontic-Caspian breeds. ENH patrilines relate to draught breeds, and to a unique haplogroup not described before. None of the 33 ENHs carried the “gait” mutation, but the mutation was found in 2 Huculs. The study demonstrates that the ENH is a genetically distinct and diverse breed of ancient origin with no notable pressure of selective breeding.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski’s horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.