Thoroughbred stallion fertility is significantly associated with FKBP6 genotype but not with inbreeding or the contribution of a leading sire.

This is a follow-up study to validate the previously detected association of the FKBP6 gene with stallion subfertility. Using a select cohort of 150 Thoroughbred stallions with detailed breeding records, we confirm significant association (P < 0.0001) between low per-cycle pregnancy rates (≤50%) and a combined A/A-A/A genotype of SNPs chr13:11 353 372G>A and chr13:11 353 436A>C in FKBP6 exon 5. We also show that stallion subfertility and the combined genotype A/A-A/A are not associated with the level of genetic diversity based on 12 autosomal microsatellite markers, or with pedigree-based inbreeding rate, or the extent of contribution of a leading Thoroughbred sire, Northern Dancer, in a stallion's pedigree. We develop a TaqMan allelic discrimination assay for the two SNPs to facilitate accurate and high-throughput genotyping. We determine allele, genotype and combined genotype frequencies of FKBP6 exon 5 SNPs in a global cohort of 518 Thoroughbreds (76% stallions or geldings and 24% mares) and show that the frequency of the A/A-A/A genotype is 4%. Because there is no similar association between the FKBP6 exon 5 genotype and stallion subfertility in Hanoverians, we suggest that the two SNPs are not causative but rather tagging a breed-specific haplotype with genetic variants unique to Thoroughbreds. Further WGS-based research is needed to identify the molecular causes underlying the observed genotype-phenotype association in Thoroughbred stallions.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era.

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite-free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high-quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Generation of an equine biobank to be used for Functional Annotation of Animal Genomes project.

The Functional Annotation of Animal Genomes (FAANG) project aims to identify genomic regulatory elements in both sexes across multiple stages of development in domesticated animals. This study represents the first stage of the FAANG project for the horse, Equus caballus. A biobank of 80 tissue samples, two cell lines and six body fluids was created from two adult Thoroughbred mares. Ante-mortem assessments included full physical examinations, lameness, ophthalmologic and neurologic evaluations. Complete blood counts and serum biochemistries were also performed. At necropsy, in addition to tissue samples, aliquots of serum, ethylenediaminetetraacetic acid (EDTA) plasma, heparinized plasma, cerebrospinal fluid, synovial fluid, urine and microbiome samples from all regions of the gastrointestinal and urogenital tracts were collected. Epidermal keratinocytes and dermal fibroblasts were cultured from skin samples. All tissues were grossly and histologically evaluated by a board-certified veterinary pathologist. The results of the clinical and pathological evaluations identified subclinical eosinophilic and lymphocytic infiltration throughout the length of the gastrointestinal tract as well as a mild clinical lameness in both animals. Each sample was cryo-preserved in multiple ways, and nuclei were extracted from selected tissues. These samples represent the first published systemically healthy equine-specific biobank with extensive clinical phenotyping ante- and post-mortem. The tissues in the biobank are intended for community-wide use in the functional annotation of the equine genome. The use of the biobank will improve the quality of the reference annotation and allow all equine researchers to elucidate unknown genomic and epigenomic causes of disease.

Skeletal variation in Tennessee Walking Horses maps to the LCORL/NCAPG gene region.

Conformation has long been a driving force in horse selection and breed creation as a predictor for performance. The Tennessee Walking Horse (TWH) ranges in size from 1.5 to 1.7 m and is often used as a trail, show, and pleasure horse. To investigate the contribution of genetics to body conformation in the TWH, we collected DNA samples, body measurements, and gait/training information from 282 individuals. We analyzed the 32 body measures with a principal component analysis. Principal component (PC)1 captured 28.5% of the trait variance, while PC2 comprised just 9.5% and PC3 6.4% of trait variance. All 32 measures correlated positively with PC1, indicating that PC1 describes overall body size. We genotyped 109 horses using the EquineSNP70 bead chip and marker association assessed the data using PC1 scores as a phenotype. Mixed-model linear analysis (EMMAX) revealed a well-documented candidate locus on ECA3 (raw P = 3.86 × 10(-9)) near the LCORL gene. A custom genotyping panel enabled fine-mapping of the PC1 body-size trait to the 3'-end of the LCORL gene (P = 7.09 × 10(-10)). This position differs from other reports suggesting single nucleotide polymorphisms (SNPs) upstream of the LCORL coding sequence regulate expression of the gene and, therefore, body size in horses. Fluorescent in situ hybridization analysis defined the position of a highly homologous 5 kb retrogene copy of LCORL (assigned to unplaced contigs of the EquCab 2.0 assembly) at ECA9 q12-q13. This is the first study to identify putative causative SNPs within the LCORL transcript itself, which are associated with skeletal size variation in horses.

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves).

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves-an asthmalike inflammatory disease in horses. Analysis of 16 RAO-susceptible (cases) and six RAO-resistant (control) horses on a custom-made whole-genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO-affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.

