Hannu Mäkinen

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.

Identifying footprints of directional and balancing selection in marine and freshwater three‐spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within‐population genetic diversity and increase among‐population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three‐spined stickleback (Gasterosteus aculeatus) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)‐associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda‐gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL‐mapping studies revealed no signatures of selection. Two novel microsatellite loci (Stn12 and Stn90) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.

Genetic relationships among marine and freshwater populations of the European three-spined stickleback (Gasterosteus aculeatus) revealed by microsatellites

To assess the population genetic structure of the three‐spined stickleback, Gasterosteus aculeatus, variability at 18 microsatellite loci was examined in 1724 individuals from 74 locations covering most of the species distribution range in Europe. The results revealed high overall degree of differentiation (FST = 0.21) but contrasting level of divergence and genetic variability between habitat types. Marine populations were genetically relatively uniform even across great geographical distances as compared to substantial differentiation among freshwater populations. Analysis of molecular variance indicated low but significant (2.7%) variation in allele frequencies between geographical regions, but a negligible effect of habitat type (0.2%). The phylogenetic pattern was not explained by habitat type, but a weak signal of populations clustering according to geographical or water system origin was found. The results support the view that three‐spined stickleback marine ancestors colonized northern European fresh waters during the postglacial marine submergence c. 10 000 years ago, whereas in the Mediterranean region colonization probably dates back to the Pleistocene. The independent origins of river and lake populations indicate that they originate from multiple colonizations rather than sharing common ancestry. In the continuous marine environment, the low degree of differentiation among populations can be explained by gene flow among subpopulations and large effective population size buffering divergence in neutral markers. In contrast, among postglacially established freshwater populations differentiation appears to be driven by genetic drift and isolation. The stepwise mutations appear to have contributed to the population differentiation in the southern part of the three‐spined stickleback distribution range.

Genetic evidence for male‐biased dispersal in the three‐spined stickleback (Gasterosteus aculeatus)

Sex‐biased dispersal is capable of generating population structure in nonisolated populations and may affect adaptation processes when selective conditions differ among populations. Intrasexual competition for local resources and/or mating opportunities predicts a male‐biased dispersal in polygynous species and a female bias in monogamous species. The patterns of sex‐biased dispersal in birds and mammals are well explained by their respective mating systems, but the picture emerging from fish studies is still mixed. Using neutral genetic markers, we investigated whether there is any evidence for sex‐biased dispersal among Baltic Sea populations of the three‐spined stickleback (Gasterosteus aculeatus). The null hypothesis of non sex‐biased dispersal was rejected in favour of male‐biased dispersal in this species. As the three‐spined stickleback has a polygynous mating system, the observed male bias in dispersal is consistent with the hypothesis that local mate competition might drive the observed pattern. Although more research both on the proximate and ultimate causes behind the observed pattern is needed, our results serve as a first step towards understanding patterns of sex‐biased dispersal in this species.

Hitchhiking Mapping Reveals a Candidate Genomic Region for Natural Selection in Three-Spined Stickleback Chromosome VIII

Identification of genes and genomic regions under directional natural selection has become one of the major goals in evolutionary genetics, but relatively little work to this end has been done by applying hitchhiking mapping to wild populations. Hitchhiking mapping starts from a genome scan using a randomly spaced set of molecular markers followed by a fine-scale analysis in the flanking regions of the candidate regions under selection. We used the hitchhiking mapping approach to narrow down a selective sweep in the genomic region flanking a candidate locus (Stn90) in chromosome VIII in the three-spined stickleback (Gasterosteus aculeatus). Twenty-four microsatellite markers were screened in an ∼800-kb region around the candidate locus in three marine and four freshwater populations. The patterns of genetic diversity and differentiation in the candidate region were compared to those of a putatively neutral set of markers. The Bayesian FST-test indicated an elevated genetic differentiation, deviating significantly from neutral expectations, at a continuous region of ∼20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of ∼90 kb flanking the candidate locus supporting the selective sweep hypothesis. Likewise, the genomewide pattern of genetic diversity differed from the candidate region in a bottleneck analysis suggesting that selection, rather than demography, explains the reduced genetic diversity at the candidate interval. The neutrality tests suggest that the selective sweep had occurred mainly in the Lake Pulmanki population, but the results from bottleneck analyses indicate that selection might have operated in other populations as well. These results suggest that the narrow interval around locus Stn90 has likely been under directional selection, but the region contains several predicted genes, each of which can be the actual targets of selection. Understanding of the functional significance of this genomic region in an ecological context will require a more detailed sequence analysis.

