Population Genetic Analysis of the Estonian Native Horse Suggests Diverse and Distinct Genetics, Ancient Origin and Contribution from Unique Patrilines

Influence of Sires on Population Substructure in Dülmen Wild Horses

The objectives of the present study were to analyze the influence of the stallions employed in the Dülmen wild horses on the genetic diversity and population substructure using Bayesian cluster analysis. The Dülmen wild horse is maintained as a unique horse population exposed to the natural conditions all year round in the Merfelder Bruch near Dülmen in Westphalia, Germany. Stallions selected for breeding have to prove their abilities to survive under this harsh environment. We used multilocus genotypic information from a set of 29 autosomal microsatellites to determine the paternity of 185 male foals sired by nine stallions. As females could not be sampled, we could not make inferences on all yearlings and test whether there are differences in the genetic population parameters between both sexes. The mean number of progeny was 19.92 with a range of 2-32, caused by the length of the service period per stallion. The average observed and unbiased expected heterozygosity was 0.688 and 0.631, the mean number of alleles was 4.448, and Wright's FIS was -0.173. Pairwise genetic distances (FST and Nei's unbiased genetic distances) were significant and varied between 0.038 to 0.091 and 0.085 to 0.290, respectively. Neighbor-joining dendrogram plots clustered a large proportion of the paternal progeny groups in different branches. Posterior Bayesian analyses using seven paternal half-sib groups with 10-74 members supported a maximum of six clusters, with two paternal progeny groups not differing, and a median of five clusters, with two groups of two sires each falling into the same clusters. When sires were employed in non-consecutive years, progeny from these different years of the same sires were grouped in the same cluster, whereas the progeny of one sire from two consecutive years were in different clusters. We were able to distinguish male progeny from Dülmen wild horse stallions and to show the effects of stallion use on the genetic substructure in the Dülmen wild horse herd. In conclusion, the analyses showed the genetic potential of the Dülmen wild horse stallions to maintain a high genetic diversity and also the effects in which breeding seasons and for how long stallions are used to sire foals. The selection of stallions may be sensitive for the further development of genetic diversity and preserve this closed population as a valuable resource for further studies on the evolution of the horse.

Genetic Diversity and Population Structure of Dülmen Wild, Liebenthal and Polish Konik Horses in Comparison with Przewalski, Sorraia, German Draught and Riding Horses

