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.     

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.     

Molecular phylogeny of Indian horse breeds with special reference to Manipuri pony based on mitochondrial D-loop

Manipuri pony is the geographically distant breed of horse from the five recognized horse breeds found in the Indian subcontinent. The phylogenetic relationship of Manipuri pony with the other breeds is unknown. The diversity in the mitochondrial (mt) DNA D-loop region is employed as an important tool to understand the origin and genetic diversification of domestic horses and to examine genetic relationships among breeds around the world. This study was carried out to understand the maternal lineages of Manipuri pony using the 247 bp region of the mtDNA D-loop. The dataset comprised of eleven numbers of self developed sequences of Manipuri pony, 59 and 35 number of retrieved sequences of Indian horse breeds and other worldwide breeds respectively. A total of 35 haplotypes was identified with a high level of genetic diversity in the Indian breeds. A total of seven major mtDNA haplogroups (A–G) was identified in the Indian horse breeds that indicated the abundance of mtDNA diversity and multiple origins of maternal lineages in them. The majority of the studied sequences of Indian breeds (33.3 %) were grouped into haplogroup D and least (3.9 %) in haplogroup E. The Manipuri breed showed the least F_ST distance (0.03866) with the most diverged Indian breeds and with Thoroughbred horse among the worldwide. This study indicated a close association between Manipuri pony and Thoroughbred.     

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.     

A genome-wide scan for tying-up syndrome in Japanese Thoroughbreds

Tying-up syndrome, also known as recurrent exertional rhabdomyolysis in Thoroughbreds, is a common muscle disorder for racehorses. In this study, we performed a multipoint linkage analysis using LOKI based on the Bayesian Markov chain Monte Carlo method using 5 half-sib families (51 affected and 277 nonaffected horses in total), and a genome-wide association study (GWAS) using microsatellites (144 affected and 144 nonaffected horses) to map candidate regions for tying-up syndrome in Japanese Thoroughbreds. The linkage analysis identified one strong L-score (82.45) between the loci UCDEQ411 and COR058 (24.9-27.9 Mb) on ECA12. The GWAS identified two suggestive genomic regions on ECA12 (24.9-27.8 Mb) and ECA20 (29.3-33.5 Mb). Based on both results, the genomic region between UCDEQ411 and TKY499 (24.9-27.8 Mb) on ECA12 was the most significant and was considered as a candidate region for tying-up syndrome in Japanese Thoroughbreds.     

Replication and fine-mapping of a QTL for recurrent airway obstruction in European Warmblood horses

Recurrent airway obstruction (RAO), or ‘heaves’, is a common performance‐limiting allergic respiratory disease of mature horses. It is related to sensitization and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. In a previous study, we detected a QTL for RAO on ECA 13 in a half‐sib family of European Warmblood horses. In this study, we genotyped additional markers in the family and narrowed the QTL down to about 1.5 Mb (23.7–25.2 Mb). We detected the strongest association with SNP BIEC2‐224511 (24 309 405 bp). We also obtained SNP genotypes in an independent cohort of 646 unrelated Warmblood horses. There was no genome‐wide significant association with RAO in these unrelated horses. However, we performed a genotypic association study of the SNPs on ECA 13 in these unrelated horses, and the SNP BIEC2‐224511 also showed the strongest association with RAO in the unrelated horses (p_raw = 0.00037). The T allele at this SNP was associated with RAO both in the family and the unrelated horses. Thus, the association study in the unrelated animals provides independent support for the previously detected QTL. The association study allows further narrowing of the QTL interval to about 0.5 Mb (24.0–24.5 Mb). We sequenced the coding regions of the genes in the critical region but did not find any associated coding variants. Therefore, the causative variant underlying this QTL is likely to be a regulatory mutation.     

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.     

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.     

猪a1-岩藻糖转移酶基因(FUT1)在26个中外猪种中的遗传变异研究

