Landscape Genetics Framework Research Articles

Landscape genetics of a sub-alpine toad: climate change predicted to induce upward range shifts via asymmetrical migration corridors.

Climate change is expected to have a major hydrological impact on the core breeding habitat and migration corridors of many amphibians in the twenty-first century. The Yosemite toad (Anaxyrus canorus) is a species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. Despite living entirely on federal lands, it has recently faced severe extirpations, yet our understanding of climatic influences on population connectivity is limited. In this study, we used a previously published double-digest RADseq dataset along with numerous remotely sensed habitat features in a landscape genetics framework to answer two primary questions in Yosemite National Park: (1) Which fine-scale climate, topographic, soil, and vegetation features most facilitate meadow connectivity? (2) How is climate change predicted to influence both the magnitude and net asymmetry of genetic migration? We developed an approach for simultaneously modeling multiple toad migration paths, akin to circuit theory, except raw environmental features can be separately considered. Our workflow identified the most likely migration corridors between meadows and used the unique cubist machine learning approach to fit and forecast environmental models of connectivity. We identified the permuted modeling importance of numerous snowpack-related features, such as runoff and groundwater recharge. Our results highlight the importance of considering phylogeographic structure, and asymmetrical migration in landscape genetics. We predict an upward elevational shift for this already high-elevation species, as measured by the net vector of anticipated genetic movement, and a north-eastward shift in species distribution via the network of genetic migration corridors across the park.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour.

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.

Modelling multilocus selection in an individual-based, spatially-explicit landscape genetics framework.

We implemented multilocus selection in a spatially-explicit, individual-based framework that enables multivariate environmental gradients to drive selection in many loci as a new module for the landscape genetics programs, CDPOP and CDMetaPOP. Our module simulates multilocus selection using a linear additive model, providing a flexible platform to evaluate a wide range of genotype-environment associations. Importantly, the module allows simulation of selection in any number of loci under the influence of any number of environmental variables. We validated the module with individual-based selection simulations under Wright-Fisher assumptions. We then evaluated results for simulations under a simple landscape selection model. Next, we simulated individual-based multilocus selection across a complex selection landscape with three loci linked to three different environmental variables. Finally, we demonstrated how the program can be used to simulate multilocus selection under varying selection strengths across different levels of gene flow in a landscape genetics framework. This new module provides a valuable addition to the study of landscape genetics, allowing for explicit evaluation of the contributions and interactions between gene flow and selection-driven processes across complex, multivariate environmental and landscape conditions.

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.

From inland to the coast: Spatial and environmental signatures on the genetic diversity in the colonization of the South Atlantic Coastal Plain

The process of colonization and range expansion to novel environments involves both demographic and selective processes that can be traceable in the genetic diversity of the organisms. Nevertheless, it can be difficult to disentangle the signatures of the demographic and selective in the current genetic diversity of the populations. In this work, we use a landscape genetics framework to investigate the association of spatial and climatic variables with the genetic diversity and differentiation patterns of wild Petunia populations involved in a process of coastal colonization during the last 400 thousand years. Over 300 individuals from 17 populations were genotyped using ten microsatellite loci. Our results suggest that the genetic diversity is higher in populations located at the center of the species range with a decline toward the edges, and that gene flow follows an inland-to-coastal and central-to-peripheral dynamic that parallels the colonization history of this coastal lineage. We identify high levels of genetic differentiation between inland and coastal populations. As part of the inland-coastal genetic differentiation, we find signals of isolation by environment associated with differences in extreme temperature regimes. The differentiation of edge populations along the coast is associated with precipitation seasonality. Our results are robust to controls for spatial and historical divergence in the observed genetic differentiation. We conclude that the ecological differentiation process during the colonization of the South Atlantic Coastal Plain (SACP) was likely facilitated by genetic enrichment resulting from gene flow from central to marginal populations, as well as by rapid allele fixation resulting from serial founder effects during the range expansion along the coast.

Elevation, Not Deforestation, Promotes Genetic Differentiation in a Pioneer Tropical Tree.

