Monitoring the effective population size of a brown bear (Ursus arctos) population using new single-sample approaches

The effective population size (N(e) ) could be the ideal parameter for monitoring populations of conservation concern as it conveniently summarizes both the evolutionary potential of the population and its sensitivity to genetic stochasticity. However, tracing its change through time is difficult in natural populations. We applied four new methods for estimating N(e) from a single sample of genotypes to trace temporal change in N(e) for bears in the Northern Dinaric Mountains. We genotyped 510 bears using 20 microsatellite loci and determined their age. The samples were organized into cohorts with regard to the year when the animals were born and yearly samples with age categories for every year when they were alive. We used the Estimator by Parentage Assignment (EPA) to directly estimate both N(e) and generation interval for each yearly sample. For cohorts, we estimated the effective number of breeders (N(b) ) using linkage disequilibrium, sibship assignment and approximate Bayesian computation methods and extrapolated these estimates to N(e) using the generation interval. The N(e) estimate by EPA is 276 (183-350 95% CI), meeting the inbreeding-avoidance criterion of N(e) > 50 but short of the long-term minimum viable population goal of N(e) > 500. The results obtained by the other methods are highly consistent with this result, and all indicate a rapid increase in N(e) probably in the late 1990s and early 2000s. The new single-sample approaches to the estimation of N(e) provide efficient means for including N(e) in monitoring frameworks and will be of great importance for future management and conservation.     

Population genomics as a new tool for wildlife management

Admixture and introgression have varied effects on population viability and fitness. Admixture might be an important source of new alleles, particularly for small, geographically isolated populations. However, admixture might also cause outbreeding depression if populations are adapted to different ecological or climatic conditions. Because of the emerging use of translocation and admixture as a conservation and wildlife management strategy to reduce genetic load (termed genetic rescue), the possible effects of admixture have practical consequences ( Bouzat et al. 2009 ; Hedrick & Fredrickson 2010 ). Importantly, genetic load and local adaptation are properties of individual loci and epistatic interactions among loci rather than properties of genomes. Likewise, the outcome and consequences of genetic rescue depend on the fitness effects of individual introduced alleles. In this issue of Molecular Ecology, Miller et al. (2012) use model‐based, population genomic analyses to document locus‐specific effects of a recent genetic rescue in the bighorn sheep population within the National Bison Range wildlife refuge (NBR; Montana, USA). They find a subset of introduced alleles associated with increased fitness in NBR bighorn sheep, some of which experienced accelerated introgression following their introduction. These loci mark regions of the genome that could constitute the genetic basis of the successful NBR bighorn sheep genetic rescue. Although population genomic analyses are frequently used to study local adaptation and selection (e.g. Hohenlohe et al. 2010 ; Lawniczak et al. 2010 ), this study constitutes a novel application of this analytical framework for wildlife management. Moreover, the detailed demographic data available for the NBR bighorn sheep population provide a rare and powerful source of information and allow more robust population genomic inference than is often possible.     

Pushing north one bottleneck at a time: site frequency spectra tell the history of Sitka spruce

Reconstructing the history of populations is a longstanding goal of molecular ecologists. In addition to a better understanding of the past, it is hoped that this knowledge would also facilitate predictions regarding species’ responses to future events such as climate change. The traditional way of doing this is through the fossil record, but these historical records are often incomplete. Inferring historical demography from patterns of nucleotide variability can help to fill these gaps. In this issue of Molecular Ecology, Holliday et al. (2010) glimpse into the demographic past of Sitka spruce, Picea sitchensis, an economically and ecologically important species native to northwestern United States and Canada, by examining the site frequency spectrum (SFS) of 153 loci in six populations covering the species entire range.     

Extinction,reintroduction, and restoration of a lizard meta‐population equilibrium in the missouri ozarks

