Performance of Swedish Warmblood fragile foal syndrome carriers and breeding prospects.

BACKGROUND: Warmblood fragile foal syndrome (WFFS) is a monogenetic defect caused by a recessive lethal missense point mutation in the procollagen-lysine, 2-oxoglutarate 5-dioxygenase 1 gene (PLOD1, c.2032G>A). The majority of homozygous WFFS horses are aborted during gestation. Clinical signs of affected horses include fragile skin, skin and mucosa lacerations, hyperextension of the articulations, and hematomas. In spite of its harmful effect, a relatively high frequency of WFFS carriers has been found in Warmblood horses, suggesting a heterozygote advantage. Thus, in this study our aims were to: (1) estimate the frequency of WFFS carriers in the Swedish Warmblood breed (SWB), (2) estimate the effect of WFFS carrier genotype on performance traits in two SWB subpopulations bred for different disciplines, and (3) simulate the potential effects of balancing selection and different selection strategies on the frequency of carriers. METHODS: In total, 2288 SWB sport horses born between 1971 and 2020 were tested for the WFFS mutation and had estimated breeding values (EBV) for ten traditional evaluating and 50 linear descriptive traits. RESULTS: The frequency of WFFS carriers calculated from a pool of 511 randomly selected SWB horses born in 2017 was equal to 7.4% and ranged from 0.0 to 12.0% among the whole set of tested SWB horses, starting from 1971 till 2020. The effect of the WFFS carrier genotype was significant for several EBV mainly related to movements and dressage traits and especially for horses not bred for the show jumping discipline. Using simulation, we showed that balancing selection can maintain a recessive lethal allele in populations such as the SWB breed over generations and that the frequency is expected to slowly decrease in absence of balancing selection. Finally, we showed that selection against carrier sires can result in a more rapid decrease of the frequency of the mutant allele over time. CONCLUSION: Further research is needed to confirm the apparent association between equine performance and the WFFS carrier genotype. Identification of such associations or new causative mutations for horse performance traits can serve as new tools in horse breeding to select for healthy, sustainable, and better performing horses.

Using droplet digital PCR for the detection of hco-acr-8b levamisole resistance marker in H. contortus.

The nematode Haemonchus contortus is one of the most prevalent and pathogenic parasites in small ruminants. Although usually controlled using anthelmintics, the development of drug resistance by the parasite has become a major issue in livestock production. While the molecular detection of benzimidazole resistance in H. contortus is well developed, the molecular tools and protocols are far less advanced for the detection of levamisole resistance. The hco-acr-8 gene encodes a critical acetylcholine susceptible subunit that confers levamisole-sensitivity to the receptor. Here, we report the development of a droplet digital PCR assay as a molecular tool to detect a 63 bp deletion in the hco-acr-8 that has been previously associated with levamisole resistance. Sanger sequencing of single adult H. contortus yielded 56 high-quality consensus sequences surrounding the region containing the deletion. Based on the sequencing data, new primers and probes were designed and validated with a novel droplet digital PCR assay for the quantification of the deletion containing "resistant" allele in genomic DNA samples. Single adult worms from six phenotypically described isolates (n = 60) and from two Swedish sheep farms (n = 30) where levamisole was effective were tested. Even though a significant difference in genotype frequencies between the resistant and susceptible reference isolates was found (p = 0.01), the homozygous "resistant" genotype was observed to be abundantly present in both the susceptible isolates as well as in some Swedish H. contortus samples. Furthermore, field larval culture samples, collected pre- (n = 7) and post- (n = 6) levamisole treatment on seven Swedish sheep farms where levamisole was fully efficacious according to Fecal Egg Count Reduction Test results, were tested to evaluate the frequency of the "resistant" allele in each. Frequencies of the deletion ranged from 35 to 80% in the pre-treatment samples, whereas no amplifiable H. contortus genomic DNA was detected in the post-treatment samples. Together, these data reveal relatively high frequencies of the 63 bp deletion in the hco-acr-8 both on individual H. contortus and field larval culture scales, and cast doubt on the utility of the deletion in the hco-acr-8 as a molecular marker for levamisole resistance detection on sheep farms.

