Putting the ‘landscape’ in landscape genetics

Spatial variation in genomic signatures of local adaptation during the cane toad invasion of Australia

Adaptive evolution can facilitate species' range expansions across environmentally heterogeneous landscapes. However, serial founder effects can limit the efficacy of selection, and the evolution of increased dispersal during range expansions may result in gene flow swamping local adaptation. Here, we study how genetic drift, gene flow and selection interact during the cane toad's (Rhinella marina) invasion across the heterogeneous landscape of Australia. Following its introduction in 1935, the cane toad colonised eastern Australia and established several stable range edges. The ongoing, more rapid range expansion in north-central Australia has occurred concomitant with an evolved increase in dispersal capacity. Using reduced representation genomic data of Australian cane toads from the expansion front and from two areas of their established range, we test the hypothesis that high gene flow constrains local adaptation at the expansion front relative to established areas. Genetic analyses indicate the three study areas are genetically distinct but show similar levels of allelic richness, heterozygosity and inbreeding. Markedly higher gene flow or recency of colonisation at the expansion front have likely hindered local adaptation at the time of sampling, as indicated by reduced slopes of genetic-environment associations (GEAs) estimated using a novel application of geographically weighted regression that accounts for allele surfing; GEA slopes are significantly steeper in established parts of the range. Our work bolsters evidence supporting adaptation of invasive species post-introduction and adds novel evidence for differing strengths of evolutionary forces among geographic areas with different invasion histories.

On the road to losing connectivity: Faecal samples provide genome-wide insights into anthropogenic impacts on two large herbivore species in central India

Humans have impacted most of the planet, and the ensuing fragmentation results in small, isolated habitat patches posing a risk of genetic diversity loss, inbreeding, and genetic load. Understanding how natural and anthropogenic landscape features affect gene flow among habitat patches is critical for maintaining connectivity. Genome-wide data are required to comprehend the impacts of recent fragmentation, which can be challenging when only non-invasive samples are available. Here, we build upon advancements in conservation genomics to address connectivity of two large herbivores, gaur (Bos gaurus) and sambar (Rusa unicolor) in central India. Given their habitat associations, we expected these species to respond similarly to habitat fragmentation. We used faecal-DNA and methylation-based host-DNA enrichment with modified ddRAD protocol to generate genome-wide single-nucleotide polymorphism (SNP) data for 124 gaur and 99 sambar individuals. Our findings reveal that gaur populations in central India are fragmented, displaying high genetic differentiation, with drift significantly affecting small populations like Umred Karhandla Wildlife Sanctuary. Although sambar shows low genetic structure, another small population, Bor Tiger Reserve is genetically differentiated. Our results suggest that although land cover change and roads restrict animal movement, the extent of this impact varies across the two species. We show that different species respond differently to landscape features, even with similar habitat associations. We highlight small and isolated populations requiring urgent conservation intervention. Such multi-species approaches enhance our understanding of cross-species connectivity patterns. We suggest shifting from single-species to multi-species holistic conservation approach in rapidly developing landscapes to better manage co-occurring endangered species.

City divided: Unveiling family ties and genetic structuring of coyotes in Seattle

Linear barriers pose significant challenges for wildlife gene flow, impacting species persistence, adaptation, and evolution. While numerous studies have examined the effects of linear barriers (e.g., fences and roadways) on partitioning urban and non-urban areas, understanding their influence on gene flow within cities remains limited. Here, we investigated the impact of linear barriers on coyote (Canis latrans) population structure in Seattle, Washington, where major barriers (i.e., interstate highways and bodies of water) divide the city into distinct quadrants. Just under 1000 scats were collected to obtain genetic data between January 2021 and December 2022, allowing us to identify 73 individual coyotes. Notably, private allele analysis underscored limited interbreeding among quadrants. When comparing one quadrant to each other, there were up to 16 private alleles within a single quadrant, representing nearly 22% of the population allelic diversity. Our analysis revealed weak isolation by distance, and despite being a highly mobile species, genetic structuring was apparent between quadrants even with extremely short geographic distance between individual coyotes, implying that Interstate 5 and the Ship Canal act as major barriers. This study uses coyotes as a model species for understanding urban gene flow and its consequences in cities, a crucial component for bolstering conservation of rarer species and developing wildlife friendly cities.

The effect of terrain on the fine-scale genetic diversity of sub-Antarctic Collembola: A landscape genetics approach

Biodiversity patterns are shaped by the interplay between geodiversity and organismal characteristics. Superimposing genetic structure onto landscape heterogeneity (i.e., landscape genetics) can help to disentangle their interactions and better understand population dynamics. Previous studies on the sub-Antarctic Prince Edward Islands (located midway between Antarctica and Africa) have highlighted the importance of landscape and climatic barriers in shaping spatial genetic patterns and have drawn attention to the value of these islands as natural laboratories for studying fundamental concepts in biology. Here, we assessed the fine-scale spatial genetic structure of the springtail, Cryptopygus antarcticus travei, which is endemic to Marion Island, in tandem with high-resolution geological data. Using a species-specific suite of microsatellite markers, a fine-scale sampling design incorporating landscape complexity and generalised linear models (GLMs), we examined genetic patterns overlaid onto high-resolution digital surface models and surface geology data across two 1-km sampling transects. The GLMs revealed that genetic patterns across the landscape closely track landscape resistance data in concert with landscape discontinuities and barriers to gene flow identified at a scale of a few metres. These results show that the island's geodiversity plays an important role in shaping biodiversity patterns and intraspecific genetic diversity. This study illustrates that fine-scale genetic patterns in soil arthropods are markedly more structured than anticipated, given that previous studies have reported high levels of genetic diversity and evidence of genetic structing linked to landscape changes for springtail species and considering the homogeneity of the vegetation complexes characteristic of the island at the scale of tens to hundreds of metres. By incorporating fine-scale and high-resolution landscape features into our study, we were able to explain much of the observed spatial genetic patterns. Our study highlights geodiversity as a driver of spatial complexity. More widely, it holds important implications for the conservation and management of the sub-Antarctic islands.

Landscape features influencing gene flow and connectivity of an endangered passerine

Dispersal of individuals and gene flow are crucial aspects to maintain genetic diversity and viability of populations, especially in the case of threatened species. Landscape composition and structure may facilitate or limit individual movement within and among populations. We used a landscape genetics approach to assess the connectivity patterns of the threatened Dupont's lark (Chersophilus duponti subsp. duponti), considering their genetic patterns and the landscape features associated with its gene flow in Spain. We analysed the genetic relatedness based on 11 species-specific polymorphic microsatellites on 416 Dupont's lark individuals sampled across peninsular Spain between 2017 and 2019, covering most of the European distribution of the species. To assess the relationship between the landscape composition and the species gene flow, we estimated genetic distance at the individual level (Dps). Next, we built a set of environmental surfaces from two time periods (years 1990 and 2018), based on factors such as land use and topography, influencing individuals' movement. We then obtained resistance surfaces from an optimization process on landscape variables. Landscape genetics analyses were done for single and composite surface models for each year separately. Our findings from both time points show that scatter or mosaic-structured vegetation composed by low agricultural and tree cover and high presence of sclerophyllous shrubs favoured Dupont's lark dispersal, while dense and continuous tree cover, as well as areas of intensive agriculture, were limiting factors. Our results suggest the importance of steppe habitat patches for the species' establishment and dispersal. In addition, our results provide key information to develop conservation measures, including conserving and restoring steppe habitats as scattered and/or mosaic-structured vegetation that could warrant the connectivity and persistence of Dupont's lark populations.

Unraveling the genetic diversity of Ceiba pubiflora (Malvaceae) in isolated limestone outcrops: Conservation strategies

