Trio-binning of a hinny refines the comparative organization of the horse and donkey X chromosomes and reveals novel species-specific features.

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids.

Production of Mare Chorionic Girdle Organoids That Secrete Equine Chorionic Gonadotropin.

The equine chorionic girdle is comprised of specialized invasive trophoblast cells that begin formation approximately 25 days after ovulation (day 0) and invade the endometrium to become endometrial cups. These specialized trophoblast cells transition from uninucleate to differentiated binucleate trophoblast cells that secrete the glycoprotein hormone equine chorionic gonadotropin (eCG; formerly known as pregnant mare serum gonadotropin or PMSG). This eCG has LH-like activity in the horse but variable LH- and FSH-like activity in other species and has been utilized for these properties both in vivo and in vitro. To produce eCG commercially, large volumes of whole blood must be collected from pregnant mares, which negatively impacts equine welfare due to repeated blood collections and the birth of an unwanted foal. Attempts to produce eCG in vitro using long-term culture of chorionic girdle explants have not been successful beyond 180 days, with peak eCG production at 30 days of culture. Organoids are three-dimensional cell clusters that self-organize and can remain genetically and phenotypically stable throughout long-term culture (i.e., months). Human trophoblast organoids have been reported to successfully produce human chorionic gonadotropin (hCG) and proliferate long-term (>1 year). The objective of this study was to evaluate whether organoids derived from equine chorionic girdle maintain physiological functionality. Here we show generation of chorionic girdle organoids for the first time and demonstrate in vitro production of eCG for up to 6 weeks in culture. Therefore, equine chorionic girdle organoids provide a physiologically representative 3D in vitro model for chorionic girdle development of early equine pregnancy.

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.

Placentation in Equids.

This chapter focuses on the early stages of placental development in horses and their relatives in the genus Equus and highlights unique features of equid reproductive biology. The equine placenta is classified as a noninvasive, epitheliochorial type. However, equids have evolved a minor component of invasive trophoblast, the chorionic girdle and endometrial cups, which links the equine placenta with the highly invasive hemochorial placentae of rodents and, particularly, with the primate placenta. Two types of fetus-to-mother signaling in equine pregnancy are mediated by the invasive equine trophoblast cells. First, endocrinological signaling mediated by equine chorionic gonadotrophin (eCG) drives maternal progesterone production to support the equine conceptus between days 40 and 100 of gestation. Only in primates and equids does the placenta produce a gonadotrophin, but the evolutionary paths taken by these two groups of mammals to produce this placental signal were very different. Second, florid expression of paternal major histocompatibility complex (MHC) class I molecules by invading chorionic girdle cells stimulates strong maternal anti-fetal antibody responses that may play a role in the development of immunological tolerance that protects the conceptus from destruction by the maternal immune system. In humans, invasive extravillous trophoblasts also express MHC class I molecules, but the loci involved, and their likely function, are different from those of the horse. Comparison of the cellular and molecular events in these disparate species provides outstanding examples of convergent evolution and co-option in mammalian pregnancy and highlights how studies of the equine placenta have produced new insights into reproductive strategies.

Cross-matching of allogeneic mesenchymal stromal cells eliminates recipient immune targeting.

Allogeneic mesenchymal stromal cells (MSCs) have been used clinically for decades, without cross-matching, on the assumption that they are immune-privileged. In the equine model, we demonstrate innate and adaptive immune responses after repeated intra-articular injection with major histocompatibility complex (MHC) mismatched allogeneic MSCs, but not MHC matched allogeneic or autologous MSCs. We document increased peri-articular edema and synovial effusion, increased synovial cytokine and chemokine concentrations, and development of donor-specific antibodies in mismatched recipients compared with recipients receiving matched allogeneic or autologous MSCs. Importantly, in matched allogeneic and autologous recipients, but not mismatched allogeneic recipients, there was increased stromal derived factor-1 along with increased MSC concentrations in synovial fluid. Until immune recognition of MSCs can be avoided, repeated clinical use of MSCs should be limited to autologous or cross-matched allogeneic MSCs. When non-cross-matched allogeneic MSCs are used in single MSC dose applications, presensitization against donor MHC should be assessed.

