Contrasting growth strategies of pond versus marine populations of nine-spined stickleback (Pungitius pungitius): a combined effect of predation and competition?

Sexual dimorphism in subterranean amphipod crustaceans covaries with subterranean habitat type

Sexual dimorphism can evolve in response to sex-specific selection pressures that vary across habitats. We studied sexual differences in subterranean amphipods Niphargus living in shallow subterranean habitats (close to the surface), cave streams (intermediate), and cave lakes (deepest and most isolated). These three habitats differ because at greater depths there is lower food availability, reduced predation, and weaker seasonality. Additionally, species near the surface have a near-even adult sex ratio (ASR), whereas species from cave lakes have a female-biased ASR. We hypothesized (a) a decrease in sexual dimorphism from shallow subterranean habitats to cave lake species because of weaker sexual selection derived from changes in the ASR and (b) an increase in female body size in cave lakes because of stronger fecundity selection on account of oligotrophy, reduced predation, and weaker seasonality. We measured body size and two sexually dimorphic abdominal appendages for all 31 species and several behaviours related to male competition (activity, risk-taking, exploration) for 12 species. Species with an equal ASR that live close to the surface exhibited sexual dimorphism in all three morphological traits, but not in behaviour. The body size of females increased from the surface to cave lakes, but no such trend was observed in males. In cave lake species, males and females differed neither morphologically nor behaviourally. Our results are consistent with the possibility that sexual and fecundity selection covary across the three habitats, which indirectly and directly, respectively, shape the degree of sexual dimorphism in Niphargus species.

Investigating effects of climate-induced changes in water temperature and diet on mercury concentrations in an Arctic freshwater forage fish

The amount of mercury (Hg) in Arctic lake food webs is, and will continue to be, affected by rapid, ongoing climate change. At warmer temperatures, fish require more energy to sustain growth; changes in their metabolic rates and consuming prey with potentially higher Hg concentrations could result in increased Hg accumulation. To examine the potential implications of climate warming on forage fish Hg accumulation in Arctic lakes, we quantified growth and Hg accumulation in Ninespine Stickleback Pungitius pungitius under different temperature and diet scenarios using bioenergetics models. Four scenarios were considered that examined the role of climate, diet, climate × diet, and climate × diet × elevated prey Hg. As expected, annual fish growth increased with warmer temperatures, but growth rates and Hg accumulation were largely diet dependent. Compared to current growth rates of 0.3 g⋅y^-1, fish growth increased at least 200% for fish consuming energy-dense benthic prey and decreased at least 40% for fish consuming pelagic prey. Compared to baseline levels, the Hg burden per kilocalorie of Ninespine Stickleback declined up to 43% with benthic consumption - indicating strong somatic growth dilution - but no more than 4% with pelagic consumption; elevated prey Hg concentrations led to moderate Hg declines in benthic-foraging fish and Hg increases in pelagic-foraging fish. Bioenergetics models demonstrated the complex interaction of water temperature, growth, prey proportions, and prey Hg concentrations that respond to climate change. Further work is needed to resolve mechanisms and rates linking climate change to Hg availability and uptake in Arctic freshwater systems.

Age-dependent genetic architecture across ontogeny of body size in sticklebacks

Heritable variation in traits under natural selection is a prerequisite for evolutionary response. While it is recognized that trait heritability may vary spatially and temporally depending on which environmental conditions traits are expressed under, less is known about the possibility that genetic variance contributing to the expected selection response in a given trait may vary at different stages of ontogeny. Specifically, whether different loci underlie the expression of a trait throughout development and thus providing an additional source of variation for selection to act on in the wild, is unclear. Here we show that body size, an important life-history trait, is heritable throughout ontogeny in the nine-spined stickleback (Pungitius pungitius). Nevertheless, both analyses of quantitative trait loci and genetic correlations across ages show that different chromosomes/loci contribute to this heritability in different ontogenic time-points. This suggests that body size can respond to selection at different stages of ontogeny but that this response is determined by different loci at different points of development. Hence, our study provides important results regarding our understanding of the genetics of ontogeny and opens an interesting avenue of research for studying age-specific genetic architecture as a source of non-parallel evolution.

Intraspecific variation in behaviour and ecology in a territorial agamid, Ctenophorus fionni

Intraspecific variation as a way to explore factors affecting the evolution of species traits in natural environments is well documented, and also important in the context of preserving biodiversity. In this study, we investigated the extent of behavioural, morphological and ecological variation in the peninsula dragon (Ctenophorus fionni), an endemic Australian agamid that displays extensive variation in colour across three allopatric populations. The aims of the study were to quantify variation across the different populations in terms of the environment, morphometric characteristics and behaviour. We found population level differences in habitat structure and encounter rates. Adult body size of C. fionni, as well as a range of morphometric traits, differed between populations, as well as the frequency of social interactions, which appears to be related to population density and abundance. Analysis of communicative signals showed differences between the southern and central populations, which appear consistent with variations in response to environmental differences between study sites. The findings of the present study, coupled with previous work examining colour variation in this species, show that the three populations of C. fionni have likely undergone substantial differentiation, and would make an interesting study system to explore trait variation in more detail.

Accurate estimates of age at maturity from the growth trajectories of fishes and other ectotherms

Age at maturity (AAM) is a key life history trait that provides insight into ecology, evolution, and population dynamics. However, maturity data can be costly to collect or may not be available. Life history theory suggests that growth is biphasic for many organisms, with a change-point in growth occurring at maturity. If so, then it should be possible to use a biphasic growth model to estimate AAM from growth data. To test this prediction, we used the Lester biphasic growth model in a likelihood profiling framework to estimate AAM from length at age data. We fit our model to simulated growth trajectories to determine minimum data requirements (in terms of sample size, precision in length at age, and the cost to somatic growth of maturity) for accurate AAM estimates. We then applied our method to a large walleye Sander vitreus data set and show that our AAM estimates are in close agreement with conventional estimates when our model fits well. Finally, we highlight the potential of our method by applying it to length at age data for a variety of ectotherms. Our method shows promise as a tool for estimating AAM and other life history traits from contemporary and historical samples.