A genome-wide scan for candidate lethal variants in Thoroughbred horses.

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

Genetic Diversity and Signatures of Selection in a Native Italian Horse Breed Based on SNP Data.

Horses are nowadays mainly used for sport and leisure activities, and several local breeds, traditionally used in agriculture, have been exposed to a dramatic loss in population size and genetic diversity. The loss of genetic diversity negatively impacts individual fitness and reduces the potential long-term survivability of a breed. Recent advances in molecular biology and bioinformatics have allowed researchers to explore biodiversity one step further. This study aimed to evaluate the loss of genetic variability and identify genomic regions under selection pressure in the Bardigiano breed based on GGP Equine70k SNP data. The effective population size based on Linkage Disequilibrium (Ne) was equal to 39 horses, and it showed a decline over time. The average inbreeding based on runs of homozygosity (ROH) was equal to 0.17 (SD = 0.03). The majority of the ROH were relatively short (91% were ≤ 2Mbp long), highlighting the occurrence of older inbreeding, rather than a more recent occurrence. A total of eight ROH islands, shared among more than 70% of the Bardigiano horses, were found. Four of them mapped to known quantitative trait loci related to morphological traits (e.g., body size and coat color) and disease susceptibility. This study provided the first genome-wide scan of genetic diversity and selection signatures in an Italian native horse breed.

Genomic Divergence in Swedish Warmblood Horses Selected for Equestrian Disciplines.

The equestrian sport horse Swedish Warmblood (SWB) originates from versatile cavalry horses. Most modern SWB breeders have specialized their breeding either towards show jumping or dressage disciplines. The aim of this study was to explore the genomic structure of SWB horses to evaluate the presence of genomic subpopulations, and to search for signatures of selection in subgroups of SWB with high or low breeding values (EBVs) for show jumping. We analyzed high density genotype information from 380 SWB horses born in the period 2010-2011, and used Principal Coordinates Analysis and Discriminant Analysis of Principal Components to detect population stratification. Fixation index and Cross Population Extended Haplotype Homozygosity scores were used to scan the genome for potential signatures of selection. In accordance with current breeding practice, this study highlights the development of two separate breed subpopulations with putative signatures of selection in eleven chromosomes. These regions involve genes with known function in, e.g., mentality, endogenous reward system, development of connective tissues and muscles, motor control, body growth and development. This study shows genetic divergence, due to specialization towards different disciplines in SWB horses. This latter evidence can be of interest for SWB and other horse studbooks encountering specialized breeding.

Inter- and intra-breed genome-wide copy number diversity in a large cohort of European equine breeds.

BACKGROUND: Copy Number Variation (CNV) is a common form of genetic variation underlying animal evolution and phenotypic diversity across a wide range of species. In the mammalian genome, high frequency of CNV differentiation between breeds may be candidates for population-specific selection. However, CNV differentiation, selection and its population genetics have been poorly explored in horses. RESULTS: We investigated the patterns, population variation and gene annotation of CNV using the Axiom® Equine Genotyping Array (670,796 SNPs) from a large cohort of individuals (N = 1755) belonging to eight European horse breeds, varying from draught horses to several warmblood populations. After quality control, 152,640 SNP CNVs (individual markers), 18,800 segment CNVs (consecutive SNP CNVs of same gain/loss state or both) and 939 CNV regions (CNVRs; overlapping segment CNVs by at least 1 bp) compared to the average signal of the reference (Belgian draught horse) were identified. Our analyses showed that Equus caballus chromosome 12 (ECA12) was the most enriched in segment CNV gains and losses (~ 3% average proportion of the genome covered), but the highest number of segment CNVs were detected on ECA1 and ECA20 (regardless of size). The Friesian horses showed private SNP CNV gains (> 20% of the samples) on ECA1 and Exmoor ponies displayed private SNP CNV losses on ECA25 (> 20% of the samples). The Warmblood cluster showed private SNP CNV gains located in ECA9 and Draught cluster showed private SNP CNV losses located in ECA7. The length of the CNVRs ranged from 1 kb to 21.3 Mb. A total of 10,612 genes were annotated within the CNVRs. The PANTHER annotation of these genes showed significantly under- and overrepresented gene ontology biological terms related to cellular processes and immunity (Bonferroni P-value < 0.05). We identified 80 CNVRs overlapping with known QTL for fertility, coat colour, conformation and temperament. We also report 67 novel CNVRs. CONCLUSIONS: This work revealed that CNV patterns, in the genome of some European horse breeds, occurred in specific genomic regions. The results provide support to the hypothesis that high frequency private CNVs residing in genes may potentially be responsible for the diverse phenotypes seen between horse breeds.

