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

Rewilded horses in European nature conservation - a genetics, ethics, and welfare perspective

In recent decades, the integration of horses (Equus ferus) in European rewilding initiatives has gained widespread popularity due to their potential for regulating vegetation and restoring natural ecosystems. However, employing horses in conservation efforts presents important challenges, which we here explore and discuss. These challenges encompass the lack of consensus on key terms inherent to conservation and rewilding, the entrenched culture and strong emotions associated with horses, low genetic diversity and high susceptibility to hereditary diseases in animals under human selection, as well as insufficient consideration for the social behaviour of horses in wild-living populations. In addition, management of wild-living horses involves intricate welfare, ethics and legislative dimensions. Anthropocentric population-control initiatives may be detrimental to horse group structures since they tend to prioritise individual welfare over the health of populations and ecosystems. To overcome these challenges, we provide comprehensive recommendations. These involve a systematic acquisition of genetic information, a focus on genetic diversity rather than breed purity and minimal veterinary intervention in wild-living populations. Further, we advise allowing for natural top-down and bottom-up control - or, if impossible, simulating this by culling or non-lethal removal of horses - instead of using fertility control for population management. We advocate for intensified collaboration between conservation biologists and practitioners and enhanced communication with the general public. Decision-making should be informed by a thorough understanding of the genetic makeup, common health issues and dynamics, and social behaviour in wild-living horse populations. Such a holistic approach is essential to reconcile human emotions associated with horses with the implementation of conservation practices that are not only effective but also sustainable for the long-term viability of functional, biodiverse ecosystems, while rehabilitating the horse as a widespread wild-living species in Europe.

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.

Unravelling of the genetic diversity, population structure, and phylogenetic relationships of the indigenous horse breeds of Pakistan

The current study aimed to evaluate the genetic diversity, level of admixture, and phylogenetic relationship, of the Pakistani horse breeds, along with their morphological characterization. Data for the body measurements showed that Morna horses had the highest values of body height, body length, chest girth, leg length, and head length, whereas the Baluchi horses had the lowest values for these traits. For the genetic diversity 64 animals, 15/breed except for Baluchi(14) and Topras(05), were genotyped by using the 17-plex equine genotyping kit. The AMOVA results showed that 13% of genetic diversity was explained by breed differences, whereas 27% and 60% came from among and within individuals respectively. The highest values of genetic diversity parameters including Na(7.29±0.29), Ne(5.73±0.28), Ho(0.74±0.05) and He(0.82±0.01) were observed for Morna, whereas their lowest values were found for Topras. However, the highest value of inbreeding coefficient (Fis) was found for Baluchi and the lowest for KB horses. Among the markers, CA425 was found as the most polymorphic and ASB23 as the least polymorphic and highly fixed marker. Results of structure analysis revealed that, except Topras, each local horse breed had very admixed genetic structure perhaps due to their continuous crossing with other breeds in the past. Moreover, the structure analysis also showed that Morna and Shiean breeds had mixing of each other which was also confirmed by the lowest value of their pairwise Fst values and likewise the phylogenetic analysis also showed their close genetic relationship with each other. The phylogenetic analysis also revealed that Shiean breed had close genetic relationship with Arabian horses. Collectively, our data showed that Morna is the largest and genetically most diverse horse breed of Pakistan, whereas Baluchi horses are the smallest in size and have the highest values of inbreeding coefficient. And the phylogeny analysis showed that Shiean breed had close genetic relationship with the Arabian.

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.

Assessment and Distribution of Runs of Homozygosity in Horse Breeds Representing Different Utility Types

The present study reports runs of homozygosity (ROH) distribution in the genomes of six horse breeds (571 horses in total) representing three horse types (primitive, light, and draft horses) based on the 65k Equine BeadChip assay. Of major interest was the length, quantity, and frequency of ROH characteristics, as well as differences between horse breeds and types. Noticeable differences in the number, length and distribution of ROH between breeds were observed, as well as in genomic inbreeding coefficients. We also identified regions of the genome characterized by high ROH coverage, known as ROH islands, which may be signals of recent selection events. Eight to fourteen ROH islands were identified per breed, which spanned multiple genes. Many were involved in important horse breed characteristics, including WFIKNN2, CACNA1G, STXBP4, NOG, FAM184B, QDPR, LCORL, and the zinc finger protein family. Regions of the genome with zero ROH occurrences were also of major interest in specific populations. Depending on the breed, we detected between 2 to 57 no-ROH regions and identified 27 genes in these regions that were common for five breeds. These genes were involved in, e.g., muscle contractility (CACNA1A) and muscle development (miR-23, miR-24, miR-27). To sum up, the obtained results can be furthered analyzed in the topic of identification of markers unique for specific horse breed characteristics.

Genetic Diversity and Structure of the Main Danubian Horse Paternal Genealogical Lineages Based on Microsatellite Genotyping

Simple Summary

The Danubian horse was created on the former Klementina stud farm near Pleven to satisfy the demands of the Bulgarian Army for light draft horses and to improve the working capacity of the local horse population. The privatization of the Klementina stud farm in the late 1990s and the lack of economic activity have led to a sharp reduction in the number of mares and stallions and their sale to private owners. At present, only six of the main paternal lines which participated in the creation of the Danubian horse breed are preserved: Zdravko, NONIUS XVII-30, Torpedo, Lider, Kalifa, and Hrabar. This is the first study on the genetic and population diversity of the Danubian horse paternal lines based on microsatellite markers (STRs). The results showed that the studied Danubian horse population was characterized by a high level of genetic diversity with a medium value of 0.84. The fixation index (F_ST) was 0.08 for all studied markers, which is indicative of the low genetic differentiation of the Danubian horse population. Our analysis also confirmed the low level of inbreeding and heterozygous deficiency among the animals selected from the six paternal lineages of the Danubian horse. The present research could be helpful for the development of breeding and conservation programs for the Danubian horse, as well as for making informed decisions on the management of paternal lines.