Genetics of Immune Disease in the Horse.

Host defenses against infection by viruses, bacteria, fungi, and parasites are critical to survival. It has been estimated that upwards of 7% of the coding genes of mammals function in immunity and inflammation. This high level of genomic investment in defense has resulted in an immune system characterized by extraordinary complexity and many levels of redundancy. Because so many genes are involved with immunity, there are many opportunities for mutations to arise that have negative effects. However, redundancy in the mammalian defense system and the adaptive nature of key immune mechanisms buffer the untoward outcomes of many such deleterious mutations.

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.

An Independent Locus Upstream of ASIP Controls Variation in the Shade of the Bay Coat Colour in Horses.

Novel coat colour phenotypes often emerge during domestication, and there is strong evidence of genetic selection for the two main genes that control base coat colour in horses-ASIP and MC1R. These genes direct the type of pigment produced, red pheomelanin (MC1R) or black eumelanin (ASIP), as well as the relative concentration and the temporal-spatial distribution of melanin pigment deposits in the skin and hair coat. Here, we describe a genome-wide association study (GWAS) to identify novel genic regions involved in the determination of the shade of bay. In total, 126 horses from five different breeds were ranked according to the extent of the distribution of eumelanin: spanning variation in phenotype from black colour restricted only to the extremities to the presence of some black pigment across nearly all the body surface. We identified a single region associated with the shade of bay ranking spanning approximately 0.5 MB on ECA22, just upstream of the ASIP gene (p = 9.76 × 10-15). This candidate region encompasses the distal 5' end of the ASIP transcript (as predicted from other species) as well as the RALY gene. Both loci are viable candidates based on the presence of similar alleles in other species. These results contribute to the growing understanding of coat colour genetics in the horse and to the mapping of genetic determinants of pigmentation on a molecular level. Given pleiotropic phenotypes in behaviour and obesity for ASIP alleles, especially those in the 5' regulatory region, improved understanding of this new Shade allele may have implications for health management in the horse.

Immunological memory and tolerance at the maternal-fetal interface: Implications for reproductive management of mares.

The development of placentation that coincided with the evolution of mammals presented new challenges to the transmission of life from one generation to the next, particularly with regard to the possibility of maternal immunological recognition and destruction of the developing conceptus. The balance between immunity and tolerance dominates the immunological relationship between mother and fetus during mammalian pregnancy, and the focal point of this relationship lies at the interface between the trophoblast cells that comprise the outermost layer of the placenta and the maternal endometrial tissues. Immune memory and tolerance are two of the cardinal characteristics of the immune system. Immune memory is essential in preventing or lessening the effect of infections to the mother or conceptus, but may also be a threat to the semi-allogeneic tissues of the fetus and placenta. The mother must develop functional immune tolerance to her fetus, but at the same time retain her ability to combat infections while pregnant. To address this imperative, mammals have developed overlapping and independent mechanisms for evading maternal anti-fetal immune responses that could result in pregnancy loss. Studies of the unusual component of equine invasive trophoblast in the epitheliochorial placenta have illuminated aspects of immune memory and tolerance that have relevance to fertility in the horse and other mammalian species.

Author Correction: Evolution of placental invasion and cancer metastasis are causally linked.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Evolution of placental invasion and cancer metastasis are causally linked.

