Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series.

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

Widespread horse-based mobility arose around 2200 BCE in Eurasia.

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and their widespread integration as a means of transport remains contentious2-4. Here we assemble a collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged around 2200 BCE, through close-kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than approximately 2700 BCE, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe around 3000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai around 3500 BCE, a settlement from central Asia associated with corrals and a subsistence economy centred on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

The Evolutionary and Historical Foundation of the Modern Horse: Lessons from Ancient Genomics.

The domestication of the horse some 5,500 years ago followed those of dogs, sheep, goats, cattle, and pigs by ∼2,500-10,000 years. By providing fast transportation and transforming warfare, the horse had an impact on human history with no equivalent in the animal kingdom. Even though the equine sport industry has considerable economic value today, the evolutionary history underlying the emergence of the modern domestic horse remains contentious. In the last decade, novel sequencing technologies have revolutionized our capacity to sequence the complete genome of organisms, including from archaeological remains. Applied to horses, these technologies have provided unprecedented levels of information and have considerably changed models of horse domestication. This review illustrates how ancient DNA, especially ancient genomes, has inspired researchers to rethink the process by which horses were first domesticated and then diversified into a variety of breeds showing a range of traits that are useful to humans.

The origins and spread of domestic horses from the Western Eurasian steppes.

Domestication of horses fundamentally transformed long-range mobility and warfare1. However, modern domesticated breeds do not descend from the earliest domestic horse lineage associated with archaeological evidence of bridling, milking and corralling2-4 at Botai, Central Asia around 3500 BC3. Other longstanding candidate regions for horse domestication, such as Iberia5 and Anatolia6, have also recently been challenged. Thus, the genetic, geographic and temporal origins of modern domestic horses have remained unknown. Here we pinpoint the Western Eurasian steppes, especially the lower Volga-Don region, as the homeland of modern domestic horses. Furthermore, we map the population changes accompanying domestication from 273 ancient horse genomes. This reveals that modern domestic horses ultimately replaced almost all other local populations as they expanded rapidly across Eurasia from about 2000 BC, synchronously with equestrian material culture, including Sintashta spoke-wheeled chariots. We find that equestrianism involved strong selection for critical locomotor and behavioural adaptations at the GSDMC and ZFPM1 genes. Our results reject the commonly held association7 between horseback riding and the massive expansion of Yamnaya steppe pastoralists into Europe around 3000 BC8,9 driving the spread of Indo-European languages10. This contrasts with the scenario in Asia where Indo-Iranian languages, chariots and horses spread together, following the early second millennium BC Sintashta culture11,12.

An endothelial regulatory module links blood pressure regulation with elite athletic performance.

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

Ancient DNA reveals an early adoption of horse culture by Native Americans.

Native Americans of the Plains and Rocky Mountains are renowned for their horsemanship. Taylor et al. recently used ancient DNA and other bioarcheological approaches to document how horses dispersed throughout America and transformed Native American societies following their introduction by the Spanish in 1519, well before the arrival of European settlers.

Functional annotation of the animal genomes: An integrated annotation resource for the horse.

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Complex N-glycans are important for interspecies transmission of H7 influenza A viruses.