Seasonally Dry Tropical Forests (SDTFs) located on limestone outcrops are vulnerable to degradation caused by timber logging and limestone extraction for cement production. Some of these forests represent the last remnants of native vegetation cover, functioning as isolated islands. Ceiba pubiflora (Malvaceae) is a tree frequently found on limestone outcrops in the central region of Brazil. This study aimed to evaluate the genetic diversity and identify suitable populations for the establishment of Management Units (MUs) for conservation. Inter-simple sequence repeat markers were employed to assess the genetic diversity in ten populations sampled from the Caatinga, Cerrado, and Atlantic Forest biomes. The species exhibited substantial genetic diversity (HT = 0.345; PLP = 97.89%). Populations SAH, JAN, and MON demonstrated elevated rates of polymorphic loci (> 84.2%) along with notable genetic diversity (He > 0.325). Additionally, these populations were the primary contributors to gene flow. The analysis of molecular variance (AMOVA) indicated that most genetic variation occurs within populations (91.5%) than between them. In the Bayesian analysis, the ten populations were clustered into five groups, revealing the presence of at least three barriers to gene flow in the landscape: 1) the Central Plateau or Paranã River valley; 2) near the Espinhaço mountains or the São Francisco River valley; and 3) around the Mantiqueira mountain range, Chapada dos Veadeiros plateau, and disturbed areas. A positive and statistically significant correlation was observed between genetic (θB) and geographic distances (r = 0.425, p = 0.008). Based on these findings, we propose the establishment of Management Units in Minas Gerais state, encompassing the (1) southern region (MIN population), (2) central region (SAH population), and (3) north region (MON population), as well as in Goiás state, covering the (4) Central Plateau region. These units can significantly contribute to preserving the genetic diversity of these trees and protecting their habitat against ongoing threats.

Genetic consequences of Last Glacial-Holocene changes in snowfall regime in Arnica mallotopus populations: A plant confined to heavy-snow areas of Japan

PREMISE

Snow is an important environmental factor affecting plant distribution. Past changes in snowfall regimes may have controlled the demographies of snow-dependent plants. However, our knowledge of changes in the distribution and demographies of such plants is limited because of the lack of fossil records.

METHODS

Population genetic and landscape genetic analyses were used to investigate the response of population dynamics of Arnica mallotopus (Asteraceae)-a plant confined to heavy-snow areas of Japan-to changes in snowfall regimes from the Last Glacial Period to the Holocene.

RESULTS

The population genetic analysis suggested that the four geographic lineages diverged during the Last Glacial Period. The interaction between reduced snowfall and lower temperatures during this period likely triggered population isolation in separate refugia. Subpopulation differentiation in the northern group was lower than in the southern group. Our ecological niche model predicted that the current distribution was patchy in the southern region; that is, the populations were isolated by topologically flat and climatically unsuitable lowlands. The landscape genetic analysis suggested that areas with little snowfall acted as barriers to the Holocene expansion of species distribution and continued limiting gene flow between local populations.

CONCLUSIONS

These findings indicate that postglacial population responses vary among regions and are controlled by environmental and geographic factors. Thus, changes in snowfall regime played a major role in shaping the distribution and genetic structure of the snow-dependent plant.

Disentangling the contemporary and historical effects of landscape on the population genomic variation of two bird species restricted to the highland forest enclaves of northeastern Brazil

Investigating the impact of landscape features on patterns of genetic variation is crucial to understand spatially dependent evolutionary processes. Here, we assess the population genomic variation of two bird species (Conopophaga cearae and Sclerurus cearensis) through the Caatinga moist forest enclaves in northeastern Brazil. To infer the evolutionary dynamics of bird populations through the Late Quaternary, we used genome-wide polymorphism data obtained from double-digestion restriction-site-associated DNA sequencing (ddRADseq), and integrated population structure analyses, historical demography models, paleodistribution modeling, and landscape genetics analyses. We found the population differentiation among enclaves to be significantly related to the geographic distance and historical resistance across the rugged landscape. The climate changes at the end of the Pleistocene to the Holocene likely triggered synchronic population decline in all enclaves for both species. Our findings revealed that both geographic distance and historical connectivity through highlands are important factors that can explain the current patterns of genetic variation. Our results further suggest that levels of population differentiation and connectivity cannot be explained purely on the basis of contemporary environmental conditions. By combining historical demographic analyses and niche modeling predictions in a historical framework, we provide strong evidence that climate fluctuations of the Quaternary promoted population differentiation and a high degree of temporal synchrony among population size changes in both species.

Individual-based landscape genomics for conservation: An analysis pipeline

Landscape genomics can harness environmental and genetic data to inform conservation decisions by providing essential insights into how landscapes shape biodiversity. The massive increase in genetic data afforded by the genomic era provides exceptional resolution for answering critical conservation genetics questions. The accessibility of genomic data for non-model systems has also enabled a shift away from population-based sampling to individual-based sampling, which now provides accurate and robust estimates of genetic variation that can be used to examine the spatial structure of genomic diversity, population connectivity and the nature of environmental adaptation. Nevertheless, the adoption of individual-based sampling in conservation genetics has been slowed due, in large part, to concerns over how to apply methods developed for population-based sampling to individual-based sampling schemes. Here, we discuss the benefits of individual-based sampling for conservation and describe how landscape genomic methods, paired with individual-based sampling, can answer fundamental conservation questions. We have curated key landscape genomic methods into a user-friendly, open-source workflow, which we provide as a new R package, A Landscape Genomics Analysis Toolkit in R (algatr). The algatr package includes novel added functionality for all of the included methods and extensive vignettes designed with the primary goal of making landscape genomic approaches more accessible and explicitly applicable to conservation biology.

Studies on genetic diversity, gene flow and landscape genetic in Avicennia marina: Spatial PCA, Random Forest, and phylogeography approaches

Mangrove forests grow in coastal areas, lagoons, estuaries, and deltas and form the main vegetation in tidal and saline wetlands. Due to the mankind activities and also changes in climate, these forests face degradations and probably extinction in some areas. Avicennia marina is one of the most distributed mangrove species throughout the world. The populations of A. marina occur in a limited region in southern parts of Iran. Very few genetic and spatial analyses are available on these plants from our country. Therefore, the present study was planned to provide detailed information on Avicennia marina populations with regard to genetic diversity, gene flow versus genetic isolation, effects of spatial variables on connectivity and structuring the genetic content of trees populations and also identifying adaptive genetic regions in respond too spatial variables. We used SCoT molecular markers for genetic analyses and utilized different computational approaches for population genetics and landscapes analyses. The results of present study showed a low to moderate genetic diversity in the studied populations and presence of significant Fst values among them. Genetic fragmentation was also observed within each province studied. A limited gene flow was noticed among neighboring populations within a particular province. One population was almost completely isolated from the gene flow with other populations and had peculiar genetic content.Spatial PCA analysis revealed both significant global and local genetic structuring in the studied populations. Spatial variables like humidity, longitude and altitude were the most important spatial features affecting genetic structure in these populations.

Move it or lose it: Predicted effects of culverts and population density on Mojave desert tortoise (Gopherus agassizii) connectivity