Genomic relatedness and diversity of Swedish native cattle breeds.

BACKGROUND: Native cattle breeds are important genetic resources given their adaptation to the local environment in which they are bred. However, the widespread use of commercial cattle breeds has resulted in a marked reduction in population size of several native cattle breeds worldwide. Therefore, conservation management of native cattle breeds requires urgent attention to avoid their extinction. To this end, we genotyped nine Swedish native cattle breeds with genome-wide 150 K single nucleotide polymorphisms (SNPs) to investigate the level of genetic diversity and relatedness between these breeds. RESULTS: We used various SNP-based approaches on this dataset to connect the demographic history with the genetic diversity and population structure of these Swedish cattle breeds. Our results suggest that the Väne and Ringamåla breeds originating from southern Sweden have experienced population isolation and have a low genetic diversity, whereas the Fjäll breed has a large founder population and a relatively high genetic diversity. Based on the shared ancestry and the constructed phylogenetic trees, we identified two major clusters in Swedish native cattle. In the first cluster, which includes Swedish mountain cattle breeds, there was little differentiation among the Fjäll, Fjällnära, Swedish Polled, and Bohus Polled breeds. The second cluster consists of breeds from southern Sweden: Väne, Ringamåla and Swedish Red. Interestingly, we also identified sub-structuring in the Fjällnära breed, which indicates different breeding practices on the farms that maintain this breed. CONCLUSIONS: This study represents the first comprehensive genome-wide analysis of the genetic relatedness and diversity in Swedish native cattle breeds. Our results show that different demographic patterns such as genetic isolation and cross-breeding have shaped the genomic diversity of Swedish native cattle breeds and that the Swedish mountain breeds have retained their authentic distinct gene pool without significant contribution from any of the other European cattle breeds that were included in this study.

Signatures of selection in the genome of Swedish warmblood horses selected for sport performance.

BACKGROUND: A growing demand for improved physical skills and mental attitude in modern sport horses has led to strong selection for performance in many warmblood studbooks. The aim of this study was to detect genomic regions with low diversity, and therefore potentially under selection, in Swedish Warmblood horses (SWB) by analysing high-density SNP data. To investigate if such signatures could be the result of selection for equestrian sport performance, we compared our SWB SNP data with those from Exmoor ponies, a horse breed not selected for sport performance traits. RESULTS: The genomic scan for homozygous regions identified long runs of homozygosity (ROH) shared by more than 85% of the genotyped SWB individuals. Such ROH were located on ECA4, ECA6, ECA7, ECA10 and ECA17. Long ROH were instead distributed evenly across the genome of Exmoor ponies in 77% of the chromosomes. Two population differentiation tests (FST and XP-EHH) revealed signatures of selection on ECA1, ECA4, and ECA6 in SWB horses. CONCLUSIONS: Genes related to behaviour, physical abilities and fertility, appear to be targets of selection in the SWB breed. This study provides a genome-wide map of selection signatures in SWB horses, and ground for further functional studies to unravel the biological mechanisms behind complex traits in horses.

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology.

The mammalian Major Histocompatibility Complex (MHC) region contains several gene families characterized by highly polymorphic loci with extensive nucleotide diversity, copy number variation of paralogous genes, and long repetitive sequences. This structural complexity has made it difficult to construct a reliable reference sequence of the horse MHC region. In this study, we used long-read single molecule, real-time (SMRT) sequencing technology from Pacific Biosciences (PacBio) to sequence eight Bacterial Artificial Chromosome (BAC) clones spanning the horse MHC class II region. The final assembly resulted in a 1,165,328 bp continuous gap free sequence with 35 manually curated genomic loci of which 23 were considered to be functional and 12 to be pseudogenes. In comparison to the MHC class II region in other mammals, the corresponding region in horse shows extraordinary copy number variation and different relative location and directionality of the Eqca-DRB, -DQA, -DQB and -DOB loci. This is the first long-read sequence assembly of the horse MHC class II region with rigorous manual gene annotation, and it will serve as an important resource for association studies of immune-mediated equine diseases and for evolutionary analysis of genetic diversity in this region.

