Takahito Shikano

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.

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.

Global analysis of genes involved in freshwater adaptation in threespine sticklebacks (Gasterosteus aculeatus).

Examples of parallel evolution of phenotypic traits have been repeatedly demonstrated in threespine sticklebacks (Gasterosteus aculeatus) across their global distribution. Using these as a model, we performed a targeted genome scan—focusing on physiologically important genes potentially related to freshwater adaptation—to identify genetic signatures of parallel physiological evolution on a global scale. To this end, 50 microsatellite loci, including 26 loci within or close to (<6 kb) physiologically important genes, were screened in paired marine and freshwater populations from six locations across the Northern Hemisphere. Signatures of directional selection were detected in 24 loci, including 17 physiologically important genes, in at least one location. Although no loci showed consistent signatures of selection in all divergent population pairs, several outliers were common in multiple locations. In particular, seven physiologically important genes, as well as reference ectodysplasin gene (EDA), showed signatures of selection in three or more locations. Hence, although these results give some evidence for consistent parallel molecular evolution in response to freshwater colonization, they suggest that different evolutionary pathways may underlie physiological adaptation to freshwater habitats within the global distribution of the threespine stickleback.

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.

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.

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.

Construction of Ultradense Linkage Maps with Lep-MAP2: Stickleback F2 Recombinant Crosses as an Example

High-density linkage maps are important tools for genome biology and evolutionary genetics by quantifying the extent of recombination, linkage disequilibrium, and chromosomal rearrangements across chromosomes, sexes, and populations. They provide one of the best ways to validate and refine de novo genome assemblies, with the power to identify errors in assemblies increasing with marker density. However, assembly of high-density linkage maps is still challenging due to software limitations. We describe Lep-MAP2, a software for ultradense genome-wide linkage map construction. Lep-MAP2 can handle various family structures and can account for achiasmatic meiosis to gain linkage map accuracy. Simulations show that Lep-MAP2 outperforms other available mapping software both in computational efficiency and accuracy. When applied to two large F2-generation recombinant crosses between two nine-spined stickleback (Pungitius pungitius) populations, it produced two high-density (∼6 markers/cM) linkage maps containing 18,691 and 20,054 single nucleotide polymorphisms. The two maps showed a high degree of synteny, but female maps were 1.5–2 times longer than male maps in all linkage groups, suggesting genome-wide recombination suppression in males. Comparison with the genome sequence of the three-spined stickleback (Gasterosteus aculeatus) revealed a high degree of interspecific synteny with a low frequency (<5%) of interchromosomal rearrangements. However, a fairly large (ca. 10 Mb) translocation from autosome to sex chromosome was detected in both maps. These results illustrate the utility and novel features of Lep-MAP2 in assembling high-density linkage maps, and their usefulness in revealing evolutionarily interesting properties of genomes, such as strong genome-wide sex bias in recombination rates.

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)

BackgroundIdentification 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.ResultsCross-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.ConclusionsOur 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.

Quantitative Genetics of Body Size and Timing of Maturation in Two Nine-Spined Stickleback (Pungitius pungitius) Populations

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.

Phylogeography and Genetic Structuring of European Nine-Spined Sticklebacks (Pungitius pungitius)—Mitochondrial DNA Evidence

As a consequence of colonisation from different glacial refugia, many northern European taxa are split into distinct western and eastern lineages. However, as for the nine-spined stickleback (Pungitius pungitius), the exact location of the contact zone between lineages often remains poorly known. We assessed the genetic differentiation and diversity in the nine-spined stickleback within Europe using 1037 base pairs of cytochrome b sequence for 320 individuals from 57 locations, including pond, lake, river, and coastal habitats. Our main aims were (i) to locate the contact zone between the previously recognized western and eastern lineages, (ii) investigate latitudinal patterns in genetic diversity, (iii) compare genetic diversity among different habitat types, and (iv) date the known split between eastern and western lineages. The data revealed the split between eastern and western to be located across the Danish Straits and roughly following the Norway/Sweden border to the North. Reference sites from Canada form their own clades, and one of the Canadian sites was found to have a haplotype common to the Eastern European lineage, possibly representing an ancestral polymorphism. The split between the two European clades was dated to approximately 1.48 million years ago (Mya), and between Canada and Europe to approximately 1.62 Mya. After controlling for habitat effects, nucleotide (but not haplotype) diversity across populations decreased with increasing latitude. Coastal populations showed significantly higher haplotype diversity (but not nucleotide diversity) than pond populations, but there were no detectable differences in haplotype diversity among different freshwater habitat types (viz. river, lake and pond populations), or between coastal and lake/river populations. Sequences were found to cluster according to their geographic proximity, rather than by habitat type, and all habitat types were found within each major clade, implying that colonisation and adaptation between the coastal and freshwater environments in different regions must have occurred in parallel.