Roadways and railways can reduce wildlife movements across landscapes, negatively impacting population connectivity. Connectivity may be improved by structures that allow safe passage across linear barriers, but connectivity could be adversely influenced by low population densities. The Mojave desert tortoise is threatened by habitat loss, fragmentation, and population declines. The tortoise continues to decline as disturbance increases across the Mojave Desert in the southwestern United States. While underground crossing structures, like hydrological culverts, have begun receiving attention, population density has not been considered in tortoise connectivity. Our work asks a novel question: How do culverts and population density affect connectivity and potentially drive genetic and demographic patterns? To explore the role of culverts and population density, we used agent-based spatially explicit forward-in-time simulations of gene flow. We constructed resistance surfaces with a range of barriers to movement and representative of tortoise habitat with anthropogenic disturbance. We predicted connectivity under variable population densities. Simulations were run for 200 non-overlapping generations (3400 years) with 30 replicates using 20 microsatellite loci. We evaluated population genetic structure and diversity and found that culverts would not entirely negate the effects of linear barriers, but gene flow improved. Our results also indicated that density is important for connectivity. Low densities resulted in declines regardless of the landscape barrier scenario (> 75% population census size, > 97% effective population size). Results from our simulation using current anthropogenic disturbance predicted decreased population connectivity over time. Genetic and demographic effects were detectable within five generations (85 years) following disturbance with estimated losses in effective population size of 69%. The pronounced declines in effective population size indicate this could be a useful monitoring metric. We suggest management strategies that improve connectivity, such as roadside fencing tied to culverts, conservation areas in a connected network, and development restricted to disturbed areas.

Incorporating ecology into gene drive modelling

Gene drive technology, in which fast-spreading engineered drive alleles are introduced into wild populations, represents a promising new tool in the fight against vector-borne diseases, agricultural pests and invasive species. Due to the risks involved, gene drives have so far only been tested in laboratory settings while their population-level behaviour is mainly studied using mathematical and computational models. The spread of a gene drive is a rapid evolutionary process that occurs over timescales similar to many ecological processes. This can potentially generate strong eco-evolutionary feedback that could profoundly affect the dynamics and outcome of a gene drive release. We, therefore, argue for the importance of incorporating ecological features into gene drive models. We describe the key ecological features that could affect gene drive behaviour, such as population structure, life-history, environmental variation and mode of selection. We review previous gene drive modelling efforts and identify areas where further research is needed. As gene drive technology approaches the level of field experimentation, it is crucial to evaluate gene drive dynamics, potential outcomes, and risks realistically by including ecological processes.

Barriers to chimpanzee gene flow at the south-east edge of their distribution

Populations on the edge of a species' distribution may represent an important source of adaptive diversity, yet these populations tend to be more fragmented and are more likely to be geographically isolated. Lack of genetic exchanges between such populations, due to barriers to animal movement, can not only compromise adaptive potential but also lead to the fixation of deleterious alleles. The south-eastern edge of chimpanzee distribution is particularly fragmented, and conflicting hypotheses have been proposed about population connectivity and viability. To address this uncertainty, we generated both mitochondrial and MiSeq-based microsatellite genotypes for 290 individuals ranging across western Tanzania. While shared mitochondrial haplotypes confirmed historical gene flow, our microsatellite analyses revealed two distinct clusters, suggesting two populations currently isolated from one another. However, we found evidence of high levels of gene flow maintained within each of these clusters, one of which covers an 18,000 km2 ecosystem. Landscape genetic analyses confirmed the presence of barriers to gene flow with rivers and bare habitats highly restricting chimpanzee movement. Our study demonstrates how advances in sequencing technologies, combined with the development of landscape genetics approaches, can resolve ambiguities in the genetic history of critical populations and better inform conservation efforts of endangered species.

An evaluation of isolation by distance and isolation by resistance on genetic structure of the Persian squirrel (Sciurus anomalus) in the Zagros forests of Iran

For the conservation of wild species, it is important to understand how landscape change and land management can affect gene flow and movement. Landscape genetic analyses provide a powerful approach to infer effects of various landscape factors on gene flow, thereby informing conservation actions. The Persian squirrel is a keystone species in the woodlands and oak forests of Western Asia, where it has experienced recent habitat loss and fragmentation. We conducted landscape genetic analyses of individuals sampled in the northern Zagros Mountains of Iran (provinces of Kurdistan, Kermanshah, and Ilam), focusing on the evaluation of isolation by distance (IBD) and isolation by resistance (IBR), using 16 microsatellite markers. The roles of geographical distance and landscape features including roads, rivers, developed areas, farming and agriculture, forests, lakes, plantation forests, rangelands, shrublands, and rocky areas of varying canopy cover, and swamp margins on genetic structure were quantified using individual-based approaches and resistance surface modeling. We found a significant pattern of IBD but only weak support for an effect of forest cover on genetic structure and gene flow. It seems that geographical distance is an important factor limiting the dispersal of the Persian squirrel in this region. The results of the current study inform ongoing conservation programs for the Persian squirrel in the Zagros oak forest.

Landscape Genetics for Brazilian Equines

Optimization of DNA collection for National gene bank and conservation programs requires information on spatial and genetic distribution of animals countrywide. The relationship between genetic and geographic distances were examined in 8 Brazilian horse breeds (Baixadeiro, Crioulo, Campeiro, Lavradeiro, Marajoara, Mangalarga Marchador, Pantaneiro and Puruca) using Single Nucleotide Polymorphism markers and collection point locations. Mantel correlations, Genetic Landscape Shape Interpolation, Allelic Aggregation Index Analyses and Spatial autocorrelation tests indicated a nonrandom distribution of horses throughout the country. Minimum collection distances for the national Gene Bank should be 530km, with clear divisions seen in genetic structure of horse populations in both North/South and East/West directions. Comparing Pantaneiro and North/Northeastern breeds, physical distance is not necessarily the defining factor for genetic differentiation. This should be considered when sampling these local breeds. These data can help optimise GenBank collection routines and conservation strategies for these breeds.

Spatial genetic structure of 2009 H1N1 pandemic influenza established as a result of interaction with human populations in mainland China

Identifying the spatial patterns of genetic structure of influenza A viruses is a key factor for understanding their spread and evolutionary dynamics. In this study, we used phylogenetic and Bayesian clustering analyses of genetic sequences of the A/H1N1pdm09 virus with district-level locations in mainland China to investigate the spatial genetic structure of the A/H1N1pdm09 virus across human population landscapes. Positive correlation between geographic and genetic distances indicates high degrees of genetic similarity among viruses within small geographic regions but broad-scale genetic differentiation, implying that local viral circulation was a more important driver in the formation of the spatial genetic structure of the A/H1N1pdm09 virus than even, countrywide viral mixing and gene flow. Geographic heterogeneity in the distribution of genetic subpopulations of A/H1N1pdm09 virus in mainland China indicates both local to local transmission as well as broad-range viral migration. This combination of both local and global structure suggests that both small-scale and large-scale population circulation in China is responsible for viral genetic structure. Our study provides implications for understanding the evolution and spread of A/H1N1pdm09 virus across the population landscape of mainland China, which can inform disease control strategies for future pandemics.

The Effect of Landscape Environmental Factors on Gene Flow of Red Deer (Cervus canadensis xanthopygus) in the Southern of the Greater Khingan Mountains, China

Red deer (Cervus canadensis xanthopygus) living in the north of China are restricted and threatened due to human activities and the changes in the natural environment, which influence the dispersal and effective gene flow between different groups of red deer. Effective gene flow plays an important role in maintaining genetic diversity and structure and ensuring population health. In order to evaluate the genetic diversity level and understand the gene flow between different red deer groups, 231 fresh fecal samples were collected from the southern part of the Greater Khingan Mountains, China. A microsatellite marker was used for genetic analysis. The results showed that the genetic diversity of red deer was intermediate in this region. Significant genetic differentiation among different groups was found in the main distribution area (p < 0.01) using F-statistics and the program STRUCTURE. Different degrees of gene flow existed in red deer groups, and the roads (importance = 40.9), elevation (importance = 38.6), and settlements (importance = 14.1) exerted main effects on gene flow between red deer groups. Human-made factors should be noticed and strictly supervised in this region to avoid excessive disturbance to the normal movement of the red deer. Further conservation and management of red deer should reduce the intensity of vehicular traffic in the concentrated distribution areas of red deer, especially during the heat season. This research helps us better understand the genetic level and health status of red deer in the southern part of the Greater Khingan Mountains and provides theoretical references for protecting and restoring the red deer populations in China.