Fertility Assessment in Sorraia Stallions by Sperm-Fish and Fkbp6 Genotyping.

The Sorraia, a critically endangered indigenous Iberian horse breed, is characterized by low genetic variability, high rate of inbreeding, bad sperm quality and subfertility. Here, we studied 11 phenotypically normal but subfertile Sorraia stallions by karyotyping, sex chromosome sperm-FISH and molecular analysis of FKBP6 - a susceptibility locus for impaired acrosome reaction (IAR). The stallions had normal sperm concentration (>300 million cells/ml), but the numbers of progressively motile sperm (21%) and morphologically normal sperm (28%) were invariably low. All stallions had a normal 64,XY karyotype. The majority of sperm (89%) had normal haploid sex chromosome content, although 11% of sperm carried various sex chromosome aneuploidies. No correlation was found between the percentage of sperm sex chromosome abnormalities and inbreeding, sperm morphology or stallion age. Direct sequencing of FKBP6 exon 4 for SNPs g.11040315G>A and g.11040379C>A revealed that none of the stallions had the susceptibility genotype (A/A-A/A) for IAR. Instead, all animals had a G/G-A/A genotype - a testimony of low genetic variability. The findings ruled out chromosomal abnormalities and genetic predisposition for IAR as contributing factors for subfertility. However, low fertility of the Sorraia stallions could be partly attributed to relatively higher rate of sex chromosome aneuploidies in the sperm.

A Non-Reciprocal Autosomal Translocation 64,XX, t(4;10)(q21;p15) in an Arabian Mare with Repeated Early Embryonic Loss.

Balanced autosomal translocations are a known cause for repeated early embryonic loss (REEL) in horses. In most cases, carriers of such translocations are phenotypically normal, but the chromosomal aberration negatively affects gametogenesis giving rise to both genetically balanced and unbalanced gametes. The latter, if involved in fertilization, result in REEL, whereas gametes with the balanced form of translocation will pass the defect into next generation. Therefore, in order to reduce the incidence of REEL, identification of translocation carriers is critical. Here, we report about a phenotypically normal 3-year-old Arabian mare that had repeated resorption of conceptuses prior to day 45 of gestation and was diagnosed with REEL. Conventional and molecular cytogenetic analyses revealed that the mare had normal chromosome number 64,XX but carried a non-mosaic and non-reciprocal autosomal translocation t(4;10)(q21;p15). This is a novel translocation described in horses with REEL and the first such report in Arabians. Previous cases of REEL due to autosomal translocations have exclusively involved Thoroughbreds. The findings underscore the importance of routine cytogenetic screening of breeding animals.

Comparative organization and gene expression profiles of the porcine pseudoautosomal region.

The pseudoautosomal region (PAR) has important biological functions in spermatogenesis, male fertility and early development. Even though pig (Sus scrofa, SSC) is an agriculturally and biomedically important species, and its genome is sequenced, current knowledge about the porcine PAR is sparse. Here we defined the PAR in SSCXp/Yp by demarcating the sequence of the pseudoautosomal boundary at X:6,743,567 bp in intron 3-4 of SHROOM2 and showed that SHROOM2 is truncated in SSCY. Cytogenetic mapping of 20 BAC clones containing 15 PAR and X-specific genes revealed that the pig PAR is largely collinear with other mammalian PARs or Xp terminal regions. The results improved the current SSCX sequence assembly and facilitated distinction between the PAR and X-specific genes to study their expression in adult and embryonic tissues. A pilot analysis showed that the PAR genes are expressed at higher levels than X-specific genes during early development, whereas the expression of PAR genes was higher at day 60 compared to day 26, and higher in embryonic tissues compared to placenta. The findings advance the knowledge about the comparative organization of the PAR in mammals and suggest that the region might have important functions in early development in pigs.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA.

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Molecular heterogeneity of XY sex reversal in horses.

