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

Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

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.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus)

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

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.

Ancient Sheep Genomes Reveal Four Millennia of North European Short-Tailed Sheep in the Baltic Sea Region

Sheep are among the earliest domesticated livestock species, with a wide variety of breeds present today. However, it remains unclear how far back this diversity goes, with formal documentation only dating back a few centuries. North European short-tailed (NEST) breeds are often assumed to be among the oldest domestic sheep populations, even thought to represent relicts of the earliest sheep expansions during the Neolithic period reaching Scandinavia <6,000 years ago. This study sequenced the genomes (up to 11.6X) of five sheep remains from the Baltic islands of Gotland and Åland, dating from the Late Neolithic (∼4,100 cal BP) to historical times (∼1,600 CE). Our findings indicate that these ancient sheep largely possessed the genetic characteristics of modern NEST breeds, suggesting a substantial degree of long-term continuity of this sheep type in the Baltic Sea region. Despite the wide temporal spread, population genetic analyses show high levels of affinity between the ancient genomes and they also exhibit relatively high genetic diversity when compared to modern NEST breeds, implying a loss of diversity in most breeds during the last centuries associated with breed formation and recent bottlenecks. Our results shed light on the development of breeds in Northern Europe specifically as well as the development of genetic diversity in sheep breeds, and their expansion from the domestication center in general.

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.

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.

Advancing responsible genomic analyses of ancient mollusc shells

The analysis of the DNA entrapped in ancient shells of molluscs has the potential to shed light on the evolution and ecology of this very diverse phylum. Ancient genomics could help reconstruct the responses of molluscs to past climate change, pollution, and human subsistence practices at unprecedented temporal resolutions. Applications are however still in their infancy, partly due to our limited knowledge of DNA preservation in calcium carbonate shells and the need for optimized methods for responsible genomic data generation. To improve ancient shell genomic analyses, we applied high-throughput DNA sequencing to 27 Mytilus mussel shells dated to ~111-6500 years Before Present, and investigated the impact, on DNA recovery, of shell imaging, DNA extraction protocols and shell sub-sampling strategies. First, we detected no quantitative or qualitative deleterious effect of micro-computed tomography for recording shell 3D morphological information prior to sub-sampling. Then, we showed that double-digestion and bleach treatment of shell powder prior to silica-based DNA extraction improves shell DNA recovery, also suggesting that DNA is protected in preservation niches within ancient shells. Finally, all layers that compose Mytilus shells, i.e., the nacreous (aragonite) and prismatic (calcite) carbonate layers, with or without the outer organic layer (periostracum) proved to be valuable DNA reservoirs, with aragonite appearing as the best substrate for genomic analyses. Our work contributes to the understanding of long-term molecular preservation in biominerals and we anticipate that resulting recommendations will be helpful for future efficient and responsible genomic analyses of ancient mollusc shells.

Eimeria leuckarti in equid coprolites from the Sassanid Era (2nd-6th century CE) excavated in Chehrabad Salt Mine archaeological site, Iran

OBJECTIVE

This study reports coccidian oocysts in an equid coprolite dated to the Sassanid Empire (2nd-6th century CE) recovered in Chehrabad Salt Mine archaeological site, Iran.

METHODS

Between 2015 and 2017, an archaeoparasitological investigation led to the discovery of an equid coprolite in the Chehrabad Salt Mine archeological site, (Douzlakh), western Iran. Samples were rehydrated using trisodium phosphate solution and were examined by light microscopy.

RESULTS

Seven oocysts of Eimeria leuckarti (Flesch, 1883) were identified; they were in various stages of sporulation.

CONCLUSION

This is the first report of ancient coccidian oocysts from equids. The importance of this observation is discussed, and current knowledge of eimeriid oocysts at archaeological sites is reviewed.

SIGNIFICANCE

The observations of E. leuckarti increases current knowledge of parasite biodiversity in ancient Iran when it rested along the Silk Road, a network of trade routes connecting the East and West that was central to economic, cultural, political, and religious interactions between these regions, and to livestock movement that could contribute to the transmission of the parasites from/to other regions.

LIMITATIONS

The contextual information about animal species present in and around the Salt Mine during its working periods, including Achaemenid dynasty (6th to 4th century BCE) and Sassanid era (2nd to 6th century CE), is very limited and does not allow secure conclusions regarding the host origin of the coprolites.

SUGGESTIONS FOR FURTHER RESEARCH

Application of molecular biology tools to identify the correct host origin of the coprolites and to detect more parasite species is advocated.

Genome-wide equine preimplantation genetic testing enabled by simultaneous haplotyping and copy number detection

