Refining the evolutionary tree of the horse Y chromosome.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

Juvenile idiopathic epilepsy in Arabian horses is not a single-gene disorder.

Valued for their temperament, beauty, athletic ability, and exhibition in the show ring, Arabian horses are an important component of the horse industry. Juvenile idiopathic epilepsy (JIE), a seizure disorder, is most often reported in Arabian foals from birth to 6 months of age. Affected foals exhibit tonic-clonic seizures lasting as long as 5 min and risking secondary complications like temporary blindness and disorientation. Some foals outgrow this condition, while others die or suffer lifelong complications if not treated. Previous work suggested a strong genetic component to JIE and proposed JIE to be a single-gene trait. In this work, we conducted a genome wide association study (GWAS) in 60 cases of JIE and 120 genetically matched controls, identifying loci suggesting JIE is not caused by a single locus. Coat color (chestnut, gray) phenotypes were used as positive control traits to assess the efficacy of GWAS in this population. Future work will attempt to future define candidate regions and explore a polygenic mode of inheritance.

Prediction of histone post-translational modification patterns based on nascent transcription data.

The role of histone modifications in transcription remains incompletely understood. Here, we examine the relationship between histone modifications and transcription using experimental perturbations combined with sensitive machine-learning tools. Transcription predicted the variation in active histone marks and complex chromatin states, like bivalent promoters, down to single-nucleosome resolution and at an accuracy that rivaled the correspondence between independent ChIP-seq experiments. Blocking transcription rapidly removed two punctate marks, H3K4me3 and H3K27ac, from chromatin indicating that transcription is required for active histone modifications. Transcription was also required for maintenance of H3K27me3, consistent with a role for RNA in recruiting PRC2. A subset of DNase-I-hypersensitive sites were refractory to prediction, precluding models where transcription initiates pervasively at any open chromatin. Our results, in combination with past literature, support a model in which active histone modifications serve a supportive, rather than an essential regulatory, role in transcription.

Y-Chromosomal Insights into Breeding History and Sire Line Genealogies of Arabian Horses.

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

Genome Diversity and the Origin of the Arabian Horse.

The Arabian horse, one of the world's oldest breeds of any domesticated animal, is characterized by natural beauty, graceful movement, athletic endurance, and, as a result of its development in the arid Middle East, the ability to thrive in a hot, dry environment. Here we studied 378 Arabian horses from 12 countries using equine single nucleotide polymorphism (SNP) arrays and whole-genome re-sequencing to examine hypotheses about genomic diversity, population structure, and the relationship of the Arabian to other horse breeds. We identified a high degree of genetic variation and complex ancestry in Arabian horses from the Middle East region. Also, contrary to popular belief, we could detect no significant genomic contribution of the Arabian breed to the Thoroughbred racehorse, including Y chromosome ancestry. However, we found strong evidence for recent interbreeding of Thoroughbreds with Arabians used for flat-racing competitions. Genetic signatures suggestive of selective sweeps across the Arabian breed contain candidate genes for combating oxidative damage during exercise, and within the "Straight Egyptian" subgroup, for facial morphology. Overall, our data support an origin of the Arabian horse in the Middle East, no evidence for reduced global genetic diversity across the breed, and unique genetic adaptations for both physiology and conformation.

Erratum: Author Correction: Improved reference genome for the domestic horse increases assembly contiguity and composition.

[This corrects the article DOI: 10.1038/s42003-018-0199-z.].

Improved reference genome for the domestic horse increases assembly contiguity and composition.

Recent advances in genomic sequencing technology and computational assembly methods have allowed scientists to improve reference genome assemblies in terms of contiguity and composition. EquCab2, a reference genome for the domestic horse, was released in 2007. Although of equal or better quality compared to other first-generation Sanger assemblies, it had many of the shortcomings common to them. In 2014, the equine genomics research community began a project to improve the reference sequence for the horse, building upon the solid foundation of EquCab2 and incorporating new short-read data, long-read data, and proximity ligation data. Here, we present EquCab3. The count of non-N bases in the incorporated chromosomes is improved from 2.33 Gb in EquCab2 to 2.41 Gb in EquCab3. Contiguity has also been improved nearly 40-fold with a contig N50 of 4.5 Mb and scaffold contiguity enhanced to where all but one of the 32 chromosomes is comprised of a single scaffold.

Ectopic Trophoblast Allografts in the Horse Resist Destruction by Secondary Immune Responses.