The regeneration of disturbed forest is an essential part of tropical forest ecology, both with respect to natural disturbance regimes and large-scale human-mediated logging, grazing, and agriculture. Pioneer tree species are critical for facilitating the transition from deforested land to secondary forest because they stabilize terrain and enhance connectivity between forest fragments by increasing matrix permeability and initiating disperser community assembly. Despite the ecological importance of early successional species, little is known about their ability to maintain gene flow across deforested landscapes. Utilizing highly polymorphic microsatellite markers, we examined patterns of genetic diversity and differentiation for the pioneer understory tree Miconia affinis across the Isthmus of Panama. Furthermore, we investigated the impact of geographic distance, forest cover, and elevation on genetic differentiation among populations using circuit theory and regression modeling within a landscape genetics framework. We report marked differences in historical and contemporary migration rates and moderately high levels of genetic differentiation in M. affinis populations across the Isthmus of Panama. Genetic differentiation increased significantly with elevation and geographic distance among populations; however, we did not find that forest cover enhanced or reduced genetic differentiation in the study region. Overall, our results reveal strong dispersal for M. affinis across human-altered landscapes, highlighting the potential use of this species for reforestation in tropical regions. Additionally, this study demonstrates the importance of considering topography when designing programs aimed at conserving genetic diversity within degraded tropical landscapes.

Beyond species distribution modeling: A landscape genetics approach to investigating range shifts under future climate change

Understanding how biodiversity will respond to future climate change is a major conservation and societal challenge. Climate change is predicted to force many species to shift their ranges in pursuit of suitable conditions. This study aims to use landscape genetics, the study of the effects of environmental heterogeneity on the spatial distribution of genetic variation, as a predictive tool to assess how species will shift their ranges to track climatic changes and inform conservation measures that will facilitate movement. The approach is based on three steps: 1) using species distribution models (SDMs) to predict suitable ranges under future climate change, 2) using the landscape genetics framework to identify landscape variables that impede or facilitate movement, and 3) extrapolating the effect of landscape connectivity on range shifts in response to future climate change. I show how this approach can be implemented using the publicly available genetic dataset of the grey long-eared bat, Plecotus austriacus, in the Iberian Peninsula. Forest cover gradient was the main landscape variable affecting genetic connectivity between colonies. Forest availability is likely to limit future range shifts in response to climate change, primarily over the central plateau, but important range shift pathways have been identified along the eastern and western coasts. I provide outputs that can be directly used by conservation managers and review the viability of the approach. Using landscape genetics as a predictive tool in combination with SDMs enables the identification of potential pathways, whose loss can affect the ability of species to shift their range into future climatically suitable areas, and the appropriate conservation management measures to increase landscape connectivity and facilitate movement.

Depth as an organizing force in Pocillopora damicornis: intra-reef genetic architecture.

Relative to terrestrial plants, and despite similarities in life history characteristics, the potential for corals to exhibit intra-reef local adaptation in the form of genetic differentiation along an environmental gradient has received little attention. The potential for natural selection to act on such small scales is likely increased by the ability of coral larval dispersal and settlement to be influenced by environmental cues. Here, we combine genetic, spatial, and environmental data for a single patch reef in Kāne‘ohe Bay, O‘ahu, Hawai‘i, USA in a landscape genetics framework to uncover environmental drivers of intra-reef genetic structuring. The genetic dataset consists of near-exhaustive sampling (n = 2352) of the coral, Pocillopora damicornis at our study site and six microsatellite genotypes. In addition, three environmental parameters – depth and two depth-independent temperature indices – were collected on a 4 m grid across 85 locations throughout the reef. We use ordinary kriging to spatially interpolate our environmental data and estimate the three environmental parameters for each colony. Partial Mantel tests indicate a significant correlation between genetic relatedness and depth while controlling for space. These results are also supported by multi-model inference. Furthermore, spatial Principle Component Analysis indicates a statistically significant genetic cline along a depth gradient. Binning the genetic dataset based on size-class revealed that the correlation between genetic relatedness and depth was significant for new recruits and increased for larger size classes, suggesting a possible role of larval habitat selection as well as selective mortality in structuring intra-reef genetic diversity. That both pre- and post-recruitment processes may be involved points to the adaptive role of larval habitat selection in increasing adult survival. The conservation importance of uncovering intra-reef patterns of genetic diversity is discussed.

A practical toolbox for design and analysis of landscape genetics studies

