Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.     

Ecological resistance surfaces predict fine‐scale genetic differentiation in a terrestrial woodland salamander

Landscape genetics has seen tremendous advances since its introduction, but parameterization and optimization of resistance surfaces still poses significant challenges. Despite increased availability and resolution of spatial data, few studies have integrated empirical data to directly represent ecological processes as genetic resistance surfaces. In our study, we determine the landscape and ecological factors affecting gene flow in the western slimy salamander (Plethodon albagula). We used field data to derive resistance surfaces representing salamander abundance and rate of water loss through combinations of canopy cover, topographic wetness, topographic position, solar exposure and distance from ravine. These ecologically explicit composite surfaces directly represent an ecological process or physiological limitation of our organism. Using generalized linear mixed‐effects models, we optimized resistance surfaces using a nonlinear optimization algorithm to minimize model AIC. We found clear support for the resistance surface representing the rate of water loss experienced by adult salamanders in the summer. Resistance was lowest at intermediate levels of water loss and higher when the rate of water loss was predicted to be low or high. This pattern may arise from the compensatory movement behaviour of salamanders through suboptimal habitat, but also reflects the physiological limitations of salamanders and their sensitivity to extreme environmental conditions. Our study demonstrates that composite representations of ecologically explicit processes can provide novel insight and can better explain genetic differentiation than ecologically implicit landscape resistance surfaces. Additionally, our study underscores the fact that spatial estimates of habitat suitability or abundance may not serve as adequate proxies for describing gene flow, as predicted abundance was a poor predictor of genetic differentiation.     

Spatial population genetic structure of a bacterial parasite in close coevolution with its host

Knowledge of a species’ population genetic structure can provide insight into fundamental ecological and evolutionary processes including gene flow, genetic drift and adaptive evolution. Such inference is of particular importance for parasites, as an understanding of their population structure can illuminate epidemiological and coevolutionary dynamics. Here, we describe the population genetic structure of the bacterium Pasteuria ramosa, a parasite that infects planktonic crustaceans of the genus Daphnia. This system has become a model for investigations of host–parasite interactions and represents an example of coevolution via negative frequency‐dependent selection (aka “Red Queen” dynamics). To sample P. ramosa, we experimentally infected a panel of Daphnia hosts with natural spore banks from the sediments of 25 ponds throughout much of the species range in Europe and western Asia. Using 12 polymorphic variable number tandem repeat loci (VNTR loci), we identified substantial genetic diversity, both within and among localities, that was structured geographically among ponds. Genetic diversity was also structured among host genotypes within ponds, although this pattern varied by locality, with P. ramosa at some localities partitioned into distinct host‐specific lineages, and other localities where recombination had shuffled genetic variation among different infection phenotypes. Across the sample range, there was a pattern of isolation by distance, and principal components analysis coupled with Procrustes rotation identified congruence between patterns of genetic variation and geography. Our findings support the hypothesis that Pasteuria is an endemic parasite coevolving closely with its host. These results provide important context for previous studies of this model system and inform hypotheses for future research.     

Detecting a hierarchical genetic population structure: the case study of the Fire Salamander (Salamandra salamandra) in Northern Italy

The multistep method here applied in studying the genetic structure of a low dispersal and philopatric species, such as the Fire Salamander Salamandra salamandra, was proved to be effective in identifying the hierarchical structure of populations living in broad‐leaved forest ecosystems in Northern Italy. In this study, 477 salamander larvae, collected in 28 sampling populations (SPs) in the Prealpine and in the foothill areas of Northern Italy, were genotyped at 16 specie‐specific microsatellites. SPs showed a significant overall genetic variation (Global F_ST = 0.032, P < 0.001). The genetic population structure was assessed by using STRUCTURE 2.3.4. We found two main genetic groups, one represented by SPs inhabiting the Prealpine belt, which maintain connections with those of the Eastern foothill lowland (PEF), and a second group with the SPs of the Western foothill lowland (WF). The two groups were significantly distinct with a Global F_ST of 0.010 (P < 0.001). While the first group showed a moderate structure, with only one divergent SP (Global F_ST = 0.006, P < 0.001), the second group proved more structured being divided in four clusters (Global F_ST = 0.017, P = 0.058). This genetic population structure should be due to the large conurbations and main roads that separate the WF group from the Prealpine belt and the Eastern foothill lowland. The adopted methods allowed the analysis of the genetic population structure of Fire Salamander from wide to local scale, identifying different degrees of genetic divergence of their populations derived from forest fragmentation induced by urban and infrastructure sprawl.     