Influenza A viruses (IAVs) can overcome species barriers by adaptation of the receptor-binding site of the hemagglutinin (HA). To initiate infection, HAs bind to glycan receptors with terminal sialic acids, which are either N-acetylneuraminic acid (NeuAc) or N-glycolylneuraminic acid (NeuGc); the latter is mainly found in horses and pigs but not in birds and humans. We investigated the influence of previously identified equine NeuGc-adapting mutations (S128T, I130V, A135E, T189A, and K193R) in avian H7 IAVs in vitro and in vivo. We observed that these mutations negatively affected viral replication in chicken cells but not in duck cells and positively affected replication in horse cells. In vivo, the mutations reduced virus virulence and mortality in chickens. Ducks excreted high viral loads longer than chickens, although they appeared clinically healthy. To elucidate why these viruses infected chickens and ducks despite the absence of NeuGc, we re-evaluated the receptor binding of H7 HAs using glycan microarray and flow cytometry studies. This re-evaluation demonstrated that mutated avian H7 HAs also bound to α2,3-linked NeuAc and sialyl-LewisX, which have an additional fucose moiety in their terminal epitope, explaining why infection of ducks and chickens was possible. Interestingly, the α2,3-linked NeuAc and sialyl-LewisX epitopes were only bound when presented on tri-antennary N-glycans, emphasizing the importance of investigating the fine receptor specificities of IAVs. In conclusion, the binding of NeuGc-adapted H7 IAV to tri-antennary N-glycans enables viral replication and shedding by chickens and ducks, potentially facilitating interspecies transmission of equine-adapted H7 IAVs.IMPORTANCEInfluenza A viruses (IAVs) cause millions of deaths and illnesses in birds and mammals each year. The viral surface protein hemagglutinin initiates infection by binding to host cell terminal sialic acids. Hemagglutinin adaptations affect the binding affinity to these sialic acids and the potential host species targeted. While avian and human IAVs tend to bind to N-acetylneuraminic acid (sialic acid), equine H7 viruses prefer binding to N-glycolylneuraminic acid (NeuGc). To better understand the function of NeuGc-specific adaptations in hemagglutinin and to elucidate interspecies transmission potential NeuGc-adapted viruses, we evaluated the effects of NeuGc-specific mutations in avian H7 viruses in chickens and ducks, important economic hosts and reservoir birds, respectively. We also examined the impact on viral replication and found a binding affinity to tri-antennary N-glycans containing different terminal epitopes. These findings are significant as they contribute to the understanding of the role of receptor binding in avian influenza infection.

A method to generate capture baits for targeted sequencing.

Hybridization capture approaches allow targeted high-throughput sequencing analysis at reduced costs compared to shotgun sequencing. Hybridization capture is particularly useful in analyses of genomic data from ancient, environmental, and forensic samples, where target content is low, DNA is fragmented and multiplex PCR or other targeted approaches often fail. Here, we describe a DNA bait synthesis approach for hybridization capture that we call Circular Nucleic acid Enrichment Reagent, or CNER (pronounced 'snare'). The CNER method uses rolling-circle amplification followed by restriction digestion to discretize microgram quantities of hybridization probes. We demonstrate the utility of the CNER method by generating probes for a panel of 23 771 known sites of single nucleotide polymorphism in the horse genome. Using these probes, we capture and sequence from a panel of ten ancient horse DNA libraries, comparing CNER capture efficiency to a commercially available approach. With about one million read pairs per sample, CNERs captured more targets (90.5% versus 66.5%) at greater mean depth than an alternative commercial approach.

Administration and detection of a multi-target rAAV gene doping vector in horses using multiple matrices and molecular techniques.

Gene doping, which includes the non-therapeutic use of genes or genetic elements that have the capacity to enhance athletic performance, is prohibited in horseracing and equestrian sports. To provide a comprehensive assessment of matrix and detection techniques, a custom adeno-associated virus serotype 8 vector was designed to include PCR binding sites for multiple target genes and assay types. The vector was injected via an intramuscular route into two Thoroughbred horses and matrices collected at defined timepoints. DNA was analysed using 3 detection methods: qPCR, digital PCR, and NGS. Overall, there was a strong correlation across the different detection methods employed, although digital PCR was less sensitive at lower concentrations. High concentrations of vector were detected at early timepoints in plasma and whole blood, which rapidly dropped after 0.5 d to trace levels by 4 d and 9 d post-administration respectively, following a similar pattern to previous studies. Vector was detected in dried blood spots at lower levels than whole blood, but with a similar detection time. Detection in hair root bulbs in one horse was observed at over a month post-administration, which opens new avenues for future gene doping testing in humans and animals.

Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies.

BACKGROUND: The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds. RESULTS: HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20). CONCLUSIONS: Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Single nucleotide polymorphisms (SNPs) in the open reading frame (ORF) of prion protein gene (PRNP) in Nigerian livestock species.

