Population genomics and the causes of local differentiation

Exactly 50 years ago, a revolution in empirical population genetics began with the introduction of methods for detecting allelic variation using protein electrophoresis (Throckmorton 1962; Hubby 1963; Lewontin & Hubby 1966). These pioneering scientists showed that populations are chock‐full of genetic variation. This variation was a surprise that required a re‐thinking of evolutionary genetic heuristics. Understanding the causes for the maintenance of this variation became and remains a major area of research. In the process of addressing the causes, this same group of scientists documented geographical genetic structure (Prakash et al. 1969), spawning the continued accumulation of what is now a huge case study catalogue of geographical differentiation (e.g. Loveless & Hamrick 1984; Linhart & Grant 1996). Geographical differentiation is clearly quite common. Yet, a truly general understanding of the patterns in and causes of spatial genetic structure across the genome remains elusive. To what extent is spatial structure driven by drift and phylogeography vs. geographical differences in environmental sources of selection? What proportion of the genome participates? A general understanding requires range‐wide data on spatial patterning of variation across the entire genome. In this issue of Molecular Ecology, Lasky et al. (2012) make important strides towards addressing these issues, taking advantage of three contemporary revolutions in evolutionary biology. Two are technological: high‐throughput sequencing and burgeoning computational power. One is cultural: open access to data from the community of scientists and especially data sets that result from large collaborative efforts. Together, these developments may at last put answers within reach.     

Comparative genomics suggests local adaptations in the invasive small hive beetle

Invasive species are a major driver of ecological and environmental changes that affect human health, food security, and natural biodiversity. The success and impact of biological invasions depend on adaptations to novel abiotic and biotic selective pressures. However, the molecular mechanisms underlying adaptations in invasive parasitic species are inadequately understood. Small hive beetles, Aethina tumida, are parasites of bee nests. Originally endemic to sub‐Saharan Africa, they are now found nearly globally. Here, we investigated the molecular bases of the adaptations to novel environments underlying their invasion routes. Genomes of historic and recent adults A. tumida from both the endemic and introduced ranges were compared. Analysis of gene–environment association identified 3049 candidate loci located in 874 genes. Functional annotation showed a significant bias toward genes linked to growth and reproduction. One of the genes from the apoptosis pathway encodes an “ecdysone‐related protein,” which is a crucial regulator in controlling body size in response to environmental cues for holometabolous insects during cell death and renewal. Genes whose proteins regulate organ size, ovary activation, and oviposition were also detected. Functions of these enriched pathways parallel behavioral differences between introduced and native A. tumida populations, which may reflect patterns of local adaptation. The results considerably improve our understanding of the underlying mechanisms and ecological factors driving adaptations of invasive species. Deep functional investigation of these identified loci will help clarify the mechanisms of local adaptation in A. tumida.     

Population genomics reveals multiple drivers of population differentiation in a sex‐role‐reversed pipefish

A major goal of molecular ecology is to identify the causes of genetic and phenotypic differentiation among populations. Population genomics is suitably poised to tackle these key questions by diagnosing the evolutionary mechanisms driving divergence in nature. Here, we set out to investigate the evolutionary processes underlying population differentiation in the Gulf pipefish, Syngnathus scovelli. We sampled approximately 50 fish from each of 12 populations distributed from the Gulf coast of Texas to the Atlantic coast of Florida and performed restriction‐site‐associated DNA sequencing to identify SNPs throughout the genome. After imposing quality and stringency filters, we selected a panel of 6348 SNPs present in all 12 populations, 1753 of which were not physically linked. We identified a genome‐wide pattern of isolation by distance, in addition to a more substantial genetic break separating populations in the Gulf of Mexico from those in the Atlantic. We also used several divergence outlier approaches and tests for genotype–environment correlations to identify 400 SNPs putatively involved in local adaptation. Patterns of phenotypic differentiation and variation diverged from the overall genomic pattern, suggesting that selection, phenotypic plasticity or demographic factors may be shaping phenotypes in distinct populations. Overall, our results suggest that population divergence is driven by a variety of factors in S. scovelli, including neutral processes and selection on multiple traits.     