Environmental selection, rather than neutral processes, best explain regional patterns of diversity in a tropical rainforest fish

To conserve the high functional and genetic variation in hotspots such as tropical rainforests, it is essential to understand the forces driving and maintaining biodiversity. We asked to what extent environmental gradients and terrain structure affect morphological and genomic variation across the wet tropical distribution of an Australian rainbowfish, Melanotaenia splendida splendida. We used an integrative riverscape genomics and morphometrics framework to assess the influence of these factors on both putative adaptive and non-adaptive spatial divergence. We found that neutral genetic population structure was largely explainable by restricted gene flow among drainages. However, environmental associations revealed that ecological variables had a similar power to explain overall genetic variation, and greater power to explain body shape variation, than the included neutral covariables. Hydrological and thermal variables were the strongest environmental predictors and were correlated with traits previously linked to heritable habitat-associated dimorphism in rainbowfishes. In addition, climate-associated genetic variation was significantly associated with morphology, supporting heritability of shape variation. These results support the inference of evolved functional differences among localities, and the importance of hydroclimate in early stages of diversification. We expect that substantial evolutionary responses will be required in tropical rainforest endemics to mitigate local fitness losses due to changing climates.

Landscape genetics identifies barriers to Natterjack toad metapopulation dispersal

Habitat fragmentation and loss reduce population size and connectivity, which imperils populations. Functional connectivity is key for species persistence in human-modified landscapes. To inform species conservation management, we investigated spatial genetic structure, gene flow and inferred dispersal between twelve breeding sites of the Natterjack toad (Bufo calamita); regionally Red-Listed as Endangered in Ireland. Spatial genetic structure was determined using both Bayesian and non-Bayesian clustering analysis of 13 polymorphic microsatellite loci genotyping 247 individuals. We tested the influence of geographic distance, climate, habitat, geographical features, and anthropogenic pressure on pairwise genetic distances between breeding sites using Isolation-by-distance and Isolation-by-resistance based on least-cost path and circuit theory models of functional connectivity. There was clear spatial structuring with genetic distances increasing with geographic distance. Gene flow was best explained by Isolation-by-resistance models with coniferous forestry plantations, bog, marsh, moor and heath, scrub, anthropogenic presence (Human Influence Index) and rivers (riparian density) identified as habitats with high resistance to gene flow while metapopulation connectivity was enhanced by coastal habitats (beaches, sand dunes and salt marshes) and coastal grassland. Despite substantial declines in census numbers over the past 15 years and its regional status as Endangered, the Natterjack toad population in Ireland retains high genetic diversity. If declines continue, maintaining habitat connectivity to prevent genetic erosion by management of coastal grasslands, pond construction and assisted migration through translocation will be increasingly important.

Whole genome demographic models indicate divergent effective population size histories shape contemporary genetic diversity gradients in a montane bumble bee

Understanding historical range shifts and population size variation provides an important context for interpreting contemporary genetic diversity. Methods to predict changes in species distributions and model changes in effective population size (N _e) using whole genomes make it feasible to examine how temporal dynamics influence diversity across populations. We investigate N _e variation and climate-associated range shifts to examine the origins of a previously observed latitudinal heterozygosity gradient in the bumble bee Bombus vancouverensis Cresson (Hymenoptera: Apidae: Bombus Latreille) in western North America. We analyze whole genomes from a latitude-elevation cline using sequentially Markovian coalescent models of N _e through time to test whether relatively low diversity in southern high-elevation populations is a result of long-term differences in N _e. We use Maxent models of the species range over the last 130,000 years to evaluate range shifts and stability. N _e fluctuates with climate across populations, but more genetically diverse northern populations have maintained greater N _e over the late Pleistocene and experienced larger expansions with climatically favorable time periods. Northern populations also experienced larger bottlenecks during the last glacial period, which matched the loss of range area near these sites; however, bottlenecks were not sufficient to erode diversity maintained during periods of large N _e. A genome sampled from an island population indicated a severe postglacial bottleneck, indicating that large recent postglacial declines are detectable if they have occurred. Genetic diversity was not related to niche stability or glacial-period bottleneck size. Instead, spatial expansions and increased connectivity during favorable climates likely maintain diversity in the north while restriction to high elevations maintains relatively low diversity despite greater stability in southern regions. Results suggest genetic diversity gradients reflect long-term differences in N _e dynamics and also emphasize the unique effects of isolation on insular habitats for bumble bees. Patterns are discussed in the context of conservation under climate change.

Biogeographic inferences across spatial and evolutionary scales

The field of biogeography unites landscape genetics and phylogeography under a common conceptual framework. Landscape genetics traditionally focuses on recent-time, population-based, spatial genetics processes at small geographical scales, while phylogeography typically investigates deep past, lineage- and species-based processes at large geographical scales. Here, we evaluate the link between landscape genetics and phylogeographical methods using the western fence lizard (Sceloporus occidentalis) as a model species. First, we conducted replicated landscape genetics studies across several geographical scales to investigate how population genetics inferences change depending on the spatial extent of the study area. Then, we carried out a phylogeographical study of population structure at two evolutionary scales informed by inferences derived from landscape genetics results to identify concordance and conflict between these sets of methods. We found significant concordance in landscape genetics processes at all but the largest geographical scale. Phylogeographical results indicate major clades are restricted to distinct river drainages or distinct hydrological regions. At a more recent timescale, we find minor clades are restricted to single river canyons in the majority of cases, while the remainder of river canyons include samples from at most two clades. Overall, the broad-scale pattern implicating stream and river valleys as key features linking populations in the landscape genetics results, and high degree of clade specificity within major topographic subdivisions in the phylogeographical results, is consistent. As landscape genetics and phylogeography share many of the same objectives, synthesizing theory, models and methods between these fields will help bring about a better understanding of ecological and evolutionary processes structuring genetic variation across space and time.

The genetic structure and connectivity in two sympatric rodent species with different life histories are similarly affected by land use disturbances

The negative impact of habitat fragmentation due to human activities may be different in different species that co-exist in the same area, with consequences on the development of environmental protection plans. Here we aim at understanding the effects produced by different natural and anthropic landscape features on gene flow patterns in two sympatric species with different specializations, one generalist and one specialist, sampled in the same locations. We collected and genotyped 194 wood mice (generalist species) and 199 bank voles (specialist species) from 15 woodlands in a fragmented landscape characterized by different potential barriers to dispersal. Genetic variation and structure were analyzed in the two species, respectively. Effective migration surfaces, isolation-by-resistance (IBR) analysis, and regression with randomization were used to investigate isolation-by-distance (IBD) and the relative importance of land cover elements on gene flow. We observed similar patterns of heterozygosity and IBD for both species, but the bank vole showed higher genetic differences among geographic areas. The IBR analysis suggests that (i) connectivity is reduced in both species by urban areas but more strongly in the specialist bank vole; (ii) cultivated areas act as dispersal corridors in both species; (iii) woodlands appear to be an important factor in increasing connectivity in the bank vole, and less so in the wood mouse. The difference in dispersal abilities between a generalist and specialist species was reflected in the difference in genetic structure, despite extensive habitat changes due to human activities. The negative effects of fragmentation due to the process of urbanization were, at least partially, mitigated by another human product, i.e., cultivated terrains subdivided by hedgerows, and this was true for both species.

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.

Forest cover and geographical distance influence fine-scale genetic structure of leaf-toed geckos in the tropical dry forests of western Mexico

