Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake–stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecological speciation. Ecological speciation is well studied in parts of the native range of the three‐spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving different stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajectories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known ecological relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.     

PHENOTYPE-DEPENDENT NATIVE HABITAT PREFERENCE FACILITATES DIVERGENCE BETWEEN PARAPATRIC LAKE AND STREAM STICKLEBACK

Adaptive divergence between adjoining populations reflects a balance between the diversifying effect of divergent selection and the potentially homogenizing effect of gene flow. In most models of migration-selection balance, gene flow is assumed to reflect individuals' inherent capacity to disperse, without regard to the match between individuals' phenotypes and the available habitats. However, habitat preferences can reduce dispersal between contrasting habitats, thereby alleviating migration load and facilitating adaptive divergence. We tested whether habitat preferences contribute to adaptive divergence in a classic example of migration-selection balance: parapatric lake and stream populations of three-spine stickleback ( Gasterosteus aculeatus ). Using a mark-transplant-recapture experiment on morphologically divergent parapatric populations, we showed that 90% of lake and stream stickleback returned to their native habitat, reducing migration between habitats by 76%. Furthermore, we found that dispersal into a nonnative habitat was phenotype dependent. Stream fish moving into the lake were morphologically more lake-like than those returning to the stream (and the converse for lake fish entering the stream). The strong native habitat preference documented here increases the extent of adaptive divergence between populations two- to fivefold relative to expectations with random movement. These results illustrate the potential importance of adaptive habitat choice in driving parapatric divergence.     

Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so‐called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral‐like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low‐dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.     

Population genetic differences along a latitudinal cline between original and recently colonized habitat in a butterfly

BACKGROUND: Past and current range or spatial expansions have important consequences on population genetic structure. Habitat-use expansion, i.e. changing habitat associations, may also influence genetic population parameters, but has been less studied. Here we examined the genetic population structure of a Palaeartic woodland butterfly Pararge aegeria (Nymphalidae) which has recently colonized agricultural landscapes in NW-Europe. Butterflies from woodland and agricultural landscapes differ in several phenotypic traits (including morphology, behavior and life history). We investigated whether phenotypic divergence is accompanied by genetic divergence between populations of different landscapes along a 700 km latitudinal gradient. METHODOLOGY/PRINCIPAL FINDINGS: Populations (23) along the latitudinal gradient in both landscape types were analyzed using microsatellite and allozyme markers. A general decrease in genetic diversity with latitude was detected, likely due to post-glacial colonization effects. Contrary to expectations, agricultural landscapes were not less diverse and no significant bottlenecks were detected. Nonetheless, a genetic signature of recent colonization is reflected in the absence of clinal genetic differentiation within the agricultural landscape, significantly lower gene flow between agricultural populations (3.494) than between woodland populations (4.183), and significantly higher genetic differentiation between agricultural (0.050) than woodland (0.034) pairwise comparisons, likely due to multiple founder events. Globally, the genetic data suggest multiple long distance dispersal/colonization events and subsequent high intra- and inter-landscape gene flow in this species. Phosphoglucomutase deviated from other enzymes and microsatellite markers, and hence may be under selection along the latitudinal gradient but not between landscape types. Phenotypic divergence was greater than genetic divergence, indicating directional selection on some flight morphology traits. MAIN CONCLUSIONS/SIGNIFICANCE: Clinal differentiation characterizes the population structure within the original woodland habitat. Genetic signatures of recent habitat expansion remain, notwithstanding high gene flow. After differentiation through drift was excluded, both latitude and landscape were significant factors inducing spatially variable phenotypic variation.     

Environmentally independent selection for hybrids between divergent freshwater stickleback lineages in semi-natural ponds

