The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560-780 thousand years before present (kyr BP). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr BP), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski's horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0-4.5 million years before present (Myr BP), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski's and domestic horse populations diverged 38-72 kyr BP, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski's horse investigated. This supports the contention that Przewalski's horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski's and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski's horse. Such regions could correspond to loci selected early during domestication.

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 – 48% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary novel centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise prior to 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.

Horses, asses, and zebras belong to a single genus, Equus, which emerged 4.0-4.5 Mya. Although the equine fossil record represents a textbook example of evolution, the succession of events that gave rise to the diversity of species existing today remains unclear. Here we present six genomes from each living species of asses and zebras. This completes the set of genomes available for all extant species in the genus, which was hitherto represented only by the horse and the domestic donkey. In addition, we used a museum specimen to characterize the genome of the quagga zebra, which was driven to extinction in the early 1900s. We scan the genomes for lineage-specific adaptations and identify 48 genes that have evolved under positive selection and are involved in olfaction, immune response, development, locomotion, and behavior. Our extensive genome dataset reveals a highly dynamic demographic history with synchronous expansions and collapses on different continents during the last 400 ky after major climatic events. We show that the earliest speciation occurred with gene flow in Northern America, and that the ancestor of present-day asses and zebras dispersed into the Old World 2.1-3.4 Mya. Strikingly, we also find evidence for gene flow involving three contemporary equine species despite chromosomal numbers varying from 16 pairs to 31 pairs. These findings challenge the claim that the accumulation of chromosomal rearrangements drive complete reproductive isolation, and promote equids as a fundamental model for understanding the interplay between chromosomal structure, gene flow, and, ultimately, speciation.equids | evolutionary genomics | speciation | admixture | chromosomal rearrangements

A reduction in the number of digits has evolved multiple times in tetrapods, particularly in cursorial mammals that travel over deserts and plains, yet the underlying developmental mechanisms have remained elusive. Here we show that digit loss can occur both during early limb patterning and at later post-patterning stages of chondrogenesis. In the “odd-toed” jerboa and horse and the “even-toed” camel, expansive cell death sculpts the tissue around the remaining toes. In contrast, digit loss in the pig is orchestrated by earlier limb patterning mechanisms including down regulation of Ptch1 expression but no increase in cell death. Together these data demonstrate remarkable plasticity in the mechanisms of vertebrate limb evolution and shed light on the complexity of morphological convergence, particularly within the artiodactyl lineage.

IntroductionThe horse is a valuable species to assess the effect of allogeneic mesenchymal stromal cells (MSCs) in regenerative treatments. No studies to date have examined recipient response to major histocompatibility complex (MHC)-mismatched equine MSCs. The purposes of this study were to immunophenotype MSCs from horses of known MHC haplotype and to compare the immunogenicity of MSCs with differing MHC class II expression.MethodsMSCs and peripheral blood leukocytes (PBLs) were obtained from Thoroughbred horses (n = 10) of known MHC haplotype (ELA-A2, -A3, and -A9 homozygotes). MSCs were cultured through P8; cells from each passage (P2 to P8) were cryopreserved until used. Immunophenotyping of MHC class I and II, CD44, CD29, CD90, LFA-1, and CD45RB was performed by using flow cytometry. Tri-lineage differentiation assays were performed to confirm MSC multipotency. Recombinant equine IFN-γ was used to stimulate MHC class II negative MSCs in culture, after which expression of MHC class II was re-examined. To assess the ability of MHC class II negative or positive MSCs to stimulate an immune response, modified one-way mixed leukocyte reactions (MLRs) were performed by using MHC-matched and mismatched responder PBLs and stimulator PBLs or MSCs. Proliferation of gated CFSE-labeled CD3+ responder T cells was evaluated via CFSE attenuation by using flow cytometry and reported as the number of cells in the proliferating T-cell gate.ResultsMSCs varied widely in MHC class II expression despite being homogenous in terms of “stemness” marker expression and ability to undergo trilineage differentiation. Stimulation of MHC class II negative MSCs with IFN-γ resulted in markedly increased expression of MHC class II. MLR results revealed that MHC-mismatched MHC class II-positive MSCs caused significantly increased responder T-cell proliferation compared with MHC-mismatched MHC class II-negative and MHC-matched MSCs, and equivalent to that of the positive control of MHC-mismatched leukocytes.ConclusionsThe results of this study suggest that MSCs should be confirmed as MHC class II negative before allogeneic application. Additionally, it must be considered that even MHC class II-negative MSCs could upregulate MHC class II expression if implanted into an area of active inflammation, as demonstrated with in vitro stimulation with IFN-γ.