In different species, embryonic aneuploidies and genome-wide errors are a major cause of developmental failure. The increasing number of equine embryos being produced worldwide provides the opportunity to characterize and rank or select embryos based on their genetic profile prior to transfer. Here, we explored the possibility of generic, genome-wide preimplantation genetic testing concurrently for aneuploidies (PGT-A) and monogenic (PGT-M) traits and diseases in the horse, meanwhile assessing the incidence and spectrum of chromosomal and genome-wide errors in in vitro-produced equine embryos. To this end, over 70,000 single nucleotide polymorphism (SNP) positions were genotyped in 14 trophectoderm biopsies and corresponding biopsied blastocysts, and in 26 individual blastomeres from six arrested cleavage-stage embryos. Subsequently, concurrent genome-wide copy number detection and haplotyping by haplarithmisis was performed and the presence of aneuploidies and genome-wide errors and the inherited parental haplotypes for four common disease-associated genes with high carrier frequency in different horse breeds (GBE1, PLOD1, B3GALNT2, MUTYH), and for one color coat-associated gene (STX17) were compared in biopsy-blastocyst combinations. The euploid (n = 12) or fully aneuploid (n = 2) state and the inherited parental haplotypes for 42/45 loci of interest of the biopsied blastocysts were predicted by the biopsy samples in all successfully analyzed biopsy-blastocyst combinations (n = 9). Two biopsies showed a loss of maternal chromosome 28 and 31, respectively, which were confirmed in the corresponding blastocysts. In one of those biopsies, additional complex aneuploidies not present in the blastocyst were found. Five out of six arrested embryos contained chromosomal and/or genome-wide errors in most of their blastomeres, demonstrating their contribution to equine embryonic arrest in vitro. The application of the described PGT strategy would allow to select equine embryos devoid of genetic errors and pathogenetic variants, and with the variants of interest, which will improve foaling rate and horse quality. We believe this approach will be a gamechanger in horse breeding.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

Mitochondrial DNA sequencing illuminates genetic diversity and origin of Hunagrian Nonius horse breed and his relatives - Danubian horse and Serbian Nonius

From a historical perspective, horse breeding in Bulgaria has been very well developed since the time of the Thracians (early Bronze Age c. 3000 BCE). Archaeological discoveries from this era present us with an extremely rich type diversity, including wild and local primitive horses, the prototype of heavy draft horses, and fine riding horses.The objective of this study was to investigate the genetic structure of unexamined populations of three closely related horse breeds - the Danubian Nonius Hungarian Nonius and Serbian Nonius horses. A 608 bp long fragment of the mtDNA D-loop region was amplified and sequenced. The obtained results showed completely different genetic profiles between the investigated breeds. We identified nine of the 17 haplogroups described in modern horses. Most of the obtained sequences fell into M, L, G, and O'P lineages, which reflects the genetic profiles of the ancestral mares that were probably used at the initial stages of formation of the breeds. The population of the Danubian horse was characterized by a high prevalence of Central Asian specific haplogroup G (45%), followed by Western Eurasian specific haplogroups L and M (both about 21%). In contrast to the Danubian horse, in the Nonius breed the highest frequency of Western Eurasian haplogroup M (43.5%) was found, followed by Middle Eastern haplogroups O'P (26.1%) Central Asian specific E (13.0%) and G (13.1%). The Serbian Nonius horse showed a completely different genetic profile with a prevalence of the rare for Europe haplogroup D (66.7%), followed by Central Asian specific G (16.7%). The high mitochondrial haplotype diversity (Hd = 0.886) found in the investigated samples is evidence for multiple maternal origins in all populations.In conclusion, the obtained results demonstrated a high percentage of haplogroup sharing especially in the Danubian and Hungarian Nonius horse breeds, which reflects the possible common origins of the two breeds. In contrast to these breeds, the Serbian Nonius, despite the small number of investigated animals, showed a specific genetic profile, which could be explained by different and independent origins.

Genetic diversity and historical demography of underutilised goat breeds in North-Western Europe

In the last decade, several studies aimed at dissecting the genetic architecture of local small ruminant breeds to discover which variations are involved in the process of adaptation to environmental conditions, a topic that has acquired priority due to climate change. Considering that traditional breeds are a reservoir of such important genetic variation, improving the current knowledge about their genetic diversity and origin is the first step forward in designing sound conservation guidelines. The genetic composition of North-Western European archetypical goat breeds is still poorly exploited. In this study we aimed to fill this gap investigating goat breeds across Ireland and Scandinavia, including also some other potential continental sources of introgression. The PCA and Admixture analyses suggest a well-defined cluster that includes Norwegian and Swedish breeds, while the crossbred Danish landrace is far apart, and there appears to be a close relationship between the Irish and Saanen goats. In addition, both graph representation of historical relationships among populations and f4-ratio statistics suggest a certain degree of gene flow between the Norse and Atlantic landraces. Furthermore, we identify signs of ancient admixture events of Scandinavian origin in the Irish and in the Icelandic goats. The time when these migrations, and consequently the introgression, of Scandinavian-like alleles occurred, can be traced back to the Viking colonisation of these two isles during the Viking Age (793-1066 CE). The demographic analysis indicates a complicated history of these traditional breeds with signatures of bottleneck, inbreeding and crossbreeding with the improved breeds. Despite these recent demographic changes and the historical genetic background shaped by centuries of human-mediated gene flow, most of them maintained their genetic identity, becoming an irreplaceable genetic resource as well as a cultural heritage.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

MTaxi: A comparative tool for taxon identification of ultra low coverage ancient genomes