Functional connectivity in replicated urban landscapes in the land snail (Cornu aspersum)

Urban areas are highly fragmented and thereby exert strong constraints on individual dispersal. Despite this, some species manage to persist in urban areas, such as the garden snail, Cornu aspersum, which is common in cityscapes despite its low mobility. Using landscape genetic approaches, we combined study area replication and multiscale analysis to determine how landscape composition, configuration and connectivity influence snail dispersal across urban areas. At the overall landscape scale, areas with a high percentage of roads decreased genetic differentiation between populations. At the population scale, genetic differentiation was positively linked with building surface, the proportion of borders where wooded patches and roads appeared side by side and the proportion of borders combining wooded patches and other impervious areas. Analyses based on pairwise genetic distances validated the isolation‐by‐distance and isolation‐by‐resistance models for this land snail, with an equal fit to least‐cost paths and circuit‐theory‐based models. Each of the 12 landscapes analysed separately yielded specific relations to environmental features, whereas analyses integrating all replicates highlighted general common effects. Our results suggest that urban transport infrastructures facilitate passive snail dispersal. At a local scale, corresponding to active dispersal, unfavourable habitats (wooded and impervious areas) isolate populations. This work upholds the use of replicated landscapes to increase the generalizability of landscape genetics results and shows how multiscale analyses provide insight into scale‐dependent processes.     

The influence of landscape configuration and environment on population genetic structure in a sedentary passerine: insights from loci located in different genomic regions

The study of the factors structuring genetic variation can help to infer the neutral and adaptive processes shaping the demographic and evolutionary trajectories of natural populations. Here, we analyse the role of isolation by distance (IBD), isolation by resistance (IBR, defined by landscape composition) and isolation by environment (IBE, estimated as habitat and elevation dissimilarity) in structuring genetic variation in 25 blue tit (Cyanistes caeruleus) populations. We typed 1385 individuals at 26 microsatellite loci classified into two groups by considering whether they are located into genomic regions that are actively (TL; 12 loci) or not (NTL; 14 loci) transcribed to RNA. Population genetic differentiation was mostly detected using the panel of NTL. Landscape genetic analyses showed a pattern of IBD for all loci and the panel of NTL, but genetic differentiation estimated at TL was only explained by IBR models considering high resistance for natural vegetation and low resistance for agricultural lands. Finally, the absence for IBE suggests a lack of divergent selection pressures associated with differences in habitat and elevation. Overall, our study shows that markers located in different genomic regions can yield contrasting inferences on landscape‐level patterns of realized gene flow in natural populations.     

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.     

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.     

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high‐latitude B. huntii habitats exhibit low ENS when compared to low‐latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low‐latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.     

Habitat use,but not gene flow,is influenced by human activities in two ecotypes of Egyptian fruit bat (Rousettus aegyptiacus)

Understanding the ecological, behavioural and evolutionary response of organisms to changing environments is of primary importance in a human‐altered world. It is crucial to elucidate how human activities alter gene flow and what are the consequences for the genetic structure of a species. We studied two lineages of the Egyptian fruit bat (Rousettus aegyptiacus) throughout the contact zone between mesic and arid Ecozones in the Middle East to evaluate the species' response to the growing proportion of human‐altered habitats in the desert. We integrated population genetics, morphometrics and movement ecology to analyse population structure, morphological variation and habitat use from GPS‐ or radio‐tagged individuals from both desert and Mediterranean areas. We classified the spatial distribution and environmental stratification by describing physical–geographical conditions and land cover. We analysed this information to estimate patch occupancy and used an isolation‐by‐resistance approach to model gene flow patterns. Our results suggest that lineages from desert and Mediterranean habitats, despite their admixture, are isolated by environment and by adaptation supporting their classification as ecotypes. We found a positive effect of human‐altered habitats on patch occupancy and habitat use of fruit bats by increasing the availability of roosting and foraging areas. While this commensalism promotes the distribution of fruit bats throughout the Middle East, gene flow between colonies has not been altered by human activities. This discrepancy between habitat use and gene flow patterns may, therefore, be explained by the breeding system of the species and modifications of natal dispersal patterns.     

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.     

Intercontinental long‐distance seed dispersal across the Mediterranean Basin explains population genetic structure of a bird‐dispersed shrub

