Evolutionary stability, landscape heterogeneity, and human land-usage shape population genetic connectivity in the Cape Floristic Region biodiversity hotspot.

As human-induced change eliminates natural habitats, it impacts genetic diversity and population connectivity for local biodiversity. The South African Cape Floristic Region (CFR) is the most diverse extratropical area for plant biodiversity, and much of its habitat is protected as a UNESCO World Heritage site. There has long been great interest in explaining the underlying factors driving this unique diversity, especially as much of the CFR is endangered by urbanization and other anthropogenic activity. Here, we use a population and landscape genetic analysis of SNP data from the CFR endemic plant Leucadendron salignum or "common sunshine conebush" as a model to address the evolutionary and environmental factors shaping the vast CFR diversity. We found that high population structure, along with relatively deeper and older genealogies, is characteristic of the southwestern CFR, whereas low population structure and more recent lineage coalescence depict the eastern CFR. Population network analyses show genetic connectivity is facilitated in areas of lower elevation and higher seasonal precipitation. These population genetic signatures corroborate CFR species-level patterns consistent with high Pleistocene biome stability and landscape heterogeneity in the southwest, but with coincident instability in the east. Finally, we also find evidence of human land-usage as a significant gene flow barrier, especially in severely threatened lowlands where genetic connectivity has been historically the highest. These results help identify areas where conservation plans can prioritize protecting high genetic diversity threatened by contemporary human activities within this unique cultural UNESCO site.

Evolutionary genomics of gypsy moth populations sampled along a latitudinal gradient.

The European gypsy moth (Lymantria dispar L.) was first introduced to Massachusetts in 1869 and within 150 years has spread throughout eastern North America. This large-scale invasion across a heterogeneous landscape allows examination of the genetic signatures of adaptation potentially associated with rapid geographical spread. We tested the hypothesis that spatially divergent natural selection has driven observed changes in three developmental traits that were measured in a common garden for 165 adult moths sampled from six populations across a latitudinal gradient covering the entirety of the range. We generated genotype data for 91,468 single nucleotide polymorphisms based on double digest restriction-site associated DNA sequencing and used these data to discover genome-wide associations for each trait, as well as to test for signatures of selection on the discovered architectures. Genetic structure across the introduced range of gypsy moth was low in magnitude (FST  = 0.069), with signatures of bottlenecks and spatial expansion apparent in the rare portion of the allele frequency spectrum. Results from applications of Bayesian sparse linear mixed models were consistent with the presumed polygenic architectures of each trait. Further analyses indicated spatially divergent natural selection acting on larval development time and pupal mass, with the linkage disequilibrium component of this test acting as the main driver of observed patterns. The populations most important for these signals were two range-edge populations established less than 30 generations ago. We discuss the importance of rapid polygenic adaptation to the ability of non-native species to invade novel environments.

Population assignment reveals low migratory connectivity in a weakly structured songbird.

Understanding migratory connectivity is essential for determining the drivers behind population dynamics and for implementing effective conservation strategies for migratory species. Genetic markers provide a means to describe migratory connectivity; however, they can be uninformative for species with weak population genetic structure, which has limited their application. Here, we demonstrated a genomic approach to describing migratory connectivity in the prothonotary warbler, Protonotaria citrea, a Neotropical songbird of conservation concern. Using 26,189 single nucleotide polymorphisms (SNPs), we revealed regional genetic structure between the Mississippi River Valley and the Atlantic Seaboard with overall weak genetic differentiation among populations (FST  = 0.0055; 95% CI: 0.0051-0.0059). Genetic variation had a stronger association with geographic rather than environmental factors, with each explaining 14.5% and 8.2% of genetic variation, respectively. By varying the numbers of genomic markers used in population assignment models with individuals of known provenance, we identified a maximum assignment accuracy (89.7% to site, 94.3% to region) using a subset of 600 highly differentiated SNPs. We then assigned samples from nonbreeding sites to breeding region and found low migratory connectivity. Our results highlight the importance of filtering markers for informative loci in models of population assignment. Quantifying migratory connectivity for weakly structured species will be useful for expanding studies to a wider range of migratory species across taxonomic groups and may contribute to a deeper understanding of the evolution of migratory strategies.