Abstract

The Danubian horse, together with the Pleven and the Eastern Bulgarian horse breeds, is one of the modern breeds in Bulgaria. The objective of this study was to compare the genetic structure and genetic diversity of six paternal genealogical lineages of the Danubian horse breed (Zdravko, NONIUS XVII-30, Torpedo, Lider, Kalifa, and Hrabar). In total, 166 individuals from the six genealogical lines were investigated, based on 15 STR markers (short tandem repeats, also known as microsatellites). In total, 184 alleles were found in the six populations, using 15 microsatellite loci. The mean number of alleles, the effective number of alleles, and the polymorphic information content (PIC) values per locus were 12.28, 9.48, and 0.73, respectively. In a comparison of the allelic diversity among sire lineages, the highest genetic diversity (Na) was observed in Lider and Kalifa (14.60 ± 0.21), while the lowest value of this parameter was observed in the Zdravko lineage 4.20 ± 0.35. The largest genetic diversity was found in loci HMS3 and HMS7, with 13 alleles, and the smallest polymorphism was noted for the locus ASB17, with 10 alleles. The level of observed heterozygosity was in the range of 0.65 ± 0.069 for the Zdravko lineage to 0.93 ± 0.01 for the Torpedo lineage. The expected heterozygosity level range was from 0.57 ± 0.048 to 0.91 ± 0.01 for all horse lineages. Structure analysis revealed three main gene pools in the study population. The first pool included the Zdravko lineage; the second had the NONIUS XVII-30, Torpedo, Lider, and Kalifa lineages; and the third defined the Hrabar lineage, which was significantly differentiated from the other genealogical lineages.

Kinetics of the Equid Herpesvirus 2 and 5 Infections among Mares and Foals from Three Polish National Studs

Equid herpesvirus 2 (EHV-2) and 5 (EHV-5) are two γ-herpesviruses that are commonly detected from horses worldwide, based on several cross-sectional molecular surveys. Comparatively few studies examined the dynamics of γ-herpesvirus infection over time in a group of horses. The aim of the current study was to investigate the dynamics of EHV-2/5 infections among mares and their foals at three Polish national studs with different breeds of horses: Arabians, Thoroughbreds and Polish Konik horses. Nasal swabs were collected from each of 38 mare-foal pairs monthly for a period of 6 to 8 months. Virus-specific quantitative PCR assays were used to determine the viral load of EHV-2 and EHV-5 in each sample. All 76 horses sampled were positive for EHV-2 or EHV-5 on at least one sampling occasion. The majority (73/76, 96%) were infected with both EHV-2 and EHV-5. In general, the mean load of viral DNA was higher in samples from foals than from mares, but similar for EHV-2 and EHV-5 at most sampling occasions. There was, however, a considerable variability in the viral DNA load between samples collected at different times from the same foal, as well as between samples from different foals. The latter was more apparent for EHV-2 than for EHV-5. All foals became infected with both viruses early in life, before weaning, and remained positive on all, or most, subsequent samplings. The virus shedding by mares was more intermittent, indicating the existence of age-related differences. Overall, the data presented extend our knowledge of EHV-2/5 epidemiology among mares and foals.

Some Genetic and Environmental Effects on Equine Asthma in Polish Konik Horses

Simple Summary

Equine Asthma (EA) is a blanket term covering inflammatory diseases of the lower airways in horses. It includes mild-to-moderate equine asthma, which affects horses of any age, and severe equine asthma, which is typically seen in horses older than 7 years of. Relationships of the disease’s occurrence with sex or breed have been proven. However, some authors consider genetic background a predisposing factor, due to the fact that in some bloodlines a clinical form of asthma is more frequently observed. This indicates serious breeding and economic consequences. This study aimed to identify the factors affecting predisposition to severe equine asthma in the population of Polish Konik horses and some environmental and inbreeding effects on the disease. Generally, in the observed population, EA is negligibly affected by the factors analysed. Individual inbreeding effects on asthma were not confirmed by various statistical approaches, but significant maternal inbreeding effects were observed. These results are very important from the perspective of the currently implemented genetic resource conservation programme.

Abstract

Current knowledge of the genetic and environmental backgrounds of equine asthma seems to be insufficient, especially for primitive horse breeds. The main objectives of this study were to estimate the effects of sex, birth period, stud, parentage line and inbreeding on asthma morbidity in Polish Konik horses. Records of 274 horses (housed in two studs) were analysed. These animals were allocated to maternal and paternal lines. Individual inbreeding coefficients were extracted from the additive relationship matrix. Horses underwent diagnosis based on observation of the basic symptoms (high frequency of coughing and excessive nasal discharge). Subsequently, some horses (28 individuals) were clinically examined to confirm the earlier observations. Generally, no significant effects of parentage line on heaves morbidity were identified by the use of logistic regression, although the Pearson’s chi-squared test had shown that individuals of some maternal and paternal lines had a predisposition to severe equine asthma. It was concluded that the individual inbreeding level is not associated with the incidence of EA, but a significant effect of the maternal inbreeding coefficient may be observed. It was also found that there is some variability in the incidence of this disease between studs.

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.