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 F(ST)-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.

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 Przewalski's 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.

Population statistics and biological traits of endangered kiso horse.

The objective of this study was to clarify the current status of endangered Kiso horse, population statistics and biological traits, in order to take a step for the conservation by scientific approach. We surveyed 125 Kiso horses (86.2% of the whole breed), analyzed the construction of the population, and calculated the coefficient of inbreeding and effective population size. Moreover, we confirmed coat color variations and the traditional traits of the Kiso horse, and measured their height at the withers and chest circumference to clarify their physical characteristics. The population pyramid of the horses was stationary or contractive, suggesting a reduction of the population in the near future. The effective population size of the horse (47.9) suggested that the diversity was much less than their census size, and the high coefficient of inbreeding, 0.11 ± 0.07 on average, suggested that the horses were surely inbred. The horses had only 4 coat colors; bay, dark bay, buckskin dun, and chestnut, and 116 horses (92.8%) were bayish color, suggesting the fixation in their coat color. Moreover, the majority of them had dorsal stripe (83 horses; 66.4%), and the average heights at withers(131.9 ± 4.4 cm) and chest circumference (167.1 ± 10.1 cm) were not significantly different between males and females.

Construction and validation of parentage testing for thoroughbred horses by 53 single nucleotide polymorphisms.

We characterized the SNP 53 JPN System for parentage verification during horse registry. The SNP 53 JPN System was constructed using 53 highly polymorphic single nucleotide polymorphisms (SNPs), which were amplified and genotyped with 2 multiplex assays. The SNP 53 JPN System showed good resolution for 95 unrelated thoroughbreds, and the exclusion probability (PE01) for each SNP ranged from 11.5 to 23.0%, resulting in a total PE01 value of 99.996%. These results indicate that the SNP 53 JPN System is useful for parentage testing of thoroughbreds. Of the 53 SNPs, 8 SNPs could be used to exclude a pseudo parent and sib combination found using the 2006 International Society for Animal Genetics (ISAG) horse comparison test, as efficiently as the parentage testing systems using short tandem repeats (STRs). Thus, we concluded that the SNP 53 JPN System could provide sufficient and reliable information for routine parentage testing of thoroughbred.

Demographic analysis of breeding structure in Japanese thoroughbred population.

To investigate the breeding structure in the Japanese Thoroughbred population, we applied a demographic analysis to the populations of foals produced from 1978 to 2005. The migration rate estimated from the proportion of foals produced by imported breeding horses was around 40% over the investigated period. After early 1990s, the migration rate through stallions imported from USA sharply increased. The average generation interval was within range of 10.5-11.5 years. The longer generation interval of Thoroughbred was considered to be a reflection of the fact that Thoroughbred horses begin breeding only after completing their performance in races. After the peak of 729 in 1993, the number of sires of foals progressively declined to 358 in 2005. Although the coefficient of variation of the progeny number of sires was within range of 1.0-1.2 until early 1990s, it gradually increased and reached the value of 1.6-1.7 in recent years. The effective number of sires consistently decreased after the peak of 302.6 in 1992, and reached 120-130 in recent years, which is 25-30% of the actual number of sires. In parallel, the demographic estimate of the effective population size declined after early 1990s. The main cause of the observed change in the breeding structure was inferred to be the intensive use of a limited number of stallions for breeding.

Gene dropping analysis of ancestral contributions and allele survival in Japanese thoroughbred population.

Genetic contributions of nine historically important ancestors and allelic diversity in the Japanese Thoroughbred population were examined by applying the gene dropping simulation to the foals produced from 1978 to 2005. Full pedigree records traced to ancestors (base animals) born around 1890 were used for the simulation. Alleles originated from some of the historically important ancestors were found to be at risk of future extinction, although their genetic contributions to the foal population have increased during the last three decades. The proportion of surviving alleles to the total alleles assigned to the base animals was 8.0% in the foal population in 2005, suggesting that a large part of genetic variability contained in the base animals is extinct in the current population.

