Quantifying the roles of ecology and geography in spatial genetic divergence

Investigating the properties of ecological landscapes that influence gene flow among populations can provide key insights into the earliest stages of biological divergence. Both ecological and geographical factors can reduce gene flow, which can lead to population divergence, but we know little of the relative strengths of these phenomena in nature. Here, we use a novel application of structural equation modelling to quantify the contributions of ecological and geographical isolation to spatial genetic divergence in 17 species of Anolis lizards. Our comparative analysis shows that although both processes contributed significantly, geographical isolation explained substantially more genetic divergence than ecological isolation (36.3 vs. 17.9% of variance respectively), suggesting that despite the proposed ubiquity of ecological divergence, non‐ecological factors play the dominant role in the evolution of spatial genetic divergence.     

Different processes shape the patterns of divergence in the nuclear and chloroplast genomes of a relict tree species in East Asia

Isolation by spatial distance (IBD), environment (IBE), and historical climatic instability (IBI) are three common processes assessed in phylogeographic and/or landscape genetic studies. However, the relative contributions of these three processes with respect to spatial genetic patterns have seldom been compared. Moreover, whether the relative contribution differs in different regions or when assessed using different genetic markers has rarely been reported. Lindera obtusiloba has been found to have two sister genetic clades of chloroplast (cpDNA) and nuclear microsatellite (nSSR), both of which show discontinuous distribution in northern and southern East Asia. In this study, we used the Mantel test and multiple matrix regression with randomization (MMRR) to determine the relative contributions of IBD, IBE, and IBI with respect to L. obtusiloba populations. Independent Mantel tests and MMRR calculations were conducted for two genetic data sets (cpDNA and nSSR) and for different regions (the overall species range, and northern and southern subregions of the range). We found a significant IBI pattern in nSSR divergence for all assessed regions, whereas no clear IBI pattern was detected with respect to cpDNA. In contrast, significant (or marginal) divergent IBD patterns were detected for cpDNA in all regions, whereas although a significant IBE was apparent with respect to the overall range, the effect was not detected in the two subregions. The differences identified in nSSR and cpDNA population divergence may be related to differences in the heredity and ploidy of the markers. Compared with the southern region, the northern region showed less significant correlation patterns, which may be related to the shorter population history and restricted population range. The findings of this study serve to illustrate that comparing between markers or regions can contribute to gaining a better understanding the population histories of different genomes or within different regions of a species' range.     

The role of geography and ecology in shaping repeated patterns of morphological and genetic differentiation between European minnows (Phoxinus phoxinus) from the Pyrenees and the Alps

Neutral and selective processes can drive repeated patterns of evolution in different groups of populations experiencing similar ecological gradients. In this paper, we used a combination of nuclear and mitochondrial DNA markers, as well as geometric morphometrics, to investigate repeated patterns of morphological and genetic divergence of European minnows in two mountain ranges: the Pyrenees and the Alps. European minnows (Phoxinus phoxinus) are cyprinid fish inhabiting most freshwater bodies in Europe, including those in different mountain ranges that could act as major geographical barriers to gene flow. We explored patterns of P. phoxinus phenotypic and genetic diversification along a gradient of altitude common to the two mountain ranges, and tested for isolation by distance (IBD), isolation by environment (IBE) and isolation by adaptation (IBA). The results indicated that populations from the Pyrenees and the Alps belong to two well differentiated, reciprocally monophyletic mtDNA lineages. Substantial genetic differentiation due to geographical isolation within and between populations from the Pyrenees and the Alps was also found using rapidly evolving AFLPs markers (isolation by distance or IBD), as well as morphological differences between mountain ranges. Also, morphology varied strongly with elevation and so did genetic differentiation to a lower extent. Despite moderate evidence for IBE and IBA, and therefore of repeated evolution, substantial population heterogeneity was found at the genetic level, suggesting that selection and population specific genetic drift act in concert to affect genetic divergence.     

Isolation by distance,resistance and/or clusters? Lessons learned from a forest‐dwelling carnivore inhabiting a heterogeneous landscape