Extensive contiguous north–south range expansion of the original population of an invasive lizard in Florida

Florida has one of the two worst non-indigenous invasive species problems in the United States, and all such species are considered by statute to be a potential ecological problem. Unfortunately, little information is available about most of Florida's invasive species. Here, we provide information on range expansion of a population of one of the lesser-known species. The Palm Beach County, Florida, population of the northern curly-tailed lizard, Leiocephalus carinatus armouri, was examined for north–south expansion of its previously documented range. Observations were conducted as contiguous expansion outward from the previously known range. We found a substantial contiguous range expansion by a minimum of 46.3 km to the south and 34.1 km to the north. This species appears to be closely associated with coastal habitat degradation in the form of human infrastructure, especially ageing pavement, sea walls, buildings, and bridge bases.     

Population genomics identifies the origin and signatures of selection of Korean weedy rice

Weedy rice is the same biological species as cultivated rice (Oryza sativa); it is also a noxious weed infesting rice fields worldwide. Its formation and population‐selective or ‐adaptive signatures are poorly understood. In this study, we investigated the phylogenetics, population structure and signatures of selection of Korean weedy rice by determining the whole genomes of 30 weedy rice, 30 landrace rice and ten wild rice samples. The phylogenetic tree and results of ancestry inference study clearly showed that the genetic distance of Korean weedy rice was far from the wild rice and near with cultivated rice. Furthermore, 537 genes showed evidence of recent positive or divergent selection, consistent with some adaptive traits. This study indicates that Korean weedy rice originated from hybridization of modern indica/indica or japonica/japonica rather than wild rice. Moreover, weedy rice is not only a notorious weed in rice fields, but also contains many untapped valuable traits or haplotypes that may be a useful genetic resource for improving cultivated rice.     

Population genetic analysis of a recent range expansion: mechanisms regulating the poleward range limit in the volcano barnacle Tetraclita rubescens

As range shifts coincident with climate change have become increasingly well documented, efforts to describe the causes of range boundaries have increased. Three mechanisms—genetic impoverishment, migration load, or a physical barrier to dispersal—are well described theoretically, but the data needed to distinguish among them have rarely been collected. We describe the distribution, abundance, genetic variation, and environment of Tetraclita rubescens, an intertidal barnacle that expanded its northern range limit by several hundreds of kilometres from San Francisco, CA, USA, since the 1970s. We compare geographic variation in abundance with abiotic and biotic patterns, including sea surface temperatures and the distributions of 387 co‐occurring species, and describe genetic variation in cytochrome c oxidase subunit I, mitochondrial noncoding region, and nine microsatellite loci from 27 locations between Bahia Magdalena (California Baja Sur, Mexico) and Cape Mendocino (CA, USA). We find very high gene flow, high genetic diversity, and a gradient in physical environmental variation coincident with the range limit. We infer that the primary cause of the northern range boundary in T. rubescens is migration load arising from flow of maladapted alleles into peripheral locations and that environmental change, which could have reduced selection against genotypes immigrating into the newly colonized portion of the range, is the most likely cause of the observed range expansion. Because environmental change could similarly affect all taxa in a region whose distributional limits are established by migration load, these mechanisms may be common causes of range boundaries and largely synchronous multi‐species range expansions.     

Genetic connectivity of lionfish (Pterois volitans) in marine protected areas of the Gulf of Mexico and Caribbean Sea