A major challenge in zooarchaeology is to morphologically distinguish closely related species' remains, especially using small bone fragments. Shotgun sequencing aDNA from archeological remains and comparative alignment to the candidate species' reference genomes will only apply when reference nuclear genomes of comparable quality are available, and may still fail when coverages are low. Here, we propose an alternative method, MTaxi, that uses highly accessible mitochondrial DNA (mtDNA) to distinguish between pairs of closely related species from ancient DNA sequences. MTaxi utilises mtDNA transversion-type substitutions between pairs of candidate species, assigns reads to either species, and performs a binomial test to determine the sample taxon. We tested MTaxi on sheep/goat and horse/donkey data, between which zooarchaeological classification can be challenging in ways that epitomise our case. The method performed efficiently on simulated ancient genomes down to 0.3x mitochondrial coverage for both sheep/goat and horse/donkey, with no false positives. Trials on n=18 ancient sheep/goat samples and n=10 horse/donkey samples of known species identity also yielded 100% accuracy. Overall, MTaxi provides a straightforward approach to classify closely related species that are difficult to distinguish through zooarchaeological methods using low coverage aDNA data, especially when similar quality reference genomes are unavailable. MTaxi is freely available at https://github.com/goztag/MTaxi.

Mitochondrial DNA sequencing of Kehilan and Hamdani horses from Saudi Arabia

The Arabian horse breed is well known for its purity and played a key role in the genetic improvement of other horses worldwide. The mitochondrial genome plays a vital role in maternal inheritance and it's helpful to evaluate its genetic diversity and conservation. It has higher mutation rates than nuclear DNA in vertebrates and therefore reveals phylogenetic relationships and haplotypes. In this study, the mitochondrial genome mutations in two Saudi horse strains, Kehilan and Hamdani demonstrated various changes in the gene and amino acid levels and included two other Saudi horses (Hadban and Seglawi) from the previous study for phylogenetic comparison. The whole mitochondrial genome sequencing resulted in intra and inter mtDNA variations between the studied horses. Interestingly, the Hamdani horse has nucleotide substitutions similar to those of the Hadban horse, which is reflected in the phylogenetic tree as a significantly close relationship. This type of study provides a better understanding of mitogenome structure and conservation of livestock species genetic data.

The Characteristics, Distribution, Function, and Origin of Alternative Lateral Horse Gaits

This article traces the characteristics, origin, distribution, and function of alternative lateral horse gaits, i.e., intermediate speed lateral-sequence gaits. Such alternative lateral gaits (running walk, rack, broken pace, hard pace, and broken trot) are prized by equestrians today for their comfort and have been found in select horse breeds for hundreds of years and even exhibited in fossil equid trackways. After exploring the evolution and development of alternative lateral gaits via fossil equid trackways, human art, and historical writings, the functional and genetic factors that led to the genesis of these gaits are discussed. Such gaited breeds were particularly favored and spread by the Scythians, Celts, Turks, and Spaniards. Fast and low-swinging hard pacing gaits are common in several horse breeds of mountainous areas of East and North Asia; high-stepping rack and running walk gaits are often displayed in European and North and South American breeds; the broken pace is found in breeds of Central Asia, Southeast Asia, West Asia, Western North America, and Brazil in South America; and the broken trot occurs in breeds of North Asia, South Asia, the Southern United States, and Brazil in South America, inhabiting desert or marshy areas.

Genetic History of the Altai Breed Horses: From Ancient Times to Modernity

This study focuses on expanding knowledge about the genetic diversity of the Altai horse native to Siberia. While studying modern horses from two Altai regions, where horses were subjected to less crossbreeding, we tested the hypothesis, formulated on the basis of morphological data, that the Altai horse is represented by two populations (Eastern and Southern) and that the Mongolian horse has a greater genetic proximity to Eastern Altai horses. Bone samples of ancient horses from different cultures of Altai were investigated to clarify the genetic history of this horse breed. As a genetic marker, we chose hypervariable region I of mitochondrial DNA. The results of the performed phylogenetic and population genetic analyses of our and previously published data confirmed the hypothesis stated above. As we found out, almost all the haplotypes of the ancient domesticated horses of Altai are widespread among modern Altai horses. The differences between the mitochondrial gene pools of the ancient horses of Altai and Mongolia are more significant than between those of modern horses of the respective regions, which is most likely due to an increase in migration processes between these regions after the Early Iron Age.

Characterization of the Caudal Ventral Tubercle in the Sixth Cervical Vertebra in Modern Equus ferus caballus

This study examined the anomalous variations of the ventral process of C6 in modern E. ferus caballus. The aim was to provide an incremental grading protocol measuring the absence of the caudal ventral tubercle (CVT) in this ventral process. The findings revealed the most prevalent absent CVT (aCVT) was left unilateral (n = 35), with bilateral (n = 29) and right unilateral (n = 12). Grading was determined in equal increments of absence 1/4, 2/4, 3/4, with 4/4 representing a complete aCVT in 56/76, with a significance of p = 0.0013. This also applied to bilateral specimens. In those C6 osseous specimens displaying a 4/4 grade aCVT, 41/56 had a partial absence of the caudal aspect of the cranial ventral tubercle (CrVT). Here, grading absent CrVTs (aCrVT) followed similarly to aCVTs, though 4/4 was not observed. The significance between 4/4 grade aCVTs and the presentation of an aCrVT was left p = 0.00001 and right p = 0.00018. In bilateral specimens, C6 morphologically resembled C5, implying a homeotic transformation that limited the attachment sites for the cranial and thoracal longus colli muscle. This potentially diminishes function and caudal cervical stability. Therefore, it is recommended that further studies examine the morphological extent of this equine complex vertebral malformation (ECVM) as well as its interrelationships and genetic code/blueprint.

