Identifying outlier loci in admixed and in continuous populations using ancestral population differentiation statistics

Genomic insights into the critically endangered King Island scrubtit

Small, fragmented, or isolated populations are at risk of population decline due to fitness costs associated with inbreeding and genetic drift. The King Island scrubtit Acanthornis magna greeniana is a critically endangered subspecies of the nominate Tasmanian scrubtit A. m. magna, with an estimated population of <100 individuals persisting in three patches of swamp forest. The Tasmanian scrubtit is widespread in wet forests on mainland Tasmania. We sequenced the scrubtit genome using PacBio HiFi and undertook a population genomic study of the King Island and Tasmanian scrubtits using a double-digest restriction site-associated DNA (ddRAD) dataset of 5,239 SNP loci. The genome was 1.48 Gb long, comprising 1,518 contigs with an N50 of 7.715 Mb. King Island scrubtits formed one of four overall genetic clusters, but separated into three distinct subpopulations when analyzed independently of the Tasmanian scrubtit. Pairwise FST values were greater among the King Island scrubtit subpopulations than among most Tasmanian scrubtit subpopulations. Genetic diversity was lower and inbreeding coefficients were higher in the King Island scrubtit than all except one of the Tasmanian scrubtit subpopulations. We observed crown baldness in 8/15 King Island scrubtits, but 0/55 Tasmanian scrubtits. Six loci were significantly associated with baldness, including one within the DOCK11 gene which is linked to early feather development. Contemporary gene flow between King Island scrubtit subpopulations is unlikely, with further field monitoring required to quantify the fitness consequences of its small population size, low genetic diversity, and high inbreeding. Evidence-based conservation actions can then be implemented before the taxon goes extinct.

Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

Identification and adaptive evolution analysis of glutaredoxin genes in Populus spp

Glutaredoxin (GRX) is a class of small redox proteins widely involved in cellular redox homeostasis and the regulation of various cellular processes. The role of GRX gene in the differentiation of Populus spp. is rarely reported. We compared the similarities and differences of GRX genes among four sections of poplar using bioinformatics, corrected the annotations of some GRX genes, and focused on analysing their transcript profiling and adaptive evolution in Populus spp. A total of 219 GRX genes were identified in four sections of poplar, among which annotations for 13 genes were corrected. Differences in GRX genes were found between sect. Turanga, represented by P. euphratica, and other poplar sections. Most notably, P. euphratica had the smallest number of duplication events for GRX genes (n = 9) and no tandem duplications, whereas there were >25 duplication events for all other poplars. Furthermore, we detected 18 pairs of GRX genes under positive selection pressure in various sections of poplar, and identified two groups of GRX genes in the Salicaceae that potentially underwent positive selection. Expression profiling results showed that the PtrGRX34 and its orthologous genes were upregulated under stress treatments. In summary, the GRX gene family underwent expansion during poplar differentiation, and some genes underwent rapid evolution during this process, which may be beneficial for Populus spp. to adapt to environmental changes. This study may provide more insights into the molecular mechanisms of Populus spp. adaptation to environmental changes and the adaptive evolution of GRX genes.

Population genomics and conservation management of the threatened black-footed tree-rat (Mesembriomys gouldii) in northern Australia

Genomic diversity is a fundamental component of Earth's total biodiversity, and requires explicit consideration in efforts to conserve biodiversity. To conserve genomic diversity, it is necessary to measure its spatial distribution, and quantify the contribution that any intraspecific evolutionary lineages make to overall genomic diversity. Here, we describe the range-wide population genomic structure of a threatened Australian rodent, the black-footed tree-rat (Mesembriomys gouldii), aiming to provide insight into the timing and extent of population declines across a large region with a dearth of long-term monitoring data. By estimating recent trajectories in effective population sizes at four localities, we confirm widespread population decline across the species' range, but find that the population in the peri-urban area of the Darwin region has been more stable. Based on current sampling, the Melville Island population made the greatest contribution to overall allelic richness of the species, and the prioritisation analysis suggested that conservation of the Darwin and Cobourg Peninsula populations would be the most cost-effective scenario to retain more than 90% of all alleles. Our results broadly confirm current sub-specific taxonomy, and provide crucial data on the spatial distribution of genomic diversity to help prioritise limited conservation resources. Along with additional sampling and genomic analysis from the far eastern and western edges of the black-footed tree-rat distribution, we suggest a range of conservation and research priorities that could help improve black-footed tree-rat population trajectories at large and fine spatial scales, including the retention and expansion of structurally complex habitat patches.

