Preferred habitat and effective population size drive landscape genetic patterns in an endangered species

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.

Assessing fine-scale pondscape connectivity with amphibian eyes: An integrative approach using genomic and capture-mark-recapture data

In the face of habitat loss, preserving functional connectivity is essential to maintain genetic diversity and the demographic dynamics required for the viability of biotic communities. This requires knowledge of the dispersal behaviour of target species, which can be modelled as kernels, or probability density functions of dispersal distances at increasing geographic distances. We present an integrative approach to investigate the relationships between genetic connectivity and demographic parameters in organisms with low vagility focusing on five syntopic pond-breeding amphibians. We genotyped 1056 individuals of two anuran and three urodele species (1732-3913 SNPs per species) from populations located in a landscape comprising 64 ponds to characterize fine-scale genetic structure in a comparative framework, and combined these genetic data with information obtained in a previous 2-year capture-mark-recapture (CMR) study. Specifically, we contrasted graphs reconstructed from genomic data with connectivity graphs based on dispersal kernels and demographic information obtained from CMR data from previous studies, and assessed the effects of population size, population density, geographical distances, inverse movement probabilities and the presence of habitat patches potentially functioning as stepping stones on genetic differentiation. Our results show a significant effect of local population sizes on patterns of genetic differentiation at small spatial scales. In addition, movement records and cluster-derived kernels provide robust inferences on most likely dispersal paths that are consistent with genomic inferences on genetic connectivity. The integration of genetic and CMR data holds great potential for understanding genetic connectivity at spatial scales relevant to individual organisms, with applications for the implementation of management actions at the landscape level.

Movement of an imperiled esocid fish in an agricultural drain

Animal movement is increasingly affected by human alterations to habitat and climate change. In wetland systems, widespread hydrologic alterations from agriculture have changed the shape, function, and stability of shallow streams and wetland habitats. These changes in habitat quality and quantity may be especially consequential for freshwater fishes such as Grass Pickerel (Esox americanus vermiculatus), a small predatory fish found in disjunct populations across southern Ontario and listed as Special Concern under Canada's Species at Risk Act. To characterize Grass Pickerel movement response to stream-channel alterations, Fisheries and Oceans Canada implemented a tracking study to monitor the movements of a Grass Pickerel population in an agricultural drain on the Niagara Peninsula (Ontario, Canada). From 2009 to 2013, 2007 Grass Pickerel were tagged and tracked in the 37.3 km2 Beaver Creek watershed using a combination of mark-recapture surveys and eight fully automated passive integrated transponder tag antennas. Most individuals moved within 500 m (i.e., stationary fish) while 16% of the fish moved > 500 m (i.e., mobile fish), with a maximum median movement distance of 1.89 km and a maximum movement distance of 13.5 km (a long-tail distribution). Most movements occurred near the largest confluence where only a few were long-distance upstream or downstream movements. Mobile fish were larger than their stationary counterparts. Grass Pickerel in sites with higher abundance had more mobile fish, implying potential density dependence. Our results highlight that, while a long-distance dispersal ability exists in extant Grass Pickerel populations, the current conditions of riverscapes may prevent these dispersals from occurring. For declining Grass Pickerel populations, limitations to their movement ecology may substantially increase the likelihood of local extirpations.

Inferring landscape resistance to gene flow when genetic drift is spatially heterogeneous

