Illumina Equine SNP50 BeadChip Research Articles

A Mechanogenetic Model of Exercise-Induced Pulmonary Haemorrhage in the Thoroughbred Horse.

Exercise-induced pulmonary haemorrhage (EIPH) occurs in horses performing high-intensity athletic activity. The application of physics principles to derive a ‘physical model’, which is coherent with existing physiology and cell biology data, shows that critical parameters for capillary rupture are cell–cell adhesion and cell stiffness (cytoskeleton organisation). Specifically, length of fracture in the capillary is a ratio between the energy involved in cell–cell adhesion and the stiffness of cells suggesting that if the adhesion diminishes and/or that the stiffness of cells increases EIPH is more likely to occur. To identify genes associated with relevant cellular or physiological phenotypes, the physical model was used in a post-genome-wide association study (GWAS) to define gene sets associated with the model parameters. The primary study was a GWAS of EIPH where the phenotype was based on weekly tracheal wash samples collected over a two-year period from 72 horses in a flat race training yard. The EIPH phenotype was determined from cytological analysis of the tracheal wash samples, by scoring for the presence of red blood cells and haemosiderophages. Genotyping was performed using the Illumina Equine SNP50 BeadChip and analysed using linear regression in PLINK. Genes within significant genome regions were selected for sets based on their GeneOntology biological process, and analysed using fastBAT. The gene set analysis showed that genes associated with cell stiffness (cytoskeleton organisation) and blood flow have the most significant impact on EIPH risk.

Selection signatures in four German warmblood horse breeds: Tracing breeding history in the modern sport horse.

The study of selection signatures helps to find genomic regions that have been under selective pressure and might host genes or variants that modulate important phenotypes. Such knowledge improves our understanding of how breeding programmes have shaped the genomes of livestock. In this study, 942 stallions were included from four, exemplarily chosen, German warmblood breeds with divergent historical and recent selection focus and different crossbreeding policies: Trakehner (N = 44), Holsteiner (N = 358), Hanoverian (N = 319) and Oldenburger (N = 221). Those breeds are nowadays bred for athletic performance and aptitude for show-jumping, dressage or eventing, with a particular focus of Holsteiner on the first discipline. Blood samples were collected during the health exams of the stallion preselections before licensing and were genotyped with the Illumina EquineSNP50 BeadChip. Autosomal markers were used for a multi-method search for signals of positive selection. Analyses within and across breeds were conducted by using the integrated Haplotype Score (iHS), cross-population Extended Haplotype Homozygosity (xpEHH) and Runs of Homozygosity (ROH). Oldenburger and Hanoverian showed very similar iHS signatures, but breed specificities were detected on multiple chromosomes with the xpEHH. The Trakehner clustered as a distinct group in a principal component analysis and also showed the highest number of ROHs, which reflects their historical bottleneck. Beside breed specific differences, we found shared selection signals in an across breed iHS analysis on chromosomes 1, 4 and 7. After investigation of these iHS signals and shared ROH for potential functional candidate genes and affected pathways including enrichment analyses, we suggest that genes affecting muscle functionality (TPM1, TMOD2-3, MYO5A, MYO5C), energy metabolism and growth (AEBP1, RALGAPA2, IGFBP1, IGFBP3-4), embryonic development (HOXB-complex) and fertility (THEGL, ZPBP1-2, TEX14, ZP1, SUN3 and CFAP61) have been targeted by selection in all breeds. Our findings also indicate selection pressure on KITLG, which is well-documented for influencing pigmentation.

Genotype imputation accuracy in multiple equine breeds from medium- to high-density genotypes.

