Mitochondrial DNA and the origins of the domestic horse

The place and date of the domestication of the horse has long been a matter for debate among archaeologists. To determine whether horses were domesticated from one or several ancestral horse populations, we sequenced the mitochondrial D-loop for 318 horses from 25 oriental and European breeds, including American mustangs. Adding these sequences to previously published data, the total comes to 652, the largest currently available database. From these sequences, a phylogenetic network was constructed that showed that most of the 93 different mitochondrial (mt)DNA types grouped into 17 distinct phylogenetic clusters. Several of the clusters correspond to breeds and/or geographic areas, notably cluster A2, which is specific to Przewalski's horses, cluster C1, which is distinctive for northern European ponies, and cluster D1, which is well represented in Iberian and northwest African breeds. A consideration of the horse mtDNA mutation rate together with the archaeological timeframe for domestication requires at least 77 successfully breeding mares recruited from the wild. The extensive genetic diversity of these 77 ancestral mares leads us to conclude that several distinct horse populations were involved in the domestication of the 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. 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.

MultiPipMaker and supporting tools: Alignments and analysis of multiple genomic DNA sequences

Analysis of multiple sequence alignments can generate important, testable hypotheses about the phylogenetic history and cellular function of genomic sequences. We describe the MultiPipMaker server, which aligns multiple, long genomic DNA sequences quickly and with good sensitivity (available at http://bio.cse.psu.edu/ since May 2001). Alignments are computed between a contiguous reference sequence and one or more secondary sequences, which can be finished or draft sequence. The outputs include a stacked set of percent identity plots, called a MultiPip, comparing the reference sequence with subsequent sequences, and a nucleotide-level multiple alignment. New tools are provided to search MultiPipMaker output for conserved matches to a user-specified pattern and for conserved matches to position weight matrices that describe transcription factor binding sites (singly and in clusters). We illustrate the use of MultiPipMaker to identify candidate regulatory regions in WNT2 and then demonstrate by transfection assays that they are functional. Analysis of the alignments also confirms the phylogenetic inference that horses are more closely related to cats than to cows.

Evolution, systematics, and phylogeography of pleistocene horses in the new world: a molecular perspective

The rich fossil record of horses has made them a classic example of evolutionary processes. However, while the overall picture of equid evolution is well known, the details are surprisingly poorly understood, especially for the later Pliocene and Pleistocene, c. 3 million to 0.01 million years (Ma) ago, and nowhere more so than in the Americas. There is no consensus on the number of equid species or even the number of lineages that existed in these continents. Likewise, the origin of the endemic South American genus Hippidion is unresolved, as is the phylogenetic position of the “stilt-legged” horses of North America. Using ancient DNA sequences, we show that, in contrast to current models based on morphology and a recent genetic study, Hippidion was phylogenetically close to the caballine (true) horses, with origins considerably more recent than the currently accepted date of c. 10 Ma. Furthermore, we show that stilt-legged horses, commonly regarded as Old World migrants related to the hemionid asses of Asia, were in fact an endemic North American lineage. Finally, our data suggest that there were fewer horse species in late Pleistocene North America than have been named on morphological grounds. Both caballine and stilt-legged lineages may each have comprised a single, wide-ranging species.

Ancient DNA simplifies the systematics of horses present in the New Word during the Pleistocene and reveals that the stilt- legged horse is native to the continent and not a Eurasian migrant.

High altitude adaptation and phylogenetic analysis of Tibetan horse based on the mitochondrial genome

To investigate genetic mechanisms of high altitude adaptations of animals living in the Tibetan Plateau, three mitochondrial genomes (mt-genome) of Tibetan horses living in Naqu (4,500 m) of Tibetan, Zhongdian (3,300 m) and Deqin (3,100 m) of Yunnan province were sequenced. The structures and lengths of these three mt-genomes are similar to the Cheju horse, which is related to Tibetan horses, but little shorter than the Swedish horse. The pair-wise identity of these three horses on nucleotide level is more than 99.3%. When the gene encoding the mitochondrial protein of Tibetan horses was analyzed, we found that NADH6 has higher non-synonymous mutation rate in all of three Tibetan horses. This implies that NADH6 may play a role in Tibetan horses' high altitude adaptation. NADH6 is one of the subunits of the complex I in the respiratory chain. Furthermore, 7 D-loop sequences of Tibetan horse from different areas were sequenced, and the phylogeny tree was constructed to study the origin and evolutionary history of Tibetan horses. The result showed that the genetic diverse was high among Tibetan horses. All of Tibetan horses from Naqu were clustered into one clade, and Tibetan horses from Zhongdian and Deqin were clustered into others clades. The first molecular evidence of Tibetan horses indicated in this study is that Tibetan horse population might have multiple origins.

Genetic diversity and relationships of Portuguese and other horse breeds based on protein and microsatellite loci variation

There are three native Portuguese horse breeds: Lusitano, Sorraia and Garrano. This study compares diversity patterns of 17 protein and 12 microsatellite markers in these three as well as 30 other breeds to infer relationships among the breeds and to compare levels of polymorphism of these breeds for use in conservation efforts. The Garrano and the Lusitano showed a high level of genetic diversity, similar to that observed for most of the other analysed breeds, while the Sorraia and Friesian breeds showed low levels of variation for both genetic marker types. The combined protein and microsatellite data produced a tree that fit historical records well and with greater confidence levels than those for either data set alone. The combined genetic diversity and relationship information provides important baseline data for future breed conservation efforts, especially for a critically endangered breed such as the Sorraia.

Evaluation of factors affecting individual assignment precision using microsatellite data from horse breeds and simulated breed crosses

Assignment tests have been utilized to investigate population classification, measure genetic diversity and to solve forensic questions. Using microsatellite data from 26 loci genotyped in eight horse breeds we examined how population differentiation, number of scored loci, number of scored animals per breed and loci variability affected individual assignment precision applying log likelihood methods. We found that both genetic differentiation and number of scored loci were highly important for recognizing the breed of origin. When comparing two and two breeds, a proportion of 95% of the most differentiated breeds (0.200 < or = FST < or = 0.259) could be identified scoring only three loci, while the corresponding number was six for the least differentiated breeds (0.080 < or = FST < or = 0.139). An identical proportion of simulated breed crosses, differentiated from their parental breeds by FST estimates in the range 0.050-0.069, was identified when scoring 12 loci. This level of source identification was not obtained for the less differentiated breed crosses. The current data further suggested that population sample size and locus variability were not critical for the assignment precision as long as moderately large sample sizes (> or = 20 animals per population) and fairly variable loci were used.

Fixed nucleotide differences on the Y chromosome indicate clear divergence between Equus przewalskii and Equus caballus

The phylogenetic relationship between Equus przewalskii and E. caballus is often a matter of debate. Although these taxa have different chromosome numbers, they do not form monophyletic clades in a phylogenetic tree based on mtDNA sequences. Here we report sequence variation from five newly identified Y chromosome regions of the horse. Two fixed nucleotide differences on the Y chromosome clearly display Przewalski's horse and domestic horse as sister taxa. At both positions the Przewalski's horse haplotype shows the ancestral state, in common with the members of the zebra/ass lineage. We discuss the factors that may have led to the differences in mtDNA and Y-chromosomal observations.

Breed demarcation and potential for breed allocation of horses assessed by microsatellite markers

Population demarcation of eight horse breeds was investigated using genotype information of 306 horses from 26 microsatellite loci. The breeds include the indigenous Norwegian breeds Fjord Horse, Nordland/Lyngen Horse, Døle Horse and Coldblooded Trotter together with Icelandic Horse, Shetland Pony, Standardbred and Thoroughbred. Both phylogenetic analysis and a maximum likelihood method were applied to examine the potential for breed allocation of individual animals. The phylogenetic analysis utilizing simple allele sharing statistics revealed clear demarcation among the breeds; 95% of the individuals clustered together with animals of the same breed in the phylogenetic tree. Even breeds with a short history of divergence like Døle Horse and Coldblooded Trotter formed distinct clusters. Implementing the maximum likelihood method allocated 96% of the individuals to their source population, applying an assignment stringency of a log of the odds ratio larger than 2. Lower allocation stringency assigned nearly all the horses. Only three individuals were wrongly allocated a breed by both methods. In conclusion, the study demonstrates clear distinction among horse breeds, and by combining the two assignment methods breed allocation could be determined for more than 99% of the individuals.

Genetic Diversity in a Feral Horse Population from Sable Island, Canada

The present-day Sable Island horse population, inhabiting an island off the eastern coast of Canada, is believed to have originated mainly from horses confiscated from the early French settlers in Nova Scotia in the latter half of the 18th century. In 1960, the Sable Island horses were given legal protected status and no human interference has since been allowed. The objective of this study was to characterize the current genetic diversity in Sable Island horses in comparison to 15 other horse breeds commonly found in Canada and 5 Spanish breeds. A total of 145 alleles from 12 microsatellite loci were detected in 1093 horses and 40 donkeys. The average number of alleles per locus ranged from 4.67 in the Sable Island horse population to 8.25 in Appaloosas, whereas the mean observed heterozygosity ranged from 0.626 in the Sable Island population to 0.787 in Asturcons. Various genetic distance estimates and clustering methods did not permit to support that the Sable Island horses originated from shipwrecked Spanish horses, according to a popular anecdote, but closely resemble light draft and multipurpose breeds commonly found in eastern Canada. Based on the Weitzman approach, the loss of the Sable Island horse population to the overall diversity in Canada is comparable or higher than any other horse breed. The Sable Island horse population has diverged enough from other breeds to deserve special attention by conservation interest groups.

Analysis of genetic diversity and the determination of relationships among western Mediterranean horse breeds using microsatellite markers

The distribution of genetic diversity and the genetic relationships among western Mediterranean horse breeds were investigated using microsatellite markers. The examined sample included seven Spanish and three Italian local horse breeds and populations, plus a Spanish Thoroughbred outgroup. The total number of animals examined was 682 (on average 62 animals per breed; range 20-122). The microsatellite marker set analysed provided 128 alleles (10.7 alleles per locus). Within-breed genetic diversity was always high (>0.70), with breeds contributing about 8% of the total genetic variability. The mean molecular coancestry of the entire population examined was 0.205, Losino being the breed that contributed most. In addition to Nei's standard and Reynolds' genetic distances, pair-wise kinship distance and molecular coancestry were estimated. Remarkably similar breed rankings were obtained with all methods. Clustering analysis provided an accurate representation of the current genetic relationships among the breeds. Determining coancestry is useful for analysing genetic diversity distribution between and within breeds and provides a good framework for jointly analysing molecular markers and pedigree information. An integrated analysis was undertaken to obtain information on the population dynamics in western Mediterranean native horse populations, and to better determine conservation priorities.

Speciation network in Laurasiatheria: retrophylogenomic signals

Rapid species radiation due to adaptive changes or occupation of new ecospaces challenges our understanding of ancestral speciation and the relationships of modern species. At the molecular level, rapid radiation with successive speciations over short time periods-too short to fix polymorphic alleles-is described as incomplete lineage sorting. Incomplete lineage sorting leads to random fixation of genetic markers and hence, random signals of relationships in phylogenetic reconstructions. The situation is further complicated when you consider that the genome is a mosaic of ancestral and modern incompletely sorted sequence blocks that leads to reconstructed affiliations to one or the other relative, depending on the fixation of their shared ancestral polymorphic alleles. The laurasiatherian relationships among Chiroptera, Perissodactyla, Cetartiodactyla, and Carnivora present a prime example for such enigmatic affiliations. We performed whole-genome screenings for phylogenetically diagnostic retrotransposon insertions involving the representatives bat (Chiroptera), horse (Perissodactyla), cow (Cetartiodactyla), and dog (Carnivora), and extracted among 162,000 preselected cases 102 virtually homoplasy-free, phylogenetically informative retroelements to draw a complete picture of the highly complex evolutionary relations within Laurasiatheria. All possible evolutionary scenarios received considerable retrotransposon support, leaving us with a network of affiliations. However, the Cetartiodactyla-Carnivora relationship as well as the basal position of Chiroptera and an ancestral laurasiatherian hybridization process did exhibit some very clear, distinct signals. The significant accordance of retrotransposon presence/absence patterns and flanking nucleotide changes suggest an important influence of mosaic genome structures in the reconstruction of species histories.

Maternal phylogenetic relationships and genetic variation among Arabian horse populations using whole mitochondrial DNA D-loop sequencing

BACKGROUND

Maternal inheritance is an essential point in Arabian horse population genetics and strains classification. The mitochondrial DNA (mtDNA) sequencing is a highly informative tool to investigate maternal lineages. We sequenced the whole mtDNA D-loop of 251 Arabian horses to study the genetic diversity and phylogenetic relationships of Arabian populations and to examine the traditional strain classification system that depends on maternal family lines using native Arabian horses from the Middle East.

RESULTS

The variability in the upstream region of the D-loop revealed additional differences among the haplotypes that had identical sequences in the hypervariable region 1 (HVR1). While the American-Arabians showed relatively low diversity, the Syrian population was the most variable and contained a very rare and old haplogroup. The Middle Eastern horses had major genetic contributions to the Western horses and there was no clear pattern of differentiation among all tested populations. Our results also showed that several individuals from different strains shared a single haplotype, and individuals from a single strain were represented in clearly separated haplogroups.

CONCLUSIONS

The whole mtDNA D-loop sequence was more powerful for analysis of the maternal genetic diversity in the Arabian horses than using just the HVR1. Native populations from the Middle East, such as Syrians, could be suggested as a hot spot of genetic diversity and may help in understanding the evolution history of the Arabian horse breed. Most importantly, there was no evidence that the Arabian horse breed has clear subdivisions depending on the traditional maternal based strain classification system.

Changes in heart rate and heart rate variability in Thoroughbreds during prolonged road transportation

OBJECTIVE

To determine whether evaluation of heart rate (HR) and HR variability (HRV) during prolonged road transportation in horses provides a sensitive index of autonomic stimulation.

ANIMALS

Five 2-year-old Thoroughbreds.

PROCEDURE

ECGs were recorded as horses were transported for 21 hours in a 9-horse van. Heart rate, high-frequency (HF) power, low-frequency (LF) power, and LF-to-HF ratio from Fourier spectral analyses of ECGs were calculated and compared with values recorded during a 24-hour period of stall rest preceding transportation.

RESULTS

HR, HF power, and LF power had diurnal rhythms during stall rest but not during road transportation. Heart rate was higher and HF power and LF power lower during road transportation than stall rest, and HR, HF power, LF power, and LF-to-HF ratio all decreased with time during road transportation. Heart rate during stall rest was weakly and inversely associated with LF power, but during road transportation was strongly associated with LF power, HF power, and LF-to-HF ratio. Neither LF power nor HF power was correlated with LF-to-HF ratio during stall rest, but LF power was strongly and HF power weakly correlated with LF-to-HF ratio during road transportation. High-frequency power and LF power were significantly correlated with each other during stall rest and road transportation. Heart rate was significantly influenced by LF power and LF-to-HF ratio during stall rest (R(2) = 0.40) and by HF power and LF-to-HF ratio during road transportation (R(2) = 0.86).

CONCLUSIONS AND CLINICAL RELEVANCE

HR is influenced by different sympathovagal mechanisms during stall rest, compared with during road transportation; HRV may be a sensitive indicator of stress in transported horses.

The distribution of ganglion cells in the equine retina and its relationship to skull morphology

It has recently been reported that a strong correlation exists between the distribution of retinal ganglion cells and nose length in the domestic dog. To determine if this phenomenon occurs in another domestic species with diverse skull morphology, the current study examined the distribution of retinal ganglion cells in 30 horses from a variety of breeds. There was a significant variation in the density of ganglion cells found across the retinae. Breed was a significant predictor for ganglion cell density within the visual streak. A strong positive correlation exists between the density of ganglion cells in the visual streak and nasal length. Significant variation was also seen in the area centralis but did not correlate with any of the recorded skull measurements. The findings of this study provide us with further understanding of the equine visual system and the level of variation that exists between individuals of the same species.

Comparison of milk oligosaccharides pattern in colostrum of different horse breeds

Colostrum oligosaccharides are known to exhibit prebiotic and immunomodulatory properties. Oligosaccharide composition is species-specific, and equine colostrum has been reported to contain unique oligosaccharides. Therefore, equine oligosaccharides (EMOS) from colostrum from different horse breeds were analyzed by CE-LIF, CE-MS(n), HILIC-MS(n), and exoglycosidase degradation. Sixteen EMOS were characterized and quantified, of which half were neutral and half were acidic. EMOS showed about 63% structural overlap with human milk oligosaccharides, known for their bioactivity. Seven EMOS were not reported before in equine oligosaccharides literature: neutral Gal(β1-4)HexNAc, Gal(β1-4)Hex-Hex, β4'-galactosyllactose, and lactose-N-hexaose, as well as acidic 6'-Sialyl-Hex-Ac-HexNAc, sialyllacto-N-tetraose-a, and disialylacto-N-tetraose (isomer not further specified). In all colostrum samples, the average oligosaccharide concentration ranged from 2.12 to 4.63 g/L; with β 6'and 3'- galactosyllactose, 3'-sialyllactose, and disialyllactose as the most abundant of all oligosaccharides (27-59, 16-37, 1-8, and 1-6%, respectively). Differences in presence and in abundance of specific EMOS were evident not only between the four breeds but also within the breed.

Analysis of pedigree and conformation data to explain genetic variability of the horse breed Franches-Montagnes

Franches-Montagnes is the only native horse breed in Switzerland, therefore special efforts should be made for ensuring its survival. The objectives of this study were to characterize the structure of this population as well as genetic variability with pedigree data, conformation traits and molecular markers. Studies were focused to clarify if this population is composed of a heavy- and a light-type subpopulation. Extended pedigree records of 3-year-old stallions (n = 68) and mares (n = 108) were available. Evaluations of body conformation traits as well as pedigree data and molecular markers did not support the two-subpopulation hypothesis. The generation interval ranged from 7.8 to 9.3 years. The complete generation equivalent was high (>12). The number of effective ancestors varied between 18.9 and 20.1, whereof 50% of the genetic variability was attributed to seven of them. Genetic contribution of Warmblood horses ranged from 36% to 42% and that of Coldblood horses from 4% to 6%. The average inbreeding coefficient reached 6%. Inbreeding effective population size was 114.5 when the average increase of the inbreeding coefficient per year since 1910 was taken. Our results suggest that bottleneck situations occurred because of selection of a small number of sire lines. Promotion of planned matings between parents that are less related is recommended in order to avoid a reduction of the genetic diversity.

Cheek tooth morphology and ancient mitochondrial DNA of late Pleistocene horses from the western interior of North America: Implications for the taxonomy of North American Late Pleistocene Equus

Horses were a dominant component of North American Pleistocene land mammal communities and their remains are well represented in the fossil record. Despite the abundant material available for study, there is still considerable disagreement over the number of species of Equus that inhabited the different regions of the continent and on their taxonomic nomenclature. In this study, we investigated cheek tooth morphology and ancient mtDNA of late Pleistocene Equus specimens from the Western Interior of North America, with the objective of clarifying the species that lived in this region prior to the end-Pleistocene extinction. Based on the morphological and molecular data analyzed, a caballine (Equus ferus) and a non-caballine (E. conversidens) species were identified from different localities across most of the Western Interior. A second non-caballine species (E. cedralensis) was recognized from southern localities based exclusively on the morphological analyses of the cheek teeth. Notably the separation into caballine and non-caballine species was observed in the Bayesian phylogenetic analysis of ancient mtDNA as well as in the geometric morphometric analyses of the upper and lower premolars. Teeth morphologically identified as E. conversidens that yielded ancient mtDNA fall within the New World stilt-legged clade recognized in previous studies and this is the name we apply to this group. Geographic variation in morphology in the caballine species is indicated by statistically different occlusal enamel patterns in the specimens from Bluefish Caves, Yukon Territory, relative to the specimens from the other geographic regions. Whether this represents ecomorphological variation and/or a certain degree of geographic and genetic isolation of these Arctic populations requires further study.

Genetic Structure of an Endangered Portuguese Semiferal Pony Breed, the Garrano

The present study intends to survey the genetic variability of an endangered semiferal Portuguese native pony breed, the Garrano. Thirteen microsatellite markers were examined in 277 animals born in 1998, belonging to eight subpopulations corresponding to eight northern Portuguese geographic regions. Mean heterozygosity (H(o)) in the Garrano breed was 0.732, ranging from 0.531 to 0.857 across subpopulations. Allelic frequencies and diversity differed significantly between regions, suggesting the existence of genetic differentiation within the breed confirmed by the population differentiation estimator F(ST). Allele sharing genetic distance (Dps) was used to determine the relationship between the analyzed subpopulations, some of which are diverging significantly from the others. Relationships among six Iberian horse breeds (including the Garrano) were assessed through the restricted maximum likelihood method, which clusters the Garrano with another Celtic pony, the Spanish Asturcon, traditionally bred in geographical continuity with the Garrano. Results reveal a statistically significant deficit of heterozygotes within the Garrano breed (F(IT) = 0.031, p < 0.05). Regular monitoring should therefore be implemented so that the effect of genetic drift within subpopulations, enhanced by inbreeding, may be successfully minimized.

Decline of genetic diversity in ancient domestic stallions in Europe

Present-day domestic horses are immensely diverse in their maternally inherited mitochondrial DNA, yet they show very little variation on their paternally inherited Y chromosome. Although it has recently been shown that Y chromosomal diversity in domestic horses was higher at least until the Iron Age, when and why this diversity disappeared remain controversial questions. We genotyped 16 recently discovered Y chromosomal single-nucleotide polymorphisms in 96 ancient Eurasian stallions spanning the early domestication stages (Copper and Bronze Age) to the Middle Ages. Using this Y chromosomal time series, which covers nearly the entire history of horse domestication, we reveal how Y chromosomal diversity changed over time. Our results also show that the lack of multiple stallion lineages in the extant domestic population is caused by neither a founder effect nor random demographic effects but instead is the result of artificial selection-initially during the Iron Age by nomadic people from the Eurasian steppes and later during the Roman period. Moreover, the modern domestic haplotype probably derived from another, already advantageous, haplotype, most likely after the beginning of the domestication. In line with recent findings indicating that the Przewalski and domestic horse lineages remained connected by gene flow after they diverged about 45,000 years ago, we present evidence for Y chromosomal introgression of Przewalski horses into the gene pool of European domestic horses at least until medieval times.

Polymorphism identification within 50 equine gene-specific sequence tagged sites

The continued discovery of polymorphisms in the equine genome will be important for future studies using genomic screens and fine mapping for the identification of disease genes. Segments of 50 equine genes were examined for variability in 10 different horse breeds using a pool-and-sequence method. We identified 11 single nucleotide polymorphisms (SNPs) in 9380 bp of sequenced exon, and 25 SNPs, six microsatellites, and one insertion/deletion in 16961 bp of sequenced intron. Of all genes studied 52% contained at least one polymorphism, and polymorphisms were found at an overall rate of 1/613 bp. Several of the putative SNPs were tested and verified by restriction enzyme analysis using natural restriction sites or ones created by primer mutagenesis. The lowest allele frequency for a SNP detected in pooled samples was 10%. Three of the SNPs verified in the diverse horse pool were further tested in six breed-specific horse pools and were found to be reasonably variable within breeds. The pool-and-sequence method allows identification of polymorphisms in horse populations and will be a valuable tool for future disease gene and comparative mapping in horses.