The discovery of genomic imprinting through studies of manipulated mouse embryos indicated that the paternal genome has a major influence on placental development. However, previous research has not demonstrated paternal bias in imprinted genes. We applied RNA sequencing to trophoblast tissue from reciprocal hybrids of horse and donkey, where genotypic differences allowed parent-of-origin identification of most expressed genes. Using this approach, we identified a core group of 15 ancient imprinted genes, of which 10 were paternally expressed. An additional 78 candidate imprinted genes identified by RNA sequencing also showed paternal bias. Pyrosequencing was used to confirm the imprinting status of six of the genes, including the insulin receptor ( INSR ), which may play a role in growth regulation with its reciprocally imprinted ligand, histone acetyltransferase-1 ( HAT1 ), a gene involved in chromatin modification, and lymphocyte antigen 6 complex, locus G6C, a newly identified imprinted gene in the major histocompatibility complex. The 78 candidate imprinted genes displayed parent-of-origin expression bias in placenta but not fetus, and most showed less than 100% silencing of the imprinted allele. Some displayed variability in imprinting status among individuals. This variability results in a unique epigenetic signature for each placenta that contributes to variation in the intrauterine environment and thus presents the opportunity for natural selection to operate on parent-of-origin differential regulation. Taken together, these features highlight the plasticity of imprinting in mammals and the central importance of the placenta as a target tissue for genomic imprinting.

A first-generation radiation hybrid (RH) map of the equine (Equus caballus) genome was assembled using 92 horse × hamster hybrid cell lines and 730 equine markers. The map is the first comprehensive framework map of the horse that (1) incorporates type I as well as type II markers, (2) integrates synteny, cytogenetic, and meiotic maps into a consensus map, and (3) provides the most detailed genome-wide information to date on the organization and comparative status of the equine genome. The 730 loci (258 type I and 472 type II) included in the final map are clustered in 101 RH groups distributed over all equine autosomes and the X chromosome. The overall marker retention frequency in the panel is ∼21%, and the possibility of adding any new marker to the map is ∼90%. On average, the mapped markers are distributed every 19 cR (4 Mb) of the equine genome-a significant improvement in resolution over previous maps. With 69 new FISH assignments, a total of 253 cytogenetically mapped loci physically anchor the RH map to various chromosomal segments. Synteny assignments of 39 gene loci complemented the RH mapping of 27 genes. The results added 12 new loci to the horse gene map. Lastly, comparison of the assembly of 447 equine genes (256 linearly ordered RH-mapped and additional 191 FISH-mapped) with the location of draft sequences of their human and mouse orthologs provides the most extensive horse-human and horse-mouse comparative map to date. We expect that the foundation established through this map will significantly facilitate rapid targeted expansion of the horse gene map and consequently, mapping and positional cloning of genes governing traits significant to the equine industry.

IntroductionThis study tested the hypothesis that Major Histocompatibility Complex (MHC) incompatible equine mesenchymal stromal cells (MSCs) would induce cytotoxic antibodies to donor MHC antigens in recipient horses after intradermal injection. No studies to date have explored recipient antibody responses to allogeneic donor MSC transplantation in the horse. This information is critical because the horse is a valuable species for assessing the safety and efficacy of MSC treatment prior to human clinical application.MethodsSix MHC heterozygote horses were identified as non-ELA-A2 haplotype by microsatellite typing and used as allogeneic MHC-mismatched MSC recipients. MHC homozygote horses of known ELA-A2 haplotype were used as MSC and peripheral blood leukocyte (PBL) donors. One MHC homozygote horse of the ELA-A2 haplotype was the recipient of ELA-A2 donor MSCs as an MHC-matched control. Donor MSCs, which were previously isolated and immunophenotyped, were thawed and culture expanded to achieve between 30x106 and 50x106 cells for intradermal injection into the recipient’s neck. Recipient serum was collected and tested for the presence of anti-donor antibodies prior to MSC injection and every 7 days after MSC injection for the duration of the 8-week study using the standard two-stage lymphocyte microcytotoxicity dye-exclusion test. In addition to anti-ELA-A2 antibodies, recipient serum was examined for the presence of cross-reactive antibodies including anti-ELA-A3 and anti-RBC antibodies.ResultsAll MHC-mismatched recipient horses produced anti-ELA-A2 antibodies following injection of ELA-A2 MSCs and developed a wheal at the injection site that persisted for the duration of the experiment. Anti-ELA-A2 antibody responses were varied both in terms of strength and timing. Four recipient horses had high-titered anti-ELA-A2 antibody responses resulting in greater than 80% donor PBL death in the microcytotoxicity assays and one of these horses also developed antibodies that cross-reacted when tested on lymphocyte targets from a horse with an unrelated MHC type.ConclusionsAllogeneic MSCs are capable of eliciting antibody responses in vivo that can be strong and also cross-reactive with MHC types other than that of the donor. Such responses could limit the effectiveness of repeated allogeneic MSC use in a single horse, and could also result in untoward inflammatory responses in recipients.Electronic supplementary materialThe online version of this article (doi:10.1186/s13287-015-0053-x) contains supplementary material, which is available to authorized users.