In connectivity models, land cover types are assigned cost values characterizing their resistance to species movements. Landscape genetic methods infer these values from the relationship between genetic differentiation and cost distances. The spatial heterogeneity of population sizes, and consequently genetic drift, is rarely included in this inference although it influences genetic differentiation. Similarly, migration rates and population spatial distributions potentially influence this inference. Here, we assessed the reliability of cost value inference under several migration rates, population spatial patterns and degrees of population size heterogeneity. Additionally, we assessed whether considering intra-population variables, here using gravity models, improved the inference when drift is spatially heterogeneous. We simulated several gene flow intensities between populations with varying local sizes and spatial distributions. We then fit gravity models of genetic distances as a function of (i) the 'true' cost distances driving simulations or alternative cost distances, and (ii) intra-population variables (population sizes, patch areas). We determined the conditions making the identification of the 'true' costs possible and assessed the contribution of intra-population variables to this objective. Overall, the inference ranked cost scenarios reliably in terms of similarity with the 'true' scenario (cost distance Mantel correlations), but this 'true' scenario rarely provided the best model goodness of fit. Ranking inaccuracies and failures to identify the 'true' scenario were more pronounced when migration was very restricted (<4 dispersal events/generation), population sizes were most heterogeneous and some populations were spatially aggregated. In these situations, considering intra-population variables helps identify cost scenarios reliably, thereby improving cost value inference from genetic data.

The Genomic Landscapes of Desert Birds Form over Multiple Time Scales

Spatial models show that genetic differentiation between populations can be explained by factors ranging from geographic distance to environmental resistance across the landscape. However, genomes exhibit a landscape of differentiation, indicating that multiple processes may mediate divergence in different portions of the genome. We tested this idea by comparing alternative geographic predctors of differentiation in ten bird species that co-occur in Sonoran and Chihuahuan Deserts of North America. Using population-level genomic data, we described the genomic landscapes across species and modeled conditions that represented historical and contemporary mechanisms. The characteristics of genomic landscapes differed across species, influenced by varying levels of population structuring and admixture between deserts, and the best-fit models contrasted between the whole genome and partitions along the genome. Both historical and contemporary mechanisms were important in explaining genetic distance, but particularly past and current environments, suggesting that genomic evolution was modulated by climate and habitat There were also different best-ftit models across genomic partitions of the data, indicating that these regions capture different evolutionary histories. These results show that the genomic landscape of differentiation can be associated with alternative geographic factors operating on different portions of the genome, which reflect how heterogeneous patterns of genetic differentiation can evolve across species and genomes.

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.

Multisource noninvasive genetics of brown bears (Ursus arctos) in Greece reveals a highly structured population and a new matrilineal contact zone in southern Europe

In human-dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244-bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos-Dinara populations. Total effective population size (N e) was estimated to be 199 individuals, and total combined population size (N C) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe's low-latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long-term persistence of the species in the Balkans.

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.

Minimum habitat thresholds required for conserving mountain lion genetic diversity

Jointly considering the ecology (e.g., habitat use) and genetics (e.g., population genetic structure and diversity) of a species can increase understanding of current conservation status and inform future management practices. Previous analyses indicate that mountain lion (Puma concolor) populations in California are genetically structured and exhibit extreme variation in population genetic diversity. Although human development may have fragmented gene flow, we hypothesized the quantity and quality of remaining habitat available would affect the genetic viability of each population. Our results indicate that area of suitable habitat, determined via a resource selection function derived using 843,500 location fixes from 263 radio‐collared mountain lions, is strongly and positively associated with population genetic diversity and viability metrics, particularly with effective population size. Our results suggested that contiguous habitat of ≥10,000 km² may be sufficient to alleviate the negative effects of genetic drift and inbreeding, allowing mountain lion populations to maintain suitable effective population sizes. Areas occupied by five of the nine geographic–genetic mountain lion populations in California fell below this habitat threshold, and two (Santa Monica Area and Santa Ana) of those five populations lack connectivity to nearby populations. Enhancing ecological conditions by protection of greater areas of suitable habitat and facilitating positive evolutionary processes by increasing connectivity (e.g., road‐crossing structures) might promote persistence of small or isolated populations. The conservation status of suitable habitat also appeared to influence genetic diversity of populations. Thus, our results demonstrate that both the area and status (i.e., protected or unprotected) of suitable habitat influence the genetic viability of mountain lion populations.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

