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.

BACKGROUND: The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds.
RESULTS: HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20).
CONCLUSIONS: Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Linking genetic diversity to extinction is a common goal in genomic studies. Recently, a debate has arisen regarding the importance of genetic variation in conservation as some studies have failed to find associations between genome-wide genetic diversity and extinction risk. However, only rarely are genetic diversity and fitness measured together in the wild, and typically demographic history and environment are ignored. It is therefore difficult to infer whether a lack of an association is real or obscured by confounding factors. To address these shortcomings, we analyzed genetic data from 7,501 individuals with extinction data from 279 meadows and mortality of 1,742 larval nests in a butterfly metapopulation. We found a strong negative association between genetic diversity and extinction when considering only heterozygosity in models. However, this association disappeared when accounting for ecological covariates, suggesting a confounding between demography and genetics and a more complex role for heterozygosity in extinction risk. Modeling interactions between heterozygosity and demographic variables revealed that associations between extinction and heterozygosity were context-dependent. For example, extinction declined with increasing heterozygosity in large, but not currently small populations, although negative associations between heterozygosity, extinction, and mortality were detected in small populations with a recent history of decline. We conclude that low genetic diversity is an important predictor of extinction, predicting >25% increase in extinction beyond ecological factors in certain contexts. These results highlight that inferences about the importance of genetic diversity for population viability should not rely on genomic data alone but require investments in obtaining demographic and environmental data from natural populations.

Cultivated potatoes are incredibly diverse, ranging from diploid to pentaploid and encompass four different species. They are adapted to disparate environments and conditions and carry unique alleles for resistance to pests and pathogens. Describing how diversity is partitioned within and among these populations is essential to understanding the potato genome and effectively utilizing landraces in breeding. This task is complicated by the difficulty of making comparisons across cytotypes and extensive admixture within section petota. We genotyped 730 accessions from the US Potato genebank including wild diploids and cultivated diploids and tetraploids using Genotype-by-sequencing. This data set allowed us to interrogate population structure and diversity as well as generate core subsets which will support breeders in efficiently screening genebank material for biotic and abiotic stress resistance alleles. We found that even controlling for ploidy, tetraploid material exhibited higher observed and expected heterozygosity than diploid accessions. In particular group chilotanum material was the most heterozygous and the only taxa not to exhibit any inbreeding. This may in part be because group chilotanum has a history of introgression not just from wild species, but landraces as well. All group chilotanum, exhibits introgression from group andigenum except clones from Southern South America near its origin, where the two groups are not highly differentiated. Moving north, we do not observe evidence for the same level of admixture back into group andigenum. This suggests that extensive history of admixture is a particular characteristic of chilotanum.

Coexistence of multiple strains of a pathogen in a host population can present significant challenges to vaccine development or treatment efficacy. Here we discuss a novel mechanism that can increase rates of long-lived strain polymorphism, rooted in the presence of social structure in a host population. We show that social preference of interaction, in conjunction with differences in immunity between host subgroups, can exert varying selection pressure on pathogen strains, creating a balancing mechanism that supports stable viral coexistence, independent of other known mechanisms. We use population genetic models to study rates of pathogen heterozygosity as a function of population size, host population composition, mutant strain fitness differences and host social preferences of interaction. We also show that even small periodic epochs of host population stratification can lead to elevated strain coexistence. These results are robust to varying social preferences of interaction, overall differences in strain fitnesses, and spatial heterogeneity in host population composition. Our results highlight the role of host population social stratification in increasing rates of pathogen strain diversity, with effects that should be considered when designing policies or treatments with a long-term view of curbing pathogen evolution.

Rails are a phenotypically diverse family of birds that includes 130 species and displays a wide distribution around the world. Here we present annotated genome assemblies for four rails from Aotearoa New Zealand: two native volant species, pūkeko Porphyrio melanotus and mioweka Gallirallus philippensis, and two endemic flightless species takahē Porphyrio hochstetteri and weka Gallirallus australis. Using the sequence read data, heterozygosity was found to be lowest in the endemic flightless species and this probably reflects their relatively small populations. The quality checks and comparison with other rallid genomes showed that the new assemblies were of good quality. This study significantly increases the number of available rallid genomes and will enable future genomic studies on the evolution of this family.

