Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long-term persistence of populations, understanding how these processes affect variation over large-scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700-km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction-site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post-glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate-mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes. K E Y W O R D S adaptive divergence, amphibians, divergent selection, genetic drift, range expansion, small populations 1 | INTRODUC TI ON Species with wide geographical distributions are powerful model systems for investigating how selection and neutral processes shape genetic and phenotypic variation across heterogeneousenvironments, which may be affected by both contemporary processes as well as the demographic history of populations. However, disentangling the relative role of these processes along spatial gradients is challenging because environmental, genetic and phenotypic variation often co-vary in the same direction. When dispersal

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene–environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.

Understanding the molecular basis of adaption is one of the central goals in evolutionary biology and when investigated across sister species it can provide detailed insight into the mechanisms of speciation. Here, we sequence the genomes of 34 individuals from three closely related grouse species in order to uncover the genomic architecture of speciation and the genes involved in adaptation. We identify 6 regions, containing 7 genes that show lineage specific signs of differential selection across the species. These genes are involved in a variety of cell processes ranging from stress response to neural, gut, olfactory and limb development. Genome wide neutrality test statistics reveal a strong signal of population expansion acting across the genomes. Additionally, we uncover a 3.5 Mb region on chromosome 20 that shows considerably lower levels of differentiation across the three grouse lineages, indicating possible action of uniform selection in this region.

Ectotherm development rates often show adaptive divergence along climatic gradients, but the genetic basis for this variation is rarely studied. Here, we investigated the genetic basis for phenotypic variation in larval development in the moor frog Rana arvalis from five regions along a latitudinal gradient from Germany to northern Sweden. We focused on the C/EBP‐1 gene, a transcription factor associated with larval development time. Allele frequencies at C/EBP‐1 varied strongly among geographical regions. Overall, the distribution of alleles along the gradient was in concordance with the dual post‐glacial colonization routes into Scandinavia, with a large number of alleles exclusively present along the southern colonization route. Only three of 38 alleles were shared between the routes. Analysis of contemporary selection on C/EBP‐1 showed divergent selection among the regions, probably reflecting adaptation to the local environmental conditions, although this was especially strong between southern and northern regions coinciding also with lineages from different colonization routes. Overall, the C/EBP‐1 gene has historically been under purifying selection, but two specific amino acid positions showed significant signals of positive selection. These positions showed divergence between southern and northern regions, and we suggest that they are functionally involved in the climatic adaptation of larval development. Using phenotypic data from a common garden experiment, we found evidence for specific C/EBP‐1 alleles being correlated with larval development time, suggesting a functional role in adaptation of larval development to large‐scale climatic variation.

In organisms in which individuals mate multiply, knowledge of the proportion of offspring sired by the last male to mate (P2) under field conditions is important for a thorough understanding of how sexual selection works in nature. In many insect groups, pronounced intraspecific variation in P2 is commonplace. Interestingly, however, in stark contrast to these observations, compilation of P2 data in dragonflies and damselflies (Odonata) indicates that a high P2, seldom below 0.95, is a feature of this taxon. Here we used double digest restriction-site associated DNA sequencing to generate a panel of single nucleotide polymorphisms (SNPs) with which we could determine paternity and estimate values of P2 in the offspring of 19 field-collected pairs of the emerald damselfly Lestes sponsa. We also estimated the relationship between P2 and male genital shape of 16 males using geometric morphometric analysis. P2 was variable (range = 0.0–1.0; mean = 0.5), and there was a marginally non-significant (P = 0.069) relationship between genital shape and P2, suggesting that males with a high P2 had an aedeagus with a broader tip. We suggest that the high P2-values reported in past studies in Odonata are partly due to the methods used to infer paternity. Use of SNPs to determine patterns of paternity and P2 in odonates is needed for a better appraisal of fitness in odonates, and would open many future avenues for use of odonates as models of sexual selection.