Genomics, environment and balancing selection in behaviourally bimodal populations: The caribou case

Selection forces that favour different phenotypes in different environments can change frequencies of genes between populations along environmental clines. Clines are also compatible with balancing forces, such as negative frequency-dependent selection (NFDS), which maintains phenotypic polymorphisms within populations. For example, NFDS is hypothesized to maintain partial migration, a dimorphic behavioural trait prominent in species where only a fraction of the population seasonally migrates. Overall, NFDS is believed to be a common phenomenon in nature, yet a scarcity of studies were published linking naturally occurring allelic variation with bimodal or multimodal phenotypes and balancing selection. We applied a Pool-seq approach and detected selection on alleles associated with environmental variables along a North-South gradient in western North American caribou, a species displaying partially migratory behaviour. On 51 loci, we found a signature of balancing selection, which could be related to NFDS and ultimately the maintenance of the phenotypic polymorphisms known within these populations. Yet, remarkably, we detected directional selection on a locus when our sample was divided into two behaviourally distinctive groups regardless of geographic provenance (a subset of GPS-collared migratory or sedentary individuals), indicating that, within populations, phenotypically homogeneous groups were genetically distinctive. Loci under selection were linked to functional genes involved in oxidative stress response, body development and taste perception. Overall, results indicated genetic differentiation along an environmental gradient of caribou populations, which we found characterized by genes potentially undergoing balancing selection. We suggest that the underlining balancing force, NFDS, plays a strong role within populations harbouring multiple haplotypes and phenotypes, as it is the norm in animals, plants and humans too.

Quantifying functional connectivity: The role of breeding habitat, abundance, and landscape features on range-wide gene flow in sage-grouse

Functional connectivity, quantified using landscape genetics, can inform conservation through the identification of factors linking genetic structure to landscape mechanisms. We used breeding habitat metrics, landscape attributes, and indices of grouse abundance, to compare fit between structural connectivity and genetic differentiation within five long-established Sage-Grouse Management Zones (MZ) I-V using microsatellite genotypes from 6,844 greater sage-grouse (Centrocercus urophasianus) collected across their 10.7 million-km2 range. We estimated structural connectivity using a circuit theory-based approach where we built resistance surfaces using thresholds dividing the landscape into "habitat" and "nonhabitat" and nodes were clusters of sage-grouse leks (where feather samples were collected using noninvasive techniques). As hypothesized, MZ-specific habitat metrics were the best predictors of differentiation. To our surprise, inclusion of grouse abundance-corrected indices did not greatly improve model fit in most MZs. Functional connectivity of breeding habitat was reduced when probability of lek occurrence dropped below 0.25 (MZs I, IV) and 0.5 (II), thresholds lower than those previously identified as required for the formation of breeding leks, which suggests that individuals are willing to travel through undesirable habitat. The individual MZ landscape results suggested terrain roughness and steepness shaped functional connectivity across all MZs. Across respective MZs, sagebrush availability (<10%-30%; II, IV, V), tree canopy cover (>10%; I, II, IV), and cultivation (>25%; I, II, IV, V) each reduced movement beyond their respective thresholds. Model validations confirmed variation in predictive ability across MZs with top resistance surfaces better predicting gene flow than geographic distance alone, especially in cases of low and high differentiation among lek groups. The resultant resistance maps we produced spatially depict the strength and redundancy of range-wide gene flow and can help direct conservation actions to maintain and restore functional connectivity for sage-grouse.

The relative contribution of natural landscapes and human-mediated factors on the connectivity of a noxious invasive weed

