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

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.

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.

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.

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.

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.

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.

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.

Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

When birds of a feather flock together: Severe genomic erosion and the implications for genetic rescue in an endangered island passerine

The Seychelles magpie-robin's (SMR) five island populations exhibit some of the lowest recorded levels of genetic diversity among endangered birds, and high levels of inbreeding. These populations collapsed during the 20th century, and the species was listed as Critically Endangered in the IUCN Red List in 1994. An assisted translocation-for-recovery program initiated in the 1990s increased the number of mature individuals, resulting in its downlisting to Endangered in 2005. Here, we explore the temporal genomic erosion of the SMR based on a dataset of 201 re-sequenced whole genomes that span the past ~150 years. Our sample set includes individuals that predate the bottleneck by up to 100 years, as well as individuals from contemporary populations established during the species recovery program. Despite the SMR's recent demographic recovery, our data reveal a marked increase in both the genetic load and realized load in the extant populations when compared to the historical samples. Conservation management may have reduced the intensity of selection by increasing juvenile survival and relaxing intraspecific competition between individuals, resulting in the accumulation of loss-of-function mutations (i.e. severely deleterious variants) in the rapidly recovering population. In addition, we found a 3-fold decrease in genetic diversity between temporal samples. While the low genetic diversity in modern populations may limit the species' adaptability to future environmental changes, future conservation efforts (including IUCN assessments) may also need to assess the threats posed by their high genetic load. Our computer simulations highlight the value of translocations for genetic rescue and show how this could halt genomic erosion in threatened species such as the SMR.

Cryptic hybridization between the ancient lineages of Natterer's bat (Myotis nattereri)

Studying hybrid zones that form between morphologically cryptic taxa offers valuable insights into the mechanisms of cryptic speciation and the evolution of reproductive barriers. Although hybrid zones have long been the focus of evolutionary studies, the awareness of cryptic hybrid zones increased recently due to rapidly growing evidence of biological diversity lacking obvious phenotypic differentiation. The characterization of cryptic hybrid zones with genome-wide analysis is in its early stages and offers new perspectives for studying population admixture and thus the impact of gene flow. In this study, we investigate the population genomics of the Myotis nattereri complex in one of its secondary contact zones, where a putative hybrid zone is formed between two of its cryptic lineages. By utilizing a whole-genome shotgun sequencing approach, we aim to characterize this cryptic hybrid zone in detail. Demographic analysis suggests that the cryptic lineages diverged during the Pliocene, c. 3.6 million years ago. Despite this ancient separation, the populations in the contact zone exhibit mitochondrial introgression and a considerable amount of mixing in nuclear genomes. The genomic structure of the populations corresponds to geographic locations and the genomic admixture changes along a geographic gradient. These findings suggest that there is no effective hybridization barrier between both lineages, nevertheless, their population structure is shaped by dispersal barriers. Our findings highlight how such deeply diverged cryptic lineages can still readily hybridize in secondary contact.

Genetic diversity and population structure of Piper nigrum (black pepper) accessions based on next-generation SNP markers

Despite the economic importance of Piper nigrum (black pepper), a highly valued crop worldwide, development and utilization of genomic resources have remained limited, with diversity assessments often relying on only a few samples or DNA markers. Here we employed restriction-site associated DNA sequencing to analyze 175 P. nigrum accessions from eight main black pepper growing regions in Sri Lanka. The sequencing effort resulted in 1,976 million raw reads, averaging 11.3 million reads per accession, revealing 150,356 high-quality single nucleotide polymorphisms (SNPs) distributed across 26 chromosomes. Population structure analysis revealed two subpopulations (K = 2): a dominant group consisting of 152 accessions sourced from both home gardens and large-scale cultivations, and a smaller group comprising 23 accessions exclusively from native collections in home gardens. This clustering was further supported by principal component analysis, with the first two principal components explaining 35.2 and 12.1% of the total variation. Genetic diversity analysis indicated substantial gene flow (Nm = 342.21) and a low fixation index (FST = 0.00073) between the two subpopulations, with no clear genetic differentiation among accessions from different agro-climatic regions. These findings demonstrate that most current black pepper genotypes grown in Sri Lanka share a common genetic background, emphasizing the necessity to broaden the genetic base to enhance resilience to biotic and abiotic stresses. This study represents the first attempt at analyzing black pepper genetic diversity using high-resolution SNP markers, laying the foundation for future genome-wide association studies for SNP-based gene discovery and breeding.

