Inferring Population Structure and Admixture Proportions in Low-Depth NGS Data

Development of SNP Panels from Low-Coverage Whole Genome Sequencing (lcWGS) to Support Indigenous Fisheries for Three Salmonid Species in Northern Canada

Single nucleotide polymorphism (SNP) panels are powerful tools for assessing the genetic population structure and dispersal of fishes and can enhance management practices for commercial, recreational and subsistence mixed-stock fisheries. Arctic Char (Salvelinus alpinus), Brook Trout (Salvelinus fontinalis) and Lake Whitefish (Coregonus clupeaformis) are among the most harvested and consumed fish species in Northern Indigenous communities in Canada, contributing significantly to food security, culture, tradition and economy. However, genetic resources supporting Indigenous fisheries have not been widely accessible to northern communities (e.g. Inuit, Cree, Dene). Here, we developed Genotyping-in-Thousands by sequencing (GT-seq) panels for population assignment and mixed-stock analyses of three salmonids, to support fisheries stewardship or co-management in Northern Canada. Using low-coverage Whole Genome Sequencing data from 418 individuals across source populations in Cambridge Bay (Nunavut), Great Slave Lake (Northwest Territories), James Bay (Québec) and Mistassini Lake (Québec), we developed a bioinformatic SNP filtering workflow to select informative SNP markers from genotype likelihoods. These markers were then used to design GT-seq panels, thus enabling high-throughput genotyping for these species. The three GT-seq panels yielded an average of 413 autosomal loci and were validated using 525 individuals with an average assignment accuracy of 83%. Thus, these GT-seq panels are powerful tools for assessing population structure and quantifying the relative contributions of populations/stocks in mixed-stock fisheries across multiple regions. Interweaving genomic data derived from these tools with Traditional Ecological Knowledge will ensure the sustainable harvest of three culturally important salmonids in Indigenous communities, contributing to food security programmes and the economy in Northern Canada.

Genomic Signatures of Adaptation Across a Precipitation Gradient From Niche Centre to Niche Edge

Evaluating the potential for species to adapt to changing climates relies on understanding current patterns of adaptive variation and selection, which might vary in intensity across a species' niche, hence affecting our inference of where adaptation might be most important in the future. Here, we investigate the genetic basis of adaptation in Lactuca serriola along a steep precipitation gradient in Israel approaching the species' arid niche limit and use candidate loci to inform predictions of its past and future adaptive evolution. Environmental association analyses combined with generalised dissimilarity models revealed 108 candidate genes showing nonlinear shifts in allele frequencies across the gradient, with 66% of these genes under strong selection near the dry niche edge. We detected selection acting on genes with separate suites of biological functions, specifically related to phenology and responses to environmental stressors, including osmotic stress, at the dry niche edge, and related to biotic interactions and defence closer to the niche centre. The adaptive genetic composition of populations, as inferred through polygenic risk scores, point to intensified selection operating towards the dry niche edge. However, inference of past and future evolutionary change predicts larger adaptive shifts occurring in the mesic part of the range, which is most affected by climate change. Our study reveals that adaptive shifts in response to climate change can be heterogeneous across a species' range and not necessarily strongest near its niche edge.

Ecological divergence of sibling allopolyploid marsh orchids is associated with species specific plasticity and distinct fungal communities

Phenotypic plasticity, the dynamic adjustment of traits to environmental variations, is crucial for enabling species to exploit broader niches and withstand suboptimal conditions. This adaptability is particularly relevant for newly formed allopolyploids, which possess redundant gene copies and must become established in diverse environments distinct from their parents and other relatives. By evaluating gene expression and root mycobiome among two ecologically divergent sibling allopolyploid marsh orchids (Dactylorhiza majalis and D. traunsteineri) in reciprocal transplants at localities where both species are native, we aimed to understand the drivers of species persistence in the face of interspecific gene flow. Despite consistent abiotic differences characterising the alternative environments at each locality, the majority of gene expression differences between the allopolyploids appears to be plastic. Ecologically relevant processes, such as photosynthesis and transmembrane transport, include some genes that are differentially expressed between the two orchids regardless of the environment, while others change their activity plastically in one species or the other. This suggests that although plasticity helps define the specific ecological range of each sibling allopolyploid, it also mediates gene flow between them, thereby preventing differentiation. Extending our investigations to the root mycobiome, we uncover more diverse fungal communities for either species when grown in the environment with nutrient-poor soils, indicating that both abiotic and biotic factors drive the distribution of sibling marsh orchids. Altogether, our results indicate that plasticity can simultaneously promote diversification and homogenisation of lineages, influencing the establishment and persistence of recurrently formed allopolyploid species.

Circumpolar Population Structure, Diversity and Recent Evolutionary History of the Bearded Seal in Relation to Past and Present Icescapes

The Arctic environment plays a critical role in the global climate system and marine biodiversity. The region's ice-covered expanses provide essential breeding and feeding grounds for a diverse assemblage of marine species, who have adapted to thrive in these harsh conditions and consequently are under threat from global warming. The bearded seal (Erignathus barbatus), including two subspecies (E. barbatus nauticus-Pacific and E. barbatus barbatus-Atlantic), is an ice-obligate Arctic species using sea ice for many aspects of its life history, rendering it particularly vulnerable to sea ice loss. It is one of the least studied and hence enigmatic of the Arctic marine mammals, with little knowledge regarding genetic structure, diversity, adaptations, and demographic history, consequently hampering management and conservation efforts. Here, we sequenced 70 whole nuclear genomes from across most of the species' circumpolar range, finding significant genetic structure between the Pacific and the Atlantic subspecies, which diverged during the Penultimate Glacial Period (~200 KYA). Remarkably, we found fine-scale genetic structure within both subspecies, with at least two distinct populations in the Pacific and three in the Atlantic. We hypothesise sea-ice dynamics and bathymetry had a prominent role in shaping bearded seal genetic structure and diversity. Our analyses of highly differentiated genomic regions can be used to complement the health, physiological, and behavioural research needed to conserve this species. In addition, we provide recommendations for management units that can be used to more specifically assess climatic and anthropogenic impacts on bearded seal populations.

Kinship clustering within an ecologically diverse killer whale metapopulation

Metapopulation dynamics can be shaped by foraging ecology, and thus be sensitive to shifts in prey availability. Genotyping 204 North Atlantic killer whales at 1346 loci, we investigated whether spatio-temporal population structuring is linked to prey type and distribution. Using population-based methods (reflecting evolutionary means), we report a widespread metapopulation connected across ecological groups based upon nuclear genome SNPs, yet spatial structuring based upon mitogenome haplotypes. These contrasting patterns of markers with maternal and biparental inheritance are consistent with matrilineal site fidelity and philopatry, and male-mediated gene flow among demes. Connectivity between fish-eating and 'mixed-diet' (eating both fish and mammal prey) killer whales, marks a deviation within a species renowned for its marked structure associated with ecology. However, relatedness estimates suggest distinct spatial clusters, and heterogeneity in recent gene flow between them. The contrasting patterns between inference of structure and inference of relatedness suggest that gene flow has been partially restricted over the past two to three generations (50-70 years). This coincides with the collapse of North Atlantic herring stocks in the late 1960s and the subsequent cessation of their seasonal connectivity. Statistically significant association between diet types and assignment of Icelandic killer whales to relatedness-based clusters indicated limited gene flow was maintained through Icelandic 'mixed-diet' whales when herring-mediated connectivity was diminished. Thus, conservation of dietary variation within this metapopulation is critical to ensure genetic health. Our study highlights the role of resource dynamics and foraging ecology in shaping population structure and emphasises the need for transnational management of this widespread migratory species and its prey.

Common Ancestry of the Id Locus: Chromosomal Rearrangement and Polygenic Possibilities

