Major inconsistencies of inferred population genetic structure estimated in a large set of domestic horse breeds using microsatellites

Population genomics provide insight into ancestral relationships and diversity of the feral horses of Theodore Roosevelt National Park

Theodore Roosevelt National Park (TRNP) manages a herd of feral horses (Equus caballus) which was present on the landscape prior to the establishment of the park. The population presents a unique scenario in that it has experienced fairly intensive and well-documented management since the park's establishment, including herd size reductions, intentional introduction of diversity, and subsequent attempts to remove introduced lineages. This provides an interesting case study on the genetic effects of diverse evolutionary forces on an isolated feral population. To explore the effects of these forces and clarify the relationship of this feral herd with other horses, we used genome-wide markers to examine the population structure of a combined dataset containing common established breeds. Using the Illumina Equine 70k BeadChip, we sampled SNPs across the genome for 118 TRNP horses and evaluated the inbreeding coefficient f and runs of homozygosity (RoH). To identify breed relationships, we compared 23 representative TRNP samples with 792 horses from 35 different breeds using genomic population structure analyses. Mean f of TRNP horses was 0.180, while the mean f for all other breeds in the dataset was 0.116 (SD 0.079). RoH analysis indicates that the TRNP population has experienced recent inbreeding in a timeframe consistent with their management. With Bayesian clustering, PCA, and maximum likelihood phylogeny, TRNP horses show genetic differentiation from other breeds, likely due to isolation, historical population bottlenecks, and genetic drift. However, maximum likelihood phylogeny places them with moderate confidence (76.8%) among draft breeds, which is consistent with the known history of breeds used on early North Dakota ranches and stallions subsequently introduced to the park herd. These findings will help resolve speculation about the origins of the herd and inform management decisions for the TRNP herd.

Defining Populations and Predicting Future Suitable Niche Space in the Geographically Disjunct, Narrowly Endemic Leafy Prairie-Clover (Dalea foliosa; Fabaceae)

Conservation actions for rare species are often based on estimates of population size and number, which are challenging to capture in natural systems. Instead, many definitions of populations rely on arbitrarily defined distances between occurrences, which is not necessarily biologically meaningful despite having utility from a conservation management perspective. Here, we introduce a case study using the narrowly endemic and highly geographically disjunct leafy prairie-clover (Dalea foliosa), for which we use nuclear microsatellite loci to assess the current delimitations of populations and management units across its entire known range. We model future potential suitable niche space for the species to assess how currently defined populations could fare under predicted changes in climate over the next 50 years. Our results indicate that genetic variation within the species is extremely limited, particularly so in the distal portions of its range (Illinois and Alabama). Within the core of its range (Tennessee), genetic structure is not consistent with populations as currently defined. Our models indicate that predicted suitable niche space may only marginally overlap with the geology associated with this species (limestone glades and dolomite prairies) by 2070. Additional studies are needed to evaluate the extent to which populations are ecologically adapted to local environments and what role this could play in future translocation efforts.

Clumppling: cluster matching and permutation program with integer linear programming

MOTIVATION

In the mixed-membership unsupervised clustering analyses commonly used in population genetics, multiple replicate data analyses can differ in their clustering solutions. Combinatorial algorithms assist in aligning clustering outputs from multiple replicates so that clustering solutions can be interpreted and combined across replicates. Although several algorithms have been introduced, challenges exist in achieving optimal alignments and performing alignments in reasonable computation time.

RESULTS

We present Clumppling, a method for aligning replicate solutions in mixed-membership unsupervised clustering. The method uses integer linear programming for finding optimal alignments, embedding the cluster alignment problem in standard combinatorial optimization frameworks. In example analyses, we find that it achieves solutions with preferred values of a desired objective function relative to those achieved by Pong and that it proceeds with less computation time than Clumpak. It is also the first method to permit alignments across replicates with multiple arbitrary values of the number of clusters K.

AVAILABILITY AND IMPLEMENTATION

Clumppling is available at https://github.com/PopGenClustering/Clumppling.

Structure and genetic variability of the Costa Rican Paso horse