Male-to-female 64,XY sex reversal is a frequently reported chromosome abnormality in horses. Despite this, the molecular causes of the condition are as yet poorly understood. This is partially because only limited molecular information is available for the horse Y chromosome (ECAY). Here, we used the recently developed ECAY map and carried out the first comprehensive study of the Y chromosome in XY mares (n=18). The integrity of the ECAY in XY females was studied by FISH and PCR using markers evenly distributed along the euchromatic region. The results showed that the XY sex reversal condition in horses has two molecularly distinct forms: (i) a Y-linked form that is characterized by Y chromosome deletions and (ii) a non-Y-linked form where the Y chromosome of affected females is molecularly the same as in normal males. Further analysis of the Y-linked form (13 cases) showed that the condition is molecularly heterogeneous: the smallest deletions spanned about 21 kb, while the largest involved the entire euchromatic region. Regardless of the size, all deletions included the SRY gene. We show that the deletions were likely caused by inter-chromatid recombination events between repeated sequences in ECAY. Further, we hypothesize that the occurrence of SRY-negative XY females in some species (horse, human) but not in others (pig, dog) is because of differences in the organization of the Y chromosome. Finally, in contrast to the Y-linked SRY-negative form of equine XY sex reversal, the molecular causes of SRY-positive XY mares (5 cases) remain as yet undefined.

Characterization of the bovine pseudoautosomal region and comparison with sheep, goat, and other mammalian pseudoautosomal regions.

The pseudoautosomal region (PAR) is a small region of sequence homology between mammalian X and Y chromosomes and is needed for sex chromosome segregation in male meiosis. The region, though studied as yet in only a few species, shows considerable variation in size and gene content. We have constructed a medium-density gene map for the cattle PAR and the adjacent X-specific region by isolating and mapping 18 BAC clones which contain 20 PAR- and 5 X-specific genes. One BAC clone containing TBL1XY and GPR143 spanned the recently demarcated bovine pseudoautosomal boundary (PAB). Comparing the gene map of cattle PAR with the high-resolution maps of human, horse, and dog PAR allowed to estimate that the size of cattle PAR is approximately 5-9 Mb. BAC end sequence analysis showed that there is a gradient of decreasing GC content from PARter towards the PAB which is consistent with findings in human, mouse, and horse. The 20 PAR- and 5 X-specific cattle genes were mapped also in goat and sheep, showing that PAR in the 3 species is similar in size, gene content, and gene order. For the first time the PAB was determined in goat sex chromosomes. Comparison of cattle, goat, and sheep PAR with homologous regions on human and horse X chromosomes showed a high degree of linkage conservation between all species. However, the most terminal human, horse, and dog PAR gene, PLCXD1, is X-specific in ruminants. Since the human/horse linkage group containing PLCXD1 is of ancestral origin, the location of PLCXD1 can be considered as a de novo event in ruminant sex chromosome evolution. The gene map of the cattle PAR adds to our knowledge about the comparative organization and evolution of the eutherian PAR and aids the sequencing, sequence assembly, and annotation of the terminal region of BTAXq.

The 450-band resolution G- and R-banded standard karyotype of the donkey (Equus asinus, 2n = 62).

Donkey chromosomes were earlier characterized separately by C-, G- and R-banding techniques. However, direct comparisons between G- and R-banding patterns have still not been carried out in this species. The present study reports this comparison at the 450-band level by using replication G- and R-banding patterns. Two sets of synchronized lymphocyte cultures were set up to obtain early (GBA+CBA-banding) and late (RBA-banding) BrdU incorporation. Slides were stained with acridine orange and observed under a fluorescence microscope. Reverse GBA+CBA- and RBA-banded karyotypes at the 450-band level were constructed. To verify G- and R-banding patterns in some acrocentric chromosomes, sequential GBA+CBA/Ag-NORs and RBA/Ag-NORs were also performed. The results of CBA-banding patterns obtained in 12 animals from 2 breeds showed a pronounced polymorphism of heterochromatin, especially in EAS1q-prox. Ideogrammatic representations of G- and R-banded karyotypes were constructed using only one common G- and R-banding nomenclature. In the present study both G- and R-banding patterns and relative ideograms are presented as standard karyotype for this species at the 450-band level.

The horse pseudoautosomal region (PAR): characterization and comparison with the human, chimp and mouse PARs.

The pseudoautosomal region (PAR) is a genomic segment on mammalian sex chromosomes where sequence homology mimics that seen between autosomal homologues. The region is essential for pairing and proper segregation of sex chromosomes during male meiosis. As yet, only human/chimp and mouse PARs have been characterized. The two groups of species differ dramatically in gene content and size of the PAR and therefore do not provide clues about the likely evolution and constitution of PAR among mammals. Here we characterize the equine PAR by i) isolating and arranging 71 BACs containing 129 markers (110 STS and 19 genes) into two contigs spanning the region, ii) precisely localizing the pseudoautosomal boundary (PAB), and iii) describing part of the contiguous X- and Y-specific regions. We also report the discovery of an approximately 200 kb region in the middle of the PAR that is present in the male-specific region of the Y (MSY) as well. Such duplication is a novel observation in mammals. Further, comparison of the equine PAR with the human counterpart shows that despite containing orthologs from an additional 1 Mb region beyond the human PAR1, the equine PAR is around 0.9 Mb smaller than the size of the human PAR. We theorize that the PAR varies in size and gene content across evolutionarily closely as well as distantly related mammals. Although striking differences like those observed between human and mouse may be rare, variations similar to those seen between horse and human may be prevalent among mammals.

