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.

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.

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.

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.

Identification of neuropathogenic Varicellovirus equidalpha1 as a potential cause of respiratory disease outbreaks among horses in North Xinjiang, China, from 2021-2023.

BACKGROUND: Varicellovirus equidalpha1 (formerly Equid alphaherpesvirus 1, EqAHV-1) is among the most important viruses responsible for respiratory disease outbreaks among horses throughout the world. No reports to date have detailed the association between EqAHV-1 and respiratory disease among horses in China. This study described one such outbreak among a population of horses in north Xinjiang that occurred from April 2021 - May 2023. RESULTS: qPCR revealed that EqAHV-1 was detectable in all samples and this virus was identified as a possible source of respiratory disease, although a limited subset of these samples were also positive for EqAHV-2, EqAHV-4, and EqAHV-5. In total, three EqAHV-1 strains responsible for causing respiratory illness in horses were isolated successfully, and full-length ORF33 sequence comparisonsand phylogenetic analyses indicated that these isolates may have originated from EqAHV-1 strains detected in Yili horse abortions. ORF30 sequence data additionally suggested that these strains were neuropathic, as evidenced by the presence of a guanine residue at nucleotide position 2254 corresponding to the aspartic acid present at position 752 in the DNA polymerase encoded by this virus. CONCLUSION: This study is the first report of an outbreak of respiratory disease among horses in China caused by EqAHV-1. ORF30 sequence characterization revealed that these EqAHV-1 strains harbored a neuropathogenic genotype. Given the detection of this virus in horses suffering from respiratory disease, concern is warranted with respect to this neuropathogenic EqAHV-1 outbreak.

Study on the origin of the Baise horse based on whole-genome resequencing.

The Baise horse, an indigenous horse breed mainly distributed in the Baise region of Guangxi province in southwest China, has a long history as draft animal. However, there is a lack of research regarding the origin and ancestral composition of the Baise horse. In this study, whole-genome resequencing data from 236 horses of seven Chinese indigenous horse breeds, five foreign horse breeds, and four Przewalski's horses were used to investigate the relationships between the Baise horse and other horse breeds. The results showed that foreign horse breeds had no significant impact on the formation of the Baise horse. The two southwestern horse populations, the Debao pony and the Jinjiang horse, exhibit the closest genetic affinity with the Baise horse. This is consistent with their adjacent geographical distribution. Analysis of the migration route revealed a gene flow from the Chakouyi horse into the Baise horse. In summary, our results confirm that the formation of the Baise horse did not involve participation from foreign breeds. Geographical distance emerges as a crucial factor in determining the genetic relationships with the Baise horse. Gene flows of indigenous horse breeds along ancient routes of trade activities had played a role in the formation of the Baise horse.

Population structure and genomic diversity of the Einsiedler horse.

The breeding history of the Einsiedler horse is closely connected with the Benedictine cloister Einsiedeln. In the mid-nineteenth century, it was decided to use European Warmblood stallions for cross-breeding and to abandon the selection of stallions. Since that time, it has only been possible to trace back the origin of Einsiedler horses using maternal ancestry information. Here, we collected high-density genotype data for European Warmblood horses (Selle Français, Swiss Warmblood and Einsiedler) and Franches-Montagnes horses, the last native Swiss horse breed, to unravel the current population structure of the Einsiedler horse. Using commonly applied methods to ascertain fine-scale population structures, it was not possible to clearly differentiate the Einsiedler from other European Warmblood horses. However, by means of runs of homozygosity (ROH) we were able to detect breed-specific ROH islands for the Einsiedler horse, including genes involved in domestication and adaptation to high altitude. Therefore, future breeding activities should involve the screening of these breed-specific ROH segments, the revival of cryopreserved sperm and the selection of Einsiedler stallions.

Evaluation of the effectiveness of single nucleotide polymorphisms compared to microsatellite markers for parentage verification in Moroccan horses.

The International Society for Animal Genetics (ISAG) currently advocates for a transition towards single nucleotide polymorphism (SNP) markers as a potential alternative for equine parentage verification. To ascertain the efficacy of this transition, it is imperative to evaluate the performance of parentage testing using SNPs in juxtaposition with short tandem repeats (STRs). As per ISAG's recommendation, we used an equine genotyping-by-sequencing panel with 144 SNPs for this purpose. Equine parentage is currently realized using 16 microsatellites (STRs) excluding the LEX3 marker. In this study, 1074 horses were genotyped using the 144 SNPs panel, including 432 foals, 414 mares, and 228 stallions, from five different breeds: 293 Arabians, 167 Barbs, 189 Thoroughbreds, 73 Anglo-Arabians, and 352 Arabian-Barbs. As a result, two SNPs markers were eliminated from the panel system due to inconsistent amplification across all examined individuals leaving 142 SNPs markers for analysis. A comparative analysis between SNPs and STRs markers revealed that the mean expected heterozygosity was 0.457 for SNPs and 0.76 for STRs, while the mean observed heterozygosity stood at 0.472 for SNPs and 0.72 for STRs. Furthermore, the probability of identity was calculated to be 5.722 × 10-57 for SNPs and 1.25 × 10-15 for STRs markers. In alignment with the Hardy-Weinberg equilibrium in polyploids test, 110 out of the total SNPs were consistent with the Hardy-Weinberg equilibrium in polyploids test (p > 0.05). Employing both SNPs and STRs markers, the mean polymorphic information content was discerned to be 0.351 for SNPs and 0.72 for STRs. The cumulative exclusion probabilities for SNP markers exceeded 99.99%, indicating that the 142 SNPs panel might be adequate for parentage testing. In contrast, when utilizing STRs markers, the combined average exclusion probabilities for one and both parents were determined to be 99.8% and 99.9%, respectively. Our comprehensive study underscores the potential of SNPs in equine parentage verification, especially when compared to STRs in terms of exclusion probabilities. As a corollary, the application of SNPs for parentage verification and identification can significantly contribute to the conservation initiative for the five Moroccan horse breeds. Nonetheless, further research is required to address and replace the deficient SNPs within the panel.