Pungitius Pungitius Research Articles

Evolution of stickleback feeding behaviour: genetics of population divergence at different ontogenetic stages

The evolutionary significance of individual consistency in a given behaviour - called animal personality - has been subject to a lot of recent research. However, the genetic underpinnings of population divergence in mean personality have rarely been studied, especially across different ontogenetic stages. Previous work has shown that marine vs. pond populations of nine-spined sticklebacks (Pungitius pungitius) have undergone adaptive divergence in a series of fitness-related traits, including behaviour. One particular behavioural trait important in this system is feeding activity: giant pond sticklebacks are more active feeders than their normal sized marine conspecifics. In a common garden experiment, we raised individuals from pure and hybrid F1 -generation crosses of a highly divergent marine - pond population pair to see if (i) feeding activity and/or its ontogenetic change was consistent between individuals, and if (ii) population divergence at different ontogenetic stages could be explained by additive genetic, nonadditive genetic or maternal effects. We found that feeding activity decreased with age, but that these changes were consistently different among both individuals and crosses. The among cross patterns were consistent with a nonadditive genetic scenario: in the early period pond sticklebacks expressed dominance for high feeding activity, while in the late period marine sticklebacks expressed dominance for low feeding activity. We conclude that nine-spined sticklebacks exhibit different feeding personalities, and that the population divergence in feeding personality is explainable by age-dependent expression of genetic dominance.

Evidence for genetic differentiation in timing of maturation among nine‐spined stickleback populations

Timing of maturation is an important life-history trait that is likely to be subjected to strong natural selection. Although population differences in timing of maturation have been frequently reported in studies of wild animal populations, little is known about the genetic basis of this differentiation. Here, we investigated population and sex differences in timing of maturation within and between two nine-spined stickleback (Pungitius pungitius) populations in a laboratory breeding experiment. We found that fish from the high-predation marine population matured earlier than fish from the low-predation pond population and males matured earlier than females. Timing of maturation in both reciprocal hybrid crosses between the two populations was similar to that in the marine population, suggesting that early timing of maturation is a dominant trait, whereas delayed timing of maturation in the pond is a recessive trait. Thus, the observed population divergence is suggestive of strong natural selection against early maturation in the piscine-predator-free pond population.

Hybridization between two cestode species and its consequences for intermediate host range

BackgroundMany parasites show an extraordinary degree of host specificity, even though a narrow range of host species reduces the likelihood of successful transmission. In this study, we evaluate the genetic basis of host specificity and transmission success of experimental F1 hybrids from two closely related tapeworm species (Schistocephalus solidus and S. pungitii), both highly specific to their respective vertebrate second intermediate hosts (three- and nine-spined sticklebacks, respectively).MethodsWe used an in vitro breeding system to hybridize Schistocephalus solidus and S. pungitii; hybridization rate was quantified using microsatellite markers. We measured several fitness relevant traits in pure lines of the parental parasite species as well as in their hybrids: hatching rates, infection rates in the copepod first host, and infection rates and growth in the two species of stickleback second hosts.ResultsWe show that the parasites can hybridize in the in vitro system, although the proportion of self-fertilized offspring was higher in the heterospecific breeding pairs than in the control pure parental species. Hybrids have a lower hatching rate, but do not show any disadvantages in infection of copepods. In fish, hybrids were able to infect both stickleback species with equal frequency, whereas the pure lines were only able to infect their normal host species.ConclusionsAlthough not yet documented in nature, our study shows that hybridization in Schistocephalus spp. is in principle possible and that, in respect to their expanded host range, the hybrids are fitter. Further studies are needed to find the reason for the maintenance of the species boundaries in wild populations.

Evolutionary Ecological Genomics

The recognition that evolution can happen on ecological timescales (Hairston et al. 2005; Pelletier et al. 2009; Ellner et al. 2011; Becks et al. 2012) has prompted the integration of ecology and evolution, while easier access to high-throughput sequencing technologies has increased the number of genetic nonmodel species entering the ‘omics’ era (e.g. Turner et al. 2010; Colbourne et al. 2011; Jones et al. 2012). We are now in a position to identify the genetic basis of adaptation and the mechanisms of adaptive responses in the wild. It nonetheless remains a challenge to go beyond descriptive measures of patterns of genetic variation and to identify the evolutionary processes driving species adaptation and evolution. This special issue represents a broad cross-section of research into evolutionary adaptation at the genetic level. The approaches used vary from classic QTL studies to RAD sequencing and RNAseq. Given the rapid advance of sequencing technology, we fully expect that ‘genomic’ as defined here will be merely ‘genetic’ in a few years, but we nonetheless hope that the results and methods described in this special issue will serve as a blueprint for future work in this field.