Camelus knoblochi genome reveals the complex evolutionary history of Old World camels

Extant Old World camels (genus Camelus) contributed to the economic and cultural exchanges between the East and West for thousands of years.1,2 Although many remains have been unearthed,3,4,5 we know neither whether the prevalent hybridization observed between extant Camelus species2,6,7 also occurred between extinct lineages and the ancestors of extant Camelus species nor why some populations became extinct while others survived. To investigate these questions, we generated paleogenomic and stable isotope data from an extinct two-humped camel species, Camelus knoblochi. We find that in the mitochondrial phylogeny, all C. knoblochi form a paraphyletic group that nests within the diversity of modern, wild two-humped camels (Camelus ferus). In contrast, they are clearly distinguished from both wild and domesticated (Camelus bactrianus) two-humped camels on the nuclear level. Moreover, the divergence pattern of the three camel species approximates a trifurcation, because the most common topology is only slightly more frequent than the two other possible topologies. This mito-nuclear phylogenetic discordance likely arose due to interspecific gene flow between all three species, suggesting that interspecific hybridization is not exclusive to modern camels but a recurrent phenomenon throughout the evolutionary history of the genus Camelus. These results suggest that the genomic complexity of Old World camels' evolutionary history is underestimated when considering data from only modern species. Finally, we find that C. knoblochi populations began declining prior to the last glacial maximum and, by integrating palaeoecological evidence and stable isotope data, suggest that this was likely due to failure to adapt to a changing environment.

Co-phylogeographic structure in a disease-causing parasite and its oyster host

With the increasing affordability of next-generation sequencing technologies, genotype-by-sequencing has become a cost-effective tool for ecologists and conservation biologists to describe a species' evolutionary history. For host–parasite interactions, genotype-by-sequencing can allow the simultaneous examination of host and parasite genomes and can yield insight into co-evolutionary processes. The eastern oyster, Crassostrea virginica, is among the most important aquacultured species in the United States. Natural and farmed oyster populations can be heavily impacted by ‘dermo’ disease caused by an alveolate protist, Perkinsus marinus. Here, we used restricted site-associated DNA sequencing (RADseq) to simultaneously examine spatial population genetic structure of host and parasite. We analysed 393 single-nucleotide polymorphisms (SNPs) for P. marinus and 52,100 SNPs for C. virginica from 36 individual oysters from the Gulf of Mexico (GOM) and mid-Atlantic coastline. All analyses revealed statistically significant genetic differentiation between the GOM and mid-Atlantic coast populations for both C. virginica and P. marinus, and genetic divergence between Chesapeake Bay and the outer coast of Virginia for C. virginica, but not for P. marinus. A co-phylogenetic analysis confirmed significant coupled evolutionary change between host and parasite across large spatial scales. The strong genetic divergence between marine basins raises the possibility that oysters from either basin would not be well adapted to parasite genotypes and phenotypes from the other, which would argue for caution with regard to both oyster and parasite transfers between the Atlantic and GOM regions. More broadly, our results demonstrate the potential of RADseq to describe spatial patterns of genetic divergence consistent with coupled evolution.

Evolution of Chromosomal Inversions across an Avian Radiation

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

Extreme High-Elevation Mammal Surveys Reveal Unexpectedly High Upper Range Limits of Andean Mice

AbstractIn the world's highest mountain ranges, uncertainty about the upper elevational range limits of alpine animals represents a critical knowledge gap regarding the environmental limits of life and presents a problem for detecting range shifts in response to climate change. Here we report results of mountaineering mammal surveys in the Central Andes, which led to the discovery of multiple species of mice living at extreme elevations that far surpass previously assumed range limits for mammals. We livetrapped small mammals from ecologically diverse sites spanning >6,700 m of vertical relief, from the desert coast of northern Chile to the summits of the highest volcanoes in the Andes. We used molecular sequence data and whole-genome sequence data to confirm the identities of species that represent new elevational records and to test hypotheses regarding species limits. These discoveries contribute to a new appreciation of the environmental limits of vertebrate life.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

Rapid vertebrate speciation via isolation, bottlenecks, and drift

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Species-specific dynamics may cause deviations from general biogeographical predictions - evidence from a population genomics study of a New Guinean endemic passerine bird family (Melampittidae)