Imputed genomes of historical horses provide insights into modern breeding

Historical genomes can provide important insights into recent genomic changes in horses, especially the development of modern breeds. In this study, we characterized 8.7 million genomic variants from a panel of 430 horses from 73 breeds, including newly sequenced genomes from 20 Clydesdales and 10 Shire horses. We used this modern genomic variation to impute the genomes of four historically important horses, consisting of publicly available genomes from 2 Przewalski's horses, 1 Thoroughbred, and a newly sequenced Clydesdale. Using these historical genomes, we identified modern horses with higher genetic similarity to those in the past and unveiled increased inbreeding in recent times. We genotyped variants associated with appearance and behavior to uncover previously unknown characteristics of these important historical horses. Overall, we provide insights into the history of Thoroughbred and Clydesdale breeds and highlight genomic changes in the endangered Przewalski's horse following a century of captive breeding.

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.

Improving the extraction of ancient Yersinia pestis genomes from the dental pulp

Ancient DNA preserved in the dental pulp offers the opportunity to characterize the genome of some of the deadliest pathogens in human history. However, while DNA capture technologies help, focus sequencing efforts, and therefore, reduce experimental costs, the recovery of ancient pathogen DNA remains challenging. Here, we tracked the kinetics of ancient Yersinia pestis DNA release in solution during a pre-digestion of the dental pulp. We found that most of the ancient Y. pestis DNA is released within 60 min at 37°C in our experimental conditions. We recommend a simple pre-digestion as an economical procedure to obtain extracts enriched in ancient pathogen DNA, as longer digestion times release other types of templates, including host DNA. Combining this procedure with DNA capture, we characterized the genome sequences of 12 ancient Y. pestis bacteria from France dating to the second pandemic outbreaks of the 17^th and 18^th centuries Common Era.

Whole-genome sequencing of cryopreserved resources from French Large White pigs at two distinct sampling times reveals strong signatures of convergent and divergent selection between the dam and sire lines

BACKGROUND

Numerous genomic scans for positive selection have been performed in livestock species within the last decade, but often a detailed characterization of the detected regions (gene or trait under selection, timing of selection events) is lacking. Cryopreserved resources stored in reproductive or DNA gene banks offer a great opportunity to improve this characterization by providing direct access to recent allele frequency dynamics, thereby differentiating between signatures from recent breeding objectives and those related to more ancient selection constraints. Improved characterization can also be achieved by using next-generation sequencing data, which helps narrowing the size of the detected regions while reducing the number of associated candidate genes.

METHODS

We estimated genetic diversity and detected signatures of recent selection in French Large White pigs by sequencing the genomes of 36 animals from three distinct cryopreserved samples: two recent samples from dam (LWD) and sire (LWS) lines, which had diverged from 1995 and were selected under partly different objectives, and an older sample from 1977 prior to the divergence.

RESULTS

French LWD and LWS lines have lost approximately 5% of the SNPs that segregated in the 1977 ancestral population. Thirty-eight genomic regions under recent selection were detected in these lines and the corresponding selection events were further classified as convergent between lines (18 regions), divergent between lines (10 regions), specific to the dam line (6 regions) or specific to the sire line (4 regions). Several biological functions were found to be significantly enriched among the genes included in these regions: body size, body weight and growth regardless of the category, early life survival and calcium metabolism more specifically in the signatures in the dam line and lipid and glycogen metabolism more specifically in the signatures in the sire line. Recent selection on IGF2 was confirmed and several other regions were linked to a single candidate gene (ARHGAP10, BMPR1B, GNA14, KATNA1, LPIN1, PKP1, PTH, SEMA3E or ZC3HAV1, among others).

CONCLUSIONS

These results illustrate that sequencing the genome of animals at several recent time points generates considerable insight into the traits, genes and variants under recent selection in a population. This approach could be applied to other livestock populations, e.g. by exploiting the rich biological resources stored in cryobanks.

DNA methylation-based profiling of horse archaeological remains for age-at-death and castration

Age profiling of archaeological bone assemblages can inform on past animal management practices, but is limited by the fragmentary nature of the fossil record and the lack of universal skeletal markers for age. DNA methylation clocks offer new, albeit challenging, alternatives for estimating the age-at-death of ancient individuals. Here, we take advantage of the availability of a DNA methylation clock based on 31,836 CpG sites and dental age markers in horses to assess age predictions in 84 ancient remains. We evaluate our approach using whole-genome sequencing data and develop a capture assay providing reliable estimates for only a fraction of the cost. We also leverage DNA methylation patterns to assess castration practice in the past. Our work opens for a deeper characterization of past husbandry and ritual practices and holds the potential to reveal age mortality profiles in ancient societies, once extended to human remains.

EWF: simulating exact paths of the Wright-Fisher diffusion

MOTIVATION

The Wright-Fisher diffusion is important in population genetics in modelling the evolution of allele frequencies over time subject to the influence of biological phenomena such as selection, mutation and genetic drift. Simulating the paths of the process is challenging due to the form of the transition density. We present EWF, a robust and efficient sampler which returns exact draws for the diffusion and diffusion bridge processes, accounting for general models of selection including those with frequency dependence.

RESULTS

Given a configuration of selection, mutation and endpoints, EWF returns draws at the requested sampling times from the law of the corresponding Wright-Fisher process. Output was validated by comparison to approximations of the transition density via the Kolmogorov-Smirnov test and QQ plots.