Wide outcrossing provides functional connectivity for new and old Banksia populations within a fragmented landscape.

Habitat fragmentation affects landscape connectivity, the extent of which is influenced by the movement capacity of the vectors of seed and pollen dispersal for plants. Negative impacts of reduced connectivity can include reduced fecundity, increased inbreeding, genetic erosion and decreased long-term viability. These are issues for not only old (remnant) populations, but also new (restored) populations. We assessed reproductive and connective functionality within and among remnant and restored populations of a common tree, Banksia menziesii R.Br. (Proteaceae), in a fragmented urban landscape, utilising a genetic and graph theoretical approach. Adult trees and seed cohorts from five remnants and two restored populations were genotyped using microsatellite markers. Genetic variation and pollen dispersal were assessed using direct (paternity assignment) and indirect (pollination graphs and mating system characterisation) methods. Restored populations had greater allelic diversity (Ar = 8.08; 8.34) than remnant populations (Ar range = 6.49-7.41). Genetic differentiation was greater between restored and adjacent remnants (FST = 0.03 and 0.10) than all other pairwise comparisons of remnant populations (mean FST = 0.01 ± 0.01; n = 16 P = 0.001). All populations displayed low correlated paternity (rp = 0.06-0.16) with wide-ranging realised pollen dispersal distances (< 1.7 km) and well-connected pollen networks. Here, we demonstrate reproductive and connective functionality of old and new populations of B. menziesii within a fragmented landscape. Due to long-distance pollination events, the physical size of these sites underestimates their effective population size. Thus, they are functionally equivalent to large populations, integrated into a larger landscape matrix.

Functional connectivity and home range inferred at a microgeographic landscape genetics scale in a desert-dwelling rodent.

Gene flow in animals is limited or facilitated by different features within the landscape matrix they inhabit. The landscape representation in landscape genetics (LG) is traditionally modeled as resistance surfaces (RS), where novel optimization approaches are needed for assigning resistance values that adequately avoid subjectivity. Also, desert ecosystems and mammals are scarcely represented in LG studies. We addressed these issues by evaluating, at a microgeographic scale, the effect of landscape features on functional connectivity of the desert-dwelling Dipodomys merriami. We characterized genetic diversity and structure with microsatellites loci, estimated home ranges and movement of individuals using telemetry-one of the first with rodents, generated a set of individual and composite environmental surfaces based on hypotheses of variables influencing movement, and assessed how these variables relate to individual-based gene flow. Genetic diversity and structure results evidenced a family-induced pattern driven by first-order-related individuals, notably determining landscape genetic inferences. The vegetation cover and soil resistance optimized surface (NDVI) were the best-supported model and a significant predictor of individual genetic distance, followed by humidity and NDVI+humidity. Based on an accurate definition of thematic resolution, we also showed that vegetation is better represented as continuously (vs. categorically) distributed. Hence, with a nonsubjective optimization framework for RS and telemetry, we were able to describe that vegetation cover, soil texture, and climatic variables influence D. merriami's functional connectivity at a microgeographic scale, patterns we could further explain based on the home range, habitat use, and activity observed between sexes. We describe the relationship between environmental features and some aspects of D. merriami's behavior and physiology.

Urban hubs of connectivity: contrasting patterns of gene flow within and among cities in the western black widow spider.