Hybridization following secondary contact of genetically divergent populations can influence the range expansion of invasive species, though specific outcomes depend on the environmental dependence of hybrid fitness. Here, using two genetically and ecologically divergent threespine stickleback lineages that differ in their history of freshwater colonization, we estimate fitness variation of parental lineages and hybrids in semi-natural freshwater ponds with contrasting histories of nutrient loading. In our experiment, we found that fish from the older freshwater lineage (Lake Geneva) and hybrids outperformed fish from the younger freshwater lineage (Lake Constance) in terms of both growth and survival, regardless of the environmental context of our ponds. Across all ponds, hybrids exhibited the highest survival. Although wild-caught adult populations differed in their functional and defence morphology, it is unclear which of these traits underlie the fitness differences observed among juveniles in our experiment. Overall, our work suggests that when hybrid fitness is insensitive to environmental conditions, as observed here, introgression may promote population expansion into unoccupied habitats and accelerate invasion success.     

Differences in rheotactic responses contribute to divergent habitat use between parapatric lake and stream threespine stickleback

Migration among populations is widely thought to undermine adaptive divergence, assuming gene flow arises from random movement of individuals. If individuals instead differ in dispersal behavior, phenotype‐dependent dispersal can reduce the effective rate of gene flow or even facilitate divergence. For example, parapatric populations of lake and stream stickleback tend to actively avoid dispersing into the adjoining habitat. However, the behavioral basis of this nonrandom dispersal was previously unknown. Here, we show that lake and stream stickleback exhibit divergent rheotactic responses (behavioral response to currents). During the breeding season, wild‐caught inlet stream stickleback were better than lake fish at maintaining position in currents, faced upstream more, and spent more time in low‐current areas. As a result, stream fish expended significantly less energy in currents than did lake fish. These divergent rheotactic responses likely contribute to divergent habitat use by lake and stream stickleback. Although rheotactic differences were absent in nonbreeding fish, divergent behavior of breeding‐season fish may suffice for assortative mating by breeding location. The resulting reproductive isolation between lake and stream fish may explain the fine‐scale evolutionary differentiation in parapatric stickleback populations.     

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.     

Parallel and nonparallel ecological,morphological and genetic divergence in lake–stream stickleback from a single catchment

Parallel phenotypic evolution in similar environments has been well studied in evolutionary biology; however, comparatively little is known about the influence of determinism and historical contingency on the nature, extent and generality of this divergence. Taking advantage of a novel system containing multiple lake–stream stickleback populations, we examined the extent of ecological, morphological and genetic divergence between three‐spined stickleback present in parapatric environments. Consistent with other lake–stream studies, we found a shift towards a deeper body and shorter gill rakers in stream fish. Morphological shifts were concurrent with changes in diet, indicated by both stable isotope and stomach contents analysis. Performing a multivariate test for shared and unique components of evolutionary response to the distance gradient from the lake, we found a strong signature of parallel adaptation. Nonparallel divergence was also present, attributable mainly to differences between river locations. We additionally found evidence of genetic substructuring across five lake–stream transitions, indicating that some level of reproductive isolation occurs between populations in these habitats. Strong correlations between pairwise measures of morphological, ecological and genetic distance between lake and stream populations supports the hypothesis that divergent natural selection between habitats drives adaptive divergence and reproductive isolation. Lake–stream stickleback divergence in Lough Neagh provides evidence for the deterministic role of selection and supports the hypothesis that parallel selection in similar environments may initiate parallel speciation.     

Niche divergence facilitated by fine‐scale ecological partitioning in a recent cichlid fish adaptive radiation

Ecomorphological differentiation is a key feature of adaptive radiations, with a general trend for specialization and niche expansion following divergence. Ecological opportunity afforded by invasion of a new habitat is thought to act as an ecological release, facilitating divergence, and speciation. Here, we investigate trophic adaptive morphology and ecology of an endemic clade of oreochromine cichlid fishes (Alcolapia) that radiated along a herbivorous trophic axis following colonization of an isolated lacustrine environment, and demonstrate phenotype‐environment correlation. Ecological and morphological divergence of the Alcolapia species flock are examined in a phylogenomic context, to infer ecological niche occupation within the radiation. Species divergence is observed in both ecology and morphology, supporting the importance of ecological speciation within the radiation. Comparison with an outgroup taxon reveals large‐scale ecomorphological divergence but shallow genomic differentiation within the Alcolapia adaptive radiation. Ancestral morphological reconstruction suggests lake colonization by a generalist oreochromine phenotype that diverged in Lake Natron to varied herbivorous morphologies akin to specialist herbivores in Lakes Tanganyika and Malawi.     