AVAILABILITY AND IMPLEMENTATION

All softwares are available at https://github.com/JaroSant/EWF.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Diversity of mitochondrial D-loop haplotypes from ancient Thracian horses in Bulgaria

The domestication of the horse began possibly more than 5000 years ago in the western part of the Eurasian steppe, and according to the leading hypothesis, horses first spread from the Steppe toward the region of the Thracian culture, starting in the second half of the 2nd millennium BCE and flourished from the fifth to first centuries BCE, mainly located in present-day Bulgaria. We analyzed 17 horse bone remains excavated from Thracian archaeological sites (fourth to first centuries BCE) in Bulgaria and successfully identified 17 sequences representing 14 different haplotypes of the mitochondrial D-loop. Compared with the mtDNA haplotypes of modern horses around the world, ancient Thracian horses in Bulgaria are thought to be more closely related to modern horses of Southern Europe and less related to those of Central Asia. In addition, the haplotypes we obtained represented 11 previously reported modern horse mtDNA haplogroups: A, B, D, E, G, H, I, L, N, P, and Q. All the haplogroups contain modern and regionally predominant haplotypes occurring in Europe, the Middle East, and Central Asia. Our results indicate that Thracian horses in Bulgaria have had relatively high genetic diversity and are closely related to modern horse breeds.

eSMC: a statistical model to infer admixture events from individual genomics data

Background

Inferring historical population admixture events yield essential insights in understanding a species demographic history. Methods are available to infer admixture events in demographic history with extant genetic data from multiple sources. Due to the deficiency in ancient population genetic data, there lacks a method for admixture inference from a single source. Pairwise Sequentially Markovian Coalescent (PSMC) estimates the historical effective population size from lineage genomes of a single individual, based on the distribution of the most recent common ancestor between the diploid’s alleles. However, PSMC does not infer the admixture event.

Results

Here, we proposed eSMC, an extended PSMC model for admixture inference from a single source. We evaluated our model’s performance on both in silico data and real data. We simulated population admixture events at an admixture time range from 5 kya to 100 kya (5 years/generation) with population admix ratio at 1:1, 2:1, 3:1, and 4:1, respectively. The root means the square error is $$\pm 7.61$$ ± 7.61 kya for all experiments. Then we implemented our method to infer the historical admixture events in human, donkey and goat populations. The estimated admixture time for both Han and Tibetan individuals range from 60 kya to 80 kya (25 years/generation), while the estimated admixture time for the domesticated donkeys and the goats ranged from 40 kya to 60 kya (8 years/generation) and 40 kya to 100 kya (6 years/generation), respectively. The estimated admixture times were concordance to the time that domestication occurred in human history.

Conclusion

Our eSMC effectively infers the time of the most recent admixture event in history from a single individual’s genomics data. The source code of eSMC is hosted at https://github.com/zachary-zzc/eSMC.

Reconstruction of the Major Maternal and Paternal Lineages in the Feral New Zealand Kaimanawa Horses

New Zealand has the fourth largest feral horse population in the world. The Kaimanawas (KHs) are feral horses descended from various domestic horse breeds released into the Kaimanawa ranges in the 19th and 20th centuries. Over time, the population size has fluctuated dramatically due to hunting, large-scale farming and forestry. Currently, the herd is managed by an annual round-up, limiting the number to 300 individuals to protect the native ecosystem. Here, we genotyped 96 KHs for uniparental markers (mitochondrial DNA, Y-chromosome) and assessed their genetic similarity with respect to other domestic horses. We show that at least six maternal and six paternal lineages contributed unequally to the KH gene pool, and today's KH population possibly represents two sub-populations. Our results indicate that three horse breeds, namely Welsh ponies, Thoroughbreds and Arabian horses had a major influence in the genetic-makeup of the extant KH population. We show that mitochondrial genetic diversity in KHs (π = 0.00687 ± 0.00355) is closer to that of the Sable Island horses (π = 0.0034 ± 0.00301), and less than other feral horse populations around the world. Our current findings, combined with ongoing genomic research, will provide insight into the population-specific genetic variation and inbreeding among KHs. This will largely advance equine research and improve the management of future breeding programs of these treasured New Zealand horse.

Unlocking Horse Y Chromosome Diversity

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

Convergence in domesticated fungi used for cheese and dry-cured meat maturation: beneficial traits, genomic mechanisms, and degeneration

Humans have domesticated genetically distant fungi for similar uses, the fermentation of lipid-rich and sugar-rich food to generate attractive aspects, odor and aroma, and to improve shelf life and product safety. Multiple independent domestication events also occurred within species. We review recent evidence of phenotypic convergence during the domestication of fungi for making cheese (Saccharomyces cerevisiae, Penicillium roqueforti, P. camemberti, and Geotrichum candidum) and for dry-cured meat making (P. nalgiovense and P. salamii). Convergence following adaptation to similar ecological niches involved colony aspect (fluffiness and color), lipolysis, proteolysis, volatile compound production, and competitive ability against food spoilers. We review evidence for convergence in genetic diversity loss in domesticated populations and in the degeneration of unused traits, such as toxin production and sexual reproduction. Phenotypic convergence sometimes occurred by similar mechanisms of genomic adaptation, in particular horizontal gene transfers and loss of genes.