Type VII and XVII Collagen mRNA Expressions in Regenerated Epidermal Laminae in Chronic Equine Laminitis.

To confirm ability forming the basement membrane of the regenerated laminar epidermis (rLE) in chronic laminitis, expression of type VII and type XVII collagen mRNAs in the rLE was studied applying sequences of two type of murine collagens. On northern blot analysis, complement DNA (cDNA) probes adjusted from the murine type VII and type XVII collagen could hybridize with the equine mRNAs, and each signal was detected as single-bands at approximately 9.5 kb and 5.6 kb, respectively. Contrasting with the expression level of equine glyceraldehyde-3-phosphate dehydrogenease mRNA, the band of type VII collagen mRNA in laminitis was stronger than normal, but the type XVII collagen mRNA in laminitis was less than normal. By in situ hybridization, positive signals in response to the murine type VII and type XVII collagen mRNA probes could be detected in the equine laminitic rLE region. From these results, it is concluded that the keratinocytes constructing the rLE in chronic stage of laminitis can express type VII and type XVII collagen mRNAs and these expression patterns were different from the normal.

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.

Improved resolution of the comparative horse-human map: investigating markers with in silico and linkage mapping approaches.

Genetic maps are extremely important tools for tracing the genes that govern economically significant traits, and microsatellites are a significant component of these. In this study, we isolated 2346 novel horse microsatellites as resources for the construction of high-density horse genetic maps. Of these 2346 markers, 339 (14.5%) horse sequences showed sequence homology to DNA sequences in the human genome, demonstrating that microsatellites as type II markers are valuable resources for developing linkage maps and that they have a potential equal to that of type I markers for developing comparative maps. Of the 339 markers, 206 (60.8%) were assigned to horse chromosomes using the Animal Health Trust (AHT) full-sib reference family, and 195 (94.6%) of these localized to the expected syntenic locations on the human genome. These results confirmed the high level of accuracy of in silico mapping. Thus, the 339 markers that exhibited homology to the human genome increased the density of markers on the horse-human comparative map. The resulting comparative map will facilitate the use of horse microsatellites as genetic markers for the identification of quantitative trait loci (QTL) that have been mapped on the human genome. In addition, although the in silico and linkage mapping data did not agree for the other 11 (5.4%) of the assigned 206 markers, these may represent new putative regions of horse-human synteny.

Polymorphism identification, RH mapping, and association analysis with the anxiety trait of the equine serotonin transporter (SLC6A4) gene.

Equine anxiety trait is considered an important temperament in various situations, including riding, training, and daily care. This study examined the polymorphism of the equine serotonin transporter (SLC6A4) gene as a candidate genetic element influencing equine anxiety trait. The sequence of the coding region of this gene was highly homologous with those of other mammals, and four single nucleotide polymorphisms were found by comparing the sequences of ten genetically unrelated thoroughbred horses. Radiation hybrid mapping revealed that this gene was located 26.92 cR from neurofibromin 1 on ECA 11. Using two-year-old thoroughbred horses (n=67), the association of these polymorphisms with the anxiety trait was examined, but no significant association was identified between each haplotype of the serotonin transporter gene and the anxiety score.

Prospects for whole genome linkage disequilibrium mapping in thoroughbreds.

Linkage disequilibrium (LD) mapping is often used in searches for genes governing economically significant traits and diseases. The D' coefficient is a commonly used measure of the extent of LD between all possible pairs of alleles at two markers. This study aimed to test the utility of the D' coefficient for LD mapping of a trait in a thoroughbred population. Microsatellite genotype data and grey coat colour as a trait model in a thoroughbred population were used to assess the extent of LD. We demonstrated that LD mapping was a reasonable approach for initial genome-wide scans in a thoroughbred population. Significant LD was demonstrated at approximately 7 cM, implying that roughly 430 appropriately spaced microsatellites were needed for systematic whole-genome LD mapping in this model. LD mapping methods using D' in a thoroughbred population were useful for identifying the chromosomal regions for diseases and economic trait loci (ETL). It was suggested that a thoroughbred population represented a population particularly suitable for LD mapping.