The diversity in dermal pigmentation and plumage color among domestic chickens is striking, with Black Bone Chickens (BBC) particularly notable for their intense melanin hyperpigmentation. This unique trait is driven by a complex chromosomal rearrangement on chromosome 20 at the Fm locus, resulting in the overexpression of the EDN3 (a gene central to melanocyte regulation). In contrast, the inhibition of dermal pigmentation is regulated by the Id locus. Although prior studies using genetic crosses, GWAS, and gene expression analysis have investigated the genetic underpinnings of the Id locus, its precise location and functional details remain elusive. Our study aims to precisely locate the Id locus, identify associated chromosomal rearrangements and candidate genes influencing dermal pigmentation, and examine the ancestral status of the Id locus in BBC breeds. Using public genomic data from BBC and non-BBC breeds, we refined the Id locus to a ~1.6 Mb region that co-localizes with Z amplicon repeat units at the distal end of the q-arm of chromosome Z within a 10.36 Mb inversion in Silkie BBC. Phylogenetic and population structure analyses reveal that the Id locus shares a common ancestry across all BBC breeds, much like the Fm locus. Selection signatures and highly differentiated BBC-specific SNPs within the MTAP gene position it as the prime candidate for the Id locus with CCDC112 and additional genes, suggesting a possible polygenic nature. Our results suggest that the Id locus is shared among BBC breeds and may function as a supergene cluster in shank and dermal pigmentation variation.

Ancient genomics and the origin, dispersal, and development of domestic sheep

The origins and prehistory of domestic sheep (Ovis aries) are incompletely understood; to address this, we generated data from 118 ancient genomes spanning 12,000 years sampled from across Eurasia. Genomes from Central Türkiye ~8000 BCE are genetically proximal to the domestic origins of sheep but do not fully explain the ancestry of later populations, suggesting a mosaic of wild ancestries. Genomic signatures indicate selection by ancient herders for pigmentation patterns, hornedness, and growth rate. Although the first European sheep flocks derive from Türkiye, in a notable parallel with ancient human genome discoveries, we detected a major influx of Western steppe-related ancestry in the Bronze Age.

Multiple, Single Trait GWAS and Supervised Machine Learning Reveal the Genetic Architecture of Fraxinus excelsior Tolerance to Ash Dieback in Europe

Common ash (Fraxinus excelsior) is under intensive attack from the invasive alien pathogenic fungus Hymenoscyphus fraxineus, causing ash dieback at epidemic levels throughout Europe. Previous studies have found significant genetic variation among genotypes in ash dieback susceptibility and that host phenology, such as autumn yellowing, is correlated with susceptibility of ash trees to H. fraxineus; however, the genomic basis of ash dieback tolerance in F. excelsior requires further investigation. Here, we integrate quantitative genetics based on multiple replicates and genome-wide association analyses with machine learning to reveal the genetic architecture of ash dieback tolerance and of phenological traits in F. excelsior populations in six European countries (Austria, Denmark, Germany, Ireland, Lithuania, Sweden). Based on phenotypic data of 486 F. excelsior replicated genotypes we observed negative genotypic correlations between crown damage caused by ash dieback and intensity of autumn leaf yellowing within multiple sampling sites. Our results suggest that the examined traits are polygenic and using genomic prediction models, with ranked single nucleotide polymorphisms (SNPs) based on GWAS associations as input, a large proportion of the variation was predicted by unlinked SNPs. Based on 100 unlinked SNPs, we can predict 55% of the variation in disease tolerance among genotypes (as phenotyped in genetic trials), increasing to a maximum of 63% when predicted from 9155 SNPs. In autumn leaf yellowing, 52% of variation is predicted by 100 unlinked SNPs, reaching a peak of 72% using 3740 SNPs. Based on feature permutations within genomic prediction models, a total of eight nonsynonymous SNPs linked to ash dieback crown damage and autumn leaf yellowing (three and five SNPs, respectively) were identified, these were located within genes related to plant defence (pattern triggered immunity, pathogen detection) and phenology (regulation of flowering and seed maturation, auxin transport). We did not find an overlap between genes associated with crown damage level and autumn leaf yellowing. Hence, our results shed light on the difference in the genomic basis of ADB tolerance and autumn leaf yellowing despite these two traits being correlated in quantitative genetic analysis. Overall, our methods show the applicability of genomic prediction models when combined with GWAS to reveal the genomic architecture of polygenic disease tolerance enabling the identification of ash dieback tolerant trees for breeding or conservation purposes.

Rapid Parallel Adaptation in Distinct Invasions of Ambrosia Artemisiifolia Is Driven by Large-Effect Structural Variants

When introduced to multiple distinct ranges, invasive species provide a compelling natural experiment for understanding the repeatability of adaptation. Ambrosia artemisiifolia is an invasive, noxious weed, and chief cause of hay fever. Leveraging over 400 whole-genome sequences spanning the native-range in North America and 2 invasions in Europe and Australia, we inferred demographically distinct invasion histories on each continent. Despite substantial differences in genetic source and effective population size changes during introduction, scans of both local climate adaptation and divergence from the native-range revealed genomic signatures of parallel adaptation between invasions. Disproportionately represented among these parallel signatures are 37 large haploblocks-indicators of structural variation-that cover almost 20% of the genome and exist as standing genetic variation in the native-range. Many of these haploblocks are associated with traits important for adaptation to local climate, like size and the timing of flowering, and have rapidly reformed native-range clines in invaded ranges. Others show extreme frequency divergence between ranges, consistent with a response to divergent selection on different continents. Our results demonstrate the key role of large-effect standing variants in rapid adaptation during range expansion, a pattern that is robust to diverse invasion histories.

Modeling Biases from Low-Pass Genome Sequencing to Enable Accurate Population Genetic Inferences

Low-pass genome sequencing is cost-effective and enables analysis of large cohorts. However, it introduces biases by reducing heterozygous genotypes and low-frequency alleles, impacting subsequent analyses such as model-based demographic history inference. Several approaches exist for inferring an unbiased allele frequency spectrum (AFS) from low-pass data, but they can introduce spurious noise into the AFS. Rather than correcting the AFS, here, we developed an approach that incorporates low-pass biases into the demographic modeling and directly analyzes the AFS from low-pass data. Our probabilistic model captures biases from the Genome Analysis Toolkit multisample calling pipeline, and we implemented it in the population genomic inference software dadi. We evaluated the model using simulated low-pass datasets and found that it alleviated low-pass biases in inferred demographic parameters. We further validated the model by downsampling 1000 Genomes Project data, demonstrating its effectiveness on real data. Our model is widely applicable and substantially improves model-based inferences from low-pass population genomic data.

Genomic Insights Into Red Squirrels in Scotland Reveal Loss of Heterozygosity Associated With Extreme Founder Effects

Remnant populations of endangered species often have complex demographic histories associated with human impact. This can present challenges for conservation as populations modified by human activity may require bespoke management. The Eurasian red squirrel, Sciurus vulgaris (L., 1758), is endangered in the UK. Scotland represents a key stronghold, but Scottish populations have been subjected to intense anthropogenic influence, including widespread extirpations, reintroductions and competition from an invasive species. This study examined the genetic legacy of these events through low coverage whole-genome resequencing of 106 red squirrels. Previously undetected patterns of population structure and gene flow were uncovered. One offshore island, four mainland Scottish populations, and a key east-coast migration corridor were observed. An abrupt historical population bottleneck, related to extreme founder effects, has led to a severe and prolonged depression in genome-wide heterozygosity, which is amongst the lowest reported for any species. Current designated red squirrel conservation stronghold locations do not encompass all existing diversity. These findings highlight the genetic legacies of past anthropogenic influence on long-term diversity in endangered taxa. Continuing management interventions and regular genetic monitoring are recommended to safeguard and improve future diversity.