DISENTANGLING THE ROLE OF PHENOTYPIC PLASTICITY AND GENETIC DIVERGENCE IN CONTEMPORARY ECOTYPE FORMATION DURING A BIOLOGICAL INVASION

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full‐sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.     

Contrasting geographic structure in evolutionarily divergent Lake Tanganyika catfishes

Geographic isolation is suggested to be among the most important processes in the generation of cichlid fish diversity in East Africa's Great Lakes, both through isolation by distance and fluctuating connectivity caused by changing lake levels. However, even broad scale phylogeographic patterns are currently unknown in many non‐cichlid littoral taxa from these systems. To begin to address this, we generated restriction‐site‐associated DNA sequence (RADseq) data to investigate phylogeographic structure throughout Lake Tanganyika (LT) in two broadly sympatric rocky shore catfish species from independent evolutionary radiations with differing behaviors: the mouthbrooding claroteine, Lophiobagrus cyclurus, and the brood‐parasite mochokid, Synodontis multipunctatus. Our results indicated contrasting patterns between these species, with strong lake‐wide phylogeographic signal in L. cyclurus including a deep divergence between the northern and southern lake basins. Further structuring of these populations was observed across a heterogeneous habitat over much smaller distances. Strong population growth was observed in L. cyclurus sampled from shallow shorelines, suggesting population growth associated with the colonization of new habitats following lake‐level rises. Conversely, S. multipunctatus, which occupies a broader depth range, showed little phylogeographic structure and lower rates of population growth. Our findings suggest that isolation by distance and/or habitat barriers may play a role in the divergence of non‐cichlid fishes in LT, but this effect varies by species.     

Parallelism and historical contingency during rapid ecotype divergence in an isopod

Recent studies on parallel evolution have focused on the relative role of selection and historical contingency during adaptive divergence. Here, we study geographically separate and genetically independent lake populations of a freshwater isopod (Asellus aquaticus) in southern Sweden. In two of these lakes, a novel habitat was rapidly colonized by isopods from a source habitat. Rapid phenotypic changes in pigmentation, size and sexual behaviour have occurred, presumably in response to different predatory regimes. We partitioned the phenotypic variation arising from habitat ('selection': 81–94%), lake ('history': 0.1–6%) and lake × habitat interaction ('unique diversification': 0.4–13%) for several traits. There was a limited role for historical contingency but a strong signature of selection. We also found higher phenotypic variation in the source populations. Phenotype sorting during colonization and strong divergent selection might have contributed to these rapid changes. Consequently, phenotypic divergence was only weakly influenced by historical contingency.     

Distinct colonization waves underlie the diversification of the freshwater sculpin (Cottus gobio) in the Central European Alpine region

Ecological speciation and adaptive radiation are key processes shaping northern temperate freshwater fish diversity. Both often involve parapatric differentiation between stream and lake populations and less often, sympatric intralacustrine diversification into habitat‐ and resource‐associated ecotypes. However, few taxa have been studied, calling for studies of others to investigate the generality of these processes. Here, we test for diversification within catchments in freshwater sculpins in a network of peri‐Alpine lakes and streams. Using 8047 and 13 182 restriction site‐associated (RADseq) SNPs, respectively, we identify three deeply divergent phylogeographic lineages associated with different major European drainages. Within the Aare catchment, we observe populations from geographically distant lakes to be genetically more similar to each other than to populations from nearby streams. This pattern is consistent with two distinct colonization waves, rather than by parapatric ecological speciation after a single colonization wave. We further find two distinct depth distribution modes in three lakes of the Aare catchment, one in very shallow and one in very deep water, and significant genomewide differentiation between these in one lake. Sculpins in the Aare catchment appear to represent an early‐stage adaptive radiation involving the evolution of a lacustrine lineage distinct from parapatric stream sculpins and the repeated onset of depth‐related intralacustrine differentiation.     

Introgressive hybridization erodes morphological divergence between lentic and lotic habitats in an endangered minnow