Landscape genetics provides a valuable framework to understand how landscape features influence gene flow and to disentangle the factors that lead to discrete and/or clinal population structure. Here, we attempt to differentiate between these processes in a forest‐dwelling small carnivore [European pine marten (Martes martes)]. Specifically, we used complementary analytical approaches to quantify the spatially explicit genetic structure and diversity and analyse patterns of gene flow for 140 individuals genotyped at 15 microsatellite loci. We first used spatially explicit and nonspatial Bayesian clustering algorithms to partition the sample into discrete clusters and evaluate hypotheses of ‘isolation by barriers’ (IBB). We further characterized the relationships between genetic distance and geographical (‘isolation by distance’, IBD) and ecological distances (‘isolation by resistance’, IBR) obtained from optimized landscape models. Using a reciprocal causal modelling approach, we competed the IBD, IBR and IBB hypotheses with each other to unravel factors driving population genetic structure. Additionally, we further assessed spatially explicit indices of genetic diversity using sGD across potentially overlapping genetic neighbourhoods that matched the inferred population structure. Our results revealed a complex spatial genetic cline that appears to be driven jointly by IBD and partial barriers to gene flow (IBB) associated with poor habitat and interspecific competition. Habitat loss and fragmentation, in synergy with past overharvesting and possible interspecific competition with sympatric stone marten (Martes foina), are likely the main factors responsible for the spatial genetic structure we observed. These results emphasize the need for a more thorough evaluation of discrete and clinal hypotheses governing gene flow in landscape genetic studies, and the potential influence of different limiting factors affecting genetic structure at different spatial scales.     

Evolutionary insights into Rhinolophus episcopus (Chiroptera,Rhinolophidae) in China: Isolation by distance,environment, or sensory system?

Evolutionary processes can be influenced by several factors, such as geographic isolation, environmental selection, and sensory variation. For most nocturnal bats, echolocation is the primary sensory system used to prey and communicate, and plays important roles in chiropteran diversification and evolution. Understanding the relative contribution of geography, the environment, and this sensory system to population genetic divergence can elucidate the processes involved in bat incipient speciation and evolution. In this study, we collected spatial and environmental information, echolocation calls, as well as the previously published genetic data (six microsatellite loci and the mitochondrial cytochrome b gene) of widely distributed Rhinolophus episcopus populations to test three hypotheses for nuclear and mitochondrial divergence (isolation by distance, isolation by environment, and isolation by sensory variation) and unveil the factors that drive intraspecific genetic differentiation. The moderate level of nuclear differentiation was correlated with geographic/spatial distance and acoustic variation, whereas the relatively high level of mitochondrial differentiation was mainly associated with acoustic divergence. No significant correlation was observed between genetic divergence and environmental variables. Among the three factors, acoustic divergence explained the highest percentage of both nuclear and mitochondrial divergence. Thus, our results indicate that sensory variation may have played important roles in driving population isolation early in bat speciation, which is consistent with the hypothesis of isolation by sensory variation. Our study emphasizes the need to consider more factors, especially sensory traits, and combine multiple statistical methods in landscape genetic studies to test their potential contributions to driving population divergence.     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.     

Evolutionary history of endangered and relict tree species Dipteronia sinensis in response to geological and climatic events in the Qinling Mountains and adjacent areas

Geological and climatic events are considered to profoundly affect the evolution and lineage divergence of plant species. However, the evolutionary histories of tree species that have responded to past geological and climate oscillations in central China''s mountainous areas remain mostly unknown. In this study, we assessed the evolutionary history of the endangered and relict tree species Dipteronia sinensis in the Qinling Mountains (QM) and adjacent areas in East Asia based on variations in the complete chloroplast genomes (cpDNA) and reduced‐genomic scale single nucleotide polymorphisms (SNPs). Population structure and phylogenetic analysis based on the cpDNA variations suggested that D. sinensis could be divided into two intraspecific genetic lineages in the eastern and western sides of the QM (EQM and WQM, respectively) in East Asia. Molecular dating suggested that the intraspecific divergence of D. sinensis occurred approximately 39.2 million years ago during the later Paleogene. It was significantly correlated with the orogeny of the QM, where the formation of this significant geographic barrier in the region may have led to the divergence of independent lineages. Bayesian clustering and demographic analysis showed that intraspecific gene flow was restricted between the EQM and WQM lineages. Isolation‐with‐migration analysis indicated that the two genetic lineages experienced significant demographic expansions after the Pleistocene ice ages. However, the genetic admixture was determined in some populations between the two lineages by the large scale of SNP variations due to DNA incompatibility, the large significant population size, and rapid gene flow of nuclear DNA markers. Our results suggest that two different conservation and management units should be constructed for D. sinensis in the EQM and WQM areas. These findings provide novel insights into the unprecedented effects of tectonic changes and climatic oscillations on lineage divergence and plant population evolution in the QM and adjacent areas in East Asia.     

