Genome-wide linkage and association analysis identifies major gene Loci for guttural pouch tympany in arabian and german warmblood horses

Equine guttural pouch tympany (GPT) is a hereditary condition affecting foals in their first months of life. Complex segregation analyses in Arabian and German warmblood horses showed the involvement of a major gene as very likely. Genome-wide linkage and association analyses including a high density marker set of single nucleotide polymorphisms (SNPs) were performed to map the genomic region harbouring the potential major gene for GPT. A total of 85 Arabian and 373 German warmblood horses were genotyped on the Illumina equine SNP50 beadchip. Non-parametric multipoint linkage analyses showed genome-wide significance on horse chromosomes (ECA) 3 for German warmblood at 16-26 Mb and 34-55 Mb and for Arabian on ECA15 at 64-65 Mb. Genome-wide association analyses confirmed the linked regions for both breeds. In Arabian, genome-wide association was detected at 64 Mb within the region with the highest linkage peak on ECA15. For German warmblood, signals for genome-wide association were close to the peak region of linkage at 52 Mb on ECA3. The odds ratio for the SNP with the highest genome-wide association was 0.12 for the Arabian. In conclusion, the refinement of the regions with the Illumina equine SNP50 beadchip is an important step to unravel the responsible mutations for GPT.     

A Genome-Wide Association Study Identifies Risk Loci to Equine Recurrent Uveitis in German Warmblood Horses

Equine recurrent uveitis (ERU) is a common eye disease affecting up to 3–15% of the horse population. A genome-wide association study (GWAS) using the Illumina equine SNP50 bead chip was performed to identify loci conferring risk to ERU. The sample included a total of 144 German warmblood horses. A GWAS showed a significant single nucleotide polymorphism (SNP) on horse chromosome (ECA) 20 at 49.3 Mb, with IL-17A and IL-17F being the closest genes. This locus explained a fraction of 23% of the phenotypic variance for ERU. A GWAS taking into account the severity of ERU, revealed a SNP on ECA18 nearby to the crystalline gene cluster CRYGA-CRYGF. For both genomic regions on ECA18 and 20, significantly associated haplotypes containing the genome-wide significant SNPs could be demonstrated. In conclusion, our results are indicative for a genetic component regulating the possible critical role of IL-17A and IL-17F in the pathogenesis of ERU. The associated SNP on ECA18 may be indicative for cataract formation in the course of ERU.     

Linked markers exclude KIT as the gene responsible for appaloosa coat colour spotting patterns in horses

The appaloosa coat colour pattern of the horse is similar to that caused by the rump-white (Rw) gene in the mouse. In the mouse Rw colour pattern is the result of an inversion involving the proto-oncogene c-kit (KIT). Therefore, we investigated KIT as a candidate gene that encodes the appaloosa coat colour gene (Lp) in horses. KIT plays a critical role in haematopoiesis, gametogenesis, and melanogenesis and encodes a transmembrane tyrosine kinase receptor that belongs to the PDGF/CSF-1/c-KIT receptor subfamily. Half-sib families segregating for Lp were uninformative for a reported polymorphism in KIT. However, KIT is located on horse chromosome 3 close to albumin (ALB), serum carboxylesterase (ES), vitamin D-binding protein (GC) and microsatellite markers ASB23, LEX007, LEX57, and UCDEQ437. Indeed, KIT and ASB23 were localized to ECA3q21-22.1 and 3q22.1-22.3, respectively, by fluorescent in situ hybridization. Family studies were conducted to investigate linkage of Lp to these markers using eight half-sib families in which Appaloosa stallions were mated to solid coloured mares. Linkage of Lp to the chromosome region containing ES, ALB, GC, ASB23, UCDEQ437, LEX57, and LEX007 was investigated by a multipoint linkage analysis using the computer program GENEHUNTER. LOD scores over the interval under investigation ranged from -4.28 to -12.48, with a score of -12.48 at the location for ASB23. Therefore, it was concluded that appaloosa (Lp) is not linked to any of the tested markers on ECA3, and thus Lp is unlikely to be the product of KIT.     