Examining how the landscape may influence gene flow is at the forefront of understanding population differentiation and adaptation. Such understanding is crucial in light of ongoing environmental changes and the elevated risk of ecosystems alteration. In particular, knowledge of how humans may influence population structure is imperative to allow for informed decisions in management and conservation as well as to gain a better understanding of anthropogenic impacts on the interplay between gene flow, genetic drift, and selection. Here, we use genome-wide molecular markers to characterize the population genetic structure and connectivity of Ipomoea purpurea (Convolvulaceae), a noxious invasive weed. We, likewise, assess the interaction between natural and human-driven influences on genetic differentiation among populations. Our analyses find that human population density is an important predictor of pairwise population differentiation, suggesting that the agricultural and/or horticultural trade may be involved in maintaining some level of connectivity across distant agricultural fields. Climatic variation appears as an additional predictor of genetic connectivity in this species. We discuss the implications of these results and highlight future research needed to disentangle the mechanistic processes underlying population connectivity of weeds.

Loss of dendritic connectivity in southern California's urban riverscape facilitates decline of an endemic freshwater fish

Life history adaptations and spatial configuration of metapopulation networks allow certain species to persist in extreme fluctuating environments, yet long-term stability within these systems relies on the maintenance of linkage habitat. Degradation of such linkages in urban riverscapes can disrupt this dynamic in aquatic species, leading to increased extinction debt in local populations experiencing environment-related demographic flux. We used microsatellites and mtDNA to examine the effects of collapsed network structure in the endemic Santa Ana sucker Catostomus santaanae of southern California, a threatened species affected by natural flood-drought cycles, "boom-and-bust" demography, hybridization and presumed artificial transplantation. Our results show a predominance of drift-mediated processes in shaping population structure and that reverse mechanisms for counterbalancing the genetic effects of these phenomena have dissipated with the collapse of dendritic connectivity. We use approximate Bayesian models to support two cases of artificial transplantation and provide evidence that one of the invaded systems better represents the historic processes that maintained genetic variation within watersheds than any remaining drainages where C. santaanae is considered native. We further show that a stable dry gap in the northern range is preventing genetic dilution of pure C. santaanae persisting upstream of a hybrid assemblage involving a non-native sucker and that local accumulation of genetic variation in the same drainage is influenced by position within the network. This work has important implications for declining species that have historically relied on dendritic metapopulation networks to maintain source-sink dynamics in phasic environments, but no longer possess this capacity in urban-converted landscapes.

Genetic footprint of population fragmentation and contemporary collapse in a freshwater cetacean

Understanding demographic trends and patterns of gene flow in an endangered species is crucial for devising conservation strategies. Here, we examined the extent of population structure and recent evolution of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis). By analysing genetic variation at the mitochondrial and nuclear microsatellite loci for 148 individuals, we identified three populations along the Yangtze River, each one connected to a group of admixed ancestry. Each population displayed extremely low genetic diversity, consistent with extremely small effective size (≤92 individuals). Habitat degradation and distribution gaps correlated with highly asymmetric gene-flow that was inefficient in maintaining connectivity between populations. Genetic inferences of historical demography revealed that the populations in the Yangtze descended from a small number of founders colonizing the river from the sea during the last Ice Age. The colonization was followed by a rapid population split during the last millennium predating the Chinese Modern Economy Development. However, genetic diversity showed a clear footprint of population contraction over the last 50 years leaving only ~2% of the pre-collapsed size, consistent with the population collapses reported from field studies. This genetic perspective provides background information for devising mitigation strategies to prevent this species from extinction.

River network architecture, genetic effective size and distributional patterns predict differences in genetic structure across species in a dryland stream fish community

Dendritic ecological network (DEN) architecture can be a strong predictor of spatial genetic patterns in theoretical and simulation studies. Yet, interspecific differences in dispersal capabilities and distribution within the network may equally affect species' genetic structuring. We characterized patterns of genetic variation from up to ten microsatellite loci for nine numerically dominant members of the upper Gila River fish community, New Mexico, USA. Using comparative landscape genetics, we evaluated the role of network architecture for structuring populations within species (pairwise F_ST ) while explicitly accounting for intraspecific demographic influences on effective population size (N_e ). Five species exhibited patterns of connectivity and/or genetic diversity gradients that were predicted by network structure. These species were generally considered to be small-bodied or habitat specialists. Spatial variation of N_e was a strong predictor of pairwise F_ST for two species, suggesting patterns of connectivity may also be influenced by genetic drift independent of network properties. Finally, two study species exhibited genetic patterns that were unexplained by network properties and appeared to be related to nonequilibrium processes. Properties of DENs shape community-wide genetic structure but effects are modified by intrinsic traits and nonequilibrium processes. Further theoretical development of the DEN framework should account for such cases.