Following the niche: the differential impact of the last glacial maximum on four European ungulates

Predicting the effects of future global changes on species requires a better understanding of the ecological niche dynamics in response to climate; the large climatic fluctuations of the last 50,000 years can be used as a natural experiment to that aim. Here we test whether the realized niche of horse, aurochs, red deer, and wild boar changed between 47,000 and 7500 years ago using paleoecological modelling over an extensive archaeological database. We show that they all changed their niche, with species-specific responses to climate fluctuations. We also suggest that they survived the climatic turnovers thanks to their flexibility and by expanding their niche in response to the extinction of competitors and predators. Irrespective of the mechanism behind such processes, the fact that species with long generation times can change their niche over thousands of years cautions against assuming it to stay constant both when reconstructing the past and predicting the future.

European megafaunal ungulates living in open habitats over the last 50,000 years showed evidence for niche change, possibly driven by climatic change and extinction of competitors and predators

The genomic history and global expansion of domestic donkeys

Donkeys transformed human history as essential beasts of burden for long-distance movement, especially across semi-arid and upland environments. They remain insufficiently studied despite globally expanding and providing key support to low- to middle-income communities. To elucidate their domestication history, we constructed a comprehensive genome panel of 207 modern and 31 ancient donkeys, as well as 15 wild equids. We found a strong phylogeographic structure in modern donkeys that supports a single domestication in Africa ~5000 BCE, followed by further expansions in this continent and Eurasia and ultimately returning to Africa. We uncover a previously unknown genetic lineage in the Levant ~200 BCE, which contributed increasing ancestry toward Asia. Donkey management involved inbreeding and the production of giant bloodlines at a time when mules were essential to the Roman economy and military.

Demographic inference provides insights into the extirpation and ecological dominance of eusocial snapping shrimps

Although eusocial animals often achieve ecological dominance in the ecosystems where they occur, many populations are unstable, resulting in local extinction. Both patterns may be linked to the characteristic demography of eusocial species-high reproductive skew and reproductive division of labor support stable effective population sizes that make eusocial groups more competitive in some species, but also lower effective population sizes that increase susceptibility to population collapse in others. Here, we examine the relationship between demography and social organization in Synalpheus snapping shrimps, a group in which eusociality has evolved recently and repeatedly. We show using coalescent demographic modelling that eusocial species have had lower but more stable effective population sizes across 100,000 generations. Our results are consistent with the idea that stable population sizes may enable competitive dominance in eusocial shrimps, but they also suggest that recent population declines are likely caused by eusocial shrimps' heightened sensitivity to environmental changes, perhaps as a result of their low effective population sizes and localized dispersal. Thus, although the unique life histories and demography of eusocial shrimps have likely contributed to their persistence and ecological dominance over evolutionary timescales, these social traits may also make them vulnerable to contemporary environmental change.

Population structure and genomic footprints of selection in five major Iranian horse breeds

The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log₁₀ rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log₁₀ rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.

mapDATAge: a ShinyR package to chart ancient DNA data through space and time

SUMMARY

Ancient DNA datasets are increasingly difficult to visualize for users lacking computational experience. Here, we describe mapDATAge, which aims to provide user-friendly automated modules for the interactive mapping of allele, haplogroup and/or ancestry distributions through space and time. mapDATAge enhances collaborative data sharing while assisting the assessment and reporting of spatiotemporal patterns of genetic changes.

AVAILABILITY AND IMPLEMENTATION

mapDATAge is a Shiny R application designed for exploring spatiotemporal patterns in ancient DNA data through a graphical user interface. It is freely available under GNU Public License in Github: https://github.com/xuefenfei712/mapDATAge.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bioarchaeological evidence of one of the earliest Islamic burials in the Levant

The Middle East plays a central role in human history harbouring a vast diversity of ethnic, cultural and religious groups. However, much remains to be understood about past and present genomic diversity in this region. Here we present a multidisciplinary bioarchaeological analysis of two individuals dated to the late 7th and early 8th centuries, the Umayyad Era, from Tell Qarassa, an open-air site in modern-day Syria. Radiocarbon dates and burial type are consistent with one of the earliest Islamic Arab burials in the Levant. Interestingly, we found genomic similarity to a genotyped group of modern-day Bedouins and Saudi rather than to most neighbouring Levantine groups. This study represents the genomic analysis of a secondary use site with characteristics consistent with an early Islamic burial in the Levant. We discuss our findings and possible historic scenarios in the light of forces such as genetic drift and their possible interaction with religious and cultural processes (including diet and subsistence practices).

Ancient genomic and archaeological data combine to identify a surprisingly early Islamic burial in modern day Syria.

Radiocarbon and genomic evidence for the survival of Equus Sussemionus until the late Holocene

The exceptionally rich fossil record available for the equid family has provided textbook examples of macroevolutionary changes. Horses, asses, and zebras represent three extant subgenera of Equus lineage, while the Sussemionus subgenus is another remarkable Equus lineage ranging from North America to Ethiopia in the Pleistocene. We sequenced 26 archaeological specimens from Northern China in the Holocene that could be assigned morphologically and genetically to Equus ovodovi, a species representative of Sussemionus. We present the first high-quality complete genome of the Sussemionus lineage, which was sequenced to 13.4× depth of coverage. Radiocarbon dating demonstrates that this lineage survived until ~3500 years ago, despite continued demographic collapse during the Last Glacial Maximum and the great human expansion in East Asia. We also confirmed the Equus phylogenetic tree and found that Sussemionus diverged from the ancestor of non-caballine equids ~2.3–2.7 million years ago and possibly remained affected by secondary gene flow post-divergence. We found that the small genetic diversity, rather than enhanced inbreeding, limited the species’ chances of survival. Our work adds to the growing literature illustrating how ancient DNA can inform on extinction dynamics and the long-term resilience of species surviving in cryptic population pockets.