The objective of the present study was to analyze the genetic diversity, individual-based assessment of population structure, and admixture in the Dülmen wild horse population in comparison to warmblood, coldblood, and primitive horse populations. The Dülmen wild horse is kept as a unique horse population in the Merfelder Bruch near Dülmen in Westphalia, Germany, and since 1856 has been managed by the Dukes of Croÿ. The Dülmen wild horse population is exposed to the natural conditions of the Merfelder Bruch all year round without human interventions for feeding and veterinary care. In the present study, genetic diversity was estimated for 101 Dülmen wild horses using multilocus genotypic information from a set of 29 autosomal microsatellites and compared with 587 horses from 17 different horse populations. Dülmen wild horses maintained a high degree of genetic diversity, with an average observed heterozygosity of 0.68, a mean number of 6.17 alleles, and heterozygote deficit of -0.035. Pairwise genetic distances (FST, Nei's standard, and Cavalli-Sforza distances) were closest to German coldblood breeds, Polish Konik, and Icelandic horses and most divergent from Sorraia and Przewalski's horses. Neighbor joining dendrogram and PCA plots showed a clear distinction of Dülmen wild horses from other populations, particularly from Przewalski horses. Posterior Bayesian analysis confirmed clear differentiation from other horse populations without an admixture pattern and a high membership index (0.92). It was possible to distinguish Dülmen wild horses from Dülmen and Polish Konik horses. In conclusion, Dülmen wild horses show a notable separation from other German horse breeds and primitive horse populations and may serve as a resource to study evolution of equine domestication.

Mitochondrial DNA and Y chromosome reveal the genetic structure of the native Polish Konik horse population

Polish Konik remains one of the most important horse breeds in Poland. The primitive, native horses with a stocky body and mouse-like coat color are protected by a conservation program, while their Polish population consists of about 3,480 individuals, representing 16 dam and six sire lines. To define the population's genetic structure, mitochondrial DNA and Y chromosome sequence variables were identified. The mtDNA whole hypervariable region analysis was carried out using the Sanger sequencing method on 233 Polish Koniks belonging to all dam lines, while the Y chromosome analysis was performed with the competitive allele-specific PCR genotyping method on 36 horses belonging to all sire lines. The analysis of the mtDNA hypervariable region detected 47 SNPs, which assigned all tested horses to 43 haplotypes. Most dam lines presented more than one haplotype; however, five dam lines were represented by only one haplotype. The haplotypes were classified into six (A, B, E, J, G, R) recognized mtDNA haplogroups, with most horses belonging to haplogroup A, common among Asian horse populations. Y chromosome analysis allocated Polish Koniks in the Crown group, condensing all modern horse breeds, and divided them into three haplotypes clustering with coldblood breeds (28 horses), warmblood breeds (two horses), and Duelmener Pony (six horses). The clustering of all Wicek sire line stallions with Duelmener horses may suggest a historical relationship between the breeds. Additionally, both mtDNA and Y chromosome sequence variability results indicate crossbreeding before the studbooks closure or irregularities in the pedigrees occurred before the DNA testing introduction.

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.

Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.

Memories, museum artefacts and excavations in resolving the history of maternal lineages in the Finnhorse

We used historical DNA samples to examine the history of a native horse breed, the Finnhorse. Samples were collected from private collections, museums, schools and excavations, representing the times prior to, during, and after the foundation of the breed; from the end of the 19th century and throughout the 20th century. We sequenced a fragment of mitochondrial DNA from these historical samples to study the history and evolution of maternal lineages of horses back to the early days of the breed, compared the mitochondrial DNA sequence diversity of different historical periods and modern day Finnhorses, estimated the effective population sizes, and searched for both temporal and geographic population genetic structure. We observed high maternal haplotype and nucleotide diversity at the time during the foundation of the breed, and a decrease in both measures during 1931-1970. In addition, we observed losses of some haplotypes present in the early stages of the breed. There was only slight evidence of geographical or temporal population structure. This study is, to our knowledge, the first to use such temporal sampling to reveal the history of a specific animal breed.

Y Chromosome Haplotypes Enlighten Origin, Influence, and Breeding History of North African Barb Horses

In horses, demographic patterns are complex due to historical migrations and eventful breeding histories. Particularly puzzling is the ancestry of the North African horse, a founding horse breed, shaped by numerous influences throughout history. A genetic marker particularly suitable to investigate the paternal demographic history of populations is the non-recombining male-specific region of the Y chromosome (MSY). Using a recently established horse MSY haplotype (HT) topology and KASP™ genotyping, we illustrate MSY HT spectra of 119 Barb and Arab-Barb males, collected from the Maghreb region and European subpopulations. All detected HTs belonged to the Crown haplogroup, and the broad MSY spectrum reflects the wide variety of influential stallions throughout the breed’s history. Distinct HTs and regional disparities were characterized and a remarkable number of early introduced lineages were observed. The data indicate recent refinement with Thoroughbred and Arabian patrilines, while 57% of the dataset supports historical migrations between North Africa and the Iberian Peninsula. In the Barb horse, we detected the HT linked to Godolphin Arabian, one of the Thoroughbred founders. Hence, we shed new light on the question of the ancestry of one Thoroughbred patriline. We show the strength of the horse Y chromosome as a genealogical tool, enlighten recent paternal history of North African horses, and set the foundation for future studies on the breed and the formation of conservation breeding programs.

Occurrence of Equine Foamy Virus Infection in Horses from Poland

Equine foamy virus (EFVeca) is a foamy virus of non-primate origin and among the least-studied members of this retroviral subfamily. By sequence comparison, EFVeca shows the highest similarity to bovine foamy virus. In contrast to simian, bovine or feline foamy viruses, knowledge about the epidemiology of EFVeca is still limited. Since preliminary studies suggested EFVeca infections among horses in Poland, we aimed to expand the diagnostics of EFVeca infections by developing specific diagnostic tools and apply them to investigate its prevalence. An ELISA test based on recombinant EFVeca Gag protein was developed for serological investigation, while semi-nested PCR for the detection of EFVeca DNA was established. 248 DNA and serum samples from purebred horses, livestock and saddle horses, Hucul horses and semi-feral Polish primitive horses were analyzed in this study. ELISA was standardized, and cut off value, sensitivity and specificity of the test were calculated using Receiver Operating Characteristic and Bayesian estimation. Based on the calculated cut off, 135 horses were seropositive to EFVeca Gag protein, while EFVeca proviral DNA was detected in 85 animals. The rate of infected individuals varied among the horse groups studied; this is the first report confirming the existence of EFVeca infections in horses from Poland using virus-specific tools.

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.

Genetic Variability and Population Structure of Polish Konik Horse Maternal Lines Based on Microsatellite Markers

Simple Summary

The aim of this study is to reveal the genetic variability and population structure of maternal lines of the Polish Konik—a native Polish horse breed of the primitive type included in the Polish conservation programme. The analysis was carried out on the basis of 17 microsatellite markers routinely used for parentage testing. The structural analysis indicated the likelihood of three genetic clusters (using the Structure and Structure Harvester software). F-statistics indicated a low level of inbreeding. All mean population coefficients were close to those obtained for representatives of male founder lines. The population seemed to be stable. However, due to a previous bottleneck event, breeding strategies should focus on avoiding inbreeding depression, preventing the decrease of genetic variability, and sustaining the active female lines.

Abstract

The aim of the conservation programme is to maintain the population size of endangered livestock breeds of less economic importance at a level that ensures the survival of the breed, the preservation of genetic diversity, and the preservation of as many pedigree lines as possible. The Polish Konik, a native Polish primitive-type horse breed and is one of the breeds included in such a programme in Poland. Presently, there are only 16 (of the 35 maternal lines known in 1962), some of which are endangered. We examined the genetic variability and structure of the Polish Konik maternal lines (176 individuals) on the basis of the pedigree data and 17 microsatellite markers (STRs) from parentage testing. The overall mean number of alleles was 7.647 (±0.411), the effective number of alleles was 3.935 (±0.271), the mean number of alleles for which the frequency was equal to or lower than 5% was 4.471 (±0.286), and the mean information index was 1.506 (±0.087). The structure of the population and admixture patterns were calculated with the Structure and Structure Harvester software. The structural analysis indicated three likely genetic clusters; as the most optimal K value was estimated as 3, with ∆K of 15.4188. The F-statistics results indicated a low level of inbreeding (average inbreeding coefficient F_IT was 0.0188, coefficient of differentiation F_ST was 0.0304, and mean inbreeding index value F_IS was −0.0119). Variability monitoring should be carried out in order to avoid inbreeding depression, while breeding strategies should be designed to prevent the decrease of genetic variability in the Polish horse breed and to sustain the active female lines.

Modern Northern Domestic Horses Carry Mitochondrial DNA Similar to Przewalski’s Horse

Several recent studies have suggested past gene flow between the Przewalski’s horse and modern domestic horse and questioned the wild origin of the Przewalski’s horse. Mitochondrial DNA has placed representatives of the Przewalski’s horse into three among the eighteen haplogroups detected from the modern horse. Of these, two haplogroups have so far been found exclusively in the Przewalski’s horse, while the one shared with the domestic horse includes captive individuals that have uncertain pedigrees. We recently found five domestic horse individuals of North European horse breeds to carry a mitochondrial haplogroup that was previously confined only to the Przewalski’s horse. These individuals were sequenced for 6039 bp of mitochondrial DNA and used, together with domestic and Przewalski’s horse sequences presenting all horse haplogroups, to examine the phylogenetic relationships and to date the divergence time between Przewalski’s and domestic horse clusters within this haplogroup. The divergence was dated to have likely occurred about 13,300–11,400 years ago, which coincides with the time of the Younger Dryas.

Genetic Diversity and Population Structure of Polish Konik Horse Based on Individuals from All the Male Founder Lines and Microsatellite Markers

The Polish Konik horse is a primitive native breed included in the genetic resource conservation program in Poland. After World War II, intensive breeding work began, aimed at rebuilding this breed. Now, the whole Polish Konik population is represented by six male founder lines (Wicek, Myszak, Glejt I, Goraj, Chochlik and Liliput). Individuals representing all six paternal lineages were selected based on their breeding documentation. We performed a fragment analysis with 17 microsatellite markers (STRs) recommended by the International Society for Animal Genetics (ISAG). The genetic diversity and structure within the paternal lineages and the whole of the studied group were investigated. The average allelic richness was 6.497 for the whole studied group. The fixation index (FST; measure of population differentiation) was low (about 3%), the mean inbreeding coefficient (FIT) was low and close to 0, and the mean inbreeding index value (FIS) was negative. The mean expected heterozygosity was established at 0.7046 and was lower than the observed heterozygosity. The power of discrimination and power of exclusion were 99.9999%. The cumulative parentage exclusion probability equaled 99.9269% when one parental genotype was known and 99.9996% with both parents' genotypic information was available. About 3% of the genetic variation was caused by differences in the breed origin and about 97% was attributed to differences among individuals. Our analysis revealed that there has been no inbreeding in the Polish Konik breed for the studied population. The genetic diversity was high, and its parameters were similar to those calculated for native breeds from other countries reported in the literature. However, due to the small number of breed founders and paternal lineages with unknown representation, the population's genetic diversity and structure should be monitored regularly.

Conservation Strategies for Local Breed Biodiversity

Breeds of domesticated animals are often overlooked as contributing to biodiversity. Their unique role at the junction of natural and human-influenced environments makes them important potential contributors to biodiversity and ecosystem function. Effective conservation of livestock breeds rests upon valid recognition of breeds as repeatable genomic packages with a reasonably high level of predictability for performance. Local or landrace breeds that lack breeder organization are especially difficult to conserve due to lack of formal recognition as breeds. Achieving success with them involves three major steps: Discover, Secure, and Sustain. Early in the process an evaluation of candidate populations for status as genetic resources is essential. This process is aided by a phenotypic matrix which can be used alongside historical investigations and genetic (DNA) studies. The goal is to include all qualifying animals and to exclude all those that do not qualify. Securing some populations depends on careful rescue protocols for maximizing the recovery of genetic variation, and this can then be followed by breeding protocols that provide for maintaining the population’s production potential along with a healthy and viable genetic structure for long-term survival and use. Sustaining breeds for the long term is also enabled by assuring market demand for the breed and its products.