Home ranges, habitat and body mass: simple correlates of home range size in ungulates

The spatial scale of animal space use, e.g. measured as individual home range size, is a key trait with important implications for ecological and evolutionary processes as well as management and conservation of populations and ecosystems. Explaining variation in home range size has therefore received great attention in ecological research. However, few studies have examined multiple hypotheses simultaneously, which is important provided the complex interactions between life history, social system and behaviour. Here, we review previous studies on home range size in ungulates, supplementing with a meta-analysis, to assess how differences in habitat use and species characteristics affect the relationship between body mass and home range size. Habitat type was the main factor explaining interspecific differences in home range size after accounting for species body mass and group size. Species using open habitats had larger home ranges for a given body mass than species using closed habitats, whereas species in open habitats showed a much weaker allometric relationship compared with species living in closed habitats. We found no support for relationships between home range size and species diet or mating system, or any sexual differences. These patterns suggest that the spatial scale of animal movement mainly is a combined effect of body mass, group size and the landscape structure. Accordingly, landscape management must acknowledge the influence of spatial distribution of habitat types on animal behaviour to ensure natural processes affecting demography and viability of ungulate populations.

Environmental and anthropogenic drivers of connectivity patterns: A basis for prioritizing conservation efforts for threatened populations

Ecosystem fragmentation and habitat loss have been the focus of landscape management due to restrictions on contemporary connectivity and dispersal of populations. Here, we used an individual approach to determine the drivers of genetic differentiation in caribou of the Canadian Rockies. We modelled the effects of isolation by distance, landscape resistance and predation risk and evaluated the consequences of individual migratory behaviour (seasonally migratory vs. sedentary) on gene flow in this threatened species. We applied distance-based and reciprocal causal modelling approaches, testing alternative hypotheses on the effects of geographic, topographic, environmental and local population-specific variables on genetic differentiation and relatedness among individuals. Overall, gene flow was restricted to neighbouring local populations, with spatial coordinates, local population size, groups and elevation explaining connectivity among individuals. Landscape resistance, geographic distances and predation risk were correlated with genetic distances, with correlations threefold higher for sedentary than for migratory caribou. As local caribou populations are increasingly isolated, our results indicate the need to address genetic connectivity, especially for populations with individuals displaying different migratory behaviours, whilst maintaining quality habitat both within and across the ranges of threatened populations.

Mapping Resource Selection Functions in Wildlife Studies: Concerns and Recommendations

Predicting the spatial distribution of animals is an important and widely used tool with applications in wildlife management, conservation, and population health. Wildlife telemetry technology coupled with the availability of spatial data and GIS software have facilitated advancements in species distribution modeling. There are also challenges related to these advancements including the accurate and appropriate implementation of species distribution modeling methodology. Resource Selection Function (RSF) modeling is a commonly used approach for understanding species distributions and habitat usage, and mapping the RSF results can enhance study findings and make them more accessible to researchers and wildlife managers. Currently, there is no consensus in the literature on the most appropriate method for mapping RSF results, methods are frequently not described, and mapping approaches are not always related to accuracy metrics. We conducted a systematic review of the RSF literature to summarize the methods used to map RSF outputs, discuss the relationship between mapping approaches and accuracy metrics, performed a case study on the implications of employing different mapping methods, and provide recommendations as to appropriate mapping techniques for RSF studies. We found extensive variability in methodology for mapping RSF results. Our case study revealed that the most commonly used approaches for mapping RSF results led to notable differences in the visual interpretation of RSF results, and there is a concerning disconnect between accuracy metrics and mapping methods. We make 5 recommendations for researchers mapping the results of RSF studies, which are focused on carefully selecting and describing the method used to map RSF studies, and relating mapping approaches to accuracy metrics.