Genotype imputation is now a key component of genomic analyses as it increases the density of available genotypes within a population. However, many factors can influence imputation accuracy. The aim of this study was to assess and compare the accuracy of imputation of high-density genotypes (Affymetrix Axiom Equine genotyping array, 670,806 SNPs) from two moderate-density genotypes (Illumina Equine SNP50 BeadChip, 54,602 SNPs and Illumina Equine SNP70 BeadChip, 65,157 SNPs), using single-breed or multiple-breed reference sets. Genotypes were available from five groups of horse breeds: Arab (AR, 1,207 horses), Trotteur Français (TF, 979 horses), Selle Français (SF, 1,979 horses), Anglo-Arab (AA, 229 horses) and various foreign sport horses (FH, 209 horses). The proportions of horses genotyped with the high-density (HD) chip in each breed group were 10% in AA, 15% in AR and FH, 30% in TF and 57% in SF. A validation set consisting of one-third of the horses genotyped with the HD chip was formed and their genotypes deleted. Two imputation strategies were compared, one in which the reference population consisted only of horses from the same breed group as in the validation set, and another with horses from all breed groups. For the first strategy, concordance rates (CRs) ranged from 97.8% (AR) to 99.0% (TF) and correlations (r²) from 0.94 (AR) to 0.99 (TF). For the second strategy, CR ranged from 97.4% (AR) to 98.9% (TF) and r² from 0.93 (AR) to 0.99 (TF). Overall, the results show a small advantage of within-breed imputation compared with multi-breed imputation. Adding horses from different breed groups to the reference population does not improve the accuracy of imputation. Imputation provides an accurate means of combining data sets from different genotyping platforms, now necessary with the increasing use of the recently developed Affymetrix Axiom Equine genotyping array.

Genetic risk factors for osteochondrosis in various horse breeds.

Osteochondrosis (OC) is an injury to cartilage canals with a following necrosis in the growth cartilage, from there it can develop to osteochondrosis dissecans (OCD). Due to its high impact in the equine industry, new insights into predisposing factors and potential high-risk genetic variants are warranted. This article reviews advancements in quantitative and molecular genetics in refining estimation of genetic parameters and identifying predisposing genetic loci. Heritabilities were highest for hock OC with estimates at 0.29-0.46 in Hanoverian warmblood and Norwegian trotters, whereas in Thoroughbreds only very low genetic variation seemed to be present in hock OC lesions. Whole genome scans using the Illumina Equine SNP50 or SNP70 Beadchip were performed in Thoroughbred, Standardbred, French and Norwegian trotter, Hanoverian and Dutch warmblood. Validation studies in Spanish Purebred and Hanoverian warmblood horses corroborated OC risk loci on ECA 3, 14, 27 and 29. Particularly, a strong association with hock-OCD was found for a single nucleotide polymorphism (SNP) on horse chromosome (ECA) 3 upstream to the LCORL gene. Gene expression and microRNA analyses may be helpful to understand pathophysiological processes in equine OC and to connect OCD-associated genomic regions with potential candidate genes. Furthermore progress in elucidating the underlying genetic variants and pathophysiological changes in OC may be expected from whole genome DNA and RNA next-generation sequencing studies.

A missense variant in the coil1A domain of the keratin 25 gene is associated with the dominant curly hair coat trait (Crd) in horse

BackgroundCurly 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.ResultsA 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.ConclusionsWe 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.

Genome-wide association study for semen quality traits in German Warmblood stallions

We performed a genome-wide association study for semen quality traits in 139 German Warmblood stallions. Stallions were genotyped using the Illumina equine SNP50 Beadchip. Traits analysed were de-regressed estimated breeding values (EBVs) for gel-free volume, sperm concentration, total number of sperm, progressive motility and the total number of progressively motile sperm. The GWAS revealed 29 SNPs on 12 different chromosomes as genome-wide significantly associated with semen quality traits. For ten genomic regions we could retrieve candidate genes influencing stallion fertility. Among the candidate genes, we could find the genes encoding cysteine-rich secretory proteins (CRISP1, CRISP2 and CRISP3). This was the first GWAS in horses performed for semen quality traits.

Genome-wide association study for jumping performances in French sport horses.