Landscape genetics integrates theory and analytical methods of population genetics and landscape ecology. Research in this area has increased in recent decades, creating a plethora of options for study design and analysis. Here we present a practical toolbox for the design and analysis of landscape genetics studies following a seven-step framework: (1) define the study objectives, (2) consider the spatial and temporal scale of the study, (3) design a sampling regime, (4) select a genetic marker, (5) generate genetic input data, (6) generate spatial input data, and (7) choose an analytical method that integrates genetic and spatial data. Study design considerations discussed include choices of spatial and temporal scale, sample size and spatial distribution, and genetic marker selection. We present analytical methods suitable for achieving different study objectives. As emerging technologies generate genetic and spatial data sets of increasing size, complexity, and resolution, landscape geneticists are challenged to execute hypothesis-driven research that combines empirical data and simulation modeling. The landscape genetics framework presented here can accommodate new design considerations and analyses, and facilitate integration of genetic and spatial data by guiding new landscape geneticists through study design, implementation, and analysis.

RIVERSCAPE GENETICS IDENTIFIES REPLICATED ECOLOGICAL DIVERGENCE ACROSS AN AMAZONIAN ECOTONE

Ecological speciation involves the evolution of reproductive isolation and niche divergence in the absence of a physical barrier to gene flow. The process is one of the most controversial topics of the speciation debate, particularly in tropical regions. Here, we investigate ecologically based divergence across an Amazonian ecotone in the electric fish, Steatogenys elegans. We combine phylogenetics, genome scans, and population genetics with a recently developed individual-based evolutionary landscape genetics approach that incorporates selection. This framework is used to assess the relative contributions of geography and divergent natural selection between environments as biodiversity drivers. We report on two closely related and sympatric lineages that exemplify how divergent selection across a major Amazonian aquatic ecotone (i.e., between rivers with markedly different hydrochemical properties) may result in replicated ecologically mediated speciation. The results link selection across an ecological gradient with reproductive isolation and we propose that assortative mating based on water color may be driving the divergence. Divergence resulting from ecologically driven selection highlights the importance of considering environmental heterogeneity in studies of speciation in tropical regions. Furthermore, we show that framing ecological speciation in a spatially explicit evolutionary landscape genetics framework provides an important first step in exploring a wide range of the potential effects of spatial dependence in natural selection.

Regional Differences in Seasonal Timing of Rainfall Discriminate between Genetically Distinct East African Giraffe Taxa

Masai (Giraffa tippelskirchi), Reticulated (G. reticulata) and Rothschild's (G. camelopardalis) giraffe lineages in East Africa are morphologically and genetically distinct, yet in Kenya their ranges abut. This raises the question of how divergence is maintained among populations of a large mammal capable of long-distance travel, and which readily hybridize in zoos. Here we test four hypotheses concerning the maintenance of the phylogeographic boundaries among the three taxa: 1) isolation-by-distance; 2) physical barriers to dispersal; 3) general habitat differences resulting in habitat segregation; or 4) regional differences in the seasonal timing of rainfall, and resultant timing of browse availability. We used satellite remotely sensed and climate data to characterize the environment at the locations of genotyped giraffes. Canonical variate analysis, random forest algorithms, and generalized dissimilarity modelling were employed in a landscape genetics framework to identify the predictor variables that best explained giraffes' genetic divergence. We found that regional differences in the timing of precipitation, and resulting green-up associated with the abundance of browse, effectively discriminate between taxa. Local habitat conditions, topographic and human-induced barriers, and geographic distance did not aid in discriminating among lineages. Our results suggest that selection associated with regional timing of events in the annual climatic cycle may help maintain genetic and phenotypic divergence in giraffes. We discuss potential mechanisms of maintaining divergence, and suggest that synchronization of reproduction with seasonal rainfall cycles that are geographically distinct may contribute to reproductive isolation. Coordination of weaning with green-up cycles could minimize the costs of lactation and predation on the young. Our findings are consistent with theory and empirical results demonstrating the efficacy of seasonal or phenologically dictated selection pressures in contributing to the reproductive isolation of parapatric populations.

Small geographic range but not panmictic: how forests structure the endangered Point Arena mountain beaver (Aplodontia rufa nigra)

The landscape genetics framework is typically applied to broad regions that occupy only small portions of a species' range. Rarely is the entire range of a taxon the subject of study. We examined the landscape genetic struc- ture of the endangered Point Arena mountain beaver (Apl- odontia rufa nigra), whose isolated geographic range is found in a restricted (85 km 2 ) but heterogenous region in California. Based on its diminutive range we may predict widespread gene flow and a relatively weak role for land- scape variation in defining genetic structure. We used skin, bone, tissue and noninvasively collected hair samples to describe genetic substructure and model gene flow. We examined spatial partitioning of multilocus DNA genotypes and mitochondrial DNA haplotypes. We identified 3 groups from microsatellite data, all of which had low estimates of effective population size consistent with significant tests for historical bottlenecks. We used least-cost-path analysis with the microsatellites to examine how vegetation type affects gene flow in a landscape genetics framework. Gene flow was best predicted by a model with ''Forest'' as the most permeable, followed by ''Riparian''. Agricultural lands demonstrated the highest resistance. MtDNA data revealed only two haplotypes: one was represented in all 57 individ- uals that occurred north of the east-west flowing Garcia River. South of the river, however, both haplotypes occurred, often at the same site suggesting that the river may have affected historical patterns of genetic divergence.