Genomic divergence between nine- and three-spined sticklebacks

BackgroundComparative genomics approaches help to shed light on evolutionary processes that shape differentiation between lineages. The nine-spined stickleback (Pungitius pungitius) is a closely related species of the ecological ‘supermodel’ three-spined stickleback (Gasterosteus aculeatus). It is an emerging model system for evolutionary biology research but has garnered less attention and lacks extensive genomic resources. To expand on these resources and aid the study of sticklebacks in a phylogenetic framework, we characterized nine-spined stickleback transcriptomes from brain and liver using deep sequencing.ResultsWe obtained nearly eight thousand assembled transcripts, of which 3,091 were assigned as putative one-to-one orthologs to genes found in the three-spined stickleback. These sequences were used for evaluating overall differentiation and substitution rates between nine- and three-spined sticklebacks, and to identify genes that are putatively evolving under positive selection. The synonymous substitution rate was estimated to be 7.1 × 10-9 per site per year between the two species, and a total of 165 genes showed patterns of adaptive evolution in one or both species. A few nine-spined stickleback contigs lacked an obvious ortholog in three-spined sticklebacks but were found to match genes in other fish species, suggesting several gene losses within 13 million years since the divergence of the two stickleback species. We identified 47 SNPs in 25 different genes that differentiate pond and marine ecotypes. We also identified 468 microsatellites that could be further developed as genetic markers in nine-spined sticklebacks.ConclusionWith deep sequencing of nine-spined stickleback cDNA libraries, our study provides a significant increase in the number of gene sequences and microsatellite markers for this species, and identifies a number of genes showing patterns of adaptive evolution between nine- and three-spined sticklebacks. We also report several candidate genes that might be involved in differential adaptation between marine and freshwater nine-spined sticklebacks. This study provides a valuable resource for future studies aiming to identify candidate genes underlying ecological adaptation in this and other stickleback species.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-14-756) contains supplementary material, which is available to authorized users.

Optimal growth strategies under divergent predation pressure

The conditions leading to gigantism in nine-spined sticklebacks Pungitius pungitius were analysed by modelling fish growth with the von Bertalanffy model searching for the optimal strategy when the model's growth constant and asymptotic fish size parameters are negatively related to each other. Predator-related mortality was modelled through the increased risk of death during active foraging. The model was parameterized with empirical growth data of fish from four different populations and analysed for optimal growth strategy at different mortality levels. The growth constant and asymptotic fish size were negatively related in most populations. Optimal fish size, fitness and life span decreased with predator-induced mortality. At low mortality, the fitness of pond populations was higher than that of sea populations. The differences disappeared at intermediate mortalities, and sea populations had slightly higher fitness at extremely high mortalities. In the scenario where all populations mature at the same age, the pond populations perform better at low mortalities and the sea populations at high mortalities. It is concluded that a trade-off between growth constant and asymptotic fish size, together with different mortality rates, can explain a significant proportion of body size differentiation between populations. In the present case, it is a sufficient explanation of gigantism in pond P. pungitius.

The evolutionary ecology of dwarfism in three‐spined sticklebacks

