Simulated validation of intron-less transgene detection using DELLY for gene-doping control in horse sports.

Gene doping is prohibited in horseracing. In a previous study, we developed a method for non-targeted transgene detection using DELLY, which is based on split-read (SR) and paired-end (PE) algorithms to detect structural variants, on WGS data. In this study, we validated the detection sensitivity of DELLY using artificially generated sequence data of 12 target genes. With DELLY, at least one intron was detected as a deletion in eight targeted genes using the 150 bp PE read WGS data, whereas all targeted genes were detected by DELLY using the 100 bp PE read data. The detection sensitivity was higher in 100 bp PE reads than in 150 bp PE reads, despite a lower total sequence coverage, probably because of mismatch tolerance between the mapped reads and reference genome. In addition, it was observed that the average intron size detected by SR alone was 293 bp and that that detected by both SR and PE was 8924 bp. Thus, we showed that transgenes with various intron-exon structures could be detected using DELLY, suggesting its application in gene-doping control in horses.

A candidate-SNP retrospective cohort study for fracture risk in Japanese Thoroughbred racehorses.

Fractures are medical conditions that compromise the athletic potential of horses and/or the safety of jockeys. Therefore, the reduction of fracture risk is an important horse and human welfare issue. The present study used molecular genetic approaches to determine the effect of genetic risk for fracture at four candidate SNPs spanning the myostatin (MSTN) gene on horse chromosome 18. Among the 3706 Japanese Thoroughbred racehorses, 1089 (29.4%) had experienced fractures in their athletic life, indicating the common occurrence of this injury in Thoroughbreds. In the case/control association study, fractures of the carpus (carpal bones and distal radius) were statistically associated with g.65809482T/C (P = 1.17 x 10-8 ), g.65868604G/T (P = 2.66 x 10-9 ), and g.66493737C/T (P = 6.41 x 10-8 ). In the retrospective cohort study using 1710 racehorses born in 2000, the relative risk (RR) was highest for male horses at g.65868604G/T, based on the dominant allele risk model (RR = 2.251, 95% confidence interval 1.407-3.604, P = 0.00041), and for female horses at g.65868604G/T, based on the recessive allele risk model (RR = 2.313, 95% confidence interval 1.380-3.877, P = 0.00163). Considering the association of these SNPs with racing performance traits such as speed, these genotypes may affect the occurrence of carpus fractures in Japanese Thoroughbred racehorses as a consequence of the non-genetic influence of the genotype on the distance and/or intensity of racing and training. The genetic information presented here may contribute to the development of strategic training programs and racing plans for racehorses that improve their health and welfare.

Genetic diversity and relationships among native Japanese horse breeds, the Japanese Thoroughbred and horses outside of Japan using genome-wide SNP data.

Eight horse breeds-Hokkaido, Kiso, Misaki, Noma, Taishu, Tokara, Miyako and Yonaguni-are native to Japan. Although Japanese native breeds are believed to have originated from ancient Mongolian horses imported from the Korean Peninsula, the phylogenetic relationships among these breeds are not well elucidated. In the present study, we compared genetic diversity among 32 international horse breeds previously evaluated by the Equine Genetic Diversity Consortium, the eight Japanese native breeds and Japanese Thoroughbreds using genome-wide SNP genotype data. The proportion of polymorphic loci and expected heterozygosity showed that the native Japanese breeds, with the exception of the Hokkaido, have relatively low diversity compared to the other breeds sampled. Phylogenetic and cluster analyses demonstrated relationships among the breeds that largely reflect their geographic distribution in Japan. Based on these data, we suggest that Japanese horses originated from Mongolian horses migrating through the Korean Peninsula. The Japanese Thoroughbreds were distinct from the native breeds, and although they maintain similar overall diversity as Thoroughbreds from outside Japan, they also show evidence of uniqueness relative to the other Thoroughbred samples. This is the first study to place the eight native Japanese breeds and Japanese Thoroughbred in context with an international sample of diverse breeds.

The evolutionary history of the DMRT3 'Gait keeper' haplotype.

A previous study revealed a strong association between the DMRT3:Ser301STOP mutation in horses and alternate gaits as well as performance in harness racing. Several follow-up studies have confirmed a high frequency of the mutation in gaited horse breeds and an effect on gait quality. The aim of this study was to determine when and where the mutation arose, to identify additional potential causal mutations and to determine the coalescence time for contemporary haplotypes carrying the stop mutation. We utilized sequences from 89 horses representing 26 breeds to identify 102 SNPs encompassing the DMRT3 gene that are in strong linkage disequilibrium with the stop mutation. These 102 SNPs were genotyped in an additional 382 horses representing 72 breeds, and we identified 14 unique haplotypes. The results provided conclusive evidence that DMRT3:Ser301STOP is causal, as no other sequence polymorphisms showed an equally strong association to locomotion traits. The low sequence diversity among mutant chromosomes demonstrated that they must have diverged from a common ancestral sequence within the last 10 000 years. Thus, the mutation occurred either just before domestication or more likely some time after domestication and then spread across the world as a result of selection on locomotion traits.

