Islands increase genetic subdivision and disrupt patterns of connectivity of intertidal snails in a complex archipelago

The view that marine species with planktonic dispersal have highly connected, demographically open populations is giving way to recognition that populations may often be largely self-recruiting, or demographically closed. This raises the question of what local conditions might favor isolation of populations. To test the importance of islands for local isolation in species with planktonic larvae, we examined allozyme variation among 35 populations of the intertidal snail Austrocochlea constricta in the Houtman Abrolhos Islands, Western Australia, spanning 60 km. Heterogeneity of allozyme frequencies among populations was high, with average F(ST) of 0.237, indicating highly localized populations. Increased subdivision was associated with islands at different scales: between island groups, separated by deep water gaps, and between disconnected sets of islands within groups. At short distances, up to two km, subdivision increased fivefold between islands compared with that between populations on the same island. Along 11 km of continuous, sheltered shore, there was isolation by distance but among a linear series of islands over similar distance, there was greater subdivision at short distances but no association with distance. These patterns had been seen previously in the direct-developing snail Bembicium vittatum, but its finding in A. constricta confirms for a planktonic disperser the importance of this complex archipelago for both retention of locally produced larvae and disruption of patterns of connectivity. Taken together, these results indicate that islands can increase both the "open" and the "closed" components of recruitment and that applicable models of genetic connectivity depend substantially on local conditions.     

Genetic and environmental components of phenotypic variation in body shape among populations of Atlantic cod (Gadus morhua L.)

A common-garden experiment was conducted on larvae to test for genetic differences in body shape among populations of Atlantic cod ( Gadus morhua ). Offspring from four north-west Atlantic regions were reared from hatching to postmetamorphosis at two temperatures (7 ± 1 °C and 11 ± 1 °C) and two food levels (1500 and 4500 prey L⁻¹). Body shape differed between populations and treatments. Population differences were greatest between south-west Scotian Shelf cod and those further north; the former were characterized by a deeper body, larger head, and longer caudal peduncle than cod from the other populations. Significant differences were also observed between two putative populations on the south-west Scotian Shelf, suggesting genetic divergence between spawning aggregations at small spatial scales (< 100 km). Temperature and food supply also influenced body shape, with the effect of the former being more pronounced. Individuals reared at the higher temperature or food level had a deeper body and a larger head than those reared at the lower temperature or food supply. Phenotypic responses to changes in the rearing environment also differed among populations, indicating genetic differences in phenotypic plasticity. Differences between populations in morphology and in phenotypic plasticity suggest genetic divergence at both large (> 1000 km) and small (< 100 km) spatial scales. The genetic differences at large spatial scales counteracted the expected effects of temperature differences in the wild, suggesting countergradient variation in morphology among these populations.  © 2006 Her Majesty the Queen in Right of Canada. Journal compilation © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 351–365.     

Recruitment limitation: A theoretical perspective

Abstract A theoretical analysis of the concept of recruitment limitation leads to the conclusion that most populations should he regarded as jointly limited by recruitment and interactions between individuals after recruitment. The open nature of local marine systems does not permit avoidance of density-dependent interactions; it simply may make them more difficult to detect. Local populations consisting of settled organisms may not experience density-dependent interactions under some circumstances, but the entire species population consisting of the collection of local populations and their planktonic larvae must have density-dependent dynamics. Any local population of settled individuals can escape density dependence if sufficient density dependence occurs among planktonic larvae or within other local populations. Common conceptions of density dependence are too narrow, leading too often to the conclusion that it is absent from a system. It is equally wrong to expect that density-dependent interactions after settlement determine local population densities independently of recruitment. Special circumstances allowing density dependence to act strongly and quickly are needed before density dependence can neutralize the effects of recruitment. Recruitment limitation and density-dependent interactions therefore should not be regarded as alternatives but as jointly acting to determine the densities of marine benthic populations. Moreover, the interaction between fluctuating recruitment and density dependence is potentially the most interesting feature of recruitment limitation. For example, this interaction may be an important diversity-maintaining mechanism for marine systems.     