Evaluating the influence of life-history characteristics on genetic structure: a comparison of small mammals inhabiting complex agricultural landscapes

Conversion of formerly continuous native habitats into highly fragmented landscapes can lead to numerous negative demographic and genetic impacts on native taxa that ultimately reduce population viability. In response to concerns over biodiversity loss, numerous investigators have proposed that traits such as body size and ecological specialization influence the sensitivity of species to habitat fragmentation. In this study, we examined how differences in body size and ecological specialization of two rodents (eastern chipmunk; Tamias striatus and white‐footed mouse; Peromyscus leucopus) impact their genetic connectivity within the highly fragmented landscape of the Upper Wabash River Basin (UWB), Indiana, and evaluated whether landscape configuration and complexity influenced patterns of genetic structure similarly between these two species. The more specialized chipmunk exhibited dramatically more genetic structure across the UWB than white‐footed mice, with genetic differentiation being correlated with geographic distance, configuration of intervening habitats, and complexity of forested habitats within sampling sites. In contrast, the generalist white‐footed mouse resembled a panmictic population across the UWB, and no landscape factors were found to influence gene flow. Despite the extensive previous work in abundance and occupancy within the UWB, no landscape factor that influenced occupancy or abundance was correlated with genetic differentiation in either species. The difference in predictors of occupancy, abundance, and gene flow suggests that species‐specific responses to fragmentation are scale dependent.

Identification of landscape features influencing gene flow: How useful are habitat selection models?

Understanding how dispersal patterns are influenced by landscape heterogeneity is critical for modeling species connectivity. Resource selection function (RSF) models are increasingly used in landscape genetics approaches. However, because the ecological factors that drive habitat selection may be different from those influencing dispersal and gene flow, it is important to consider explicit assumptions and spatial scales of measurement. We calculated pairwise genetic distance among 301 Dall's sheep (Ovis dalli dalli) in southcentral Alaska using an intensive noninvasive sampling effort and 15 microsatellite loci. We used multiple regression of distance matrices to assess the correlation of pairwise genetic distance and landscape resistance derived from an RSF, and combinations of landscape features hypothesized to influence dispersal. Dall's sheep gene flow was positively correlated with steep slopes, moderate peak normalized difference vegetation indices (NDVI), and open land cover. Whereas RSF covariates were significant in predicting genetic distance, the RSF model itself was not significantly correlated with Dall's sheep gene flow, suggesting that certain habitat features important during summer (rugged terrain, mid-range elevation) were not influential to effective dispersal. This work underscores that consideration of both habitat selection and landscape genetics models may be useful in developing management strategies to both meet the immediate survival of a species and allow for long-term genetic connectivity.

Repeated Habitat Disturbances by Fire Decrease Local Effective Population Size

Effective population size is a fundamental parameter in population genetics, and factors that alter effective population size will shape the genetic characteristics of populations. Habitat disturbance may have a large effect on genetic characteristics of populations by influencing immigration and gene flow, particularly in fragmented habitats. We used the Florida Sand Skink (Plestiodon reynoldsi) to investigate the effect of fire-based habitat disturbances on the effective population size in the highly threatened, severely fragmented, and fire dependent Florida scrub habitat. We screened 7 microsatellite loci in 604 individuals collected from 12 locations at Archbold Biological Station. Archbold Biological Station has an active fire management plan and detailed records of fires dating to 1967. Our objective was to determine how the timing, number, and intervals between fires affect effective population size, focusing on multiple fires in the same location. Effective population size was higher in areas that had not been burned for more than 10 years and decreased with number of fires and shorter time between fires. A similar pattern was observed in abundance: increasing abundance with time-since-fire and decreasing abundance with number of fires. The ratio of effective population size to census size was higher at sites with more recent fires and tended to decrease with time-since-last-fire. These results suggest that habitat disturbances, such as fire, may have a large effect in the genetic characteristics of local populations and that Florida Sand Skinks are well adapted to the natural fire dynamics required to maintain Florida scrub.

