Consequences of extensive habitat fragmentation in landscape-level patterns of genetic diversity and structure in the Mediterranean esparto grasshopper

Niche modelling and landscape genetics of the yellow-legged hornet (Vespa velutina): An integrative approach for evaluating central-marginal population dynamics in Europe

Genetic diversity is an important biological trait for a successful invasion. During the expansion across a new territory, an invasive species may face unprecedented ecological conditions that will determine its demography and genetic diversity. The first record of the yellow-legged hornet (Vespa velutina) in Europe dates back to 2004 in France, from where it has successfully spread through a large territory in the continent, including Italy, Spain and Portugal. Integrative approaches offer a powerful strategy to detect and understand patterns of genetic variation in central and marginal populations. Here, we have analysed the relationship between genetic diversity parameters inferred from 15 V. velutina nuclear DNA microsatellite loci, and geographical and environmental drivers, such as the distance to the introduction focus, environmental suitability and distance to native and invasive niche centroids. Our results revealed a central-marginal dynamic, where allelic richness decreased towards the edge of the expansion range. The low environmental suitability of the territories invaded by marginal populations could prevent a diverse population from establishing and reducing the genetic diversity in populations at the expansion edge. Moreover, Markov chain Monte Carlo analysis showed both geographical and environmental distances were influencing population genetic differentiation. This study highlights the importance of combining genetic analysis with geographical and environmental drivers to understand genetic trends of invasive species to new environment.

Analysis of Spatiotemporal Variation in Habitat Suitability for Oedaleus decorus asiaticus Bei-Bienko on the Mongolian Plateau Using Maxent and Multi-Source Remote Sensing Data

O. decorus asiaticus is a major grasshopper species that harms the development of agriculture on the Mongolian Plateau. Therefore, it is important to enhance the monitoring of O. decorus asiaticus. In this study, the spatiotemporal variation in the habitat suitability for O. decorus asiaticus on the Mongolian Plateau was assessed using maximum entropy (Maxent) modeling along with multi-source remote sensing data (meteorology, vegetation, soil, and topography). The predictions of the Maxent model were accurate (AUC = 0.910). The key environmental variables affecting the distribution of grasshoppers and their contribution were grass type (51.3%), accumulated precipitation (24.9%), altitude (13.0%), vegetation coverage (6.6%), and land surface temperature (4.2%). Based on the assessment results of suitability by Maxent model, the model threshold settings, and the formula for calculating the inhabitability index, the 2000s, 2010s, and 2020s inhabitable areas were calculated. The results show that the distribution of suitable habitat for O. decorus asiaticus in 2000 was similar to that in 2010. From 2010 to 2020, the suitability of the habitat for O. decorus asiaticus in the central region of the Mongolian Plateau changed from moderate to high. The main factor contributing to this change was accumulated precipitation. Few changes in the areas of the habitat with low suitability were observed across the study period. The results of this study enhance our understanding of the vulnerability of different regions on the Mongolian Plateau to plagues of O. decorus asiaticus and will aid the monitoring of grasshopper plagues in this region.

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.

Impact of Limited Dispersion Capacity and Natural Barriers on the Population Structure of the Grasshopper Ommexecha virens (Orthoptera: Ommexechidae)

The grasshopper Ommexecha virens Serville has low dispersion capacity, and it is regarded as a specialist, only being found in sandy, dry environments with high incidence of sunlight. Considering these aspects, we evaluated the diversity and genetic structure of O. virens natural populations using ISSR (Inter Simple Sequence Repeat) markers. The data pointed to low expected heterozygosity for some populations (H_E = 0.06-0.09), probably a consequence of positive inbreeding, which is typical of species showing low or null dispersion indices. Moreover, significant genetic differentiation was observed (F_ST = 0.50 and G_ST = 0.51), as well as low number of migrants (N_m = 0.47), indicating that the populations are genetically differentiated. This is likely related to the limitation in dispersing and fragmentation of suitable environment localities colonized by O. virens. The populations of O. virens were structured in three genetic groups associated to different landscapes, revealing the presence of a secondary contact zone, possibly arisen from isolation followed by genetic divergence among populations and subsequent gene flow of divergent individuals of O. virens. At last, we found positive isolation by distance (IBD; r: 0.427; P: 0.025) which is an important factor, since it may be adding to the emergence of reproductive barriers among individuals of O. virens that have been experiencing isolation.