Long‐distance dispersal (LDD) is a pivotal process for plants determining their range of distribution and promoting gene flow among distant populations. Most fleshy‐fruited species rely on frugivorous vertebrates to disperse their seeds across the landscape. While LDD events are difficult to record, a few ecological studies have shown that birds move a sizeable number of ingested seeds across geographic barriers, such as sea straits. The foraging movements of migrant frugivores across distant populations, including those separated by geographic barriers, creates a constant flow of propagules that in turn shapes the spatial distributions of the genetic variation in populations. Here, we have analysed the genetic diversity and structure of 74 populations of Pistacia lentiscus, a fleshy‐fruited shrub widely distributed in the Mediterranean Basin, to elucidate whether the Mediterranean Sea acts as a geographic barrier or alternatively whether migratory frugivorous birds promote gene flow among populations located on both sides of the sea. Our results show reduced genetic distances among populations, including intercontinental populations, and they show a significant genetic structure across an eastern‐western axis. These findings are consistent with known bird migratory routes that connect the European and African continents following a north‐southwards direction during the fruiting season of many fleshy‐fruited plants. Further, approximate Bayesian analysis failed to explain the observed patterns as a result of historical population migrations at the end of Last Glacial Maximum. Therefore, anthropic and/or climatic changes that would disrupt the migratory routes of frugivorous birds might have genetic consequences for the plant species they feed upon.     

Living in two worlds: Evolutionary mechanisms act differently in the native and introduced ranges of an invasive plant

Identifying the factors that influence spatial genetic structure among populations can provide insights into the evolution of invasive plants. In this study, we used the common reed (Phragmites australis), a grass native in Europe and invading North America, to examine the relative importance of geographic, environmental (represented by climate here), and human effects on population genetic structure and its changes during invasion. We collected samples of P. australis from both the invaded North American and native European ranges and used molecular markers to investigate the population genetic structure within and between ranges. We used path analysis to identify the contributions of each of the three factors—geographic, environmental, and human‐related—to the formation of spatial genetic patterns. Genetic differentiation was observed between the introduced and native populations, and their genetic structure in the native and introduced ranges was different. There were strong effects of geography and environment on the genetic structure of populations in the native range, but the human‐related factors manifested through colonization of anthropogenic habitats in the introduced range counteracted the effects of environment. The between‐range genetic differences among populations were mainly explained by the heterogeneous environment between the ranges, with the coefficient 2.6 times higher for the environment than that explained by the geographic distance. Human activities were the primary contributor to the genetic structure of the introduced populations. The significant environmental divergence between ranges and the strong contribution of human activities to the genetic structure in the introduced range suggest that invasive populations of P. australis have evolved to adapt to a different climate and to human‐made habitats in North America.     

Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes

High‐throughput Plasmodium genomic data is increasingly useful in assessing prevalence of clinically important mutations and malaria transmission patterns. Understanding parasite diversity is important for identification of specific human or parasite populations that can be targeted by control programmes, and to monitor the spread of mutations associated with drug resistance. An up‐to‐date understanding of regional parasite population dynamics is also critical to monitor the impact of control efforts. However, this data is largely absent from high‐burden nations in Africa, and to date, no such analysis has been conducted for malaria parasites in Tanzania countrywide. To this end, over 1,000 P. falciparum clinical isolates were collected in 2017 from 13 sites in seven administrative regions across Tanzania, and parasites were genotyped at 1,800 variable positions genome‐wide using molecular inversion probes. Population structure was detectable among Tanzanian P. falciparum parasites, approximately separating parasites from the northern and southern districts and identifying genetically admixed populations in the north. Isolates from nearby districts were more likely to be genetically related compared to parasites sampled from more distant districts. Known drug resistance mutations were seen at increased frequency in northern districts (including two infections carrying pfk13‐R561H), and additional variants with undetermined significance for antimalarial resistance also varied by geography. Malaria Indicator Survey (2017) data corresponded with genetic findings, including average region‐level complexity‐of‐infection and malaria prevalence estimates. The parasite populations identified here provide important information on extant spatial patterns of genetic diversity of Tanzanian parasites, to which future surveys of genetic relatedness can be compared.     

Optimizing the trade‐off between spatial and genetic sampling efforts in patchy populations: towards a better assessment of functional connectivity using an individual‐based sampling scheme