Genomic exploration of the journey of Plasmodium vivax in Latin America

Plasmodium vivax is the predominant malaria parasite in Latin America. Its colonization history in the region is rich and complex, and is still highly debated, especially about its origin(s). Our study employed cutting-edge population genomic techniques to analyze whole genome variation from 620 P. vivax isolates, including 107 newly sequenced samples from West Africa, Middle East, and Latin America. This sampling represents nearly all potential source populations worldwide currently available. Analyses of the genetic structure, diversity, ancestry, coalescent-based inferences, including demographic scenario testing using Approximate Bayesian Computation, have revealed a more complex evolutionary history than previously envisioned. Indeed, our analyses suggest that the current American P. vivax populations predominantly stemmed from a now-extinct European lineage, with the potential contribution also from unsampled populations, most likely of West African origin. We also found evidence that P. vivax arrived in Latin America in multiple waves, initially during early European contact and later through post-colonial human migration waves in the late 19th-century. This study provides a fresh perspective on P. vivax's intricate evolutionary journey and brings insights into the possible contribution of West African P. vivax populations to the colonization history of Latin America.

Phenotypic and Genomic Signatures of Latitudinal Local Adaptation Along With Prevailing Ocean Current in a Coastal Goby

In the marine realm, unidirectional ocean currents often lead to high migration rates of marine organisms and, therefore, inhibit the formation of their latitudinal genetic structure. In contrast, cryptic latitudinal structures associated with local adaptation may frequently exist in widespread species generally exposed to a strong environmental heterogeneity. However, our understanding of the evolvability of locally adapted populations in open marine environments still needs to be completed. The coastal area along the Sea of Japan, where the Tsushima Warm Current flows from south to north in the Japanese Archipelago, provides a good model system for exploring this question. This study explored evidence for latitudinal local adaptation along with the prevailing ocean current in the ice goby Leucopsarion petersii at the phenotypic and genomic levels. Common garden experiments clearly showed genetically based clinal variation in growth rate, strongly suggesting local adaptation through conutergradient selection of this fitness-related trait. Analyses based on reduced-representation sequencing revealed a slight signal of genetic differentiation between the southern and northern populations, although continuous historical gene flow between them was supported by demographic modelling. Also, whole-genome resequencing showed their independent demographic history during the last glacial period. Thus, these results suggest that gene flow along with the prevailing ocean current is somewhat limited, and the populations are not completely panmictic. Furthermore, the selection scan based on low-coverage genome-wide sequencing detected putative genomic signatures of latitudinal adaptation of growth-related genes. Thus, our integrative study provided a novel example of marine local adaptation under a large ocean current.

Whole Genome Sequencing Reveals Substantial Genetic Structure and Evidence of Local Adaptation in Alaskan Red King Crab

High-latitude ocean basins are the most productive on earth, supporting high diversity and biomass of economically and socially important species. A long tradition of responsible fisheries management has sustained these species for generations, but modern threats from climate change, habitat loss, and new fishing technologies threaten their ecosystems and the human communities that depend on them. Among these species, Alaska's most charismatic megafaunal invertebrate, the red king crab, faces all three of these threats and has declined substantially in many parts of its distribution. Managers have identified stock structure and local adaptation as crucial information to help understand biomass declines and how to potentially reverse them, with regulation and possible stock enhancement. We generated low-coverage whole genome sequencing (lcWGS) data on red king crabs from five regions: The Aleutian Islands, eastern Bering Sea, northern Bering Sea, Gulf of Alaska, and Southeast Alaska. We used data from millions of genetic markers generated from lcWGS to build on previous studies of population structure in Alaska that used < 100 markers and to investigate local adaptation. We found each of the regions formed their own distinct genetic clusters, some containing subpopulation structure. Most notably, we found that the Gulf of Alaska and eastern Bering Sea were significantly differentiated, something that had not been previously documented. Inbreeding in each region was low and not a concern for fisheries management. We found genetic patterns consistent with local adaptation on several chromosomes and one particularly strong signal on chromosome 100. At this locus, the Gulf of Alaska harbors distinct genetic variation that could facilitate local adaptation to their environment. Our findings support the current practice of managing red king crab at a regional scale, and they strongly favor sourcing broodstock from the target population if stock enhancement is considered to avoid genetic mismatch.

Widespread admixture blurs population structure and confounds Lake Trout (Salvelinus namaycush) conservation even in the genomic era

Intraspecific variation is important for species' long-term persistence in changing environments. Conservation strategies targeting intraspecific variation often rely on the identification of management or policy units below the species level based on biological differences among populations. To identify management units, this paper examines intraspecific divergence of Lake Trout (Salvelinus namaycush) in Great Slave Lake (GSL), Canada, using low-coverage whole-genome sequencing data. Specifically, we evaluate genetic differentiation and assess the relationship with morphological, mitochondrial, and putatively adaptive divergence. We show that at least three and up to five genetically distinct Lake Trout populations co-occur in GSL and exhibit differences in spatial distribution and body size, with signatures of selection. However, admixture was widespread (60% of the fish) and population structure was shallow (average FST = 0.022). These findings highlight that, even in the era of whole genome sequencing, identifying discrete management units to implement conservation efforts and policy can remain challenging in systems where gene flow among genetically and ecologically distinct populations is ubiquitous. To give more recognition to this complexity, management efforts should also focus on the area where adaptive variation is evident, and evolutionary acts are at play, to better protect species' resilience and adaptive potential in some natural systems.

Population genomics of the southern giraffe

Studying wildlife taxonomic diversity and identifying distinct populations has traditionally been largely based on morphology and geographic origin. More recently, this method has been supplemented by genetic data from the mitochondrial genome. However, this is limited as only maternally inherited and may not reflect the true nature of a population's genetics. Within the giraffe (Giraffa spp.), subspecies and unique populations were successfully characterized using both mitochondrial and genomic DNA studies, which led to new insights and, in some cases, unexpected results that required further verification. Here, we sequenced the genomes of 85 southern giraffe (G. giraffa) individuals from ten populations across southern Africa for a detailed investigation into the genetic diversity and history of its two subspecies, the Angolan (G. g. angolensis) and the South African (G. g. giraffa) giraffe. While the overall genotypes show low levels of runs of homozygosity compared to other mammals, the degree of heterozygosity is limited despite the large population size of South African giraffe. The nuclear genotype is largely congruent with the mitochondrial genotype. However, we have identified that the distribution of the Angolan giraffe is not as far east as indicated in an earlier mitochondrial DNA study. Botswana's Central Kalahari Game Reserve giraffe are unique, with a clear admixture of Angolan and South African giraffe populations. However, the enigmatic desert-dwelling giraffe of northwest Namibia is locally distinct from other Angolan giraffe yet exhibits intra-subspecies signs of admixture resulting from a recent introduction of individuals from Namibia's Etosha National Park. Whole genome sequencing is an invaluable and nearly indispensable tool for wildlife management to uncover genetic diversity that is undetectable through mitogenomic, geographical, and morphological means.

Substygophily in Dinaric Karst: A Model Case of Locally Endemic Minnows Phoxinellus (Leuciscinae)

