Tuomas Leinonen

Comparative studies of quantitative trait and neutral marker divergence: a meta-analysis

Comparative studies of quantitative genetic and neutral marker differentiation have provided means for assessing the relative roles of natural selection and random genetic drift in explaining among‐population divergence. This information can be useful for our fundamental understanding of population differentiation, as well as for identifying management units in conservation biology. Here, we provide comprehensive review and meta‐analysis of the empirical studies that have compared quantitative genetic (QST) and neutral marker (FST) differentiation among natural populations. Our analyses confirm the conclusion from previous reviews – based on ca. 100% more data – that the QST values are on average higher than FST values [mean difference 0.12 (SD 0.27)] suggesting a predominant role for natural selection as a cause of differentiation in quantitative traits. However, although the influence of trait (life history, morphological and behavioural) and marker type (e.g. microsatellites and allozymes) on the variance of the difference between QST and FST is small, there is much heterogeneity in the data attributable to variation between specific studies and traits. The latter is understandable as there is no reason to expect that natural selection would be acting in similar fashion on all populations and traits (except for fitness itself). We also found evidence to suggest that QST and FST values across studies are positively correlated, but the significance of this finding remains unclear. We discuss these results in the context of utility of the QST–FST comparisons as a tool for inferring natural selection, as well as associated methodological and interpretational problems involved with individual and meta‐analytic studies.

Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.

Female-Biased Expression on the X Chromosome as a Key Step in Sex Chromosome Evolution in Threespine Sticklebacks

Given that the genome of males and females are almost identical with the exception of genes on the Y (or W) chromosome or sex-determining alleles (in organisms without sex chromosomes), it is likely that many downstream processes resulting in sexual dimorphism are produced by changes in regulation. In early stages of sex chromosome evolution, as the Y-chromosome degenerates, gene expression should be significantly impacted for genes residing on the sex chromosome pair as regulatory mutations accumulate. However, this has rarely been examined because most model organisms have clearly diverged sex chromosomes. Fish provide a unique opportunity to examine the evolution of sex chromosomes because genetic sex determination has evolved quite recently in some groups of fish. We compared sex-specific transcription in threespine stickleback (Gasterosteus aculeatus) liver tissue using a long-oligo microarray. Of the 1,268 genes that were differentially expressed between sexes, a highly significant proportion (23%) was concentrated on chromosome 19, corresponding to the recently described nascent sex chromosomes. The sex-biased genes are enriched for different functional categories in males and females, although there is no specific functional enrichment on the sex chromosomes. Female-biased genes are concentrated at one end of the sex chromosome, corresponding to a deletion in the Y, suggesting a lack of global dosage compensation. Prior research on threespine sticklebacks has demonstrated various degrees of dissimilarity in upstream regions of genes on the Y providing a potential mechanism for the observed patterns of female-biased expression. We hypothesize that degeneration of the Y chromosome results in regulatory mutations that create a sex-specific expression pattern and that this physical concentration of sex-biased expression on the nascent sex chromosome may be a key feature characterizing intermediate phases of sex chromosome evolution.

Genetics of body shape and armour variation in threespine sticklebacks

Patterns of genetic variation and covariation can influence the rate and direction of phenotypic evolution. We explored the possibility that the parallel morphological evolution seen in threespine stickleback (Gasterosteus aculeatus) populations colonizing freshwater environments is facilitated by patterns of genetic variation and covariation in the ancestral (marine) population. We estimated the genetic (G) and phenotypic (P) covariance matrices and directions of maximum additive genetic (gmax) and phenotypic (pmax) covariances of body shape and armour traits. Our results suggest a role for the ancestral G in explaining parallel morphological evolution in freshwater populations. We also found evidence of genetic constraints owing to the lack of variance in the ancestral G. Furthermore, strong genetic covariances and correlations among traits revealed that selective factors responsible for threespine stickleback body shape and armour divergence may be difficult to disentangle. The directions of gmax and pmax were correlated, but the correlations were not high enough to imply that phenotypic patterns of trait variation and covariation within populations are very informative of underlying genetic patterns.

The Evolution and Adaptive Potential of Transcriptional Variation in Sticklebacks—Signatures of Selection and Widespread Heritability

Evidence implicating differential gene expression as a significant driver of evolutionary novelty continues to accumulate, but our understanding of the underlying sources of variation in expression, both environmental and genetic, is wanting. Heritability in particular may be underestimated when inferred from genetic mapping studies, the predominant “genetical genomics” approach to the study of expression variation. Such uncertainty represents a fundamental limitation to testing for adaptive evolution at the transcriptomic level. By studying the inheritance of expression levels in 10,495 genes (10,527 splice variants) in a threespine stickleback pedigree consisting of 563 individuals, half of which were subjected to a thermal treatment, we show that 74–98% of transcripts exhibit significant additive genetic variance. Dominance variance is also prevalent (41–99% of transcripts), and genetic sources of variation seem to play a more significant role in expression variance in the liver than a key environmental variable, temperature. Among-population comparisons suggest that the majority of differential expression in the liver is likely due to neutral divergence; however, we also show that signatures of directional selection may be more prevalent than those of stabilizing selection. This predominantly aligns with the neutral model of evolution for gene expression but also suggests that natural selection may still act on transcriptional variation in the wild. As genetic variation both within- and among-populations ultimately defines adaptive potential, these results indicate that broad adaptive potential may be found within the transcriptome.