A BAC contig map over the proximal approximately 3.3 Mb region of horse chromosome 21.

A total of 207 BAC clones containing 155 loci were isolated and arranged into a map of linearly ordered overlapping clones over the proximal part of horse chromosome 21 (ECA21), which corresponds to the proximal half of the short arm of human chromosome 19 (HSA19p) and part of HSA5. The clones form two contigs - each corresponding to the respective human chromosomes - that are estimated to be separated by a gap of approximately 200 kb. Of the 155 markers present in the two contigs, 141 (33 genes and 108 STS) were generated and mapped in this study. The BACs provide a 4-5x coverage of the region and span an estimated length of approximately 3.3 Mb. The region presently contains one mapped marker per 22 kb on average, which represents a major improvement over the previous resolution of one marker per 380 kb obtained through the generation of a dense RH map for this segment. Dual color fluorescence in situ hybridization on metaphase and interphase chromosomes verified the relative order of some of the BACs and helped to orient them accurately in the contigs. Despite having similar gene order and content, the equine region covered by the contigs appears to be distinctly smaller than the corresponding region in human (3.3 Mb vs. 5.5-6 Mb) because the latter harbors a host of repetitive elements and gene families unique to humans/primates. Considering limited representation of the region in the latest version of the horse whole genome sequence EquCab2, the dense map developed in this study will prove useful for the assembly and annotation of the sequence data on ECA21 and will be instrumental in rapid search and isolation of candidate genes for traits mapped to this region.

A 4,103 marker integrated physical and comparative map of the horse genome.

A comprehensive second-generation whole genome radiation hybrid (RH II), cytogenetic and comparative map of the horse genome (2n = 64) has been developed using the 5000rad horse x hamster radiation hybrid panel and fluorescence in situ hybridization (FISH). The map contains 4,103 markers (3,816 RH; 1,144 FISH) assigned to all 31 pairs of autosomes and the X chromosome. The RH maps of individual chromosomes are anchored and oriented using 857 cytogenetic markers. The overall resolution of the map is one marker per 775 kilobase pairs (kb), which represents a more than five-fold improvement over the first-generation map. The RH II incorporates 920 markers shared jointly with the two recently reported meiotic maps. Consequently the two maps were aligned with the RH II maps of individual autosomes and the X chromosome. Additionally, a comparative map of the horse genome was generated by connecting 1,904 loci on the horse map with genome sequences available for eight diverse vertebrates to highlight regions of evolutionarily conserved syntenies, linkages, and chromosomal breakpoints. The integrated map thus obtained presents the most comprehensive information on the physical and comparative organization of the equine genome and will assist future assemblies of whole genome BAC fingerprint maps and the genome sequence. It will also serve as a tool to identify genes governing health, disease and performance traits in horses and assist us in understanding the evolution of the equine genome in relation to other species.

Assignment of six genes to bovine chromosomes 5 and 16 by fluorescence in situ hybridization, radiation hybrid mapping and genetic linkage analysis.

Several quantitative trait loci for beef carcass traits have been mapped to bovine chromosome 5. The objective of this study was to map six candidate genes for these traits by fluoresence in situ hybridization, genetic linkage analysis and radiation hybrid mapping. MYF5 and MYF6 were assigned to 5q13, WIF1 to 5q23 and MMP19 to 5q25. A paralog of MYF5 (putatively MYOG) was assigned to 16q12. A novel microsatellite placed MYF5 and MYF6 10.4 cM from BM6026 and 19.1 cM from BL23 on the genetic linkage map. MYF5 (62.6 cR), WNT10B (319.5 cR), WIF1 (500.8 cR) and MMP19 (701.2 cR) were also integrated into the 5000(Rad) radiation hybrid map.

A 1.3-Mb interval map of equine homologs of HSA2.