Among mammals, placental invasion is correlated with vulnerability to malignancy. Animals with more invasive placentation (for example, humans) are more vulnerable to malignancy. To explain this correlation, we propose the hypothesis of 'Evolved Levels of Invasibility' proposing that the evolution of invasibility of stromal tissue affects both placental and cancer invasion. We provide evidence for this using an in vitro model. We find that bovine endometrial and skin fibroblasts are more resistant to invasion than are their human counterparts. Gene expression profiling identified genes with high expression in human but not in bovine fibroblasts. Knocking down a subset of them in human fibroblasts leads to stronger resistance to cancer cell invasion. Identifying the evolutionary determinants of stromal invasibility can provide important insights to develop rational antimetastatic therapeutics.

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

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites.

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

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.

In vitro MSC function is related to clinical reaction in vivo.

BACKGROUND: We recently demonstrated that intracellular xenogen-contaminated autologous MSCs (FBS) and non-xenogen-contaminated allogeneic (ALLO) MSCs caused an adverse clinical response after repeated intra-articular injection in horses, whereas autologous (AUTO) MSCs did not. Our current objective was to use clinical data from the previous study to compare MSC stemness against adverse response indicated by synovial total nucleated cell count (TNCC) following intra-articular MSC injection. METHODS: Stemness, quantified by a trilineage differentiation (TLD) score; immunomodulation, quantified by mixed lymphocyte reactions (MLRs); and degree of MHCI expression, quantified by mean fluorescent intensity (MFI); were correlated to the synovial TNCC 24 h after naïve and primed injection. RESULTS: There was a trend of a negative correlation (p = 0.21, r = - 0.44) between TLD score and TNCC after primed injection in the ALLO group. Within the ALLO group only, there was a significant positive correlation (p = 0.05, r = 0.77) between MHCI MFI and TNCC after naïve injection and a trend (p = 0.16, r = 0.49) of a positive association of MHCI MFI to TNCC after primed injection. Within the FBS group only, there was a positive correlation (p = 0.04, r = 1) between TNCC and lymphocyte proliferation after both injections. CONCLUSIONS: The trend of a negative correlation of TLD score and TNCC in the ALLO, but not the FBS group, together with the association of MHCI expression and TNCC in the ALLO group, indicates that improved stemness is associated with reduced MSC immunogenicity. When inflammation was incited by xenogen, there was a strong correlation of lymphocyte activation in vitro to adverse response in vivo, confirming that MLRs in vitro reflect MSC immunomodulatory activity in vivo. The relationship of stemness in vitro, suppression of lymphocyte activation in vitro, MHCI expression in vitro, and clinical response in vivo should be further investigated.

Correction: miRNA independent hepacivirus variants suggest a strong evolutionary pressure to maintain miR-122 dependence.

[This corrects the article DOI: 10.1371/journal.ppat.1006694.].

Birth, evolution, and transmission of satellite-free mammalian centromeric domains.

Mammalian centromeres are associated with highly repetitive DNA (satellite DNA), which has so far hindered molecular analysis of this chromatin domain. Centromeres are epigenetically specified, and binding of the CENPA protein is their main determinant. In previous work, we described the first example of a natural satellite-free centromere on Equus caballus Chromosome 11. Here, we investigated the satellite-free centromeres of Equus asinus by using ChIP-seq with anti-CENPA antibodies. We identified an extraordinarily high number of centromeres lacking satellite DNA (16 of 31). All of them lay in LINE- and AT-rich regions. A subset of these centromeres is associated with DNA amplification. The location of CENPA binding domains can vary in different individuals, giving rise to epialleles. The analysis of epiallele transmission in hybrids (three mules and one hinny) showed that centromeric domains are inherited as Mendelian traits, but their position can slide in one generation. Conversely, centromere location is stable during mitotic propagation of cultured cells. Our results demonstrate that the presence of more than half of centromeres void of satellite DNA is compatible with genome stability and species survival. The presence of amplified DNA at some centromeres suggests that these arrays may represent an intermediate stage toward satellite DNA formation during evolution. The fact that CENPA binding domains can move within relatively restricted regions (a few hundred kilobases) suggests that the centromeric function is physically limited by epigenetic boundaries.