In this issue of Molecular Ecology, Neuwald & Templeton (2013) report on a 22‐year study of natural populations of Collared Lizards (Crotaphytus collaris) that evolved on isolated on rock outcrops (‘glades’) in the Ozark Mountains in eastern Missouri. This ecosystem was originally maintained by frequent fires that kept the forest understory open, but fire‐suppression was adopted as official policy in about 1945, which led to a loss of native biodiversity, including local extinctions of some lizard populations. Policies aimed at restoring biodiversity included controlled burns and re‐introductions of lizards to some glades, which began in 1984. Populations were monitored from 1984–2006, and demographic and genetic data collected from 1 679 lizards were used to documents shifts in meta‐population dynamics over four distinct phases of lizard recovery: 1–an initial translocation of lizards drawn from the same source populations onto three glades that were likely part of one meta‐population; 2–a period of isolation and genetic drift associated with the absence of fires; 3–a period of rapid colonization and population increase following restoration of fire; and 4–stabilization of the meta‐population under regular prescribed burning. This study system thus provides a rare opportunity to characterize the dynamics of a landscape‐scale management strategy on the restoration of the meta‐population of a reintroduced species; long‐term case studies of the extinction, founding, increase, and stabilization of a well‐defined meta‐population, based on both demographic and population genetic data, are rare in the conservation, ecological, and evolutionary literature.     

A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour‐intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra 2008 ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.     

Multiple estimates of effective population size for monitoring a long‐lived vertebrate: an application to Yellowstone grizzly bears

Effective population size (N_e) is a key parameter for monitoring the genetic health of threatened populations because it reflects a population's evolutionary potential and risk of extinction due to genetic stochasticity. However, its application to wildlife monitoring has been limited because it is difficult to measure in natural populations. The isolated and well‐studied population of grizzly bears (Ursus arctos) in the Greater Yellowstone Ecosystem provides a rare opportunity to examine the usefulness of different N_e estimators for monitoring. We genotyped 729 Yellowstone grizzly bears using 20 microsatellites and applied three single‐sample estimators to examine contemporary trends in generation interval (GI), effective number of breeders (N_b) and N_e during 1982–2007. We also used multisample methods to estimate variance (N_eV) and inbreeding N_e (N_eI). Single‐sample estimates revealed positive trajectories, with over a fourfold increase in N_e (≈100 to 450) and near doubling of the GI (≈8 to 14) from the 1980s to 2000s. N_eV (240–319) and N_eI (256) were comparable with the harmonic mean single‐sample N_e (213) over the time period. Reanalysing historical data, we found N_eV increased from ≈80 in the 1910s–1960s to ≈280 in the contemporary population. The estimated ratio of effective to total census size (N_e/N_c) was stable and high (0.42–0.66) compared to previous brown bear studies. These results support independent demographic evidence for Yellowstone grizzly bear population growth since the 1980s. They further demonstrate how genetic monitoring of N_e can complement demographic‐based monitoring of N_c and vital rates, providing a valuable tool for wildlife managers.     

Let's talk about sex: A rigorous statistical framework to assign the sex of individuals from reduced‐representation sequencing data

Molecular markers have been used to identify the sex of sampled individuals for several decades, but the time‐consuming development phase prevented their application in many systems. Recently, a growing number of papers have applied reduced‐representation sequencing (RRS) protocols to the identification of sex‐specific markers without the use of test crosses or prior genomic information. While such an approach has great advantages in terms of versatility and ease of use, the “shotgun sequencing” nature of RRS data sets leads to a high amount of missing data, which results in statistical challenges to the confident assignment of sex to individuals. In this issue of Molecular Ecology Resources, Stovall et al. (Molecular Ecology Resources, 18, 2018) provide a statistical framework to answer two questions: (1) how many individuals of one sex only must possess a genotype for this locus to be considered significantly sex‐specific? and (2) How many sex‐specific loci must an individual of unknown sex possess (in a given data set) to be confidently assigned a sex? The statistical pipeline introduced, and applied to samples of New Zealand fur seal (Arctocephalus forsteri) to identify 90 sex‐specific loci, should be broadly applicable to a large number of species and constitutes a nice addition to the molecular ecology toolkit in the genomics era.     

Widespread yet heterogeneous genomic divergence