A Nonsense Variant in the ST14 Gene in Akhal-Teke Horses with Naked Foal Syndrome.

Naked foal syndrome (NFS) is a genodermatosis in the Akhal-Teke horse breed. We provide the first scientific description of this phenotype. Affected horses have almost no hair and show a mild ichthyosis. So far, all known NFS affected horses died between a few weeks and 3 yr of age. It is not clear whether a specific pathology caused the premature deaths. NFS is inherited as a monogenic autosomal recessive trait. We mapped the disease causing genetic variant to two segments on chromosomes 7 and 27 in the equine genome. Whole genome sequencing of two affected horses, two obligate carriers, and 75 control horses from other breeds revealed a single nonsynonymous genetic variant on the chromosome 7 segment that was perfectly associated with NFS. The affected horses were homozygous for ST14:c.388G>T, a nonsense variant that truncates >80% of the open reading frame of the ST14 gene (p.Glu130*). The variant leads to partial nonsense-mediated decay of the mutant transcript. Genetic variants in the ST14 gene are responsible for autosomal recessive congenital ichthyosis 11 in humans. Thus, the identified equine ST14:c.388G>T variant is an excellent candidate causative variant for NFS, and the affected horses represent a large animal model for a known human genodermatosis. Our findings will enable genetic testing to avoid the nonintentional breeding of NFS-affected foals.

Large Deletions at the SHOX Locus in the Pseudoautosomal Region Are Associated with Skeletal Atavism in Shetland Ponies.

Skeletal atavism in Shetland ponies is a heritable disorder characterized by abnormal growth of the ulna and fibula that extend the carpal and tarsal joints, respectively. This causes abnormal skeletal structure and impaired movements, and affected foals are usually killed. In order to identify the causal mutation we subjected six confirmed Swedish cases and a DNA pool consisting of 21 control individuals to whole genome resequencing. We screened for polymorphisms where the cases and the control pool were fixed for opposite alleles and observed this signature for only 25 SNPs, most of which were scattered on genome assembly unassigned scaffolds. Read depth analysis at these loci revealed homozygosity or compound heterozygosity for two partially overlapping large deletions in the pseudoautosomal region (PAR) of chromosome X/Y in cases but not in the control pool. One of these deletions removes the entire coding region of the SHOX gene and both deletions remove parts of the CRLF2 gene located downstream of SHOX. The horse reference assembly of the PAR is highly fragmented, and in order to characterize this region we sequenced bacterial artificial chromosome (BAC) clones by single-molecule real-time (SMRT) sequencing technology. This considerably improved the assembly and enabled size estimations of the two deletions to 160-180 kb and 60-80 kb, respectively. Complete association between the presence of these deletions and disease status was verified in eight other affected horses. The result of the present study is consistent with previous studies in humans showing crucial importance of SHOX for normal skeletal development.

Conformation Traits and Gaits in the Icelandic Horse are Associated with Genetic Variants in Myostatin (MSTN).

Many genes are known to have an influence on conformation and performance traits; however, the role of one gene, Myostatin (MSTN), has been highlighted in recent studies on horses. Myostatin acts as a repressor in the development and regulation of differentiation and proliferative growth of skeletal muscle. Several studies have examined the link between MSTN, conformation, and performance in racing breeds, but no studies have investigated the relationship in Icelandic horses. Icelandic horses, a highly unique breed, are known both for their robust and compact conformation as well as their additional gaits tölt and pace. Three SNPs (g.65868604G>T [PR8604], g.66493737C>T [PR3737], and g.66495826A>G [PR5826]) flanking or within equine MSTN were genotyped in 195 Icelandic horses. The SNPs and haplotypes were analyzed for association with official estimated breeding values (EBV) for conformation traits (n = 11) and gaits (n = 5). The EBV for neck, withers, and shoulders was significantly associated with both PR8604 and PR3737 (P < 0.05). PR8604 was also associated with EBV for total conformation (P = 0.05). These associations were all supported by the haplotype analysis. However, while SNP PR5826 showed a significant association with EBVs for leg stance and hooves (P < 0.05), haplotype analyses for these traits failed to fully support these associations. This study demonstrates the possible role of MSTN on both the form and function of horses from non-racing breeds. Further analysis of Icelandic horses as well as other non-racing breeds would be beneficial and likely help to completely understand the influence of MSTN on conformation and performance in horses.