Invasive trophoblast from Day 34 horse conceptuses survives in extrauterine sites in allogeneic recipients that are immunologically naive to donor major histocompatibility complex class I antigens. The ectopic trophoblast retains its in utero characteristics, including similar lifespan, physiologic effect of its secreted product (equine chorionic gonadotropin) upon the recipient's ovaries, and induction of host immune responses. Immunologic memory has not been considered previously in this experimental system. We hypothesized that primary exposure to ectopic trophoblast would affect the recipient's immune status such that the survival time of subsequent transplants would be altered. Secondary transplant lifespans could be shortened by destructive memory responses, as has been observed in ectopic trophoblast studies in rodents, or lengthened, as occurs when male skin grafts follow multiple syngeneic pregnancies in mice. Eight mares received two closely spaced trophoblast transplants. Both grafts for each recipient were obtained from conceptuses sired by the same stallion to provide consistency in histocompatibility antigen exposure. Donor stallions were major histocompatibility complex class I homozygotes. Cytotoxic antibody production was tracked to monitor recipients' immune responses to the transplants. Detection of serum equine chorionic gonadotropin was used as a proxy for transplant lifespan. There was no significant difference between the distributions of primary and secondary transplant lifespans, despite evidence of immunologic memory. These data demonstrate that secondary ectopic trophoblast transplants in horses do not experience earlier destruction or prolonged survival following immune priming of recipients. Mechanisms responsible for the eventual demise of the transplants remain unperturbed by secondary immune responses or chronic antigenic exposure.

Speciation with gene flow in equids despite extensive chromosomal plasticity.

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.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse.

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.

Paternally expressed genes predominate in the placenta.

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.

Random X inactivation in the mule and horse placenta.

In eutherian mammals, dosage compensation of X-linked genes is achieved by X chromosome inactivation. X inactivation is random in embryonic and adult tissues, but imprinted X inactivation (paternal X silencing) has been identified in the extra-embryonic membranes of the mouse, rat, and cow. Few other species have been studied for this trait, and the data from studies of the human placenta have been discordant or inconclusive. Here, we quantify X inactivation using RNA sequencing of placental tissue from reciprocal hybrids of horse and donkey (mule and hinny). In placental tissue from the equid hybrids and the horse parent, the allelic expression pattern was consistent with random X inactivation, and imprinted X inactivation can clearly be excluded. We characterized horse and donkey XIST gene and demonstrated that XIST allelic expression in female hybrid placental and fetal tissues is negatively correlated with the other X-linked genes chromosome-wide, which is consistent with the XIST-mediated mechanism of X inactivation discovered previously in mice. As the most structurally and morphologically diverse organ in mammals, the placenta also appears to show diverse mechanisms for dosage compensation that may result in differences in conceptus development across species.

IL-22 is expressed by the invasive trophoblast of the equine (Equus caballus) chorionic girdle.

The invasive trophoblast cells of the equine placenta migrate into the endometrium to form endometrial cups, dense accumulations of trophoblast cells that produce equine chorionic gonadotropin between days 40 and 120 of normal pregnancy. The mechanisms by which the trophoblast cells invade the endometrium while evading maternal immune destruction are poorly defined. A gene expression microarray analysis performed on placental tissues obtained at day 34 of gestation revealed a >900-fold upregulation of mRNA encoding the cytokine IL-22 in chorionic girdle relative to noninvasive chorion. Quantitative RT-PCR assays were used to verify high expression of IL-22 in chorionic girdle. Additional quantitative RT-PCR analysis showed a striking increase in IL-22 mRNA expression in chorionic girdle from days 32 to 35 and an absence of IL-22 expression in other conceptus tissues. Bioinformatic analysis and cDNA sequencing confirmed the predicted length of horse IL-22, which carries a 3' extension absent in IL-22 genes of humans and mice, but present in the cow and pig. Our discovery of IL-22 in the chorionic girdle is a novel finding, as this cytokine has been previously reported in immune cells only. IL-22 has immunoregulatory functions, with primary action on epithelial cells. mRNA of IL-22R1 was detected in pregnant endometrium at levels similar to other equine epithelia. Based upon these findings, we hypothesize that IL-22 cytokine produced by the chorionic girdle binds IL-22R1 on endometrium, serving as a mechanism of fetal-maternal communication by modulating endometrial responses to trophoblast invasion.

Analysis of MHC class I genes across horse MHC haplotypes.

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.

The complete map of the Ig heavy chain constant gene region reveals evidence for seven IgG isotypes and for IgD in the horse.

This report contains the first map of the complete Ig H chain constant (IGHC) gene region of the horse (Equus caballus), represented by 34 overlapping clones from a new bacterial artificial chromosome library. The different bacterial artificial chromosome inserts containing IGHC genes were identified and arranged by hybridization using overgo probes specific for individual equine IGHC genes. The analysis of these IGHC clones identified two previously undetected IGHC genes of the horse. The newly found IGHG7 gene, which has a high homology to the equine IGHG4 gene, is located between the IGHG3 and IGHG4 genes. The high degree of conservation shared between the nucleotide sequences of the IGHG7 and IGHG4 genes is unusual for the IGHG genes of the horse and suggests that these two genes duplicated most recently during evolution of the equine IGHG genes. Second, we present the genomic nucleotide sequence of the equine IGHD gene, which is located downstream of the IGHM gene. Both the IGHG7 and IGHD genes were found to be expressed at the mRNA level. The order of the 11 IGHC genes in the IGH-locus of the horse was determined to be 5'-M-D-G1-G2-G3-G7-G4-G6-G5-E-A-3', confirming previous studies using lambda phage clones, with the exception that the IGHG5 gene was found to be the most downstream-located IGHG gene. Fluorescence in situ hybridization was used to localize the IGHC region to Equus caballus (ECA) 24qter, the horse chromosome corresponding to human chromosome 14, where the human IGH locus is found.