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

Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

Relaxed risk of predation drives parallel evolution of stickleback behavior

The occurrence of similar phenotypes in multiple independent populations derived from common ancestral conditions (viz. parallel evolution) is a testimony of evolution by natural selection. Parallel evolution implies that populations share a common phenotypic response to a common selection pressure associated with habitat similarity. Examples of parallel evolution at genetic and phenotypic levels are fairly common, but the driving selective agents often remain elusive. Similarly, the role of phenotypic plasticity in facilitating early stages of parallel evolution is unclear. We investigated whether the relaxation of predation pressure associated with the colonization of freshwater ponds by nine-spined sticklebacks (Pungitius pungitius) likely explains the divergence in complex behaviors between marine and pond populations, and whether this divergence is parallel. Using laboratory-raised individuals exposed to different levels of perceived predation risk, we calculated vectors of phenotypic divergence for four behavioral traits between habitats and predation risk treatments. We found a significant correlation between the directions of evolutionary divergence and phenotypic plasticity, suggesting that divergence in behavior between habitats is aligned with the response to relaxation of predation pressure. Finally, we show alignment across multiple pairs of populations, and that relaxation of predation pressure has likely driven parallel evolution of behavior in this species.

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.

Patterns, Predictors, and Consequences of Dominance in Hybrids

AbstractCompared to those of their parents, are the traits of first-generation (F₁) hybrids typically intermediate, biased toward one parent, or mismatched for alternative parental phenotypes? To address this empirical gap, we compiled data from 233 crosses in which traits were measured in a common environment for two parent taxa and their F₁ hybrids. We find that individual traits in F₁s are halfway between the parental midpoint and one parental value. Considering pairs of traits together, a hybrid's bivariate phenotype tends to resemble one parent (parent bias) about 50% more than the other, while also exhibiting a similar magnitude of mismatch due to different traits having dominance in conflicting directions. Using data from an experimental field planting of recombinant hybrid sunflowers, we illustrate that parent bias improves fitness, whereas mismatch reduces fitness. Our study has three major conclusions. First, hybrids are not phenotypically intermediate but rather exhibit substantial mismatch. Second, dominance is likely determined by the idiosyncratic evolutionary trajectories of individual traits and populations. Finally, selection against hybrids likely results from selection against both intermediate and mismatched phenotypes.

Egg size and the adaptive capacity of early life history traits in Chinook salmon (Oncorhynchus tshawytscha)

Offspring traits are greatly influenced by maternal effects, and these maternal effects may provide an important pathway through which populations can adapt to changing thermal environments. We investigated the effect of egg size on the among- and within-population variation in early life history traits among introduced Great Lakes Chinook salmon (Oncorhynchus tshawytscha) populations under varying thermal conditions. We reared Chinook salmon from three populations in a common-garden hatchery study at 6.5, 9.4, and 15.2°C and measured a variety of fitness-related traits during development. We found that most of the among-population variation in early life history traits was explained by egg size. However, the contribution of egg size to the among-population variation decreased with an increase in temperature suggesting that other effects, such as genetic, contribute at high temperature. Within populations, egg size explained much of the dam variance and maternal effect for traits in every temperature, whereas egg size generally had little to no influence on the sire variance and heritability. Overall, our results demonstrate the significant contribution egg size makes to shaping early life history phenotypes among and within populations, and suggest that egg size is an important pathway through which offspring phenotypes can evolve on contemporary timescales.

Female fecundity traits in wild populations of African annual fish: the role of the aridity gradient

The evolution of life history is shaped by life expectancy. Life-history traits coevolve, and optimal states for particular traits are constrained by trade-offs with other life-history traits. Life histories contrast among species, but may also diverge intraspecifically, at the level of populations. We studied the evolution of female reproductive allocation strategy, using natural populations of two sympatric species of African annual fishes, Nothobranchius furzeri and Nothobranchius orthonotus. These species inhabit pools in the Mozambican savanna that are formed in the rainy season and persist for only 2-10 months. Using 207 female N. furzeri from 11 populations and 243 female N. orthonotus from 14 populations, we tested the effects of genetic background (intraspecific lineage) and life expectancy (position on the aridity gradient determining maximum duration of their temporary habitat) on female fecundity traits. First, we found that variation in female body mass was small within populations, but varied considerably among populations. Second, we found that fecundity was largely defined by female body mass and that females spawned most of their eggs in the morning. Third, we found that the trade-off between egg size and egg number varied among lineages of N. furzeri and this outcome has been confirmed by data from two separate years. Overall, we demonstrate that local conditions were important determinants for Nothobranchius growth and fecundity and that eggs size in arid region was less limited by female fecundity than in humid region.

Quantitative trait locus analysis of body shape divergence in nine-spined sticklebacks based on high-density SNP-panel