Genomic data support multiple introductions and explosive demographic expansions in a highly invasive aquatic insect

The study of the genetic makeup and demographic fate of alien species is essential to understand their capacity to recover from founder effects, adapt to new environmental conditions and, ultimately, become invasive and potentially damaging. Here, we employ genomic data to gain insights into key demographic processes that might help to explain the extraordinarily successful invasion of the Western Mediterranean region by the North American boatman Trichocorixa verticalis (Hemiptera: Corixidae). Our analyses revealed the genetic distinctiveness of populations from the main areas comprising the invasive range and coalescent-based simulations supported that they originated from independent introductions events probably involving different source populations. Testing of alternative demographic models indicated that all populations experienced a strong bottleneck followed by a recent and instantaneous demographic expansion that restored a large portion (>30%) of their ancestral effective population sizes shortly after introductions took place (<60 years ago). Considerable genetic admixture of some populations suggest that hypothetical barriers to dispersal (i.e., land and sea water) are permeable to gene flow and/or that they originated from introductions involving multiple lineages. This study demonstrates the repeated arrival of propagules with different origins and short time lags between arrival and establishment, emphasizing the extraordinary capacity of the species to recover from founder effects and genetically admix in invaded areas. This can explain the demonstrated capacity of this aquatic insect to spread and outcompete native species once it colonizes new suitable regions. Future genomic analyses of native range populations could help to infer the genetic makeup of introduced populations and track invasion routes.

Genetic structure across urban and agricultural landscapes reveals evidence of resource specialization and philopatry in the Eastern carpenter bee, Xylocopa virginica L

Human activity continues to impact global ecosystems, often by altering the habitat suitability, persistence, and movement of native species. It is thus critical to examine the population genetic structure of key ecosystemservice providers across human-altered landscapes to provide insight into the forces that limit wildlife persistence and movement across multiple spatial scales. While some studies have documented declines of bee pollinators as a result of human-mediated habitat alteration, others suggest that some bee species may benefit from altered land use due to increased food or nesting resource availability; however, detailed population and dispersal studies have been lacking. We investigated the population genetic structure of the Eastern carpenter bee, Xylocopa virginica, across 14 sites spanning more than 450 km, including dense urban areas and intensive agricultural habitat. X. virginica is a large bee which constructs nests in natural and human-associated wooden substrates, and is hypothesized to disperse broadly across urbanizing areas. Using 10 microsatellite loci, we detected significant genetic isolation by geographic distance and significant isolation by land use, where urban and cultivated landscapes were most conducive to gene flow. This is one of the first population genetic analyses to provide evidence of enhanced insect dispersal in human-altered areas as compared to semi-natural landscapes. We found moderate levels of regional-scale population structure across the study system (G'ST = 0.146) and substantial co-ancestry between the sampling regions, where co-ancestry patterns align with major human transportation corridors, suggesting that human-mediated movement may be influencing regional dispersal processes. Additionally, we found a signature of strong site-level philopatry where our analyses revealed significant levels of high genetic relatedness at very fine scales (<1 km), surprising given X. virginica's large body size. These results provide unique genetic evidence that insects can simultaneously exhibit substantial regional dispersal as well as high local nesting fidelity in landscapes dominated by human activity.

Comparison of Genetic Structure between Endangered and Common Butterflies, Ypthima multistriata and Y. argus (Lepidoptera; Nymphalidae), Inhabiting Japan

Genetic structures of two closely related butterflies, Ypthima multistriata and Y. argus, inhabiting Japan were compared based on the mitochondrial cytochrome c oxidase subunit I (COI) sequences. The former species is classified as a vulnerable species and exhibits a characteristic pattern of voltinism: univoltine and bivoltine populations are distributed in a scattered manner. The latter species is common and has a normal geographical pattern of voltinism: the number of annual generations is correlated with latitude. Our genetic analyses of these two species yielded contrasting results: a spatial analysis of molecular variance (SAMOVA) and F_ST between each pair of populations revealed a locally fragmented genetic structure for Y. multistriata, compared to three distinct geographic groups of Y. argus within which range-wide gene flow occurs. Although Y. argus is a common species, only the southernmost populations in Japan had higher genetic diversity, while the other populations had the same or lower levels of genetic diversity, compared to Y. multistriata. These results indicate that: 1) the degree of fragmentation of Y. multistriata populations was higher; however, markedly lower genetic diversity was not found, and 2) although Y. argus is a common species, its populations may not be genetically robust. In addition, AMOVA revealed a relationship between voltinism and genetic variation in Y. multistriata. This result suggests a phylogenetic constraint of voltinism in this butterfly.