肠毒素大肠杆菌F18(ECF18)是引起仔猪断奶后水肿和腹泻病的主要病原菌,a1—岩藻糖转移酶基因(FUT1)是ECF18侵染猪小肠的受体蛋白候选基因。通过采用PCR—RFLP方法检测了5个西方商业猪种以及21个中国地方猪种(群)1458个个体在FUT1基因开放阅读框架的307核苷酸位点的G-A点突变(M307^G-A)遗传变异。结果表明：5个外来猪种以及中国地方猪种中的临高猪在该FUT1基因位点存在多态性,其他中国地方猪种均表现为极端的单态分布,只有易感的GG基因型,没有多态性。由此提示：1)如果猪FUT1 M307^G-A点突变是决定猪小肠上ECF18受体表达与否的关键因素,则绝大部分中国地方猪种均不具备抵抗ECF18的遗传基础,这除了表明ECF18抗性基因有可能起源于西方猪种外,同时也表明对中国地方猪种中在这个位点惟一存在多态性的海南临高猪的品种资源保存具有非常重要的意义。2)一般而言,在中国的养猪生产实践中,中国地方猪种的仔猪抗水肿与腹泻病能力普遍强于外来猪种,研究的结果提示有必要对中国地方猪种所具备的上述遗传抗性做更深入的研究,寻找、定位其相应的QTL或／和抗性基因。     

Canine cystinuria: polymorphism in the canine<Emphasis Type="Italic"> SLC3A1</Emphasis> gene and identification of a nonsense mutation in cystinuric Newfoundland dogs

Cystinuria is an inherited renal and intestinal disease characterized by defective amino acid reabsorption and cystine urolithiasis. Different forms of the disease, designated type I and non-type I in cystinuric humans, can be distinguished clinically and biochemically, and have been associated with mutations in the SLC3A1 (rBAT) and SLC7A9 genes, respectively. Type I cystinuria is the most common form and is inherited as an autosomal recessive trait in humans. Cystinuria has been recognized in more than 60 breeds of dogs and a severe form, resembling type I cystinuria, has been characterized in the Newfoundland breed. Here we report the cloning and sequencing of the canine SLC3A1 cDNA and gene, and the identification of a nonsense mutation in exon 2 of the gene in cystinuric Newfoundland dogs. A mutation-specific test was developed for the diagnosis and control of cystinuria in Newfoundland dogs. In cystinuric dogs of six other breeds, either heterozygosity at the SLC3A1 locus or lack of mutations in the coding region of the SLC3A1 gene were observed, indicating that cystinuria is genetically heterogeneous in dogs, as it is in humans. The canine homologue of human type I cystinuria provides the opportunity to use a large animal model to investigate molecular approaches for the treatment of cystinuria and other renal tubular diseases.     

Genetic variability in the Skyros pony and its relationship with other Greek and foreign horse breeds

In Greece, seven native horse breeds have been identified so far. Among these, the Skyros pony is outstanding through having a distinct phenotype. In the present study, the aim was to assess genetic diversity in this breed, by using different types of genetic loci and available genealogical information. Its relationships with the other Greek, as well as foreign, domestic breeds were also investigated. Through microsatellite and pedigree analysis it appeared that the Skyros presented a similar level of genetic diversity to the other European breeds. Nevertheless, comparisons between DNA-based and pedigree-based results revealed that a loss of genetic diversity had probably already occurred before the beginning of breed registration. Tests indicated the possible existence of a recent bottleneck in two of the three main herds of Skyros pony. Nonetheless, relatively high levels of heterozygosity and Polymorphism Information Content indicated sufficient residual genetic variability, probably useful in planning future strategies for breed conservation. Three other Greek breeds were also analyzed. A comparison of these with domestic breeds elsewhere, revealed the closest relationships to be with the Middle Eastern types, whereas the Skyros itself remained isolated, without any close relationship, whatsoever.     

An imputation-based genome-wide association study for growth and fatness traits in Sujiang pigs

Sujiang pigs are a synthetic breed derived from Jiangquhai, Fengjing, and Duroc pigs. In this study, we sequenced the genome of 62 pigs with a coverage depth of 10× to 20×, including 27 Sujiang and 35 founder breed pigs, and we collected 360 global pigs’ genome sequence data from public databases including 39 Duroc pigs. We obtained a high-quality variant dataset of 365 Sujiang pigs by imputing the porcine 80 K single nucleotide polymorphism (SNP) Beadchip to the whole-genome scale with a total of 422 pigs as a reference panel. A dataset of 365 imputated Sujiang pigs was used to perform single-trait genome-wide association study (GWAS) and meta-analyses for growth and fatness traits. Single-trait GWAS identified 1 907, 18, and 14 SNPs surpassing the suggestively significant threshold for backfat thickness, chest circumference, and chest width, respectively. Meta-analyses identified 2 400 genome-wide significant SNPs and 520 suggestively significant SNPs for backfat thickness and chest circumference, and 719 genome-wide significant SNPs and 1 225 suggestively significant SNPs for all seven traits. According to the meta-analysis of backfat thickness and chest circumference, a remarkable region of 2.69 Mb on Sus scrofa chromosome 4 containing FAM110B, IMPAD1, LYN, MOS, PENK, PLAG1, SDR16C5 and XKR4 was identified as a candidate region. The haplotype heat map of the 2.69 Mb region verified that Sujiang pigs were derived from Duroc and Chinese indigenous pigs, especially Jiangquhai pigs. The Kruskal-Wallis test showed that haplotypes of the 2.69 Mb region significantly affected backfat thickness and chest circumference traits. We then focused on PLAG1, an important growth-related gene, and identified two synonymous SNPs with obvious differences among different breeds in the PLAG1 gene. We then performed genotyping of 365 Sujiang, 150 Duroc, 95 Jiangquhai, and 100 Fengjing pigs to confirm the above result and verified that the two variants significantly affected phenotypes of growth and fatness traits. Our findings not only provide insights into the genetic architecture of porcine growth and fatness traits but also provide potential markers for selective breeding of these traits in Sujiang pigs.     