Genomic signals of local adaptation and hybridization in Asian white birch

Disentangling the numerous processes that affect patterns of genome-wide diversity in widespread tree species has important implications for taxonomy, conservation, and forestry. Here, we investigate the population genomic structure of Asian white birch (Betula platyphylla) in China and seek to explain it in terms of hybridization, demography and adaptation. We generate whole genome sequence data from 83 individuals across the species range in China. Combining this with an existing data set for 79 European and Russian white birches, we show a clear distinction between B. pendula and B. platyphylla, which have sometimes been lumped taxonomically. Genomic diversity of B. platyphylla in north-western China and Central Russia is affected greatly by hybridization with B. pendula. Excluding these hybridized populations, B. platyphylla in China has a linear distribution from north-eastern to south-western China, along the edge of the inland mountainous region. Within this distribution, three genetic clusters are found, which we model as long diverged with subsequent episodes of gene flow. Patterns of covariation between allele frequencies and environmental variables in B. platyphylla suggest the role of natural selection in the distribution of diversity at 7609 SNPs of which 3767 were significantly differentiated among the genetic clusters. The putative adaptive SNPs are distributed throughout the genome and span 1633 genic regions. Of these genic regions, 87 were previously identified as candidates for selective sweeps in Eurasian B. pendula. We use the 7609 environmentally associated SNPs to estimate the risk of nonadaptedness for each sequenced B. platyphylla individual under a scenario of future climate change, highlighting areas where populations may be under future threat from rising temperatures.

Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units.

Using genomics to guide seed-sourcing at the right taxonomical level for ecological restoration projects: The complex case of Carex bigelowii s.lat. in Norway

There is a growing demand for ecological restoration using suitable seeds following international standards or national legal demands for local seed-sourcing. However, before selecting the appropriate geographic origin of seeds, it is vital to explore taxonomic complexity related to the focal taxa. We used ddRAD-seq to screen genomic diversity within Carex bigelowii s.lat. focussing on Norway. This species complex is considered a candidate for seeding, but presents considerable morphological, ecological, and genetic variation. The genetic structure of 132 individuals of C. bigelowii s.lat., including Carex nigra as an outgroup, was explored using ordinations, clustering analyses, and a genetic barrier algorithm. Two highly divergent clusters were evident, supporting the recognition of two taxonomic units "C. dacica" and C. bigelowii "subsp. bigelowii". Previously defined seed-sourcing regions for C. bigelowii s.lat. did not consider the known taxonomic complexity, and therefore interpreted the overall genetic structure as seed-sourcing regions, not taxa. We estimated genetic neighborhood sizes within each taxon to be 100-150 km and 300 km, respectively, indicating species-specific delimitations of local seed-sourcing regions. Frequent hybrids, local genetic distinctiveness, and suggested ecotypes add complexity to the discussed seed-sourcing regions. Our results show how genomic screening of diversity and structure in a species complex can alleviate the taxonomic impediment, inform practical questions, and legal requirements related to seed-sourcing, and together with traditional taxonomic work provide necessary information for a sound management of biodiversity.

Assessing signals of selection and historical demography to develop conservation strategies in the Chilean emblematic Araucaria araucana

Loss of genetic diversity reduces the ability of species to evolve and respond to environmental change. Araucaria araucana is an emblematic conifer species from southern South America, with important ethnic value for the Mapuche people (Pehuenche); the Chilean Government has catalogued its conservation status as vulnerable. Climatic fluctuations were potentially a major impact in the genetic variation within many tree species. In this context, the restricted geographic distribution of A. araucana in Chile appears to be a consequence of the Last Glacial Maximum (LGM). During the past two centuries, strong human intervention has also affected the geographical distribution and population sizes of A. araucana. Reduction of population size may cause loss of genetic diversity, which could affect frequency of adaptive loci. The aims of this study were to know the existence of potential loci under selection and populations with genetic, demographic disequilibrium in the Chilean distribution of A. araucana. Based on 268 polymorphic AFLP loci, we have investigated potential loci under selection and genetic, demographic disequilibrium within seven Chilean populations of Araucaria araucana. Correlation of 41 outlier loci with the environmental variables of precipitation and temperature reveals signatures of selection, whereas 227 neutral loci provide estimates of demographic equilibrium and genetic population structure. Three populations are recommended as priorities for conservation.