BACKGROUND: Prion diseases, also known as transmissible spongiform encephalopathies (TSEs) remain one of the deleterious disorders, which have affected several animal species. Polymorphism of the prion protein (PRNP) gene majorly determines the susceptibility of animals to TSEs. However, only limited studies have examined the variation in PRNP gene in different Nigerian livestock species. Thus, this study aimed to identify the polymorphism of PRNP gene in Nigerian livestock species (including camel, dog, horse, goat, and sheep). We sequenced the open reading frame (ORF) of 65 camels, 31 village dogs and 12 horses from Nigeria and compared with PRNP sequences of 886 individuals retrieved from public databases. RESULTS: All the 994 individuals were assigned into 162 haplotypes. The sheep had the highest number of haplotypes (n = 54), and the camel had the lowest (n = 7). Phylogenetic tree further confirmed clustering of Nigerian individuals into their various species. We detected five non-synonymous SNPs of PRNP comprising of G9A, G10A, C11G, G12C, and T669C shared by all Nigerian livestock species and were in Hardy-Weinberg Equilibrium (HWE). The amino acid changes in these five non-synonymous SNP were all "benign" via Polyphen-2 program. Three SNPs G34C, T699C, and C738G occurred only in Nigerian dogs while C16G, G502A, G503A, and C681A in Nigerian horse. In addition, C50T was detected only in goats and sheep. CONCLUSION: Our study serves as the first to simultaneously investigate the polymorphism of PRNP gene in Nigerian livestock species and provides relevant information that could be adopted in programs targeted at breeding for prion diseases resistance.

Estimating Temporally Variable Selection Intensity from Ancient DNA Data.

Novel technologies for recovering DNA information from archaeological and historical specimens have made available an ever-increasing amount of temporally spaced genetic samples from natural populations. These genetic time series permit the direct assessment of patterns of temporal changes in allele frequencies and hold the promise of improving power for the inference of selection. Increased time resolution can further facilitate testing hypotheses regarding the drivers of past selection events such as the incidence of plant and animal domestication. However, studying past selection processes through ancient DNA (aDNA) still involves considerable obstacles such as postmortem damage, high fragmentation, low coverage, and small samples. To circumvent these challenges, we introduce a novel Bayesian framework for the inference of temporally variable selection based on genotype likelihoods instead of allele frequencies, thereby enabling us to model sample uncertainties resulting from the damage and fragmentation of aDNA molecules. Also, our approach permits the reconstruction of the underlying allele frequency trajectories of the population through time, which allows for a better understanding of the drivers of selection. We evaluate its performance through extensive simulations and demonstrate its utility with an application to the ancient horse samples genotyped at the loci for coat coloration. Our results reveal that incorporating sample uncertainties can further improve the inference of selection.

Unveiling the population genetic structure of Iranian horses breeds by whole-genome resequencing analysis.

Preserving genetic diversity is pivotal for enhancing genetic improvement and facilitating adaptive responses to selection. This study focuses on identifying key genetic variants, including single nucleotide polymorphisms (SNPs), insertion/deletion polymorphisms (INDELs), and copy number variants (CNVs), while exploring the genomic evolutionary connectedness among seven Iranian horses representing five indigenous breeds: Caspian, Turkemen, DareShuri, Kurdish, and Asil. Using whole-genome resequencing, we generated 2.7 Gb of sequence data, with raw reads ranging from 1.2 Gb for Caspian horses to 0.38 Gb for Turkoman horses. Post-filtering, approximately 1.9 Gb of reads remained, with ~ 1.5 Gb successfully mapped to the horse reference genome (EquCab3.0), achieving mapping rates between 76.4% (Caspian) and 98.35% (Turkoman). We identified 2,909,816 SNPs in Caspian horses, constituting around 0.1% of the genome. Notably, 71% of these SNPs were situated in intergenic regions, while 8.5 and 6.8% were located upstream and downstream, respectively. A comparative analysis of SNPs between Iranian and non-Iranian horse breeds showed that Caspian horses had the lowest number of shared SNPs with Turkoman horses. Instead, they showed a closer genetic relationship with DareShuri, Quarter, Arabian, Standardbred, and Asil breeds. Hierarchical clustering highlighted Caspian horses as a distinct cluster, underscoring their distinctive genomic signature. Caspian horses exhibit a unique genetic profile marked by an enrichment of private mutations in neurological genes, influencing sensory perception and awareness. This distinct genetic makeup shapes mating preferences and signifies a separate evolutionary trajectory. Additionally, significant non-synonymous single nucleotide polymorphisms (nsSNPs) in reproductive genes offer intervention opportunities for managing Caspian horses. These findings reveal the population genetic structure of Iranian horse breeds, contributing to the advancement of knowledge in areas such as conservation, performance traits, climate adaptation, reproduction, and resistance to diseases in equine science.