Polymorphism in developmental mode and its effect on population genetic structure of a spionid polychaete, Pygospio elegans

Population genetic structure of sedentary marine species is expected to be shaped mainly by the dispersal ability of their larvae. Long-lived planktonic larvae can connect populations through migration and gene flow, whereas species with nondispersive benthic or direct-developing larvae are expected to have genetically differentiated populations. Poecilogonous species producing different larval types are ideal when studying the effect of developmental mode on population genetic structure and connectivity. In the spionid polychaete Pygospio elegans, different larval types have been observed between, and sometimes also within, populations. We used microsatellite markers to study population structure of European P. elegans from the Baltic Sea (BS) and North Sea (NS). We found that populations with planktonic larvae had higher genetic diversity than did populations with benthic larvae. However, this pattern may not be related to developmental mode, since in P. elegans, developmental mode may be associated with geography. Benthic larvae were more commonly seen in the brackish BS and planktonic larvae were predominant in the NS, although both larval types also are found from both areas. Significant isolation-by-distance (IBD) was found overall and within regions. Most of the pair-wise F(ST) comparisons among populations were significant, although some geographically close populations with planktonic larvae were found to be genetically similar. However, these results, together with the pattern of IBD, autocorrelation within populations, as well as high estimated local recruitment, suggest that dispersal is limited in populations with planktonic larvae as well as in those with benthic larvae. The decrease in salinity between the NS and BS causes a barrier to gene flow in many marine species. In P. elegans, low, but significant, differentiation was detected between the NS and BS (3.34% in AMOVA), but no clear transition zone was observed, indicating that larvae are not hampered by the change in salinity.     

Within-Population Developmental and Morphological Plasticity is Mirrored in Between-Population Differences: Linking Plasticity and Diversity

It has been suggested that phenotypic plasticity can facilitate evolutionary diversification of organisms. If life-history and morphological diversification across a lineage is mirrored in diversification in the same traits due to phenotypic plasticity within a lineage it fulfils one of the expectations that are needed to support this diversification hypothesis. We carried out a laboratory study to examine development rate and morphology between and within populations of the parsley frog, Pelodytes punctatus. We found that frogs reared in the laboratory had a longer development time, relatively longer hind legs and relatively narrower heads under constant water level compared to those under decreasing water level simulating pool drying. This adaptive phenotypic plasticity response to pool drying was mirrored across populations because frogs from permanent waters had longer development times, relatively longer hind legs and relatively narrower heads compared to frogs from temporary waters. Hence the developmental and morphological plasticity observed within populations was also observed between populations as constitutive expressed traits. We suggest that the morphology pattern observed was driven by a common developmental process (time to metamorphosis), indicating that plasticity may contribute to evolutionary change, ultimately resulting in genetic accommodation of the morphological traits.     

Plasticity matches phenotype to local conditions despite genetic homogeneity across 13 snake populations

In a widespread species, a matching of phenotypic traits to local environmental optima is generally attributed to site-specific adaptation. However, the same matching can occur via adaptive plasticity, without requiring genetic differences among populations. Adult sea kraits (Laticauda saintgironsi) are highly philopatric to small islands, but the entire population within the Neo-Caledonian Lagoon is genetically homogeneous because females migrate to the mainland to lay their eggs at communal sites; recruits disperse before settling, mixing up alleles. Consequently, any matching between local environments (e.g. prey sizes) and snake phenotypes (e.g. body sizes and relative jaw sizes (RJSs)) must be achieved via phenotypic plasticity rather than spatial heterogeneity in gene frequencies. We sampled 13 snake colonies spread along an approximately 200 km northwest–southeast gradient (n > 4500 individuals) to measure two morphological features that affect maximum ingestible prey size in gape-limited predators: body size and RJS. As proxies of habitat quality (HQ), we used protection status, fishing pressure and lagoon characteristics (lagoon width and distance of islands to the barrier reef). In both sexes, spatial variation in body sizes and RJSs was linked to HQ; albeit in different ways, consistent with sex-based divergences in foraging ecology. Strong spatial divergence in morphology among snake colonies, despite genetic homogeneity, supports the idea that phenotypic plasticity can facilitate speciation by creating multiple phenotypically distinct subpopulations shaped by their environment.     