Genetic differentiation during adaptive divergence and speciation is heterogeneous among genomic regions. Some regions can be highly differentiated between populations, for example, because they harbour genes under divergent selection or those causing reproductive isolation and thus are resistant to gene flow. Other regions might be homogenized by gene flow and thus weakly differentiated. Debates persist about the number of differentiated regions expected under divergence with gene flow, and their causes, size, and genomic distribution. In this issue of Molecular Ecology, a study of freshwater stickleback used next-generation sequencing to shed novel insight into these issues (Roesti et al. 2012). Many genomic regions distributed across the genome were strongly differentiated, indicating divergence with gene flow can involve a greater number of loci than often thought. Nonetheless, differentiation of some regions, such as those near the centre of chromosomes where recombination is reduced, was strongly accentuated over others. Thus, divergence was widespread yet highly heterogeneous across the genome. Moreover, different population pairs varied in patterns of differentiation, illustrating how genomic divergence builds up across stages of the speciation process. The study demonstrates how variation in different evolutionary processes, such as selection and recombination rate, can combine to result in similar genomic patterns. Future work could focus on teasing apart the contributions of different processes for causing differentiation, a task facilitated by experimental manipulations.     

The chimerical genome of Isla del Coco feral pigs (Costa Rica), an isolated population since 1793 but with remarkable levels of diversity

The history of domestic species and of their wild ancestors is not a simple one, and feral processes can clarify key aspects of this history, including the adaptive processes triggered by new environments. Here, we provide a comprehensive genomic study of Isla del Coco (Costa Rica) feral pigs, a unique population that was allegedly founded by two individuals and has remained isolated since 1793. Using SNP arrays and genome sequencing, we show that Cocos pigs are hybrids between Asian and European pigs, as are modern international pig breeds. This conclusively shows that, as early as the 18th century, British vessels were loading crossbred pigs in Great Britain and transporting them overseas. We find that the Y chromosome has Asian origin, which has not been reported in any international pig breed. Chinese haplotypes seem to have been transmitted independently between Cocos and other pig breeds, suggesting independent introgression events and a complex pattern of admixing. Although data are compatible with a founder population of N = 2, variability levels are as high in Cocos pigs as in international pig breeds (~1.9 SNPs/kb) and higher than in European wild boars or local breeds (~1.7 SNPs/kb). Nevertheless, we also report a 10‐Mb region with a marked decrease in variability across all samples that contains four genes (CPE, H3F3C, SC4MOL and KHL2) previously identified as highly differentiated between wild and domestic pigs. This work therefore illustrates how feral population genomic studies can help to resolve the history of domestic species and associated admixture events.     

Genomic evidence for gene flow between monarchs with divergent migratory phenotypes and flight performance

Monarch butterflies are known for their spectacular annual migration in eastern North America, with millions of monarchs flying up to 4,500 km to overwintering sites in central Mexico. Monarchs also live west of the Rocky Mountains, where they travel shorter distances to overwinter along the Pacific Coast. It is often assumed that eastern and western monarchs form distinct evolutionary units, but genomic studies to support this notion are lacking. We used a tethered flight mill to show that migratory eastern monarchs have greater flight performance than western monarchs, consistent with their greater migratory distances. However, analysing more than 20 million SNPs in 43 monarch genomes, we found no evidence for genomic differentiation between eastern and western monarchs. Genomic analysis also showed identical and low levels of genetic diversity, and demographic analyses indicated similar effective population sizes and ongoing gene flow between eastern and western monarchs. Gene expression analysis of a subset of candidate genes during active flight revealed differential gene expression related to nonmuscular motor activity. Our results demonstrate that eastern and western monarchs maintain migratory differences despite ongoing gene flow, and suggest that migratory differences between eastern and western monarchs are not driven by select major‐effects alleles. Instead, variation in migratory distance and destination may be driven by environmentally induced differential gene expression or by many alleles of small effect.     

Patterns of genetic variation in US federal bison herds

Like many wide‐ranging mammals, American bison (Bison bison) have experienced significant range contraction over the past two centuries and are maintained in artificially isolated populations. A basic understanding of the distribution of genetic variation among populations is necessary to facilitate long‐term germplasm preservation and species conservation. The 11 herds maintained within the US federal system are a critically important source of germplasm for bison conservation, as they include many of the oldest herds in the USA and have served as a primary resource for the establishment of private and public herds worldwide. In this study, we used a panel of 51 nuclear markers to investigate patterns of neutral genetic variation among these herds. Most of these herds have maintained remarkably high levels of variation despite the severe bottleneck suffered in the late 1800s. However, differences were noted in the patterns of variation and levels of differentiation among herds, which were compared with historical records of establishment, supplementation, herd size, and culling practices. Although some lineages have been replicated across multiple herds within the US federal system, other lineages with high levels of genetic variation exist in isolated herds and should be considered targets for the establishment of satellite herds. From this and other studies, it is clear that the genetic variation represented in the US federal system is unevenly distributed among National Park Service and Fish and Wildlife Service herds, and that these resources must be carefully managed to ensure long‐term species conservation.     