Analysis of genetic variation contributing to measured speed in Thoroughbreds identifies genomic regions involved in the transcriptional response to exercise

Despite strong selection for athletic traits in Thoroughbred horses, there is marked variation in speed and aptitude for racing performance within the breed. Using global positioning system monitoring during exercise training, we measured speed variables and temporal changes in speed with age to derive phenotypes for GWAS. The aim of the study was to test the hypothesis that genetic variation contributes to variation in end‐point physiological traits, in this case galloping speed measured during field exercise tests. Standardisation of field‐measured phenotypes was attempted by assessing horses exercised on the same gallop track and managed under similar conditions by a single trainer. PCA of six key speed indices captured 73.9% of the variation with principal component 1 (PC1). Verifying the utility of the phenotype, we observed that PC1 (median) in 2‐year‐old horses was significantly different among elite, non‐elite and unraced horses (P < 0.001) and the temporal change with age in PC1 varied among horses with different myostatin (MSTN) g.66493737C>T SNP genotypes. A GWAS for PC1 in 2‐year‐old horses (n = 122) identified four SNPs reaching the suggestive threshold for association (P < 4.80 × 10^?5), defining a 1.09 Mb candidate region on ECA8 containing the myosin XVIIIB (MYO18B) gene. In a GWAS for temporal change in PC1 with age (n = 168), five SNPs reached the suggestive threshold for association and defined candidate regions on ECA2 and ECA11. Both regions contained genes that are significantly differentially expressed in equine skeletal muscle in response to acute exercise and training stimuli, including MYO18A. As MYO18A plays a regulatory role in the skeletal muscle response to exercise, the identified genomic variation proximal to the myosin family genes may be important for the regulation of the response to exercise and training.     

Physical anchorage and orientation of equine linkage groups by FISH mapping BAC clones containing microsatellite markers

A horse bacterial artificial chromosome (BAC) library was screened for 19 microsatellite markers from unassigned or non-oriented linkage groups. Clones containing 11 (AHT20, EB2E8, HMS45, LEX005, LEX014, LEX023, LEX044, TKY111, UCDEQ425, UCDEQ464 and VIASH21) of these were found, which were from eight different linkage groups. The BAC clones were used as probes in dual colour FISH to identify their precise chromosomal origin. The microsatellite markers are located on nine different horse chromosomes, four of which (ECA6, ECA25, ECA27 and ECA28) had no previously in situ assigned markers.     

A genome-wide association study for racing performances in Thoroughbreds clarifies a candidate region near the MSTN gene

Using 1400 microsatellites, a genome-wide association study (GWAS) was performed to identify genomic regions associated with lifetime earnings and performance ranks, as determined by the Japan Racing Association (JRA). The minimum heritability (h(2) ) was estimated at 7-8% based on the quantitative trait model, suggesting that the racing performance is heritable. Following GWAS with microsatellites, fine mapping led to identification of three SNPs on ECA18, namely, g.65809482T>C (P=1.05E-18), g.65868604G>T (P=6.47E-17), and g.66539967A>G (P=3.35E-14) associated with these performance measures. The haplotype of these SNPs, together with a recently published nearby SNP, g.66493737C>T (P=9.06E-16) in strong linkage disequilibrium, also showed a very clear association with the performance (P<1E-05). The candidate genomic region contained eight genes annotated by ENSEMBL, including the myostatin gene (MSTN). These findings suggest the presence of a gene affecting the racing performance in Thoroughbred racehorses in this region on ECA18.     

Refinement of quantitative trait loci on equine chromosome 10 for radiological signs of navicular disease in Hanoverian warmblood horses