Evolutionarily Young African Rhinoceros Gammaretroviruses.

High-throughput sequences were generated from DNA and cDNA from four Southern white rhinoceros (Ceratotherium simum simum) located in the Taronga Western Plain Zoo in Australia. Virome analysis identified reads that were similar to Mus caroli endogenous gammaretrovirus (McERV). Previous analysis of perissodactyl genomes did not recover gammaretroviruses. Our analysis, including the screening of the updated white rhinoceros (Ceratotherium simum) and black rhinoceros (Diceros bicornis) draft genomes identified high-copy orthologous gammaretroviral ERVs. Screening of Asian rhinoceros, extinct rhinoceros, domestic horse, and tapir genomes did not identify related gammaretroviral sequences in these species. The newly identified proviral sequences were designated SimumERV and DicerosERV for the white and black rhinoceros retroviruses, respectively. Two long terminal repeat (LTR) variants (LTR-A and LTR-B) were identified in the black rhinoceros, with different copy numbers associated with each (n = 101 and 373, respectively). Only the LTR-A lineage (n = 467) was found in the white rhinoceros. The African and Asian rhinoceros lineages diverged approximately 16 million years ago. Divergence age estimation of the identified proviruses suggests that the exogenous retroviral ancestor of the African rhinoceros ERVs colonized their genomes within the last 8 million years, a result consistent with the absence of these gammaretroviruses from Asian rhinoceros and other perissodactyls. The black rhinoceros germ line was colonized by two lineages of closely related retroviruses and white rhinoceros by one. Phylogenetic analysis indicates a close evolutionary relationship with ERVs of rodents including sympatric African rats, suggesting a possible African origin of the identified rhinoceros gammaretroviruses. IMPORTANCE Rhinoceros genomes were thought to be devoid of gammaretroviruses, as has been determined for other perissodactyls (horses, tapirs, and rhinoceros). While this may be true of most rhinoceros, the African white and black rhinoceros genomes have been colonized by evolutionarily young gammaretroviruses (SimumERV and DicerosERV for the white and black rhinoceros, respectively). These high-copy endogenous retroviruses (ERVs) may have expanded in multiple waves. The closest relative of SimumERV and DicerosERV is found in rodents, including African endemic species. Restriction of the ERVs to African rhinoceros suggests an African origin for the rhinoceros gammaretroviruses.

Long-term persistence of viral RNA and inflammation in the CNS of macaques exposed to aerosolized Venezuelan equine encephalitis virus.

Venezuelan equine encephalitis virus (VEEV) is a positively-stranded RNA arbovirus of the genus Alphavirus that causes encephalitis in humans. Cynomolgus macaques are a relevant model of the human disease caused by VEEV and are useful in exploring pathogenic mechanisms and the host response to VEEV infection. Macaques were exposed to small-particle aerosols containing virus derived from an infectious clone of VEEV strain INH-9813, a subtype IC strain isolated from a human infection. VEEV-exposed macaques developed a biphasic fever after infection similar to that seen in humans. Maximum temperature deviation correlated with the inhaled dose, but fever duration did not. Neurological signs, suggestive of virus penetration into the central nervous system (CNS), were predominantly seen in the second febrile period. Electroencephalography data indicated a statistically significant decrease in all power bands and circadian index during the second febrile period that returned to normal after fever resolved. Intracranial pressure increased late in the second febrile period. On day 6 post-infection macaques had high levels of MCP-1 and IP-10 chemokines in the CNS, as well as a marked increase of T lymphocytes and activated microglia. More than four weeks after infection, VEEV genomic RNA was found in the brain, cerebrospinal fluid and cervical lymph nodes. Pro-inflammatory cytokines & chemokines, infiltrating leukocytes and pathological changes were seen in the CNS tissues of macaques euthanized at these times. These data are consistent with persistence of virus replication and/or genomic RNA and potentially, inflammatory sequelae in the central nervous system after resolution of acute VEEV disease.