The utility of QTL-linked markers to detect selective sweeps in natural populations — a case study of the EDA gene and a linked marker in threespine stickleback

Sequence polymorphisms in coding genes and variability in quantitative trait loci (QTL)‐linked markers can be used to uncover the evolutionary mechanisms of traits involved in adaptive processes. We studied sequence variation in the EDA gene and allelic variation in 18 microsatellites — one of which (Gac4174) is linked with the EDA QTL — in low, partially and completely plated morphs from eight threespine stickleback European populations. The results agree with previous studies in that EDA polymorphism is closely related to plate number variation: EDA sequences grouped populations into low and completely plated morphs, whereas microsatellites failed to do so. Furthermore, partially plated fish were heterozygous with respect to the distinctive EDA alleles for completely and low plated morphs, indicating that completely plated morph alleles are not entirely dominant in controlling the expression of lateral plate number. An examination of population differentiation in plate number with quantitative genetic methods revealed that the degree of differentiation exceeded that expected from genetic drift alone (QST > FST). Our results support the adaptive genetic differentiation of plate morphs and the view that distinctive EDA gene polymorphism occurs in similar sites across the distribution range of this species. Yet, allele frequency differentiation in the Gac4174 microsatellite locus, informative in experimental crosses for plate number variation, did not differ from that of neutral markers and, was therefore unable to detect the signature of natural selection responsible for population divergence.

Population genomic analyses of early‐phase Atlantic Salmon (Salmo salar) domestication/captive breeding

Domestication can have adverse genetic consequences, which may reduce the fitness of individuals once released back into the wild. Many wild Atlantic salmon (Salmo salar L.) populations are threatened by anthropogenic influences, and they are supplemented with captively bred fish. The Atlantic salmon is also widely used in selective breeding programs to increase the mean trait values for desired phenotypic traits. We analyzed a genomewide set of SNPs in three domesticated Atlantic salmon strains and their wild conspecifics to identify loci underlying domestication. The genetic differentiation between domesticated strains and wild populations was low (FST < 0.03), and domesticated strains harbored similar levels of genetic diversity compared to their wild conspecifics. Only a few loci showed footprints of selection, and these loci were located in different linkage groups among the different wild population/hatchery strain comparisons. Simulated scenarios indicated that differentiation in quantitative trait loci exceeded that in neutral markers during the early phases of divergence only when the difference in the phenotypic optimum between populations was large. This study indicates that detecting selection using standard approaches in the early phases of domestication might be challenging unless selection is strong and the traits under selection show simple inheritance patterns.

Plastic and Evolutionary Gene Expression Responses Are Correlated in European Grayling (Thymallus thymallus) Subpopulations Adapted to Different Thermal Environments

Understanding how populations adapt to changing environmental conditions is a long-standing theme in evolutionary biology. Gene expression changes have been recognized as an important driver of local adaptation, but relatively little is known regarding the direction of change and in particular, about the interplay between plastic and evolutionary gene expression. We have previously shown that the gene expression profiles of European grayling (Thymallus thymallus) populations inhabiting different thermal environments include both plastic and evolutionary components. However, whether the plastic and evolutionary responses were in the same direction was not investigated in detail, nor was the identity of the specific genes involved. In this study, we show that the plastic changes in protein expression in response to different temperatures are highly correlated with the evolutionary response in grayling subpopulations adapted to different thermal environments. This finding provides preliminary evidence that the plastic response most likely facilitates adaptation during the early phases of colonization of thermal environments. The proteins that showed significant changes in expression level between warm and cold temperature treatments were mostly related to muscle development, which is consistent with earlier findings demonstrating muscle mass differentiation between cold and warm grayling populations.