Genetic mapping of recurrent exertional rhabdomyolysis in a population of North American Thoroughbreds

Recurrent exertional rhabdomyolysis is a heritable disorder that results in painful skeletal muscle cramping with exercise in up to 10% of all Thoroughbred racehorses. Here, we report a genome‐wide association study with 48 282 SNPs analyzed among 48 case and 37 control Thoroughbreds. The most significant SNPs spanned approximately 13 Mb on ECA16, and the P‐value of the most significant SNP after correcting for population structure was 8.0 × 10^?6. This region on ECA16 was further evaluated by genotyping 247 SNPs in both the initial population and a second population of 34 case and 98 control Thoroughbreds. Several SNPs across the 13‐Mb region on ECA16 showed significance when each population was analyzed separately; however, the exact positions of the most significant SNPs within this region on ECA16 varied between populations. This variability in location may be attributed to lack of power owing to insufficient sample sizes within each population individually, or to the relative distribution of long, conserved haplotypes, characteristic of the Thoroughbred breed. Future genome‐wide association studies with additional horses would likely improve the power to resolve casual loci located on ECA16 and increase the likelihood of detecting any additional loci on other chromosomes contributing to disease susceptibility.     

Localization of Canine Brachycephaly Using an Across Breed Mapping Approach

The domestic dog, Canis familiaris, exhibits profound phenotypic diversity and is an ideal model organism for the genetic dissection of simple and complex traits. However, some of the most interesting phenotypes are fixed in particular breeds and are therefore less tractable to genetic analysis using classical segregation-based mapping approaches. We implemented an across breed mapping approach using a moderately dense SNP array, a low number of animals and breeds carefully selected for the phenotypes of interest to identify genetic variants responsible for breed-defining characteristics. Using a modest number of affected (10–30) and control (20–60) samples from multiple breeds, the correct chromosomal assignment was identified in a proof of concept experiment using three previously defined loci; hyperuricosuria, white spotting and chondrodysplasia. Genome-wide association was performed in a similar manner for one of the most striking morphological traits in dogs: brachycephalic head type. Although candidate gene approaches based on comparable phenotypes in mice and humans have been utilized for this trait, the causative gene has remained elusive using this method. Samples from nine affected breeds and thirteen control breeds identified strong genome-wide associations for brachycephalic head type on Cfa 1. Two independent datasets identified the same genomic region. Levels of relative heterozygosity in the associated region indicate that it has been subjected to a selective sweep, consistent with it being a breed defining morphological characteristic. Genotyping additional dogs in the region confirmed the association. To date, the genetic structure of dog breeds has primarily been exploited for genome wide association for segregating traits. These results demonstrate that non-segregating traits under strong selection are equally tractable to genetic analysis using small sample numbers.     

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.     

Selection signatures in Shetland ponies

Shetland ponies were selected for numerous traits including small stature, strength, hardiness and longevity. Despite the different selection criteria, Shetland ponies are well known for their small stature. We performed a selection signature analysis including genome‐wide SNPs of 75 Shetland ponies and 76 large‐sized horses. Based upon this dataset, we identified a selection signature on equine chromosome (ECA) 1 between 103.8 Mb and 108.5 Mb. A total of 33 annotated genes are located within this interval including the IGF1R gene at 104.2 Mb and the ADAMTS17 gene at 105.4 Mb. These two genes are well known to have a major impact on body height in numerous species including humans. Homozygosity mapping in the Shetland ponies identified a region with increased homozygosity between 107.4 Mb and 108.5 Mb. None of the annotated genes in this region have so far been associated with height. Thus, we cannot exclude the possibility that the identified selection signature on ECA1 is associated with some trait other than height, for which Shetland ponies were selected.     