1. Body size is a defining phenotypic trait, but the ecological causes of its evolution are poorly understood. Most studies have considered only a single putative causal agent and have failed to recognise that different environmental agents are often correlated. 2. Darwin suggested that although trait variation across populations is often associated with abiotic variation, evolution is more likely to be driven by biotic factors correlated with the abiotic variation. This hypothesis has received little explicit attention. 3. We use structural equation modelling to quantify the relative importance of abiotic (pH, metal concentrations) and biotic (competition, predation) factors in the evolution of body size in three-spined sticklebacks Gasterosteus aculeatus on the island of North Uist, Scotland. We combine phenotypic data from multiple isolated populations, detailed characterisation of their environment and a common garden experiment that establishes the genetic basis of size differences. 4. Three-spined sticklebacks on North Uist show almost unprecedented intraspecific evolution of body size that has taken place rapidly (<16,000 years). The smallest fish mature at only 7% of the mass of ancestral, anadromous fish. Dwarfism is associated with reduced abundance of a smaller competitor species, the nine-spined stickleback Pungitius pungitius, and with low pH indicative of poor resource conditions. Dwarfism also tends to occur where an important predator, the brown trout Salmo trutta, is also small. The abundance of P. pungitius and the size of S. trutta are themselves related to underlying abiotic environmental variation. 5. Despite the close association between abiotic and biotic factors across populations, our results support Darwin's hypothesis that biotic factors, associated with variation in the abiotic environment, are more important in explaining evolution than is abiotic variation per se. This study demonstrates the importance of considering the relationships between environmental variables before conclusions can be drawn about the causes of (body size) evolution on islands.

Factors influencing nine-spined stickleback (Pungitus pungitus)trapping success

Information about factors influencing catch per unit effort (CPUE) are not only interesting in commercial and sport fisheries contexts, but also for scientists fishing for research purposes. Here, the influence of trap type, baiting and attractors on CPUE in pond nine-spined sticklebacks (Pungitius pungitius) were investigated. The results show that baiting improved CPUE considerably, but trap type had a smaller effect on CPUE. Visual attractors (i.e. aluminium foil) did not improve CPUE. Apart from showing that nine-spined sticklebacks are likely to use olfactory cues for locating food, the results should provide useful information for those interested in catching nine-spined sticklebacks with the aid of minnow traps.

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, ••, ••–••.

Estimating population-level coancestry coefficients by an admixture F model.

In this article, we develop an admixture F model (AFM) for the estimation of population-level coancestry coefficients from neutral molecular markers. In contrast to the previously published F model, the AFM enables disentangling small population size and lack of migration as causes of genetic differentiation behind a given level of F(ST). We develop a Bayesian estimation scheme for fitting the AFM to multiallelic data acquired from a number of local populations. We demonstrate the performance of the AFM, using simulated data sets and real data on ninespine sticklebacks (Pungitius pungitius) and common shrews (Sorex araneus). The results show that the parameterization of the AFM conveys more information about the evolutionary history than a simple summary parameter such as F(ST). The methods are implemented in the R package RAFM.

Molecular evolutionary and population genomic analysis of the nine‐spined stickleback using a modified restriction‐site‐associated DNA tag approach

In recent years, the explosion of affordable next generation sequencing technology has provided an unprecedented opportunity to conduct genome-wide studies of adaptive evolution in organisms previously lacking extensive genomic resources. Here, we characterize genome-wide patterns of variability and differentiation using pooled DNA from eight populations of the nine-spined stickleback (Pungitius pungitius L.) from marine, lake and pond environments. We developed a novel genome complexity reduction protocol, defined as paired-end double restriction-site-associated DNA (PE dRAD), to maximize read coverage at sequenced locations. This allowed us to identify over 114,000 short consensus sequences and 15,000 SNPs throughout the genome. A total of 6834 SNPs mapped to a single position on the related three-spined stickleback genome, allowing the detection of genomic regions affected by divergent and balancing selection, both between species and between freshwater and marine populations of the nine-spined stickleback. Gene ontology analysis revealed 15 genomic regions with elevated diversity, enriched for genes involved in functions including immunity, chemical stimulus response, lipid metabolism and signalling pathways. Comparisons of marine and freshwater populations identified nine regions with elevated differentiation related to kidney development, immunity and MAP kinase pathways. In addition, our analysis revealed that a large proportion of the identified SNPs mapping to LG XII is likely to represent alternative alleles from divergent X and Y chromosomes, rather than true autosomal markers following Mendelian segregation. Our work demonstrates how population-wide sequencing and combining inter- and intra-specific RAD analysis can uncover genome-wide patterns of differentiation and adaptations in a non-model species.

Heterozygosity–behaviour correlations in nine‐spined stickleback (Pungitius pungitius) populations: contrasting effects at random and functional loci