Characterizing the genomic variation and population dynamics of Plasmodium falciparum malaria parasites in and around Lake Victoria, Kenya

Characterising the genomic variation and population dynamics of Plasmodium falciparum parasites in high transmission regions of Sub-Saharan Africa is crucial to the long-term efficacy of regional malaria elimination campaigns and eradication. Whole-genome sequencing (WGS) technologies can contribute towards understanding the epidemiology and structural variation landscape of P. falciparum populations, including those within the Lake Victoria basin, a region of intense transmission. Here we provide a baseline assessment of the genomic diversity of P. falciparum isolates in the Lake region of Kenya, which has sparse genetic data. Lake region isolates are placed within the context of African-wide populations using Illumina WGS data and population genomic analyses. Our analysis revealed that P. falciparum isolates from Lake Victoria form a cluster within the East African parasite population. These isolates also appear to have distinct ancestral origins, containing genome-wide signatures from both Central and East African lineages. Known drug resistance biomarkers were observed at similar frequencies to those of East African parasite populations, including the S160N/T mutation in the pfap2mu gene, which has been associated with delayed clearance by artemisinin-based combination therapy. Overall, our work provides a first assessment of P. falciparum genetic diversity within the Lake Victoria basin, a region targeting malaria elimination.

Landscape-scale population connectivity in two parasitoid species associated with the spruce budworm: Testing the birdfeeder effect using genetic data

Periodic and spatially synchronous outbreaks of insect pests have dramatic consequences for boreal and sub-boreal forests. Within these multitrophic systems, parasitoids can be stabilizing agents by dispersing toward patches containing higher host density (the so-called birdfeeder effect). However, we know little about the dispersal abilities of parasitoids in continuous forested landscapes, limiting our understanding of the spatiotemporal dynamics of host-parasitoid systems, and constraining our ability to predict forest resilience in the context of global changes. In this study, we investigate the spatial genetic structure and spatial variation in genetic diversity of two important species of spruce budworm larval parasitoids during outbreaks: Apanteles fumiferanae Viereck (Braconidae) and Glypta fumiferanae (Viereck) (Ichneumonidae). Using parasitoids sampled in 2014 from 26 and 29 locations across a study area of 350,000 km² , we identified 1,012 and 992 neutral SNP loci for A. fumiferanae (N = 279 individuals) and G. fumiferanae (N = 382), respectively. Using DAPC, PCA, AMOVA, and IBD analyses, we found evidence for panmixia and high genetic connectivity for both species, matching the previously described genetic structure of the spruce budworm within the same context, suggesting similar effective dispersal during outbreaks and high parasitoid population densities between outbreaks. We also found a significant negative relationship between genetic diversity and latitude for A. fumiferanae but not for G. fumiferanae, suggesting that northern range limits may vary by species within the spruce budworm parasitoid community. These spatial dynamics should be considered when predicting future insect outbreak severities in boreal landscapes.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

A spectral theory for Wright's inbreeding coefficients and related quantities

Wright’s inbreeding coefficient, F_ST, is a fundamental measure in population genetics. Assuming a predefined population subdivision, this statistic is classically used to evaluate population structure at a given genomic locus. With large numbers of loci, unsupervised approaches such as principal component analysis (PCA) have, however, become prominent in recent analyses of population structure. In this study, we describe the relationships between Wright’s inbreeding coefficients and PCA for a model of K discrete populations. Our theory provides an equivalent definition of F_ST based on the decomposition of the genotype matrix into between and within-population matrices. The average value of Wright’s F_ST over all loci included in the genotype matrix can be obtained from the PCA of the between-population matrix. Assuming that a separation condition is fulfilled and for reasonably large data sets, this value of F_ST approximates the proportion of genetic variation explained by the first (K − 1) principal components accurately. The new definition of F_ST is useful for computing inbreeding coefficients from surrogate genotypes, for example, obtained after correction of experimental artifacts or after removing adaptive genetic variation associated with environmental variables. The relationships between inbreeding coefficients and the spectrum of the genotype matrix not only allow interpretations of PCA results in terms of population genetic concepts but extend those concepts to population genetic analyses accounting for temporal, geographical and environmental contexts.