Anthropogenic Habitats Facilitate Dispersal of an Early Successional Obligate: Implications for Restoration of an Endangered Ecosystem

Landscape modification and habitat fragmentation disrupt the connectivity of natural landscapes, with major consequences for biodiversity. Species that require patchily distributed habitats, such as those that specialize on early successional ecosystems, must disperse through a landscape matrix with unsuitable habitat types. We evaluated landscape effects on dispersal of an early successional obligate, the New England cottontail (Sylvilagus transitionalis). Using a landscape genetics approach, we identified barriers and facilitators of gene flow and connectivity corridors for a population of cottontails in the northeastern United States. We modeled dispersal in relation to landscape structure and composition and tested hypotheses about the influence of habitat fragmentation on gene flow. Anthropogenic and natural shrubland habitats facilitated gene flow, while the remainder of the matrix, particularly development and forest, impeded gene flow. The relative influence of matrix habitats differed between study areas in relation to a fragmentation gradient. Barrier features had higher explanatory power in the more fragmented site, while facilitating features were important in the less fragmented site. Landscape models that included a simultaneous barrier and facilitating effect of roads had higher explanatory power than models that considered either effect separately, supporting the hypothesis that roads act as both barriers and facilitators at all spatial scales. The inclusion of LiDAR-identified shrubland habitat improved the fit of our facilitator models. Corridor analyses using circuit and least cost path approaches revealed the importance of anthropogenic, linear features for restoring connectivity between the study areas. In fragmented landscapes, human-modified habitats may enhance functional connectivity by providing suitable dispersal conduits for early successional specialists.

Forecasting Ecological Genomics: High-Tech Animal Instrumentation Meets High-Throughput Sequencing

Recent advancements in animal tracking technology and high-throughput sequencing are rapidly changing the questions and scope of research in the biological sciences. The integration of genomic data with high-tech animal instrumentation comes as a natural progression of traditional work in ecological genetics, and we provide a framework for linking the separate data streams from these technologies. Such a merger will elucidate the genetic basis of adaptive behaviors like migration and hibernation and advance our understanding of fundamental ecological and evolutionary processes such as pathogen transmission, population responses to environmental change, and communication in natural populations.

Spatial Niche Segregation of Sympatric Stone Marten and Pine Marten--Avoidance of Competition or Selection of Optimal Habitat?

Coexistence of ecologically similar species relies on differences in one or more dimensions of their ecological niches, such as space, time and resources in diel and/or seasonal scales. However, niche differentiation may result from other mechanisms such as avoidance of high predation pressure, different adaptations or requirements of ecologically similar species. Stone marten (Martes foina) and pine marten (Martes martes) occur sympatrically over a large area in Central Europe and utilize similar habitats and food, therefore it is expected that their coexistence requires differentiation in at least one of their niche dimensions or the mechanisms through which these dimensions are used. To test this hypothesis, we used differences in the species activity patterns and habitat selection, estimated with a resource selection function (RSF), to predict the relative probability of occurrence of the two species within a large forest complex in the northern geographic range of the stone marten. Stone martens were significantly heavier, have a longer body and a better body condition than pine martens. We found weak evidence for temporal niche segregation between the species. Stone and pine martens were both primarily nocturnal, but pine martens were active more frequently during the day and significantly reduced the duration of activity during autumn-winter. Stone and pine martens utilized different habitats and almost completely separated their habitat niches. Stone marten strongly preferred developed areas and avoided meadows and coniferous or deciduous forests. Pine marten preferred deciduous forest and small patches covered by trees, and avoided developed areas and meadows. We conclude that complete habitat segregation of the two marten species facilitates sympatric coexistence in this area. However, spatial niche segregation between these species was more likely due to differences in adaptation to cold climate, avoidance of high predator pressure and/or food preferences by both species than competitive interaction between them.