The study of heterozygosity-fitness correlations (HFCs) has a long history in the fields of ecology and evolutionary biology but remains controversial. Recently, it has been shown that the genetic distance of markers from functional loci can be an important factor to be considered in addition to marker numbers and variability. In this study, we investigated the correlation between individual heterozygosity and behaviour (aggression, boldness and feeding activity) in nine-spined stickleback (Pungitius pungitius) individuals originating from four populations in two contrasting environments. Offspring of full-sib families raised in a common garden setting were assessed for behaviour and genotyped using 84 microsatellite markers that were either located within or near behaviourally or physiologically important genes (termed 'functional') or were randomly selected. No associations were detected with any behavioural trait in any population or over all populations when genetic variability was measured using all 84 markers combined. However, when the markers were separated into three functional categories (behavioural, physiological and random), several significant associations were observed both with functional markers and random markers in one of the four populations assessed. Interestingly, contrasting correlations with behaviour were observed when using physiological gene (negative) and random (positive) markers. Upon dividing the physiological gene markers into further subcategories based on their specific physiological functions, a strong relationship between the heterozygosity of markers linked to osmoregulation-related genes, and behaviour was revealed in the brackish water population. Our results indicate that both local (physiological) and general (neutral) effects are important in shaping behaviour and that heterozygosity-behaviour correlations are population dependent.

Body size divergence in nine-spined sticklebacks: disentangling additive genetic and maternal effects

Interpopulation differences in body size are of common occurrence in vertebrates, but the relative importance of genetic, maternal, and environmental effects as causes of observed differentiation have seldom been assessed in the wild. Gigantism in pond nine-spined sticklebacks (Pungitius pungitius Linnaeus, 1758) has been repeatedly observed, but the quantitative genetic basis of population divergence in size has remained unstudied. We conducted a common garden experiment – using ‘pure’ and reciprocal crosses between two populations (‘giant’ pond versus ‘normal’ marine) – to test for the relative importance of additive genetic, non-additive genetic, and maternal effects on body size after 11 months of growth in the laboratory. We found that body size difference between the two populations in laboratory conditions owed mainly to additive genetic effects, and only to a minor degree to maternal effects. Furthermore, the weak maternal effects were seen only in the offspring of ‘giant’ mothers, and appeared to be mediated through differences in egg size. Thus, the results suggest that gigantism in pond populations of P. pungitius is based on the effects of additively acting genes, rather than to direct environmental induction, or maternal or non-additive gene action.

Diversifying and correlational selection on behavior toward conspecific and heterospecific competitors in brook stickleback (Culaea inconstans)

Behaviors toward heterospecifics and conspecifics may be correlated because of shared mechanisms of expression in both social contexts (nonadaptive covariation) or because correlational selection favors adaptive covariation. We evaluated these hypotheses by comparing behavior toward conspecifics and heterospecifics in brook stickleback (Culaea inconstans) from three populations sympatric with and three allopatric from a competitor, the ninespine stickleback (Pungitius pungitius). Behavioral traits were classified into three multivariate components: overt aggression, sociability, and activity. The correlation of behavior between social contexts for both overt aggression and activity varied among populations in a way unrelated to sympatry with ninespine stickleback, while mean aggression was reduced in sympatry. Correlations in allopatric populations suggest that overt aggression and activity may genetically covary between social contexts for nonadaptive reasons. Sociability was rarely correlated in allopatry but was consistently correlated in sympatry despite reduced mean sociability, suggesting that correlational selection may favor a sociability syndrome in brook stickleback when they coexist with ninespine stickleback. Thus, interspecific competition may impose diversifying selection on behavior among populations, although the causes of correlated behavior toward conspecifics and heterospecifics and whether it can evolve in one social context independent of the other may depend on the type of behavior.

Spectral tuning by selective chromophore uptake in rods and cones of eight populations of nine-spined stickleback (Pungitius pungitius)