As urbanization drastically alters the natural landscape and generates novel habitats within cities, the potential for changes to gene flow for urban-dwelling species increases. The western black widow spider (Latrodectus hesperus) is a medically relevant urban adapter pest species, for which we have previously identified population genetic signatures consistent with urbanization facilitating gene flow, likely due to human-mediated transport. Here, in an analysis of 1.9 million genome-wide SNPs, we contrast broad-scale geographical analyses of 10 urban and 11 non-urban locales with fine-scale within-city analyses including 30 urban locales across the western USA. These hierarchical datasets enable us to test hypotheses of how urbanization impacts multiple urban cities and their genetic connectivity at different spatial scales. Coupled fine-scale and broad-scale analyses reveal contrasting patterns of high and low genetic differentiation among locales within cities as a result of low and high genetic connectivity, respectively, of these cities to the overall population network. We discuss these results as they challenge the use of cities as replicates of urban eco-evolution, and have implications for conservation and human health in a rapidly growing urban habitat.

Urbanization as a facilitator of gene flow in a human health pest.

Urban fragmentation can reduce gene flow that isolates populations, reduces genetic diversity and increases population differentiation, all of which have negative conservation implications. Alternatively, gene flow may actually be increased among urban areas consistent with an urban facilitation model. In fact, urban adapter pests are able to thrive in the urban environment and may be experiencing human-mediated transport. Here, we used social network theory with a population genetic approach to investigate the impact of urbanization on genetic connectivity in the Western black widow spider, as an urban pest model of human health concern. We collected genomewide single nucleotide polymorphism variation from mitochondrial and nuclear double-digest RAD (ddRAD) sequence data sets from 210 individuals sampled from 11 urban and 10 nonurban locales across its distribution of the Western United States. From urban and nonurban contrasts of population, phylogenetic, and network analyses, urban locales have higher within-population genetic diversity, lower between-population genetic differentiation and higher estimates of genetic connectivity. Social network analyses show that urban locales not only have more connections, but can act as hubs that drive connectivity among nonurban locales, which show signatures of historical isolation. These results are consistent with an urban facilitation model of gene flow and demonstrate the importance of sampling multiple cities and markers to identify the role that urbanization has had on larger spatial scales. As the urban landscape continues to grow, this approach will help determine what factors influence the spread and adaptation of pests, like the venomous black widow spider, in building policies for human and biodiversity health.

Sex-specific graphs: Relating group-specific topology to demographic and landscape data.

Sex-specific genetic structure is a commonly observed pattern among vertebrate species. Facing differential selective pressures, individuals may adopt sex-specific life history traits that ultimately shape genetic variation among populations. Although differential dispersal dynamics are commonly detected in the literature, few studies have used genetic structure to investigate sex-specific functional connectivity. The recent use of graph theoretic approaches in landscape genetics has demonstrated network capacities to describe complex system behaviours where network topology represents genetic interaction among subunits. Here, we partition the overall genetic structure into sex-specific graphs, revealing different male and female dispersal dynamics of a fisher (Pekania [Martes] pennanti) metapopulation in southern Ontario. Our analyses based on network topologies supported the hypothesis of male-biased dispersal. Furthermore, we demonstrated that the effect of the landscape, identified at the population level, could be partitioned among sex-specific strata. We found that female connectivity was negatively correlated with snow depth, whereas connectivity among males was not. Our findings underscore the potential of conducting sex-specific analysis by identifying landscape elements or configuration that differentially promotes or impedes functional connectivity between sexes, revealing processes that may otherwise remain cryptic. We propose that the sex-specific graph approach would be applicable to other vagile species where differential sex-specific processes are expected to occur.

Model Comparison for Abiotic versus Biotic Pollen Dispersal.

An agent-based model with a correlated random walk is used to explore pollination within a forest. For abiotic dispersal, say via the wind, we use a purely random walk where there is no correlation between consecutive steps and for biotic dispersal, say via insect, we use a moderate or highly correlated random walk. In particular, we examine the differences in a number of biological measurement between a purely random walk and a correlated random walk in terms of gene dispersal in low and high plant densities.

Identification of Eastern United States Reticulitermes Termite Species via PCR-RFLP, Assessed Using Training and Test Data.