Prediction of harmful variants on mitochondrial genes: Test of habitat‐dependent and demographic effects in a euryhaline fish

Abstract Both effective population size and life history may influence the efficacy of purifying selection, but it remains unclear if the environment affects the accumulation of weakly deleterious nonsynonymous polymorphisms. We hypothesize that the reduced energetic cost of osmoregulation in brackish water habitat may cause relaxation of selective constraints at mitochondrial oxidative phosphorylation (OXPHOS) genes. To test this hypothesis, we analyzed 57 complete mitochondrial genomes of Pungitius pungitius collected from brackish and freshwater habitats. Based on inter‐ and intraspecific comparisons, we estimated that 84% and 68% of the nonsynonymous polymorphisms in the freshwater and brackish water populations, respectively, are weakly or moderately deleterious. Using in silico prediction tools (MutPred, SNAP2), we subsequently identified nonsynonymous polymorphisms with potentially harmful effect. Both prediction methods indicated that the functional effects of the fixed nonsynonymous substitutions between nine‐ and three‐spined stickleback were weaker than for polymorphisms within species, indicating that harmful nonsynonymous polymorphisms within populations rarely become fixed between species. No significant differences in mean estimated functional effects were identified between freshwater and brackish water nine‐spined stickleback to support the hypothesis that reduced osmoregulatory energy demand in the brackish water environment reduces the strength of purifying selection at OXPHOS genes. Instead, elevated frequency of nonsynonymous polymorphisms in the freshwater environment (P n /P s = 0.549 vs. 0.283; Fishers exact test p = .032) suggested that purifying selection is less efficient in small freshwater populations. This study shows the utility of in silico functional prediction tools in population genetic and evolutionary research in a nonmammalian vertebrate and demonstrates that mitochondrial energy production genes represent a promising system to characterize the demographic, life history and potential habitat‐dependent effects of segregating amino acid variants.

Displaced phylogeographic signals from Gyrodactylus arcuatus, a parasite of the three-spined stickleback Gasterosteus aculeatus, suggest freshwater glacial refugia in Europe

We examined the global mitochondrial phylogeography of Gyrodactylus arcuatus, a flatworm ectoparasite of three-spined stickleback Gasterosteus aculeatus. In accordance with the suggested high divergence rate of 13%/million years, the genetic variation of the parasite was high: haplotype diversity h=0.985 and nucleotide diversity π=0.0161. The differentiation among the parasite populations was substantial (Φst=0.759), with two main allopatric clades (here termed Euro and North) accounting for 54% of the total genetic variation. The diversity center of the Euro clade was in the Baltic Sea, while the North clade was spread across the Barents and White Seas. A single haplotype within the North clade was found in the western and eastern Pacific Ocean. Divergence of main clades was estimated to be circa 200 thousand years ago. Each main clade was further divided into six distinct subclades, estimated to have diverged in isolation since 135 thousand years ago. This second division corresponds approximately to the Eemian interglacial predating the last glacial maximum. A demographic expansion of the subclades is associated with colonisation of northern Europe since the last glacial maximum, circa 15-40 thousand years ago. The parasite phylogeny is most likely explained by sequential isolated bottlenecks and expansions in numerous allopatric refugia. The postglacial intermingling and high variation in the marine parasite populations, separately in the Baltic and Barents Seas, suggest low competition of divergent parasite matrilines, coupled with a large population size and high rate of dispersal of hosts. The genetic contribution of the assumed refugial fish populations maintaining the parasite during the last glacial maximum was not detected among the marine sticklebacks, which perhaps were infected after range expansion.