Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

A genomic selection component analysis characterizes migration‐selection balance

The genetic differentiation of populations in response to local selection pressures has long been studied by evolutionary biologists, but key details about the process remain obscure. How rapidly can local adaptation evolve, how extensive is the process across the genome, and how strong are the opposing forces of natural selection and gene flow? Here, we combine direct measurement of survival and reproduction with whole‐genome genotyping of a plant species (Mimulus guttatus) that has recently invaded a novel habitat (the Quarry population). We renovate the classic selection component method to accommodate genomic data and observe selection at SNPs throughout the genome. SNPs showing viability selection in Quarry exhibit elevated divergence from neighboring populations relative to neutral SNPs. We also find that nonsignificant SNPs exhibit a subtle, but still significant, change in allele frequency toward neighboring populations, a predicted effect of gene flow. Given that the Quarry population is most probably only 30–40 generations old, the alleles conferring local advantage are almost certainly older than the population itself. Thus, local adaptation owes to the recruitment of standing genetic variation.     

Population dynamics and evolutionary processes: the manifold roles of habitat selection

Summary Any character that has a substantial effect on a species' distribution and abundance can exert a variety of indirect effects on evolutionary processes. It is suggested that an organism's capacity for habitat selection is just such a character. Habitat selection can constrain the selective environment experienced by a population. Habitat selection can also indirectly influence the relative importance of natural selection, drift, and gene flow, through its effect on population size and growth rate. In many circumstances (but not all), habitat selection increases population size and growth rate, and thereby makes selection in a local environment more effective than drift and gene flow.     

Local adaptation and ecological differentiation under selection,migration, and drift in Arabidopsis lyrata*

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole‐genome‐based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4N_em < 11) and have small effective population sizes (N_e < 6000), suggesting that both migration and drift can counteract local selection. However, the three‐year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.     

MAINTENANCE OF ECOLOGICALLY SIGNIFICANT GENETIC VARIATION IN THE TIGER SWALLOWTAIL BUTTERFLY THROUGH DIFFERENTIAL SELECTION AND GENE FLOW

Differential selection in a heterogeneous environment is thought to promote the maintenance of ecologically significant genetic variation. Variation is maintained when selection is counterbalanced by the homogenizing effects of gene flow and random mating. In this study, we examine the relative importance of differential selection and gene flow in maintaining genetic variation in Papilio glaucus. Differential selection on traits contributing to successful use of host plants (oviposition preference and larval performance) was assessed by comparing the responses of southern Ohio, north central Georgia, and southern Florida populations of P. glaucus to three hosts: Liriodendron tulipifera, Magnolia virginiana, and Prunus serotina. Gene flow among populations was estimated using allozyme frequencies from nine polymorphic loci. Significant genetic differentiation was observed among populations for both oviposition preference and larval performance. This differentiation was interpreted to be the result of selection acting on Florida P. glaucus for enhanced use of Magnolia, the prevalent host in Florida. In contrast, no evidence of population differentiation was revealed by allozyme frequencies. F_ST-values were very small and Nm, an estimate of the relative strengths of gene flow and genetic drift, was large, indicating that genetic exchange among P. glaucus populations is relatively unrestricted. The contrasting patterns of spatial differentiation for host-use traits and lack of differentiation for electrophoretically detectable variation implies that differential selection among populations will be counterbalanced by gene flow, thereby maintaining genetic variation for host-use traits.     

Ecological adaptation drives wood frog population divergence in life history traits

Phenotypic variation among populations is thought to be generated from spatial heterogeneity in environments that exert selection pressures that overcome the effects of gene flow and genetic drift. Here, we tested for evidence of isolation by distance or by ecology (i.e., ecological adaptation) to generate variation in early life history traits and phenotypic plasticity among 13 wood frog populations spanning 1200 km and 7° latitude. We conducted a common garden experiment and related trait variation to an ecological gradient derived from an ecological niche model (ENM) validated to account for population density variation. Shorter larval periods, smaller body weight, and relative leg lengths were exhibited by populations with colder mean annual temperatures, greater precipitation, and less seasonality in precipitation and higher population density (high-suitability ENM values). After accounting for neutral genetic variation, the Q_ST–F_ST analysis supported ecological selection as the key process generating population divergence. Further, the relationship between ecology and traits was dependent upon larval density. Specifically, high-suitability/high-density populations in the northern part of the range were better at coping with greater conspecific competition, evidenced by greater postmetamorphic survival and no difference in body weight when reared under stressful conditions of high larval density. Our results support that both climate and competition selection pressures drive clinal variation in larval and metamorphic traits in this species. Range-wide studies like this one are essential for accurate predictions of population’s responses to ongoing ecological change.Subject terms: Biogeography, Ecological modelling, Evolutionary ecology, Evolution     