A genome-wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo-Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54,602 SNPs available on this chip, 44,424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single-marker mixed model and a haplotype-based mixed model (significant: P < 1E-05; suggestive: P < 1E-04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.

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.

A genome-wide association study indicatesLCORL/NCAPGas a candidate locus for withers height in German Warmblood horses

A genome-wide association scan for loci affecting withers height was conducted in 782 German Warmblood stallions, which were genotyped using the Illumina EquineSNP50 Bead Chip. A principal components approach was applied to correct for population structure. The analysis revealed a single major QTL on ECA3 explaining ~18 per cent of the phenotypic variance, which is in concordance with recent reports from other horse populations. The LCORL/NCAPG locus represents a strong candidate gene for this QTL. This locus is among a small number that have consistently been identified to influence human height in several large meta-analyses. Furthermore, a mutation within the NCAPG gene was found to affect growth and body frame size in cattle. Together with the results of this study in German Warmbloods, these findings strongly indicate LCORL/NCAPG as a candidate locus for withers height in horses. Further studies are, however, needed to confirm this.

Expression Levels of LCORL Are Associated with Body Size in Horses

Body size is an important characteristic for horses of various breeds and essential for the classification of ponies concerning the limit value of 148 cm (58.27 inches) height at the withers. Genome-wide association analyses revealed the highest associated quantitative trait locus for height at the withers on horse chromosome (ECA) 3 upstream of the candidate gene LCORL. Using 214 Hanoverian horses genotyped on the Illumina equine SNP50 BeadChip and 42 different horse breeds across all size ranges, we confirmed the highly associated single nucleotide polymorphism BIEC2-808543 (−log10P = 8.3) and the adjacent gene LCORL as the most promising candidate for body size. We investigated the relative expression levels of LCORL and its two neighbouring genes NCAPG and DCAF16 using quantitative real-time PCR (RT-qPCR). We could demonstrate a significant association of the relative LCORL expression levels with the size of the horses and the BIEC2-808543 genotypes within and across horse breeds. In heterozygous C/T-horses expression levels of LCORL were significantly decreased by 40% and in homozygous C/C-horses by 56% relative to the smaller T/T-horses. Bioinformatic analyses indicated that this SNP T>C mutation is disrupting a putative binding site of the transcription factor TFIID which is important for the transcription process of genes involved in skeletal bone development. Thus, our findings suggest that expression levels of LCORL play a key role for body size within and across horse breeds and regulation of the expression of LCORL is associated with genetic variants of BIEC2-808543. This is the first functional study for a body size regulating polymorphism in horses and a further step to unravel the mechanisms for understanding the genetic regulation of body size in horses.

Computation of deregressed proofs for genomic selection when own phenotypes exist with an application in French show-jumping horses1

Genomic evaluations often use as pseudo-phenotypes corrected means of progeny performances, like daughter yield deviations (DYD) in dairy species. In horse breeding, own performances are also available and performances from other relatives (as half sibs) may play an important part in the EBV because the number of progeny remains low, even for stallions. The first step for genomic selection in horses is therefore to generate pseudo-phenotypes for genomic analysis when parental or own information is considered. This work presents an easy method to compute deregressed EBV from regular pedigree-based genetic evaluations (EBV, reliabilities) to be used in genomic evaluations. The proposed methodology builds deregressed proofs so that they combine own performances (from genotyped individuals) and performances of relatives (outside of the genotyped sample). An application to show jumping horse data is presented. A sample of 908 stallions specialized in show jumping [71% Selle Français (SF), 17% foreign sport horses (FH), 13% Anglo Arab (AA)] were genotyped. Genotyping was performed using the Illumina Equine SNP50 BeadChip, and after quality tests, 44,444 SNP were retained. Two methods were used for genomic evaluation: GBLUP and BayesCπ, and 6 validation data sets were compared, chosen according to breeds SF + FH + AA or SF + FH, family structure (more than 3 half sibs), reliability of sires (>0.97) or sons (>0.72). In spite of a favorable genetic structure [linkage disequilibrium equal to 0.24 at 50 kb pairs], results showed low advantage of genomic evaluation. On the validation sample SF + FH + AA, the correlation between deregressed proofs and GBLUP or BayesCπ predictions was 0.39, 0.37, 0.51 according to the different validation data sets, compared with 0.36, 0.33, 0.53 obtained with BLUP predictions. Correlations were much lower on the SF + FH sample. Research is pursued to understand this low advantage of genomic selection and to improve the methodology for genomic evaluation in this context, which is less favorable than dairy cattle breeding.

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.

Detection of copy number variants in the horse genome and examination of their association with recurrent laryngeal neuropathy

We used the data from a recently performed genome-wide association study using the Illumina Equine SNP50 beadchip for the detection of copy number variants (CNVs) and examined their association with recurrent laryngeal neuropathy (RLN), an important equine upper airway disease compromising performance. A total of 2797 CNVs were detected for 477 horses, covering 229kb and seven SNPs on average. Overlapping CNVs were merged to define 478 CNV regions (CNVRs). CNVRs, particularly deletions, were shown to be significantly depleted in genes. Fifty-two of the 67 common CNVRs (frequency≥1%) were validated by association mapping, Mendelian inheritance, and/or Mendelian inconsistencies. None of the 67 common CNVRs were significantly associated with RLN when accounting for multiple testing. However, a duplication on chromosome 10 was detected in 10 cases (representing three breeds) and two unphenotyped parents but in none of the controls. The duplication was embedded in an 8-Mb haplotype shared across breeds.

Inbreeding in the Thoroughbred horse

Changes in the inbreeding coefficient, F, in the Thoroughbred horse over the past 45 years have been investigated by genotyping 467 Thoroughbred horses (born between 1961 and 2006) using the Illumina Equine SNP50 bead chip, which comprises 54,602 SNPs uniformly distributed across the equine genome. The Spearman rank correlation coefficient, r, between the year of birth and F was estimated. The results indicate that inbreeding in Thoroughbreds has increased over the past 40 years, with r = 0.24, P < 0.001 demonstrating that there is a highly significant, though relatively weak correlation between the year of birth and inbreeding coefficients. Interestingly, the majority of the increase in inbreeding is post-1996 and coincides with the introduction of stallions covering larger numbers of mares.

A genome-wide association study of osteochondritis dissecans in the Thoroughbred.

Osteochondrosis is a developmental orthopaedic disease that occurs in horses, other livestock species, companion animal species, and humans. The principal aim of this study was to identify quantitative trait loci (QTL) associated with osteochondritis dissecans (OCD) in the Thoroughbred using a genome-wide association study. A secondary objective was to test the effect of previously identified QTL in the current population. Over 300 horses, classified as cases or controls according to clinical findings, were genotyped for the Illumina Equine SNP50 BeadChip. An animal model was first implemented in order to adjust each horse's phenotypic status for average relatedness among horses and other potentially confounding factors which were present in the data. The genome-wide association test was then conducted on the residuals from the animal model. A single SNP on chromosome 3 was found to be associated with OCD at a genome-wide level of significance, as determined by permutation. According to the current sequence annotation, the SNP is located in an intergenic region of the genome. The effects of 24 SNPs, representing QTL previously identified in a sample of Hanoverian Warmblood horses, were tested directly in the animal model. When fitted alongside the significant SNP on ECA3, two of these SNPs were found to be associated with OCD. Confirmation of the putative QTL identified on ECA3 requires validation in an independent sample. The results of this study suggest that a significant challenge faced by equine researchers is the generation of sufficiently large data sets to effectively study complex diseases such as osteochondrosis.

Genome-wide association study among four horse breeds identifies a common haplotype associated with in vitro CD3+ T cell susceptibility/resistance to equine arteritis virus infection.

Previously, we have shown that horses could be divided into susceptible and resistant groups based on an in vitro assay using dual-color flow cytometric analysis of CD3+ T cells infected with equine arteritis virus (EAV). Here, we demonstrate that the differences in in vitro susceptibility of equine CD3+ T lymphocytes to EAV infection have a genetic basis. To investigate the possible hereditary basis for this trait, we conducted a genome-wide association study (GWAS) to compare susceptible and resistant phenotypes. Testing of 267 DNA samples from four horse breeds that had a susceptible or a resistant CD3+ T lymphocyte phenotype using both Illumina Equine SNP50 BeadChip and Sequenom's MassARRAY system identified a common, genetically dominant haplotype associated with the susceptible phenotype in a region of equine chromosome 11 (ECA11), positions 49572804 to 49643932. The presence of a common haplotype indicates that the trait occurred in a common ancestor of all four breeds, suggesting that it may be segregated among other modern horse breeds. Biological pathway analysis revealed several cellular genes within this region of ECA11 encoding proteins associated with virus attachment and entry, cytoskeletal organization, and NF-κB pathways that may be associated with the trait responsible for the in vitro susceptibility/resistance of CD3+ T lymphocytes to EAV infection. The data presented in this study demonstrated a strong association of genetic markers with the trait, representing de facto proof that the trait is under genetic control. To our knowledge, this is the first GWAS of an equine infectious disease and the first GWAS of equine viral arteritis.

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.

The same ELA class II risk factors confer equine insect bite hypersensitivity in two distinct populations

Insect bite hypersensitivity (IBH) is a chronic allergic dermatitis common in horses. Affected horses mainly react against antigens present in the saliva from the biting midges, Culicoides ssp, and occasionally black flies, Simulium ssp. Because of this insect dependency, the disease is clearly seasonal and prevalence varies between geographical locations. For two distinct horse breeds, we genotyped four microsatellite markers positioned within the MHC class II region and sequenced the highly polymorphic exons two from DRA and DRB3, respectively. Initially, 94 IBH-affected and 93 unaffected Swedish born Icelandic horses were tested for genetic association. These horses had previously been genotyped on the Illumina Equine SNP50 BeadChip, which made it possible to ensure that our study did not suffer from the effects of stratification. The second population consisted of 106 unaffected and 80 IBH-affected Exmoor ponies. We show that variants in the MHC class II region are associated with disease susceptibility (p raw = 2.34 × 10−5), with the same allele (COR112:274) associated in two separate populations. In addition, we combined microsatellite and sequencing data in order to investigate the pattern of homozygosity and show that homozygosity across the entire MHC class II region is associated with a higher risk of developing IBH (p = 0.0013). To our knowledge this is the first time in any atopic dermatitis suffering species, including man, where the same risk allele has been identified in two distinct populations.Electronic supplementary materialThe online version of this article (doi:10.1007/s00251-011-0573-1) contains supplementary material, which is available to authorized users.

The identification of SNPs with indeterminate positions using the Equine SNP50 BeadChip

We have used linkage disequilibrium (LD) to identify single nucleotide polymorphisms (SNPs) on the Illumina Equine SNP50 BeadChip, which may be incorrectly positioned on the genome map. A total of 1201 Thoroughbred horses were genotyped using the Illumina Equine SNP50 BeadChip. LD was evaluated in a pairwise fashion between all autosomal SNPs, both within and across chromosomes. Filters were then applied to the data, firstly to identify SNPs that may have been mapped to the wrong chromosome and secondly to identify SNPs that may have been incorrectly positioned within chromosomes. We identified a single SNP on ECA28, which showed low LD with neighbouring SNPs but considerable LD with a group of SNPs on ECA10. Furthermore, a cluster of SNPs on ECA5 showed unusually low LD with surrounding SNPs. A total of 39 SNPs met the criteria for unusual within-chromosome LD. The results of this study indicate that some SNPs may be misplaced. This finding is significant, as misplaced SNPs may lead to difficulties in the application of genomic methods, such as homozygosity mapping, for which SNP order is important.

Results of a haplotype-based GWAS for recurrent laryngeal neuropathy in the horse

Recurrent laryngeal neuropathy (RLN) is a major upper-airway disease of horses that causes abnormal respiratory noise during exercise and can impair performance. Etiopathogenesis remains unclear but genetic factors have been suspected for many decades. The objective of this study was to identify risk loci associated with RLN. To that end we genotyped 234 cases (196 Warmbloods, 20 Trotters, 14 Thoroughbreds, and 4 Draft horses), 228 breed-matched controls, and 69 parents with the Illumina Equine SNP50 BeadChip. Using these data, we quantified population structure and performed single-marker and haplotype-based association studies, as well as family-based linkage analyses. We accounted for population stratification by modeling a random polygenic background effect with covariance structure estimated from genome-wide SNP data. Using the haplotype-based approach, we identified two genome-wide suggestive loci in Warmbloods, respectively on chromosomes 21 (p = 1.62 × 10(-6)) and 31 (p = 1.69 × 10(-5)). The two signals were driven by the enrichment of a "protective" haplotype in controls compared to cases.