Nobushige Ishida

Mitochondrial DNA Sequences of Various Species of the Genus Equus with Special Reference to the Phylogenetic Relationship Between Przewalskii's Wild Horse and Domestic Horse

The noncoding region between tRNAPro and the large conserved sequence block is the most variable region in the mammalian mitochondrial DNA D-loop region. This variable region (ca. 270 bp) of four species of Equus, including Mongolian and Japanese native domestic horses as well as Przewalskiis (or Mongolian) wild horse, were sequenced. These data were compared with our recently published Thoroughbred horse mitochondrial DNA sequences. The evolutionary rate of this region among the four species of Equus was estimated to be 2–4 × 10−8 per site per year. Phylogenetic trees of Equus species demonstrate that Przewalskiis wild horse is within the genetic variation among the domestic horse. This suggests that the chromosome number change (probably increase) of the Przewalskiis wild horse occurred rather recently.

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 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.

Immunohistochemical Localization of Chromogranin A in the Acinar Cells of Equine Salivary Glands Contrasts with Rodent Glands

We investigated the existence of chromogranin A (CgA) in salivary glands of the horse by Western blotting and enzyme immunoassay (EIA) using an antiserum against a peptide sequence of equine CgA. We also compared its cellular distribution between the horse and rat salivary glands with a tyramide signal amplification immunofluorescence technique. Western blotting gave three significant immunoreactive bands (74, 56 and 48 kDa) in adrenal medulla and three major salivary glands of horses. Immunoreactivities for CgA measured by EIA in horses were 154.05 ± 41.46, 20.32 ± 5.59 and 4.43 ± 2.23 pmol/g wet weight in the parotid gland, submandibular gland and sublingual gland, respectively, and 1.03 ± 0.407 pmol/mg protein in the saliva. Immunohistochemically, the positive reactivity was mainly recognized at acinar cells in equine salivary glands. This exhibits a contrast to the finding in the rat salivary glands that the CgA immunoreactivity is localized at the duct cells of the submandibular gland. These results provide novel evidence that in the horse, CgA is stored in the acinar cells of salivary glands, and secreted into saliva.

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.

Whole-genome linkage disequilibrium screening for complex traits in horses

The identification of candidate genes for significant traits is crucial. In this study, we developed and tested effective and systematic methods based on linkage disequilibrium (LD) for the identification of candidate regions for genes with Mendelian inheritance and those associated with complex traits. Our approach entailed the combination of primary screening using pooled DNA samples based on ΔTAC, secondary screening using an individual typing method and tertiary screening using a permutation test based on the differences in the haplotype frequency between two neighbouring microsatellites. This series of methods was evaluated using horse coat colour traits (chestnut/non-chestnut) as a simple Mendelian inheritance model. In addition, the methods were evaluated using a complex trait model constructed by mixing samples from chestnut and non-chestnut horses. Using both models, the methods could detect the expected regions for the horse coat colour trait. The results revealed that LD extends up to several centimorgans in horses, indicating that whole-genome LD screening in horses could be performed systematically and efficiently by combining the above-mentioned methods. Since genetic maps based on microsatellites have been constructed for many other species, the approaches present here could have wide applicability.

Expression of Equine Chromogranin A

In conclusion, we could confirm the expression of CGA in the equine submandibular gland. We consider that the salivary level as well as the plasma level of the CGA-LI can be a quantitative index for monitoring the activity of the sympathetic nervous system in physical and psychological stress of race horses.

Estimation models for the morbidity of the horses infected with equine influenza virus.

Estimation formulas for the morbidity of horses infected with equine influenza virus by linear regression, logistic regression and probit transformation were developed, using data from the outbreak at the Sha Tin Racing Track in Hong Kong in 1992. Using these formulas, we estimated the equine influenza virus morbidity rates at training centers belonging to the Japan Racing Association (JRA) in October 1997 and in October 1998. In 1998 JRA started a new vaccination program, and every horse must now be vaccinated twice per year. At that time, the vaccine included two US lineage virus strains, the A/equine/Kentucky/81 strain and the A/equine/La Plata/93 (LP93) strain, against equine type-2 influenza viruses; it did not include any EU lineage virus strains, such as A/equine/Suffolk/89 (SF89). Comparing the geometric mean (GM) values of hemagglutination inhibition (HI) titers between the LP93 strain and the SF89 strain in 1997 and in 1998, they both rose significantly at every age (p<0.05) by Wilcoxon test. Calculations by the simulation models show the morbidity rates for LP93 diminished from 0.439 (linear), 0.423 (logistic) and 0.431 (probit) to 0.276 (linear), 0.265 (logistic) and 0.271 (probit), respectively. On the other hand, the estimated morbidity rates for SF89 diminished only slightly from 0.954 (linear), 0.932 (logistic) and 0.944 (probit) to 0.946 (linear), 0.914 (logistic) and 0.927 (probit), respectively. Our simulation models could estimate the effect of the vaccine on each of the equine virus strains represented by the morbidity of infected horses. Thus, they are useful for vaccine evaluation.