The Dinaric Karst extends along the Adriatic coast of the Western Balkan Peninsula and is home to a group of "karst minnows" of the genera Delminichthys, Phoxinellus, and Telestes, which have adapted to the highly variable water conditions in the karst by spending up to several months underground, but require surface habitats for spawning, defining them as substygophiles. The three species of the genus Phoxinellus, P. alepidotus, P. pseudalepidotus, and P. dalmaticus, are defined by restricted ranges, making them vulnerable to pollution and extended draughts caused by the climate change. In this study, the phylogeny of Leusciscinae was reconstructed using 15 Phoxinellus and one Delminichthys adspersus, one Pelasgus epiroticus, and one Telestes polylepis complete mitochondrial genomes and the position of the genus Phoxinellus within the subfamily as sister species to the Chondrostoma clade was confirmed. The inter- and intrapopulation structure of the genus Phoxinellus was inferred using molecular (nuclear and mitochondrial data) and morphological analyses. For the molecular analysis, more than 150 historical specimens were analyzed for a short fragment of the cytochrome oxidase I (COI) barcoding region and 15 Phoxinellus specimens were subjected to single nucleotide polymorphism analysis. For morphological analysis, 121 Phoxinellus specimens were analyzed for 51 measurements and 8 counts. All analyses confirmed the clear delimitation of the three Phoxinellus species, but were insufficient to fully resolve the intrapopulation structure within the species. This study also included data from field surveys of Phoxinellus collected over the past 20 years, which showed that abundance is declining and ranges are shrinking. Phoxinellus are also threatened by invasive/introduced species. Based on cave observations/occurrence and morphological analysis, P. dalmaticus was classified as an advanced substygophile and P. alepidotus and P. pseudalepidotus were classified as basic stygophiles.

Whole-Genome Resequencing Reveals Polygenic Signatures of Directional and Balancing Selection on Alternative Migratory Life Histories

Migration in animals and associated adaptations to contrasting environments are underpinned by complex genetic architecture. Here, we explore the genomic basis of facultative anadromy in brown trout (Salmo trutta), wherein some individuals migrate to sea while others remain resident in natal rivers, to better understand how alternative migratory tactics (AMTs) are maintained evolutionarily. To identify genomic variants associated with AMTs, we sequenced whole genomes for 194 individual trout from five anadromous-resident population pairs, situated above and below waterfalls, in five different Irish rivers. These waterfalls act as natural barriers to upstream migration and hence we predicted that loci underpinning AMTs should be under similar divergent selection across these replicate pairs. A sliding windows based analysis revealed a highly polygenic adaptive divergence between anadromous and resident populations, encompassing 329 differentiated genomic regions. These regions were associated with 292 genes involved in various processes crucial for AMTs, including energy homeostasis, reproduction, osmoregulation, immunity, circadian rhythm and neural function. Furthermore, examining patterns of diversity we were able to link specific genes and biological processes to putative AMT trait classes: migratory-propensity, migratory-lifestyle and residency. Importantly, AMT outlier regions possessed higher genetic diversity than the background genome, particularly in the anadromous group, suggesting balancing selection may play a role in maintaining genetic variation. Overall, the results from this study provide important insights into the genetic architecture of migration and the evolutionary mechanisms shaping genomic diversity within and across populations.

Restricted Dispersal and Phenotypic Response to Water Depth in a Foundation Seagrass

Species conservation and management benefit from precise understanding of natural patterns of dispersal and genetic variation. Using recent advances in indirect genetic methods applied to both adult plants and dispersed seeds, we find that the mean seed dispersal in a threatened marine foundation plant (the eelgrass Zostera marina) is approximately 100-200 m. This distance is surprisingly more similar to that of wind-dispersed terrestrial seeds (~10s to 100s of meters) than the passive dispersal of marine propagules via currents (~10s to 100s of kilometres). Because nearshore marine plants like Zostera are commonly distributed across strong selective gradients driven by bathymetry (depth) even within these restricted spatial scales, seeds are capable of dispersing to novel water depths and experiencing profound shifts in light availability, temperature and wave exposure. We documented strong phenotypic variation and genome-wide differentiation among plants separated by approximately the spatial scale of mean realised dispersal. This result suggests genetic isolation by environment in response to depth-related environmental gradients as one plausible explanation for this pattern. The ratio of effective to census size (or Ne/Nc) approximated 0.1%, indicating that a fraction of existing plants provides the genetic variation to allow adaptation to environmental change. Our results suggest that successful conservation of seagrass meadows that can adapt to microspatial and temporal variation in environmental conditions will be low without direct and persistent intervention using large numbers of individuals or a targeted selection of genotypes.

Persistent Gene Flow Suggests an Absence of Reproductive Isolation in an African Antelope Speciation Model

African antelope diversity is a globally unique vestige of a much richer world-wide Pleistocene megafauna. Despite this, the evolutionary processes leading to the prolific radiation of African antelopes are not well understood. Here, we sequenced 145 whole genomes from both subspecies of the waterbuck (Kobus ellipsiprymnus), an African antelope believed to be in the process of speciation. We investigated genetic structure and population divergence and found evidence of a mid-Pleistocene separation on either side of the eastern Great Rift Valley, consistent with vicariance caused by a rain shadow along the so-called "Kingdon's Line." However, we also found pervasive evidence of both recent and widespread historical gene flow across the Rift Valley barrier. By inferring the genome-wide landscape of variation among subspecies, we found 14 genomic regions of elevated differentiation, including a locus that may be related to each subspecies' distinctive coat pigmentation pattern. We investigated these regions as candidate speciation islands. However, we observed no significant reduction in gene flow in these regions, nor any indications of selection against hybrids. Altogether, these results suggest a pattern whereby climatically driven vicariance is the most important process driving the African antelope radiation and suggest that reproductive isolation may not set in until very late in the divergence process. This has a significant impact on taxonomic inference, as many taxa will be in a gray area of ambiguous systematic status, possibly explaining why it has been hard to achieve consensus regarding the species status of many African antelopes. Our analyses demonstrate how population genetics based on low-depth whole genome sequencing can provide new insights that can help resolve how far lineages have gone along the path to speciation.

High Crimean-Congo hemorrhagic fever incidence linked to greater genetic diversity and differentiation in Hyalomma marginatum populations in Türkiye

BACKGROUND

Ticks are crucial vectors of a wide range of pathogens, posing significant threats to human and animal health globally. Understanding the genetic basis of tick biology and host-parasite interactions is essential for developing effective control programs. This study investigates the fine-scale genetic structure of Hyalomma marginatum Koch, 1844, the primary vector of Crimean-Congo hemorrhagic fever (CCHF) in Türkiye. Despite its significant public health importance, information regarding its population structure and genetic diversity is quite limited.

METHODS

We used restriction site-associated DNA sequencing (RAD-Seq) to obtain genome-wide sequence data from 10 tick populations in Türkiye, collected from regions with low, moderate, and high incidence rates of CCHF. Based on these data, we determined population structure and diversity of populations using principal component analysis (PCA) and admixture analysis. Furthermore, we calculated pairwise FST and utilized discriminant analysis of principal components (DAPC) to understand genetic differentiation between populations.

RESULTS

PCA and admixture analysis indicated minimal genetic structure between populations, but we detected notable genetic differentiation and high genetic diversity from regions with high CCHF rates. Furthermore, our DAPC identified 31 significant single-nucleotide polymorphisms (SNPs) associated with regions with high CCHF incidence, with 25 SNPs located near genes involved in critical biological functions such as nucleic acid binding, transmembrane transport, and proteolysis. These findings suggest that genetic variations in these regions may confer adaptive advantages in environments with high pathogen loads.

CONCLUSIONS

This study provides the first comprehensive analysis of H. marginatum genetic diversity in Türkiye, revealing significant differentiation in populations from CCHF-endemic regions. These results underscore the importance of considering fine-scale genetic diversity to fully understand the drivers of genetic variation in ticks and their implications for vectorial capacity.

Prioritizing Conservation Areas for the Hyacinth Macaw (Anodorhynchus hyacinthinus) in Brazil From Low-Coverage Genomic Data