Molecular cloning, nucleotide sequence and presence of multiple functional polyadenylation signals in the 3'-untranslated region of equine dopamine beta-hydroxylase cDNA.

Complementary DNA (cDNA) encoding equine dopamine beta-hydroxylase (DBH) was amplified with a combination of reverse transcription-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE) method, and their nucleotide sequences (Accession No. AB029430: the DDBJ nucleotide sequence database) was determined. A total of 3842 bp cDNA sequence was consisted with 5 bp of 5' flanking untranslated sequence, 1833 bp of open reading frame encoding 610 amino acids, and 2004 bp of 3' flanking untranslated sequence. The deduced amino acid sequence of equine DBH was very similar to the known mammalian DBH sequences. The similarity between amino acid sequence of equine DBH to sequences of bovine, human, rat and mouse DBH were 86.3, 84.6, 82.2 and 81.2%, respectively. Northern blot analysis and in situ hybridization revealed three different sizes of mRNA expressions in equine adrenal medulla tissue. Then we found three putative polyadenylation signal sites in the 3' flanking untranslated sequence. These results indicate that alternative use of three polyadenylation sites generates the equine DBH mRNA that have different sizes of 3' flanking untranslated region. These results may provide further evidence for understanding DBH molecule and clues for the equine DBH gene analysis.

Determination of nitrotyrosine by HPLC-ECD and its application.

3-nitrotyrosine, a product of tyrosine nitration, is a useful indicator of oxidative damage. We modified the previously reported HPLC-electrochemical detection (ECD) method: specifically, a through-type porous carbon electrode was used as a reducing electrode instead of the mercury-gold amalgam electrode, because the response of the latter changes over time. A combination of reverse-phase HPLC and electrochemical detector passed through -800 mV reduction potential and subsequently under +250 mV oxidation potential allows measurement of 3-nitrotyrosine. The detection limit of this assay was less than 10 fmol. In mice to which lipopolysaccharide (LPS) was administered intraperitoneally, plasma 3-nitrotyrosine levels were elevated, corresponding to LPS dosage. These findings suggest that the improved HPLC-ECD method can be used as a specific and sensitive assay of biological 3-nitrotyrosine and can be applied clinically.

Nucleotide sequence of equine erythropoietin and characterization of region-specific antibodies.

OBJECTIVE: To determine the full-length complementary DNA (cDNA) sequence of equine erythropoietin (EPO) and to develop region-specific antibodies to differentiate equine EPO (eEPO) and human EPO (hEPO). SAMPLE POPULATION: RNA and lysate extracted from renal tissues of an adult Thoroughbred. PROCEDURE: Full-length cDNA was determined by use of a reverse transcriptase-polymerase chain reaction assay and a rapid amplification of cDNA ends method. The deduced amino acid sequence was compared with sequences of EPO reported for other species. Furthermore, 4 synthetic peptides were designed in 2 distinctive parts of the eEPO and hEPO amino acid sequences to obtain antibodies specific for eEPO and hEPO. Specificity of the antibodies was tested against supernatant of homogenized equine kidney and recombinant hEPO (rhEPO) by use of western immunoblotting techniques. RESULTS: Analysis of the 1,181 bp in the nucleotide sequence revealed that eEPO was a residue of 192 amino acids. Similarity of eEPO with amino acid sequences of EPO from other species was 81.0% to 90.6%. Antibodies were specifically recognized by eEPO or rhEPO molecules. Anti-hEPO (161 to 165) antibody specifically recognized rhEPO. In contrast, anti-eEPO (133 to 144) antibody reacted with the equine kidney lysate. CONCLUSIONS AND CLINICAL RELEVANCE: We determined the cDNA and amino acid sequence of eEPO and developed region-specific antibodies that specifically recognized eEPO or rhEPO. These antibodies may be useful in distinguishing rhEPO from eEPO in a test to detect the misuse of rhEPO in racehorses.

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. F(ST) 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.