Lionfish (Pterois volitans) have rapidly invaded the tropical Atlantic and spread across the wider Caribbean in a relatively short period of time. Because of its high invasion capacity, we used it as a model to identify the connectivity among nine marine protected areas (MPAs) situated in four countries in the Gulf of Mexico and the Caribbean Sea. This study provides evidence of local genetic differentiation of P. volitans in the Gulf of Mexico and the Caribbean Sea. A total of 475 lionfish samples were characterized with 12 microsatellites, with 6–20 alleles per locus. Departures from Hardy–Weinberg equilibrium (HWE) were found in 10 of the 12 loci, all caused by heterozygous excess. Moderate genetic differentiation was observed between Chiriviche, Venezuela and Xcalak, México localities (F_ST = 0.012), and between the Los Roques and the Veracruz (F_ST = 0.074) sites. STRUCTURE analysis found that four genetic entities best fit our data. A unique genetic group in the Gulf of Mexico may imply that the lionfish invasion unfolded both in a counterclockwise manner in the Gulf of Mexico. In spite of the notable dispersion of P. volitans, our results show some genetic structure, as do other noninvasive Caribbean fish species, suggesting that the connectivity in some MPAs analyzed in the Caribbean is limited and caused by only a few source individuals with subsequent genetic drift leading to local genetic differentiation. This indicates that P. volitans dispersion could be caused by mesoscale phenomena, which produce stochastic connectivity pulses. Due to the isolation of some MPAs from others, these findings may hold a promise for local short‐term control of by means of intensive fishing, even in MPAs, and may have regional long‐term effects.     

Gene flow in the anemone Anthopleura elegantissima limits signatures of local adaptation across an extensive geographic range

Species inhabiting marine environments face a wide range of environmental conditions that vary spatially across several orders of magnitude. The selective pressures that these conditions impose on marine organisms, in combination with potentially high rates of gene flow between distant populations, make it difficult to predict the extent to which these populations can locally adapt. Here, I identify how selection and gene flow influence the population genetic structure of the anemone Anthopleura elegantissima along the Pacific coast of North America. Isolation by distance is the dominant pattern across the range of this species, with a genetic break near Pt. Conception, CA. Furthermore, demographic modelling suggests that this species was historically confined to southerly latitudes before expanding northward. Outlier analyses identify 24 loci under selection (out of ~1,100), but the same analysis on simulated genetic data generated using the most likely demographic model erroneously identified the same number of loci under selection, if not more. Taken together, these results suggest that demographic processes are the dominant force shaping population genetic patterns in A. elegantissima along the Pacific coast of North America. I discuss these patterns in terms of the evolutionary history of A. elegantissima, the potential for local adaptation, and their consequences with respect to interactions with the endosymbiont Breviolum muscatinei across their geographic range.     

Changes in composition of cuticular biochemicals of the facultatively polygynous ant Petalomyrmex phylax during range expansion in Cameroon with respect to social, spatial and genetic variation

In social insects, biochemicals found at the surface of the cuticle are involved in the recognition process and in protection against desiccation and pathogens. However, the relative contribution of evolutionary forces in shaping diversity of these biochemicals remains largely unresolved in ants. We determined the composition of epicuticular biochemicals for workers sampled in 12 populations of the ant Petalomyrmex phylax from Cameroon. Genetic variation at 12 microsatellite markers was used to infer population history and to provide null expectations under the neutrality hypothesis. Genetic data suggest a recent southward range expansion of this ant species. Furthermore, there is a decline southward in the numbers of queens present in mature colonies. Here, we contrast the pattern of biochemical variation against genetic, social and spatial parameters. We thus provide the first estimates of the relative contribution of neutral and selective processes on variation of ant cuticular profile. Populations in migration-drift disequilibrium showed reduction of within-population variation for genetic markers as well as for cuticular profiles. In these populations, the cuticular profile became biased towards a limited number of high molecular weight molecules. Within- and among-population biochemical variation was explained by both genetic and social variation and by the spatial distribution of populations. We therefore propose that during range expansion of P. phylax, the composition of epicuticular compounds has been affected by a combination of neutral processes - genetic drift and spatially limited dispersal - and spatially varying selection, social organization and environmental effects.     

