Telhisa Hasegawa

Genome Sequence, Comparative Analysis, and Population Genetics of the Domestic Horse

A Horse Is a Horse, of Course The history of horse domestication is closely tied to the history of the human society. Wade et al. (p. 865) report on the sequencing and provide a single nucleotide polymorphism map of the horse (Equus caballus) genome. Horses are a member of the order perissodactyla (odd-toed animals with hooves). The analysis reveals an evolutionarily new centromere on equine chromosome 11 that displays properties of an immature but fully functioning centromere and is devoid of centromeric satellite sequence. The findings clarify the nature of genetic diversity within and across horse breeds and suggest that the horse was domesticated from a relatively large number of females, but few males. The horse genome reveals an evolutionary new centromere and conserved chromosomal sequences relative to other mammals. We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

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

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

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

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

Equine inhibin/activin βA-subunit mRNA is expressed in the endometrial gland, but not in the trophoblast, during pregnancy

The expression of both inhibin α‐ and inhibin/activin βA‐subunit mRNA was examined in equine uteroplacental tissues collected during pregnancy (days 90 to 300). Northern blot analysis revealed that 5 transcripts (7.0, 4.1, 3.4, 2.6, 1.5 kb) of βA‐subunit were present, and the most abundantly expressed transcript was the 1.5 kb one. Relatively high levels of the 1.5 kb transcript were seen in the second trimester of pregnancy compared to what was found in the third trimester. To identify the tissue localization of βA‐subunit mRNA, in situ hybridization was performed, and the positive signal was observed exclusively in the endometrial glands, but not in the fetal placental tissue (trophoblast) at days 150, 210, and 300 of pregnancy. On the other hand, inhibin α‐subunit transcript could not be detected at any stage of pregnancy examined either by Northern blot analysis or in situ hybridization. Although the factor(s) regulating the gene expression of βA‐subunit in this equine tissue is currently unknown, these results suggest that activin, but not inhibin, is predominantly produced in the endometrial glands of the pregnant mare, and thus produced activin may play a paracrine or endocrine role during pregnancy in this species. Mol. Reprod. Dev. 47:363–369, 1997.

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.