Using an Inbred Horse Breed in a High Density Genome-Wide Scan for Genetic Risk Factors of Insect Bite Hypersensitivity (IBH).

While susceptibility to hypersensitive reactions is a common problem amongst humans and animals alike, the population structure of certain animal species and breeds provides a more advantageous route to better understanding the biology underpinning these conditions. The current study uses Exmoor ponies, a highly inbred breed of horse known to frequently suffer from insect bite hypersensitivity, to identify genomic regions associated with a type I and type IV hypersensitive reaction. A total of 110 cases and 170 controls were genotyped on the 670K Axiom Equine Genotyping Array. Quality control resulted in 452,457 SNPs and 268 individuals being tested for association. Genome-wide association analyses were performed using the GenABEL package in R and resulted in the identification of two regions of interest on Chromosome 8. The first region contained the most significant SNP identified, which was located in an intron of the DCC netrin 1 receptor gene. The second region identified contained multiple top SNPs and encompassed the PIGN, KIAA1468, TNFRSF11A, ZCCHC2, and PHLPP1 genes. Although additional studies will be needed to validate the importance of these regions in horses and the relevance of these regions in other species, the knowledge gained from the current study has the potential to be a step forward in unraveling the complex nature of hypersensitive reactions.

Regulatory mutations in TBX3 disrupt asymmetric hair pigmentation that underlies Dun camouflage color in horses.

Dun is a wild-type coat color in horses characterized by pigment dilution with a striking pattern of dark areas termed primitive markings. Here we show that pigment dilution in Dun horses is due to radially asymmetric deposition of pigment in the growing hair caused by localized expression of the T-box 3 (TBX3) transcription factor in hair follicles, which in turn determines the distribution of hair follicle melanocytes. Most domestic horses are non-dun, a more intensely pigmented phenotype caused by regulatory mutations impairing TBX3 expression in the hair follicle, resulting in a more circumferential distribution of melanocytes and pigment granules in individual hairs. We identified two different alleles (non-dun1 and non-dun2) causing non-dun color. non-dun2 is a recently derived allele, whereas the Dun and non-dun1 alleles are found in ancient horse DNA, demonstrating that this polymorphism predates horse domestication. These findings uncover a new developmental role for T-box genes and new aspects of hair follicle biology and pigmentation.

Evaluation of whole-genome sequencing of four Chinese crested dogs for variant detection using the ion proton system.

BACKGROUND: Next generation sequencing (NGS) has traditionally been performed by large genome centers, but in recent years, the costs for whole-genome sequencing (WGS) have decreased substantially. With the introduction of smaller and less expensive "desktop" systems, NGS is now moving into the general laboratory. To evaluate the Ion Proton system for WGS we sequenced four Chinese Crested dogs and analyzed the data quality in terms of genome and exome coverage, the number of detected single nucleotide variants (SNVs) and insertions and deletions (INDELs) and the genotype concordance with the Illumina HD canine SNP array. For each of the four dogs, a 200 bp fragment library was constructed from genomic DNA and sequenced on two Ion PI chips per dog to reach mean coverage of 6-8x of the canine genome (genome size ≈ 2.4 Gb). RESULTS: On average, each Ion PI chip yielded approximately 73.3 million reads with a mean read length of 130 bp (~9.5 Gb sequence data) of which 98.5 % could be aligned to the canine reference genome (CanFam3.1). By sequencing a single dog using one fragment library and two Ion PI chips, on average 80 % of the genome and 77 % exome was covered by at least four reads. After removing duplicate reads (20.7 %) the mean coverage across the whole genome was 6x. Using sequence data from all four individuals (four fragment libraries and eight Ion PI chips) the genome and exome coverage could be further increased to 97.2 and 94.3 %, respectively. We detected 4.83 million unique SNPs and 6.10 million unique INDEL positions across all individuals. A comparison between SNP genotypes detected with the WGS and the 170 K Illumina HD canine SNP array showed 90 % concordance. CONCLUSIONS: We have evaluated whole-genome sequencing on the Ion Proton system for genetic variant detection in four Chinese crested dogs. Even though INDEL calling with Ion Proton data is challenging due to specific platform errors, in case of SNP calling it can serve as an alternative to other next-generation sequencing platforms and SNP genotyping arrays, in studies aiming to identify causative mutations for rare monogenic diseases. In addition, we have identified new genetic variants of the Chinese Crested dog that will contribute to further whole-genome sequencing studies aimed to identify mutations associated with monogenic diseases with autosomal recessive inheritance.