Wide dispersion of B chromosomes in Rhammatocerus brasiliensis (Orthoptera, Acrididae)

The grasshopper Rhammatocerus brasiliensis shows polymorphism of B chromosomes, but the magnitude of B-chromosome occurrence and the factors that may contribute to their dispersion in the species remain unknown thus far. The present study analyzed the occurrence and dispersion of B chromosomes in R. brasiliensis individuals from 21 populations widely distributed in the Brazilian Northeast. The genetic connectivity between 10 populations was verified through analysis of ISSR markers from 200 individuals. Of the 21 populations, 19 presented individuals with one B chromosome, three with two, and one with three B chromosomes. The B chromosome is of medium size and constitutive heterochromatin (CH) located in the pericentromeric region. A variant B chromosome was observed in three populations, similar in size to that of chromosome X, gap and CH, and located in the terminal region. B chromosome frequencies in different populations varied from 0% to 18,8%, mean 8,5%. The wide distribution of the B chromosome is likely a consequence of the positive gene flow among the analyzed populations. B-chromosome occurrence in populations of R. brasiliensis possibly follows the population genetic structure of the species and, owing to the existence of a variant, its origin may not be recent.

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus)

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site-associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.

Spatiotemporal Differentiation of Alpine Butterfly Parnassius glacialis (Papilionidae: Parnassiinae) in China: Evidence from Mitochondrial DNA and Nuclear Single Nucleotide Polymorphisms

The Apollo butterfly, Parnassius glacialis, is one of the most charming members of its genus and includes two subspecies locally distributed in montane areas of south-central China and Japan. In this study, we investigated the genetic structure and demographic history of P. glacialis by analyzing partial sequences of four mitochondrial genes and nuclear single nucleotide polymorphisms (SNPs) via genotyping-by-sequencing (GBS) of samples from nearly the entire known distributional range in China. The mitochondrial DNA (mtDNA) data demonstrated that a total of 39 haplotypes were present, and the species was estimated to have diverged about 0.95 million years ago during the middle Pleistocene transition into two main clades that likely formed during the Kunlun-Huanghe tectonic movement. The two clades then dispersed independently in distinct geographic areas alongside the mountainous routes in central and southern China, most likely driven by the Pleistocene glacial-interglacial cycles. Nuclear SNP analysis was generally congruent with mtDNA results at the individual level. A minor incongruence of genetic structures that was detected between mtDNA and nuclear SNP data from the Laojunshan and Tiantangzhai populations was likely due to secondary contact and male-biased dispersal. Our work demonstrates that complicated dispersal-vicariance evolutionary processes likely led to the current geographic distribution of P. glacialis in China, particularly the uplift of the Qinghai-Tibet Plateau and related climatic oscillations during the Quaternary period.

Landscape genetics reveals unique and shared effects of urbanization for two sympatric pool-breeding amphibians

