Structure and stability of genetic variance-covariance matrices: A Bayesian sparse factor analysis of transcriptional variation in the three-spined stickleback.

The genetic variance-covariance matrix (G) is a quantity of central importance in evolutionary biology due to its influence on the rate and direction of multivariate evolution. However, the predictive power of empirically estimated G-matrices is limited for two reasons. First, phenotypes are high-dimensional, whereas traditional statistical methods are tuned to estimate and analyse low-dimensional matrices. Second, the stability of G to environmental effects and over time remains poorly understood. Using Bayesian sparse factor analysis (BSFG) designed to estimate high-dimensional G-matrices, we analysed levels variation and covariation in 10,527 expressed genes in a large (n = 563) half-sib breeding design of three-spined sticklebacks subject to two temperature treatments. We found significant differences in the structure of G between the treatments: heritabilities and evolvabilities were higher in the warm than in the low-temperature treatment, suggesting more and faster opportunity to evolve in warm (stressful) conditions. Furthermore, comparison of G and its phenotypic equivalent P revealed the latter is a poor substitute of the former. Most strikingly, the results suggest that the expected impact of G on evolvability-as well as the similarity among G-matrices-may depend strongly on the number of traits included into analyses. In our results, the inclusion of only few traits in the analyses leads to underestimation in the differences between the G-matrices and their predicted impacts on evolution. While the results highlight the challenges involved in estimating G, they also illustrate that by enabling the estimation of large G-matrices, the BSFG method can improve predicted evolutionary responses to selection.

A universal and reliable assay for molecular sex identification of three-spined sticklebacks (Gasterosteus aculeatus).

In heterogametic species, biological differences between the two sexes are ubiquitous, and hence, errors in sex identification can be a significant source of noise and bias in studies where sex-related sources of variation are of interest or need to be controlled for. We developed and validated a universal multimarker assay for reliable sex identification of three-spined sticklebacks (Gasterosteus aculeatus). The assay makes use of genotype scores from three sex-linked loci and utilizes Bayesian probabilistic inference to identify sex of the genotyped individuals. The results, validated with 286 phenotypically sexed individuals from six populations of sticklebacks representing all major genetic lineages (cf. Pacific, Atlantic and Japan Sea), indicate that in contrast to commonly used single-marker-based sex identification assays, the developed multimarker assay should be 100% accurate. As the markers in the assay can be scored from agarose gels, it provides a quick and cost-efficient tool for universal sex identification of three-spined sticklebacks. The general principle of combining information from multiple markers to improve the reliability of sex identification is transferable and can be utilized to develop and validate similar assays for other species.

The genetic contribution to sex determination and number of sex chromosomes vary among populations of common frogs (Rana temporaria).

The patterns of sex determination and sex differentiation have been shown to differ among geographic populations of common frogs. Notably, the association between phenotypic sex and linkage group 2 (LG2) has been found to be perfect in a northern Swedish population, but weak and variable among families in a southern one. By analyzing these populations with markers from other linkage groups, we bring two new insights: (1) the variance in phenotypic sex not accounted for by LG2 in the southern population could not be assigned to genetic factors on other linkage groups, suggesting an epigenetic component to sex determination; (2) a second linkage group (LG7) was found to co-segregate with sex and LG2 in the northern population. Given the very short timeframe since post-glacial colonization (in the order of 1000 generations) and its seemingly localized distribution, this neo-sex chromosome system might be the youngest one described so far. It does not result from a fusion, but more likely from a reciprocal translocation between the original Y chromosome (LG2) and an autosome (LG7), causing their co-segregation during male meiosis. By generating a strict linkage between several important genes from the sex-determination cascade (Dmrt1, Amh and Amhr2), this neo-sex chromosome possibly contributes to the 'differentiated sex race' syndrome (strictly genetic sex determination and early gonadal development) that characterizes this northern population.

Evolution of anuran brains: disentangling ecological and phylogenetic sources of variation.

