When environmental variation is spatially continuous, dispersing individuals move among nearby sites with similar habitat conditions. But as an environmental gradient becomes steeper, gene flow may connect more divergent habitats, and this is predicted to reduce the slope of the adaptive cline that evolves. We compared quantitative genetic divergence of Rana temporaria frog populations along a 2,000-m elevational gradient in eastern Switzerland (new experimental results) with divergence along a 1,550-km latitudinal gradient in Fennoscandia (previously published results). Both studies found significant countergradient variation in larval development rate (i.e., animals from cold climates developed more rapidly). The cline was weaker with elevation than with latitude. Animals collected on both gradients were genotyped at 2,000 single-nucleotide polymorphism markers, revealing that dispersal distance was 30% farther on the latitudinal gradient but 3.9 times greater with respect to environmental conditions on the elevational gradient. A meta-analysis of 19 experimental studies of anuran populations spanning temperature gradients revealed that countergradient variation in larval development, while significant overall, was weaker when measured on steeper gradients. These findings support the prediction that adaptive population divergence is less pronounced, and maladaptation more pervasive, on steep environmental gradients. abstract: When environmental variation is spatially continuous, dispersing individuals move among nearby sites with similar habitat conditions. But as an environmental gradient becomes steeper, gene flow may connect more divergent habitats, and this is predicted to reduce the slope of the adaptive cline that evolves. We compared quantitative genetic divergence of Rana temporaria frog populations along a 2,000-m elevational gradient in eastern Switzerland (new experimental results) with divergence along a 1,550-km latitudinal gradient in Fennoscandia (previously published results). Both studies found significant countergradient variation in larval development rate (i.e., animals from cold climates developed more rapidly). The cline was weaker with elevation than with latitude. Animals collected on both gradients were genotyped at ∼2,000 singlenucleotide polymorphism markers, revealing that dispersal distance was 30% farther on the latitudinal gradient but 3.9 times greater with respect to environmental conditions on the elevational gradient. A meta-analysis of 19 experimental studies of anuran populations spanning temperature gradients revealed that countergradient variation in larval development, while significant overall, was weaker when measured on steeper gradients. These findings support the prediction that adaptive population divergence is less pronounced, and maladaptation more pervasive, on steep environmental gradients.

A pattern of population structure called isolation-by-environment (IBE) evolves when gene flow connecting populations in different habitats is lower than expected. Although IBE is widespread, there is limited information on its magnitude compared with other factors influencing gene flow. We estimated the relative importance of IBE in the frog Rana temporaria in the Swiss Alps, a geographic context in which IBE should be relatively pronounced. The environmental factor potentially causing IBE was the length of the growing season, which is highly correlated with elevation. A sample of 992 individuals from 82 breeding sites were genotyped at 1827 single-nucleotide polymorphism markers; gene flow was estimated in four ways related to F ST , genetic distance, allele sharing, and distance on a population graph. Gravity modeling and random forest regression evaluated the importance of six at-site covariates, 10 between-site covariates, and geographic distance. There was broad agreement among analysis methods and measures of gene flow: isolation-by-distance (IBD) and habitat quality between sites were of highest importance, the elevation and ruggedness of the dispersal path were about half as important, and IBE was about 10-20% as important as IBD. These results combine with other evidence to suggest that population divergence across elevational gradients is underway in amphibians. K E Y W O R D S : Altitude, connectivity, ddRAD, isolation-by-adaptation, topography. Clearly, the alternative habitats must impose divergent selection on traits that influence gene flow (Widmer et al. 2009; Wang and Bradburd 2014; Butlin and Smadja 2018). Selection imposed by phenological differences between habitats has attracted particular 9 6 2

Sex chromosomes in vertebrates range from highly heteromorphic (as in most birds and mammals) to strictly homomorphic (as in many fishes, amphibians, and nonavian reptiles). Reasons for these contrasted evolutionary trajectories remain unclear, but species such as common frogs with polymorphism in the extent of sex chromosome differentiation may potentially deliver important clues. By investigating 92 common frog populations from a wide range of elevations throughout Switzerland, we show that sex chromosome differentiation strongly correlates with alleles at the candidate sex‐determining gene Dmrt1. Y‐specific Dmrt1 haplotypes cluster into two main haplogroups, YA and YB, with a phylogeographic signal that parallels mtDNA haplotypes: YA populations, with mostly well‐differentiated sex chromosomes, occur primarily south of the main alpine ridge that bisects Switzerland, whereas YB populations, with mostly undifferentiated (proto‐)sex chromosomes, occur north of this ridge. Elevation has only a marginal effect, opposing previous suggestions of a major role for climate on sex chromosome differentiation. The Y‐haplotype effect might result from differences in the penetrance of alleles at the sex‐determining locus (such that sex reversal and ensuing X‐Y recombination are more frequent in YB populations), and/or fixation of an inversion on YA (as supported by the empirical observation that YA haplotypes might not recombine in XYA females).

Molecular methods are commonly used to investigate cryptic populations that are difficult to locate or observe directly. The population dynamics of many subterranean organisms have been overlooked, at least in part, as a result of the absence of appropriate molecular markers. Recent studies in African mole-rats have raised questions about the modes of dispersal and mate acquisition. In the present study, we apply a suite of 25 microsatellite markers to test the overground/underground dispersal hypotheses. Using these data, we also apply an approach to estimate population size and look for signal of demographic expansion or contraction. The genetic data suggest that the same breeding population extends between locations (approximately 50 km), with elevated inbreeding coefficients suggestive of some degree of isolation of the urban location. Low genetic differentiation between study sites supports the proposed high levels of vagility of dispersing individuals overground. We find a signal of long-term population decline of Bathyergus suillus in this region. Their adherence to mesic conditions potentially recommends B. suillus to be of utility in monitoring the proposed climate-induced desiccation of the Western Cape. Of potential interest is the discovery of a second divergent population at the rural location, with microsatellite data suggesting contemporary reproductive isolation and a mitochondrial divergence putatively dated at approximately 0.6 Mya.