The Costa Rican Paso Horse (CPC) is a breed developed in Costa Rica. The objectives were to estimate the genetic structure and evaluate the levels of genetic variability of the population. The genotypes of 14 microsatellites in 3654 records (2052 females and 1602 males) were analyzed. Expected (He) and observed (Ho) heterozygosity, polymorphic information content (CIP), fixation index (FIS), Shannon index, as well as Hardy-Weinberg disequilibrium (DHW) were evaluated. Kinship relationships (Rij) were estimated throughout the entire population. The effective population size (Ne) was calculated, alternating allele frequencies less than 0.05, 0.02 and 0.01. The Bayesian clustering study was carried out to infer how many lines are appropriate from the analysis of genotypes using multiple loci. The number of alleles per locus ranged from 7 to 17, with an average value of 9.6; nine loci presented DHW (P < 0.05); two loci presented negative FIS values, the same as Ho > He; the average of CIP, Ho and He was 0.254, 0.756 and 0.785, respectively. At the 12 loci where He > Ho, the differences ranged from 0.002 to 0.341 (0.036 on average). For Ne, the estimates were 201.9, 230.1, and 241.5. In the Rij, 54.86% of the estimates were in the interval of 0.01 to 77.7%. The number of lines that define the population corresponds to three, with an approximate composition of 33.1%, 32.4% and 34.5%, respectively. The CPC, as a subdivided population with DHW and a reduction in heterozygotes may be associated with possible Wahlund effects. Keywords: Wahlund effect, equines, genetic markers, synthetic breed, Hardy Weinberg.

Are rabid raccoons (Procyon lotor) ready for the rapture? Determining the geographic origin of rabies virus-infected raccoons using RADcapture and microhaplotypes

North America is recognized for the exceptional richness of rabies virus (RV) wildlife reservoir species. Management of RV is accomplished through vaccination targeting mesocarnivore reservoir populations, such as the raccoon (Procyon lotor) in Eastern North America. Raccoons are a common generalist species, and populations may reach high densities in developed areas, which can result in contact with humans and pets with potential exposures to the raccoon variant of RV throughout the eastern United States. Understanding the spatial movement of RV by raccoon populations is important for monitoring and refining strategies supporting the landscape-level control and local elimination of this lethal zoonosis. We developed a high-throughput genotyping panel for raccoons based on hundreds of microhaplotypes to identify population structure and genetic diversity relevant to rabies management programs. Throughout the eastern United States, we identified hierarchical population genetic structure with clusters that were connected through isolation-by-distance. We also illustrate that this genotyping approach can be used to support real-time management priorities by identifying the geographic origin of a rabid raccoon that was collected in an area of the United States that had been raccoon RV-free for 8 years. The results from this study and the utility of the microhaplotype panel and genotyping method will provide managers with information on raccoon ecology that can be incorporated into future management decisions.

Genetic Diversity of Kazakhstani Equus caballus (Linnaeus, 1758) Horse Breeds Inferred from Microsatellite Markers

Understanding the genetic diversity and structure of domesticated horse (Equus caballus) populations is critical for long-term herd management and breeding programs. This study examines 435 horses from Kazakhstan, covering seven groups in three geographic areas using 11 STR markers. Identified are 136 alleles, with the mean number of alleles per locus ranging from 9 to 19. VHL20 is the most variable locus across groups, while loci HTG4, AHT4, AHT5, HTG7, and HMS3 are variable in most populations. The locus AHT5 in the Emba population shows the highest frequency of rare alleles, while the lowest frequency, 0.005, is observed in the Kulandy population. All loci were highly informative for the Kazakhstani populations of E. caballus, with PIC values higher than 0.5. Pairwise variations in Wright's FST distances show that the examined varieties have little genetic differentiation (0.05%), indicating a high degree of admixture and a continuing lineage sorting process. Phylogenetic and population structure analyses reveal three major clusters of Kazakh horses, representing (I) the Uralsk population of the Kushum breed and the monophyly of two groups: (II) the Kozhamberdy population of the Mugalzhar breed, and (III) the Mugalzhar-Kushum breed populations. Kazakhstani horse populations, while being regionally isolated, were recently in contact with each other.

Spatial and temporal patterns in the population genomics of the European cockchafer Melolontha melolontha in the Alpine region