Signatures of polygenic adaptation associated with climate across the range of a threatened fish species with high genetic connectivity

Adaptive differences across species’ ranges can have important implications for population persistence and conservation management decisions. Despite advances in genomic technologies, detecting adaptive variation in natural populations remains challenging. Key challenges in gene–environment association studies involve distinguishing the effects of drift from those of selection and identifying subtle signatures of polygenic adaptation. We used paired‐end restriction site‐associated DNA sequencing data (6,605 biallelic single nucleotide polymorphisms; SNPs) to examine population structure and test for signatures of adaptation across the geographic range of an iconic Australian endemic freshwater fish species, the Murray cod Maccullochella peelii. Two univariate gene–association methods identified 61 genomic regions associated with climate variation. We also tested for subtle signatures of polygenic adaptation using a multivariate method (redundancy analysis; RDA). The RDA analysis suggested that climate (temperature‐ and precipitation‐related variables) and geography had similar magnitudes of effect in shaping the distribution of SNP genotypes across the sampled range of Murray cod. Although there was poor agreement among the candidate SNPs identified by the univariate methods, the top 5% of SNPs contributing to significant RDA axes included 67% of the SNPs identified by univariate methods. We discuss the potential implications of our findings for the management of Murray cod and other species generally, particularly in relation to informing conservation actions such as translocations to improve evolutionary resilience of natural populations. Our results highlight the value of using a combination of different approaches, including polygenic methods, when testing for signatures of adaptation in landscape genomic studies.     

Population fragmentation and major histocompatibility complex variation in the spotted suslik,Spermophilus suslicus

The fragmentation of populations typically enhances depletion of genetic variation, but highly polymorphic major histocompatibility complex (MHC) genes are thought to be under balancing selection and therefore retain polymorphism despite population bottlenecks. In this study, we investigate MHC DRB (class II) exon 2 variation in 14 spotted suslik populations from two regions differing in their degree of habitat fragmentation and gene flow. We found 16 alleles that segregated in a sample of 248 individuals. The alleles were highly divergent and revealed the hallmark signs of positive selection acting on them in the past, showing a significant excess of nonsynonymous substitutions. This excess was concentrated in putative antigen‐binding sites, which suggests that past selection was driven by pathogens. MHC diversity was significantly lower in fragmented western populations than in the eastern populations, characterized by significant gene flow. In contrast to neutral variation, amova did not reveal genetic differentiation between the two regions. This may indicate similar selective pressures shaping MHC variation in both regions until the recent past. However, MHC allelic richness within a population was correlated with that for microsatellites. F_ST outlier analyses have shown that population differentiation at DRB was neither higher nor lower than expected under neutrality. The results suggest that selection on MHC is not strong enough to counteract drift that results from recent fragmentation of spotted suslik populations.     

Do genetic structure and landscape heterogeneity impact color morph frequency in a polymorphic salamander?

Landscape heterogeneity plays an important role in population structure and divergence, particularly for species with limited vagility. Here, we used a landscape genetic approach to identify how landscape and environmental variables affect genetic structure and color morph frequency in a polymorphic salamander. The eastern red‐backed salamander, Plethodon cinereus, is widely distributed in northeastern North America and contains two common color morphs, striped and unstriped, that are divergent in ecology, behavior, and physiology. To quantify population structure, rates of gene flow, and genetic drift, we amplified 10 microsatellite loci from 648 individuals across 28 sampling localities. This study was conducted in northern Ohio, where populations of P. cinereus exhibit an unusually wide range of morph frequency variation. To test whether genetic distance was more correlated with morph frequency, elevation, canopy cover, waterways, ecological niche or geographic distance, we used resistance distance and least cost path analyses. We then examined whether landscape and environmental variables, genetic distance or geographic distance were correlated with variation in morph frequency. Tests for population structure revealed three genetic clusters across our sampling range, with one cluster monomorphic for the striped morph. Rates of gene flow and genetic drift were low to moderate across sites. Genetic distance was most correlated with ecological niche, elevation and a combination of landscape and environmental variables. In contrast, morph frequency variation was correlated with waterways and geographic distance. Thus, our results suggest that selection is also an important evolutionary force across our sites, and a balance between gene flow, genetic drift and selection interact to maintain the two color morphs.     

Landscape genomic signatures indicate reduced gene flow and forest‐associated adaptive divergence in an endangered neotropical turtle