The family Melampittidae is endemic to New Guinea and consists of two monotypic genera: Melampitta lugubris (Lesser Melampitta) and Megalampitta gigantea (Greater Melampitta). Both Melampitta species have scattered and disconnected distributions across New Guinea in the central mountain range and in some of the outlying ranges. While M. lugubris is common and found in most montane regions of the island, M. gigantaea is elusive and known from only six localities in isolated pockets on New Guinea with very specific habitats of limestone and sinkholes. In this project, we apply museomics to determine the population structure and demographic history of these two species. We re-sequenced the genomes of all seven known M. gigantaea samples housed in museum collections as well as 24 M. lugubris samples from across its distribution. By comparing population structure between the two species, we investigate to what extent habitat dependence, such as in M. gigantaea, may affect population connectivity. Phylogenetic and population genomic analyses, as well as acoustic variation revealed that M. gigantaea consists of a single population in contrast to M. lugubris that shows much stronger population structure across the island. We suggest a recent collapse of M. gigantaea into its fragmented habitats as an explanation to its unexpected low diversity and lack of population structure. The deep genetic divergences between the M. lugubris populations on the Vogelkop region, in the western central range and the eastern central range, respectively, suggests that these three populations should be elevated to full species level. This work sheds new light on the mechanisms that have shaped the intriguing distribution of the two species within this family and is a prime example of the importance of museum collections for genomic studies of poorly known and rare species.

Patterns of genetic divergence in the Rio Grande cooter (Pseudemys gorzugi), a riverine turtle inhabiting an arid and anthropogenically modified system

The lower Rio Grande and Pecos River of the southwest United States have been heavily modified by human activities, profoundly impacting the integrity of their aquatic wildlife. In this context, we focused our study on the population genomics of the Rio Grande Cooter (Pseudemys gorzugi), a freshwater turtle of increasing conservation concern, residing in these two rivers and their tributaries. The genetic data revealed two distinct populations: one in the Pecos and Black Rivers of New Mexico and another in the Rio Grande and Devils River of Texas, with admixed individuals identified at the confluence of the Rio Grande and Pecos River. In addition to having a smaller geographic range, we found lower observed heterozygosity, reduced nucleotide diversity, and a smaller effective population size (Ne) in New Mexico population. Our results depict a significant isolation-by-distance pattern across their distribution, with migration being notably infrequent at river confluences. These findings are pivotal for future conservation and restoration strategies, emphasizing the need to recognize the unique needs of each population.

Unlocking the secret life of blue mussels: Exploring connectivity in the Skagerrak through biophysical modeling and population genomics

Knowledge of functional dispersal barriers in the marine environment can be used to inform a wide variety of management actions, such as marine spatial planning, restoration efforts, fisheries regulations, and invasive species management. Locations and causes of dispersal barriers can be studied through various methods, including movement tracking, biophysical modeling, demographic models, and genetics. Combining methods illustrating potential dispersal, such as biophysical modeling, with realized dispersal through, e.g., genetic connectivity estimates, provides particularly useful information for teasing apart potential causes of observed barriers. In this study, we focus on blue mussels (Mytilus edulis) in the Skagerrak-a marginal sea connected to the North Sea in Northern Europe-and combine biophysical models of larval dispersal with genomic data to infer locations and causes of dispersal barriers in the area. Results from both methods agree; patterns of ocean currents are a major structuring factor in the area. We find a complex pattern of source-sink dynamics with several dispersal barriers and show that some areas can be isolated despite an overall high dispersal capability. Finally, we translate our finding into management advice that can be used to sustainably manage this ecologically and economically important species in the future.

Evolution and genetic architecture of sex-limited polymorphism in cuckoos

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

Parallel evolution despite low genetic diversity in three-spined sticklebacks

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south.

Giraffe lineages are shaped by major ancient admixture events

Strong genetic structure has prompted discussion regarding giraffe taxonomy,1,2,3 including a suggestion to split the giraffe into four species: Northern (Giraffa c. camelopardalis), Reticulated (G. c. reticulata), Masai (G. c. tippelskirchi), and Southern giraffes (G. c. giraffa).4,5,6 However, their evolutionary history is not yet fully resolved, as previous studies used a simple bifurcating model and did not explore the presence or extent of gene flow between lineages. We therefore inferred a model that incorporates various evolutionary processes to assess the drivers of contemporary giraffe diversity. We analyzed whole-genome sequencing data from 90 wild giraffes from 29 localities across their current distribution. The most basal divergence was dated to 280 kya. Genetic differentiation, FST, among major lineages ranged between 0.28 and 0.62, and we found significant levels of ancient gene flow between them. In particular, several analyses suggested that the Reticulated lineage evolved through admixture, with almost equal contribution from the Northern lineage and an ancestral lineage related to Masai and Southern giraffes. These new results highlight a scenario of strong differentiation despite gene flow, providing further context for the interpretation of giraffe diversity and the process of speciation in general. They also illustrate that conservation measures need to target various lineages and sublineages and that separate management strategies are needed to conserve giraffe diversity effectively. Given local extinctions and recent dramatic declines in many giraffe populations, this improved understanding of giraffe evolutionary history is relevant for conservation interventions, including reintroductions and reinforcements of existing populations.