The biodiversity within tropical dry forests (TDFs) is astounding and yet poorly catalogued due to inadequate sampling and the presence of cryptic species. In the Mexican TDF, endemic species are common, and the landscape has been continually altered by geological and anthropogenic changes. To understand how landscape and environmental variables have shaped the population structure of endemic species, we studied the recently described species of leaf-toed gecko, Phyllodactylus benedettii, in coastal western Mexico. Using double-digest restriction site-associated DNA sequencing data, we first explore population structure and estimate the number of ancestral populations. The results indicate a high degree of genetic structure with little admixture, and patterns corresponding to both latitudinal and altitudinal gradients. We find that genetic structure cannot be explained purely by geographical distance, and that ecological corridors may facilitate dispersal and gene flow. We then model the spatial distribution of P. benedettii in the TDF through time and find that the coastline has been climatically suitable for the species since the Last Glacial Maximum. Landscape genetic analyses suggest that the combination of isolation by distance (IBD) and isolation by resistance (IBR; forest cover) has influenced the spatial genetic structure of the species. Overall, our genomic data demonstrate fine-scale population structure in TDF habitat, a complex colonization history, and spatial patterns consistent with both IBD and other ecological factors. These results further highlight the Mexican TDF as a diversity hotspot and suggest that continued anthropogenic changes are likely to affect native fauna.

Integrating data types to estimate spatial patterns of avian migration across the Western Hemisphere

For many avian species, spatial migration patterns remain largely undescribed, especially across hemispheric extents. Recent advancements in tracking technologies and high-resolution species distribution models (i.e., eBird Status and Trends products) provide new insights into migratory bird movements and offer a promising opportunity for integrating independent data sources to describe avian migration. Here, we present a three-stage modeling framework for estimating spatial patterns of avian migration. First, we integrate tracking and band re-encounter data to quantify migratory connectivity, defined as the relative proportions of individuals migrating between breeding and nonbreeding regions. Next, we use estimated connectivity proportions along with eBird occurrence probabilities to produce probabilistic least-cost path (LCP) indices. In a final step, we use generalized additive mixed models (GAMMs) both to evaluate the ability of LCP indices to accurately predict (i.e., as a covariate) observed locations derived from tracking and band re-encounter data sets versus pseudo-absence locations during migratory periods and to create a fully integrated (i.e., eBird occurrence, LCP, and tracking/band re-encounter data) spatial prediction index for mapping species-specific seasonal migrations. To illustrate this approach, we apply this framework to describe seasonal migrations of 12 bird species across the Western Hemisphere during pre- and postbreeding migratory periods (i.e., spring and fall, respectively). We found that including LCP indices with eBird occurrence in GAMMs generally improved the ability to accurately predict observed migratory locations compared to models with eBird occurrence alone. Using three performance metrics, the eBird + LCP model demonstrated equivalent or superior fit relative to the eBird-only model for 22 of 24 species-season GAMMs. In particular, the integrated index filled in spatial gaps for species with over-water movements and those that migrated over land where there were few eBird sightings and, thus, low predictive ability of eBird occurrence probabilities (e.g., Amazonian rainforest in South America). This methodology of combining individual-based seasonal movement data with temporally dynamic species distribution models provides a comprehensive approach to integrating multiple data types to describe broad-scale spatial patterns of animal movement. Further development and customization of this approach will continue to advance knowledge about the full annual cycle and conservation of migratory birds.

Main roads and land cover shaped the genetic structure of a Mediterranean island wild boar population

Patterns of genetic differentiation within and among animal populations might vary due to the simple effect of distance or landscape features hindering gene flow. An assessment of how landscape connectivity affects gene flow can help guide management, especially in fragmented landscapes. Our objective was to analyze population genetic structure and landscape genetics of the native wild boar (Sus scrofa meridionalis) population inhabiting the island of Sardinia (Italy), and test for the existence of Isolation-by-Distance (IBD), Isolation-by-Barrier (IBB), and Isolation-by-Resistance (IBR). A total of 393 Sardinian wild boar samples were analyzed using a set of 16 microsatellite loci. Signals of genetic introgression from introduced non-native wild boars or from domestic pigs were revealed by a Bayesian cluster analysis including 250 reference individuals belonging to European wild populations and domestic breeds. After removal of introgressed individuals, genetic structure in the population was investigated by different statistical approaches, supporting a partition into five discrete subpopulations, corresponding to five geographic areas on the island: north-west (NW), central west (CW), south-west (SW), north-central east (NCE), and south-east (SE). To test the IBD, IBB, and IBR hypotheses, we optimized resistance surfaces using genetic algorithms and linear mixed-effects models with a maximum likelihood population effects parameterization. Landscape genetics analyses revealed that genetic discontinuities between subpopulations can be explained by landscape elements, suggesting that main roads, urban settings, and intensively cultivated areas are hampering gene flow (and thus individual movements) within the Sardinian wild boar population. Our results reveal how human-transformed landscapes can affect genetic connectivity even in a large-sized and highly mobile mammal such as the wild boar, and provide crucial information to manage the spread of pathogens, including the African Swine Fever virus, endemic in Sardinia.

Riverscape Genomics Clarifies Neutral and Adaptive Evolution in an Amazonian Characin Fish (Triportheus albus)

Understanding the role of natural selection in the evolution of wild populations is challenging due to the spatial complexity of natural systems. The richest diversity of freshwater fishes in the world is found in the Amazon Basin, a system where marked hydrochemical differences exist at the interface of major rivers with distinct “water colors” (i.e., black, white, and clear water). We hypothesize that divergent natural selection associated with these “aquatic ecotones” influences population-level adaptive divergence in the non-migratory Amazonian fish fauna. This hypothesis was tested using a landscape genomics framework to compare the relative contribution of environmental and spatial factors to the evolutionary divergence of the Amazonian characin fish Triportheus albus. The framework was based on spatial data, in situ hydrochemical measurements, and 15,251 filtered SNPs (single nucleotide polymorphisms) for T. albus sampled from three major Amazonian rivers. Gradient Forest, redundancy analysis (RDA) and BayPass analyses were used to test for signals of natural selection, and model-based and model-free approaches were used to evaluate neutral population differentiation. After controlling for a signal of neutral hierarchical structure which was consistent with the expectations for a dendritic system, variation in turbidity and pH were key factors contributing to adaptive divergence. Variation in genes involved in acid-sensitive ion transport pathways and light-sensitive photoreceptor pathways was strongly associated with pH and turbidity variability. This study improves our understanding of how natural selection and neutral evolution impact on the distribution of aquatic biodiversity from the understudied and ecologically complex Amazonia.

Environmental differences explain subtle yet detectable genetic structure in a widespread pollinator

BACKGROUND

The environment is a strong driver of genetic structure in many natural populations, yet often neglected in population genetic studies. This may be a particular problem in vagile species, where subtle structure cannot be explained by limitations to dispersal. Consequently, these species might falsely be considered quasi-panmictic and hence potentially mismanaged. A species this might apply to, is the buff-tailed bumble bee (Bombus terrestris), an economically important and widespread pollinator, which is considered to be quasi-panmictic at mainland continental scales. Here we aimed to (i) quantify genetic structure in 21+ populations of the buff-tailed bumble bee, sampled throughout two Eastern European countries, and (ii) analyse the degree to which structure is explained by environmental differences, habitat permeability and geographic distance. Using 12 microsatellite loci, we characterised populations of this species with Fst analyses, complemented by discriminant analysis of principal components and Bayesian clustering approaches. We then applied generalized dissimilarity modelling to simultaneously assess the informativeness of geographic distance, habitat permeability and environmental differences among populations in explaining divergence.

RESULTS