Individual diet specialisation in sparrows is driven by phenotypic plasticity in traits related to trade‐offs in animal performance

Individual diet specialisation (IS) is frequent in many animal taxa and affects population and community dynamics. The niche variation hypothesis (NVH) predicts that broader population niches should exhibit greater IS than populations with narrower niches, and most studies that examine the ecological factors driving IS focus on intraspecific competition. We show that phenotypic plasticity of traits associated with functional trade‐offs is an important, but unrecognised mechanism that promotes and maintains IS. We measured nitrogen isotope (δ¹⁵N) and digestive enzyme plasticity in four populations of sparrows (Zonotrichia capensis) to explore the relationship between IS and digestive plasticity. Our results show that phenotypic plasticity associated with functional trade‐offs is related in a nonlinear fashion with the degree of IS and positively with population niche width. These findings are opposite to the NVH and suggest that among individual differences in diet can be maintained via acclimatisation and not necessarily require a genetic component.     

Morphological anti‐predator defences in the nine‐spined stickleback: constitutive,induced or both?

Environmental differences among populations are expected to lead to local adaptation, while spatial or temporal environmental variation within a population will favour evolution of phenotypic plasticity. As plasticity itself can be under selection, locally adapted populations can vary in levels of plasticity. Nine‐spined stickleback (Pungitius pungitius) originating from isolated ponds (low piscine predation risk, high competition) vs. lake and marine populations (high piscine predation risk, low competition) are known to be morphologically adapted to their respective environments. However, nothing is known about their ability to express phenotypic plasticity in morphology in response to perceived predation risk or food availability/competition. We studied predator‐induced phenotypic plasticity in body shape and armour of marine and pond nine‐spined stickleback in a factorial common garden experiment with two predator treatments (present vs. absent) and two feeding regimes (low vs. high). The predation treatment did not induce any morphological shifts in fish from either habitat or food regime. However, strong habitat‐dependent differences between populations as well as strong sexual dimorphism in both body shape and armour were found. The lack of predator‐induced plasticity in development of the defence traits (viz. body armour and body depth) suggests that morphological anti‐predator traits in nine‐spined stickleback are strictly constitutive, rather than inducible. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

The origins of ecotypic variation of rainbow trout: a test of environmental vs. genetically based differences in morphology

Although morphological plasticity has been observed in a variety of taxa, few experimental studies have compared the relative proportion of morphological variability that is accounted for by environmentally induced plasticity, and how much is because of genetically based differences among populations. We compared the morphology of six rainbow trout (Oncorhynchus mykiss) populations from different ecotypic categories that were raised under flowing vs. standing-water conditions. Our data indicate that both environmental conditions and ecotypic differences account for a significant proportion of variation in morphology. Among ecotype effects, however, accounted for a much larger proportion of morphological variability than environmental conditions. Rainbow trout from stream populations had deeper caudal peduncles, and longer fins than lake populations, and rainbow trout from a piscivorous population had larger mouth and head lengths than all other ecotypes. Environmentally induced differences in morphology were primarily related to differences in mouth and head lengths, as well as fin length. Relative to morphometric differences from natural rainbow trout populations, most characteristics deviated in the same direction in our experimental populations. Our data indicate that morphological differences across rainbow trout populations have a genetic basis and may represent locally adaptive characteristics and highlight the role of ecology in promoting phenotypic divergence.     

Variable salinity tolerance in ascidian larvae is primarily a plastic response to the parental environment