Correlates of genetic differentiation and isolation by distance in 17 congeneric Silene species

The contemporary pattern of intraspecific genetic variation can indicate the relative role of gene flow and local differentiation in shaping the evolutionary history and future trajectory of a species. To assess the recent influence of contrasting life history and demographic characteristics on genetic structure within a group of closely related species, patterns of genetic differentiation (F(ST) and related statistics) and isolation by distance (IBD) were compared among 17 congeneric herbaceous plant species. Data came from 35 published studies of 16 species, and a previously unpublished analysis of chloroplast genetic variation in the rare endemic Silene rotundifolia. Among-population genetic variance was most strongly influenced by the type of genetic marker used; cytoplasmic markers showed larger values than allozyme and anonymous nuclear markers. Other independently significant factors were geographical range size and, for allozyme studies, reproductive system; in particular, endemism and hermaphroditism were associated with higher among-population genetic variance, whereas large native geographical range and dioecy were associated with lower among-population variance. Over equivalent spatial scales, dioecious populations also showed weaker IBD than hermaphrodites, perhaps because increased population transience and/or variance in the spatial pattern of gene flow are more closely associated with dioecy in this genus. Invasive populations had both highly variable among-population genetic variance, and no evidence for IBD, consistent with nonequilibrium conditions. Other analysed factors including predominant pollinator had no discernable influence on genetic structure or patterns of IBD. In general, this comparative approach appears to be valuable for synthesizing the complementary information provided by F-statistics and IBD, and for indicating the relative importance of particular biological factors in shaping genetic variation within different species of a closely related plant group.     

Divergent selection along elevational gradients promotes genetic and phenotypic disparities among small mammal populations

Species distributed along mountain slopes, facing contrasting habitats in short geographic scale, are of particular interest to test how ecologically based divergent selection promotes phenotypic and genetic disparities as well as to assess isolation‐by‐environment mechanisms. Here, we conduct the first broad comparative study of phenotypic variation along elevational gradients, integrating a large array of ecological predictors and disentangling population genetic driver processes. The skull form of nine ecologically distinct species distributed over a large altitudinal range (100–4200 m) was compared to assess whether phenotypic divergence is a common phenomenon in small mammals and whether it shows parallel patterns. We also investigated the relative contribution of biotic (competition and predation) and abiotic parameters on phenotypic divergence via mixed models. Finally, we assessed the population genetic structure of a rodent species (Niviventer confucianus) via analysis of molecular variance and F_ST along three mountain slopes and tested the isolation‐by‐environment hypothesis using Mantel test and redundancy analysis. We found a consistent phenotypic divergence and marked genetic structure along elevational gradients; however, the species showed mixed patterns of size and skull shape trends across mountain zones. Individuals living at lower altitudes differed greatly in both phenotype and genotype from those living at high elevations, while middle‐elevation individuals showed more intermediate forms. The ecological parameters associated with phenotypic divergence along elevation gradients are partly related to species' ecological and evolutionary constraints. Fossorial and solitary animals are mainly affected by climatic factors, while terrestrial and more gregarious species are influenced by biotic and abiotic parameters. A novel finding of our study is that predator richness emerged as an important factor associated with the intraspecific diversification of the mammalian skull along elevational gradients, a previously overlooked parameter. Population genetic structure was mainly driven by environmental heterogeneity along mountain slopes, with no or a week spatial effect, fitting the isolation‐by‐environment scenario. Our study highlights the strong and multifaceted effects of heterogeneous steep habitats and ecologically based divergent selective forces in small mammal populations.     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Evolution along the Great Rift Valley: phenotypic and genetic differentiation of East African white‐eyes (Aves,Zosteropidae)