Temporally dynamic habitat suitability predicts genetic relatedness among caribou

Landscape heterogeneity plays a central role in shaping ecological and evolutionary processes. While species utilization of the landscape is usually viewed as constant within a year, the spatial distribution of individuals is likely to vary in time in relation to particular seasonal needs. Understanding temporal variation in landscape use and genetic connectivity has direct conservation implications. Here, we modelled the daily use of the landscape by caribou in Quebec and Labrador, Canada and tested its ability to explain the genetic relatedness among individuals. We assessed habitat selection using locations of collared individuals in migratory herds and static occurrences from sedentary groups. Connectivity models based on habitat use outperformed a baseline isolation-by-distance model in explaining genetic relatedness, suggesting that variations in landscape features such as snow, vegetation productivity and land use modulate connectivity among populations. Connectivity surfaces derived from habitat use were the best predictors of genetic relatedness. The relationship between connectivity surface and genetic relatedness varied in time and peaked during the rutting period. Landscape permeability in the period of mate searching is especially important to allow gene flow among populations. Our study highlights the importance of considering temporal variations in habitat selection for optimizing connectivity across heterogeneous landscape and counter habitat fragmentation.

Predator experience overrides learned aversion to heterospecifics in stickleback species pairs

Predation risk can alter female mating decisions because the costs of mate searching and selecting attractive mates increase when predators are present. In response to predators, females have been found to plastically adjust mate preference within species, but little is known about how predators alter sexual isolation and hybridization among species. We tested the effects of predator exposure on sexual isolation between benthic and limnetic threespine sticklebacks (Gasterosteus spp.). Female discrimination against heterospecific mates was measured before and after females experienced a simulated attack by a trout predator or a control exposure to a harmless object. In the absence of predators, females showed increased aversion to heterospecifics over time. We found that predator exposure made females less discriminating and precluded this learned aversion to heterospecifics. Benthic and limnetic males differ in coloration, and predator exposure also affected sexual isolation by weakening female preferences for colourful males. Predator effects on sexual selection were also tested but predators had few effects on female choosiness among conspecific mates. Our results suggest that predation risk may disrupt the cognitive processes associated with mate choice and lead to fluctuations in the strength of sexual isolation between species.

The influence of habitat structure on genetic differentiation in red fox populations in north-eastern Poland

The red fox (Vulpes vulpes) has the widest global distribution among terrestrial carnivore species, occupying most of the Northern Hemisphere in its native range. Because it carries diseases that can be transmitted to humans and domestic animals, it is important to gather information about their movements and dispersal in their natural habitat but it is difficult to do so at a broad scale with trapping and telemetry. In this study, we have described the genetic diversity and structure of red fox populations in six areas of north-eastern Poland, based on samples collected from 2002–2003. We tested 22 microsatellite loci isolated from the dog and the red fox genome to select a panel of nine polymorphic loci suitable for this study. Genetic differentiation between the six studied populations was low to moderate and analysis in Structure revealed a panmictic population in the region. Spatial autocorrelation among all individuals showed a pattern of decreasing relatedness with increasing distance and this was not significantly negative until 93 km, indicating a pattern of isolation-by-distance over a large area. However, there was no correlation between genetic distance and either Euclidean distance or least-cost path distance at the population level. There was a significant relationship between genetic distance and the proportion of large forests and water along the Euclidean distances. These types of habitats may influence dispersal paths taken by red foxes, which is useful information in terms of wildlife disease management.