Introgressive hybridization may erode phenotypic divergence along environmental gradients, collapsing locally adapted populations into a hybrid swarm. Alternatively, introgression may promote phenotypic divergence by providing variation on which natural selection can act. In freshwater fishes, water flow often selects for divergent morphological traits in lake versus stream habitats. We tested the effects of introgression on lake–stream morphological divergence in the minnow Owens Tui Chub (Siphateles bicolor snyderi), which has been rendered endangered by introgession from the introduced Lahontan Tui Chub (Siphateles bicolor obesa). Using geometric morphometric analysis of 457 individual Tui Chub from thirteen populations, we found that both native and introgressing parent taxa exhibited divergent body and caudal fin shapes in lake versus stream habitats, but their trajectories of divergence were distinct. In contrast, introgressed populations exhibited intermediate body and caudal fin shapes that were not differentiated by habitat type, indicating that introgression has eroded phenotypic divergence along the lentic–lotic gradient throughout the historic range of the Owens Tui Chub. Individuals within hybrid populations were less morphologically variable than those within parent populations, suggesting hybrid adaptation to selective agents other than water flow or loss of variance by drift.     

Contemporary Parallel Diversification,Antipredator Adaptations and Phenotypic Integration in an Aquatic Isopod

It is increasingly being recognized that predation can be a strong diversifying agent promoting ecological divergence. Adaptations against different predatory regimes can emerge over short periods of time and include many different traits. We studied antipredator adaptations in two ecotypes of an isopod (Asellus aquaticus) that have, diverged in parallel in two Swedish lakes over the last two decades. We quantified differences in escape speed, morphology and behavior for isopods from different ecotypes present in these lakes. Isopods from the source habitat (reed) coexist with mainly invertebrate predators. They are more stream-profiled and have higher escape speeds than isopods in the newly colonized stonewort habitat, which has higher density of fish predators. Stonewort isopods also show more cautious behaviors and had higher levels of phenotypic integration between coloration and morphological traits than the reed isopods. Colonization of a novel habitat with a different predation regime has thus strengthened the correlations between pigmentation and morphology and weakened escape performance. The strong signature of parallelism for these phenotypic traits indicates that divergence is likely to be adaptive and is likely to have been driven by differences in predatory regimes. Furthermore, our results indicate that physical performance, behavior and morphology can change rapidly and in concert as new habitats are colonized.     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.     

Recent local adaptation of sockeye salmon to glacial spawning habitats

Salmonids spawn in highly diverse habitats, exhibit strong genetic population structuring, and can quickly colonize newly created habitats with few founders. Spawning traits often differ among populations, but it is largely unknown if these differences are adaptive or due to genetic drift. To test if sockeye salmon (Oncorhynchus nerka) populations are adapted to glacial, beach, and tributary spawning habitats, we examined variation in heritable phenotypic traits associated with spawning in 13 populations of wild sockeye salmon in Lake Clark, Alaska. These populations were commonly founded between 100 and 400 hundred sockeye salmon generations ago and exhibit low genetic divergence at 11 microsatellite loci (F _ST < 0.024) that is uncorrelated with spawning habitat type. We found that mean P _ST (phenotypic divergence among populations) exceeded neutral F _ST for most phenotypic traits measured, indicating that phenotypic differences among populations could not be explained by genetic drift alone. Phenotypic divergence among populations was associated with spawning habitat differences, but not with neutral genetic divergence. For example, female body color was lighter and egg color was darker in glacial than non-glacial habitats. This may be due to reduced sexual selection for red spawning color in glacial habitats and an apparent trade-off in carotenoid allocation to body and egg color in females. Phenotypic plasticity is an unlikely source of phenotypic differences because Lake Clark sockeye salmon spend nearly all their lives in a common environment. Our data suggest that Lake Clark sockeye salmon populations are adapted to spawning in glacial, beach and tributary habitats and provide the first evidence of a glacial spawning ecotype in salmonids. Glacial spawning habitats are often young (i.e., <200 years old) and ephemeral. Thus, local adaptation of sockeye salmon to glacial habitats appears to have occurred recently.     

Quantifying the constraining influence of gene flow on adaptive divergence in the lake-stream threespine stickleback system