Spatial Genetic Structure of a Symbiotic Beetle-Fungal System: Toward Multi-Taxa Integrated Landscape Genetics

Spatial patterns of genetic variation in interacting species can identify shared features that are important to gene flow and can elucidate co-evolutionary relationships. We assessed concordance in spatial genetic variation between the mountain pine beetle (Dendroctonus ponderosae) and one of its fungal symbionts, Grosmanniaclavigera, in western Canada using neutral genetic markers. We examined how spatial heterogeneity affects genetic variation within beetles and fungi and developed a novel integrated landscape genetics approach to assess reciprocal genetic influences between species using constrained ordination. We also compared landscape genetic models built using Euclidean distances based on allele frequencies to traditional pair-wise Fst. Both beetles and fungi exhibited moderate levels of genetic structure over the total study area, low levels of structure in the south, and more pronounced fungal structure in the north. Beetle genetic variation was associated with geographic location while that of the fungus was not. Pinevolume and climate explained beetle genetic variation in the northern region of recent outbreak expansion. Reciprocal genetic relationships were only detectedin the south where there has been alonger history of beetle infestations. The Euclidean distance and Fst-based analyses resulted in similar models in the north and over the entire study area, but differences between methods in the south suggest that genetic distances measures should be selected based on ecological and evolutionary contexts. The integrated landscape genetics framework we present is powerful, general, and can be applied to other systems to quantify the biotic and abiotic determinants of spatial genetic variation within and among taxa.

Landscape genetics of a specialized grasshopper inhabiting highly fragmented habitats: a role for spatial scale

AbstractAim The study of geographical discontinuities in the distribution of genetic variability in natural populations is a central topic in both evolutionary and conservation research. In this study, we aimed to analyse (1) the factors associated with genetic diversity at the landscape spatial scale in the highly specialized grasshopper Mioscirtus wagneri and (2) to identify the relative contribution of alternative factors to the observed patterns of genetic structure in this species.Location La Mancha region, Central Spain.Methods We sampled 28 populations of the grasshopper M. wagneri and genotyped 648 individuals at seven microsatellite loci. We employed a causal modelling approach to identify the most influential variables associated with genetic differentiation within a multiple hypothesis‐testing framework.Results We found that genetic diversity differs among populations located in different river basins and decreases with population isolation. Causal modelling analyses showed variability in the relative influence of the studied landscape features across different spatial scales. When a highly isolated population is considered, the analyses suggested that geographical distance is the only factor explaining the genetic differentiation between populations. When that population is excluded, the causal modelling analysis revealed that elevation and river basins are also relevant factors contributing to explaining genetic differentiation between the studied populations.Main conclusions These results indicate that the spatial scale considered and the inclusion of outlier populations may have important consequences on the inferred contribution of alternative landscape factors on the patterns of genetic differentiation even when all populations are expected to similarly respond to landscape structure. Thus, a multiscale perspective should also be incorporated into the landscape genetics framework to avoid biased conclusions derived from the spatial scale analysed and/or the geographical distribution of the studied populations.

This month in Evolutionary Applications.

This third issue of Evolutionary Applications features a wide range of topics, including evolutionary causes of invasiveness (Lee and Gelembiuk 2008), human sexual dimorphism (Crespi 2008), estimation of effective population size (Nomura 2008), genetic connectivity and habitat restoration (Raeymaekers et al. 2008), maize adaptation (Mercer et al. 2008), fitness consequences of wild-farm salmon hybridization (Fraser et al. 2008), adaptation to environmental stress (Gilchrist et al. 2008), and Darwinian perspectives on fisheries management (Olsen et al. 2008). Similar to the first two issues in this inaugural year, once again we see that evolutionary biology greatly contributes to our understanding of a wide array of practical issues. Below we highlight just a few of the articles representing the breath of topics covered in this month’s issue.