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

Complex and Dynamic Gene-by-Age and Gene-by-Environment Interactions Underlie Functional Morphological Variation in Adaptive Divergence in Arctic Charr (Salvelinus alpinus)

The evolution of adaptive phenotypic divergence requires heritable genetic variation. However, it is underappreciated that trait heritability is molded by developmental processes interacting with the environment. We hypothesized that the genetic architecture of divergent functional traits was dependent on age and foraging environment. Thus, we induced plasticity in full-sib families of Arctic charr (Salvelinus alpinus) morphs from two Icelandic lakes by mimicking prey variation in the wild. We characterized variation in body shape and size at two ages and investigated their genetic architecture with quantitative trait locus (QTL) analysis. Age had a greater effect on body shape than diet in most families, suggesting that development strongly influences phenotypic variation available for selection. Consistent with our hypothesis, multiple QTL were detected for all traits and their location depended on age and diet. Many of the genome-wide QTL were located within a subset of duplicated chromosomal regions suggesting that ancestral whole genome duplication events have played a role in the genetic control of functional morphological variation in the species. Moreover, the detection of two body shape QTL after controlling for the effects of age provides additional evidence for genetic variation in the plastic response of morphological traits to environmental variation. Thus, functional morphological traits involved in phenotypic divergence are molded by complex genetic interactions with development and environment.

Heterogeneous genomic architecture of skeletal armour traits in sticklebacks

Whether populations adapt to similar selection pressures using the same underlying genetic variants depends on population history and the distribution of standing genetic variation at the metapopulation level. Studies of sticklebacks provide a case in point: when colonizing and adapting to freshwater habitats, three-spined sticklebacks (Gasterosteus aculeatus) with high gene flow tend to fix the same adaptive alleles in the same major loci, whereas nine-spined sticklebacks (Pungitius pungitius) with limited gene flow tend to utilize a more heterogeneous set of loci. In accordance with this, we report results of quantitative trait locus (QTL) analyses using a backcross design showing that lateral plate number variation in the western European nine-spined sticklebacks mapped to 3 moderate-effect QTL, contrary to the major-effect QTL in three-spined sticklebacks and different from the 4 QTL previously identified in the eastern European nine-spined sticklebacks. Furthermore, several QTL were identified associated with variation in lateral plate size, and 3 moderate-effect QTL with body size. Together, these findings indicate more heterogenous and polygenic genetic underpinnings of skeletal armour variation in nine-spined than three-spined sticklebacks, indicating limited genetic parallelism underlying armour trait evolution in the family Gasterostidae.

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.