Genetic structure of the buff-tailed bumble bee quantified by means of Fst was subtle and not detected by Bayesian clustering. Discriminant analysis of principal components suggested insignificant but still noticeable structure that slightly exceeded estimates obtained through Fst analyses. As expected, geographic distance and habitat permeability were not informative in explaining the spatial pattern of genetic divergence. Yet, environmental variables related to temperature, vegetation and topography were highly informative, explaining between 33 and 39% of the genetic variation observed.

CONCLUSIONS

In contrast to previous studies reporting quasi-panmixia in continental populations of this species, we demonstrated the presence of subtle population structure related to environmental heterogeneity. Environmental data proved to be highly useful in unravelling the drivers of genetic structure in this vagile and opportunistic species. We highlight the potential of including these data to obtain a better understanding of population structure and the processes driving it in species considered to be quasi-panmictic.

Host, microbiome, and complex space: applying population and landscape genetic approaches to gut microbiome research in wild populations

In recent years, emerging sequencing technologies and computational tools have driven a tidal wave of research on host-associated microbiomes, particularly the gut microbiome. These studies demonstrate numerous connections between the gut microbiome and vital host functions, primarily in humans, model organisms, and domestic animals. As the adaptive importance of the gut microbiome becomes clearer, interest in studying the gut microbiomes of wild populations has increased, in part due to the potential for discovering conservation applications. The study of wildlife gut microbiomes holds many new challenges and opportunities due to the complex genetic, spatial, and environmental structure of wild host populations, and the potential for these factors to interact with the microbiome. The emerging picture of adaptive coevolution in host-microbiome relationships highlights the importance of understanding microbiome variation in the context of host population genetics and landscape heterogeneity across a wide range of host populations. We propose a conceptual framework for understanding wildlife gut microbiomes in relation to landscape variables and host population genetics, including the potential of approaches derived from landscape genetics. We use this framework to review current research, synthesize important trends, highlight implications for conservation, and recommend future directions for research. Specifically, we focus on how spatial structure and environmental variation interact with host population genetics and microbiome variation in natural populations, and what we can learn from how these patterns of covariation differ depending on host ecological and evolutionary traits.

Seascape Genetics of the Atlantic Spotted Dolphin (Stenella frontalis) Based on Mitochondrial DNA

The Atlantic spotted dolphin (Stenella frontalis) is endemic to tropical, subtropical, and warm temperate waters of the Atlantic Ocean. Throughout its distribution, both geographic distance and environmental variation may contribute to population structure of the species. In this study, we follow a seascape genetics approach to investigate population differentiation of Atlantic spotted dolphins based on a large worldwide dataset and the relationship with marine environmental variables. The results revealed that the Atlantic spotted dolphin exhibits population genetic structure across its distribution based on mitochondrial DNA control region (mtDNA-CR) data. Analyses based on the contemporary landscape suggested, at both the individual and population level, that the population genetic structure is consistent with the isolation-by-distance model. However, because geography and environmental matrices were correlated, and because in some, but not all analyses, we found a significant effect for the environment, we cannot rule out the addition contribution of environmental factors in structuring genetic variation. Future analyses based on nuclear data are needed to evaluate whether local processes, such as social structure and some level of philopatry within populations, may be contributing to the associations among genetic structure, geographic, and environmental distance.

Environmentally associated chromosomal structural variation influences fine-scale population structure of Atlantic Salmon (Salmo salar)

Chromosomal rearrangements (e.g., inversions, fusions, and translocations) have long been associated with environmental variation in wild populations. New genomic tools provide the opportunity to examine the role of these structural variants in shaping adaptive differences within and among wild populations of non-model organisms. In Atlantic Salmon (Salmo salar), variations in chromosomal rearrangements exist across the species natural range, yet the role and importance of these structural variants in maintaining adaptive differences among wild populations remains poorly understood. We genotyped Atlantic Salmon (n = 1429) from 26 populations within a highly genetically structured region of southern Newfoundland, Canada with a 220K SNP array. Multivariate analysis, across two independent years, consistently identified variation in a structural variant (translocation between chromosomes Ssa01 and Ssa23), previously associated with evidence of trans-Atlantic secondary contact, as the dominant factor influencing population structure in the region. Redundancy analysis suggested that variation in the Ssa01/Ssa23 chromosomal translocation is strongly correlated with temperature. Our analyses suggest environmentally mediated selection acting on standing genetic variation in genomic architecture introduced through secondary contact may underpin fine-scale local adaptation in Placentia Bay, Newfoundland, Canada, a large and deep embayment, highlighting the importance of chromosomal structural variation as a driver of contemporary adaptive divergence.

Population genetic structure of the Natterjack toad (Epidalea calamita) in Ireland: implications for conservation management

Molecular methods can play a crucial role in species management and conservation. Despite the usefulness of genetic approaches, they are often not explicitly included as part of species recovery plans and conservation practises. The Natterjack toad (Epidalea calamita) is regionally Red-Listed as Endangered in Ireland. The species is declining and is now present at just seven sites within a highly restricted range. This study used 13 highly polymorphic microsatellite markers to analyse the population genetic diversity and structure. Genetic diversity was high with expected heterozygosity between 0.55 and 0.61 and allelic richness between 4.77 and 5.92. Effective population sizes were small (Ne < 100 individuals), but not abnormal for pond breeding amphibians. However, there was no evidence of historical or contemporary genetic bottlenecks or high levels of inbreeding. We identified a positive relationship between Ne and breeding pond surface area, suggesting that environmental factors are a key determinant of population size. Significant genetic structuring was detected throughout the species’ range, and we identified four genetic entities that should be considered in the species’ conservation strategies. Management should focus on preventing further population declines and future loss of genetic diversity overall and within genetic entities while maintaining adequate local effective population size through site-specific protection, human-mediated translocations and head-start programs. The apparent high levels of genetic variation give hope for the conservation of Ireland’s rarest amphibian if appropriately protected and managed.

Agricultural intensification alters marbled newt genetic diversity and gene flow through density and dispersal reduction

Recent agricultural intensification threatens global biodiversity with amphibians being one of the most impacted groups. Because of their biphasic life cycle, amphibians are particularly vulnerable to habitat loss and fragmentation that often result in small, isolated populations and loss of genetic diversity. Here, we studied how landscape heterogeneity affects genetic diversity, gene flow and demographic parameters in the marbled newt, Triturus marmoratus, over a hedgerow network landscape in Western France. While the northern part of the study area consists of preserved hedged farmland, the southern part was more profoundly converted for intensive arable crops production after WWII. Based on 67 sampled ponds and 10 microsatellite loci, we characterized regional population genetic structure and evaluated the correlation between landscape variables and (i) local genetic diversity using mixed models and (ii) genetic distance using multiple regression methods and commonality analysis. We identified a single genetic population characterized by a spatially heterogeneous isolation-by-distance pattern. Pond density in the surrounding landscape positively affected local genetic diversity while arable crop land cover negatively affected gene flow and connectivity. We used demographic inferences to quantitatively assess differences in effective population density and dispersal between the contrasted landscapes characterizing the northern and southern parts of the study area. Altogether, results suggest recent land conversion affected T. marmoratus through reduction in both effective population density and dispersal due to habitat loss and reduced connectivity.

EST-SSR-based landscape genetics of Pseudotaxus chienii, a tertiary relict conifer endemic to China

