Spatial Genetic Variation Research Articles

Predicting rapid adaptation in time from adaptation in space: A 30-year field experiment in marine snails

Predicting the outcomes of adaptation is a major goal of evolutionary biology. When temporal changes in the environment mirror spatial gradients, it opens up the potential for predicting the course of adaptive evolution over time based on patterns of spatial genetic and phenotypic variation. We assessed this approach in a 30-year transplant experiment in the intertidal snail Littorina saxatilis . In 1992, snails were transplanted from a predation-dominated environment to one dominated by wave action. On the basis of spatial patterns, we predicted transitions in shell size and morphology, allele frequencies at positions throughout the genome, and chromosomal rearrangement frequencies. Observed changes closely agreed with predictions and transformation was both dramatic and rapid. Hence, adaptation can be predicted from knowledge of the phenotypic and genetic variation among populations.

Environmental suitability throughout the late quaternary explains population genetic diversity

Genetic variation among populations is reflected in biogeographic patterns for many species, but general rules of spatial genetic variation have not been established. In this paper, we establish a theoretical framework based on projecting environmental Grinellian niches back through time to relate the present geographic distribution of population genetic structure to a given species' historical evolutionary context. Thanks to advances in next‐generation sequencing technologies, as well as more accurate climate models and the amassing of information stored in biological collections, it is possible to implement this theoretical framework directly. We develop a case study of the tassel‐eared squirrel Sciurus aberti to jointly analyze spatial, environmental, and genetic data to predict the historical endemic area of this species. Our results reveal that in cases of genetic isolation by geographic distance, the prevalence of environmental suitability over time corresponds to the genetic fixation index (Fst) of populations with respect to a source population. Populations closer to the historical endemic area show higher genetic diversity and a lower Fst value. This empirical example relates back to the theoretical framework, allowing two further advances: 1) a layer of biogeographic explanation for the results obtained from population genomic methods; and 2) predictive maps of this genetic structure to support biodiversity conservation efforts. Overall, this work advances a perspective that integrates population genetics with historical patterns of species distribution. The limitations posed in the theoretical framework should be considered before implementing the suitability prevalence area (SPA) in a general way over different taxa. Otherwise, the predictability of the genetic diversity of populations as a product of environmental stability over time may not be adequate.

Oceanographical-driven dispersal and environmental variation explain genetic structure in an upwelling coastal ecosystem

The seascape comprises multiple environmental variables that interact with species biology to determine patterns of spatial genetic variation. The environment imposes spatially variable selective forces together with homogenizing and diverging drivers that facilitate or restrict dispersal, which is a complex, time-dependent process. Understanding how the seascape influences spatial patterns of genetic variation remains elusive, particularly in coastal upwelling systems. Here, we combine genome-wide SNP data, Lagrangian larval dispersal simulated over a hydrodynamic model, and ocean environmental information to quantify the relative contribution of ocean circulation and environmental heterogeneity as drivers of the spatial genetic structure of two congeneric intertidal limpets, Scurria scurra and S. araucana, along the central coast of Chile. We find that a genetic break observed in both limpet species coincides with a break in connectivity shown by the Lagrangian dispersal, suggesting that mean ocean circulation is an important seascape feature, in particular for S. scurra. For S. araucana, environmental variation appears as a better predictor of genetic structure than ocean circulation. Overall, our study shows broad patterns of seascape forcing on genetic diversity and contributes to our understanding of the complex ecological and evolutionary interactions along coastal upwelling systems.

Sex and landscape influence spatial genetic variation in a large fossorial mammal, the Bare-nosed Wombat (Vombatus ursinus).

Dispersal is an important process that is widely studied across species, and it can be influenced by intrinsic and extrinsic factors. Intrinsic factors commonly assessed include the sex and age of individuals, while landscape features are frequently-tested extrinsic factors. Here, we investigated the effects of both sex and landscape composition and configuration on genetic distances among bare-nosed wombats (Vombatus ursinus)-one of the largest fossorial mammals in the world and subject to habitat fragmentation, threats from disease, and human persecution including culling as an agricultural pest. We analyzed a data set comprising 74 Tasmanian individuals (30 males and 44 females), genotyped for 9,064 single-nucleotide polymorphisms. We tested for sex-biased dispersal and the influence of landscape features on genetic distances including land use, water, vegetation, elevation, and topographic ruggedness. We detected significant female-biased dispersal, which may be related to females donating burrows to their offspring due to the energetic cost of excavation, given their large body sizes. Land use, waterbodies, and elevation appeared to be significant landscape predictors of genetic distance. Land use potentially reflects land clearing and persecution over the last 200 years. If our findings based on a limited sample size are valid, retention and restoration of nonanthropogenic landscapes in which wombats can move and burrow may be important for gene flow and maintenance of genetic diversity.