Both phenotypic plasticity and local genetic adaptation may contribute to a species’ ability to inhabit different environmental conditions. While phenotypic plasticity is usually considered costly, local adaptation takes generations to respond to environmental change and may be constrained by strong gene flow. The majority of marine species have complex life-cycles with pelagic stages that might be expected to promote gene flow and plastic responses, and yet several notable examples of local adaptation have been found in species with broadcast larvae. In the ascidian, Ciona intestinalis (Linnaeus, 1767),—a common marine species with broadcast spawning and a short larval stage—previous studies have found marked differences in salinity tolerance of early life-history stages among populations from different salinity regimes. We used common-garden experiments to test whether observed differences in salinity tolerance could be explained by phenotypic plasticity. Adult ascidians from two low salinity populations [2–5 m depth, ~25 practical salinity units (PSU)], and two full salinity populations (25–27 m depth, ~31 PSU) were acclimated for 2–4 weeks at both 25 and 31 PSU. Gametes were fertilized at the acclimation salinities, and the newly formed embryos were transferred to 10 different salinities (21–39 PSU) and cultured to metamorphosis. Adult acclimation salinity had an overriding and significant effect on larval metamorphic success: tolerance norms for larvae almost fully matched the acclimation salinity of the parents, independent of parental origin (deep or shallow). However we also detected minor population differences that could be attributed to either local adaptation or persistent environmental effects. We conclude that differences in salinity tolerance of C. intestinalis larvae from different populations are driven primarily by transgenerational phenotypic plasticity, a strategy that seems particularly favourable for an organism living in coastal waters where salinity is less readily predicted than in the open oceans.     

Temporal dynamics of NMDA receptor-induced changes in spine morphology and AMPA receptor recruitment to spines

Recent studies have shown that the activation of NMDA receptors can induce rapid changes in dendritic morphology and synaptic recruitment of AMPA receptors in dendritic spines. Here, we analyze the time course of NMDA receptor-induced changes in dendrite morphology and recruitment of AMPA receptors to synapses in cultured neurons. Activation of NMDA receptors causes a rapid transient increase in the size of preexisting spines and then the gradual formation of new dendritic protrusions and spines. NMDA receptor activation also induced GFP-tagged AMPA receptors to cluster in dendrites and to be inserted into the surface of dendritic spines. These results indicate that NMDA receptor activation induces several phases of dendritic plasticity, initial expansion of dendritic spines, followed by the de novo formation of spines and AMPA receptor dendritic clustering and surface expression on spines. Each of these forms of plasticity may have significant effects on the efficacy of synaptic transmission.     

Natural selection, larval dispersal, and the geography of phenotype in the sea

Populations evolve generalist, specialist, and plastic strategies in response to environmental heterogeneity. Describing such within-species variation in phenotype and how it arises is central to understanding a variety of ecological and evolutionary topics. The literature on phenotypic differences among populations is highly biased; for every one article published on a marine species, at least 10 articles are published on a terrestrial species and eight focus on terrestrial plants. Here, I outline what we know from the marine literature about geographic variation in phenotype in the sea, with a principal focus on local adaptation. The theory of environmental "grain" predicts that the most likely evolutionary response (e.g., local adaptation, phenotypic plasticity, generalism, and balanced polymorphism) depends on the spatial scale of environmental variation relative to the distance that an organism disperses. Consistent with these predictions, phenotypic plasticity is stronger among invertebrates with geographically broad dispersal versus restricted dispersal (i.e., planktonic-dispersers versus direct-developers). However, contrary to predictions, the relative frequency, and spatial scale of local adaptation is not consistently greater among direct-developers relative to planktonic disperers. This indicates that the likelihood of local adaptation depends on other organismal or environmental traits. Two of the most vexing issues that remain include (1) predicting the extent to which barriers to dispersal are a cause versus consequence of phenotypic differentiation and (2) delineating the relative importance of evolutionary forces that favor or impede local adaptation. Understanding the mechanistic basis of the geography of phenotypic differences, or phenogeography, has gained recent momentum because of a need to predict impacts of global climatic change, anthropogenic disturbances, and dispersal of organisms to non-native habitats.     

Microgeographic adaptations of spotted salamander morphological defenses in response to a predaceous salamander and beetle