Evidence for the effect of serotonin receptor 1A gene (HTR1A) polymorphism on tractability in Thoroughbred horses.

Tractability, or how easily animals can be trained and controlled, is an important behavioural trait for the management and training of domestic animals, but its genetic basis remains unclear. Polymorphisms in the serotonin receptor 1A gene (HTR1A) have been associated with individual variability in anxiety-related traits in several species. In this study, we examined the association between HTR1A polymorphisms and tractability in Thoroughbred horses. We assessed the tractability of 167 one-year-old horses reared at a training centre for racehorses using a questionnaire consisting of 17 items. A principal components analysis of answers contracted the data to five principal component (PC) scores. We genotyped two non-synonymous single nucleotide polymorphisms (SNPs) in the horse HTR1A coding region. We found that one of the two SNPs, c.709G>A, which causes an amino acid change at the intracellular region of the receptor, was significantly associated with scores of four of five PCs in fillies (all Ps < 0.05) and one PC in colts (P < 0.01). Horses carrying an A allele at c.709G>A showed lower tractability. This result provides the first evidence that a polymorphism in a serotonin-related gene may affect tractability in horses with the effect partially different depending on sex.

Worldwide frequency distribution of the 'Gait keeper' mutation in the DMRT3 gene.

For centuries, domestic horses have represented an important means of transport and served as working and companion animals. Although their role in transportation is less important today, many horse breeds are still subject to intense selection based on their pattern of locomotion. A striking example of such a selected trait is the ability of a horse to perform additional gaits other than the common walk, trot and gallop. Those could be four-beat ambling gaits, which are particularly smooth and comfortable for the rider, or pace, used mainly in racing. Gaited horse breeds occur around the globe, suggesting that gaitedness is an old trait, selected for in many breeds. A recent study discovered that a nonsense mutation in DMRT3 has a major impact on gaitedness in horses and is present at a high frequency in gaited breeds and in horses bred for harness racing. Here, we report a study of the worldwide distribution of this mutation. We genotyped 4396 horses representing 141 horse breeds for the DMRT3 stop mutation. More than half (2749) of these horses also were genotyped for a SNP situated 32 kb upstream of the DMRT3 nonsense mutation because these two SNPs are in very strong linkage disequilibrium. We show that the DMRT3 mutation is present in 68 of the 141 genotyped horse breeds at a frequency ranging from 1% to 100%. We also show that the mutation is not limited to a geographical area, but is found worldwide. The breeds with a high frequency of the stop mutation (>50%) are either classified as gaited or bred for harness racing.

Empirical growth curve estimation considering multiple seasonal compensatory growths of body weights in Japanese Thoroughbred colts and fillies.

