Heritability of Recurrent Exertional Rhabdomyolysis in Standardbred and Thoroughbred Racehorses Derived From SNP Genotyping Data.

Recurrent exertional rhabdomyolysis (RER) in Thoroughbred and Standardbred racehorses is characterized by episodes of muscle rigidity and cell damage that often recur upon strenuous exercise. The objective was to evaluate the importance of genetic factors in RER by obtaining an unbiased estimate of heritability in cohorts of unrelated Thoroughbred and Standardbred racehorses. Four hundred ninety-one Thoroughbred and 196 Standardbred racehorses were genotyped with the 54K or 74K SNP genotyping arrays. Heritability was calculated from genome-wide SNP data with a mixed linear and Bayesian model, utilizing the standard genetic relationship matrix (GRM). Both the mixed linear and Bayesian models estimated heritability of RER in Thoroughbreds to be approximately 0.34 and in Standardbred racehorses to be approximately 0.45 after adjusting for disease prevalence and sex. To account for potential differences in the genetic architecture of the underlying causal variants, heritability estimates were adjusted based on linkage disequilibrium weighted kinship matrix, minor allele frequency and variant effect size, yielding heritability estimates that ranged between 0.41-0.46 (Thoroughbreds) and 0.39-0.49 (Standardbreds). In conclusion, between 34-46% and 39-49% of the variance in RER susceptibility in Thoroughbred and Standardbred racehorses, respectively, can be explained by the SNPs present on these 2 genotyping arrays, indicating that RER is moderately heritable. These data provide further rationale for the investigation of genetic mutations associated with RER susceptibility.

Accumulating mutations in series of haplotypes at the KIT and MITF loci are major determinants of white markings in Franches-Montagnes horses.

Coat color and pattern variations in domestic animals are frequently inherited as simple monogenic traits, but a number are known to have a complex genetic basis. While the analysis of complex trait data remains a challenge in all species, we can use the reduced haplotypic diversity in domestic animal populations to gain insight into the genomic interactions underlying complex phenotypes. White face and leg markings are examples of complex traits in horses where little is known of the underlying genetics. In this study, Franches-Montagnes (FM) horses were scored for the occurrence of white facial and leg markings using a standardized scoring system. A genome-wide association study (GWAS) was performed for several white patterning traits in 1,077 FM horses. Seven quantitative trait loci (QTL) affecting the white marking score with p-values p≤10(-4) were identified. Three loci, MC1R and the known white spotting genes, KIT and MITF, were identified as the major loci underlying the extent of white patterning in this breed. Together, the seven loci explain 54% of the genetic variance in total white marking score, while MITF and KIT alone account for 26%. Although MITF and KIT are the major loci controlling white patterning, their influence varies according to the basic coat color of the horse and the specific body location of the white patterning. Fine mapping across the MITF and KIT loci was used to characterize haplotypes present. Phylogenetic relationships among haplotypes were calculated to assess their selective and evolutionary influences on the extent of white patterning. This novel approach shows that KIT and MITF act in an additive manner and that accumulating mutations at these loci progressively increase the extent of white markings.

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.

A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies.

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50-100 kb and reached background levels within 1-2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species.

Accurate determination of phenotypic information from historic thoroughbred horses by single base extension.

Historic DNA data have the potential to identify phenotypic information otherwise invisible in the historical, archaeological and palaeontological record. In order to determine whether a single nucleotide polymorphism typing protocol based on single based extension (SNaPshot™) could produce reliable phenotypic data from historic samples, we genotyped three coat colour markers for a sample of historic Thoroughbred horses for which both phenotypic and correct genotypic information were known from pedigree information in the General Stud Book. Experimental results were consistent with the pedigrees in all cases. Thus we demonstrate that historic DNA techniques can produce reliable phenotypic information from museum specimens.

L-2-hydroxyglutaric aciduria: characterisation of the molecular defect in a spontaneous canine model.

l-2-hydroxyglutaric aciduria (l-2-HGA) is a neurometabolic disorder that produces a variety of clinical neurological deficits, including psychomotor retardation, seizures and ataxia. The biochemical hallmark of l-2-HGA is the accumulation of l-2-hydroxyglutaric acid (l-2-HG) in cerebrospinal fluid, plasma and urine. Mutations within the gene L2HGDH (Entrez Gene ID 79944) on chromosome 14q22 encoding L-2-hydroxyglutaric acid dehydrogenase have recently been shown to cause l-2-HGA in humans. Using a candidate gene approach in an outbred pet dog population segregating l-2-HGA, the causal molecular defect was identified in the canine homologue of L2HGDH and characterised. DNA sequencing and pedigree analysis indicate a common founder effect in the canine model. The canine model shares many of the clinical and MRI features of the disease in humans and represents a valuable resource as a spontaneous model of l-2-HGA.

Exclusion of linkage of the RYR1, CACNA1S, and ATP2A1 genes to recurrent exertional rhabdomyolysis in Thoroughbreds.