The geno-geo-climate nexus: contributions of geographic and ecological factors in shaping the genomic divergence of two closely related threatened rainforest species of Fontainea Heckel (Euphorbiaceae)

ContextProcesses that shape genomic and ecological divergence can reveal important evolutionary dynamics to inform the conservation of threatened species. Fontainea is a genus of rainforest shrubs and small trees including critically endangered and threatened species restricted to narrow, but complex geographic and ecological regions. Several species of Fontainea are subject to spatially explicit conditions and experience limited intra-specific gene flow, likely generating genetic differentiation and local adaptation.ObjectivesHere, we explored the genetic and ecological mechanisms underlying patterns of diversification in two, closely related threatened Fontainea species. Our aim was to compare spatial patterns of genetic variation between the vulnerable Fontainea australis (Southern Fontainea) and critically endangered F. oraria (Coastal Fontainea), endemic to the heterogeneous subtropical region of central, eastern Australia, where large-scale clearing has severely reduced rainforest habitat to a fraction (< 1%) of its pre-European settlement extent.MethodsWe used a set of 10,000 reduced-representation markers to infer genetic relationships and the drivers of spatial genetic variation across the two species. In addition, we employed a combination of univariate and multivariate genome-environment association analysis using a set of topo-climatic variables to explore potential patterns of local adaptation as a factor impacting genomic divergence.ResultsOur study revealed that Coastal Fontainea have a close genetic relationship with Southern Fontainea. We showed that isolation by distance has played a key role in their genetic variation, indicating that vicariance can explain the spatial genetic distribution of the two species. Genotype-environment analyses showed a strong association with temperature and topographic features, suggesting adaptation to localised thermal environments. We used a multivariate redundancy analysis to identify a range of putatively adapted loci associated with local environmental conditions.ConclusionsDivergent selection at the local-habitat scale as a result of dispersal limitations and environmental heterogeneity (including physical barriers) are likely contributors to adaptive divergence between the two Fontainea species. Our findings have presented evidence to indicate that Southern and Coastal Fontainea were comprised of distinct genetic groups and ecotypes, that together may form a single species continuum, with further phenotype research suggested to confirm the current species boundaries. Proactive conservation actions, including assisted migration to enhance the resilience of populations lacking stress-tolerant single nucleotide polymorphisms (SNPs) may be required to secure the long-term future of both taxa. This is especially vital for the critically endangered Coastal Fontainea given projections of habitat decline for the species under future climate scenarios.

DisperseNN2: a neural network for estimating dispersal distance from georeferenced polymorphism data

Spatial genetic variation is shaped in part by an organism’s dispersal ability. We present a deep learning tool, disperseNN2, for estimating the mean per-generation dispersal distance from georeferenced polymorphism data. Our neural network performs feature extraction on pairs of genotypes, and uses the geographic information that comes with each sample. These attributes led disperseNN2 to outperform a state-of-the-art deep learning method that does not use explicit spatial information: the mean relative absolute error was reduced by 33% and 48% using sample sizes of 10 and 100 individuals, respectively. disperseNN2 is particularly useful for non-model organisms or systems with sparse genomic resources, as it uses unphased, single nucleotide polymorphisms as its input. The software is open source and available from https://github.com/kr-colab/disperseNN2, with documentation located at https://dispersenn2.readthedocs.io/en/latest/.

A metapopulation model reveals connectivity-driven hotspots in treatment resistance evolution in a marine parasite

Abstract In salmon aquaculture, the sustainable management of salmon lice (Lepeophtheirus salmonis) is limited by the adaptive capacity of the parasite. This is evident in the repeated evolution of pesticide resistance in the salmon louse population. To better prepare for resistance, we constructed a numerical metapopulation model that predicts the evolutionary dynamics of lice across an interconnected farm network. This model integrates within-farm population dynamics and between-farm louse dispersal, the latter using outputs from a state-of-the-art particle-tracking model. Distinct from previous metapopulation models, it also simulates spatial and temporal genetic variation arising from selection. The model was parameterized to simulate the evolution of resistance to the pesticide azamethiphos on farms in southern Norway. It successfully reproduced the rapid (within 10 years) evolution of azamethiphos resistance following extensive delousing treatments. It also identified strong spatial patterns in resistance, with regions of high farm connectivity being potential hotspots of louse adaptation. Rates of infestation and evolution were significantly reduced when highly connected farms were excluded from the simulation, compared to when low-connectivity or random sites were excluded. This model can be a valuable tool for coordinating pest management at a regional scale, in a way that slows or prevents the spread of resistance.