Thoroughbred horses are seasonal mating animals, and their foals are born yearly in spring seasons. In northern regions or countries, the foals generally show a typical seasonal compensatory growth pattern, where their growth rate declines in winter and increases in the next spring. In this study, a new empirical approach is proposed to adjust for this compensatory growth when growth curve equations are estimated, by using BW of Japanese Thoroughbred colts and fillies raised in Hidaka, Hokkaido. Based on the traditional Richards growth curve equation, new growth curve equations were developed and fit to the weight-age data. The foals generally experience 2 major winter seasons before their debut in horseracing. The new equations had sigmoid subfunctions that can empirically adjust the first and second year compensatory growths, combined with the Richards biological parameter responsible for the maturity of animals. The unknown parameters included in the equations were estimated by SAS NLMIXED procedure. The goodness-of-fit was examined by using several indices of goodness-of-fit (i.e., Akaike's information criterion, Bayesian information criterion, -2 log likelihood, and residual sum of squares) for the multiple applications of the subfunctions. The indices indicated the best fit of the new equations including both subfunctions for the first and second compensatory growths to the weight-age data. The shapes of the growth curves were improved during the periods of compensatory growth. The proposed method is one of the useful approaches for adjusting multiple seasonal compensatory growths in growth curve estimations of Thoroughbreds and for the management of young horses during the compensatory periods.

Heritability estimates for racing performance in Japanese Thoroughbred racehorses using linear and non-linear model analyses.

This study evaluated the differences between linear and non-linear modelled heritability estimates of racing performance based on lifetime earnings (LE) and lifetime ranking (LR) in Japanese Thoroughbred racehorses. The heritability estimate (h(2) = 0.25) obtained from a non-linear model based on formal Japan Racing Association ranking was much higher than that obtained from a linear model based on the original trait phenotype (h(2) = 0.11). The linear models showed slightly higher heritability estimates under the trait categorizations than under the original phenotypes, while the non-linear categorical trait models showed much higher heritability estimates than the linear models, especially for binary trait categorizations (h(2) = 0.34) with non-winning and winning horses. The binary trait categorizations were consistent with the case and control classifications in the previous genome-wide association study (GWAS), which identified possible sequence variants on ECA18 that affect racing performance in Japanese Thoroughbred racehorses. Those findings suggested that the different heritability estimates obtained from several trait categorizations would reflect the possible presence of susceptibility gene segregations in the analyzed population, indicating that heritability estimates from non-linear models are useful for the selection of case and control populations in GWAS.

A cohort study of racing performance in Japanese Thoroughbred racehorses using genome information on ECA18.

Using 1710 Thoroughbred racehorses in Japan, a cohort study was performed to evaluate the influence of genotypes at four single nucleotide polymorphisms (SNPs) on equine chromosome 18 (ECA18), which were associated in a previous genome-wide association study for racing performance with lifetime earnings and performance rank. In males, both g.65809482T>C and g.65868604G>T were related to performance rank (P= 0.005). In females, g.65809482T>C (P = 1.76E-6), g.65868604G>T (P=6.81E-6) and g.66493737C>T (P=4.42E-5) were strongly related to performance rank and also to lifetime earnings (P < 0.05). When win-race distance (WRD) among all winning racehorses and best race distance (BRD) among elite racehorses were considered as the phenotypes, significant associations (P<0.001) were observed for all four SNPs. The favourable race distance of both elite (BRD) and novice racehorses (WRD) was also associated with genotypes in the ECA18 region, indicating the presence of a gene in this region influencing optimum race distance in Thoroughbred racehorses. Therefore, the association with performance rank is likely due to the bias in the race distances. The location of the SNPs within and proximal to the gene encoding myostatin (MSTN) strongly suggests that regulation of the MSTN gene affects racing performance. In particular, the g.65809482T>C, g.65868604G>T and g.66493737C>T SNPs, or their combinations, may be genetic diagnostic markers for racing performance indicators such as WRD and BRD.

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

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

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.

Development of a method for simultaneously genotyping multiple horse coat colour loci and genetic investigation of basic colour variation in Thoroughbred and Misaki horses in Japan.

In order to develop a genotyping method that can be used in the registration procedure for Thoroughbreds, we developed a method for simultaneously genotyping multiple coat colour genes on the basis of single nucleotide polymorphism typing by using the SNaPshot(TM) technique. This method enabled precise and reasonable detection of causal mutations; it was effective for genotyping of MC1R, ASIP, and SLC45A2 at the Extension (E), Agouti (A), Cream dilution (C) loci, and the possibility of identification of rare variants of MC1R, EDNRB and KIT at the E, Overo (O) and Sabino 1 (SB1) loci, respectively, was also indicated. It was considered that this genotyping method would provide information not only for the registration of Thoroughbreds but also for the preservation of phenotypic characters, such as coat colour, of endangered Misaki native horses in Japan. Therefore, genetic variations at the five coat colour loci were investigated in 1111 Thoroughbred and 99 Misaki native horses. Allele frequencies at the polymorphic E and A loci were estimated, and the proportions of basic coat colours that could be expected in the Thoroughbred population were bay, 0.662; black, 0.070; chestnut, 0.268. In the Misaki population, they were bay, 0.792; black, 0.129; chestnut, 0.080. The data presented were the first of its kind on genetic coat colour variation, and will be important with regard to the registration of Thoroughbreds and the management of Misaki horses.

Seven novel KIT mutations in horses with white coat colour phenotypes.

White coat colour in horses is inherited as a monogenic autosomal dominant trait showing a variable expression of coat depigmentation. Mutations in the KIT gene have previously been shown to cause white coat colour phenotypes in pigs, mice and humans. We recently also demonstrated that four independent mutations in the equine KIT gene are responsible for the dominant white coat colour phenotype in various horse breeds. We have now analysed additional horse families segregating for white coat colour phenotypes and report seven new KIT mutations in independent Thoroughbred, Icelandic Horse, German Holstein, Quarter Horse and South German Draft Horse families. In four of the seven families, only one single white horse, presumably representing the founder for each of the four respective mutations, was available for genotyping. The newly reported mutations comprise two frameshift mutations (c.1126_1129delGAAC; c.2193delG), two missense mutations (c.856G>A; c.1789G>A) and three splice site mutations (c.338-1G>C; c.2222-1G>A; c.2684+1G>A). White phenotypes in horses show a remarkable allelic heterogeneity. In fact, a higher number of alleles are molecularly characterized at the equine KIT gene than for any other known gene in livestock species.

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.

International Equine Gene Mapping Workshop Report: a comprehensive linkage map constructed with data from new markers and by merging four mapping resources.

A comprehensive male linkage map was generated by adding 359 new, informative microsatellites to the International Equine Gene Map half-sibling reference families and by combining genotype data from three independent mapping resources: a full sibling family created at the Animal Health Trust in Newmarket, United Kingdom, eight half-sibling families from Sweden and two half-sibling families from the University of California, Davis. Because the combined data were derived primarily from half-sibling families, only autosomal markers were analyzed. The map was constructed from a total of 766 markers distributed on the 31 equine chromosomes. It has a higher marker density than that of previously reported maps, with 626 markers linearly ordered and 140 other markers assigned to a chromosomal region. Fifty-nine markers (7%) failed to meet the criteria for statistical evidence of linkage and remain unassigned. The map spans 3,740 cM with an average distance of 6.3 cM between markers. Fifty-five percent of the intervals are < or = 5 cM and only 3% > or = 20 cM. The present map demonstrates the cohesiveness of the different data sets and provides a single resource for genome scan analyses and integration with the radiation hybrid map.