Navicular disease is characterized by a progressive degenerative alteration of the equine podotrochlea. In this study, we refined a previously identified quantitative trait locus (QTL) on horse chromosome 10 for the abnormal development of canales sesamoidales (DCS) of the navicular bone in Hanoverian warmblood horses. Genotyping was done in 192 Hanoverian warmblood horses from 17 paternal half-sib groups. The whole marker set comprised 45 markers including seven newly developed microsatellites and 13 single nucleotide polymorphisms (SNPs) within positional candidate genes. Chromosome-wide significant QTL were confirmed and refined for DCS on horse chromosome (ECA) 10 at 0.16-2.70 Mb and at 14.45-36.37 Mb. Nine microsatellites and three SNP markers reached the highest multipoint Zmeans and LOD scores at 19.34-20.38 Mb and at 23.17-30.73 Mb with genome-wide error probabilities of P<0.05. In addition, a significant association of a SNP within VSTM1 and a significant haplotype-trait association within IRF3 could be shown. These results support a possible role of the candidate genes VSTM1 and IRF3 within the QTL on ECA10 for DCS. This study is a further step towards the identification of the genes responsible for navicular disease in Hanoverian warmblood horses.     

Whole-genome scan identifies quantitative trait loci for chronic pastern dermatitis in German draft horses

Chronic pastern dermatitis (CPD), also known as chronic progressive lymphedema (CPL), is a skin disease that affects draft horses. This disease causes painful lower-leg swelling, nodule formation, and skin ulceration, interfering with movement. The aim of this whole-genome scan was to identify quantitative trait loci (QTL) for CPD in German draft horses. We recorded clinical data for CPD in 917 German draft horses and collected blood samples from these horses. Of these 917 horses, 31 paternal half-sib families comprising 378 horses from the breeds Rhenish German, Schleswig, Saxon-Thuringian, and South German were chosen for genotyping. Each half-sib family was constituted by only one draft horse breed. Genotyping was done for 318 polymorphic microsatellites evenly distributed on all equine autosomes and the X chromosome with a mean distance of 7.5 Mb. An across-breed multipoint linkage analysis revealed chromosome-wide significant QTL on horse chromosomes (ECA) 1, 9, 16, and 17. Analyses by breed confirmed the QTL on ECA1 in South German and the QTL on ECA9, 16, and 17 in Saxon-Thuringian draft horses. For the Rhenish German and Schleswig draft horses, additional QTL on ECA4 and 10 and for the South German draft horses an additional QTL on ECA7 were found. This is the first whole-genome scan for CPD in draft horses and it is an important step toward the identification of candidate genes.     

Refinement of a quantitative gene locus on equine chromosome 16 responsible for osteochondrosis in Hanoverian warmblood horses

Osteochondrosis (OC) is an inherited developmental disease in young horses most frequently observed in thoroughbreds, trotters, warmblood and coldblood horses. Quantitative trait loci (QTL) for equine OC have been identified in Hanoverian warmblood horses employing a whole genome scan with microsatellites. A QTL on ECA16 reached the genome-wide significance level for hock osteochondrosis dissecans (OCD). The aim of this study was to refine this QTL on ECA16 using an extended marker set of 34 newly developed microsatellites and 15 single nucleotide polymorphisms (SNPs). We used the same 14 paternal half-sib groups as in the above-mentioned whole genome scan. The QTL for OCD in hock joints on ECA16 could be delimited at an interval between 17.60 and 45.18 Mb using multipoint non-parametric linkage analyses. In addition, six microsatellites and one SNP were significantly associated with hock OCD in the QTL region between 24.26 and 42.41 Mb. Furthermore, our analysis revealed a second QTL for fetlock OC between 6.55 and 24.26 Mb on ECA16. This report is a further step towards unravelling the genes underlying QTL for equine OC and towards the development of a marker test for OC in Hanoverian warmblood horses.     

Genome-wide association analysis of osteochondrosis of the tibiotarsal joint in Norwegian Standardbred trotters