Metapopulation-structured species can be negatively affected when landscape fragmentation impairs connectivity. We investigated the effects of urbanization on genetic diversity and gene flow for two sympatric amphibian species, spotted salamanders (Ambystoma maculatum) and wood frogs (Lithobates sylvaticus), across a large (>35,000 km2) landscape in Maine, USA, containing numerous natural and anthropogenic gradients. Isolation-by-distance (IBD) patterns differed between the species. Spotted salamanders showed a linear and relatively high variance relationship between genetic and geographic distances (r = .057, p < .001), whereas wood frogs exhibited a strongly nonlinear and lower variance relationship (r = 0.429, p < .001). Scale dependence analysis of IBD found gene flow has its most predictable influence (strongest IBD correlations) at distances up to 9 km for spotted salamanders and up to 6 km for wood frogs. Estimated effective migration surfaces revealed contrasting patterns of high and low genetic diversity and gene flow between the two species. Population isolation, quantified as the mean IBD residuals for each population, was associated with local urbanization and less genetic diversity in both species. The influence of geographic proximity and urbanization on population connectivity was further supported by distance-based redundancy analysis and multiple matrix regression with randomization. Resistance surface modeling found interpopulation connectivity to be influenced by developed land cover, light roads, interstates, and topography for both species, plus secondary roads and rivers for wood frogs. Our results highlight the influence of anthropogenic landscape features within the context of natural features and broad spatial genetic patterns, in turn supporting the premise that while urbanization significantly restricts interpopulation connectivity for wood frogs and spotted salamanders, specific landscape elements have unique effects on these two sympatric species.

Using high-throughput sequencing to investigate the factors structuring genomic variation of a Mediterranean grasshopper of great conservation concern

Inferring the demographic history of species is fundamental for understanding their responses to past climate/landscape alterations and improving our predictions about the future impacts of the different components of ongoing global change. Estimating the time-frame at which population fragmentation took place is also critical to determine whether such process was shaped by ancient events (e.g. past climate/geological changes) or if, conversely, it was driven by recent human activities (e.g. habitat loss). We employed genomic data (ddRAD-Seq) to determine the factors shaping contemporary patterns of genetic variation in the Iberian cross-backed grasshopper Dociostaurus crassiusculus, an endangered species with limited dispersal capacity and narrow habitat requirements. Our analyses indicate the presence of two ancient lineages and three genetic clusters resulted from historical processes of population fragmentation (~18-126 ka) that predate the Anthropocene. Landscape genetic analyses indicate that the limits of major river basins are the main geographical feature explaining large-scale patterns of genomic differentiation, with no apparent effect of human-driven habitat fragmentation. Overall, our study highlights the importance of detailed phylogeographic, demographic and spatially-explicit landscape analyses to identify evolutionary significant units and determine the relative impact of historical vs. anthropogenic factors on processes of genetic fragmentation in taxa of great conservation concern.

The influence of latitude, geographic distance, and habitat discontinuities on genetic variation in a high latitude montane species

Examining the factors that influence contemporary genetic patterns is important given the alarming rate at which natural environments are changing. In particular habitat fragmentation and climate change are expected to influence the distribution and diversity of natural populations. In this study we used both mitochondrial control region (mtDNA) and microsatellite data to answer the following questions about genetic diversity and divergence in mountain chickadees (Poecile gambeli) a resident bird species in western North America: (1) Do populations exhibit similar levels of genetic diversity across the range? (2) What is the genetic affinity of western populations in Oregon and Washington? (3) Do genetic patterns exhibit isolation by distance, or are genetic patterns more heavily influenced by habitat discontinuity? We tested the effects of isolation by distance and habitat distribution on genetic structure by analyzing 266 samples from 17 sites across western Canada and the United States. We found a near significant relationship between genetic diversity and latitude, however, our results indicate that overall, latitude is not a strong predictor of genetic diversity. Our analyses of populations in Oregon and Washington revealed a mismatch between patterns detected with mtDNA and microsatellite data. In particular, Washington clustered with the Coast Range/Cascades/Rocky Mountain mtDNA group, but with populations in southern Oregon/California based on microsatellite data. These results suggest the presence of a contact zone in Washington between the two mtDNA clades Coast Range/Cascades/Rocky Mountain and southern Oregon/California clades. Finally, our study revealed a greater effect of isolation by distance than isolation by habitat for both mtDNA and microsatellite data. Overall the isolation by distance signal was greater for mtDNA than microsatellite patterns. The greater signal of isolation by distance on mtDNA patterns likely reflects the strong effects of Pleistocene glaciations in shaping genetic patterns in western North America.