Genetic diversity of five local Swedish chicken breeds detected by microsatellite markers.

This study aimed at investigating the genetic diversity, relationship and population structure of 110 local Swedish chickens derived from five breeds (Gotlandshöna, Hedemorahöna, Öländsk dvärghöna, Skånsk blommehöna, and Bohuslän- Dals svarthöna, in the rest of the paper the shorter name Svarthöna is used) using 24 microsatellite markers. In total, one hundred thirteen alleles were detected in all populations, with a mean of 4.7 alleles per locus. For the five chicken breeds, the observed and expected heterozygosity ranged from 0.225 to 0.408 and from 0.231 to 0.515, with the lowest scores for the Svarthöna and the highest scores for the Skånsk blommehöna breeds, respectively. Similarly, the average within breed molecular kinship varied from 0.496 to 0.745, showing high coancestry, with Skånsk blommehöna having the lowest and Svarthöna the highest coancestry. Furthermore, all breeds showed significant deviations from Hardy-Weinberg expectations. Across the five breeds, the global heterozygosity deficit (FIT) was 0.545, population differentiation index (FST) was 0.440, and the global inbreeding of individuals within breed (FIS) was 0.187. The phylogenetic relationships of chickens were examined using neighbor-joining trees constructed at the level of breeds and individual samples. The neighbor-joining tree constructed at breed level revealed two main clusters, with Hedemorahöna and Öländsk dvärghöna breeds in one cluster, and Gotlandshöna and Svarthöna breeds in the second cluster leaving the Skånsk blommehöna in the middle. Based on the results of the STRUCTURE analysis, the most likely number of clustering of the five breeds was at K = 4, with Hedemorahöna, Gotlandshöna and Svarthöna breeds forming their own distinct clusters, while Öländsk dvärghöna and Skånsk blommehöna breeds clustered together. Losses in the overall genetic diversity of local Swedish chickens due to breeds extinction varied from -1.46% to -6.723%. The results of the current study can be used as baseline genetic information for genetic conservation program, for instance, to control inbreeding and to implement further genetic studies in local Swedish chickens.

Genetic diversity in the modern horse illustrated from genome-wide SNP data.

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. F(ST) 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.

Genome-wide analysis reveals selection for important traits in domestic horse breeds.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.

Mutations in DMRT3 affect locomotion in horses and spinal circuit function in mice.

Locomotion in mammals relies on a central pattern-generating circuitry of spinal interneurons established during development that coordinates limb movement. These networks produce left-right alternation of limbs as well as coordinated activation of flexor and extensor muscles. Here we show that a premature stop codon in the DMRT3 gene has a major effect on the pattern of locomotion in horses. The mutation is permissive for the ability to perform alternate gaits and has a favourable effect on harness racing performance. Examination of wild-type and Dmrt3-null mice demonstrates that Dmrt3 is expressed in the dI6 subdivision of spinal cord neurons, takes part in neuronal specification within this subdivision, and is critical for the normal development of a coordinated locomotor network controlling limb movements. Our discovery positions Dmrt3 in a pivotal role for configuring the spinal circuits controlling stride in vertebrates. The DMRT3 mutation has had a major effect on the diversification of the domestic horse, as the altered gait characteristics of a number of breeds apparently require this mutation.