A comparative approach that utilizes information from more densely mapped or sequenced genomes is a proven and efficient means to increase our knowledge of the structure of the horse genome. Human chromosome 2 (HSA2), the second largest human chromosome, comprising 243 Mb, and containing 1246 known genes, corresponds to all or parts of three equine chromosomes. This report describes the assignment of 140 new markers (78 genes and 62 microsatellites) to the equine radiation hybrid (RH) map, and the anchoring of 24 of these markers to horse chromosomes by FISH. The updated equine RH maps for ECA6p, ECA15, and ECA18 resulting from this work have one, two, and three RH linkage groups, respectively, per chromosome/chromosome-arm. These maps have a three-fold increase in the number of mapped markers compared to previous maps of these chromosomes, and an increase in the average marker density to one marker per 1.3 Mb. Comparative maps of ECA6p, ECA15, and ECA18 with human, chimpanzee, dog, mouse, rat, and chicken genomes reveal blocks of conserved synteny across mammals and vertebrates.

Linkage and comparative mapping of the locus controlling susceptibility towards E. COLI F4ab/ac diarrhoea in pigs.

In 1995, Edfors-Lilja and coworkers mapped the locus for the E. COLI K88ab (F4ab) and K88ac (F4ac) intestinal receptor to pig chromosome 13 (SSC13). Using the same family material we have refined the map position to a region between the microsatellite markers Sw207 and Sw225. Primers from these markers were used to screen a pig BAC library and the positive clones were used for fluorescent in situ hybridization (FISH) analysis. The results of the FISH analysis helped to propose a candidate gene region in the SSC13q41-->q44 interval. Shotgun sequencing of the FISH-mapped BAC clones revealed that the candidate region contains an evolutionary breakpoint between human and pig. In order to further characterise the rearrangements between SSC13 and human chromosome 3 (HSA3), detailed gene mapping of SSC13 was carried out. Based on this mapping data we have constructed a detailed comparative map between SSC13 and HSA3. Two candidate regions on human chromosome 3 have been identified that are likely to harbour the human homologue of the gene responsible for susceptibility towards E. COLI F4ab/ac diarrhoea in pigs.

An ordered BAC contig map of the equine major histocompatibility complex.

A physical map of ordered bacterial artificial chromosome (BAC) clones was constructed to determine the genetic organization of the horse major histocompatibility complex. Human, cattle, pig, mouse, and rat MHC gene sequences were compared to identify highly conserved regions which served as source templates for the design of overgo primers. Thirty-five overgo probes were designed from 24 genes and used for hybridization screening of the equine USDA CHORI 241 BAC library. Two hundred thirty-eight BAC clones were assembled into two contigs spanning the horse MHC region. The first contig contains the MHC class II region and was reduced to a minimum tiling path of nine BAC clones that span approximately 800 kb and contain at least 20 genes. A minimum tiling path of a second contig containing the class III/I region is comprised of 14 BAC clones that span approximately 1.6 Mb and contain at least 34 genes. Fluorescence in situ hybridization (FISH) using representative clones from each of the three regions of the MHC localized the contigs onto ECA20q21 and oriented the regions relative to one another and the centromere. Dual-colored FISH revealed that the class I region is proximal to the centromere, the class II region is distal, and the class III region is located between class I and II. These data indicate that the equine MHC is a single gene-dense region similar in structure and organization to the human MHC and is not disrupted as in ruminants and pigs.

Genetic mapping of GBE1 and its association with glycogen storage disease IV in American Quarter horses.

Comparative biochemical and histopathological data suggest that a deficiency in the glycogen branching enzyme (GBE) is responsible for a fatal neonatal disease in Quarter Horse foals that closely resembles human glycogen storage disease type IV (GSD IV). Identification of DNA markers closely linked to the equine GBE1 gene would assist us in determining whether a mutation in this gene leads to the GSD IV-like condition. FISH using BAC clones as probes assigned the equine GBE1 gene to a marker deficient region of ECA26q12-->q13. Four other genes, ROBO2, ROBO1, POU1F1, and HTR1F, that flank GBE1 within a 10-Mb segment of HSA3p12-->p11, were tightly linked to equine GBE1 when analyzed on the Texas A&M University 5000 rad equine radiation hybrid panel, while the GLB1, MITF, RYBP, and PROS1 genes that flank this 10-Mb interval were not linked with markers in the GBE1 group. A polymorphic microsatellite (GBEms1) in a GBE1 BAC clone was then identified and genetically mapped to ECA26 on the Animal Health Trust full-sibling equine reference family. All Quarter Horse foals affected with GSD IV were homozygous for an allele of GBEms1, as well as an allele of the most closely linked microsatellite marker, while a control horse population showed significant allelic variation with these markers. This data provides strong molecular genetic support for the candidacy of the GBE1 locus in equine GSD IV.