Pseudotaxus chienii, belonging to the monotypic genus Pseudotaxus (Taxaceae), is a relict conifer endemic to China. Its populations are usually small and patchily distributed, having a low capacity of natural regeneration. To gain a clearer understanding of how landscape variables affect the local adaptation of P. chienii, we applied EST-SSR markers in conjunction with landscape genetics methods: (a) to examine the population genetic pattern and spatial genetic structure; (b) to perform genome scan and selection scan to identify outlier loci and the associated landscape variables; and (c) to model the ecological niche under climate change. As a result, P. chienii was found to have a moderate level of genetic variation and a high level of genetic differentiation. Its populations displayed a significant positive relationship between the genetic and geographical distance (i.e., "isolation by distance" pattern) and a strong fine-scale spatial genetic structure within 2 km. A putatively adaptive locus EMS6 (functionally annotated to cellulose synthase A catalytic subunit 7) was identified, which was found significantly associated with soil Cu, K, and Pb content and the combined effects of temperature and precipitation. Moreover, P. chienii was predicted to experience significant range contractions in future climate change scenarios. Our results highlight the potential of specific soil metal content and climate variables as the driving force of adaptive genetic differentiation in P. chienii. The data would also be useful to develop a conservation action plan for P. chienii.

Can We Compare Effect Size of Spatial Genetic Structure Between Studies and Species Using Moran Eigenvector Maps?

As the field of landscape genetics is progressing toward comparative empirical studies and meta-analysis, it is important to know how best to compare the strength of spatial genetic structure between studies and species. Moran’s Eigenvector Maps are a promising method that does not make an assumption of isolation-by-distance in a homogeneous environment but can discern cryptic structure that may result from multiple processes operating in heterogeneous landscapes. MEMgene uses spatial filters from Moran’s Eigenvector Maps as predictor variables to explain variation in a genetic distance matrix, and it returns adjusted R2 as a measure of the amount of genetic variation that is spatially structured. However, it is unclear whether, and under which conditions, this value can be used to compare the degree of spatial genetic structure (effect size) between studies. This study addresses the fundamental question of comparability at two levels: between independent studies (meta-analysis mode) and between species sampled at the same locations (comparative mode). We used published datasets containing 9,900 haploid, biallelic, neutral loci simulated on a quasi-continuous, square landscape under four demographic scenarios (island model, isolation-by-distance, expansion from one or two refugia). We varied the genetic resolution (number of individuals and loci) and the number of random sampling locations. We considered two measures of effect size, the MEMgene adjusted R2 and multivariate Moran’s I, which is related to Moran’s Eigenvector Maps. Both metrics were highly sensitive to the number of locations, even when using standardized effect sizes, SES, and the number of individuals sampled per location, but not to the number of loci. In comparative mode, using the same Moran Eigenvector Maps for all species, even those with missing values at some sampling locations, reduced bias due to the number of locations under isolation-by-distance (stationary process) but increased it under expansion from one or two refugia (non-stationary process). More robust measures of effect size need to be developed before the strength of spatial genetic structure can be accurately compared, either in a meta-analysis of independent empirical studies or within a comparative, multispecies landscape genetic study.

Primate landscape genetics: A review and practical guide

Landscape genetics is an emerging field that integrates population genetics, landscape ecology, and spatial statistics to investigate how geographical and environmental features and evolutionary processes such as gene flow, genetic drift, and selection structure genetic variation at both the population and individual levels, with implications for ecology, evolution, and conservation biology. Despite being particularly well suited for primatologists, this method is currently underutilized. Here, we synthesize the current state of research on landscape genetics in primates. We begin by outlining how landscape genetics has been used to disentangle the drivers of diversity, followed by a review of how landscape genetic methods have been applied to primates. This is followed by a section highlighting special considerations when applying the methods to primates, and a practical guide to facilitate further landscape genetics studies using both existing and de novo datasets. We conclude by exploring future avenues of inquiry that could be facilitated by recent developments as well as underdeveloped applications of landscape genetics to primates.

Connectivity Predicts Presence but Not Population Density in the Habitat-Specific Mountain Lizard Iberolacerta martinezricai

The Batuecan lizard Iberolacerta martinezricai is a critically endangered species due to its significantly reduced distribution, which is restricted to the scree slopes (SS) of a few mountain peaks within the Batuecas-Sierra de Francia Natural Park (western Spain). Given its high specialisation in this type of discontinuous habitat, the long-term conservation of the species requires maintaining the connectivity between populations. This study analyses the contribution of connectivity, as well as other patch-related factors, in the distribution and density patterns of the species. With this aim, 67 SS were sampled by line transects from May to October 2018. Each SS was characterised using variables indicative of the microhabitat conditions for the lizard. Inter-SS connectivity was quantified using graph theory for seven distances. Generalised linear models (GLMs) were performed for both presence and density. Model results showed that while connectivity was a relevant factor in the presence of lizards, density only involved patch-related variables. Discrepancies probably occurred because the factors influencing presence operate on a wider scale than those of abundance. In view of the results, the best-connected SS, but also those where the lizard is most abundant and from which more dispersed individuals are likely to depart, seem to be the essential patches in any conservation strategy. The results may also be relevant to other species with habitat-specific requirements.

Adjacent Terrestrial Landscapes Impact the Biogeographical Pattern of Soil Escherichia coli Strains in Produce Fields by Modifying the Importance of Environmental Selection and Dispersal

High-quality habitats for wildlife (e.g., forest) provide essential ecosystem services while increasing species diversity and habitat connectivity. Unfortunately, the presence of such habitats adjacent to produce fields may increase risk for contamination of fruits and vegetables by enteric bacteria, including Escherichia coliE. coli survives in extrahost environments (e.g., soil) and could be dispersed across landscapes by wildlife. Understanding how terrestrial landscapes impact the distribution of soil E. coli strains is of importance in assessing the contamination risk of agricultural products. Here, using multilocus sequence typing, we characterized 938 E. coli soil isolates collected from two watersheds with different landscape patterns in New York State, USA, and compared the distribution of E. coli and the influence that environmental selection and dispersal have on the distribution between these two watersheds. Results showed that for the watershed with widespread produce fields, sparse forests, and limited interaction between the two land use types, E. coli composition was significantly different between produce field sites and forest sites; this distribution appears to be shaped by relatively strong environmental selection, likely from soil phosphorus, and slight dispersal limitation. For the watershed with more forested areas and stronger interaction between produce field sites and forest sites, E. coli composition between these two land use types was relatively homogeneous; this distribution appeared to be a consequence of wildlife-driven dispersal, inferred by competing models. Collectively, our results suggest that terrestrial landscape attributes could impact the biogeographic pattern of enteric bacteria by adjusting the importance of environmental selection and dispersal.IMPORTANCE Understanding the ecology of enteric bacteria in extrahost environments is important for the development and implementation of strategies to minimize preharvest contamination of produce with enteric pathogens. Our findings suggest that watershed landscape is an important factor influencing the importance of ecological drivers and dispersal patterns of E. coli Agricultural areas in such watersheds may have a higher risk of produce contamination due to fewer environmental constraints and higher potential of dispersal of enteric bacteria between locations. Thus, there is a perceived trade-off between priorities of environmental conservation and public health in on-farm food safety, with limited ecological data supporting or refuting the role of wildlife in dispersing pathogens under normal operating conditions. By combining field sampling and spatial modeling, we explored ecological principles underlying the biogeographic pattern of enteric bacteria at the regional level, which can benefit agricultural, environmental, and public health scientists who aim to reduce the risk of food contamination by enteric bacteria while minimizing negative impacts on wildlife habitats.

Analysing landscape effects on dispersal networks and gene flow with genetic graphs