Heritable phenotypic differences between populations, caused by the selective effects of distinct environmental conditions, are of commonplace occurrence in nature. However, the actual genomic targets of this kind of selection are still poorly understood. We conducted a quantitative trait locus (QTL) mapping study to identify genomic regions responsible for morphometric differentiation between genetically and phenotypically divergent marine and freshwater nine-spined stickleback (Pungitius pungitius) populations. Using a dense panel of SNP-markers obtained by restriction site associated DNA sequencing of an F2 recombinant cross, we found 22 QTL that explained 3.5-12.9% of phenotypic variance in the traits under investigation. We detected one fairly large-effect (PVE = 9.6%) QTL for caudal peduncle length-a trait with a well-established adaptive function showing clear differentiation among marine and freshwater populations. We also identified two large-effect QTL for lateral plate numbers, which are different from the lateral plate QTL reported in earlier studies of this and related species. Hence, apart from identifying several large-effect QTL in shape traits showing adaptive differentiation in response to different environmental conditions, the results suggest intra- and interspecific heterogeneity in the genomic basis of lateral plate number variation.

Genetic basis for body size variation between an anadromous and two derived lacustrine populations of threespine stickleback Gasterosteus aculeatus in southwest Alaska

Body size is a highly variable trait among geographically separated populations. Size-assortative reproductive isolation has been linked to recent adaptive radiations of threespine stickleback (Gasterosteus aculeatus) into freshwater, but the genetic basis of the commonly found size difference between anadromous and derived lacustrine sticklebacks has not been tested. We studied the genetic basis of size differences between recently diverging stickleback lineages in southwest Alaska using a common environment experiment. We crossed stickleback within one anadromous (Naknek River) and one lake (Pringle Lake) population and between the anadromous and two lake populations (Pringle and JoJo Lakes), and raised them in a salinity of 4–6 ppt. The F1 anadromous and freshwater forms differed significantly in size, whereas hybrids were intermediate or exhibited dominance toward the anadromous form. Additionally, the size of freshwater F1s differed from their wild counterparts, with within-population F1s from Pringle Lake growing larger than their wild counterparts, while there was no size difference between lab-raised and wild anadromous fish. Sexual dimorphism was always present in anadromous fish, but not in freshwater, and not always in the hybrid crosses. These results, along with parallel changes among anadromous and freshwater forms in other regions, suggest that this heritable trait is both plastic and may be under divergent and/or sexual selection.

Effects of Large-Scale Releases on the Genetic Structure of Red Sea Bream (Pagrus major, Temminck et Schlegel) Populations in Japan

Large-scale hatchery releases are carried out for many marine fish species worldwide; nevertheless, the long-term effects of this practice on the genetic structure of natural populations remains unclear. The lack of knowledge is especially evident when independent stock enhancement programs are conducted simultaneously on the same species at different geographical locations, as occurs with red sea bream (Pagrus major, Temminck et Schlegel) in Japan. In this study, we examined the putative effects of intensive offspring releases on the genetic structure of red sea bream populations along the Japanese archipelago by genotyping 848 fish at fifteen microsatellite loci. Our results suggests weak but consistent patterns of genetic divergence (F_ST = 0.002, p < 0.001). Red sea bream in Japan appeared spatially structured with several patches of distinct allelic composition, which corresponded to areas receiving an important influx of fish of hatchery origin, either released intentionally or from unintentional escapees from aquaculture operations. In addition to impacts upon local populations inhabiting semi-enclosed embayments, large-scale releases (either intentionally or from unintentional escapes) appeared also to have perturbed genetic structure in open areas. Hence, results of the present study suggest that independent large-scale marine stock enhancement programs conducted simultaneously on one species at different geographical locations may compromise native genetic structure and lead to patchy patterns in population genetic structure.

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.

The association between parental life history and offspring phenotype

In many taxa there is considerable intraspecific variation in life-history strategies from within a single population, reflecting alternative routes through which organisms can achieve successful reproduction. Atlantic salmon Salmo salar (Linnaeus) show some of the greatest within-population variability in life history strategies amongst vertebrates, with multiple discrete male and female life histories co-existing and interbreeding on many spawning grounds, although the effect of the various combinations of life histories on offspring traits remains unknown. Using crosses of wild fish we show here that the life history strategy of both parents was significantly associated with a range of offspring traits. Mothers that had spent longer at sea (two versus one year) produced offspring which were heavier, longer and in better condition at the time of first feeding. However, these relationships disappeared shortly after fry had begun feeding exogenously. At this stage, the juvenile rearing environment (i.e. time spent in fresh water as juveniles) of the mother was a better predictor of offspring traits, with mothers that were faster to develop in fresh water (migrating to sea after two rather than three years of age) producing offspring that had higher maximal metabolic rates, aerobic scopes, and that grew faster. Faster developing fathers (1 year old sneaker males) tended to produce offspring that had higher maximal metabolic rates, were in better body condition and grew faster. The results suggest that both genetic effects and those related to parental early and late life history contribute to offspring traits.

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.

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.

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.