There are two primary measures of the amount of genetic variation in a population at a locus: heterozygosity and the number of alleles. Effective population size (N e) provides both an expectation of the amount of heterozygosity in a population at drift-mutation equilibrium and the rate of loss of heterozygosity because of genetic drift. In contrast, the number of alleles in a population at drift-mutation equilibrium is a function of both N e and census size (N C). In addition, populations with the same N e can lose allelic variation at very different rates. Allelic variation is generally much more sensitive to bottlenecks than heterozygosity. Expressions used to adjust for the effects of violations of the ideal population on N e do not provide good predictions of the loss of allelic variation. These effects are much greater for loci with many alleles, which are often important for adaptation. We show that there is a linear relationship between the reduction of N C and the corresponding reduction of the expected number of alleles at drift-mutation equilibrium. This makes it possible to predict the expected effect of a bottleneck on allelic variation. Heterozygosity provides good estimates of the rate of adaptive change in the short-term, but allelic variation provides important information about long-term adaptive change. The guideline of long-term N e being greater than 500 is often used as a primary genetic metric for evaluating conservation status. We recommend that this guideline be expanded to take into account allelic variation as well as heterozygosity.

Despite their critical roles in genetic sex determination, sex chromosomes remain unknown in many non-model organisms, especially those having recently evolved sex-linked regions (SLRs). These evolutionarily young and labile sex chromosomes are important for understanding early sex chromosome evolution but are difficult to identify due to the lack of Y/W degeneration and SLRs limited to small genomic regions. Here, we present SLRfinder, a method to identify candidate SLRs using linkage disequilibrium (LD) clustering, heterozygosity and genetic divergence. SLRfinder does not rely on specific sequencing methods or a specific type of reference genome (e.g., from the homomorphic sex). In addition, the input of SLRfinder does not require phenotypic sexes, which may be unknown from population sampling, but sex information can be incorporated and is necessary to validate candidate SLRs. We tested SLRfinder using various published datasets and compared it to the local principal component analysis (PCA) method and the depth-based method Sex Assignment Through Coverage (SATC). As expected, the local PCA method could not be used to identify unknown SLRs. SATC works better on conserved sex chromosomes, whereas SLRfinder outperforms SATC in analysing labile sex chromosomes, especially when SLRs harbour inversions. Power analyses showed that SLRfinder worked better when sampling more populations that share the same SLR. If analysing one population, a relatively larger sample size (around 50) is needed for sufficient statistical power to detect significant SLR candidates, although true SLRs are likely always top-ranked. SLRfinder provides a novel and complementary approach for identifying SLRs and uncovering additional sex chromosome diversity in nature.

vonHoldt et al. ((2024), Molecular Ecology, 33, e17231) (vH24) used low-coverage (average ~ 7X read depth) restriction site-associated DNA sequence data to estimate individual inbreeding and heterozygosity, and recent effective population size (Ne), in Great Lakes (GL) and Northern Rocky Mountain (RM) wolves. They concluded that RM heterozygosity rapidly declined between 1991 and 2020, and that Ne declined substantially in GL and RM over the last 50 generations. Here, we evaluate the effects of low sequence coverage and sampling strategy on vH24\'s findings and provide general recommendations for using sequence data to evaluate inbreeding, heterozygosity and Ne. Low-coverage sequencing resulted in downwardly biased estimates of individual inbreeding and heterozygosity, and an erroneous large temporal decline in RM heterozygosity due to declining read depth through time. Additionally, vH24\'s sampling strategy-which combined individuals from several genetically differentiated populations and across at least eight wolf generations-is expected to have resulted in severe downward bias in estimates of recent Ne for RM. We recommend using high-coverage sequence data ( ≥ $$ \\ge $$ 15-20X) to estimate inbreeding and heterozygosity. Carefully filtering individuals, loci and genotypes, and using genotype imputation or likelihoods can help to minimise bias when low-coverage sequence data must be used. For estimation of contemporary Ne, the marginal benefits of using more than 103-104 loci are small, so aggressive filtering of loci with low average read depth potentially can retain most individuals without sacrificing much precision. Individuals are relatively more valuable than loci because analyses of contemporary Ne should focus on roughly single-generation samples from local breeding populations.

Factors that increase reproductive variance among individuals act to reduce effective population size (Ne), which accelerates the loss of genetic diversity and decreases the efficacy of purifying selection. These factors include sexual cannibalism, offspring investment and mating system. Pre-copulatory sexual cannibalism, where the female consumes the male prior to mating, exacerbates this effect. We performed comparative transcriptomics in two spider species, the cannibalistic Trechaleoides biocellata and the non-cannibalistic T. keyserlingi, to generate genomic evidence to support these predictions. First, we estimated heterozygosity and found that genetic diversity is relatively lower in the cannibalistic species. Second, we calculated dN/dS ratios as a measure of purifying selection; a higher dN/dS ratio indicated relaxed purifying selection in the cannibalistic species. These results are consistent with the hypothesis that sexual cannibalism impacts operational sex ratio and demographic processes, which interact with evolutionary forces to shape the genetic structure of populations. However, other factors such as the mating system and life-history traits contribute to shaping Ne. Comparative analyses across multiple contrasting species pairs would be required to disentangle these effects. Our study highlights that extreme behaviours such as pre-copulatory cannibalism may have profound eco-evolutionary effects.