Graph-theoretic approaches have relevant applications in landscape genetic analyses. When species form populations in discrete habitat patches, genetic graphs can be used (a) to identify direct dispersal paths followed by propagules or (b) to quantify landscape effects on multi-generational gene flow. However, the influence of their construction parameters remains to be explored. Using a simulation approach, we constructed genetic graphs using several pruning methods (geographical distance thresholds, topological constraints, statistical inference) and genetic distances to weight graph links (F_ST , D_PS , Euclidean genetic distances). We then compared the capacity of these different graphs to (a) identify the precise topology of the dispersal network and (b) to infer landscape resistance to gene flow from the relationship between cost-distances and genetic distances. Although not always clear-cut, our results showed that methods based on geographical distance thresholds seem to better identify dispersal networks in most cases. More interestingly, our study demonstrates that a sub-selection of pairwise distances through graph pruning (thereby reducing the number of data points) can counter-intuitively lead to improved inferences of landscape effects on dispersal. Finally, we showed that genetic distances such as the D_PS or Euclidean genetic distances should be preferred over the F_ST for landscape effect inference as they respond faster to landscape changes.

Conservation Genomics of the Threatened Western Spadefoot, Spea hammondii, in Urbanized Southern California

Populations of the western spadefoot (Spea hammondii) in southern California occur in one of the most urbanized and fragmented landscapes on the planet and have lost up to 80% of their native habitat. Orange County is one of the last strongholds for this pond-breeding amphibian in the region, and ongoing restoration efforts targeting S. hammondii have involved habitat protection and the construction of artificial breeding ponds. These efforts have successfully increased breeding activity, but genetic characterization of the populations, including estimates of effective population size and admixture between the gene pools of constructed artificial and natural ponds, has never been undertaken. Using thousands of genome-wide single-nucleotide polymorphisms, we characterized the population structure, genetic diversity, and genetic connectivity of spadefoots in Orange County to guide ongoing and future management efforts. We identified at least 2, and possibly 3 major genetic clusters, with additional substructure within clusters indicating that individual ponds are often genetically distinct. Estimates of landscape resistance suggest that ponds on either side of the Los Angeles Basin were likely interconnected historically, but intense urban development has rendered them essentially isolated, and the resulting risk of interruption to natural metapopulation dynamics appears to be high. Resistance surfaces show that the existing artificial ponds were well-placed and connected to natural populations by low-resistance corridors. Toad samples from all ponds (natural and artificial) returned extremely low estimates of effective population size, possibly due to a bottleneck caused by a recent multi-year drought. Management efforts should focus on maintaining gene flow among natural and artificial ponds by both assisted migration and construction of new ponds to bolster the existing pond network in the region.

Integration of georeferenced and genetic data for the management of biodiversity in sheep genetic resources in Brazil

There are few animal germplasm/gene bank collections in Brazil, and basic studies are needed to attend the future internal and external demands from international partners. The aim of this work was to validate a "proof of concept" that integrates spatial (georeferenced data) and genetic data regarding the local of origin from 3518 DNA samples from 17 different genetic groups or breeds of sheep in the Brazilian Germplasm bank. Spatialisation shows that not all genetic groups have samples in the bank, and collection is concentrated in the conservation nuclei spread nationwide. Only 21% of states with a specific breed have samples in the gene bank. The mean number of animals sampled per collection was 32, while the mean distance travelled to collect samples was 262 km from the conservation nuclei. For example, the Brazilian Somali were only collected in the conservation nucleus in Ceará State. No samples were collected to date for the Cariri breed, which is recognised by the Brazilian Ministry of Agriculture. Only two farms and one breed in the bank are from the northern region. Of the 27 states, there are samples in the gene bank of sheep from 13, so several states have no samples, requiring collection from herds outside the official system of conservation to make sure that studies using this germplasm realised are not biased. Significant genetic differences are seen above 332 km, which should guide future sampling efforts. Suggestions are given for improving the quantity, quality and diversity of samples in the gene bank.

Transcriptome characterization and generation of marker resource for Himalayan vulnerable species, Ulmus wallichiana

Ulmus wallichiana is a traditional medicinal plant listed as a vulnerable in the IUCN red list data. Genomic and transcriptomic resources for this species are lacking, hindering its genetic exploration. Further, no polymorphic marker resource is available for this species, thus limiting the elucidation of its underlying genetic diversity, which is a pre-requisite for its conservation. This study was therefore aimed to generate a functionally annotated transcriptomic resource and screen it for SSR regions. We used paired-end Illumina based RNAseq technology and trinity based de novo assembly approach to generate full length transcripts, which were screened for SSR regions and functionally annotated. Around 6.6 million raw reads were de novo assembled transcripts, which were clustered into 146,083 unigenes. 19,909 transcripts were provided with 3986 unique KEGG ids, 70,519 transcripts with 6621 unique Pfam domains, and 45,125 transcripts with 7302 unique INTERPRO domains. 1456 transcripts were identified as transcriptions factors (TFs). Further, 8868 unique GO terms were obtained for the unigenes. The transcripts mapped to 23,056 known pre-determined orthology clusters in the eggNOG database. A total of 16,570 SSRs were identified from the unigenes. Out of the 90 SSRs selected for characterization on 20 genotypes, 28 were polymorphic. Mean effective alleles (Ne) of 2.53, mean observed heterozygosity (Ho) of 0.77, and average polymorphic information content (PIC) of 0.57 were found. This study may facilitate the genetic exploration of this species. The polymorphic SSRs would prove useful to explore its genetic diversity patterns, required for its conservation.

Type and intensity of surrounding human land use, not local environment, shape genetic structure of a native grassland plant

Landscape heterogeneity can shape genetic structure and functional connectivity of populations. When this heterogeneity imposes variable costs of moving across the landscape, populations can be structured according to a pattern of "isolation by resistance" (IBR). At the same time, divergent local environmental filters can limit gene flow, creating an alternative pattern of "isolation by environment" (IBE). Here, we evaluate IBR and IBE in the insect-pollinated, biennial plant Sabatia angularis (L.) Pursh (Gentianaceae) across serpentine grasslands in the fragmented landscape of SE Pennsylvania, USA using ~4500 neutral SNP loci. Specifically, we test the extent to which radical alteration of the landscape matrix by humans has fundamentally altered the cost of movement, imprinting a pattern of IBR dictated by land use type and intensity, and the potential for IBE in relation to a gradient of heavy metal concentrations found in serpentine soil. We reveal a strong signal of IBR and a weak signal of IBE across sites, indicating the greater importance of the landscape matrix in shaping genetic structure of S. angularis populations in the study region. Based on Circuitscape and least cost path approaches, we find that both low- and high-intensity urbanization resist gene flow by orders of magnitude greater than "natural" habitats, although resistance to low-intensity urbanization weakens at larger spatial scales. While cropland presents a substantially lower barrier than urban development, cumulative human land use surrounding populations predicts within-population genetic diversity and inbreeding in S. angularis. Our results emphasize the role of forest buffers and corridors in facilitating gene flow between serpentine grassland patches and averting local extinction of plant populations.

Fine-scale genetic structure in the critically endangered red-fronted macaw in the absence of geographic and ecological barriers

Behavioural and socio-cultural traits are recognized in the restriction of gene flow in species with high cognitive capacity and complex societies. This isolation by social barriers has been generally overlooked in threatened species by assuming disrupted gene flow due to population fragmentation and decline. We examine the genetic structure and ecology of the global population of the Critically Endangered red-fronted macaw (Ara rubrogenys), an endemic species to the inter-Andean valleys of Bolivia. We found a fine-scale genetic structuring in four genetic clusters. Genetic diversity was higher in wild compared to captive-bred macaws, but similar to that of captive wild-caught macaws. We found no clear evidence of severe genetic erosion in the population in recent decades, but it was patent in historic times, overlapping with drastic human habitat transformation and macaw persecution over millennia. We found no evidence of geographical and ecological barriers, owing to the high dispersal ability, nesting and foraging habits between genetic clusters. The lack of genetic intermixing despite long-distance foraging and seasonal movements suggests recruitment in natal colonies and other social factors reinforcing philopatry-related genetic structure. Conservation efforts should be specifically focussed on major threats in each genetic cluster as independent conservation units, and also considered in ex-situ management.