Variation in ecological selection pressures has been implicated to explain variation in brain size and architecture in fishes, birds and mammals, but little is known in this respect about amphibians. Likewise, the relative importance of constraint vs. mosaic hypotheses of brain evolution in explaining variation in brain size and architecture remains contentious. Using phylogenetic comparative methods, we studied interspecific variation in brain size and size of different brain parts among 43 Chinese anuran frogs and explored how much of this variation was explainable by variation in ecological factors (viz. habitat type, diet and predation risk). We also evaluated which of the two above-mentioned hypotheses best explains the observed patterns. Although variation in brain size explained on average 80.5% of the variation in size of different brain parts (supporting the constraint hypothesis), none of the three ecological factors were found to explain variation in overall brain size. However, habitat and diet type explained a significant amount of variation in telencephalon size, as well in three composite measures of brain architecture. Likewise, predation risk explained a significant amount of variation in bulbus olfactorius and optic tecta size. Our results show that evolution of anuran brain accommodates features compatible with both constraint (viz. strong allometry among brain parts) and mosaic (viz. independent size changes in response to ecological factors in certain brain parts) models of brain size evolution.

Evidence for sex-specific selection in brain: a case study of the nine-spined stickleback.

Theory predicts that the sex making greater investments into reproductive behaviours demands higher cognitive ability, and as a consequence, larger brains or brain parts. Further, the resulting sexual dimorphism can differ between populations adapted to different environments, or among individuals developing under different environmental conditions. In the nine-spine stickleback (Pungitius pungitius), males perform nest building, courtship, territory defence and parental care, whereas females perform mate choice and produce eggs. Also, predation-adapted marine and competition-adapted pond populations have diverged in a series of ecologically relevant traits, including the level of phenotypic plasticity. Here, we studied sexual dimorphism in brain size and architecture in nine-spined stickleback from marine and pond populations reared in a factorial experiment with predation and food treatments in a common garden experiment. Males had relatively larger brains, larger telencephala, cerebella and hypothalami (6-16% divergence) than females, irrespective of habitat. Females tended to have larger bulbi olfactorii than males (13%) in the high food treatment, whereas no such difference was found in the low food treatment. The strong sexual dimorphism in brain architecture implies that the different reproductive allocation strategies (behaviour vs. egg production) select for different investments into the costly brains between males and females. The lack of habitat dependence in brain sexual dimorphism suggests that the sex-specific selection forces on brains differ only negligibly between habitats. Although significance of the observed sex-specific brain plasticity in the size of bulbus olfactorius remains unclear, it demonstrates the potential for sex-specific neural plasticity.

Disentangling plastic and genetic changes in body mass of Siberian jays.

Spatial and temporal phenotypic differentiation in mean body size is of commonplace occurrence, but the underlying causes remain often unclear: both genetic differentiation in response to selection (or drift) and environmentally induced plasticity can create similar phenotypic patterns. Studying changes in body mass in Siberian jays (Perisoreus infaustus) over three decades, we discovered that mean body mass declined drastically (ca. 10%) over the first two decades, but increased markedly thereafter back to almost the initial level. Quantitative genetic analyses revealed that although body mass was heritable (h(2) = 0.46), the pronounced temporal decrease in body mass was mainly a product of phenotypic plasticity. However, a concomitant and statistically significant decrease in predicted breeding values suggests a genetic component to this change. The subsequent increase in mean body mass was indicated to be entirely due to plasticity. Selection on body mass was estimated to be too weak to fully account for the observed genetic decline in body mass, but bias in selection differential estimates due to environmental covariance between body mass and fitness is possible. Hence, the observed body mass changes appear to be driven mainly by phenotypic plasticity. Although we were not able to identify the ecological driver of the observed plastic changes, the results highlight the utility of quantitative genetic approaches in disentangling genetic and phenotypic changes in natural populations.

Mechanism of hybridization between bream Abramis brama and roach Rutilus rutilus in their native range.

Mechanisms of hybridization between bream Abramis brama and roach Rutilus rutilus were studied within the native range of the species in a lake in southern Finland. Through the genetic analysis of A. brama, R. rutilus and putative hybrids, hybridization is shown to have occurred between female A. brama and male R. rutilus. These results match with previous findings from introduced habitats, suggesting that mating between female A. brama and male R. rutilus is the predominant mechanism through which the two species hybridize.

Bringing habitat information into statistical tests of local adaptation in quantitative traits: a case study of nine-spined sticklebacks.

Detection of footprints of historical natural selection on quantitative traits in cross-sectional data sets is challenging, especially when the number of populations to be compared is small and the populations are subject to strong random genetic drift. We extend a recent Bayesian multivariate approach to differentiate between selective and neutral causes of population differentiation by the inclusion of habitat information. The extended framework allows one to test for signals of selection in two ways: by comparing the patterns of population differentiation in quantitative traits and in neutral loci, and by comparing the similarity of habitats and phenotypes. We illustrate the framework using data on variation of eight morphological and behavioral traits among four populations of nine-spined sticklebacks (Pungitius pungitius). In spite of the strong signal of genetic drift in the study system (average FST = 0.35 in neutral markers), strong footprints of adaptive population differentiation were uncovered both in morphological and behavioral traits. The results give quantitative support for earlier qualitative assessments, which have attributed the observed differentiation to adaptive divergence in response to differing ecological conditions in pond and marine habitats.