The visual pigments of rods and cones were studied in eight Fennoscandian populations of nine-spined stickleback (Pungitius pungitius). The wavelength of maximum absorbance of the rod pigment (λ(max)) varied between populations from 504 to 530 nm. Gene sequencing showed that the rod opsins of all populations were identical in amino acid composition, implying that the differences were due to varying proportions of chromophores A1 and A2. Four spectral classes of cones were found (two S-cones, M-cones and L-cones), correlating with the four classes of vertebrate cone pigments. For quantitative estimation of chromophore proportions, we considered mainly rods and M-cones. In four populations, spectra of both photoreceptor types indicated A2 dominance (population mean λ(max)=525-530 nm for rods and 535-544 nm for M-cones). In the four remaining populations, however, rod spectra (mean λ(max)=504-511 nm) indicated strong A1 dominance, whereas M-cone spectra (mean λ(max)=519-534 nm) suggested substantial fractions of A2. Quantitative analysis of spectra by three methods confirmed that rods and cones in these populations use significantly different chromophore proportions. The outcome is a shift of M-cone spectra towards longer wavelengths and a better match to the photic environment (light spectra peaking >560 nm in all the habitats) than would result from the chromophore proportions of the rods. Chromophore content was also observed to vary partly independently in M- and L-cones with potential consequences for colour discrimination. This is the first demonstration that selective processing of chromophore in rods and cones, and in different cone types, may be ecologically relevant.

High levels of fluctuating asymmetry in isolated stickleback populations

BackgroundFluctuating asymmetry (FA), defined as small random deviations from the ideal bilateral symmetry, has been hypothesized to increase in response to both genetic and environmental stress experienced by a population. We compared levels of FA in 12 bilateral meristic traits (viz. lateral-line system neuromasts and lateral plates), and heterozygosity in 23 microsatellite loci, among four marine (high piscine predation risk) and four pond (zero piscine predation risk) populations of nine-spined sticklebacks (Pungitius pungitius).ResultsPond sticklebacks had on average three times higher levels of FA than marine fish and this difference was highly significant. Heterozygosity in microsatellite markers was on average two times lower in pond (HE ≈ 0.3) than in marine (HE ≈ 0.6) populations, and levels of FA and heterozygosity were negatively correlated across populations. However, after controlling for habitat effect on heterozygosity, levels of FA and heterozygosity were uncorrelated.ConclusionsThe fact that levels of FA in traits likely to be important in the context of predator evasion were elevated in ponds compared to marine populations suggests that relaxed selection for homeostasis in ponds lacking predatory fish may be responsible for the observed habitat difference in levels of FA. This inference also aligns with the observation that the levels of genetic variability across the populations did not explain population differences in levels of FA after correcting for habitat effect. Hence, while differences in strength of selection, rather than in the degree of genetic stress could be argued to explain habitat differences in levels of FA, the hypothesis that increased FA in ponds is caused by genetic stress cannot be rejected.

Habitat use by fishes of Lake Superior. I. Diel patterns of habitat use in nearshore and offshore waters of the Apostle Islands region

Diel patterns of distribution of fishes in nearshore (15–80 m depth) and offshore (&amp;gt;80 m) waters of the Apostle Islands region of Lake Superior were described using bottom trawls, mid-water trawls, and acoustic gear during day and night sampling. These data revealed three types of diel migration: diel vertical migration (DVM), diel bank migration (DBM), and no migration. DVM was expressed by fishes migrating from benthopelagic to pelagic strata and DBM was expressed by fishes migrating horizontally from deeper waters in the day to shallower waters at night while remaining within the benthopelagic stratum. Most fishes that did not exhibit diel migration showed increased nighttime densities as a result of increased activity and movement from benthic to benthopelagic strata. Rainbow Smelt (Osmerus mordax), Cisco (Coregonus artedi), Bloater (C. hoyi), Kiyi (C. kiyi), juvenile Trout-Perch (Percopsis omiscomaycus), and adult siscowet (Salvelinus namaycush siscowet) exhibited DVM. Lake Whitefish (C. clupeaformis), lean Lake Trout (Salvelinus namaycush namaycush), and juvenile siscowet exhibited DBM. Adult Trout-Perch and adult Pygmy Whitefish (Prosopium coulteri) exhibited a mixture of DBM and DVM. Burbot (Lota lota), Slimy Sculpin (Cottus cognatus), Spoonhead Sculpin (C. ricei), and Deepwater Sculpin (Myoxocephalus thompsonii) did not exhibit diel migration, but showed evidence of increased nocturnal activity. Ninespine Stickleback (Pungitius pungitius) exhibited a mixture of DVM and non-migration. Juvenile Pygmy Whitefish did not show a diel change in density or depth distribution. Species showing ontogenetic shifts in depth distribution with larger, adult life stages occupying deeper waters included, Rainbow Smelt, lean and siscowet Lake Trout, Lake Whitefish, Pygmy Whitefish, Ninespine Stickleback and Trout-Perch. Of these species, siscowet also showed an ontogenetic shift from primarily DBM as juveniles to primarily DVM as adults. Across all depths, fishes expressing DVM accounted for 73% of the total estimated community areal biomass (kg ha−1) while those expressing DBM accounted for 25% and non-migratory species represented 2% of the biomass. The proportion of total community biomass exhibiting DVM increased with depth, from 59% to 95% across ≤30 m to &amp;gt;90 m depth zones. Along the same depth gradient, the proportion of total community biomass exhibiting DBM declined from 40% to 1%, while non-migrators increased from 1% to 4%. These results indicate that DVM and DBM behaviors are pervasive in the Lake Superior fish community and potentially provide strong linkages that effect coupling of benthic and pelagic and nearshore and offshore habitats.