Estimates of current genetic diversity and population connectivity are especially important for endangered species that are subject to illegal harvesting and trafficking. Genetic monitoring can also ensure that management units are sustaining viable populations, while estimating genetic structure and population dynamics can influence genetic rescue efforts and reintroduction from captive breeding and confiscated animals. The Hyacinth Macaw (Anodorhynchus hyacinthinus) is a charismatic endangered species with a fragmented (allopatric) distribution. Using low coverage genomes, we aimed to investigate the dynamics across the remaining three large disjunct populations of Hyacinth Macaws in Brazil to inform conservation strategies. We obtained low coverage DNA data for 54 individuals from seven sampling sites. Our results showed that Hyacinth Macaws have four genetically structured clusters with relatively high levels of diversity. The Pantanal biome had two genetically distinct populations, with no obvious physical barriers that might explain this differentiation. We detected signs of gene flow between populations, with some geographical regions being more connected than others. Estimates of effective population size in the past million years of the species' evolutionary history showed a decline trend with the lowest Ne in all populations reached within the last few thousand years. Our findings suggest that populations from the Pantanal biome are key to connecting sites across its distribution, and maintaining the integrity of this habitat is important for protecting the species. Given the genetic structure found, we also highlight the need of conserving all wild populations to ensure the protection of the species' evolutionary potential.

Genomic signatures of inbreeding and mutation load in tree ferns

Ferns (Pteridophyta), as the second largest group of vascular plants, play important roles in ecosystem functioning. Homosporous ferns exhibit a remarkable range of mating systems, from extreme inbreeding to obligate outcrossing, which may have significant evolutionary and ecological implications. Despite their significance, the impact of genome-wide inbreeding on genetic diversity and mutation load within the fern lineage remain largely unexplored. In this study, we utilized whole-genome sequencing to investigate the genomic signatures of inbreeding and genetic load in three Alsophila tree fern species. Our analysis revealed extremely high inbreeding in A. spinulosa, in contrast to the predominantly outcrossing observed in A. costularis and A. latebrosa. This difference likely reflects divergent mating systems and demographic histories. Consistent with its extreme inbreeding propensity, A. spinulosa exhibits reduced genetic diversity and a pronounced decline in effective population size. Comparison of genetic load revealed an overall reduction in deleterious mutations in the highly inbred A. spinulosa, highlighting that long-term inbreeding may have contributed to the purging of strongly deleterious mutations, thereby prolonging the survival of A. spinulosa. Despite this, however, A. spinulosa carries a substantive realized genetic load that may potentially instigate future fitness decline. Our findings illuminate the complex evolutionary interplay between inbreeding and mutation load in homosporous ferns, yielding insights with important implications for the conservation and management of these species.

Extensive Population Structure Highlights an Apparent Paradox of Stasis in the Impala (Aepyceros melampus)

Impalas are unusual among bovids because they have remained morphologically similar over millions of years-a phenomenon referred to as evolutionary stasis. Here, we sequenced 119 whole genomes from the two extant subspecies of impala, the common (Aepyceros melampus melampus) and black-faced (A. m. petersi) impala. We investigated the evolutionary forces working within the species to explore how they might be associated with its evolutionary stasis as a taxon. Despite being one of the most abundant bovid species, we found low genetic diversity overall, and a phylogeographic signal of spatial expansion from southern to eastern Africa. Contrary to expectations under a scenario of evolutionary stasis, we found pronounced genetic structure between and within the two subspecies with indications of ancient, but not recent, gene flow. Black-faced impala and eastern African common impala populations had more runs of homozygosity than common impala in southern Africa, and, using a proxy for genetic load, we found that natural selection is working less efficiently in these populations compared to the southern African populations. Together with the fossil record, our results are consistent with a fixed-optimum model of evolutionary stasis, in which impalas in the southern African core of the range are able to stay near their evolutionary fitness optimum as a generalist ecotone species, whereas eastern African impalas may struggle to do so due to the effects of genetic drift and reduced adaptation to the local habitat, leading to recurrent local extinction in eastern Africa and re-colonisation from the South.

Crossing the Pacific: Genomics Reveals the Presence of Japanese Sardine (Sardinops melanosticta) in the California Current Large Marine Ecosystem

Recent increases in frequency and intensity of warm water anomalies and marine heatwaves have led to shifts in species ranges and assemblages. Genomic tools can be instrumental in detecting such shifts. In the early stages of a project assessing population genetic structure in Pacific Sardine (Sardinops sagax), we detected the presence of Japanese Sardine (Sardinops melanosticta) along the west coast of North America for the first time. We assembled a high quality, chromosome-scale reference genome of the Pacific Sardine and generated low coverage, whole genome sequence (lcWGS) data for 345 sardine collected in the California Current Large Marine Ecosystem (CCLME) in 2021 and 2022. Fifty individuals sampled in 2022 were identified as Japanese Sardine based on strong differentiation observed in lcWGS SNP and full mitogenome data. Although we detected a single case of mitochondrial introgression, we did not observe evidence for recent hybridization events. These findings change our understanding of Sardinops spp. distribution and dispersal in the Pacific and highlight the importance of long-term monitoring programs.

Whole-genome resequencing improves the utility of otoliths as a critical source of DNA for fish stock research and monitoring

Fish ear bones, known as otoliths, are often collected in fisheries to assist in management, and are a common sample type in museum and national archives. Beyond their utility for ageing, morphological and trace element analysis, otoliths are a repository of valuable genomic information. Previous work has shown that DNA can be extracted from the trace quantities of tissue remaining on the surface of otoliths, despite the fact that they are often stored dry at room temperature. However, much of this work has used reduced representation sequencing methods in clean lab conditions, to achieve adequate yields of DNA, libraries and ultimately single-nucleotide polymorphisms (SNPs). Here, we pioneer the use of small-scale (spike-in) sequencing to screen contemporary otolith samples prepared in regular molecular biology (in contrast to clean) laboratories for contamination and quality levels, submitting for whole-genome resequencing only samples above a defined endogenous DNA threshold. Despite the typically low quality and quantity of DNA extracted from otoliths, we are able to produce whole-genome libraries and ultimately sets of filtered, unlinked and even putatively adaptive SNPs of ample numbers for downstream uses in population, climate and conservation genomics. By comparing with a set of tissue samples from the same species, we are able to highlight the quality and efficacy of otolith samples from DNA extraction and library preparation, to bioinformatic preprocessing and SNP calling. We provide detailed schematics, protocols and scripts of our approach, such that it can be adopted widely by the community, improving the use of otoliths as a source of valuable genomic data.

The genomic natural history of the aurochs

Now extinct, the aurochs (Bos primigenius) was a keystone species in prehistoric Eurasian and North African ecosystems, and the progenitor of cattle (Bos taurus), domesticates that have provided people with food and labour for millennia1. Here we analysed 38 ancient genomes and found 4 distinct population ancestries in the aurochs-European, Southwest Asian, North Asian and South Asian-each of which has dynamic trajectories that have responded to changes in climate and human influence. Similarly to Homo heidelbergensis, aurochsen first entered Europe around 650 thousand years ago2, but early populations left only trace ancestry, with both North Asian and European B. primigenius genomes coalescing during the most recent glaciation. North Asian and European populations then appear separated until mixing after the climate amelioration of the early Holocene. European aurochsen endured the more severe bottleneck during the Last Glacial Maximum, retreating to southern refugia before recolonizing from Iberia. Domestication involved the capture of a small number of individuals from the Southwest Asian aurochs population, followed by early and pervasive male-mediated admixture involving each ancestral strain of aurochs after domestic stocks dispersed beyond their cradle of origin.

Whole genomes from the extinct Xerces Blue butterfly can help identify declining insect species