Local adaptation to salinity in the three-spined stickleback?

Different lines of evidence suggest that the occurrence and extent of local adaptation in high gene flow marine environments - even in mobile and long-lived vertebrates with complex life cycles - may be more widespread than earlier thought. We conducted a common garden experiment to test for local adaptation to salinity in Baltic Sea sticklebacks (Gasterosteus aculeatus). Fish from three different native salinity regimes (high, mid and low) were subjected to three salinity treatments (high, mid and low) in a full-factorial experimental design. Irrespective of their origin, fish subjected to low (and mid) salinity treatments exhibited higher juvenile survival, grew to largest sizes and were in better condition than fish subjected to the high salinity treatment. However, a significant interaction between native and treatment salinities - resulting mainly from the poor performance of fish native to low salinity in the high salinity treatment - provided clear cut evidence for adaptation to local variation in salinity. Additional support for this inference was provided by the fact that the results concur with an earlier demonstration of significant differentiation in a number of genes with osmoregulatory functions across the same populations and that the population-specific responses to salinity treatments exceeded that to be expected by random genetic drift.

Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks.

The evidence for adaptive phenotypic differentiation in mobile marine species remains scarce, partly due to the difficulty of obtaining quantitative genetic data to demonstrate the genetic basis of the observed phenotypic differentiation. Using a combination of phenotypic and molecular genetic approaches, we elucidated the relative roles of natural selection and genetic drift in explaining lateral plate number differentiation in threespine sticklebacks (Gasterosteus aculeatus) across the entire Baltic Sea basin (approximately 392 000 km(2) ). We found that phenotypic differentiation (PST  = 0.213) in plate number exceeded that in neutral markers (FST  = 0.008), suggesting an adaptive basis for the observed differentiation. Because a close correspondence was found between plate phenotype and genotype at a quantitative trait loci (QTL; STN381) tightly linked to the gene (Ectodysplasin) underlying plate variation, the evidence for adaptive differentiation was confirmed by comparison of FST at the QTL (FSTQ  = 0.089) with FST at neutral marker loci. Hence, the results provide a comprehensive demonstration of adaptive phenotypic differentiation in a high-gene-flow marine environment with direct, rather than inferred, verification for the genetic basis of this differentiation. In general, the results illustrate the utility of PST -FST -FSTQ comparisons in uncovering footprints of natural selection and evolution and add to the growing evidence for adaptive genetic differentiation in high-gene-flow marine environments, including that of the relatively young Baltic Sea.

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.

Evolution of stickleback feeding behaviour: genetics of population divergence at different ontogenetic stages.

The evolutionary significance of individual consistency in a given behaviour - called animal personality - has been subject to a lot of recent research. However, the genetic underpinnings of population divergence in mean personality have rarely been studied, especially across different ontogenetic stages. Previous work has shown that marine vs. pond populations of nine-spined sticklebacks (Pungitius pungitius) have undergone adaptive divergence in a series of fitness-related traits, including behaviour. One particular behavioural trait important in this system is feeding activity: giant pond sticklebacks are more active feeders than their normal sized marine conspecifics. In a common garden experiment, we raised individuals from pure and hybrid F1 -generation crosses of a highly divergent marine - pond population pair to see if (i) feeding activity and/or its ontogenetic change was consistent between individuals, and if (ii) population divergence at different ontogenetic stages could be explained by additive genetic, nonadditive genetic or maternal effects. We found that feeding activity decreased with age, but that these changes were consistently different among both individuals and crosses. The among cross patterns were consistent with a nonadditive genetic scenario: in the early period pond sticklebacks expressed dominance for high feeding activity, while in the late period marine sticklebacks expressed dominance for low feeding activity. We conclude that nine-spined sticklebacks exhibit different feeding personalities, and that the population divergence in feeding personality is explainable by age-dependent expression of genetic dominance.

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.

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.

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin.

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.

Population differences in levels of linkage disequilibrium in the wild.