Osteochondrosis (OC), a disturbance in the process of endochondral ossification, is by far the most important equine developmental orthopaedic disease and is also common in other domestic animals and humans. The purpose of this study was to identify quantitative trait loci (QTL) associated with osteochondrosis dissecans (OCD) at the intermediate ridge of the distal tibia in Norwegian Standardbred (SB) using the Illumina Equine SNP50 BeadChip whole-genome single-nucleotide polymorphism (SNP) assay. Radiographic data and blood samples were obtained from 464 SB yearlings. Based on the radiographic examination, 162 horses were selected for genotyping; 80 of these were cases with an OCD at the intermediate ridge of the distal tibia, and 82 were controls without any developmental lesions in the joints examined. Genotyped horses descended from 22 sires, and the number of horses in each half-sib group ranged from 3 to 14. The population structure necessitated statistical correction for stratification. When conducting a case-control genome-wide association study (GWAS), mixed-model analyses displayed regions on chromosomes (Equus callabus chromosome - ECA) 5, 10, 27 and 28 that showed moderate evidence of association (P ≤ 5 × 10(-5); this P-value is uncorrected i.e. not adjusted for multiple comparisons) with OCD in the tibiotarsal joint. Two SNPs on ECA10 represent the most significant hits (uncorrected P=1.19 × 10(-5) in the mixed-model). In the basic association (chi-square) test, these SNPs achieved statistical significance with the Bonferroni correction (P=0.038) and were close in the permuted logistic regression test (P=0.054). Putative QTL on ECA 5, 10, 27 and 28 represent interesting areas for future research, validation studies and fine mapping of candidate regions. Results presented here represent the first GWAS of OC in horses using the recently released Illumina Equine SNP50 BeadChip.     

Fine mapping of a quantitative trait locus for osteochondrosis on horse chromosome 2

In this study, we refine a quantitative trait locus for equine osteochondrosis (OC) on horse chromosome (ECA) 2 to a genome-wide significant interval at 20.08-30.94 Mb. The marker set contained 27 newly developed microsatellites equidistantly distributed over ECA2 and 44 nucleotide polymorphisms, located in 16 positional candidate genes for OC. Genotyping was performed in 211 Hanoverian horses from 14 paternal half-sib groups. A NCDN-associated SNP and haplotype were significantly associated with OC in fetlock and/or hock joints. This study is a further step towards the identification of genes responsible for OC in horses.     

A genome-wide association study for quantitative trait loci of show-jumping in Hanoverian warmblood horses

Show-jumping is an economically important breeding goal in Hanoverian warmblood horses. The aim of this study was a genome-wide association study (GWAS) for quantitative trait loci (QTL) for show-jumping in Hanoverian warmblood horses, employing the Illumina equine SNP50 Beadchip. For our analyses, we genotyped 115 stallions of the National State stud of Lower Saxony. The show-jumping talent of a horse includes style and ability in free-jumping. To control spurious associations based on population stratification, two different mixed linear animal model (MLM) approaches were employed, besides linear models with fixed effects only and adaptive permutations for correcting multiple testing. Population stratification was explained best in the MLM considering Hanoverian, Thoroughbred, Trakehner and Holsteiner genes and the marker identity-by-state relationship matrix. We identified six QTL for show-jumping on horse chromosomes (ECA) 1, 8, 9 and 26 (-log(10) P-value >5) and further putative QTL with -log(10) P-values of 3-5 on ECA1, 3, 11, 17 and 21. Within six QTL regions, we identified human performance-related genes including PAPSS2 on ECA1, MYL2 on ECA8, TRHR on ECA9 and GABPA on ECA26 and within the putative QTL regions NRAP on ECA1, and TBX4 on ECA11. The results of our GWAS suggest that genes involved in muscle structure, development and metabolism are crucial for elite show-jumping performance. Further studies are required to validate these QTL in larger data sets and further horse populations.     

Refinement of a quantitative trait locus on equine chromosome 5 responsible for fetlock osteochondrosis in Hanoverian warmblood horses

In this report, we provide 29 new informative microsatellites distributed over a region of 21 Mb on horse chromosome (ECA) 5 and refine a quantitative trait locus (QTL) for fetlock osteochondrosis dissecans (OCD) to a genome-wide significant interval between 78.03 and 90.23 Mb on ECA5. Genotyping was performed in 211 Hanoverian warmblood horses from 14 paternal half-sib groups. Within this OCD-QTL, collagen type XXIV alpha 1 was identified as a potential functional candidate gene for equine osteochondrosis. This report is a further step towards unravelling the genes that cause equine osteochondrosis.     