The Xerces Blue (Glaucopsyche xerces) is considered to be the first butterfly to become extinct in historical times. It was notable for its chalky lavender wings with conspicuous white spots on the ventral wings. The last individuals were collected in their restricted habitat, in the dunes near the Presidio military base in San Francisco, in 1941. We sequenced the genomes of four 80- to 100-year-old Xerces Blue, and seven historical and one modern specimens of its closest relative, the Silvery Blue (Glaucopsyche lygdamus). We compared these to a novel annotated genome of the Green-Underside Blue (Glaucopsyche alexis). Phylogenetic relationships inferred from complete mitochondrial genomes indicate that Xerces Blue was a distinct species that diverged from the Silvery Blue lineage at least 850,000 years ago. Using nuclear genomes, both species experienced population growth during the Eemian interglacial period, but the Xerces Blue decreased to a very low effective population size subsequently, a trend opposite to that observed in the Silvery Blue. Runs of homozygosity and deleterious load in the former were significantly greater than in the later, suggesting a higher incidence of inbreeding. These signals of population decline observed in Xerces Blue could be used to identify and monitor other insects threatened by human activities, whose extinction patterns are still not well known.

The last days of Aporia crataegi (L.) in Britain: Evaluating genomic erosion in an extirpated butterfly

Current rates of habitat degradation and climate change are causing unprecedented declines in global biodiversity. Studies on vertebrates highlight how conservation genomics can be effective in identifying and managing threatened populations, but it is unclear how vertebrate-derived metrics of genomic erosion translate to invertebrates, with their markedly different population sizes and life histories. The Black-veined White butterfly (Aporia crataegi) was extirpated from Britain in the 1920s. Here, we sequenced historical DNA from 17 specimens collected between 1854 and 1924 to reconstruct demography and compare levels of genomic erosion between extirpated British and extant European mainland populations. We contrast these results using modern samples of the Common Blue butterfly (Polyommatus icarus); a species with relatively stable demographic trends in Great Britain. We provide evidence for bottlenecks in both these species around the period of post-glacial colonization of the British Isles. Our results reveal different demographic histories and Ne for both species, consistent with their fates in Britain, likely driven by differences in life history, ecology and genome size. Despite a difference, by an order of magnitude, in historical effective population sizes (Ne), reduction in genome-wide heterozygosity in A. crataegi was comparable to that in P. icarus. Symptomatic of A. crataegi's disappearance were marked increases in runs-of-homozygosity (RoH), potentially indicative of recent inbreeding, and accumulation of putatively mildly and weakly deleterious variants. Our results provide a rare glimpse of genomic erosion in a regionally extinct insect and support the potential use of genomic erosion metrics in identifying invertebrate populations or species in decline.

Species limits and hybridization in Andean leaf-eared mice (Phyllotis)

Leaf-eared mice (genus Phyllotis) are among the most widespread and abundant small mammals in the Andean Altiplano, but species boundaries and distributional limits are often poorly delineated due to sparse survey data from remote mountains and high-elevation deserts. Here we report a combined analysis of mitochondrial DNA variation and whole-genome sequence (WGS) variation in Phyllotis mice to delimit species boundaries, to assess the timescale of diversification of the group, and to examine evidence for interspecific hybridization. Estimates of divergence dates suggest that most diversification of Phyllotis occurred during the past 3 million years. Consistent with the Pleistocene Aridification hypothesis, our results suggest that diversification of Phyllotis largely coincided with climatically induced environmental changes in the mid- to late Pleistocene. Contrary to the Montane Uplift hypothesis, most diversification in the group occurred well after the major phase of uplift of the Central Andean Plateau. Species delimitation analyses revealed surprising patterns of cryptic diversity within several nominal forms, suggesting the presence of much undescribed alpha diversity in the genus. Results of genomic analyses revealed evidence of ongoing hybridization between the sister species Phyllotis limatus and P. vaccarum and suggest that the contemporary zone of range overlap between the two species represents an active hybrid zone.

Local adaptation to climate facilitates a global invasion

Local adaptation may facilitate range expansion during invasions, but the mechanisms promoting destructive invasions remain unclear. Cheatgrass (Bromus tectorum), native to Eurasia and Africa, has invaded globally, with particularly severe impacts in western North America. We sequenced 307 genotypes and conducted controlled experiments. We found that diverse lineages invaded North America, where long-distance gene flow is common. Ancestry and phenotypic clines in the native range predicted those in the invaded range, indicating pre-adapted genotypes colonized different regions. Common gardens showed directional selection on flowering time that reversed between warm and cold sites, potentially maintaining clines. In the Great Basin, genomic predictions of strong local adaptation identified sites where cheatgrass is most dominant. Preventing new introductions that may fuel adaptation is critical for managing ongoing invasions.

Late Pleistocene polar bear genomes reveal the timing of allele fixation in key genes associated with Arctic adaptation

The polar bear (Ursus maritimus) occupies a relatively narrow ecological niche, with many traits adapted for cold temperatures, movement across snow, ice and open water, and for consuming highly lipid-dense prey species. The divergence of polar bears from brown bears (Ursus arctos) and their adaptation to their Arctic lifestyle is a well-known example of rapid evolution. Previous research investigating whole genomes uncovered twelve key genes that are highly differentiated between polar and brown bears, show signatures of selection in the polar bear lineage, and are associated with polar bear adaptation to the Arctic environment. Further research suggested fixed derived alleles in these genes arose from selection on both standing variation and de novo mutations in the evolution of polar bears. Here, we reevaluate these findings based on a larger and geographically more representative dataset of 119 polar bears and 135 brown bears, and assess the timing of derived allele fixation in polar bears by incorporating the genomes of two Late Pleistocene individuals (aged 130-100,000 years old and 100-70,000 years old). In contrast with previous results, we found no evidence of derived alleles fixed in present-day polar bears within the key genes arising from de novo mutation. Most derived alleles fixed in present-day polar bears were also fixed in the Late Pleistocene polar bears, suggesting selection occurred prior to 70,000 years ago. However, some derived alleles fixed in present-day polar bears were not fixed in the two Late Pleistocene polar bears, including at sites within APOB, LYST, and TTN. These three genes are associated with cardiovascular function, metabolism, and pigmentation, suggesting selection may have acted on different loci at different times.

Recent speciation and adaptation to aridity in the ecologically diverse Pilbara region of Australia enabled the native tobaccos (Nicotiana; Solanaceae) to colonize all Australian deserts

Over the last 6 million years, the arid Australian Eremaean Zone (EZ) has remained as dry as it is today. A widely accepted hypothesis suggests that the flora and fauna of arid regions were more broadly distributed before aridification began. In Australia, this process started around 20 million years ago (Ma), leading to gradual speciation as the climate became increasingly arid. Here, we use genomic data to investigate the biogeography and timing of divergence of native allotetraploid tobaccos, Nicotiana section Suaveolentes (Solanaceae). The original allotetraploid migrants from South America were adapted to mesic areas of Australia and recently radiated in the EZ, including in sandy dune fields (only 1.2 Ma old), after developing drought adaptations. Coalescent and maximum likelihood analyses suggest that Nicotiana section Suaveolentes arrived on the continent around 6 Ma, with the ancestors of the Pilbara (Western Australian) lineages radiating there at the onset of extreme aridity 5 Ma by locally adapting to these various ancient, highly stable habitats. The Pilbara thus served as both a mesic refugium and cradle for adaptations to harsher conditions, due to its high topographical diversity, providing microhabitats with varying moisture levels and its proximity to the ocean, which buffers against extreme aridity. This enabled species like Nicotiana to survive in mesic refugia and subsequently adapt to more arid conditions. These results demonstrate that initially poorly adapted plant groups can develop novel adaptations in situ, permitting extensive and rapid dispersal despite the highly variable and unpredictable extreme conditions of the EZ.