Predation-imposed selection on threespine stickleback (Gasterosteus aculeatus) morphology: a test of the refuge use hypothesis.

The transition from marine to freshwater life in the threespine stickleback (Gasterosteus aculeatus) is accompanied by complex morphological changes—including reduction in bony armor and change in body shape—but experimental evidence for the selective agents behind these evolutionary transitions is sparse. We investigated whether selection by predatory fish affects threespine stickleback morphology differentially when refuge is absent (pelagic lifestyle—ancestral condition) or present (benthic lifestyle—derived condition). Our results show that selection favors low numbers of lateral plates in habitats with refuge, whereas fully plated individuals have a selective advantage in habitats without refuge. We also found that a decrease in the length of the caudal peduncle increased survival probability, irrespective of habitat. The effect of spine lengths on survival was evident only in a multivariate analysis of selection, implying that it is essential to account for phenotypic and genetic correlations between traits before drawing conclusions about the effects of selection on single traits. Apart from uncovering targets and patterns of predator‐induced selection on threespine stickleback morphology, our results provide direct evidence to support the hypothesis that differences in antipredator strategies in pelagic versus benthic sticklebacks could play a role in the repeated, independent cases of plate number reduction following freshwater colonization in this species.

Genetic basis of sexual dimorphism in the threespine stickleback Gasterosteus aculeatus

Sexual dimorphism (SD) in morphological, behavioural and physiological features is common, but the genetics of SD in the wild has seldom been studied in detail. We investigated the genetic basis of SD in morphological traits of threespine stickleback (Gasterosteus aculeatus) by conducting a large breeding experiment with fish from an ancestral marine population that acts as a source of morphological variation. We also examined the patterns of SD in a set of 38 wild populations from different habitats to investigate the relationship between the genetic architecture of SD of the marine ancestral population in relation to variation within and among natural populations. The results show that genetic architecture in terms of heritabilities, additive genetic variances and covariances (as well as correlations) is very similar in the two sexes in spite of the fact that many of the traits express significant SD. Furthermore, population differences in threespine stickleback body shape and armour SD appear to have evolved despite constraints imposed by genetic architecture. This implies that constraints for the evolution of SD imposed by strong genetic correlations are not as severe and absolute as commonly thought.

Quantitative Genetics of Body Size and Timing of Maturation in Two Nine-Spined Stickleback (Pungitius pungitius) Populations

Due to its influence on body size, timing of maturation is an important life-history trait in ectotherms with indeterminate growth. Comparison of patterns of growth and maturation within and between two populations (giant vs. normal sized) of nine-spined sticklebacks (Pungitius pungitius) in a breeding experiment revealed that the difference in mean adult body size between the populations is caused by differences in timing of maturation, and not by differential growth rates. The fish in small-sized population matured earlier than those from large-sized population, and maturation was accompanied by a reduction in growth rate in the small-sized population. Males matured earlier and at smaller size than females, and the fish that were immature at the end of the experiment were larger than those that had already matured. Throughout the experimental period, body size in both populations was heritable (h2 = 0.10–0.64), as was the timing of maturation in the small-sized population (h2 = 0.13–0.16). There was a significant positive genetic correlation between body size and timing of maturation at 140 DAH, but not earlier (at 80 or 110 DAH). Comparison of observed body size divergence between the populations revealed that Q ST exceeded F ST at older ages, indicating adaptive basis for the observed divergence. Hence, the results suggest that the body size differences within and between populations reflect heritable genetic differences in the timing of maturation, and that the observed body size divergence is adaptive.

driftsel: an R package for detecting signals of natural selection in quantitative traits

Approaches and tools to differentiate between natural selection and genetic drift as causes of population differentiation are of frequent demand in evolutionary biology. Based on the approach of Ovaskainen et al. (2011), we have developed an R package (driftsel) that can be used to differentiate between stabilizing selection, diversifying selection and random genetic drift as causes of population differentiation in quantitative traits when neutral marker and quantitative genetic data are available. Apart from illustrating the use of this method and the interpretation of results using simulated data, we apply the package on data from three‐spined sticklebacks (Gasterosteus aculeatus) to highlight its virtues. driftsel can also be used to perform usual quantitative genetic analyses in common‐garden study designs.

Transcription and redox enzyme activities: comparison of equilibrium and disequilibrium levels in the three-spined stickleback

Evolutionary and acclimatory responses require functional variability, but in contrast with mRNA and protein abundance data, most physiological measurements cannot be obtained in a high-throughput manner. Consequently, one must either rely on high-throughput transcriptomic or proteomic data with only predicted functional information, or accept the limitation that most physiological measurements can give fewer data than those provided by transcriptomics or proteomics. We evaluated how transcriptional and redox enzyme activity data agreed with regard to population differentiation (i.e. a system in steady state in which any time lag between transcription, translation and post-translational effects would be irrelevant) and in response to an acute 6°C increase in temperature (i.e. a disequilibrium state wherein translation could not have caught up with transcription) in the three-spined stickleback (Gasterosteus aculeatus). Transcriptional and enzyme activity data corresponded well with regard to population differentiation, but less so with regard to acute temperature increase. The data thus suggest that transcriptional and functional measurements can lead to similar conclusions when a biological system is in a steady state. The responses to acute changes must, as has been demonstrated earlier, be based on changes in cellular conditions or properties of existing proteins without significant de novo synthesis of new gene products.