Parallel Evolution of Polydactyly Traits in Chinese and European Chickens

Polydactyly is one of the most common hereditary congenital limb malformations in chickens and other vertebrates. The zone of polarizing activity regulatory sequence (ZRS) is critical for the development of polydactyly. The causative mutation of polydactyly in the Silkie chicken has been mapped to the ZRS; however, the causative mutations of other chicken breeds are yet to be established. To understand whether the same mutation decides the polydactyly phenotype in other chicken breeds, we detected the single-nucleotide polymorphism in 26 different chicken breeds, specifically, 24 Chinese indigenous breeds and 2 European breeds. The mutation was found to have fully penetrated chickens with polydactyly in China, indicating that it is causative for polydactyly in Chinese indigenous chickens. In comparison, the mutation showed no association with polydactyly in Houdan chickens, which originate from France, Europe. Based on the different morphology of polydactyly in Chinese and European breeds, we assumed that the trait might be attributable to different genetic foundations. Therefore, we subsequently performed genome-wide association analysis (GWAS) to locate the region associated with polydactyly. As a result, a ~0.39 Mb genomic region on GGA2p was identified. The region contains six candidate genes, with the causative mutation found in Chinese indigenous breeds also being located in this region. Our results demonstrate that polydactyly in chickens from China and Europe is caused by two independent mutation events that are closely located in the chicken genome.     

Selection signature analysis implicates the PC1/PCSK1 region for chicken abdominal fat content

We conducted a selection signature analysis using the chicken 60k SNP chip in two chicken lines that had been divergently selected for abdominal fat content (AFC) for 11 generations. The selection signature analysis used multiple signals of selection, including long-range allele frequency differences between the lean and fat lines, long-range heterozygosity changes, linkage disequilibrium, haplotype frequencies, and extended haplotype homozygosity. Multiple signals of selection identified ten signatures on chromosomes 1, 2, 4, 5, 11, 15, 20, 26 and Z. The 0.73 Mb PC1/PCSK1 region of the Z chromosome at 55.43-56.16 Mb was the most heavily selected region. This region had 26 SNP markers and seven genes, Mar-03, SLC12A2, FBN2, ERAP1, CAST, PC1/PCSK1 and ELL2, where PC1/PCSK1 are the chicken/human names for the same gene. The lean and fat lines had two main haplotypes with completely opposite SNP alleles for the 26 SNP markers and were virtually line-specific, and had a recombinant haplotype with nearly equal frequency (0.193 and 0.196) in both lines. Other haplotypes in this region had negligible frequencies. Nine other regions with selection signatures were PAH-IGF1, TRPC4, GJD4-CCNY, NDST4, NOVA1, GALNT9, the ESRP2-GALR1 region with five genes, the SYCP2-CADH4 with six genes, and the TULP1-KIF21B with 14 genes. Genome-wide association analysis showed that nearly all regions with evidence of selection signature had SNP effects with genome-wide significance (P<10(-6)) on abdominal fat weight and percentage. The results of this study provide specific gene targets for the control of chicken AFC and a potential model of AFC in human obesity.     

Whole-genome scan for guttural pouch tympany in Arabian and German warmblood horses

Equine guttural pouch tympany (GPT) is a hereditary disease in foals of several breeds, including thoroughbreds, Arabian, Quarter and warmblood horses. We performed a whole-genome scan for GPT in 143 horses from five Arabian and five German warmblood families and genotyped 257 microsatellites. Chromosome-wide significant linkage was detected on ECA2 and ECA15 using multipoint non-parametric linkage analyses. Analyses stratified by sex revealed chromosome-wide significant linkage on ECA2 for fillies and chromosome-wide significant linkage on ECA15 for colts. For Arabian colts, the quantitative trait locus (QTL) on ECA15 was genome-wide significant. Haplotypes including two to four microsatellites within the QTL on ECA2 and 15 in fillies and colts, respectively, were significantly associated with GPT for both breeds. Thus, our analysis indicated sex-specific QTL, a fact which is in agreement with a two- to fourfold higher incidence of GPT in females. This is the first report of QTL for equine GPT and a first step towards identifying genes responsible for GPT.