Ice Age megafauna rock art in the Colombian Amazon?

Megafauna paintings have accompanied the earliest archaeological contexts across the continents, revealing a fundamental inter-relationship between early humans and megafauna during the global human expansion as unfamiliar landscapes were humanized and identities built into new territories. However, the identification of extinct megafauna from rock art is controversial. Here, we examine potential megafauna depictions in the rock art of Serranía de la Lindosa, Colombian Amazon, that includes a giant sloth, a gomphothere, a camelid, horses and three-toed ungulates with trunks. We argue that they are Ice Age rock art based on the (i) naturalistic appearance and diagnostic morphological features of the animal images, (ii) late Pleistocene archaeological dates from La Lindosa confirming the contemporaneity of humans and megafauna, (iii) recovery of ochre pigments in late Pleistocene archaeological strata, (iv) the presence of most megafauna identified in the region during the late Pleistocene as attested by archaeological and palaeontological records, and (v) widespread depiction of extinct megafauna in rock art across the Americas. Our findings contribute to the emerging picture of considerable geographical and stylistic variation of geometric and figurative rock art from early human occupations across South America. Lastly, we discuss the implications of our findings for understanding the early human history of tropical South America.

This article is part of the theme issue ‘Tropical forests in the deep human past’.

Metagenomic investigation of the equine faecal microbiome reveals extensive taxonomic diversity

Background

The horse plays crucial roles across the globe, including in horseracing, as a working and companion animal and as a food animal. The horse hindgut microbiome makes a key contribution in turning a high fibre diet into body mass and horsepower. However, despite its importance, the horse hindgut microbiome remains largely undefined. Here, we applied culture-independent shotgun metagenomics to thoroughbred equine faecal samples to deliver novel insights into this complex microbial community.

Results

We performed metagenomic sequencing on five equine faecal samples to construct 123 high- or medium-quality metagenome-assembled genomes from Bacteria and Archaea. In addition, we recovered nearly 200 bacteriophage genomes. We document surprising taxonomic diversity, encompassing dozens of novel or unnamed bacterial genera and species, to which we have assigned new Candidatus names. Many of these genera are conserved across a range of mammalian gut microbiomes.

Conclusions

Our metagenomic analyses provide new insights into the bacterial, archaeal and bacteriophage components of the horse gut microbiome. The resulting datasets provide a key resource for future high-resolution taxonomic and functional studies on the equine gut microbiome.

The Diversity Decline in Wild and Managed Honey Bee Populations Urges for an Integrated Conservation Approach

Many parts of the globe experience severe losses and fragmentation of habitats, affecting the self-sustainability of pollinator populations. A number of bee species coexist as wild and managed populations. Using honey bees as an example, we argue that several management practices in beekeeping threaten genetic diversity in both wild and managed populations, and drive population decline. Large-scale movement of hive stocks, introductions into new areas, breeding programs and trading of queens contribute to reducing genetic diversity, as recent research demonstrated for wild and managed honey bees within a few decades. Examples of the effects of domestication in other organisms show losses of both genetic diversity and fitness functions. Cases of natural selection and feralization resulted in maintenance of a higher genetic diversity, including in a Varroa destructor surviving population of honey bees. To protect the genetic diversity of honey bee populations, exchange between regions should be avoided. The proposed solution to selectively breed all local subspecies for a use in beekeeping would reduce the genetic diversity of each, and not address the value of the genetic diversity present in hybridized populations. The protection of Apis mellifera’s, Apis cerana’s and Apis koschevnikovi’s genetic diversities could be based on natural selection. In beekeeping, it implies to not selectively breed but to leave the choice of the next generation of queens to the colonies, as in nature. Wild populations surrounded by beekeeping activity could be preserved by allowing Darwinian beekeeping in a buffer zone between the wild and regular beekeeping area.

Morphological Characteristics of a Horse Discovered in an Avar-Period Grave from Sâncraiu de Mureș (Alba County, Romania)

Simple Summary

The present study deals with the morphological description of the animal remains in an Avar grave that contained the remnants of a horse burial, along with its owner, in a cemetery situated in Sâncraiu de Mures, Alba County, Romania, located geographically in the Transylvanian Plateau. The investigation of the grave and the data about the remains were published by archaeologists from the Institute of Archaeology of Cluj Napoca, Romania. For analysis and study, the animal bone remains were sent to the Faculty of Veterinary Medicine Cluj-Napoca. After a careful analysis of the morphological characteristics of the bones, the conclusions showed that the bones belonged to one young horse of 2–2.5 years, with an average recalculated shoulder height as being situated in the 1300–1400 mm interval. The study provides a set of morphological data regarding the characteristics of the horses from this scantily represented history in Romania, and the results show that the data obtained are consistent with the similar Avar horse burials in the region.

Abstract

