Genetic population structure constrains local adaptation in sticklebacks.

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

A phylogenomic perspective on diversity, hybridization and evolutionary affinities in the stickleback genus Pungitius.

Hybridization and convergent evolution are phenomena of broad interest in evolutionary biology, but their occurrence poses challenges for reconstructing evolutionary affinities among affected taxa. Sticklebacks in the genus Pungitius are a case in point: evolutionary relationships and taxonomic validity of different species and populations in this circumpolarly distributed species complex remain contentious due to convergent evolution of traits regarded as diagnostic in their taxonomy, and possibly also due to frequent hybridization among taxa. To clarify the evolutionary relationships among different Pungitius species and populations globally, as well as to study the prevalence and extent of introgression among recognized species, genomic data sets of both reference genome-anchored single nucleotide polymorphisms and de novo assembled RAD-tag loci were constructed with RAD-seq data. Both data sets yielded topologically identical and well-supported species trees. Incongruence between nuclear and mitochondrial DNA-based trees was found and suggested possibly frequent hybridization and mitogenome capture during the evolution of Pungitius sticklebacks. Further analyses revealed evidence for frequent nuclear genetic introgression among Pungitius species, although the estimated proportions of autosomal introgression were low. Apart from providing evidence for frequent hybridization, the results challenge earlier mitochondrial and morphology-based hypotheses regarding the number of species and their affinities in this genus: at least seven extant species can be recognized on the basis of genetic data. The results also shed new light on the biogeographical history of the Pungitius-complex, including suggestion of several trans-Arctic invasions of Europe from the Northern Pacific. The well-resolved phylogeny should facilitate the utility of this genus as a model system for future comparative evolutionary studies.

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene-mapping approach.

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small-effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine-spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL-mapping approaches based on a de-biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL-mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single-locus analyses of an F2 interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1-6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small-effect loci. Apart from suggesting moderately heritable (h2  ≈ 0.15-0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL-mapping approach developed here has distinctive advantages over the traditional QTL-mapping methods in analyses of dense marker panels.

Progressive recombination suppression and differentiation in recently evolved neo-sex chromosomes.

Recombination suppression leads to the structural and functional differentiation of sex chromosomes and is thus a crucial step in the process of sex chromosome evolution. Despite extensive theoretical work, the exact processes and mechanisms of recombination suppression and differentiation are not well understood. In threespine sticklebacks (Gasterosteus aculeatus), a different sex chromosome system has recently evolved by a fusion between the Y chromosome and an autosome in the Japan Sea lineage, which diverged from the ancestor of other lineages approximately 2 Ma. We investigated the evolutionary dynamics and differentiation processes of sex chromosomes based on comparative analyses of these divergent lineages using 63 microsatellite loci. Both chromosome-wide differentiation patterns and phylogenetic inferences with X and Y alleles indicated that the ancestral sex chromosomes were extensively differentiated before the divergence of these lineages. In contrast, genetic differentiation appeared to have proceeded only in a small region of the neo-sex chromosomes. The recombination maps constructed for the Japan Sea lineage indicated that recombination has been suppressed or reduced over a large region spanning the ancestral and neo-sex chromosomes. Chromosomal regions exhibiting genetic differentiation and suppressed or reduced recombination were detected continuously and sequentially in the neo-sex chromosomes, suggesting that differentiation has gradually spread from the fusion point following the extension of recombination suppression. Our study illustrates an ongoing process of sex chromosome differentiation, providing empirical support for the theoretical model postulating that recombination suppression and differentiation proceed in a gradual manner in the very early stage of sex chromosome evolution.

QTL analysis of behavior in nine-spined sticklebacks (Pungitius pungitius).

The genetic architecture of behavioral traits is yet relatively poorly understood in most non-model organisms. Using an F2-intercross (n = 283 offspring) between behaviorally divergent nine-spined stickleback (Pungitius pungitius) populations, we tested for and explored the genetic basis of different behavioral traits with the aid of quantitative trait locus (QTL) analyses based on 226 microsatellite markers. The behaviors were analyzed both separately (viz. feeding activity, risk-taking and exploration) and combined in order to map composite behavioral type. Two significant QTL-explaining on average 6 % of the phenotypic variance-were detected for composite behavioral type on the experiment-wide level, located on linkage groups 3 and 8. In addition, several suggestive QTL located on six other linkage groups were detected on the chromosome-wide level. Apart from providing evidence for the genetic basis of behavioral variation, the results provide a good starting point for finer-scale analyses of genetic factors influencing behavioral variation in the nine-spined stickleback.

High degree of genetic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus).

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene-based genome-scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection-as determined with several outlier detection methods-was low (F(ST) = 0.033 or 0.011, respectively), whereas average FST for directionally selected markers was significantly higher (F(ST) = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.

Quantitative trait loci for growth and body size in the nine-spined stickleback Pungitius pungitius L.

Body size is an ecologically important trait shown to be genetically variable both within and among different animal populations as revealed by quantitative genetic studies. However, few studies have looked into underlying genetic architecture of body size variability in the wild using genetic mapping methods. With the aid of quantitative trait loci (QTL) analyses based on 226 microsatellite markers, we mapped body size and growth rate traits in the nine-spined stickleback (Pungitius pungitius) using an F2 -intercross (n = 283 offspring) between size-divergent populations. In total, 17 QTL locations were detected. The proportion of phenotypic variation explained by individual body size-related QTL ranged from 3% to 12% and those related to growth parameters and increments from 3% to 10%. Several of the detected QTL affected either early or late growth. These results provide a solid starting point for more in depth investigations of structure and function of genomic regions involved in determination of body size in this popular model of ecological and evolutionary research.

High levels of fluctuating asymmetry in isolated stickleback populations.

BACKGROUND: Fluctuating asymmetry (FA), defined as small random deviations from the ideal bilateral symmetry, has been hypothesized to increase in response to both genetic and environmental stress experienced by a population. We compared levels of FA in 12 bilateral meristic traits (viz. lateral-line system neuromasts and lateral plates), and heterozygosity in 23 microsatellite loci, among four marine (high piscine predation risk) and four pond (zero piscine predation risk) populations of nine-spined sticklebacks (Pungitius pungitius). RESULTS: Pond sticklebacks had on average three times higher levels of FA than marine fish and this difference was highly significant. Heterozygosity in microsatellite markers was on average two times lower in pond (HE ≈ 0.3) than in marine (HE ≈ 0.6) populations, and levels of FA and heterozygosity were negatively correlated across populations. However, after controlling for habitat effect on heterozygosity, levels of FA and heterozygosity were uncorrelated. CONCLUSIONS: The fact that levels of FA in traits likely to be important in the context of predator evasion were elevated in ponds compared to marine populations suggests that relaxed selection for homeostasis in ponds lacking predatory fish may be responsible for the observed habitat difference in levels of FA. This inference also aligns with the observation that the levels of genetic variability across the populations did not explain population differences in levels of FA after correcting for habitat effect. Hence, while differences in strength of selection, rather than in the degree of genetic stress could be argued to explain habitat differences in levels of FA, the hypothesis that increased FA in ponds is caused by genetic stress cannot be rejected.

A high incidence of selection on physiologically important genes in the three-spined stickleback, Gasterosteus aculeatus.

Genome scan approaches to detect footprints of directional selection in the genomes of wild animal and plant populations have become popular tools to study local adaptation and speciation at the molecular level. Most studies thus far have used random molecular markers and found footprints of directional selection at, on average, 5% (range: 1-15%) of the examined loci. We focused on physiologically important genes that exhibit transcriptional responses to specific environmental or developmental conditions and assessed if these genes have been subject to directional selection and are responsible for local adaptation in the three-spined stickleback (Gasterosteus aculeatus). Using microsatellite markers located within or closely linked to (<6 kb) target genes, we investigated footprints of directional selection for 157 genes with known physiological functions in three marine and six freshwater populations. A high incidence (16.6%) of footprints of directional selection for these genes was revealed by four different outlier tests. In a subset of four populations screened with both physiologically important and random genes, footprints of directional selection were more frequent in physiologically important genes (13.4%) as compared with random genes (2.4%). In general, our findings indicate strong selective pressures on physiologically important genes, suggesting that these genes have significant functions in evolutionary adaptation to environmental heterogeneity.

Heterozygosity-behaviour correlations in nine-spined stickleback (Pungitius pungitius) populations: contrasting effects at random and functional loci.

The study of heterozygosity-fitness correlations (HFCs) has a long history in the fields of ecology and evolutionary biology but remains controversial. Recently, it has been shown that the genetic distance of markers from functional loci can be an important factor to be considered in addition to marker numbers and variability. In this study, we investigated the correlation between individual heterozygosity and behaviour (aggression, boldness and feeding activity) in nine-spined stickleback (Pungitius pungitius) individuals originating from four populations in two contrasting environments. Offspring of full-sib families raised in a common garden setting were assessed for behaviour and genotyped using 84 microsatellite markers that were either located within or near behaviourally or physiologically important genes (termed 'functional') or were randomly selected. No associations were detected with any behavioural trait in any population or over all populations when genetic variability was measured using all 84 markers combined. However, when the markers were separated into three functional categories (behavioural, physiological and random), several significant associations were observed both with functional markers and random markers in one of the four populations assessed. Interestingly, contrasting correlations with behaviour were observed when using physiological gene (negative) and random (positive) markers. Upon dividing the physiological gene markers into further subcategories based on their specific physiological functions, a strong relationship between the heterozygosity of markers linked to osmoregulation-related genes, and behaviour was revealed in the brackish water population. Our results indicate that both local (physiological) and general (neutral) effects are important in shaping behaviour and that heterozygosity-behaviour correlations are population dependent.

Spectral tuning by selective chromophore uptake in rods and cones of eight populations of nine-spined stickleback (Pungitius pungitius).

The visual pigments of rods and cones were studied in eight Fennoscandian populations of nine-spined stickleback (Pungitius pungitius). The wavelength of maximum absorbance of the rod pigment (λ(max)) varied between populations from 504 to 530 nm. Gene sequencing showed that the rod opsins of all populations were identical in amino acid composition, implying that the differences were due to varying proportions of chromophores A1 and A2. Four spectral classes of cones were found (two S-cones, M-cones and L-cones), correlating with the four classes of vertebrate cone pigments. For quantitative estimation of chromophore proportions, we considered mainly rods and M-cones. In four populations, spectra of both photoreceptor types indicated A2 dominance (population mean λ(max)=525-530 nm for rods and 535-544 nm for M-cones). In the four remaining populations, however, rod spectra (mean λ(max)=504-511 nm) indicated strong A1 dominance, whereas M-cone spectra (mean λ(max)=519-534 nm) suggested substantial fractions of A2. Quantitative analysis of spectra by three methods confirmed that rods and cones in these populations use significantly different chromophore proportions. The outcome is a shift of M-cone spectra towards longer wavelengths and a better match to the photic environment (light spectra peaking >560 nm in all the habitats) than would result from the chromophore proportions of the rods. Chromophore content was also observed to vary partly independently in M- and L-cones with potential consequences for colour discrimination. This is the first demonstration that selective processing of chromophore in rods and cones, and in different cone types, may be ecologically relevant.

High degree of sex chromosome differentiation in stickleback fishes.

BACKGROUND: Studies of closely related species with different sex chromosome systems can provide insights into the processes of sex chromosome differentiation and evolution. To investigate the potential utility of molecular markers in studying sex chromosome differentiation at early stages of their divergence, we examined the levels and patterns of genetic differentiation between sex chromosomes in nine-spined (Pungitius pungitius) and three-spined sticklebacks (Gasterosteus aculeatus) using microsatellite markers. RESULTS: A set of novel microsatellite markers spanning the entire length of the sex chromosomes were developed for nine-spined sticklebacks using the sequenced genomes of other fish species. Sex-specific patterns of genetic variability and male-specific alleles were identified at most of these loci, indicating a high degree of differentiation between the X and Y chromosomes in nine-spined sticklebacks. In three-spined sticklebacks, male-specific alleles were detected at some loci confined to two chromosomal regions. In addition, male-specific null alleles were identified at several other loci, implying the absence of Y chromosomal alleles at these loci. Overall, male-specific alleles and null alleles were found over a region spanning 81% of the sex chromosomes in three-spined sticklebacks. CONCLUSIONS: High levels but distinct patterns of sex chromosome differentiation were uncovered in the stickleback species that diverged 13 million years ago. Our results suggest that the Y chromosome is highly degenerate in three-spined sticklebacks, but not in nine-spined sticklebacks. In general, the results demonstrate that microsatellites can be useful in identifying the degree and patterns of sex chromosome differentiation in species at initial stages of sex chromosome evolution.

Identification of local- and habitat-dependent selection: scanning functionally important genes in nine-spined sticklebacks (Pungitius pungitius).

Understanding the selective forces promoting adaptive population divergence is a central issue in evolutionary biology. The role of environmental salinity in driving adaptation and evolution in aquatic organisms is still poorly understood. We investigated the relative impacts of habitat type (cf. saltwater vs. freshwater) and geographic area in shaping adaptive population divergence, as well as genes responsible for adaptation to different salinities in nine-spined sticklebacks (Pungitius pungitius). To this end, we employed a hitchhiking mapping approach with 111 microsatellite loci and one insertion/deletion locus including 63 loci situated within or close to genes with important physiological functions such as osmoregulation, growth, and thermal response. Using three pairs of marine and freshwater populations from different geographic areas, we identified several loci showing consistent evidence of being under directional selection in different outlier tests. Analyses of molecular variance at the loci under selection indicated that geographic area rather than habitat type has been acting as a central force in shaping adaptive population divergence. Nevertheless, both outlier tests and a spatial analysis method indicated that two loci (growth hormone receptor 2 and DEAD box polypeptide 56) are involved in adaptation to different habitats, implying that environmental salinity has been affecting them as a selective force. These loci are promising candidates for further investigations focusing on the molecular mechanisms of adaptation to marine and freshwater environments.

Species introduction promotes hybridization and introgression in Coregonus: is there sign of selection against hybrids?

Species introductions are considered one of the major drivers of biodiversity loss via ecological interactions and genetic admixture with local fauna. We examined two well-recognized fish species, native whitefish (Coregonus lavaretus) and introduced vendace (Coregonus albula), as well as their morphological hybrids in a single lake to test for selection against hybrids and backcrosses in the wild. A representative random subsample of 693 individuals (27.8%) was taken from the total catch of coregonids. This subsample was examined with the aim to select c. 50 individuals of pure whitefish (n = 52), pure vendace (n = 55) and putative hybrid (n = 19) for genetic analyses. The subsequent microsatellites and mitochondrial (mt) DNA analyses provided compelling evidence of hybridization and introgression. Of the 126 fish examined, four were found to be F(1) , 14 backcrosses to whitefish and seven backcrosses to vendace. The estimates of historical gene flow suggested higher rates from introduced vendace into native whitefish than vice versa, whereas estimates of contemporary gene flow were equal. Mitochondrial introgression was skewed, with 18 backcrosses having vendace mtDNA and only three with whitefish mtDNA. Hybrids and backcrosses had intermediate morphology and niche utilization compared with parental species. No evidence of selection against hybrids or backcrosses was apparent, as both hybrid and backcross growth rates and fecundities were high. Hybrids (F(1) ) were only detected in 2 year-classes, suggesting temporal variability in mating between vendace and whitefish. However, our data show that hybrids reached sexual maturity and reproduced actively, with backcrosses recorded from six consecutive year-classes, whereas no F(2) individuals were found. The results indicate widespread introgression, as 10.8% of coregonids were estimated to be backcrosses.

Utility of sequenced genomes for microsatellite marker development in non-model organisms: a case study of functionally important genes in nine-spined sticklebacks (Pungitius pungitius).

BACKGROUND: Identification of genes involved in adaptation and speciation by targeting specific genes of interest has become a plausible strategy also for non-model organisms. We investigated the potential utility of available sequenced fish genomes to develop microsatellite (cf. simple sequence repeat, SSR) markers for functionally important genes in nine-spined sticklebacks (Pungitius pungitius), as well as cross-species transferability of SSR primers from three-spined (Gasterosteus aculeatus) to nine-spined sticklebacks. In addition, we examined the patterns and degree of SSR conservation between these species using their aligned sequences. RESULTS: Cross-species amplification success was lower for SSR markers located in or around functionally important genes (27 out of 158) than for those randomly derived from genomic (35 out of 101) and cDNA (35 out of 87) libraries. Polymorphism was observed at a large proportion (65%) of the cross-amplified loci independently of SSR type. To develop SSR markers for functionally important genes in nine-spined sticklebacks, SSR locations were surveyed in or around 67 target genes based on the three-spined stickleback genome and these regions were sequenced with primers designed from conserved sequences in sequenced fish genomes. Out of the 81 SSRs identified in the sequenced regions (44,084 bp), 57 exhibited the same motifs at the same locations as in the three-spined stickleback. Di- and trinucleotide SSRs appeared to be highly conserved whereas mononucleotide SSRs were less so. Species-specific primers were designed to amplify 58 SSRs using the sequences of nine-spined sticklebacks. CONCLUSIONS: Our results demonstrated that a large proportion of SSRs are conserved in the species that have diverged more than 10 million years ago. Therefore, the three-spined stickleback genome can be used to predict SSR locations in the nine-spined stickleback genome. While cross-species utility of SSR primers is limited due to low amplification success, SSR markers can be developed for target genes and genomic regions using our approach, which should be also applicable to other non-model organisms. The SSR markers developed in this study should be useful for identification of genes responsible for phenotypic variation and adaptive divergence of nine-spined stickleback populations, as well as for constructing comparative gene maps of nine-spined and three-spined sticklebacks.

History vs. habitat type: explaining the genetic structure of European nine-spined stickleback (Pungitius pungitius) populations.

The genetic structure of contemporary populations can be shaped by both their history and current ecological conditions. We assessed the relative importance of postglacial colonization history and habitat type in the patterns and degree of genetic diversity and differentiation in northern European nine-spined sticklebacks (Pungitius pungitius), using mitochondrial DNA (mtDNA) sequences and 12 nuclear microsatellite and insertion/deletion loci. The mtDNA analyses identified - and microsatellite analyses supported - the existence of two historically distinct lineages (eastern and western). The analyses of nuclear loci among 51 European sites revealed clear historically influenced and to minor degree habitat dependent, patterns of genetic diversity and differentiation. While the effect of habitat type on the levels of genetic variation (coastal > freshwater) and differentiation (freshwater > coastal) was clear, the levels of genetic variability and differentiation in the freshwater sites were independent of habitat type (viz. river, lake and pond). However, levels of genetic variability, together with estimates of historical effective population sizes, decreased dramatically and linearly with increasing latitude. These geographical patterns of genetic variability and differentiation suggest that the contemporary genetic structure of freshwater nine-spined sticklebacks has been strongly impacted by the founder events associated with postglacial colonization and less by current ecological conditions (cf. habitat type). In general, the results highlight the strong and persistent effects of postglacial colonization history on genetic structuring of northern European fauna and provide an unparalleled example of latitudinal trends in levels of genetic diversity.

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 approximately 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 approximately 20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of approximately 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 evolution of sex determination associated with a chromosomal inversion.

Sex determination is a fundamentally important and highly diversified biological process, yet the mechanisms behind the origin of this diversity are mostly unknown. Here we suggest that the evolution of sex determination systems can be driven by a chromosomal inversion. We show that an XY system evolved recently in particular nine-spined stickleback (Pungitius pungitius) populations, which arose from ancient hybridization between two divergent lineages. Our phylogenetic and genetic mapping analyses indicate that the XY system is formed in a large inversion that is associated with hybrid sterility between the divergent lineages. We suggest that a new male-determining gene evolved in the inversion in response to selection against impaired male fertility in a hybridized population. Given that inversions are often associated with hybrid incompatibility in animals and plants, they might frequently contribute to the diversification of sex determination systems.

A four-questions perspective on public information use in sticklebacks (Gasterosteidae).

Whether learning primarily reflects general processes or species-specific challenges is a long-standing matter of dispute. Here, we present a comprehensive analysis of public information use (PI-use) in sticklebacks (Gasterosteidae). PI-use is a form of social learning by which animals are able to assess the relative quality of resources, here prey patches, by observing the behaviour of others. PI-use was highly specific with only Pungitius and their closest relative Culaea inconstans showing evidence of PI-use. We saw no effects of ontogenetic experience upon PI-use in Pungitius pungitius. Experiments with live demonstrators and animated fish revealed that heightened activity and feeding strikes by foraging conspecifics are important cues in the transmission of PI. Finally, PI-use was the only form of learning in which P. pungitius and another stickleback, Gasterosteus aculeatus differed. PI-use in sticklebacks is species-specific and may represent an 'ecological specialization' for social foraging. Whether this reflects selection on perception, attentional or cognitive processes remains to be determined.

Estimation of quantitative genetic parameters using marker-inferred relatedness in Japanese flounder: a case study of upward bias.

Marker-based methods for estimating heritability have been proposed as an effective means to study quantitative traits in long-lived organisms and natural populations. However, practical examinations to evaluate the usefulness and robustness of a regression method are limited. Using several quantitative traits of Japanese flounder Paralichthys olivaceus, the present study examined the influence of relatedness estimator and population structure on the estimation of heritability and genetic correlation under a regression method with 7 microsatellite loci. Significant heritability and genetic correlation were detected for several quantitative traits in 2 laboratory populations but not in a natural population. In the laboratory populations, upward bias in heritability appeared depending on the relatedness estimators and the populations. Upward bias in heritability increased with decreasing the actual variance of relatedness, suggesting that the estimates of heritability under the regression method tend to be overestimated due to the underestimation of the actual variance of relatedness. Therefore, relationship structure and precise estimation of relatedness are critical for applying this method.