The European cockchafer Melolontha melolontha is an agricultural pest in many European countries. Populations have a synchronized 3 or 4 years life cycle, leading to temporally isolated populations. Despite the economic importance and availability of comprehensive historical as well as current records on cockchafer occurrence, population genomic analyses of M. melolontha are missing. For example, the effects of geographic separation caused by the mountainous terrain of the Alps and of temporal isolation on the genomic structure of M. melolontha still remain unknown. To address this gap, we genotyped 475 M. melolontha adults collected during 3 years from 35 sites in a central Alpine region. Subsequent population structure analyses discriminated two main genetic clusters, i.e., the South Tyrol cluster including collections located southeast of the Alpine mountain range, and a northwestern alpine cluster with all the other collections, reflecting distinct evolutionary history and geographic barriers. The "passo di Resia" linking South and North Tyrol represented a regional contact zone of the two genetic clusters, highlighting genomic differentiation between the collections from the northern and southern regions. Although the collections from northwestern Italy were assigned to the northwestern alpine genetic cluster, they displayed evidence of admixture with the South Tyrolean genetic cluster, suggesting shared ancestry. A linear mixed model confirmed that both geographic distance and, to a lower extent, also temporal isolation had a significant effect on the genetic distance among M. melolontha populations. These effects may be attributed to limited dispersal capacity and reproductive isolation resulting from synchronized and non-synchronized swarming flights, respectively. This study contributes to the understanding of the phylogeography of an organism that is recognized as an agricultural problem and provides significant information on the population genomics of insect species with prolonged temporally shifted and locally synchronized life cycles.

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.

Crimp: An efficient tool for summarizing multiple clusterings in population structure analysis and beyond

When a data set is repeatedly clustered using unsupervised techniques, the resulting clusterings, even if highly similar, may list their clusters in different orders. This so-called 'label-switching' phenomenon obscures meaningful differences between clusterings, complicating their comparison and summary. The problem often arises in the context of population structure analysis based on multilocus genotype data. In this field, a variety of popular tools apply model-based clustering, assigning individuals to a prespecified number of ancestral populations. Since such methods often involve stochastic components, it is a common practice to perform multiple replicate analyses based on the same input data and parameter settings. Available postprocessing tools allow to mitigate label switching, but leave room for improvements, in particular, regarding large input data sets. In this work, I present Crimp, a lightweight command-line tool, which offers a relatively fast and scalable heuristic to align clusters across replicate clusterings consisting of the same number of clusters. For small problem sizes, an exact algorithm can be used as an alternative. Additional features include row-specific weights, input and output files similar to those of CLUMPP (Jakobsson & Rosenberg, 2007) and the evaluation of a given solution in terms of CLUMPP as well as its own objective functions. Benchmark analyses show that Crimp, especially when applied to larger data sets, tends to outperform alternative tools considering runtime requirements and various quality measures. While primarily targeting population structure analysis, Crimp can be used as a generic tool to correct multiple clusterings for label switching. This facilitates their comparison and allows to generate an averaged clustering. Crimp's computational efficiency makes it even applicable to relatively large data sets while offering competitive solution quality.

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.

Genetic Differentiation among Livestock Breeds—Values for Fst

Simple Summary

The degree of relationship among livestock breeds can be quantified by the F_st statistic, which measures the extent of genetic differentiation between them. An F_st value of 0.1 has often been taken as indicating that two breeds are indeed genetically distinct, but this concept has not been evaluated critically. Here, F_st values have been collated for the six major livestock species: cattle, sheep, goats, pigs, horses, and chickens. These values are remarkably variable both within and between species, demonstrating that F_st > 0.1 is not a reliable criterion for breed distinctiveness. However, the large body of F_st data accumulated in the last 20–30 years represents an untapped database that could contribute to the development of interdisciplinary research involving livestock breeds.

Abstract

(1) Background: The F_st statistic is widely used to characterize between-breed relationships. F_st = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding F_st values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed F_st were as follows, cattle: microsatellite 0.06–0.12, SNP 0.08–0.15; sheep: microsatellite 0.06–0.10, SNP 0.06–0.17; horses: microsatellite 0.04–0.11, SNP 0.08–0.12; goats: microsatellite 0.04–0.14, SNP 0.08–0.16; pigs: microsatellite 0.06–0.27, SNP 0.15–0.22; chickens: microsatellite 0.05–0.28, SNP 0.08–0.26. (4) Conclusions: (1) Large amounts of F_st data are available for a substantial proportion of the world’s livestock breeds, (2) the value for between-breed F_st of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated F_st data may have value for interdisciplinary research.