Development and characterization of microsatellite markers for endangered species Stipa pennata (Poaceae) and their usefulness in intraspecific delimitation

Stipa pennata (Poaceae), has become a rare and endangered species in Central Europe due habitat loss and fragmentation. This species is characterized by high morphological variability, which has resulted in the description of numerous intraspecific taxa. The aim of present work is to develop microsatellite markers useful in population genetics studies as well as in intraspecific taxonomy of S. pennata s.l. We developed ten microsatellite markers using Illumina high-throughput. Polymorphism at each marker was evaluated using 4–15 individuals from four morphotypes of S. pennata s.l. Seven markers showed polymorphism while three were monomorphic. The number of alleles per locus ranged from 7 to 12, and the observed and expected heterozygosity varied from 0.000 to 1.000 and 0.000 to 0.8670, respectively. Our results confirm that three of four studied morphotypes are genetically distinct. The microsatellite markers developed here will be useful for evaluating levels of genetic diversity and differentiation, to study gene flow, population dynamics and in future conservation studies as well as for intraspecific delimitation of morphologically similar taxa within S. pennata s.l.

Spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus revealed by museum genomics

Analyzing genetic variation through time and space is important to identify key evolutionary and ecological processes in populations. However, using contemporary genetic data to infer the dynamics of genetic diversity may be at risk of a bias, as inferences are performed from a set of extant populations, setting aside unavailable, rare, or now extinct lineages. Here, we took advantage of new developments in next-generation sequencing to analyze the spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus, a steppic Southwestern-Palearctic species. We applied a recently developed hybridization capture (hyRAD) protocol that allows retrieving orthologous sequences even from degraded DNA characteristic of museum specimens. We identified single nucleotide polymorphisms in 68 historical and 51 modern samples in order to (i) unravel the spatial genetic structure across part of the species distribution and (ii) assess the loss of genetic diversity over the past century in Swiss populations. Our results revealed (i) the presence of three potential glacial refugia spread across the European continent and converging spatially in the Alpine area. In addition, and despite a limited population sample size, our results indicate (ii) a loss of allelic richness in contemporary Swiss populations compared to historical populations, whereas levels of expected heterozygosities were not significantly different. This observation is compatible with an increase in the bottleneck magnitude experienced by central European populations of O. decorus following human-mediated land-use change impacting steppic habitats. Our results confirm that application of hyRAD to museum samples produces valuable information to study genetic processes across time and space.

Inferring the demographic history of an oligophagous grasshopper: Effects of climatic niche stability and host-plant distribution

Understanding the consequences of past environmental changes on the abiotic and biotic components of the landscape and deciphering their impacts on the demographic trajectories of species is a major issue in evolutionary biogeography. In this study, we combine nuclear and mitochondrial genetic data to study the phylogeographical structure and lineage-specific demographic histories of the scrub-legume grasshopper (Chorthippus binotatus binotatus), a montane taxon distributed in the Iberian Peninsula and France that exclusively feeds on certain scrub-legume species. Genetic data and paleo-distribution modelling indicate the presence of four main lineages that seem to have diverged in allopatry and long-term persisted in Iberian and French refugia since the Mid Pleistocene. Comparisons of different demographic hypotheses in an Approximate Bayesian Computation (ABC) framework supported a population bottleneck in the northwestern French clade and paleo-distribution modelling indicate that the populations of this lineage have experienced more severe environmental fluctuations during the last 21 000 years than those from the Iberian Peninsula. Accordingly, we found that nuclear genetic diversity of the populations of scrub-legume grasshopper is positively associated with local stability of suitable habitats defined by both Pleistocene climate changes and historical distributional shifts of host-plant species. Overall, our study highlights the importance of integrating the potential effects of abiotic (i.e. climate and geography) and biotic components (i.e. inter-specific interactions) into the study of the evolutionary and demographic history of specialist taxa with narrow ecological requirements.

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.

Impacts of human-induced environmental disturbances on hybridization between two ecologically differentiated Californian oak species