Spatial and temporal variation of genetic diversity and genetic differentiation in Daphnia galeata populations in four large reservoirs in southern China.

Daphnia galeata is a common and dominant species in warmer waters, and has a strong top-down effect on both phytoplankton and bacteria. The knowledge of its temporal and spatial patterns of genetic diversity is fundamental in understanding its population dynamics and potential ecological function in ecosystems. Its population genetics have been investigated at regional scales but few within regions or at smaller spatial scales. Here, we examined the fine-scale spatial genetic variation of D. galeata within four large, deep reservoirs in wet and dry seasons and the six-year variation of genetic diversity in one of the reservoirs by using cytochrome c oxidase subunit I and microsatellites (simple sequence repeat). Our study shows that fine-scale spatial genetic variation commonly occurred within the reservoirs, indicating strong environmental selection at least in the two of reservoirs with strong longitudinal gradients. Since the environmental gradients established in the dry season was largely reduced in the wet season, the fine-scale spatial genetic variation was much higher in the dry season. The dynamics of local genetic diversity did not follow the theoretical pattern of rapid erosion but peaked in mid or mid-late growth season. The local genetic diversity of D. galeata appears to be shaped and maintained not only by recruitment from resting egg banks but also by gene flow within reservoirs. The temporal and fine-scale genetic variation within a water body suggests that it is necessary to pay attention to sampling periods and locations of a given water body in regional studies.

Population genetics of a recent range expansion and subsequent loss of migration in monarch butterflies.

Range expansions—whether permanent or transient—strongly influence the distribution of genetic variation in space. Monarch butterflies are best known for long‐distance seasonal migration within North America but are also established as nonmigratory populations around the world, including on Pacific Islands. Previous research has highlighted stepwise expansion across the Pacific, though questions remain about expansion timing and the population genetic consequences of migration loss. Here, we present reduced‐representation sequencing data for 275 monarchs from North America (n = 85), 12 Pacific Islands (n = 136) and three locations in Australia (n = 54), with the goal of understanding (i) how the monarch's Pacific expansion has shaped patterns of population genetic variation and (ii) how loss of migration has influenced spatial patterns of differentiation. We find support for previously described stepwise dispersal across the Pacific and document an additional expansion from Hawaii into the Mariana Islands. Nonmigratory monarchs within the Mariana Islands show strong patterns of differentiation, despite their proximity; by contrast, migratory North American samples form a single genetically panmictic population across the continent. Estimates of Pacific establishment timing are highly uncertain (~100–1,000,000 years ago) but overlap with historical records that indicate a recent expansion. Our data support (i) a recent expansion across the Pacific whose timing overlaps with available historical records of establishment and (ii) a strong role for seasonal migration in determining patterns of spatial genetic variation. Our results are noteworthy because they demonstrate how the evolution of partial migration can drive population differentiation over contemporary timescales.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef.

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A.tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Temporal and spatial genetic variation of Engraulis encrasicolus in the Adriatic Sea

Small pelagic fish play a key role in the marine ecosystem, controlling predator abundance and regulating primary production rates by foraging on plankton. Alterations induced at different ecological levels by fishing activities and/or environmental changes are affecting the reproductive success of several small pelagics, including European anchovy (Engraulis encrasicolus), which is a major target of Adriatic mid-water trawl and purse-seine fisheries. In this study, we evaluated short-term genetic changes of the species in the Adriatic Sea by applying molecular markers in samples of three generations of European anchovy. Thirteen polymorphic microsatellite markers and a mitochondrial gene were used in cohorts of adults and larvae, collected at three sites in the north-east, central, and southern Adriatic Sea in 2015. Furthermore, temporal, and spatial genetic variation was assessed by comparing the above dataset with a sample of adult anchovy collected in 2012 at three sites close to those sampled in 2015. Expected heterozygosity was higher in adults than in larvae, suggesting a loss of genetic diversity and uneven reproduction. In addition, a comparison of the two datasets demonstrated a change in the anchovy population structure from 2012 to 2015. In the reproductive event of 2015, this change led the two main genetic stocks described in the Adriatic Sea to merge into one. We suggest that the population structure of European anchovy in the north-eastern Adriatic may be influenced by changes in environmental parameters and by periodic alterations in the temporal pattern of population connectivity.