Utility of pooled sequencing for association mapping in nonmodel organisms

High‐density genome‐wide sequencing increases the likelihood of discovering genes of major effect and genomic structural variation in organisms. While there is an increasing availability of reference genomes across broad taxa, the greatest limitation to whole‐genome sequencing of multiple individuals continues to be the costs associated with sequencing. To alleviate excessive costs, pooling multiple individuals with similar phenotypes and sequencing the homogenized DNA (Pool‐Seq) can achieve high genome coverage, but at the loss of individual genotypes. Although Pool‐Seq has been an effective method for association mapping in model organisms, it has not been frequently utilized in natural populations. To extend bioinformatic tools for rapid implementation of Pool‐Seq data in nonmodel organisms, we developed a pipeline called PoolParty and illustrate its effectiveness in genetic association mapping. Alignment expectations based on five pooled Chinook salmon (Oncorhynchus tshawytscha) libraries showed that approximately 48% genome coverage per library could be achieved with reasonable sequencing effort. We additionally examined male and female O. tshawytscha libraries to illustrate how Pool‐Seq techniques can successfully map known genes associated with functional differences among sexes such as growth hormone 2. Finally, we compared pools of individuals of different spawning ages for each sex to discover novel genes involved with age at maturity in O. tshawytscha such as opsin4 and transmembrane protein19. While not appropriate for every system, Pool‐Seq data processed by the PoolParty pipeline is a practical method for identifying genes of major effect in nonmodel organisms when high genome coverage is necessary and cost is a limiting factor.     

Divergence with gene flow in a population of thermophilic bacteria: a potential role for spatially varying selection

A fundamental goal of evolutionary biology is to understand how ecological diversity arises and is maintained in natural populations. We have investigated the contributions of gene flow and divergent selection to the distribution of genetic variation in an ecologically differentiated population of a thermophilic cyanobacterium (Mastigocladus laminosus) found along the temperature gradient of a nitrogen‐limited stream in Yellowstone National Park. For most loci sampled, gene flow appears to be sufficient to prevent substantial genetic divergence. However, one locus (rfbC) exhibited a comparatively low migration rate as well as other signatures expected for a gene experiencing spatially varying selection, including an excess of common variants, an elevated level of polymorphism and extreme genetic differentiation along the gradient. rfbC is part of an expression island involved in the production of the polysaccharide component of the protective envelope of the heterocyst, the specialized nitrogen‐fixing cell of these bacteria. SNP genotyping in the vicinity of rfbC revealed a ~5‐kbp region including a gene content polymorphism that is tightly associated with environmental temperature and therefore likely contains the target of selection. Two genes have been deleted both in the predominant haplotype found in the downstream region of White Creek and in strains from other Yellowstone populations of M. laminosus, which may result in the production of heterocysts with different envelope properties. This study implicates spatially varying selection in the maintenance of variation related to thermal performance at White Creek despite on‐going or recent gene flow.     

Critical review of NGS analyses for de novo genotyping multigene families

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next‐generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra‐deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the ‘degree of change’ (DOC) method avoids using arbitrary cut‐off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third‐generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second‐generation sequencing approaches.     

Use of genotyping by sequencing data to develop a high‐throughput and multifunctional SNP panel for conservation applications in Pacific lamprey

Next‐generation sequencing data can be mined for highly informative single nucleotide polymorphisms (SNPs) to develop high‐throughput genomic assays for nonmodel organisms. However, choosing a set of SNPs to address a variety of objectives can be difficult because SNPs are often not equally informative. We developed an optimal combination of 96 high‐throughput SNP assays from a total of 4439 SNPs identified in a previous study of Pacific lamprey (Entosphenus tridentatus) and used them to address four disparate objectives: parentage analysis, species identification and characterization of neutral and adaptive variation. Nine of these SNPs are F_ST outliers, and five of these outliers are localized within genes and significantly associated with geography, run‐timing and dwarf life history. Two of the 96 SNPs were diagnostic for two other lamprey species that were morphologically indistinguishable at early larval stages and were sympatric in the Pacific Northwest. The majority (85) of SNPs in the panel were highly informative for parentage analysis, that is, putatively neutral with high minor allele frequency across the species’ range. Results from three case studies are presented to demonstrate the broad utility of this panel of SNP markers in this species. As Pacific lamprey populations are undergoing rapid decline, these SNPs provide an important resource to address critical uncertainties associated with the conservation and recovery of this imperiled species.     