Reticulitermes termites play key roles in dead wood decomposition and nutrient cycling in forests. They also damage man-made structures, resulting in considerable economic loss. In the eastern United States, five species (R. flavipes, R. virginicus, R. nelsonae, R. hageni and R. malletei) have overlapping ranges and are difficult to distinguish morphologically. Here we present a molecular tool for species identification. It is based on polymerase chain reaction (PCR) amplification of a section of the mitochondrial cytochrome oxidase subunit II gene, followed by a three-enzyme restriction fragment length polymorphism (RFLP) assay, with banding patterns resolved via agarose gel electrophoresis. The assay was designed using a large set of training data obtained from a public DNA sequence database, then evaluated using an independent test panel of Reticulitermes from the Southern Appalachian Mountains, for which species assignments were determined via phylogenetic comparison to reference sequences. After refining the interpretive framework, the PCR-RFLP assay was shown to provide accurate identification of four co-occurring species (the fifth species, R. hageni, was absent from the test panel, so accuracy cannot yet be extended to training data). The assay is cost- and time-efficient, and will help improve knowledge of Reticulitermes species distributions.

Is there such a thing as landscape genetics?

For a scientific discipline to be interdisciplinary, it must satisfy two conditions; it must consist of contributions from at least two existing disciplines, and it must be able to provide insights, through this interaction, that neither progenitor discipline could address. In this study, I examine the complete body of peer-reviewed literature self-identified as landscape genetics (LG) using the statistical approaches of text mining and natural language processing. The goal here was to quantify the kinds of questions being addressed in LG studies, the ways in which questions are evaluated mechanistically, and how they are differentiated from the progenitor disciplines of landscape ecology and population genetics. I then circumscribe the main factions within published LG studies examining the extent to which emergent questions are being addressed and highlighting a deep bifurcation between existing individual- and population-based approaches. I close by providing some suggestions on where theoretical and analytical work is needed if LGs is to serve as a real bridge connecting evolution and ecology sensu lato.

The gravity of pollination: integrating at-site features into spatial analysis of contemporary pollen movement.

Pollen-mediated gene flow is a major driver of spatial genetic structure in plant populations. Both individual plant characteristics and site-specific features of the landscape can modify the perceived attractiveness of plants to their pollinators and thus play an important role in shaping spatial genetic variation. Most studies of landscape-level genetic connectivity in plants have focused on the effects of interindividual distance using spatial and increasingly ecological separation, yet have not incorporated individual plant characteristics or other at-site ecological variables. Using spatially explicit simulations, we first tested the extent to which the inclusion of at-site variables influencing local pollination success improved the statistical characterization of genetic connectivity based upon examination of pollen pool genetic structure. The addition of at-site characteristics provided better models than those that only considered interindividual spatial distance (e.g. IBD). Models parameterized using conditional genetic covariance (e.g. population graphs) also outperformed those assuming panmixia. In a natural population of Cornus florida L. (Cornaceae), we showed that the addition of at-site characteristics (clumping of primary canopy opening above each maternal tree and maternal tree floral output) provided significantly better models describing gene flow than models including only between-site spatial (IBD) and ecological (isolation by resistance) variables. Overall, our results show that including interindividual and local ecological variation greatly aids in characterizing landscape-level measures of contemporary gene flow.

Ecological coassociations influence species' responses to past climatic change: an example from a Sonoran Desert bark beetle.