A population founded by a single pair of individuals: establishment,expansion, and evolution

Events occurring at the founding of a population, and in the next few generations, are potentially of great importance for the future evolution of the population. This study reports demographic, genetic, and morphological changes that took place during and after the colonization of the small Galápagos island of Daphne Major by three male and two female large ground finches, Geospiza magnirostris, at the end of 1982. Using assignment tests with microsatellite DNA data we demonstrate heterogeneity among the immigrants. Their sources included both a near island (Santa Cruz) and a far island (Marchena). However, almost all immigrants that stayed to breed were from an intermediate island (Santiago) and its satellites. Song may have been responsible for this selectivity. Mean heterozygosity stayed roughly constant over the next 15 years while allelic diversity almost doubled, after an initial decline, as the breeding population increased to a maximum of 30 pairs. Although close inbreeding occurred, with a reduction in heterozygosity, an expected net decline in heterozygosity did not occur, for two reasons: it was counteracted by continuing gene flow from immigrants at a low rate, and inbred birds (in one cohort) were at a selective disadvantage. An abrupt step-function shift in beak shape occurred after 9 years. Thus the study provides evidence of drift and selection causing morphological and genetic divergence in the establishment of a new population and in the first few generations.     

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment‐specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2‐5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.     

Herbarium specimens reveal a historical shift in phylogeographic structure of common ragweed during native range disturbance

Invasive plants provide ample opportunity to study evolutionary shifts that occur after introduction to novel environments. However, although genetic characters pre‐dating introduction can be important determinants of later success, large‐scale investigations of historical genetic structure have not been feasible. Common ragweed (Ambrosia artemisiifolia L.) is an invasive weed native to North America that is known for its allergenic pollen. Palynological records from sediment cores indicate that this species was uncommon before European colonization of North America, and ragweed populations expanded rapidly as settlers deforested the landscape on a massive scale, later becoming an aggressive invasive with populations established globally. Towards a direct comparison of genetic structure now and during intense anthropogenic disturbance of the late 19th century, we sampled 45 natural populations of common ragweed across its native range as well as historical herbarium specimens collected up to 140 years ago. Bayesian clustering analyses of 453 modern and 473 historical samples genotyped at three chloroplast spacer regions and six nuclear microsatellite loci reveal that historical ragweed's spatial genetic structure mirrors both the palaeo‐record of Ambrosia pollen deposition and the historical pattern of agricultural density across the landscape. Furthermore, for unknown reasons, this spatial genetic pattern has changed substantially in the intervening years. Following on previous work relating morphology and genetic expression between plants collected from eastern North America and Western Europe, we speculate that the cluster associated with humans’ rapid transformation of the landscape is a likely source of these aggressive invasive populations.     

Population genomic analyses reveal a history of range expansion and trait evolution across the native and invaded range of yellow starthistle (Centaurea solstitialis)

Identifying sources of genetic variation and reconstructing invasion routes for non‐native introduced species is central to understanding the circumstances under which they may evolve increased invasiveness. In this study, we used genome‐wide single nucleotide polymorphisms to study the colonization history of Centaurea solstitialis in its native range in Eurasia and invasions into the Americas. We leveraged this information to pinpoint key evolutionary shifts in plant size, a focal trait associated with invasiveness in this species. Our analyses revealed clear population genomic structure of potential source populations in Eurasia, including deep differentiation of a lineage found in the southern Apennine and Balkan Peninsulas and divergence among populations in Asia, eastern Europe and western Europe. We found strongest support for an evolutionary scenario in which western European populations were derived from an ancient admixture event between populations from eastern Europe and Asia, and subsequently served as the main genetic ‘bridgehead’ for introductions to the Americas. Introductions to California appear to be from a single source region, and multiple, independent introductions of divergent genotypes likely occurred into the Pacific Northwest. Plant size has evolved significantly at three points during range expansion, including a large size increase in the lineage responsible for the aggressive invasion of the California interior. These results reveal a long history of colonization, admixture and trait evolution in C. solstitialis, and suggest routes for improving evidence‐based management decisions for one of the most ecologically and economically damaging invasive species in the western United States.     