Human‐induced transformations of ecosystems usually result in fragmented populations subject to increased extinction risk. Fragmentation is also often associated with novel environmental heterogeneity, which in combination with restricted gene flow may increase the opportunity for local adaptation. To manage at‐risk populations in these landscapes, it is important to understand how gene flow is changing, and how populations respond to habitat loss. We conducted a landscape genomics analysis using Restriction‐site Associated DNA sequencing to investigate the evolutionary response of the critically endangered Dahl's Toad‐headed turtle (Mesoclemmys dahli) to severe habitat modification. The species has lost almost all of its natural habitat in the southwestern part of its range and about 70% in the northeast. Based on least cost path analysis across different resistance surfaces for 3,211 SNPs, we found that the landscape matrix is restricting gene flow, causing the fragmentation of the species into at least six populations. Genome scans and allele‐environment association analyses indicate that the population fragments in the deforested grasslands of the southwest are adaptively different from those in the more forested northeast. Populations in areas with no forest had low levels of adaptive genetic diversity and the fixation of ancestrally‐polymorphic SNPs, consistent with directional selection in this novel environment. Our results suggest that this forest‐stream specialist is adapting to pond‐grassland conditions, but it is also suffering from negative consequences of habitat loss, including genetic erosion, isolation, small effective population sizes, and inbreeding. We recommend gene flow restoration via genetic rescue to counteract these threats, and provide guidance for this strategy.     

Genetic and geographic differentiation in the Rio Negro tuco-tuco (Ctenomys rionegrensis): inferring the roles of migration and drift from multiple genetic markers

Among tuco-tucos, Ctenomys rionegrensis is especially amenable to the study of the forces driving population differentiation because of the restricted geographic range it occupies in Uruguay. Within this limited area, the Rio Negro tuco-tuco is limited to sandy soils. It nonetheless exhibits remarkable variation in pelage color, including melanic, agouti, and dark-backed individuals. Two hypotheses have been put forth to explain this pattern: (1) local differentiation and fixation of alternative pelage types by genetic drift under limited gene flow; or (2) fixation by natural selection that may take place even in the presence of gene flow. A previous allozyme study rejected the genetic drift hypothesis on the basis of high inferred levels of migration. New estimates of gene flow from microsatellites and mitochondrial cytochrome b sequences were obtained for C. rionegrensis populations to further test these hypotheses. Much lower levels of gene flow were estimated with these more sensitive markers. Microsatellite-based estimates of gene flow are close to zero and may come closest to estimating current levels of migration. A lack of equilibrium between migration and genetic drift is also strongly suggested by the absence of an isolation-by-distance pattern found in all three genetic datasets. The microsatellite genotype data show that the species is strongly structured geographically, with subpopulations constituting distinct genetic entities. If current levels of gene flow are very low, as indicated by the new data, the local fixation of alternative alleles, including those responsible for pelage color polymorphism, is possible by drift alone. A scenario is thus proposed in which the species expanded in the recent past from a more restricted geographic range and has subsequently differentiated in near isolation, with genetic drift possibly playing a primary role in overall genetic differentiation. The local fixation of pelage color types could also be due to drift, but selection on this trait cannot be ruled out without direct analysis.     

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.     

Contrasting patterns of quantitative and neutral genetic variation in locally adapted populations of the natterjack toad, Bufo calamita

The relative importance of natural selection and genetic drift in determining patterns of phenotypic diversity observed in nature is still unclear. The natterjack toad (Bufo calamita) is one of a few amphibian species capable of breeding in saline ponds, even though water salinity represents a considerable stress for them. Results from two common-garden experiments showed a pattern of geographic variation in embryonic salinity tolerance among populations from either fresh or brackish environments, consistent with the hypothesis of local adaptation. Full-sib analysis showed increased variation in survival among sibships within population for all populations as osmotic stress was increased (broad-sense heritability increased as salinity raised). Nevertheless, toads native to the brackish water environment had the highest overall survival under brackish conditions. Levels of population genetic differentiation for salinity tolerance were higher than those of neutral genetic differentiation, the latter obtained through the analysis of eight microsatellite loci. Microsatellite markers also revealed little population differentiation, lack of an isolation-by-distance pattern, and moderate gene flow connecting the populations. Therefore, environmental stress tolerance appears to have evolved in absence of geographic isolation, and consequently we reject the null hypothesis of neutral differentiation.     

Genetic patchiness of the shore crab Pachygrapsus marmoratus along the Portuguese coast