Ecologically interacting species may have phylogeographical histories that are shaped both by features of their abiotic landscape and by biotic constraints imposed by their coassociation. The Baja California peninsula provides an excellent opportunity to examine the influence of abiotic vs. biotic factors on patterns of diversity in plant-insect species.This is because past climatic and geological changes impacted the genetic structure of plants quite differently to that of codistributed free-living animals (e.g. herpetofauna and small mammals). Thus, 'plant-like' patterns should be discernible in host-specific insect herbivores. Here, we investigate the population history of a monophagous bark beetle, Araptus attenuatus, and consider drivers of phylogeographical patterns in the light of previous work on its host plant, Euphorbia lomelii. Using a combination of phylogenetic, coalescent-simulation-based and exploratory analyses of mitochondrial DNA sequences and nuclear genotypic data, we found that the evolutionary history of A. attenuatus exhibits similarities to its host plant that are attributable to both biotic and abiotic processes. Southward range expansion and recent colonization of continental Sonora from the Baja peninsula appear to be unique to this taxon pair and probably reflect influences of the host plant. On the other hand, abiotic factors with landscape level influences on a diverse suite of codistributed arid-adapted taxa, such as Plio- and Pleistocene-aged marine incursions in the region, also left genetic signatures in beetle and host plant populations. Superimposed on these similarities, bark beetle-specific patterns and processes were also evident: our data revealed two secondarily sympatric,reproductively isolated genetic lineages, as well as a previously unrecognized mid peninsular warm desert refuge. Taken together, this work illustrates that the evolutionary history of species-specific insect herbivores may represent a mosaic of influences,including­but not limited to­those imposed by the host plant.

Defining the landscape of adaptive genetic diversity.

Whether they are used to describe fitness, genome architecture or the spatial distribution of environmental variables, the concept of a landscape has figured prominently in our collective reasoning. The tradition of landscapes in evolutionary biology is one of fitness mapped onto axes defined by phenotypes or molecular sequence states. The characteristics of these landscapes depend on natural selection, which is structured across both genomic and environmental landscapes, and thus, the bridge among differing uses of the landscape concept (i.e. metaphorically or literally) is that of an adaptive phenotype and its distribution across geographical landscapes in relation to selective pressures. One of the ultimate goals of evolutionary biology should thus be to construct fitness landscapes in geographical space. Natural plant populations are ideal systems with which to explore the feasibility of attaining this goal, because much is known about the quantitative genetic architecture of complex traits for many different plant species. What is less known are the molecular components of this architecture. In this issue of Molecular Ecology, Parchman et al. (2012) pioneer one of the first truly genome-wide association studies in a tree that moves us closer to this form of mechanistic understanding for an adaptive phenotype in natural populations of lodgepole pine (Pinus contorta Dougl. ex Loud.).

Landscape modelling of gene flow: improved power using conditional genetic distance derived from the topology of population networks.

Landscape genetics is a burgeoning field of interest that focuses on how site-specific factors influence the distribution of genetic variation and the genetic connectivity of individuals and populations. In this manuscript, we focus on two methodological extensions for landscape genetic analyses: the use of conditional genetic distance (cGD) derived from population networks and the utility of extracting potentially confounding effects caused by correlations between phylogeographic history and contemporary ecological factors. Individual-based simulations show that when describing the spatial distribution of genetic variation, cGD consistently outperforms the traditional genetic distance measure of linearized F(ST) under both 1- and 2-dimensional stepping stone models and Cavalli-Sforza and Edward's chord distance D(c) in 1-dimensional landscapes. To show how to identify and extract the effects of phylogeographic history prior to embarking on landscape genetic analyses, we use nuclear genotypic data from the Sonoran desert succulent Euphorbia lomelii (Euphrobiaceae), for which a detailed phylogeographic history has previously been determined. For E. lomelii, removing the effect of phylogeographic history significantly influences our ability to infer both the identity and the relative importance of spatial and bio-climatic variables in subsequent landscape genetic analyses. We close by discussing the utility of cGD in landscape genetic analyses.

Nuclear gene phylogeography using PHASE: dealing with unresolved genotypes, lost alleles, and systematic bias in parameter estimation.