Natural hybridization, which can be involved in local adaptation and in speciation processes, has been linked to different sources of anthropogenic disturbance. Here, we use genotypic data to study range-wide patterns of genetic admixture between the serpentine-soil specialist leather oak (Quercus durata) and the widespread Californian scrub oak (Quercus berberidifolia). First, we estimated hybridization rates and the direction of gene flow. Second, we tested the hypothesis that genetic admixture increases with different sources of environmental disturbance, namely anthropogenic destruction of natural habitats and wildfire frequency estimated from long-term records of fire occurrence. Our analyses indicate considerable rates of hybridization (> 25%), asymmetric gene flow from Q. durata into Q. berberidifolia, and a higher occurrence of hybrids in areas where both species live in close parapatry. In accordance with the environmental disturbance hypothesis, we found that genetic admixture increases with wildfire frequency, but we did not find a significant effect of other sources of human-induced habitat alteration (urbanization, land clearing for agriculture) or a suite of ecological factors (climate, elevation, soil type). Our findings highlight that wildfires constitute an important source of environmental disturbance, promoting hybridization between two ecologically well-differentiated native species.

High gene flow in the American badger overrides habitat preferences and limits broadscale genetic structure

Habitat associations are a function of habitat preferences and dispersal capabilities, both of which can influence how species responded to Quaternary climatic changes and contemporary habitat heterogeneity. Predicting resultant genetic structure is not always straightforward, especially in species where high dispersal potential and habitat preferences yield opposing predictions. The American badger has high dispersal capabilities that predict widespread panmixia, but avoids closed-canopy forests and clay soils, which could restrict gene flow and create ecologically based population genetic structure. We used mitochondrial sequence and microsatellite data sets to characterize how these opposing forces contribute to genetic structure in badgers at a continent-wide scale. Our data revealed an overall lack of ecologically based population genetic structure, suggesting that high dispersal capabilities were sufficiently realized to overcome most habitat-based genetic structure. At a broadscale, badger gene flow is limited only by geographic distance (isolation by distance) and large water barriers (Lake Michigan and the Mississippi River). The absence of genetic structure in a species with strong avoidance of unsuitable habitats advances our understanding of when and how genetic structure emerges in widespread, highly mobile species.

Hierarchical genetic structure shaped by topography in a narrow-endemic montane grasshopper

BACKGROUND

Understanding the underlying processes shaping spatial patterns of genetic structure in free-ranging organisms is a central topic in evolutionary biology. Here, we aim to disentangle the relative importance of neutral (i.e. genetic drift) and local adaptation (i.e. ecological divergence) processes in the evolution of spatial genetic structure of the Morales grasshopper (Chorthippus saulcyi moralesi), a narrow-endemic taxon restricted to the Central Pyrenees. More specifically, we analysed range-wide patterns of genetic structure and tested whether they were shaped by geography (isolation-by-distance, IBD), topographic complexity and present and past habitat suitability models (isolation-by-resistance, IBR), and environmental dissimilarity (isolation-by-environment, IBE).

RESULTS

Different clustering analyses revealed a deep genetic structure that was best explained by IBR based on topographic complexity. Our analyses did not reveal a significant role of IBE, a fact that may be due to low environmental variation among populations and/or consequence of other ecological factors not considered in this study are involved in local adaptation processes. IBR scenarios informed by current and past climate distribution models did not show either a significant impact on genetic differentiation after controlling for the effects of topographic complexity, which may indicate that they are not capturing well microhabitat structure in the present or the genetic signal left by dispersal routes defined by habitat corridors in the past.

CONCLUSIONS

Overall, these results indicate that spatial patterns of genetic variation in our study system are primarily explained by neutral divergence and migration-drift equilibrium due to limited dispersal across abrupt reliefs, whereas environmental variation or spatial heterogeneity in habitat suitability associated with the complex topography of the region had no significant effect on genetic discontinuities after controlling for geography. Our study highlights the importance of considering a comprehensive suite of potential isolating mechanisms and analytical approaches in order to get robust inferences on the processes promoting genetic divergence of natural populations.