Spatial heterogeneity in the selection imposed by different predator species could promote the adaptive diversification of local prey populations. However, high gene flow might swamp local adaptations at limited spatial scales or generalized phenotypic plasticity might evolve in place of local diversification. Spotted salamander larvae Ambystoma maculatum face strongly varying risks from gape‐limited marbled salamander larvae Ambystoma opacum and gape‐unconstrained diving beetle larvae Dytiscus spp. across natural landscapes. To evaluate if A. maculatum adapts to these predation risk across micro‐geographic scales, I measured selection gradients in response to the two focal predators and then assayed the defensive morphologies of ten populations in a common garden experiment. I found that A. opacum induced selection on A. maculatum for larger tailfins and bodies whereas beetles induced selection for larger tail muscles and smaller bodies. In accordance with the local adaptation hypothesis, A. maculatum populations inhabiting ponds with high beetle densities grew larger tail muscles relative to other populations when raised in a common environment. However, populations exposed to strong A. opacum selection did not evolve larger tailfins as predicted. High gene flow or morphological plasticity could explain the absence of this morphological response to A. opacum. Overall, results suggest that populations can sometimes evolve adaptive traits in response to locally variable selection regimes even across the very limited distances that separate populations in this study. If prey populations often differ in their defenses against local predators, then this variation could affect the outcome of species interactions in local communities.     

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.     

Phenotypic divergence between seasnake (Emydocephalus annulatus) populations from adjacent bays of the New Caledonian Lagoon

Populations of widespread species often differ in phenotypic traits, although rarely in such a dramatic fashion as revealed by research on turtle‐headed seasnakes (Emydocephalus annulatus). These snakes are highly philopatric, with mark–recapture studies showing that the interchange of individuals rarely occurs even between two adjacent bays (separated by < 1.2 km) in Noumea, New Caledonia. Data on > 500 field‐captured snakes from these two bays reveal significant differences between these two locations in snake morphology (mean body length, relative tail length, head shape), colour, ecology (body condition, growth rate, incidence of algal fouling), behaviour (antipredator tactics), and locomotor performance. For some traits, the disparity was very marked (e.g. mean swimming speeds differed by > 30%). The causal bases for these phenotypic divergences may involve founder effects, local adaptation, and phenotypic plasticity. The spatial divergence in phenotypic traits offers a cautionary tale both for researchers (sampling of only a few populations may fail to provide a valid overview of the morphology, performance, and behaviour of a species) and managers (loss of local populations may eliminate distinctive genetic variation). © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Isolation and characterization of 11 polymorphic microsatellite markers for the marine gastropod Concholepas concholepas (Brugière, 1789)

Because of its long‐lived planktonic stage, the marine gastropod Concholepas concholepas is expected to exchange larvae over large distances. However, discrepancies between expected and realized dispersal have been documented in marine invertebrates. To investigate relationships between potential and effective (i.e. gene flow) dispersal, we developed 11 microsatellite markers and investigate their usefulness by analysing two populations distant by c. 4000 km. The 11 loci were found to be highly polymorphic in both populations, with 12–51 alleles according to the locus. This polymorphism is strong enough to allow fine‐scale population analyses including larval studies and paternity analyses.     

Range margin populations show high climate adaptation lags in European trees

How populations of long‐living species respond to climate change depends on phenotypic plasticity and local adaptation processes. Marginal populations are expected to have lags in adaptation (i.e. differences between the climatic optimum that maximizes population fitness and the local climate) because they receive pre‐adapted alleles from core populations preventing them from reaching a local optimum in their climatically marginal habitat. Yet, whether adaptation lags in marginal populations are a common feature across phylogenetically and ecologically different species and how lags can change with climate change remain unexplored. To test for range‐wide patterns of phenotypic variation and adaptation lags of populations to climate, we (a) built model ensembles of tree height accounting for the climate of population origin and the climate of the site for 706 populations monitored in 97 common garden experiments covering the range of six European forest tree species; (b) estimated populations' adaptation lags as the differences between the climatic optimum that maximizes tree height and the climate of the origin of each population; (c) identified adaptation lag patterns for populations coming from the warm/dry and cold/wet margins and from the distribution core of each species range. We found that (a) phenotypic variation is driven by either temperature or precipitation; (b) adaptation lags are consistently higher in climatic margin populations (cold/warm, dry/wet) than in core populations; (c) predictions for future warmer climates suggest adaptation lags would decrease in cold margin populations, slightly increasing tree height, while adaptation lags would increase in core and warm margin populations, sharply decreasing tree height. Our results suggest that warm margin populations are the most vulnerable to climate change, but understanding how these populations can cope with future climates depend on whether other fitness‐related traits could show similar adaptation lag patterns.     