The moist and cool cloud forests of East Africa represent a network of isolated habitats that are separated by dry and warm lowland savannah, offering an opportunity to investigate how strikingly different selective regimes affect species diversification. Here, we used the passerine genus Zosterops (white‐eyes) from this region as our model system. Species of the genus occur in contrasting distribution settings, with geographical mountain isolation driving diversification, and savannah interconnectivity preventing differentiation. We analyze (1) patterns of phenotypic and genetic differentiation in high‐ and lowland species (different distribution settings), (2) investigate the potential effects of natural selection and temporal and spatial isolation (evolutionary drivers), and (3) critically review the taxonomy of this species complex. We found strong phenotypic and genetic differentiation among and within the three focal species, both in the highland species complex and in the lowland taxa. Altitude was a stronger predictor of phenotypic patterns than the current taxonomic classification. We found longitudinal and latitudinal phenotypic gradients for all three species. Furthermore, wing length and body weight were significantly correlated with altitude and habitat type in the highland species Z. poliogaster. Genetic and phenotypic divergence showed contrasting inter‐ and intraspecific structures. We suggest that the evolution of phenotypic characters is mainly driven by natural selection due to differences in the two macro‐habitats, cloud forest and savannah. In contrast, patterns of neutral genetic variation appear to be rather driven by geographical isolation of the respective mountain massifs. Populations of the Z. poliogaster complex, as well as Z. senegalensis and Z. abyssinicus, are not monophyletic based on microsatellite data and have higher levels of intraspecific differentiation compared to the currently accepted species.     

Genomics of the divergence continuum in an African plant biodiversity hotspot,I: drivers of population divergence in Restio capensis (Restionaceae)

Understanding the drivers of population divergence, speciation and species persistence is of great interest to molecular ecology, especially for species‐rich radiations inhabiting the world's biodiversity hotspots. The toolbox of population genomics holds great promise for addressing these key issues, especially if genomic data are analysed within a spatially and ecologically explicit context. We have studied the earliest stages of the divergence continuum in the Restionaceae, a species‐rich and ecologically important plant family of the Cape Floristic Region (CFR) of South Africa, using the widespread CFR endemic Restio capensis (L.) H.P. Linder & C.R. Hardy as an example. We studied diverging populations of this morphotaxon for plastid DNA sequences and >14 400 nuclear DNA polymorphisms from Restriction site Associated DNA (RAD) sequencing and analysed the results jointly with spatial, climatic and phytogeographic data, using a Bayesian generalized linear mixed modelling (GLMM) approach. The results indicate that population divergence across the extreme environmental mosaic of the CFR is mostly driven by isolation by environment (IBE) rather than isolation by distance (IBD) for both neutral and non‐neutral markers, consistent with genome hitchhiking or coupling effects during early stages of divergence. Mixed modelling of plastid DNA and single divergent outlier loci from a Bayesian genome scan confirmed the predominant role of climate and pointed to additional drivers of divergence, such as drift and ecological agents of selection captured by phytogeographic zones. Our study demonstrates the usefulness of population genomics for disentangling the effects of IBD and IBE along the divergence continuum often found in species radiations across heterogeneous ecological landscapes.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Isolation by instability: Historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present‐day tropical biota. Using genome‐wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post‐Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as “isolation by instability,” highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.     

Drivers of population genetic differentiation in the wild: isolation by dispersal limitation,isolation by adaptation and isolation by colonization