Information about the levels of linkage disequilibrium (LD) in wild animal populations is still limited, and this is true particularly with respect to possible interpopulation variation in the levels of LD. We compared the levels and extent of LD at the genome-wide scale in three Siberian jay (Perisoreus infaustus) populations, two of which (Kuusamo and Ylläs) represented outbred populations within the main distribution area of the species, whereas the third (Suupohja) was a semi-isolated, partially inbred population at the margin of the species' distribution area. Although extensive long-range LD (>20 cM) was observed in all three populations, LD generally decayed to background levels at a distance of 1-5 cM or c. 200-600 kb. The degree and extent of LD differed markedly between populations but aligned closely with both observed levels of within-population genetic variation and expectations based on population history. The levels of LD were highest in the most inbred population with strong population substructure (Suupohja), compared with the two outbred populations. Furthermore, the decay of LD with increasing distance was slower in Suupohja, compared with the other two populations. By demonstrating that levels of LD can vary greatly over relatively short geographical distances within a species, these results suggest that prospects for association mapping differ from population to population. In this example, the prospects are best in the Suupohja population, given that minimized marker genotyping and a minimum marker spacing of 1-5 cM (c. 200-600 kb) would be sufficient for a whole genome scan for detecting QTL.

Intraspecific divergence in the lateral line system in the nine-spined stickleback (Pungitius pungitius).

The mechanosensory lateral line system of fishes is an important organ system conveying information crucial to individual fitness. Yet, our knowledge of lateral line diversity is almost exclusively based on interspecific studies, whereas intraspecific variability and possible population divergence have remained largely unexplored. We investigated lateral line system variability in four marine and five pond populations of nine-spined stickleback (Pungitius pungitius). We found significant differences in neuromast number between pond and marine fish. In particular, three of seventeen lateral line regions (viz. caudal peduncle superficial neuromasts; canal neuromasts from the anterior trunk and caudal peduncle) showed strong divergence between habitats. Similar results were obtained with laboratory-reared individuals from a subset of populations, suggesting that the patterns found in nature likely have a genetic basis. Interestingly, we also found habitat-dependent population divergence in neuromast variability, with pond populations showing greater heterogeneity than marine populations, although only in wild-caught fish. A comparison of neutral genetic (F(ST)) and phenotypic (P(ST)) differentiation suggested that natural selection is likely associated with habitat-dependent divergence in neuromast counts. Hence, the results align with the conclusion that the mechanosensory lateral line system divergence among marine and pond nine-spined sticklebacks is adaptive.

Brain plasticity over the metamorphic boundary: carry-over effect of larval environment on froglet brain development.

Brain development shows high plasticity in response to environmental heterogeneity. However, it is unknown how environmental variation during development may affect brain architecture across life history switch points in species with complex life cycles. Previously, we showed that predation and competition affect brain development in common frog (Rana temporaria) tadpoles. Here, we studied whether larval environment had carry-over effects in brains of metamorphs. Tadpoles grown at high density had large optic tecta at metamorphosis, whereas tadpoles grown under predation risk had small diencephala. We found that larval density had a carry-over effect on froglet optic tectum size, whereas the effect of larval predation risk had vanished by metamorphosis. We discuss the possibility that the observed changes may be adaptive, reflecting the needs of an organism in given environmental and developmental contexts.

First-generation linkage map for the common frog Rana temporaria reveals sex-linkage group.

The common frog (Rana temporaria) has become a model species in the fields of ecology and evolutionary biology. However, lack of genomic resources has been limiting utility of this species for detailed evolutionary genetic studies. Using a set of 107 informative microsatellite markers genotyped in a large full-sib family (800 F1 offspring), we created the first linkage map for this species. This partial map-distributed over 15 linkage groups-has a total length of 1698.8 cM. In line with the fact that males are the heterogametic sex in this species and a reduction of recombination is expected, we observed a lower recombination rate in the males (map length: 1371.5 cM) as compared with females (2089.8 cM). Furthermore, three loci previously documented to be sex-linked (that is, carrying male-specific alleles) in adults from the wild mapped to the same linkage group. The linkage map described in this study is one of the densest ones available for amphibians. The discovery of a sex linkage group in Rana temporaria, as well as other regions with strongly reduced male recombination rates, should help to uncover the genetic underpinnings of the sex-determination system in this species. As the number of linkage groups found (n=15) is quite close to the actual number of chromosomes (n=13), the map should provide a useful resource for further evolutionary, ecological and conservation genetic work in this and other closely related species.