Genetic Variability Trend of Lusitano Horse Breed Reared in Italy

The Lusitano Horse (LH) originates from Portugal, but is reared worldwide. Since 1994, the University of Milan has routinely tested the LHs bred in Italy for parentage control. This study aims to assess the genetic variability of the LH reared in Italy using 16 microsatellites markers. Moreover, the genetic variability changes over the years in the total population (n.384) and in unrelated horses (n.47) were evaluated. Horses were grouped according to their date of birth (1975–1990, 1991–2000, 2001–2010, 2010–2019). Standard genetic diversity parameters, including observed (Ho) and expected (He) heterozygosity, Hardy-Weinberg equilibrium (HWE; P-Val), allelic richness, and inbreeding coefficient (Fis) were estimated. In the whole period, the total population showed Ho as high as 0.69, low Fis (0.057), and imbalance for HWE. When considering the unrelated horses, Ho was seen to increase over time (from 0.594 in 1975–1990 to 0.68 in 2010–2019) and frequencies were in HWE, again having low and decreasing values of Fis (from 0.208 in 1975–1990 to 0.019 in 2010–2019). Bottleneck analysis excluded a recent population decline. Principal Coordinate Analysis at the individual level defined two clusters, the major cluster including all the most recent horses. An increasing number of dams (156% more from 2001–2010 to 2011–2019) supports the good variability recorded in the population so far. However, the high number of foals (77.2%) sired by only four stallions in recent years suggests caution in the choice of the sires for the future.

The impact of estimator choice: Disagreement in clustering solutions across K estimators for Bayesian analysis of population genetic structure across a wide range of empirical data sets

The software program STRUCTURE is one of the most cited tools for determining population structure. To infer the optimal number of clusters from STRUCTURE output, the ΔK method is often applied. However, a recent study relying on simulated microsatellite data suggested that this method has a downward bias in its estimation of K and is sensitive to uneven sampling. If this finding holds for empirical data sets, conclusions about the scale of gene flow may have to be revised for a large number of studies. To determine the impact of method choice, we applied recently described estimators of K to re-estimate genetic structure in 41 empirical microsatellite data sets; 15 from a broad range of taxa and 26 from one phylogenetic group, coral. We compared alternative estimates of K (Puechmaille statistics) with traditional (ΔK and posterior probability) estimates and found widespread disagreement of estimators across data sets. Thus, one estimator alone is insufficient for determining the optimal number of clusters; this was regardless of study organism or evenness of sampling scheme. Subsequent analysis of molecular variance (AMOVA) did not necessarily clarify which clustering solution was best. To better infer population structure, we suggest a combination of visual inspection of STRUCTURE plots and calculation of the alternative estimators at various thresholds in addition to ΔK. Disagreement between traditional and recent estimators may have important biological implications, such as previously unrecognized population structure, as was the case for many studies reanalysed here.

The future of endangered crayfish in light of protected areas and habitat fragmentation

The long-term survival of a species requires, among other things, gene flow between populations. Approaches for the evaluation of fragmentation in the frame of freshwater habitats consider only a small amount of the information that combined demography and geography are currently able to provide. This study addresses two species of Austropotamobius crayfish in the light of population genetics, spatial ecology and protected areas of the Carpathians. Advancing the classical approaches, we defined ecological distances upon the rasterised river network as a surrogate of habitat resistance to migration, quantifying the deviations from the species´ suitability range for a set of relevant geospatial variables in each cell of the network. Molecular analyses revealed the populations of the two Austropotamobius crayfish species are clearly distinct, lacking hybridisation. Comparing pairs of populations, we found, in some cases, a strong disagreement regarding genetic and ecological distances, potentially due to human-mediated translocations or the geophysical phenomena of regressive erosion, which may have led to unexpected colonisation routes. Protected areas were found to offer appropriate local habitat conditions but failed to ensure connectivity. The methodology applied in this study allowed us to quantify the contribution of each geospatial (environmental) variable to the overall effect of fragmentation, and we found that water quality was the most important variable. A multilevel approach proved to reveal a better understanding of drivers behind the distribution patterns, which can lead to more adequate conservation measures.