Genetic characterization of the endangered Kiso horse using 31 microsatellite DNAs.

In order to contribute to conservation of the endangered Kiso horse, we clarified their genetic information using 31 microsatellite DNAs, and genotyped 125 horses, 83% of the existing breed. First, we clarified the current status of the horses. The horses were confirmed to have experienced rapid loss of population causing a bottleneck, and their effective population size was much smaller than their census size. Moreover, the number of alleles (6.3), observed heterozygosity (0.674), and expected heterozygosity (0.662) were in the same range as other endangered horses all over the world. Therefore, although their inbreeding level was not so severe (F(is): -0.017), the Kiso horse is surely one of the endangered. Second, we obtained genetic information of individuals. This information allowed us to understand the genetic distance of individuals, and might help in development of a reproductive strategy concerning the genetic distance between the mating pairs. Moreover, there appeared to be 4 subpopulations of Kiso horse, and this result was in good agreement with their historical background. Third, we confirmed that the parentage test for identification using the 31 microsatellite DNAs was highly reliable (probability of exclusion: 0.999999993). This identification increases the reliability of stud certification, and is also helpful for effective management. Understanding the genetic diversity within the population and the relationships among individuals is important to ensuring effective management for maintenance of genetic variation, and this study may help in conservation of the endangered Kiso horse.

In vitro development of equine oocytes from preserved ovaries after intracytoplasmic sperm injection.

In this study, we evaluated the meiotic competence of equine oocytes from ovaries preserved for one day. We also investigated fertilization, cleavage rate, developmental competence and freezability of equine embryos after intracytoplasmic sperm injection (ICSI). After collection from ovaries, the oocytes were classified into two groups comprised of those having compact cumulus layers (Cp) or those having expanded cumulus layers (Ex). Oocytes with a first polar body were subjected to fertilization by ICSI using frozen-thawed stallion spermatozoa and were then cultured in CR1aa medium. The rates of metaphase II-stage oocytes, normal fertilization and cleavage were not significantly different between the two oocyte categories (38.5, 70.0 and 48.7% for CP and 43.5, 60.0 and 58.8% for Ex, respectively). However, the blastocyst development rate of Ex was significantly (P<0.05) higher than that of Cp (25.5 vs. 7.7%). Three Cp-derived and 12 Ex-derived early blastocysts were cryopreserved using the slow cooling protocol, and all of them developed to hatching blastocysts after thawing. These results suggest that equine oocytes fertilized by ICSI can develop to the preimplantation stage in culture conditions similar to those used in the bovine. Furthermore, the Ex oocytes had higher developmental competence than the Cp oocytes, and the in vitro-produced blastocysts had high viability after freezing and thawing.

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

Single linkage group per chromosome genetic linkage map for the horse, based on two three-generation, full-sibling, crossbred horse reference families.

A genetic linkage map of the horse consisting of 742 markers, which comprises a single linkage group for each of the autosomes and the X chromosome, is presented. The map has been generated from two three-generation full-sibling reference families, sired by the same stallion, in which there are 61 individuals in the F2 generation. Each linkage group has been assigned to a chromosome and oriented with reference to markers mapped by fluorescence in situ hybridization. The average interval between markers is 3.7 cM and the linkage groups collectively span 2772 cM. The 742 markers comprise 734 microsatellite and 8 gene-based markers. The utility of the microsatellite markers for comparative mapping has been significantly enhanced by comparing their flanking sequences with the human genome sequence; this enabled conserved segments between human and horse to be identified. The new map provides a valuable resource for genetically mapping traits of interest in the horse.

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.