Population divergence in compensatory growth responses and their costs in sticklebacks

Compensatory growth (CG) may be an adaptive mechanism that helps to restore an organisms’ growth trajectory and adult size from deviations caused by early life resource limitation. Yet, few studies have investigated the genetic basis of CG potential and existence of genetically based population differentiation in CG potential. We studied population differentiation, genetic basis, and costs of CG potential in nine-spined sticklebacks (Pungitius pungitius) differing in their normal growth patterns. As selection favors large body size in pond and small body size in marine populations, we expected CG to occur in the pond but not in the marine population. By manipulating feeding conditions (viz. high, low and recovery feeding treatments), we found clear evidence for CG in the pond but not in the marine population, as well as evidence for catch-up growth (i.e., size compensation without growth acceleration) in both populations. In the marine population, overcompensation occurred individuals from the recovery treatment grew eventually larger than those from the high feeding treatment. In both populations, the recovery feeding treatment reduced maturation probability. The recovery feeding treatment also reduced survival probability in the marine but not in the pond population. Analysis of interpopulation hybrids further suggested that both genetic and maternal effects contributed to the population differences in CG. Hence, apart from demonstrating intrinsic costs for recovery growth, both genetic and maternal effects were identified to be important modulators of CG responses. The results provide an evidence for adaptive differentiation in recovery growth potential.

Geographic variation in age structure and longevity in the nine-spined stickleback (Pungitius pungitius)

Variation in age and size of mature nine-spined sticklebacks (Pungitius pungitius) within and among 16 Fennoscandian populations were assessed using skeletochronology. The average age of individuals in a given population varied from 1.7 to 4.7 years. Fish from pond populations were on average older than those from lake and marine populations, and females tended to be older than males. Reproduction in marine and lake populations commenced typically at an age of two years, whereas that in ponds at an age of three years. The maximum life span of the fish varied from 3 to 7 years. Mean body size within and among populations increased with increasing age, but the habitat and population differences in body size persisted even after accounting for variation in population age (and sex) structure. Hence, the population differences in mean body size are not explainable by age differences alone. As such, much of the pronounced intraspecific variation in population age structure can be attributed to delayed maturation and extended longevity of the pond fish. The results are contrasted and discussed in the context of similar data from the three-spined stickleback (Gasterosteus aculeatus) occupying the same geographic area.

Extensive behavioural divergence following colonisation of the freshwater environment in threespine sticklebacks

Colonisation of novel environments means facing new ecological challenges often resulting in the evolution of striking divergence in phenotypes. However, little is known about behavioural divergence following colonisation, despite the predicted importance of the role of behavioural phenotype-environment associations in adaptive divergence. We studied the threespine stickleback (Gasterosteus aculeatus), a model system for postglacial colonisation of freshwater habitats largely differing in ecological conditions from the ones faced by the descendants of the marine ancestor. We found that common-environment reared freshwater juveniles were less social, more active and more aggressive than their marine counterparts. This behavioural divergence could represent the result of natural selection that acted on individuals following freshwater colonisation, with predation as a key selection agent. Alternatively, the behavioural profile of freshwater juveniles could represent the characteristics of individuals that preferentially invaded freshwater after the glacial retreat, drawn from the standing variation present in the marine population.

QTL analysis of behavior in nine-spined sticklebacks (Pungitius pungitius)

The genetic architecture of behavioral traits is yet relatively poorly understood in most non-model organisms. Using an F2-intercross (n = 283 offspring) between behaviorally divergent nine-spined stickleback (Pungitius pungitius) populations, we tested for and explored the genetic basis of different behavioral traits with the aid of quantitative trait locus (QTL) analyses based on 226 microsatellite markers. The behaviors were analyzed both separately (viz. feeding activity, risk-taking and exploration) and combined in order to map composite behavioral type. Two significant QTL-explaining on average 6 % of the phenotypic variance-were detected for composite behavioral type on the experiment-wide level, located on linkage groups 3 and 8. In addition, several suggestive QTL located on six other linkage groups were detected on the chromosome-wide level. Apart from providing evidence for the genetic basis of behavioral variation, the results provide a good starting point for finer-scale analyses of genetic factors influencing behavioral variation in the nine-spined stickleback.

Ecological causes of morphological evolution in the three-spined stickleback

The central assumption of evolutionary theory is that natural selection drives the adaptation of populations to local environmental conditions, resulting in the evolution of adaptive phenotypes. The three-spined stickleback (Gasterosteus aculeatus) displays remarkable phenotypic variation, offering an unusually tractable model for understanding the ecological mechanisms underpinning adaptive evolutionary change. Using populations on North Uist, Scotland we investigated the role of predation pressure and calcium limitation on the adaptive evolution of stickleback morphology and behavior. Dissolved calcium was a significant predictor of plate and spine morph, while predator abundance was not. Stickleback latency to emerge from a refuge varied with morph, with populations with highly reduced plates and spines and high predation risk less bold. Our findings support strong directional selection in three-spined stickleback evolution, driven by multiple selective agents.

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.

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.

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.

High levels of fluctuating asymmetry in isolated stickleback populations

Background

Fluctuating 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).

Results

Pond 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 (H_E ≈ 0.3) than in marine (H_E ≈ 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.

Conclusions

The 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.