Genomics of local adaptation with gene flow

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher Fst values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher F_st values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

The genetics of adaptation to discrete heterogeneous environments: frequent mutation or large-effect alleles can allow range expansion

Range expansions are complex evolutionary and ecological processes. From an evolutionary standpoint, a populations' adaptive capacity can determine the success or failure of expansion. Using individual-based simulations, we model range expansion over a two-dimensional, approximately continuous landscape. We investigate the ability of populations to adapt across patchy environmental gradients and examine how the effect sizes of mutations influence the ability to adapt to novel environments during range expansion. We find that genetic architecture and landscape patchiness both have the ability to change the outcome of adaptation and expansion over the landscape. Adaptation to new environments succeeds via many mutations of small effect or few of large effect, but not via the intermediate between these cases. Higher genetic variance contributes to increased ability to adapt, but an alternative route of successful adaptation can proceed from low genetic variance scenarios with alleles of sufficiently large effect. Steeper environmental gradients can prevent adaptation and range expansion on both linear and patchy landscapes. When the landscape is partitioned into local patches with sharp changes in phenotypic optimum, the local magnitude of change between subsequent patches in the environment determines the success of adaptation to new patches during expansion.

Size selection by a gape-limited predator of a marine snail: Insights into magic traits for speciation

The intertidal snail Littorina saxatilis has repeatedly evolved two parallel ecotypes assumed to be wave adapted and predatory shore crab adapted, but the magnitude and targets of predator‐driven selection are unknown. In Spain, a small, wave ecotype with a large aperture from the lower shore and a large, thick‐shelled crab ecotype from the upper shore meet in the mid‐shore and show partial size‐assortative mating. We performed complementary field tethering and laboratory predation experiments; the first set compared the survival of two different size‐classes of the crab ecotype while the second compared the same size‐class of the two ecotypes. In the first set, the large size‐class of the crab ecotype survived significantly better than the small size‐class both on the upper shore and in the laboratory. In the second set, the small size‐class of the crab ecotype survived substantially better than that of the wave ecotype both on the upper shore and in the laboratory. Shell‐breaking predation on tethered snails was almost absent within the lower shore. In the laboratory shore crabs (Pachygrapsus marmoratus) with larger claw heights selected most strongly against the small size‐class of the crab ecotype, whereas those with medium claw heights selected most strongly against the thin‐shelled wave ecotype. Sexual maturity occurred at a much larger size in the crab ecotype than in the wave ecotype. Our results showed that selection on the upper shore for rapid attainment of a size refuge from this gape‐limited predator favors large size, thick shells, and late maturity. Model parameterization showed that size‐selective predation restricted to the upper shore resulted in the evolution of the crab ecotype despite gene flow from the wave ecotype snails living on the lower shore. These results on gape‐limited predation and previous ones showing size‐assortative mating between ecotypes suggest that size may represent a magic trait for the thick‐shelled ecotype.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

Strong neutral genetic differentiation in a host, but not in its parasite

The genetic diversity and population structure of a parasite with a complex life cycle generally depends on the dispersal by its most motile host. Given that high gene flow is assumed to hinder local adaptation, this can impose significant constraints on a parasite's potential to adapt to local environmental conditions, intermediate host populations, and ultimately to host-parasite coevolution. Here, we aimed to examine the population genetic basis for local host-parasite interactions between the eye fluke Diplostomum lineage 6, a digenean trematode with a multi-host life cycle (including a snail, a fish, and a bird) and its second intermediate host, the three-spined stickleback Gasterosteus aculeatus L. We developed the first microsatellite primers for D. lineage 6 and used them together with published stickleback markers to analyse host and parasite population structures in 19 freshwater lakes, which differ in their local environmental characteristics regarding water chemistry and Diplostomum abundance. Our analyses suggest that one parasite population successfully infects a range of genetically differentiated stickleback populations. The lack of neutral genetic differentiation in D. lineage 6, which could be attributed to the motility of the parasite's definitive host as well as its life cycle characteristics, makes local host-parasite co-adaptations seem more likely on a larger geographical scale than among the lakes of our study site. Our study provides a suitable background for future studies in this system and the first microsatellite primers for a widespread fish parasite.