Right out of the gate: the genomics of Lessepsian invaders in the vicinity of the Suez Canal

Marine organisms that enter the Mediterranean from the Red Sea via the Suez Canal are known as Lessepsian bioinvaders. In general, genetic studies of Lessepsian fishes have shown little structure between Red Sea and Mediterranean populations. Yet notable exceptions suggest the importance of life-history factors that may influence patterns of spatial genetic variation. In this study, by sampling two invasive fishes with different life histories (the rabbitfish Siganus rivulatus and the filefish Stephanolepis diaspros), we looked at evidence of population structure and selection at the boundary between the Red Sea and the Mediterranean (the Suez Canal), using thousands of molecular markers. Results illustrate two divergent patterns of genetic patterns, with little genetic structure in S. rivulatus and strong population structure in S. diaspros, even at such small spatial scale. We discuss differences in ecological characteristics between the two species to account for such differences. In addition, we report that in the face of both high (S. rivulatus) and low (S. diaspros) gene flow, loci under selection were uncovered, and some protein coding genes were identified as being involved with osmoregulation, which seems to be an important feature of individuals crossing the salinity-variable Suez Canal. The presence of genes under selection in populations near the Suez Canal supports the idea that selection may be active and essential for successful invasions right out of the gate.

Assessment of drivers of spatial genetic variation of a ground-dwelling bird species and its implications for conservation.

In modern wildlife ecology, spatial population genetic methods are becoming increasingly applied. Especially for animal species in fragmented landscapes, preservation of gene flow becomes a high priority target in order to restore genetic diversity and prevent local extinction. Within Central Europe, the Alps represent the core distribution area of the black grouse, Lyrurus tetrix. At its easternmost Alpine range, events of subpopulation extinction have already been documented in the past decades. Molecular data combined with spatial analyses can help to assess landscape effects on genetic variation and therefore can be informative for conservation management. Here, we addressed whether the genetic pattern of the easternmost Alpine black grouse metapopulation system is driven by isolation by distance or isolation by resistance. Correlative ecological niche modeling was used to assess geographic distances and landscape resistances. We then applied regression‐based approaches combined with population genetic analyses based on microsatellite data to disentangle effects of isolation by distance and isolation by resistance among individuals and subpopulations. Although population genetic analyses revealed overall low levels of genetic differentiation, the ecological niche modeling showed subpopulations to be clearly delimited by habitat structures. Spatial genetic variation could be attributed to effects of isolation by distance among individuals and isolation by resistance among subpopulations, yet unknown effects might factor in. The easternmost subpopulation was the most differentiated, and at the same time, immigration was not detected; hence, its long‐term survival might be threatened. Our study provides valuable insights into the spatial genetic variation of this small‐scale metapopulation system of Alpine black grouse.

Spatial and temporal genetic variation in Ethiopian barley (Hordeum vulgare L.) landraces as revealed by simple sequence repeat (SSR) markers

BackgroundEthiopia is considered as center of diversity for barley (Hordeum vulgare L.) and it is grown across different agro-ecologies of the country. Unraveling population structure and gene flow status on temporal scales assists an evaluation of the consequences of physical, demographic and overall environmental changes on the stability and persistence of populations. This study was to examine spatial and temporal genetic variation within and among barley landrace samples collected over a period of four decades, using simple sequence repeat markers.ResultsResults from STRUCTURE, neighbor joining tree and discriminant analysis of principal component (DAPC) analysis revealed presence low-to-high genetic diversity among the landraces and grouped the landraces into three clusters. The cluster analysis revealed a close relationship between landraces along geographic proximity with genetic distance increases along with geographic distance. From analysis of molecular variance (AMOVA) in terms of collection year, it was observed that within-population genetic diversity much higher than between population and that the temporal differentiation is considerably smaller. The low-to-high genetic differentiation between landraces could be attributed to gene flow across the region as a consequence of seed exchange among farmers.ConclusionThe results demonstrate that this set of SSRs was highly informative and useful in generating a meaningful classification of barley germplasms. Furthermore, results obtained from this study also suggest that landraces are a source of valuable germplasm for sustainable agriculture in the context of future climate change and in situ conservation strategies following adaptation to local environments.