The population genetic structure of the shore crab Pachygrapsus marmoratus was studied along the Portuguese coast based on six variable microsatellite loci. Genetic differentiation among populations according to a geographic gradient was not detected. This lack of genetic structure reflects the continuous distribution of the species along the Portuguese coast and suggests that gene flow occurs within the studied distribution range. Gene flow is probably maintained by the planktonic larvae of P. marmoratus that can last up to 31 days in the plankton. Tests for population differentiation demonstrated that the Praia das Avencas population is genetically more separated from all other populations, and Bayesian methodologies tend to form 4 groups that clustered together populations that are several hundred kilometres apart. This grouping pattern could be due to coastal hydrological events that are apparently influencing larval flux. Other hypotheses to explain the significant genetic heterogeneity among populations on a local scale and the absence of geographic variation are pre- and post-settlement natural selection events. Results suggest that the forces causing genetic differentiation may be acting on a local scale and that the larval pool is possibly not always mixed homogeneously.     

Environmental Heterogeneity Explains the Genetic Structure of Continental and Mediterranean Populations of Fraxinus angustifolia Vahl

Tree species with wide distributions often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in a wind-pollinated Mediterranean tree species, Fraxinus angustifolia Vahl, within a recognised glacial refugium in Croatia. We sampled 11 populations from environmentally divergent habitats within the Continental and Mediterranean biogeographical regions. We combined genetic data analyses based on nuclear microsatellite loci, multivariate statistics on environmental data and ecological niche modelling (ENM). We identified a geographic structure with a high genetic diversity and low differentiation in the Continental region, which contrasted with the significantly lower genetic diversity and higher population divergence in the Mediterranean region. The positive and significant correlation between environmental and genetic distances after controlling for geographic distance suggests an important influence of ecological divergence of the sites in shaping genetic variation. The ENM provided support for niche differentiation between the populations from the Continental and Mediterranean regions, suggesting that contemporary populations may represent two divergent ecotypes. Ecotype differentiation was also supported by multivariate environmental and genetic distance analyses. Our results suggest that despite extensive gene flow in continental areas, long-term stability of heterogeneous environments have likely promoted genetic divergence of ashes in this region and can explain the present-day genetic variation patterns of these ancient populations.     

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80‐year period. Screening of >1000 gene‐associated single‐nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.     

Simulating natural selection in landscape genetics

Linking landscape effects to key evolutionary processes through individual organism movement and natural selection is essential to provide a foundation for evolutionary landscape genetics. Of particular importance is determining how spatially-explicit, individual-based models differ from classic population genetics and evolutionary ecology models based on ideal panmictic populations in an allopatric setting in their predictions of population structure and frequency of fixation of adaptive alleles. We explore initial applications of a spatially-explicit, individual-based evolutionary landscape genetics program that incorporates all factors--mutation, gene flow, genetic drift and selection--that affect the frequency of an allele in a population. We incorporate natural selection by imposing differential survival rates defined by local relative fitness values on a landscape. Selection coefficients thus can vary not only for genotypes, but also in space as functions of local environmental variability. This simulator enables coupling of gene flow (governed by resistance surfaces), with natural selection (governed by selection surfaces). We validate the individual-based simulations under Wright-Fisher assumptions. We show that under isolation-by-distance processes, there are deviations in the rate of change and equilibrium values of allele frequency. The program provides a valuable tool (cdpop v1.0; http://cel.dbs.umt.edu/software/CDPOP/) for the study of evolutionary landscape genetics that allows explicit evaluation of the interactions between gene flow and selection in complex landscapes.     

The genomic footprint of climate adaptation in Chironomus riparius

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool‐Seq data and population genetic modelling. Common‐garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an F_ST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.     

Genetic connectivity and diversity in inselberg populations of Acacia woodmaniorum,a rare endemic of the Yilgarn Craton banded iron formations

Historically rare plant species with disjunct population distributions and small population sizes might be expected to show significant genetic structure and low levels of genetic diversity because of the effects of inbreeding and genetic drift. Across the globe, terrestrial inselbergs are habitat for rich, often rare and endemic flora and are valuable systems for investigating evolutionary processes that shape patterns of genetic structure and levels of genetic diversity at the landscape scale. We assessed genetic structure and levels of genetic diversity across the range of the historically rare inselberg endemic Acacia woodmaniorum. Phylogeographic and genetic structure indicates that connectivity is not sufficient to produce a panmictic population across the limited geographic range of the species. However, historical levels of gene flow are sufficient to maintain a high degree of adaptive connectivity across the landscape. Genetic diversity indicates gene flow is sufficient to largely counteract any negative genetic effects of inbreeding and random genetic drift in even the most disjunct or smallest populations. Phylogeographic and genetic structure, a signal of isolation by distance and a lack of evidence of recent genetic bottlenecks suggest long-term stability of contemporary population distributions and population sizes. There is some evidence that genetic connectivity among disjunct outcrops may be facilitated by the occasional long distance dispersal of Acacia polyads carried by insect pollinators moved by prevailing winds.