15The high mountain ranges of western Europe have had a profound effect on 16 the recolonisation of Europe from glacial refugia. The Alps present a 17 particularly interesting case, because they present an absolute barrier to 18 dispersal to most lineages, obstructing recolonisation from multiple refugia in 19 the Italian Alps. Here we investigate the effect of the European Alps on the 20 the phylogeographic history of Rana temporaria across its range in 21Switzerland. Based on partial cytochrome b and COX1 sequences we find two 22 mitochondrial lineages that occur roughly north and south of the alpine ridge 23 bisecting Switzerland, with contact zones between them in the east and west. 24The northern haplogroup falls within the previously identified common western 25European haplogroup, while the southern haplogroup is unique to 26Switzerland. We find that the lineages diverged ~110 kya, approximately the 27 onset of the last glacial maximum, indicative of origins in separate refugia. 28Phylogenetic analyses suggest that the lineages originate from two refugia in 29 northern Italy, and colonised Switzerland via trans-and curcum-alpine routes. 30Our results show that the European Alps is a semi-permeable barrier to 31 dispersal for R. temporaria, and have contributed to the complex 32

BackgroundThe major histocompatibility complex (MHC) is an important component of the vertebrate immune system and is frequently used to characterise adaptive variation in wild populations due to its co-evolution with pathogens. Passerine birds have an exceptionally diverse MHC with multiple gene copies and large numbers of alleles compared to other avian taxa. The Nesospiza bunting species complex (two species on Nightingale Island; one species with three sub-species on Inaccessible Island) represents a rapid adaptive radiation at a small, isolated archipelago, and is thus an excellent model for the study of adaptation and speciation. In this first study of MHC in Nesospiza buntings, we aim to characterize MHCIIß variation, determine the strength of selection acting at this gene region and assess the level of shared polymorphism between the Nesospiza species complex and its putative sister taxon, Rowettia goughensis, from Gough Island.ResultsIn total, 23 unique alleles were found in 14 Nesospiza and 2 R. goughensis individuals encoding at least four presumably functional loci and two pseudogenes. There was no evidence of ongoing selection on the peptide binding region (PBR). Of the 23 alleles, 15 were found on both the islands inhabited by Nesospiza species, and seven in both Nesospiza and Rowettia; indications of shared, ancestral polymorphism. A gene tree of Nesospiza MHCIIß alleles with several other passerine birds shows three highly supported Nesospiza-specific groups. All R. goughensis alleles were shared with Nesospiza, and these alleles were found in all three Nesospiza sequence groups in the gene tree, suggesting that most of the observed variation predates their phylogenetic split.ConclusionsLack of evidence of selection on the PBR, together with shared polymorphism across the gene tree, suggests that population variation of MHCIIß among Nesospiza and Rowettia is due to ancestral polymorphism rather than local selective forces. Weak or no selection pressure could be attributed to low parasite load at these isolated Atlantic islands. The deep divergence between the highly supported Nesospiza-specific sequence Groups 2 and 3, and the clustering of Group 3 close to the distantly related passerines, provide strong support for preserved ancestral polymorphism, and present evidence of one of the rare cases of extensive ancestral polymorphism in birds.

High mountain ranges such as the European Alps and Pyrenees were the last to deglaciate in mainland Europe during glacial minima over the last ~700 ky (Darnault et al., 2012;Ehlers et al., 2018).For taxa adapted to warm climates, the persistence and extent of Alpine ice sheets and the east-west orientation of the mountain ranges hindered recolonization from more southerly refugia in Italy and the Iberian Peninsula into Europe (Taberlet et al., 1998). Coldtolerant organisms exhibit a more complex history. Many expanded their range during glacial periods, some persisted in nunataks (e.g., ice-free rocky protrusions) within the Alps, and trans-Alpine recolonization has occurred in several species (e.g., Oak species, the common shrew, and barred grass snake, and the European stag beetle;

Understanding the rate and extent to which populations can adapt to novel environments at their ecological margins is fundamental to predicting the persistence of biological communities during ongoing and rapid global change. Recent range expansion in response to climate change in the UK butterfly Aricia agestis is associated with the evolution of novel interactions with a larval food plant, and the loss of its ability to use its ancestral larval host species. Using ddRAD analysis of 61210 variable SNPs from 261 females from throughout the UK range of this species, we identify genomic regions at multiple chromosomes that are associated with these evolutionary responses, and their association with demographic history and ecological variation. Gene flow appears widespread throughout the range, despite the apparently fragmented nature of the habitats used by this species. Patterns of haplotype variation between selected and neutral genomic regions suggest that evolution associated with climate adaptation is polygenic, resulting from the independent spread of existing alleles throughout the established range of this species, rather than the colonisation of pre-adapted genotypes from coastal populations. These data suggest that rapid responses to climate change do not depend on the availability of pre-adapted genotypes. Instead, the evolution of novel forms of biotic interaction in Aricia agestis has occurred during range expansion, through the assembly of novel genotypes from alleles from multiple localities.