Speciation in savanna birds in South America: The case of the Least Nighthawk Chordeiles pusillus (Aves: Caprimulgidae) in and out of the Amazon

The Least Nighthawk Chordeiles pusillus is widespread wherever there are savannas in the South American tropics, often in isolated patches, such as white-sands savannas in the Amazon rainforest realm. Here, we investigate genetic relationships between populations of the Least Nighthawk to understand historical processes leading to its diversification and to determine dispersal routes between northern and southern savannas by way of three hypothesized dispersal corridors by comparing samples from white-sand savannas to samples from other savannas outside of the Amazon rainforest region. We use 32 mtDNA samples from the range of C. pusillus to infer a dated phylogeny. In a subset of 17 samples, we use shotgun sequences to infer a distance-based phylogeny and to estimate individual admixture proportions. We calculate gene flow and shared alleles between white-sand and non-Amazonian populations using the ABBA-BABA test (D statistics), and Principal Component Analysis (PCA) to examine genetic structure within and between lineages. Finally, we use species distribution modelling (SDM) of conditions during the Last Glacial Maximum (LGM), currently, and in the future (2050-2080) to predict potential species occurrence under a climate change scenario. Two main clades (estimated to have diverged around 1.07 million years ago) were recovered with mtDNA sequences and Single Nucleotide Polymorphism (SNPs) and were supported by NGSadmix and PCA: one in the Amazon basin white-sand savannas, the other in the non-Amazonian savannas. Possible allele sharing between these clades was indicated by the D-statistics between northern non-Amazonian populations and the white-sand savanna population, but this was not corroborated by the admixture analyses. Dispersal among northern non-Amazonian populations may have occurred in a dry corridor between the Guianan and the Brazilian Shield, which has since moved eastward. Our data suggest that the lineages separated well before the Last Glacial Maximum, consequently dispersal could have happened at any earlier time during similar climatic conditions. Subsequently, non-Amazonian lineages became more divergent among themselves, possibly connecting and dispersing across the mouth of the Amazon River across Marajó island during favourable climatic conditions in the Pleistocene.

Contemporary intergeneric hybridization and backcrossing among birds-of-paradise

Despite large differences in morphology, behavior and lek-mating strategies the birds-of-paradise are known to hybridize occasionally, even across different genera. Many of these bird-of-paradise hybrids were originally described as distinct species based on large morphological differences when compared to recognized species. Nowadays, these specimens are generally recognized as hybrids based on morphological assessments. Having fascinated naturalists for centuries, hybrid specimens of birds-of-paradise have been collected and the specimens kept in Natural History Collections. In the present study, we utilize this remarkable resource in a museomics framework and evaluate the genomic composition of most described intergeneric hybrids and some intrageneric hybrids. We show that the majority of investigated specimens are first-generation hybrids and that the parental species, in most cases, are in line with prior morphological assessments. We also identify two specimens that are the result of introgressive hybridization between different genera. Additionally, two specimens exhibit hybrid morphologies but have no identifiable signals of hybridization, which may indicate that minor levels of introgression can have large morphological effects. Our findings provide direct evidence of contemporary introgressive hybridization taking place between genera of birds-of-paradise in nature, despite markedly different morphologies and lek-mating behaviors.

Cryptic diversity of shallow and mesophotic Stephanocoenia intersepta corals across Florida Keys National Marine Sanctuary

Population genetic analyses can provide useful data on species' regional connectivity and diversity which can inform conservation and restoration efforts. In this study, we quantified the genetic connectivity and diversity of Stephanocoenia intersepta corals from shallow (<30 m) to mesophotic (30-45 m) depths across Florida Keys National Marine Sanctuary. We generated single nucleotide polymorphism (SNP) markers to identify genetic structuring of shallow and mesophotic S. intersepta corals. We uncovered four distinct, cryptic genetic lineages with varying levels of depth-specificity. Shallow-specific lineages exhibited lower heterozygosity and higher inbreeding relative to depth-generalist lineages found across both shallow and mesophotic reefs. Estimation of recent genetic migration rates demonstrated that mesophotic sites are more prolific sources than shallow sites, particularly in the Lower Keys and Upper Keys. Additionally, we compared endosymbiotic Symbiodiniaceae among sampled S. intersepta using the ITS2 region and SYMPORTAL analysis framework, identifying symbionts from the genera Symbiodinium, Breviolum, and Cladocopium. Symbiodiniaceae varied significantly across depth and location and exhibited significant, but weak correlation with host lineage and genotype. Together, these data demonstrate that despite population genetic structuring across depth, some mesophotic populations may provide refuge for shallow populations moving forward and remain important contributors to the overall genetic diversity of this species throughout the region. This study highlights the importance of including mesophotic as well as shallow corals in population genetic assessments and informs future science-based management, conservation, and restoration efforts within Florida Keys National Marine Sanctuary.

Forest genomics in the Caucasus through the lens of its dominant tree species - Fagus orientalis

The last glacial period is known to have greatly influenced the demographic history of temperate forest trees, with important range contractions and post-glacial expansions that led to the formation of multiple genetic lineages and secondary contact zones in the Northern Hemisphere. These dynamics have been extensively studied for European and North American species but are still poorly understood in other temperate regions of rich biodiversity such as the Caucasus. Our study helps filling that gap by deciphering the genomic landscapes of F. orientalis across the South Caucasus. The use of genome-wide data confirmed a past demographic history strongly influenced by the Last Glacial Maximum, revealing two disjunct glacial refugia in the Colchis and Hyrcanian regions. The resulting patterns of genetic diversity, load and differentiation are not always concordant across the region, with genetic load pinpointing the location of the glacial refugia more efficiently than genetic diversity alone. The Hyrcanian forests show depleted genetic diversity and substantial isolation, even if long-distance gene flow is still present with the main centre of diversity in the Greater Caucasus. Finally, we characterize a strong heterogeneity of genetic diversity and differentiation along the species chromosomes, with noticeably a first chromosome showing low diversity and weak differentiation.

Scots pine - panmixia and the elusive signal of genetic adaptation

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Discordance between taxonomy and population genomic data: An avian example relevant to the United States Endangered Species Act

Population genomics can reveal cryptic biological diversity that may impact fitness while simultaneously serving to delineate relevant conservation units. Here, we leverage the power of whole-genome resequencing for conservation by studying 433 individual lesser prairie-chicken (Tympanuchus pallidicinctus; LEPC, a federally endangered species of conservation concern in the United States) and greater prairie-chicken (Tympanuchus cupido; GRPC, a legally huntable species throughout much of its range). The genomic diversity of two formally recognized distinct population segments (DPSs) of LEPCs is similar, but they are genetically distinct. Neither DPS is depleted of its genomic diversity, neither is especially inbred, and temporal diversity is relatively stable in both conservation units. Interspecific differentiation between the two species was only slightly higher than that observed between LEPC DPSs, due largely to bidirectional introgression. The high resolution provided by our dataset identified a genomic continuum between the two species such that individuals sampled from the hybrid zone were imperfectly assigned to their presumptive species when considering only their physical characteristics. The admixture between the two species is reflected in the spectrum of individual ancestry coefficients, which has legal implications for the "take" of individuals under the Endangered Species Act. Overall, our data highlight the recurring dissonance between static policies and dynamic species boundaries that are increasingly obvious in the population genomic era.

Population structure in Mixornis tit-babblers across Sunda Shelf matches interfluvia of paleo-rivers