Principal component analysis (PCA) is the most-frequently used approach to describe population genetic structure from large population genomic data sets. In this study, we show that PCA not only estimates ancestries of sampled individuals, but also computes the average value of Wright’s inbreeding coefficient over the loci included in the genotype matrix. Our result shows that inbreeding coefficients and PCA eigenvalues provide equivalent descriptions of population structure. As a consequence, PCA extends the definition of those coefficients beyond the framework of allelic frequencies. We give examples on how F_ST can be computed from ancient DNA samples for which genotypes are corrected for coverage, and in an ecological genomic example where a proportion of genetic variation is explained by environmental variables.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (F_ST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

Scale-dependent signatures of local adaptation in a foundation tree species

Understanding local adaptation is critical for conservation management under rapidly changing environmental conditions. Local adaptation inferred from genotype-environment associations may show different genomic patterns depending on the spatial scale of sampling, due to differences in the slope of environmental gradients and the level of gene flow. We compared signatures of local adaptation across the genome of mountain ash (Eucalyptus regnans) at two spatial scales: A species-wide data set and a topographically-complex subregional data set. We genotyped 367 individual trees at over 3700 single-nucleotide polymorphisms (SNPs), quantified patterns of spatial genetic structure among populations, and used two analytical methods to identify loci associated with at least one of three environmental variables at each spatial scale. Together, the analyses identified 549 potentially adaptive SNPs at the subregion scale, and 435 SNPs at the range-wide scale. A total of 39 genic or near-genic SNPs, associated with 28 genes, were identified at both spatial scales, although no SNP was identified by both methods at both scales. We observed that nongenic regions had significantly higher homozygote excess than genic regions, possibly due to selective elimination of inbred genotypes during stand development. Our results suggest that strong environmental selection occurs in mountain ash, and that the identification of putatively adaptive loci can differ substantially depending on the spatial scale of analyses. We also highlight the importance of multiple adaptive genetic architectures for understanding patterns of local adaptation across large heterogenous landscapes, with comparison of putatively adaptive loci among spatial scales providing crucial insights into the process of adaptation.

LEA 3: Factor models in population genetics and ecological genomics with R

A major objective of evolutionary biology is to understand the processes by which organisms have adapted to various environments, and to predict the response of organisms to new or future conditions. The availability of large genomic and environmental data sets provides an opportunity to address those questions, and the R package LEA has been introduced to facilitate population and ecological genomic analyses in this context. By using latent factor models, the program computes ancestry coefficients from population genetic data and performs genotype-environment association analyses with correction for unobserved confounding variables. In this study, we present new functionalities of LEA, which include imputation of missing genotypes, fast algorithms for latent factor mixed models using multivariate predictors for genotype-environment association studies, population differentiation tests for admixed or continuous populations, and estimation of genetic offset based on climate models. The new functionalities are implemented in version 3.1 and higher releases of the package. Using simulated and real data sets, our study provides evaluations and examples of applications, outlining important practical considerations when analysing ecological genomic data in R.

Landscape genomics of Quercus lobata reveals genes involved in local climate adaptation at multiple spatial scales

Understanding how the environment shapes genetic variation provides critical insight about the evolution of local adaptation in natural populations. At multiple spatial scales and multiple geographic contexts within a single species, such information could address a number of fundamental questions about the scale of local adaptation and whether or not the same loci are involved at different spatial scales or geographic contexts. We used landscape genomic approaches from three local elevational transects and rangewide sampling to (a) identify genetic variation underlying local adaptation to environmental gradients in the California endemic oak, Quercus lobata; (b) examine whether putatively adaptive SNPs show signatures of selection at multiple spatial scales; and (c) map putatively adaptive variation to assess the scale and pattern of local adaptation. Of over 10 k single-nucleotide polymorphisms (SNPs) generated with genotyping-by-sequencing, we found signatures of natural selection by climate or local environment at over 600 SNPs (536 loci), some at multiple spatial scales across multiple analyses. Candidate SNPs identified with gene-environment tests (LFMM) at the rangewide scale also showed elevated associations with climate variables compared to the background at both rangewide and elevational transect scales with gradient forest analysis. Some loci overlap with those detected in other oak species, raising the question of whether the same loci might be involved in local climate adaptation in different congeneric species that inhabit different geographic contexts. Mapping landscape patterns of adaptive versus background genetic variation identified regions of marked local adaptation and suggests nonlinear association of candidate SNPs and environmental variables. Taken together, our results offer robust evidence for novel candidate genes for local climate adaptation at multiple spatial scales.

Functional genetic diversity of domestic and wild American mink (Neovison vison)

The release of domestic organisms to the wild threatens biodiversity because the introduction of domestic genes through interbreeding can negatively impact wild conspecifics via outbreeding depression. In North America, farmed American mink (Neovison vison) frequently escape captivity, yet the impact of these events on functional genetic diversity of wild mink populations is unclear. We characterized domestic and wild mink in Ontario at 17 trinucleotide microsatellites located in functional genes thought to be associated with traits affected by domestication. We found low functional genetic diversity in both mink types, as only four of 17 genes were variable, yet allele frequencies varied widely between captive and wild populations. To determine whether allele frequencies of wild populations were affected by geographic location, we performed redundancy analysis and spatial analysis of principal components on three polymorphic loci (AR, ATN1 and IGF-1). We found evidence to suggest domestic release events are affecting the functional genetic diversity of wild mink, as sPCA showed clear distinctions between wild individuals near mink farms and those located in areas without mink farms. This is further substantiated through RDA, where spatial location was associated with genetic variation of AR, ATN1 and IGF1.

Genomic signatures of seed mass adaptation to global precipitation gradients in sorghum

Seed mass is a key component of adaptation in plants and a determinant of yield in crops. The climatic drivers and genomic basis of seed mass variation remain poorly understood. In the cereal crop Sorghum bicolor, globally-distributed landraces harbor abundant variation in seed mass, which is associated with precipitation in their agroclimatic zones of origin. This study aimed to test the hypothesis that diversifying selection across precipitation gradients, acting on ancestral cereal grain size regulators, underlies seed mass variation in global sorghum germplasm. We tested this hypothesis in a set of 1901 georeferenced and genotyped sorghum landraces, 100-seed mass from common gardens, and bioclimatic precipitation variables. As predicted, 100-seed mass in global germplasm varies significantly among botanical races and is correlated to proxies of the precipitation gradients. With general and mixed linear model genome-wide associations, we identified 29 and 56 of 100 a priori candidate seed size genes with polymorphisms in the top 1% of seed mass association, respectively. Eleven of these genes harbor polymorphisms associated with the precipitation gradient, including orthologs of genes that regulate seed size in other cereals. With FarmCPU, 13 significant SNPs were identified, including one at an a priori candidate gene. Finally, we identified eleven colocalized outlier SNPs associated with seed mass and precipitation that also carry signatures of selection based on F_ST scans and PCAdapt, which represents a significant enrichment. Our findings suggest that seed mass in sorghum was shaped by diversifying selection on drought stress, and can inform genomics-enabled breeding for climate-resilient cereals.

The Effect of Neutral Recombination Variation on Genome Scans for Selection

Recently, there has been an increasing interest in identifying the role that regions of low recombination or inversion play in adaptation of species to local environments. Many examples of groups of adapted genes located within inversions are arising in the literature, in part inspired by theory that predicts the evolution of these so-called “supergenes.” We still, however, have a poor understanding of how genomic heterogeneity, such as varying rates of recombination, may confound signals of selection. Here, I evaluate the effect of neutral inversions and recombination variation on genome scans for selection, including tests for selective sweeps, differentiation outlier tests, and association tests. There is considerable variation among methods in their performance, with some methods being unaffected and some showing elevated false positive signals within a neutral inversion or region of low recombination. In some cases the false positive signal can be dampened or removed, if it is possible to use a quasi-independent set of SNPs to parameterize the model before performing the test. These results will be helpful to those seeking to understand the importance of regions of low recombination in adaptation.

Genomics of end-Pleistocene population replacement in a small mammal

Current species distributions at high latitudes are the product of expansion from glacial refugia into previously uninhabitable areas at the end of the last glaciation. The traditional view of postglacial colonization is that southern populations expanded their ranges into unoccupied northern territories. Recent findings on mitochondrial DNA (mtDNA) of British small mammals have challenged this simple colonization scenario by demonstrating a more complex genetic turnover in Britain during the Pleistocene-Holocene transition where one mtDNA clade of each species was replaced by another mtDNA clade of the same species. Here, we provide evidence from one of those small mammals, the bank vole (Clethrionomys glareolus), that the replacement was genome-wide. Using more than 10 000 autosomal SNPs we found that similar to mtDNA, bank vole genomes in Britain form two (north and south) clusters which admix. Therefore, the genome of the original postglacial colonists (the northern cluster) was probably replaced by another wave of migration from a different continental European population (the southern cluster), and we gained support for this by modelling with approximate Bayesian computation. This finding emphasizes the importance of analysis of genome-wide diversity within species under changing climate in creating opportunities for sophisticated testing of population history scenarios.