Tumors Are Evolutionary Island-Like Ecosystems.

Integration of ecological and evolutionary features has begun to understand the interplay of tumor heterogeneity, microenvironment, and metastatic potential. Developing a theoretical framework is intrinsic to deciphering tumors’ tremendous spatial and longitudinal genetic variation patterns in patients. Here, we propose that tumors can be considered evolutionary island-like ecosystems, that is, isolated systems that undergo evolutionary and spatiotemporal dynamic processes that shape tumor microenvironments and drive the migration of cancer cells. We examine attributes of insular systems and causes of insularity, such as physical distance and connectivity. These properties modulate migration rates of cancer cells through processes causing spatial and temporal isolation of the organs and tissues functioning as a supply of cancer cells for new colonizations. We discuss hypotheses, predictions, and limitations of tumors as islands analogy. We present emerging evidence of tumor insularity in different cancer types and discuss their relevance to the islands model. We suggest that the engagement of tumor insularity into conceptual and mathematical models holds promise to illuminate cancer evolution, tumor heterogeneity, and metastatic potential of cells.

Spatio-temporal patterns of genetic variation of the silverside Odontesthes regia in the highly productive Humboldt Current System

The Humboldt Current System (HCS) is characterised by latitudinal upwelling differences and interannual El Niño events that bring equatorial warm waters into the system. Within the HCS, the silverside Odontesthes regia is a coastal fish that faces intense harvesting, and physical and oceanographic variability. Here, we evaluated how different processes in the HCS have shaped the spatial and temporal genetic variation of O. regia populations by using 650-bp alignment of the mitochondrial control region and nine microsatellite loci. We found that both DNA markers indicated high genetic diversity, presence of at least two co-distributed genetic groups, and a slight genetic structure shaped by a biogeographic break within the HCS rather than the latitudinal upwelling differences or geographical distances. Demographic history of these groups suggests an allopatric origin and posterior reconnection as the result of historical changes in the HCS. No genetic interannual divergences were found in subsystems with intermittent presence of warm waters, suggesting no significant influence of El Niño on the genetic variation of O. regia from 2013 to 2018. Despite the inherent disturbances of the HCS, past and recent estimates of effective population sizes showed high values in most populations. Our results provide essential information for the management of O. regia in the HCS. Our study also helps to understand how the spatio-temporal dynamics of the HCS influences the genetic variation of coastal marine fishes.

Association between spatial genetic variation and potential distribution in tree fern Alsophila gigantea (Cyatheaceae) in Hainan Island, China

Spatial genetic variation involves spatial genetic structure (SGS) and genetic diversity is important genetic features of plants. We first evaluated spatial genetic structure (SGS) and genetic diversity among four populations of Alsophila gigantea from Hainan Island, China, using inter-simple sequence repeat (ISSR) markers. Significant but weak FSGS was found in A. gigantea. High genetic diversity was identified at the species level and the population level. AMOVA analysis revealed a low level of genetic differentiation among the four populations with high gene flow. Mantel test showed no significant correlation between genetic distance and geographic distance. It was found that association between annual mean temperature and annual precipitation with FSGS. Combined with these spatial genetic variation, abundant precipitation and suitable temperature create a stable environment for A. gigantea in Hainan Island, which allows the fern to expand rapidly during the LGM. These results further emphasized the role of outcrossing, and history and environmental factors in the evolution of A. gigantea. This study also provided new insights on in local adaptation of A. gigantea to environmental fluctuations, and available genetic data to enhance the conservation for relict tree ferns.

Coupled Source-Sink Habitats Produce Spatial and Temporal Variation of Cancer Cell Molecular Properties as an Alternative to Branched Clonal Evolution and Stem Cell Paradigms