Genome-wide association study of insect bite hypersensitivity in two horse populations in the Netherlands

Background

Insect bite hypersensitivity is a common allergic disease in horse populations worldwide. Insect bite hypersensitivity is affected by both environmental and genetic factors. However, little is known about genes contributing to the genetic variance associated with insect bite hypersensitivity. Therefore, the aim of our study was to identify and quantify genomic associations with insect bite hypersensitivity in Shetland pony mares and Icelandic horses in the Netherlands.

Methods

Data on 200 Shetland pony mares and 146 Icelandic horses were collected according to a matched case–control design. Cases and controls were matched on various factors (e.g. region, sire) to minimize effects of population stratification. Breed-specific genome-wide association studies were performed using 70 k single nucleotide polymorphisms genotypes. Bayesian variable selection method Bayes-C with a threshold model implemented in GenSel software was applied. A 1 Mb non-overlapping window approach that accumulated contributions of adjacent single nucleotide polymorphisms was used to identify associated genomic regions.

Results

The percentage of variance explained by all single nucleotide polymorphisms was 13% in Shetland pony mares and 28% in Icelandic horses. The 20 non-overlapping windows explaining the largest percentages of genetic variance were found on nine chromosomes in Shetland pony mares and on 14 chromosomes in Icelandic horses. Overlap in identified associated genomic regions between breeds would suggest interesting candidate regions to follow-up on. Such regions common to both breeds (within 15 Mb) were found on chromosomes 3, 7, 11, 20 and 23. Positional candidate genes within 2 Mb from the associated windows were identified on chromosome 20 in both breeds. Candidate genes are within the equine lymphocyte antigen class II region, which evokes an immune response by recognizing many foreign molecules.

Conclusions

The genome-wide association study identified several genomic regions associated with insect bite hypersensitivity in Shetland pony mares and Icelandic horses. On chromosome 20, associated genomic regions in both breeds were within 2 Mb from the equine lymphocyte antigen class II region. Increased knowledge on insect bite hypersensitivity associated genes will contribute to our understanding of its biology, enabling more efficient selection, therapy and prevention to decrease insect bite hypersensitivity prevalence.     

Identification of equine major histocompatibility complex haplotypes using polymorphic microsatellites

A system for identifying equine major histocompatibility complex (MHC) haplotypes was developed based on five polymorphic microsatellites located within the MHC region on ECA 20. Molecular signatures for 50 microsatellite haplotypes were recognized from typing 353 horses. Of these, 23 microsatellite haplotypes were associated with 12 established equine leucocyte antigen (ELA) haplotypes in Thoroughbreds and Standardbreds. Five ELA serotypes were associated with multiple microsatellite subhaplotypes, expanding the estimates of diversity in the equine MHC. The strong correlations between serological and microsatellite typing demonstrated a linkage to known MHC class I protein polymorphisms and validated this assay as a useful supplement to ELA serotyping, and in some applications, a feasible alternative method for MHC genotyping in horse families and in population studies.     

Homozygosity mapping approach identifies a missense mutation in equine cyclophilin B (PPIB) associated with HERDA in the American Quarter Horse

Hereditary equine regional dermal asthenia (HERDA), a degenerative skin disease that affects the Quarter Horse breed, was localized to ECA1 by homozygosity mapping. Comparative genomics allowed the development of equine gene-specific markers which were used with a set of affected horses to detect a homozygous, identical-by-descent block spanning approximately 2.5 Mb, suggesting a recent origin for the HERDA mutation. We report a mutation in cyclophilin B (PPIB) as a novel, causal candidate gene for HERDA. A c.115G>A missense mutation in PPIB alters a glycine residue that has been conserved across vertebrates. The mutation was homozygous in 64 affected horses and segregates concordant with inbreeding loops apparent in the genealogy of 11 affected horses. Screening of control Quarter Horses indicates a 3.5% carrier frequency. The development of a test that can detect affected horses prior to development of clinical signs and carriers of HERDA will allow Quarter Horse breeders to eliminate this debilitating disease.     