OBJECTIVE: To determine whether there was genetic linkage between the recurrent exertional rhabdomyolysis (RER) trait in Thoroughbred horse pedigrees and DNA markers in genes (the sarcoplasmic reticulum calcium release channel [RYR1] gene, the sarcoplasmic reticulum calcium ATPase [ATP2A1] gene, and the transverse tubule dihydropyridine receptor-voltage sensor [CACNA1S] gene) that are important in myoplasmic calcium regulation. ANIMALS: 34 horses in the University of Minnesota RER resource herd and 62 Thoroughbreds from 3 families of Thoroughbreds outside of the university in which RER-affected status was assigned after 2 or more episodes of ER had been observed. PROCEDURES: Microsatellite DNA markers from the RYR1, ATP2A1, and CACNA1S gene loci on equine chromosomes 10, 13, and 30 were identified. Genotypes were obtained for all horses in the 4 families affected by RER, and data were used to test for linkage of these 3 loci to the RER phenotype. RESULTS: Analysis of the RYR1, CACNA1S, and ATP2A1 microsatellites excluded a link between those markers and the RER trait. CONCLUSIONS AND CLINICAL RELEVANCE: It is likely that the heritable alterations in muscle contractility that are characteristic of RER are caused by a gene that is not yet known to cause related muscle disease in other species.

A human-horse comparative map based on equine BAC end sequences.

In an effort to increase the density of sequence-based markers for the horse genome we generated 9473 BAC end sequences (BESs) from the CHORI-241 BAC library with an average read length of 677 bp. BLASTN searches with the BESs revealed 4036 meaningful hits (E

Genomic sequence of the class II region of the canine MHC: comparison with the MHC of other mammalian species.

The domestic dog, Canis familiaris, is an excellent model species in which to study complex inherited diseases, having over 200 recognized breeds, each of which represents a closed gene pool. Overlapping canine genomic BAC clones were sequenced to obtain 711,521 bp of the canine classical and extended MHC class II regions. Analysis and annotation of this sequence reveals that it contains 45 loci, of which 29 are predicted to be functionally expressed. Comparison of the DLA class II sequence with those of the cat, human, and mouse highlights regions of syntenic conservation and species-specific gene rearrangement and duplication and gives an insight into the evolution of the DR region in the order Carnivora. Elucidation of functionally important dog class II genes and the identification of 23 microsatellite markers spanning this region will contribute significantly to the study of canine diseases that have an immune component.

Glycogen branching enzyme (GBE1) mutation causing equine glycogen storage disease IV.

Comparative biochemical and histopathological evidence suggests that a deficiency in the glycogen branching enzyme, encoded by the GBE1 gene, is responsible for a recently identified recessive fatal fetal and neonatal glycogen storage disease (GSD) in American Quarter Horses termed GSD IV. We have now derived the complete GBE1 cDNA sequences for control horses and affected foals, and identified a C to A substitution at base 102 that results in a tyrosine (Y) to stop (X) mutation in codon 34 of exon 1. All 11 affected foals were homozygous for the X34 allele, their 11 available dams and sires were heterozygous, and all 16 control horses were homozygous for the Y34 allele. The previous findings of poorly branched glycogen, abnormal polysaccharide accumulation, lack of measurable GBE1 enzyme activity and immunodetectable GBE1 protein, coupled with the present observation of abundant GBE1 mRNA in affected foals, are all consistent with the nonsense mutation in the 699 amino acid GBE1 protein. The affected foal pedigrees have a common ancestor and contain prolific stallions that are likely carriers of the recessive X34 allele. Defining the molecular basis of equine GSD IV will allow for accurate DNA testing and the ability to prevent occurrence of this devastating disease affecting American Quarter Horses and related breeds.

Construction and integration of radiation-hybrid and cytogenetic maps of dog Chromosome X.

Chromosome (chr) X is under-represented in current maps of the genome of the domestic dog ( Canis familiaris). To address this problem, we have constructed a small-insert, genomic DNA library in pBluescript from flow-sorted canine Chr X DNA. Fluorescence in situ hybridization (FISH) studies confirmed that the library was highly enriched for Chr X. Clones containing microsatellites were identified and sequenced. Database searches detected significant sequence identity between four X-derived clones and genes previously characterized in other species. Thirty-seven markers derived from these clones were mapped on Chr X by FISH, and of these, 28 were mapped by using the female-derived T72 whole-genome radiation hybrid (RH) panel (Research Genetics). Four X-linked canine genes from publicly available data were also mapped. Eight RH linkage groups with LOD >4.0 were identified, and FISH data were used to locate the groups on the chromosome; four groups could be unambiguously orientated by FISH data. In each case, the FISH and RH data were mutually consistent. The data suggest strongly conserved synteny between canine and human X Chrs. The pseudoautosomal region has been further characterized, and the putative or actual locations of nine genes of clinical relevance have been suggested.