Intratumoral molecular cancer cell heterogeneity is conventionally ascribed to the accumulation of random mutations that occasionally generate fitter phenotypes. This model is built upon the “mutation-selection” paradigm in which mutations drive ever-fitter cancer cells independent of environmental circumstances. An alternative model posits spatio-temporal variation (e.g., blood flow heterogeneity) drives speciation by selecting for cancer cells adapted to each different environment. Here, spatial genetic variation is the consequence rather than the cause of intratumoral evolution. In nature, spatially heterogenous environments are frequently coupled through migration. Drawing from ecological models, we investigate adjacent well-perfused and poorly-perfused tumor regions as “source” and “sink” habitats, respectively. The source habitat has a high carrying capacity resulting in more emigration than immigration. Sink habitats may support a small (“soft-sink”) or no (“hard-sink”) local population. Ecologically, sink habitats can reduce the population size of the source habitat so that, for example, the density of cancer cells directly around blood vessels may be lower than expected. Evolutionarily, sink habitats can exert a selective pressure favoring traits different from those in the source habitat so that, for example, cancer cells adjacent to blood vessels may be suboptimally adapted for that habitat. Soft sinks favor a generalist cancer cell type that moves between the environment but can, under some circumstances, produce speciation events forming source and sink habitat specialists resulting in significant molecular variation in cancer cells separated by small distances. Finally, sink habitats, with limited blood supply, may receive reduced concentrations of systemic drug treatments; and local hypoxia and acidosis may further decrease drug efficacy allowing cells to survive treatment and evolve resistance. In such cases, the sink transforms into the source habitat for resistant cancer cells, leading to treatment failure and tumor progression. We note these dynamics will result in spatial variations in molecular properties as an alternative to the conventional branched evolution model and will result in cellular migration as well as variation in cancer cell phenotype and proliferation currently described by the stem cell paradigm.

Multiscale patterns of isolation by ecology and fine-scale population structure in Texas bobcats.

Patterns of spatial genetic variation can be generated by a variety of ecological processes, including individual preferences based on habitat. These ecological processes act at multiple spatial and temporal scales, generating scale-dependent effects on gene flow. In this study, we focused on bobcats (Lynx rufus), a highly mobile, generalist felid that exhibits ecological and behavioral plasticity, high abundance, and broad connectivity across much of their range. However, bobcats also show genetic differentiation along habitat breaks, a pattern typically observed in cases of isolation-by-ecology (IBE). The IBE observed in bobcats is hypothesized to occur due to habitat-biased dispersal, but it is unknown if this occurs at other habitat breaks across their range or at what spatial scale IBE becomes most apparent. Thus, we used a multiscale approach to examine isolation by ecology (IBE) patterns in bobcats (Lynx rufus) at both fine and broad spatial scales in western Texas. We genotyped 102 individuals at nine microsatellite loci and used partial redundancy analysis (pRDA) to test if a suite of landscape variables influenced genetic variation in bobcats. Bobcats exhibited a latitudinal cline in population structure with a spatial signature of male-biased dispersal, and no clear barriers to gene flow. Our pRDA tests revealed high genetic similarity in similar habitats, and results differed by spatial scale. At the fine spatial scale, herbaceous rangeland was an important influence on gene flow whereas mixed rangeland and agriculture were significant at the broad spatial scale. Taken together, our results suggests that complex interactions between spatial-use behavior and landscape heterogeneity can create non-random gene flow in highly mobile species like bobcats. Furthermore, our results add to the growing body of data highlighting the importance of multiscale study designs when assessing spatial genetic structure.

Distinct phylogeographic structures and evolutionary histories of wild medicinal Salvia miltiorrhiza Bunge plant populations in the mountains of Central China

Geographic features play an important role in shaping the phylogeographic structures of species. Salvia miltiorrhiza Bunge is a widely studied medicinal crop. In this study, we used three cpDNA sequences (trnL–trnF, atpB–rbcL, and rpl16) to investigate the phylogeographic structure and evolutionary history of 190 S. miltiorrhiza endemic plants taken from 11 locations in the Qinling, Dabie, and Mengshan Mountains of Central China. Fifty-three total haplotypes were identified, and 50 of these were population-specific. In addition, moderate-to-high genetic differentiation was observed among the populations (mean Fst = 0.325), while phylogeographic structure was also observed based on STRUCTURE and spatial analysis of molecular variance (SAMOVA) analyses. In addition, BEAST and DIYABC analyses suggested that divergence times of the haplotypes in the AHS and ALQ2 populations of the Dabie Mountains were earlier than those in the Qinling and Mengshan Mountains. Consequently, an evolutionary scenario of S. miltiorrhiza populations in the Qinling and Mengshan Mountains diverging from populations in the Dabie Mountains was highly supported. The results reported here provide a useful framework for understanding the role of geographic features in shaping the spatial genetic variation and evolutionary history of S. miltiorrhiza plants.