Rivers constitute an important biogeographic divide in vast areas of tropical rainforest, such as the Amazon and Congo Basins. Southeast Asia's rainforests are currently fragmented across islands divided by sea, which has long obscured their extensive history of terrestrial connectivity as part of a vast (but now submerged) subcontinent - Sundaland - during most of the Quaternary. The role of paleo-rivers in determining population structure in Sundaic rainforests at a time when these forests were connected remains little understood. We examined the coloration of museum skins and used the genomic DNA of museum samples and freshly-collected blood tissue of a pair of Sundaic songbird species, the pin-striped and bold-striped tit-babblers (Mixornis gularis and M. bornensis, respectively), to assess the genetic affinity of populations on small Sundaic islands that have largely been ignored by modern research. Our genomic and morphological results place the populations from the Anambas and Natuna Islands firmly within M. gularis from the Malay Peninsula in western Sundaland, even though some of these islands are geographically much closer to Borneo, where M. bornensis resides. Our results reveal genetic structure consistent with the course of Sundaic paleo-rivers and the location of the interfluvia they formed, and add to a small but growing body of evidence that rivers would have been of equal biogeographic importance in Sundaland's former connected forest landscape as they are in Amazonia and the Congo Basin today.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

Site-directed genotype screening for elimination of antinutritional saponins in quinoa seeds identifies TSARL1 as a master controller of saponin biosynthesis selectively in seeds

Climate change may result in a drier climate and increased salinization, threatening agricultural productivity worldwide. Quinoa (Chenopodium quinoa) produces highly nutritious seeds and tolerates abiotic stresses such as drought and high salinity, making it a promising future food source. However, the presence of antinutritional saponins in their seeds is an undesirable trait. We mapped genes controlling seed saponin content to a genomic region that includes TSARL1. We isolated desired genetic variation in this gene by producing a large mutant library of a commercial quinoa cultivar and screening the library for specific nucleotide substitutions using droplet digital PCR. We were able to rapidly isolate two independent tsarl1 mutants, which retained saponins in the leaves and roots for defence, but saponins were undetectable in the seed coat. We further could show that TSARL1 specifically controls seed saponin biosynthesis in the committed step after 2,3-oxidosqualene. Our work provides new important knowledge on the function of TSARL1 and represents a breakthrough for quinoa breeding.

Selection against domestication alleles in introduced rabbit populations

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

Modeling biases from low-pass genome sequencing to enable accurate population genetic inferences

Low-pass genome sequencing is cost-effective and enables analysis of large cohorts. However, it introduces biases by reducing heterozygous genotypes and low-frequency alleles, impacting subsequent analyses such as demographic history inference. We developed a probabilistic model of low-pass biases from the Genome Analysis Toolkit (GATK) multi-sample calling pipeline, and we implemented it in the population genomic inference software dadi. We evaluated the model using simulated low-pass datasets and found that it alleviated low-pass biases in inferred demographic parameters. We further validated the model by downsampling 1000 Genomes Project data, demonstrating its effectiveness on real data. Our model is widely applicable and substantially improves model-based inferences from low-pass population genomic data.

Hybridization in birds-of-paradise: Widespread ancestral gene flow despite strong sexual selection in a lek-mating system

Sexual selection can directly contribute to reproductive isolation and is an important mechanism that can lead to speciation. Lek-mating is one of the most extreme forms of sexual selection, but surprisingly does not seem to preclude occasional hybridization in nature. However, hybridization among lekking species may still be trivial if selection against offspring with intermediate phenotypes prohibits introgression. Here we investigate this further by sequencing the genomes of nearly all bird-of-paradise (Paradisaeidae) species and 10 museum specimens of putative hybrid origin. We find that intergeneric hybridization indeed still takes place despite extreme differentiation in form, plumage, and behavior. In parallel, the genomes of contemporary species contain widespread signatures of past introgression, demonstrating that hybridization has repeatedly resulted in shared genetic variation despite strong sexual isolation. Our study raises important questions about extrinsic factors that modulate hybridization probability and the evolutionary consequences of introgressive hybridization between lekking species.

loco-pipe: an automated pipeline for population genomics with low-coverage whole-genome sequencing

Summary

We developed loco-pipe, a Snakemake pipeline that seamlessly streamlines a set of essential population genomic analyses for low-coverage whole genome sequencing (lcWGS) data. loco-pipe is highly automated, easily customizable, massively parallelized, and thus is a valuable tool for both new and experienced users of lcWGS.

Availability and implementation

loco-pipe is published under the GPLv3. It is freely available on GitHub (github.com/sudmantlab/loco-pipe) and archived on Zenodo (doi.org/10.5281/zenodo.10425920).

Genomic analysis of wolves from Pakistan clarifies boundaries among three divergent wolf lineages

Among the three main divergent lineages of gray wolf (Canis lupus), the Holarctic lineage is the most widespread and best studied, particularly in North America and Europe. Less is known about Tibetan (also called Himalayan) and Indian wolf lineages in southern Asia, especially in areas surrounding Pakistan where all three lineages are thought to meet. Given the endangered status of the Indian wolf in neighboring India and unclear southwestern boundary of the Tibetan wolf range, we conducted mitochondrial and genome-wide sequencing of wolves from Pakistan and Kyrgyzstan. Sequences of the mitochondrial D-loop region of 81 wolves from Pakistan indicated contact zones between Holarctic and Indian lineages across the northern and western mountains of Pakistan. Reduced-representation genome sequencing of eight wolves indicated an east-to-west cline of Indian to Holarctic ancestry, consistent with a contact zone between these two lineages in Pakistan. The western boundary of the Tibetan lineage corresponded to the Ladakh region of India's Himalayas with a narrow zone of admixture spanning this boundary from the Karakoram Mountains of northern Pakistan into Ladakh, India. Our results highlight the conservation significance of Pakistan's wolf populations, especially the remaining populations in Sindh and Southern Punjab that represent the highly endangered Indian lineage.

Local genetic adaptation to habitat in wild chimpanzees

How populations adapt to their environment is a fundamental question in biology. Yet we know surprisingly little about this process, especially for endangered species such as non-human great apes. Chimpanzees, our closest living relatives, are particularly interesting because they inhabit diverse habitats, from rainforest to woodland-savannah. Whether genetic adaptation facilitates such habitat diversity remains unknown, despite having wide implications for evolutionary biology and conservation. Using 828 newly generated exomes from wild chimpanzees, we find evidence of fine-scale genetic adaptation to habitat. Notably, adaptation to malaria in forest chimpanzees is mediated by the same genes underlying adaptation to malaria in humans. This work demonstrates the power of non-invasive samples to reveal genetic adaptations in endangered populations and highlights the importance of adaptive genetic diversity for chimpanzees.

Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy

South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.

Ancient and Recent Hybridization in the Oreochromis Cichlid Fishes

Cichlid fishes of the genus Oreochromis (tilapia) are among the most important fish for inland capture fisheries and global aquaculture. Deliberate introductions of non-native species for fisheries improvement and accidental escapees from farms have resulted in admixture with indigenous species. Such hybridization may be detrimental to native biodiversity, potentially leading to genomic homogenization of populations and the loss of important genetic material associated with local adaptation. By contrast, introgression may fuel diversification when combined with ecological opportunity, by supplying novel genetic combinations. To date, the role of introgression in the evolutionary history of tilapia has not been explored. Here we studied both ancient and recent hybridization in tilapia, using whole genome resequencing of 575 individuals from 23 species. We focused on Tanzania, a natural hotspot of tilapia diversity, and a country where hybridization between exotic and native species in the natural environment has been previously reported. We reconstruct the first genome-scale phylogeny of the genus and reveal prevalent ancient gene flow across the Oreochromis phylogeny. This has likely resulted in the hybrid speciation of one species, O. chungruruensis. We identify multiple cases of recent hybridization between native and introduced species in the wild, linked to the use of non-native species in both capture fisheries improvement and aquaculture. This has potential implications for both conservation of wild populations and the development of the global tilapia aquaculture industry.