Comparative analyses of population structure in two subspecies of Nigella degenii: evidence for diversifying selection on pollen-color dimorphisms

Flower color can be a major determinant of plant fitness, not only because of preferential visitation by pollinators but also because of pleiotropic relationships between the expression of floral pigments and biochemically related compounds that influence vegetative performance variables. Different environments may therefore favor different pigmentation phenotypes. We examined whether spatially varying selection has played a major role in shaping large-scale patterns of differentiation in two subspecies of Nigella degenii (Ranunculaceae), with particular emphasis on pollen color. The two subspecies appear to have been genetically isolated for substantial periods of time and, therefore, provide a "replicated test" for the effect of natural selection. Estimates of population structure based on a suite of floral and vegetative characters were compared with the corresponding data for amplified fragment length polymorphism (AFLP) markers, which were assumed to be selectively neutral. We found low levels of genetic structure within the subspecies using both the AFLP markers (FST < or = 0.05) and quantitative characters (QST < or = 0.15), with no statistically significant differences between the two measures. There is, therefore, no evidence of diversifying selection being important in structuring variation in quantitative characters within each of the two subspecies. In contrast, estimates of differentiation in pollen color (FST > or = 0.10) significantly exceeded the neutral expectations (estimated from AFLP data), suggesting that local adaptation has played a key role in the evolution of this monogenic character.     

Testing the molecular and evolutionary causes of a 'leapfrog' pattern of geographical variation in coloration

Understanding the mechanisms accounting for the evolution of phenotypic diversity is central to evolutionary biology. We use molecular and phenotypic data to test hypotheses for 'leapfrog' patterns of geographical variation, in which phenotypically similar, disjunct populations are separated by distinct populations of the same species. Phylogenetic reconstructions revealed independent evolution of melanic plumage characters in different populations in the Neotropical avian genus Arremon. Thus, phenotypic similarities between distant populations cannot be explained by close phylogenetic affinity. Nor can they be attributed to recurring mutations in the MC1R gene, a locus involved in melanic pigmentation. A coalescent analysis indicates that plumage traits have become fixed at a faster rate than expected under genetic drift, suggesting that selection underlies their repeated evolution. In contrast to views that genetic drift drives phenotypic differentiation in Neotropical montane birds, our results imply that geographical variation preceding speciation may reflect the action of deterministic selective processes.     

Mixed population genomics support for the central marginal hypothesis across the invasive range of the cane toad (Rhinella marina) in Australia

Understanding factors that cause species' geographic range limits is a major focus in ecology and evolution. The central marginal hypothesis (CMH) predicts that species cannot adapt to conditions beyond current geographic range edges because genetic diversity decreases from core to edge due to smaller, more isolated edge populations. We employed a population genomics framework using 24 235–33 112 SNP loci to test major predictions of the CMH in the ongoing invasion of the cane toad (Rhinella marina) in Australia. Cane toad tissue samples were collected along broad‐scale, core‐to‐edge transects across their invasive range. Geographic and ecological core areas were identified using GIS and habitat suitability indices from ecological niche modelling. Bayesian clustering analyses revealed three genetic clusters, in the northwest invasion‐front region, northeast precipitation‐limited region and southeast cold temperature‐limited region. Core‐to‐edge patterns of genetic diversity and differentiation were consistent with the CMH in the southeast, but were not supported in the northeast and showed mixed support in the northwest. Results suggest cold temperatures are a likely contributor to southeastern range limits, consistent with CMH predictions. In the northeast and northwest, ecological processes consisting of a steep physiological barrier and ongoing invasion dynamics, respectively, are more likely explanations for population genomic patterns than the CMH.