QTL‐seq for rapid identification of candidate genes for 100‐seed weight and root/total plant dry weight ratio under rainfed conditions in chickpea

Terminal drought is a major constraint to chickpea productivity. Two component traits responsible for reduction in yield under drought stress include reduction in seeds size and root length/root density. QTL‐seq approach, therefore, was used to identify candidate genomic regions for 100‐seed weight (100SDW) and total dry root weight to total plant dry weight ratio (RTR) under rainfed conditions. Genomewide SNP profiling of extreme phenotypic bulks from the ICC 4958 × ICC 1882 population identified two significant genomic regions, one on CaLG01 (1.08 Mb) and another on CaLG04 (2.7 Mb) linkage groups for 100SDW. Similarly, one significant genomic region on CaLG04 (1.10 Mb) was identified for RTR. Comprehensive analysis revealed four and five putative candidate genes associated with 100SDW and RTR, respectively. Subsequently, two genes (Ca_04364 and Ca_04607) for 100SDW and one gene (Ca_04586) for RTR were validated using CAPS/dCAPS markers. Identified candidate genomic regions and genes may be useful for molecular breeding for chickpea improvement.     

A genome-wide association study reveals loci influencing height and other conformation traits in horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ～200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ～40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ～18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.     

Genome-wide search for microsatellite markers associated with radiologic alterations in the navicular bone of Hanoverian warmblood horses

The aim of this study was to identify quantitative trait loci (QTLs) for pathologic changes in the navicular bone in Hanoverian warmblood horses. Seventeen paternal half-sib groups comprising 192 individuals were analyzed in a whole-genome scan. These families included 144 progeny and grandchildren, which were randomly chosen from the Hanoverian warmblood. Three different traits were considered: deformed canales sesamoidales and radiographic changes in the contour and in the structure of the navicular bone. The genome scan included in total 214 highly polymorphic microsatellite markers. The putatively linked genomic regions on equine chromosomes (ECA) 2, 3, 10, and 15 were refined using 53 additional microsatellites. Chromosome-wide significant QTLs were located on five different equine chromosomes (ECA2, 3, 4, 10, and 26). Genome-wide significant QTLs were on ECA2 at 48 cM and on ECA10 from 45.5 to 49.8 cM. This study was a first step to get more insight into the molecular genetic determination of radiologic changes in the equine navicular bone. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetics of swayback in American Saddlebred horses

Extreme lordosis, also called swayback, lowback or softback, can occur as a congenital trait or as a degenerative trait associated with ageing. In this study, the hereditary aspect of congenital swayback was investigated using whole genome association studies of 20 affected and 20 unaffected American Saddlebred (ASB) Horses for 48,165 single-nucleotide polymorphisms (SNPs). A statistically significant association was identified on ECA20 (corrected P=0.017) for SNP BIEC2-532523. Of the 20 affected horses, 17 were homozygous for this SNP when compared to seven homozygotes among the unaffected horses, suggesting a major gene with a recessive mode of inheritance. The result was confirmed by testing an additional 13 affected horses and 166 unaffected horses using 35 SNPs in this region of ECA20 (corrected P=0.036). Combined results for 33 affected horses and 287 non-affected horses allowed identification of a region of homozygosity defined by four SNPs in the region. Based on the haplotype defined by these SNPs, 80% of the 33 affected horses were homozygous, 21% heterozygous and 9% did not possess the haplotype. Among the non-affected horses, 15% were homozygous, 47% heterozygous and 38% did not possess the haplotype. The differences between the two groups were highly significant (P<0.00001). The region defined by this haplotype includes 53 known and predicted genes. Exons from three candidate genes, TRERF1, RUNX2 and CNPY3 were sequenced without finding distinguishing SNPs. The mutation responsible for swayback may lie in other genes or in regulatory regions outside exons. This information can be used by breeders to reduce the occurrence of swayback among their livestock. This condition may serve as a model for investigation of congenital skeletal deformities in other species.