Multiple scales of genetic connectivity in a brooding coral on isolated reefs following catastrophic bleaching

Understanding the pattern of connectivity among populations is crucial for the development of realistic and spatially explicit population models in marine systems. Here we analysed variation at eight microsatellite loci to assess the genetic structure and to infer patterns of larval dispersal for a brooding coral, Seriatopora hystrix, at an isolated system of reefs in northern Western Australia. Spatial autocorrelation analyses show that populations are locally subdivided, and that the majority of larvae recruit to within 100 m of their natal colony. Further, a combination of F- and R- statistics showed significant differentiation at larger spatial scales (2-60 km) between sites, and this pattern was clearly not associated with distance. However, Bayesian analysis demonstrated that recruitment has been supplemented by less frequent but recent input of larvae from outside the local area; 2-6% of colonies were excluded from the site at which they were sampled. Individual assignments of these migrants to the most likely populations suggest that the majority of migrants were produced at the only site that was not decimated by a recent and catastrophic coral bleaching event. Furthermore, the only site that recovered to prebleaching levels received most of these immigrants. We conclude that the genetic structure of this brooding coral reflects its highly opportunistic life history, in which prolific, philopatric recruitment is occasionally supplemented by exogenously produced larvae.     

LOCAL ADAPTATION AND THE EVOLUTION OF PHENOTYPIC PLASTICITY IN TRINIDADIAN GUPPIES (POECILIA RETICULATA)

Divergent selection pressures across environments can result in phenotypic differentiation that is due to local adaptation, phenotypic plasticity, or both. Trinidadian guppies exhibit local adaptation to the presence or absence of predators, but the degree to which predator‐induced plasticity contributes to population differentiation is less clear. We conducted common garden experiments on guppies obtained from two drainages containing populations adapted to high‐ and low‐predation environments. We reared full‐siblings from all populations in treatments simulating the presumed ancestral (predator cues present) and derived (predator cues absent) conditions and measured water column use, head morphology, and size at maturity. When reared in presence of predator cues, all populations had phenotypes that were typical of a high‐predation ecotype. However, when reared in the absence of predator cues, guppies from high‐ and low‐predation regimes differed in head morphology and size at maturity; the qualitative nature of these differences corresponded to those that characterize adaptive phenotypes in high‐ versus low‐predation environments. Thus, divergence in plasticity is due to phenotypic differences between high‐ and low‐predation populations when reared in the absence of predator cues. These results suggest that plasticity might initially play an important role during colonization of novel environments, and then evolve as a by‐product of adaptation to the derived environment.     

Spatial phenotypic and genetic structure of threespine stickleback (Gasterosteus aculeatus) in a heterogeneous natural system,Lake Mývatn,Iceland

Eco‐evolutionary responses of natural populations to spatial environmental variation strongly depend on the relative strength of environmental differences/natural selection and dispersal/gene flow. In absence of geographic barriers, as often is the case in lake ecosystems, gene flow is expected to constrain adaptive divergence between environments – favoring phenotypic plasticity or high trait variability. However, if divergent natural selection is sufficiently strong, adaptive divergence can occur in face of gene flow. The extent of divergence is most often studied between two contrasting environments, whereas potential for multimodal divergence is little explored. We investigated phenotypic (body size, defensive structures, and feeding morphology) and genetic (microsatellites) structure in threespine stickleback (Gasterosteus aculeatus) across five habitat types and two basins (North and South) within the geologically young and highly heterogeneous Lake Mývatn, North East Iceland. We found that (1) North basin stickleback were, on average, larger and had relatively longer spines than South basin stickleback, whereas (2) feeding morphology (gill raker number and gill raker gap width) differed among three of five habitat types, and (3) there was only subtle genetic differentiation across the lake. Overall, our results indicate predator and prey mediated phenotypic divergence across multiple habitats in the lake, in face of gene flow.