Contrasting population structures in two sympatric fishes in the Baltic Sea basin

Detailed multispecies studies on the patterns of genetic variability and differentiation in marine environments are still rare. Using mitochondrial and nuclear genetic markers, we compared genetic variability and population structuring of threespine (Gasterosteus aculeatus) and ninespine (Pungitius pungitius) sticklebacks from the same eleven marine and six freshwater locations within the Baltic Sea basin. Analyses of both marker types revealed a significantly lower degree of genetic structuring in both marine and freshwater populations of threespine than those ninespine sticklebacks. Isolation-by-distance (IBD) was detected across the marine populations in both species, suggesting spatially limited gene flow. However, the levels of genetic diversity and differentiation across the localities were uncorrelated between the two species in both marine and freshwater environments. Accordingly, estimates of effective population sizes were larger and migration rates were higher for three- than for ninespine sticklebacks. Hence, ninespine stickleback populations from the Baltic Sea basin appear to be subject to stronger genetic drift than sympatric threespine sticklebacks, and the proximate reason for this difference is likely to be found from autecological differences between the two species. In accordance with the earlier studies, genetic variability was higher and the degree of genetic differentiation was lower in marine than in freshwater populations in both species.

On nine-spined sticklebacks of the genus Pungitius (Gasterosteiformes) of Lake Bol’shoe (Bol’shoi Shantar Island, Shantar Islands)

Specific features of morphology and biology of two species of nine-spined sticklebacks—Pungitius sinensis and P. pungitius—dwelling in brackish Lake Bol’shoe (Bol’shoi Shantar Island) have been studied. Differences in the number of dorsal spines, vertebrae, and coloration, as well as the absence of specimens with an intermediate phenotype indicate the reproductive isolation of two species supported by the temporal separation of spawning periods. Species P. pungitius was found in Lake Bol’shoe for the first time.

Niche segregation in two closely related species of stickleback along a physiological axis: explaining multidecadal changes in fish distribution from iron-induced respiratory impairment.

Acute exposure to iron can be lethal to fish, but long-term sublethal impacts of iron require further study. Here we investigated whether the spatial and temporal distribution (1967–2004) of two closely related species of stickleback matched the spatial distribution of iron concentrations in the groundwater. We used the ‘Northern Peel region’, a historically iron-rich peat landscape in The Netherlands as a case study. This allowed us to test the hypothesis that niche segregation in two closely related species of stickleback occurred along a physiological axis. Patterns in stickleback occurrence were strongly associated with spatial patterns in iron concentrations before 1979: iron-rich grid cells were avoided by three-spined stickleback (Gasterosteus aculeatus, Linnaeus 1758) and preferred by nine-spined stickleback (Pungitius pungitius, [Linnaeus, 1758]). After 1979, the separation between both sticklebacks became weaker, corresponding to a decreased influence of local groundwater on stream water quality. The way both species changed their distribution in the field provides a strong indication that they differ in their susceptibility to iron-rich conditions. These observed differences correspond with differences in their respiration physiology, tolerance of poor oxygen conditions and overall life-history strategy documented in the literature. Our results exemplify how species can partition niche along a non-structural niche axis, such as sublethal iron-rich conditions. Other fish species may similarly segregate along concentration gradients in iron, while sublethal concentrations of other metals such as copper may similarly impact fish via respiratory impairment and reduced aerobic scope.Electronic supplementary materialThe online version of this article (doi:10.1007/s10452-012-9395-y) contains supplementary material, which is available to authorized users.