The constraining effect of gene flow on adaptive divergence is often inferred but rarely quantified. We illustrate ways of doing so using stream populations of threespine stickleback (Gasterosteus aculeatus) that experience different levels of gene flow from a parapatric lake population. In the Misty Lake watershed (British Columbia, Canada), the inlet stream population is morphologically divergent from the lake population, and presumably experiences little gene flow from the lake. The outlet stream population, however, shows an intermediate phenotype and may experience more gene flow from the lake. We first used microsatellite data to demonstrate that gene flow from the lake is low into the inlet but high into the outlet, and that gene flow from the lake remains relatively constant with distance along the outlet. We next combined gene flow data with morphological and habitat data to quantify the effect of gene flow on morphological divergence. In one approach, we assumed that inlet stickleback manifest well-adapted phenotypic trait values not constrained by gene flow. We then calculated the deviation between the observed and expected phenotypes for a given habitat in the outlet. In a second approach, we parameterized a quantitative genetic model of adaptive divergence. Both approaches suggest a large impact of gene flow, constraining adaptation by 80-86% in the outlet (i.e., only 14-20% of the expected morphological divergence in the absence of gene flow was observed). Such approaches may be useful in other taxa to estimate how important gene flow is in constraining adaptive divergence in nature.     

Molecular and morphological data suggest weak phylogeographic structure in the fairy shrimp <Emphasis Type="Italic">Streptocephalus torvicornis</Emphasis> (Branchiopoda,Anostraca)

Lake Tanganyika has undergone substantial climate-driven lake level fluctuations that have repeatedly changed the distribution and extent of habitat for endemic fishes. Here we consider whether patterns of population genetic structure and phenotypic divergence within the cichlid fish Telmatochromis temporalis have been affected by changing lake levels. The species has a large-bodied rock-living ecomorph and a small-bodied shell-living ecomorph, and both are found in close proximity in littoral habitats. Using mtDNA sequences we found that geographically distant (>50 km) populations within the southern lake region diverged approximately 130,000–230,000 years ago, suggesting that the regional genetic structure persisted through a low stand of over 400 m ~106,000 years ago that ended with a rise to present levels ~100,000 years ago. We also found signatures of large population expansions since this rise across the study region, suggesting that the populations positively responded to new habitat as lake levels rose to present levels. Finally, we found that geographically adjacent (<10 km) ecomorphs exhibit both significant genetic differentiation and signatures of gene flow after the lake level rise. The results suggest that local ecomorph divergence progressed with gene flow after the last major rise in lake level ~100,000, potentially facilitated by new ecological opportunities.     

Incipient sympatric speciation in Midas cichlid fish from the youngest and one of the smallest crater lakes in Nicaragua due to differential use of the benthic and limnetic habitats?

Understanding how speciation can occur without geographic isolation remains a central objective in evolutionary biology. Generally, some form of disruptive selection and assortative mating are necessary for sympatric speciation to occur. Disruptive selection can arise from intraspecific competition for resources. If this competition leads to the differential use of habitats and variation in relevant traits is genetically determined, then assortative mating can be an automatic consequence (i.e., habitat isolation). In this study, we caught Midas cichlid fish from the limnetic (middle of the lake) and benthic (shore) habitats of Crater Lake Asososca Managua to test whether some of the necessary conditions for sympatric speciation due to intraspecific competition and habitat isolation are given. Lake As. Managua is very small (<900 m in diameter), extremely young (maximally 1245 years of age), and completely isolated. It is inhabited by, probably, only a single endemic species of Midas cichlids, Amphilophus tolteca. We found that fish from the limnetic habitat were more elongated than fish collected from the benthic habitat, as would be predicted from ecomorphological considerations. Stable isotope analyses confirmed that the former also exhibit a more limnetic lifestyle than the latter. Furthermore, split‐brood design experiments in the laboratory suggest that phenotypic plasticity is unlikely to explain much of the observed differences in body elongation that we observed in the field. Yet, neutral markers (microsatellites) did not reveal any genetic clustering in the population. Interestingly, demographic inferences based on RAD‐seq data suggest that the apparent lack of genetic differentiation at neutral markers could simply be due to a lack of time, as intraspecific competition may only have begun a few hundred generations ago.