Genetic data are increasingly used in landscape ecology for the indirect assessment of functional connectivity, that is, the permeability of landscape to movements of organisms. Among available tools, matrix correlation analyses (e.g. Mantel tests or mixed models) are commonly used to test for the relationship between pairwise genetic distances and movement costs incurred by dispersing individuals. When organisms are spatially clustered, a population‐based sampling scheme (PSS) is usually performed, so that a large number of genotypes can be used to compute pairwise genetic distances on the basis of allelic frequencies. Because of financial constraints, this kind of sampling scheme implies a drastic reduction in the number of sampled aggregates, thereby reducing sampling coverage at the landscape level. We used matrix correlation analyses on simulated and empirical genetic data sets to investigate the efficiency of an individual‐based sampling scheme (ISS) in detecting isolation‐by‐distance and isolation‐by‐barrier patterns. Provided that pseudo‐replication issues are taken into account (e.g. through restricted permutations in Mantel tests), we showed that the use of interindividual measures of genotypic dissimilarity may efficiently replace interpopulation measures of genetic differentiation: the sampling of only three or four individuals per aggregate may be sufficient to efficiently detect specific genetic patterns in most situations. The ISS proved to be a promising methodological alternative to the more conventional PSS, offering much flexibility in the spatial design of sampling schemes and ensuring an optimal representativeness of landscape heterogeneity in data, with few aggregates left unsampled. Each strategy offering specific advantages, a combined use of both sampling schemes is discussed.     

Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east‐west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.     

Postglacial dispersal patterns and mitochondrial genetic structure of the Pyrenean desman (Galemys pyrenaicus) in the northwestern region of the Iberian Peninsula

The genetic structure of small semiaquatic animals may be influenced by dispersal across both rivers and land. The relative importance of these two modes of dispersal may vary across different species and with ecological conditions and evolutionary periods. The Pyrenean desman (Galemys pyrenaicus) is an endemic mammal of the Iberian Peninsula with a strong phylogeographic structure and semiaquatic habits, thus making it an ideal model to study the effects of river and overland dispersal on its genetic structure. Thanks to different types of noninvasive samples, we obtained an extensive sampling of the Pyrenean desman from the northwestern region of the Iberian Peninsula and sequenced two mitochondrial DNA fragments. We then analyzed, using an isolation‐by‐distance approach, the correlation between phylogenetic distances and geographical distances measured along both river networks and land to infer the relative importance of river and overland dispersal. We found that the correlations in the whole area and in a large basin were consistent with an effect of overland dispersal, which may be due to the postglacial colonization of new territories using terrestrial corridors and, possibly, a more extensive fluvial network that may have been present during the Holocene. However, in a small basin, likely to be less influenced by the impact of ancient postglacial dispersal, the correlations suggested significant overall effects of both overland and river dispersal, as expected for a semiaquatic mammal. Therefore, different scales and geographical regions reflect different aspects of the evolutionary history and ecology of this semiaquatic species using this isolation‐by‐distance method. The results we obtained may have crucial implications for the conservation of the Pyrenean desman because they reinforce the importance of interbasin dispersal for this species in the studied area and the need to protect the whole riverine ecosystem, including rivers, upland streams and terrestrial corridors between basins.     

Comparative landscape genetics reveals the evolution of viviparity reduces genetic connectivity in fire salamanders

Evolutionary changes in reproductive mode may affect co‐evolving traits, such as dispersal, although this subject remains largely underexplored. The shift from aquatic oviparous or larviparous reproduction to terrestrial viviparous reproduction in some amphibians entails skipping the aquatic larval stage and, thus, greater independence from water. Accordingly, amphibians exhibiting terrestrial viviparous reproduction may potentially disperse across a wider variety of suboptimal habitats and increase population connectivity in fragmented landscapes compared to aquatic‐breeding species. We investigated this hypothesis in the fire salamander (Salamandra salamandra), which exhibits both aquatic‐ (larviparity) and terrestrial‐breeding (viviparity) strategies. We genotyped 426 larviparous and 360 viviparous adult salamanders for 13 microsatellite loci and sequenced a mitochondrial marker for 133 larviparous and 119 viviparous individuals to compare population connectivity and landscape resistance to gene flow within a landscape genetics framework. Contrary to our predictions, viviparous populations exhibited greater differentiation and reduced genetic connectivity compared to larviparous populations. Landscape genetic analyses indicate viviparity may be partially responsible for this pattern, as water courses comprised a significant barrier only in viviparous salamanders, probably due to their fully terrestrial life cycle. Agricultural areas and, to a lesser extent, topography also decreased genetic connectivity in both larviparous and viviparous populations. This study is one of very few to explicitly demonstrate the evolution of a derived reproductive mode affects patterns of genetic connectivity. Our findings open avenues for future research to better understand the eco‐evolutionary implications underlying the emergence of terrestrial reproduction in amphibians.