Non-Synonymous Substitutions in Cadherin 13, Solute Carrier Family 6 Member 4, and Monoamine Oxidase A Genes are Associated with Personality Traits in Thoroughbred Horses.

Retraining retired racehorses for various purposes can help correct behavioral issues. However, ensuring efficiency and preventing accidents present global challenges. Based on the hypothesis that a simple personality assessment could help address these challenges, the present study aimed to identify genetic markers associated with personality. Eight genes were selected from 18 personality-related candidate genes that are orthologs of human personality genes, and their association with personality was verified based on actual behavior. A total of 169 Thoroughbred horses were assessed for their tractability (questionnaire concerning tractability in 14 types of situations and 3 types of impressions) during the training process. Personality factors were extracted from the data using principal component analysis and analyzed for their association with single nucleotide variants as non-synonymous substitutions in the target genes. Three genes, CDH13, SLC6A4, and MAOA, demonstrated significant associations based on simple linear regression, marking the identification of these genes for the first time as contributors to temperament in Thoroughbred horses. All these genes, as well as the previously identified HTR1A, are involved in the serotonin neurotransmitter system, suggesting that the tractability of horses may be correlated with their social personality. Assessing the genotypes of these genes before retraining is expected to prevent problems in the development of a racehorse's second career and shorten the training period through individual customization of training methods, thereby improving racehorse welfare.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed.

BACKGROUND: The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity. RESULTS: To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders. CONCLUSIONS: The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes.

BACKGROUND: Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected. RESULTS: A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3. CONCLUSIONS: The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Intrauterine devices influence prostaglandin secretion by equine uterus: in vitro and in vivo studies.

BACKGROUND: Intrauterine devices (IUD) are used in the veterinary practice as the non-pharmacological method of oestrus suppression in mares. When placed in the uterus, IUD create a physical contact with the endometrium that mimics the presence of an equine embryo. However, the mechanism of their action has not been fully elucidated. The objective of the present study was to examine the effect of mechanical stimulation of IUD on mare`s endometrium in both in vitro and in vivo study. For this purpose, we demonstrated the effect of IUD on prostaglandin (PG) F2α and PGE2 secretion, and mRNA transcription of genes involved in PG synthesis pathway in equine endometrial cells in vitro. In the in vivo study, we aimed to compare short-term effect of IUD inserted on day 0 (oestrus) with day 5-6 post-ovulation (the specific time when embryo reaches uterus after fertilization) on PG secretion from equine endometrium. To determine the long-term effect on PG synthase mRNA transcription, a single endometrial biopsy was taken only once within each group of mares at certain time points of the estrous cycle from mares placement with IUD on days 0 or 5-6 post-ovualtion. RESULTS: We showed for the first time that the incubation of the endometrial cells with the presence of IUD altered the pattern of PG synthase mRNA transcription in equine epithelial and stromal endometrial cells. In vivo, in mares placement with IUD on day 0, PGE2 concentrations in blood plasma were upregulated between 1 and 6, and at 10 h after the IUD insertion, compared with the control mares (P < 0.05). Moreover, the decrease of PTGFS mRNA transcription on day 16- 18, associated with an elevation in PTGES mRNA transcription on day 20 -21 of the estrous cycle in endometrial biopsies collected from mares placement with IUD on days 5-6 suggest an antiluteolytic action of IUD during the estrous cycle. CONCLUSION: We conclude that the application of IUD may mimic the equine conceptus presence through the physical contact with the endometrium altering PG synthase transcription, and act as a potent modulator of endometrial PG secretion both in vitro and in vivo.