Comparative studies of the divergence of quantitative traits and neutral molecular markers, known as Q(ST)-F(ST) comparisons, provide a means for researchers to distinguish between natural selection and genetic drift as causes of population differentiation in complex polygenic traits. The use of Q(ST)-F(ST) comparisons has increased rapidly in the last few years, highlighting the utility of this approach for addressing a wide range of questions that are relevant to evolutionary and ecological genetics. These studies have also provided lessons for the design of future Q(ST)-F(ST) comparisons. Methods based on the Q(ST)-F(ST) approach could also be used to analyse various types of 'omics' data in new and revealing ways.

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 (Q(ST)) and neutral marker (F(ST)) differentiation among natural populations. Our analyses confirm the conclusion from previous reviews - based on ca. 100% more data - that the Q(ST) values are on average higher than F(ST) 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 Q(ST) and F(ST) 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 Q(ST) and F(ST) 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 Q(ST)-F(ST) comparisons as a tool for inferring natural selection, as well as associated methodological and interpretational problems involved with individual and meta-analytic studies.

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.

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.

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.

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.

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.

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.