Horse body size varies greatly due to intense selection within each breed. American Miniatures are less than one meter tall at the withers while Shires and Percherons can exceed two meters. The genetic basis for this variation is not known. We hypothesize that the breed population structure of the horse should simplify efforts to identify genes controlling size. In support of this, here we show with genome-wide association scans (GWAS) that genetic variation at just four loci can explain the great majority of horse size variation. Unlike humans, which are naturally reproducing and possess many genetic variants with weak effects on size, we show that horses, like other domestic mammals, carry just a small number of size loci with alleles of large effect. Furthermore, three of our horse size loci contain the LCORL, HMGA2 and ZFAT genes that have previously been found to control human height. The LCORL/NCAPG locus is also implicated in cattle growth and HMGA2 is associated with dog size. Extreme size diversification is a hallmark of domestication. Our results in the horse, complemented by the prior work in cattle and dog, serve to pinpoint those very few genes that have played major roles in the rapid evolution of size during domestication.

The appaloosa coat spotting pattern in horses is caused by a single incomplete dominant gene (LP). Homozygosity for LP (LP/LP) is directly associated with congenital stationary night blindness (CSNB) in Appaloosa horses. LP maps to a 6-cM region on ECA1. We investigated the relative expression of two functional candidate genes located in this LP candidate region (TRPM1 and OCA2), as well as three other linked loci (TJP1, MTMR10, and OTUD7A) by quantitative real-time RT-PCR. No large differences were found for expression levels of TJP1, MTMR10, OTUD7A, and OCA2. However, TRPM1 (Transient Receptor Potential Cation Channel, Subfamily M, Member 1) expression in the retina of homozygous appaloosa horses was 0.05% the level found in non-appaloosa horses (R ¼ 0.0005). This constitutes a .1800-fold change (FC) decrease in TRPM1 gene expression in the retina (FC ¼ À1870.637, P ¼ 0.001) of CSNB-affected (LP/LP) horses. TRPM1 was also downregulated in LP/LP pigmented skin (R ¼ 0.005, FC ¼ À193.963, P ¼ 0.001) and in LP/LP unpigmented skin (R ¼ 0.003, FC ¼ À288.686, P ¼ 0.001) and was downregulated to a lesser extent in LP/lp unpigmented skin (R ¼ 0.027, FC ¼ À36.583, P ¼ 0.001). TRP proteins are thought to have a role in controlling intracellular Ca 21 concentration. Decreased expression of TRPM1 in the eye and the skin may alter bipolar cell signaling as well as melanocyte function, thus causing both CSNB and LP in horses.

Recent advances in genomic sequencing technology and computational assembly methods have allowed scientists to improve reference genome assemblies in terms of contiguity and composition. EquCab2, a reference genome for the domestic horse, was released in 2007. Although of equal or better quality compared to other first-generation Sanger assemblies, it had many of the shortcomings common to them. In 2014, the equine genomics research community began a project to improve the reference sequence for the horse, building upon the solid foundation of EquCab2 and incorporating new short-read data, long-read data, and proximity ligation data. Here, we present EquCab3. The count of non-N bases in the incorporated chromosomes is improved from 2.33 Gb in EquCab2 to 2.41 Gb in EquCab3. Contiguity has also been improved nearly 40-fold with a contig N50 of 4.5 Mb and scaffold contiguity enhanced to where all but one of the 32 chromosomes is comprised of a single scaffold.

During fetal development neural-crest-derived melanoblasts migrate across the entire body surface and differentiate into melanocytes, the pigment-producing cells. Alterations in this precisely regulated process can lead to white spotting patterns. White spotting patterns in horses are a complex trait with a large phenotypic variance ranging from minimal white markings up to completely white horses. The “splashed white” pattern is primarily characterized by an extremely large blaze, often accompanied by extended white markings at the distal limbs and blue eyes. Some, but not all, splashed white horses are deaf. We analyzed a Quarter Horse family segregating for the splashed white coat color. Genome-wide linkage analysis in 31 horses gave a positive LOD score of 1.6 in a region on chromosome 6 containing the PAX3 gene. However, the linkage data were not in agreement with a monogenic inheritance of a single fully penetrant mutation. We sequenced the PAX3 gene and identified a missense mutation in some, but not all, splashed white Quarter Horses. Genome-wide association analysis indicated a potential second signal near MITF. We therefore sequenced the MITF gene and found a 10 bp insertion in the melanocyte-specific promoter. The MITF promoter variant was present in some splashed white Quarter Horses from the studied family, but also in splashed white horses from other horse breeds. Finally, we identified two additional non-synonymous mutations in the MITF gene in unrelated horses with white spotting phenotypes. Thus, several independent mutations in MITF and PAX3 together with known variants in the EDNRB and KIT genes explain a large proportion of horses with the more extreme white spotting phenotypes.

BackgroundTo date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array.ResultsUsing whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation.ConclusionsHere, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3943-8) contains supplementary material, which is available to authorized users.