Resilience to periodic disturbances and the long-term genetic stability in Acropora coral

Climate change is restructuring natural ecosystems. The direct impacts of these events on biodiversity and community structure are widely documented, but the impacts on the genetic variation of populations remains largely unknown. We monitored populations of Acropora coral on a remote coral reef system in northwest Australia for two decades and through multiple cycles of impact and recovery. We combined these demographic data with a temporal genetic dataset of a common broadcast spawning corymbose Acropora to explore the spatial and temporal patterns of connectivity underlying recovery. Our data show that broad-scale dispersal and post-recruitment survival drive recovery from recurrent disturbances, including mass bleaching and mortality. Consequently, genetic diversity and associated patterns of connectivity are maintained through time in the broader metapopulation. The results highlight an inherent resilience in these globally threatened species of coral and showcase their ability to cope with multiple disturbances, given enough time to recover is permitted.

Calling Structural Variants with Confidence from Short-Read Data in Wild Bird Populations

Comprehensive characterization of structural variation in natural populations has only become feasible in the last decade. To investigate the population genomic nature of structural variation, reproducible and high-confidence structural variation callsets are first required. We created a population-scale reference of the genome-wide landscape of structural variation across 33 Nordic house sparrows (Passer domesticus). To produce a consensus callset across all samples using short-read data, we compare heuristic-based quality filtering and visual curation (Samplot/PlotCritic and Samplot-ML) approaches. We demonstrate that curation of structural variants is important for reducing putative false positives and that the time invested in this step outweighs the potential costs of analyzing short-read-discovered structural variation data sets that include many potential false positives. We find that even a lenient manual curation strategy (e.g. applied by a single curator) can reduce the proportion of putative false positives by up to 80%, thus enriching the proportion of high-confidence variants. Crucially, in applying a lenient manual curation strategy with a single curator, nearly all (>99%) variants rejected as putative false positives were also classified as such by a more stringent curation strategy using three additional curators. Furthermore, variants rejected by manual curation failed to reflect the expected population structure from SNPs, whereas variants passing curation did. Combining heuristic-based quality filtering with rapid manual curation of structural variants in short-read data can therefore become a time- and cost-effective first step for functional and population genomic studies requiring high-confidence structural variation callsets.

Variable parallelism in the genomic basis of age at maturity across spatial scales in Atlantic Salmon

Complex traits often exhibit complex underlying genetic architectures resulting from a combination of evolution from standing variation, hard and soft sweeps, and alleles of varying effect size. Increasingly, studies implicate both large-effect loci and polygenic patterns underpinning adaptation, but the extent that common genetic architectures are utilized during repeated adaptation is not well understood. Sea age or age at maturation represents a significant life history trait in Atlantic Salmon (Salmo salar), the genetic basis of which has been studied extensively in European Atlantic populations, with repeated identification of large-effect loci. However, the genetic basis of sea age within North American Atlantic Salmon populations remains unclear, as does the potential for a parallel trans-Atlantic genomic basis to sea age. Here, we used a large single-nucleotide polymorphism (SNP) array and low-coverage whole-genome resequencing to explore the genomic basis of sea age variation in North American Atlantic Salmon. We found significant associations at the gene and SNP level with a large-effect locus (vgll3) previously identified in European populations, indicating genetic parallelism, but found that this pattern varied based on both sex and geographic region. We also identified nonrepeated sets of highly predictive loci associated with sea age among populations and sexes within North America, indicating polygenicity and low rates of genomic parallelism. Despite low genome-wide parallelism, we uncovered a set of conserved molecular pathways associated with sea age that were consistently enriched among comparisons, including calcium signaling, MapK signaling, focal adhesion, and phosphatidylinositol signaling. Together, our results indicate parallelism of the molecular basis of sea age in North American Atlantic Salmon across large-effect genes and molecular pathways despite population-specific patterns of polygenicity. These findings reveal roles for both contingency and repeated adaptation at the molecular level in the evolution of life history variation.