Archeozoological studies provide an insight into human–environment relations, bringing important information on the morphology of the animal and the role of the animal and its functions. The purpose of this study was to reveal the morphological characteristics of the horse identified in an 8th century BC Avar cemetery dated, as it resulted from the investigation carried out on the materials presented to the Anatomy Lab of Department of Comparative Anatomy. The cleaning and conditioning of the materials were performed in the lab, followed by anatomical and osteometrical study. The identification of the species, the osteometrical interpretation and the assessment of age at death constituted the basis of the main conclusions. Based on the morphological and metrical data, we concluded that the fragments originated from a single young horse individual (Equus caballus) no older than 2.5 (2–2.5 years age span) included into class 5 of height (Vitt scale) with slender extremities. An overview of the available data (horse morphology) from similar sites in Romania and neighboring territories (Hungary and Croatia) is presented, with the intention of a general framing for the characteristics of the horses used by the Avar populations.

One Hundred Years of Coat Colour Influences on Genetic Diversity in the Process of Development of a Composite Horse Breed

Genetic diversity and demographic parameters were computed to evaluate the historic effects of coat colour segregation in the process of configuration of the Hispano-Arabian horse (Há). Pedigree records from 207,100 individuals born between 1884 and 2019 were used. Although coat colour is not a determinant for the admission of Hispano-Arabian individuals as apt for breeding, it may provide a representative visual insight into the gene contribution of Spanish Purebred horses (PRE), given many of the dilution genes described in Há are not present in the Arabian Purebred breed (PRá). The lack of consideration of coat colour inheritance patterns by the entities in charge of individual registration and the dodging behaviour of breeders towards the historic banning policies, may have acted as a buffer for diversity loss (lower than 8%). Inbreeding levels ranged from 1.81% in smokey cream horses to 8.80 for white horses. Contextually, crossbred breeding may increase the likelihood for double dilute combinations to occur as denoted by the increased number of Há horses displaying Pearl coats (53 Há against 3 PRE and 0 PRá). Bans against certain coat colours and patterns may have prevented an appropriate registration of genealogical information from the 4th generation onwards for decades. This may have brought about the elongation of generation intervals. Breeder tastes may have returned to the formerly officially-recognised coat colours (Grey and Bay) and Chestnut/Sorrel. However, coat colour conditioning effects must be evaluated timely for relatively short specific periods, as these may describe cyclic patterns already described in owners’ and breeders’ tastes over the centuries.

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.

A single-nucleotide mutation within the TBX3 enhancer increased body size in Chinese horses

Chinese ponies are endemic to the mountainous areas of southwestern China and were first reported in the archaeological record at the Royal Tomb of Zhongshan King, Mancheng, dated to approximately ∼2,100 YBP.¹ Previous work has started uncovering the genetic basis of size variation in western ponies and horses, revealing a limited number of loci, including HMGA2,²LCORL/NCAPG,³ZFAT, and LASP1.^4,5 Whether the same genetic pathways also drive the small body size of Chinese ponies, which show striking anatomical differences to Shetland ponies,⁶ remains unclear.^2,7 To test this, we combined whole-genome sequences of 187 horses across China. Statistical analyses revealed top association between genetic variation at the T-box transcription factor 3 (TBX3) and the body size. Fine-scale analysis across an extended population of 189 ponies and 574 horses narrowed down the association to one A/G SNP at an enhancer region upstream of the TBX3 (ECA8:20,644,555, p = 2.34e−39). Luciferase assays confirmed the single-nucleotide G mutation upregulating TBX3 expression, and enhancer-knockout mice exhibited shorter limbs than wild-type littermates (p < 0.01). Re-analysis of ancient DNA data showed that the G allele, which is most frequent in modern horses, first occurred some ∼2,300 years ago and rose in frequency since. This supports selection for larger size in Asia from approximately the beginning of the Chinese Empire. Overall, this study characterized the causal regulatory mutation underlying small body size in Chinese ponies and revealed size as one of the main selection targets of past Chinese breeders.

•

One single A/G SNP in TBX3 enhancer region drives size variation in Chinese horses

•

The frequency of the G variant correlates positively with size in 763 horses

•

Cellular and mice models confirm it affects TBX3 transcription and the limb length

•

The G variant first occurred ∼2,300 years ago and rose in frequency since

The Sequence Analysis of Mitochondrial DNA Revealed Some Major Centers of Horse Domestications: The Archaeologist's Cut

The question about the time and the place of horse domestication, a process which had a profound impact on the progress of mankind, is disputable. According to the most widely accepted hypothesis, the earliest domestication of the horse happened in the western parts of the Eurasian steppes, between the Northern Black Sea region and present-day Kazakhstan and Turkmenistan. It seems that it occurred not earlier than the first half and most probably during the middle (even the last third) of the fourth millennium BC (from ∼ 5.5 kya). The next steps of large-scale horse breeding occurred almost simultaneously in Eurasia and North Africa due to the development of the social structure of human communities. On the other hand, the morphological differences between wild and domestic animals are rather vague and the genetic introgression between them is speculative. In this review, we have tried to gather all available scientific data on the existing possible hypotheses for the earliest domestication of the horse, as well as to highlight some data on the most plausible ones. This is due to the frequency of some significant data on the frequency of strictly defined mitotypes in different historical periods of human civilizations existing in the same periods.