Illuminating the mysteries of wolf history

One of the most enduring surprises about the genetic history of Late Pleistocene populations is that continuity is often disturbed by upheaval. In fact, studies that support population continuity are increasingly rare in humans, a variety of vertebrate taxa, and vascular plants (Hofreiter & Stewart 2009; Burbrink et al. 2016). Perhaps such continuity should not be expected as the Pleistocene is marked by episodes of climate change, glaciation and the invasions of humans into previously isolated areas. Although fossils are one of the primary sources for inferring population continuity, a problem with fossil material is that, even if similar morphological forms might exist in a place over time, they may not be from the same genetic lineage. There are now readily available methods to assess genetic continuity solely from DNA found in fossil material, provided the record is fairly continuous. In a From the Cover article in this issue of Molecular Ecology, Loog et al. (2020) apply some of these readily available methods to analyse mitochondrial genomes and model the demography of wolves over the last 50,000 years.     

Is Niagara Falls a barrier to gene flow in riverine fishes? A test using genome‐wide SNP data from seven native species

Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post‐glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water‐flow drive the downstream accumulation of genetic diversity, whereas post‐glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m high Niagara Falls – a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls’ role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome‐wide SNPs and coalescent models to test whether populations are: (a) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (b) more genetically diverse upstream, consistent with post‐glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and (c) have migrated either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to migration, but two species have probably dispersed upstream past the falls after glacial retreat yet before opening of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species’ genetic diversity and highlight the need to better understand post‐glacial dispersal pathways.     

Structure_threader: An improved method for automation and parallelization of programs structure,fastStructure and MavericK on multicore CPU systems

Structure_threader is a program to parallelize multiple runs of genetic clustering software that does not make use of multithreading technology (structure , fastStructure and MavericK) on multicore computers. Our approach was benchmarked across multiple systems and displayed great speed improvements relative to the single‐threaded implementation, scaling very close to linearly with the number of physical cores used. Structure_threader was compared to previous software written for the same task—ParallelStructure and StrAuto and was proven to be the faster (up to 25% faster) wrapper under all tested scenarios. Furthermore, Structure_threader can perform several automatic and convenient operations, assisting the user in assessing the most biologically likely value of ‘K’ via implementations such as the “Evanno,” or “Thermodynamic Integration” tests and automatically draw the “meanQ” plots (static or interactive) for each value of K (or even combined plots). Structure_threader is written in python 3 and licensed under the GPLv3. It can be downloaded free of charge at https://github.com/StuntsPT/Structure_threader .     

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.     

The genomic and physiological basis of life history variation in a butterfly metapopulation

Unravelling the mechanisms underlying variation in life history traits is of fundamental importance for our understanding of adaptation by natural selection. While progress has been made in mapping fitness-related phenotypes to genotypes, mainly in a handful of model organisms, functional genomic studies of life history adaptations are still in their infancy. In particular, despite a few notable exceptions, the genomic basis of life history variation in natural populations remains poorly understood. This is especially true for the genetic underpinnings of life history phenotypes subject to diversifying selection driven by ecological dynamics in patchy environments--as opposed to adaptations involving strong directional selection owing to major environmental changes, such as latitudinal gradients, extreme climatic events or transitions from salt to freshwater. In this issue of Molecular Ecology,Wheat et al. (2011) now make a significant leap forward by applying the tools of functional genomics to dispersal-related life history variation in a butterfly metapopulation. Using a combination of microarrays, quantitative PCR and physiological measurements, the authors uncover several metabolic and endocrine factors that likely contribute to the observed life history phenotypes. By identifying molecular candidate mechanisms of fitness variation maintained by dispersal dynamics in a heterogeneous environment,they also begin to address fascinating interactions between the levels of physiology, ecology and evolution.