Empirical population genetic studies have been dominated by a neutralist view, according to which gene flow and drift are the main forces driving population genetic structure in nature. The neutralist view in essence describes a process of isolation by dispersal limitation (IBDL) that generally leads to a pattern of isolation by distance (IBD). Recently, however, conceptual frameworks have been put forward that view local genetic adaptation as an important driver of population genetic structure. Isolation by adaptation (IBA) and monopolization (M) posit that gene flow among natural populations is reduced as a consequence of local genetic adaptation. IBA stresses that effective gene flow is reduced among habitats that show dissimilar ecological characteristics, leading to a pattern of isolation by environment. In monopolization, local genetic adaptation of initial colonizing genotypes results in a reduction in gene flow that fosters the persistence of founder effects. Here, we relate these different processes driving landscape genetic structure to patterns of IBD and isolation by environment (IBE). We propose a method to detect whether IBDL, IBA and M shape genetic differentiation in natural landscapes by studying patterns of variation at neutral and non‐neutral markers as well as at ecologically relevant traits. Finally, we reinterpret a representative number of studies from the recent literature by associating patterns to processes and identify patterns associated with local genetic adaptation to be as common as IBDL in structuring regional genetic variation of populations in the wild. Our results point to the importance of quantifying environmental gradients and incorporating ecology in the analysis of population genetics.     

Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close‐by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (F_ST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south‐western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub‐populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.     

Environmental selection is a main driver of divergence in house sparrows (Passer domesticus) in Romania and Bulgaria

Both neutral and adaptive evolutionary processes can cause population divergence, but their relative contributions remain unclear. We investigated the roles of these processes in population divergence in house sparrows (Passer domesticus) from Romania and Bulgaria, regions characterized by high landscape heterogeneity compared to Western Europe. We asked whether morphological divergence, complemented with genetic data in this human commensal species, was best explained by environmental variation, geographic distance, or landscape resistance—the effort it takes for an individual to disperse from one location to the other—caused by either natural or anthropogenic barriers. Using generalized dissimilarity modeling, a matrix regression technique that fits biotic beta diversity to both environmental predictors and geographic distance, we found that a small set of climate and vegetation variables explained up to ~30% of the observed divergence, whereas geographic and resistance distances played much lesser roles. Our results are consistent with signals of selection on morphological traits and of isolation by adaptation in genetic markers, suggesting that selection by natural environmental conditions shapes population divergence in house sparrows. Our study thus contributes to a growing body of evidence that adaptive evolution may be a major driver of diversification.     

Isolation by environment in White‐breasted Nuthatches (Sitta carolinensis) of the Madrean Archipelago sky islands: a landscape genomics approach

Understanding landscape processes driving patterns of population genetic differentiation and diversity has been a long‐standing focus of ecology and evolutionary biology. Gene flow may be reduced by historical, ecological or geographic factors, resulting in patterns of isolation by distance (IBD) or isolation by environment (IBE). Although IBE has been found in many natural systems, most studies investigating patterns of IBD and IBE in nature have used anonymous neutral genetic markers, precluding inference of selection mechanisms or identification of genes potentially under selection. Using landscape genomics, the simultaneous study of genomic and ecological landscapes, we investigated the processes driving population genetic patterns of White‐breasted Nuthatches (Sitta carolinensis) in sky islands (montane forest habitat islands) of the Madrean Archipelago. Using more than 4000 single nucleotide polymorphisms and multiple tests to investigate the relationship between genetic differentiation and geographic or ecological distance, we identified IBE, and a lack of IBD, among sky island populations of S. carolinensis. Using three tests to identify selection, we found 79 loci putatively under selection; of these, seven matched CDS regions in the Zebra Finch. The loci under selection were highly associated with climate extremes (maximum temperature of warmest month and minimum precipitation of driest month). These results provide evidence for IBE – disentangled from IBD – in sky island vertebrates and identify potential adaptive genetic variation.     

History,geography and host use shape genomewide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei)

Divergent host use has long been suspected to drive population differentiation and speciation in plant‐feeding insects. Evaluating the contribution of divergent host use to genetic differentiation can be difficult, however, as dispersal limitation and population structure may also influence patterns of genetic variation. In this study, we use double‐digest restriction‐associated DNA (ddRAD) sequencing to test the hypothesis that divergent host use contributes to genetic differentiation among populations of the redheaded pine sawfly (Neodiprion lecontei), a widespread pest that uses multiple Pinus hosts throughout its range in eastern North America. Because this species has a broad range and specializes on host plants known to have migrated extensively during the Pleistocene, we first assess overall genetic structure using model‐based and model‐free clustering methods and identify three geographically distinct genetic clusters. Next, using a composite‐likelihood approach based on the site frequency spectrum and a novel strategy for maximizing the utility of linked RAD markers, we infer the population topology and date divergence to the Pleistocene. Based on existing knowledge of Pinus refugia, estimated demographic parameters and patterns of diversity among sawfly populations, we propose a Pleistocene divergence scenario for N. lecontei. Finally, using Mantel and partial Mantel tests, we identify a significant relationship between genetic distance and geography in all clusters, and between genetic distance and host use in two of three clusters. Overall, our results indicate that Pleistocene isolation, dispersal limitation and ecological divergence all contribute to genomewide differentiation in this species and support the hypothesis that host use is a common driver of population divergence in host‐specialized insects.