BACKGROUND: A widely-used approach for screening nuclear DNA markers is to obtain sequence data and use bioinformatic algorithms to estimate which two alleles are present in heterozygous individuals. It is common practice to omit unresolved genotypes from downstream analyses, but the implications of this have not been investigated. We evaluated the haplotype reconstruction method implemented by PHASE in the context of phylogeographic applications. Empirical sequence datasets from five non-coding nuclear loci with gametic phase ascribed by molecular approaches were coupled with simulated datasets to investigate three key issues: (1) haplotype reconstruction error rates and the nature of inference errors, (2) dataset features and genotypic configurations that drive haplotype reconstruction uncertainty, and (3) impacts of omitting unresolved genotypes on levels of observed phylogenetic diversity and the accuracy of downstream phylogeographic analyses. RESULTS: We found that PHASE usually had very low false-positives (i.e., a low rate of confidently inferring haplotype pairs that were incorrect). The majority of genotypes that could not be resolved with high confidence included an allele occurring only once in a dataset, and genotypic configurations involving two low-frequency alleles were disproportionately represented in the pool of unresolved genotypes. The standard practice of omitting unresolved genotypes from downstream analyses can lead to considerable reductions in overall phylogenetic diversity that is skewed towards the loss of alleles with larger-than-average pairwise sequence divergences, and in turn, this causes systematic bias in estimates of important population genetic parameters. CONCLUSIONS: A combination of experimental and computational approaches for resolving phase of segregating sites in phylogeographic applications is essential. We outline practical approaches to mitigating potential impacts of computational haplotype reconstruction on phylogeographic inferences. With targeted application of laboratory procedures that enable unambiguous phase determination via physical isolation of alleles from diploid PCR products, relatively little investment of time and effort is needed to overcome the observed biases.

A systems biology approach to invasive behavior: comparing cancer metastasis and suburban sprawl development.

BACKGROUND: Despite constant progress, cancer remains the second leading cause of death in the United States. The ability of tumors to metastasize is central to this dilemma, as many studies demonstrate successful treatment correlating to diagnosis prior to cancer spread. Hence a better understanding of cancer invasiveness and metastasis could provide critical insight. PRESENTATION OF THE HYPOTHESIS: We hypothesize that a systems biology-based comparison of cancer invasiveness and suburban sprawl will reveal similarities that are instructive. TESTING THE HYPOTHESIS: We compare the structure and behavior of invasive cancer to suburban sprawl development. While these two systems differ vastly in dimension, they appear to adhere to scale-invariant laws consistent with invasive behavior in general. We demonstrate that cancer and sprawl have striking similarities in their natural history, initiating factors, patterns of invasion, vessel distribution and even methods of causing death. IMPLICATIONS OF THE HYPOTHESIS: We propose that metastatic cancer and suburban sprawl provide striking analogs in invasive behavior, to the extent that conclusions from one system could be predictive of behavior in the other. We suggest ways in which this model could be used to advance our understanding of cancer biology and treatment.

GeneticStudio: a suite of programs for spatial analysis of genetic-marker data.

The analysis of genetic marker data is increasingly being conducted in the context of the spatial arrangement of strata (e.g. populations) necessitating a more flexible set of analysis tools. GeneticStudio consists of four interacting programs: (i) Geno a spreadsheet-like interface for the analysis of spatially explicit marker-based genetic variation; (ii) Graph software for the analysis of Population Graph and network topologies, (iii) Manteller, a general purpose for matrix analysis program; and (iv) SNPFinder, a program for identifying single nucleotide polymorphisms. The GeneticStudio suite is available as source code as well as binaries for OSX and Windows and is distributed under the GNU General Public License.

Powers of discerning: challenges to understanding dispersal processes in natural populations.

In this issue of Molecular Ecology, authors Robledo-Arnuncio & Garcia present a compelling approach for quantifying seed dispersal in plant populations. Building upon methods previously used for quantification of pollen dispersal, the authors not only examine the behaviour of the model with respect to sample sizes, dispersal distance, and the kurtosis of the dispersal function but also provide an empirical example using Prunus mahaleb.