|Home | About | Journals | Submit | Contact Us | Français|
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.
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.
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.
Recent advances in our understanding of the physiological and molecular functions of genes have paved the road for investigating functional genomic variation associated with adaptation and speciation in the wild [1,2]. Consequently, targeting specific genes and genomic regions of interest - rather than random genomic regions - holds a great promise as a shortcut to identify genes involved in phenotypic variation and adaptive divergence [3-5]. Despite a steadily increasing number of completed genome sequences, genomic resources and tools are still very limited for the vast majority of non-model organisms. Therefore, ability to develop molecular markers in or around target genes is essential for application of this approach for non-model organisms. In addition, molecular markers associated with functionally important genes are useful in construction of comparative genetic maps, in which they can be exploited as comparative anchor tagged sequence loci [6,7].
Microsatellites or simple sequence repeats (SSRs) are highly abundant in eukaryotic genomes, accounting for 3-5% of the mammalian genomes [8,9]. Owing to their wide genomic distribution, codominant inheritance and hypervariability, they are widely recognized as one of the most powerful molecular markers in the field of genetics. As a result of the widespread use of SSRs, substantial efforts have been made to devise procedures for developing SSR markers [10,11]. In addition, cross-species transfer of SSR primers is commonly attempted in many taxa . However, SSR markers developed with conventional approaches are derived from the genome more or less in a random manner. Expressed sequence tags (ESTs) are commonly used as an alternative to genomic libraries as a source of SSR markers . SSR markers derived from ESTs have some advantages over those developed from genomic libraries because EST-derived markers can associate with genes of known or putative function, and they exhibit relatively high transferability between closely related species [13-17]. However, SSRs are generally much less abundant in transcribed regions than in non-transcribed regions [18-21] and found typically only in a few percentage of ESTs [15,22-26]. Besides, designing primers requires sufficient flanking sequences, resulting in a considerable reduction in number of ESTs available to develop SSR markers [23,24,26,27]. Therefore, even if a large EST database is available for a target species, ESTs have limitations as a material for development of SSR markers for specific genes.
One way to obtain SSR markers for specific genes and genomic regions in a given species is to use SSR primers developed for the closest relative with a sequenced genome. However, an obvious limitation of this approach is that mutations in SSR flanking sequences will inhibit cross-species amplification success - a problem that is likely to attenuate with an increasing divergence time . In general, success of cross-species transfer is a negative function of the evolutionary distance separating the source and focal species [28-31]. Another crucial issue is related to evolution and persistence of SSRs among different species. Investigations of SSR conservation have demonstrated that several SSRs are retained not only in closely related species, but also in species that have diverged more than 100 million years ago [32-35]. Nevertheless, comprehensive surveys of SSR conservation using aligned sequences of different species have rarely been reported [36,37], making it difficult to estimate the patterns and degree of SSR conservation in different taxa.
For the reasons elaborated above, development of SSR markers for target genes and genomic regions in non-model organisms is challenging. Yet, while the closest relative with a sequenced genome is too distantly related to the focal species, one can take an advantage of the increasing number of completed genomes for different species. For instance, as often used in species for which no direct species-specific sequence information is available, conserved sequences in specific genes and genomic regions of interest can be used to design primer sequences applicable to a wide variety of organisms [e.g. ].
Teleosts consist of approximately 28,000 species , which correspond to more than half of all living vertebrates. Despite a number of features of evolutionary interest and economical importance, genomic resources and tools are still lacking for most teleost taxa. Currently, genome sequences are available for five species - zebrafish (Danio rerio), three-spined stickleback (Gasterosteus aculeatus), medaka (Oryzias latipes), spotted green pufferfish (Tetraodon nigroviridis) and fugu (Takifugu rubripes) . The development of genome sequences for three-spined sticklebacks has made great contribution to an understanding of the genetic architecture of several phenotypic traits [41-44]. Because three-spined and nine-spined (Pungitius pungitius) sticklebacks exhibit similar ecological and morphological characteristics , these species provide an opportunity to study whether the same genes or genomic regions are responsible for phenotypic variation of certain traits and adaptive divergence in different lineages. This would facilitate a molecular understanding of the parallel evolution of these species, which have diverged more than 10 million years ago - equivalent to 5-10 millions of generations [46,47]. A potentially effective strategy to this end would be to develop SSR markers targeting functionally important genes.
The main objective of this study was to develop a large set of SSR markers targeting specific genes and genomic regions for a non-model organism - the nine-spined stickleback - in which genome sequences and ESTs are not yet available. To this end, two strategies were adopted. First, we tested cross-species utility of 158 SSR primer sets for functionally important genes originally developed in three-spined sticklebacks together with 188 SSR markers derived from genomic libraries and ESTs. Secondly, we investigated the potential utility of available sequenced fish genomes to develop SSR markers for functionally important genes in nine-spined sticklebacks. To address prospects for this approach, the patterns and degree of SSR conservation were examined in three-spined and nine-spined sticklebacks using their aligned sequences.
Out of the 158 SSR markers for functionally important genes (gene-based SSRs), 27 showed robust and specific amplification within the expected size range in nine-spined sticklebacks (Table (Table1,1, see also Additional files 1 and 2), resulting in a low level (17.1%) of cross-species amplification. In contrast, amplification success was 34.7% (35 out of 101) and 40.2% (35 out of 87) in the SSR markers derived from genomic libraries (genomic SSRs) and ESTs (EST-derived SSRs), respectively (Table (Table1).1). The tendency for higher amplification success with the EST-derived SSRs than with the genomic SSRs is in agreement with the results of previous studies [48-51] - finding which has been explained by high sequence conservation in coding regions [13-15].
Factors affecting cross-species amplification success were assessed using the 388 SSR markers. A hierarchical generalized linear model (GLM) revealed a significant influence of SSR type (cf. gene-based, genomic vs. EST-derived SSRs) on amplification success (F2,331 = 8.28, P = 0.016). In addition, amplification success was significantly affected by primer site (cf. exonic, intronic, intergenic vs. other combinations; F3,331 = 16.43, P < 0.001). Across the three SSR types, amplification success was high for the SSR markers in which both forward and reverse primers were located in exonic regions (62.5%, 15 out of 24), whereas it was lower if primers were located either in intronic (20.4%, 20 out of 98) or intergenic regions (26.3%, 42 out of 160; Table Table1).1). This effect was particularly obvious for the gene-based SSRs with intronic primers, in which case the amplification success was very low (3.8%, two out of 53; Table Table1).1). As for EST-derived SSRs, trinucleotide SSRs are the most abundant repeat motif in ESTs and tend to be found in coding regions, whereas dinucleotide SSRs are often found in untranslated regions . The fact that most of the EST-derived SSRs used in our study are dinucleotide repeats (85 out of 87) suggests that a number of the EST-derived SSRs might be located in untranslated regions. While in theory EST-derived SSRs should be located in exonic regions, 54 (out of 87) SSRs were located in intergenic regions according to the Ensembl genebuild. This inconsistency could be due to artifacts such as prediction errors and contamination of cDNA libraries with genomic DNA. Nevertheless, the result that amplification success tended to be higher for the EST-derived SSRs with exonic primers (50.0%, seven out of 14) than for those with intergenic primers (36.2%, 17 out of 47) might, at least in part, result from the fact that sequence homology is less in untranslated regions and increase toward the start codon of the coding regions in related species .
While the effects of SSR type and primer site were significant, amplification success was not significantly associated with average primer length (F1,331 = 1.09, P = 0.297) or average GC content (F1,331 = 3.12, P = 0.077). Similarly, amplification success was independent of differences in GC content (F1,331 = 0.00, P = 0.993) and melting temperature (F1,331 = 0.53, P = 0.465) between primers within a given primer pair. In addition, there was no association between amplification success and expected PCR product size (F1,331 = 0.393, P = 0.531). While the positive effect of average melting temperature appeared to be significant (F1,331 = 5.82, P = 0.016), no clear difference of melting temperature was found among the SSR types. Based on these results, it is unlikely that differential amplification success among the three sets of SSR markers stemmed from different primer conditions. In fact, all of the SSR markers were successfully amplified under the same PCR and DNA conditions in three-spined sticklebacks. Multiple gene copies are known to exist in several functionally important genes [52-55]. Rather than primer conditions, divergence of functionally important genes might be the cause of low cross-species amplification success of the gene-based SSRs.
In the gene-based SSRs, polymorphism was found at 16 out of the 27 amplifying loci (59.3%; Table Table1).1). This rate was similar to that observed in the genomic SSRs (74.3%, 26 out of 35) and EST-derived SSRs (60.0%, 21 out of 35; Table Table1).1). In total, 63 out of the 97 amplified loci exhibited polymorphism in Fennoscandian populations (Additional file 1). For the amplified loci, incidence of polymorphism was independent of SSR type (GLM, F2,90 = 3.37, P = 0.185), SSR location (cf. exonic, intronic vs. intergenic regions; F2,90 = 5.11, P = 0.078) and SSR repeat motif (cf. di- vs. trinucleotide repeats; F1,90 = 0.02, P = 0.876). The relatively high proportion of polymorphic loci across the different SSR types and SSR locations suggests that several SSRs are conserved in three-spined and nine-spined sticklebacks. For the polymorphic loci of gene-based SSRs, an average of 8.8 alleles per locus (range = 2-38) were identified in the three populations (Additional file 1). This value was equivalent to that obtained in the genomic SSRs (9.3) and EST-derived SSRs (8.3; Additional file 1). Average heterozygosity varied from 0.19 to 0.55 in the gene-based SSRs, from 0.07 to 0.64 in the genomic SSRs and from 0.14 to 0.63 in the EST-derived SSRs among the three populations (Additional file 1). MICRO-CHECKER analyses did not indicate the presence of null alleles, with the possible exceptions of the CLCN7, GS1, Gac7080P and Stn18 in the Baltic Sea, the Stn127 and GAest41 in the Lake 1 and the GAest16 in the Pyöreälampi. There was no evidence for deviations from Hardy-Weinberg equilibrium at any locus in any of the populations.
In general, our results demonstrate that cross-species utility of SSR primers for functionally important genes is less efficient as compared to that of genomic and EST-derived SSR markers. This is attributed to limited amplification success rather than a low incidence of polymorphism. Therefore, the development of species-specific primers would be necessary for obtaining SSR markers for functionally important genes.
To investigate the potential utility of available sequenced fish genomes for SSR marker development in nine-spined sticklebacks, we surveyed SSRs within and around 67 functionally important genes in the three-spined stickleback genome and designed 70 primer sets for amplification and sequencing of these SSR regions using the conserved sequences determined by sequenced fish genomes (Additional file 3). The PCR product size of respective genomic regions obtained in nine-spined sticklebacks was concordant with that estimated from the three-spined stickleback genome (Additional file 4). All of the sequences of nine-spined sticklebacks for the 70 regions exhibited the highest BLAST hit scores in the target regions and high homologies to the sequences of the three-spined stickleback (Additional file 4). Out of the 70 genomic regions representing 44,084 bp, 49 contained at least one SSR in nine-spined sticklebacks (Additional file 4). The total number of SSRs observed in these regions was 81, including 9 mono-, 52 di-, 18 tri- and two tetranucleotide motifs (Figure (Figure1).1). In the three-spined stickleback genome, 96 SSRs were found in the 70 homologous regions, including 12 mono-, 57 di-, 23 tri- and four tetranucleotide motifs (Figure (Figure1).1). Out of the 81 SSRs found in nine-spined sticklebacks, 64 were identified at the same locations as those of the three-spined stickleback (Additional file 4). In addition, 57 out of the 64 SSRs exhibited the same motifs as those of the three-spined stickleback (Figure (Figure1,1, see also Additional file 4), indicating that a large proportion of SSRs are conserved in the genomes of these species. Our results also demonstrated that SSRs with di- and trinucleotide repeat motifs are highly conserved but those with mononucleotide repeat motifs are less so (Figure (Figure1).1). Hence, the level of SSR conservation may differ among SSRs differing in repeat motif type.
To further address SSR conservation in stickleback species, we investigated if SSRs randomly derived from genomic libraries of Pungitius species are found at the homologous genomic locations of three-spined sticklebacks. For this analysis, we used publicly available SSR and flanking sequences of Pungitius pungitius (i.e. nine-spined stickleback)  and Pungitius sp.  - so called the Omono-type, which has been regarded as an independent species from Pungitius pungitius based on the biological species concept . The 13 Pungitius pungitius sequences (5,310 bp) contained one mono- and 16 dinucleotide motif SSRs. In the three-spined stickleback genome, 18 SSRs were identified in the homologous regions, including one mono- and 17 dinucleotide motifs. Out of the 18 SSRs identified in the three-spined stickleback, 15 (83.3%) exhibited the same motifs at the same locations as in Pungitius pungitius (Figure (Figure2).2). In the 19 Pungitius sp. sequences (4,117 bp) containing 20 dinucleotide motif SSRs, 17 SSRs were identified in the homologous regions of the three-spined stickleback genome, including one mono- and 16 dinucleotide motifs (Figure (Figure2).2). Out of the 17 SSRs identified in the three-spined stickleback, 15 (88.2%) exhibited the same motifs at the same locations as in Pungitius sp. (Figure (Figure2).2). The comparative analyses of randomly selected Pungitius SSRs in the three-spined stickleback genome further indicated a high degree of SSR conservation in stickleback species.
While several studies have reported conservation of single SSRs between different taxa [32-35,59], a comprehensive survey of SSR conservation is limited to a comparison of human (Homo sapiens) and chimpanzee (Pan troglodytes) [36,37], which have diverged six million years ago . According to Vowles and Amos , 70% of human SSRs are homologues in chimpanzees. Our results demonstrated that a similar proportion (70% for randomly selected motifs) of SSRs are retained in three-spined sticklebacks and Pungitius species despite longer divergence time (cf. more than 10 million years)  and much shorter generation times (cf. one or two years) .
Based on the sequences obtained in nine-spined sticklebacks, species-specific primer sets were designed to amplify 58 SSRs targeting 57 functionally important genes (Table (Table2).2). Among them, polymorphism was identified at 41 loci (Table (Table2)2) in Fennoscandian populations. On average, 7.7 alleles per locus (range = 2-27) were identified across the three populations (Table (Table3).3). Average heterozygosity was 0.57 in the Baltic Sea, 0.37 in the Lake 1 and 0.06 in the Pyöreälampi (Table (Table3).3). There was no indication for the presence of null alleles, with the possible exceptions of the Ppgm40 and Ppgm50 in the Baltic Sea and the Ppgm52 and Ppgm56 in the Lake 1. Deviations from Hardy-Weinberg equilibrium were not observed at any locus in any of the populations.
The level of SSR variability is known to be associated with repeat motifs due to their different mutation rates [61,62]. In addition, cross-species transfer of SSR primers often results in a lower level of SSR variability in a focal species relative to a source species because of ascertainment bias [63,64]. We investigated the patterns and degree of SSR variability using three Fennoscandian populations. Across the 104 polymorphic loci identified in this study, an average number of alleles per locus and average heterozygosity were 7.4 and 0.60 in the Baltic Sea, 2.6 and 0.32 in the Lake 1, and 1.6 and 0.10 in the Pyöreälampi, respectively. As expected, the levels of SSR variability were significantly dependent on population. The genome-wide survey indicated that genetic variation of the Pyöreälampi is very low, as also shown in a previous study with 11 SSR and one insertion/deletion loci . In our data set, the levels of SSR variability were not dependent on marker origin, SSR type and SSR repeat motif (Table (Table4).4). However, a significant influence of SSR location on the levels of allele number and heterozygosity was apparent (Table (Table4).4). Across the three populations, average allele number and heterozygosity were 4.7 and 0.41 in exonic regions, 3.3 and 0.30 in intronic regions and 4.2 and 0.36 in intergenic regions, respectively. While the level of SSR variability is known to differ between coding and untranslated regions [66,67], EST-derived SSRs tend to show lower variability than genomic SSRs [50,68]. These differences are thought to arise due to heterogeneous distributions of SSR repeat motifs. However, the higher variability in exonic SSRs than in other SSRs is not explainable by an artifact stemming from different repeat motifs because a majority of the polymorphic SSRs were dimeric repeats independently of their location. Several lines of evidence suggest that SSR variation may affect various traits and be subject to natural selection [21,69,70]. While the potential effect of variable mutation rates can not be ruled out, the heterogeneous distribution of SSR variability observed in this study might be ascribable to natural selection.
Our study demonstrated that a large proportion of SSRs are conserved in the stickleback species which have diverged from a common ancestor more than 10 million years ago . Therefore, the three-spined stickleback genome can be used to predict SSR locations in Pungitius species. Our results also suggest that the main limitation of cross-species utility of SSR markers lies in the failure of amplification success probably due to mutations in SSR flanking sequences. While it is possible to predict to some degree the likelihood of amplification success based on the information of primer binding sites, cross-species transferability of SSR primers for functionally important genes is particularly low as compared to that of genomic and EST-derived SSR primers. Yet, SSR markers can be developed for functionally important genes and target genomic regions using the approach outlined in this paper. This approach should be applicable also to other non-model organisms. The SSR markers developed for functionally important genes should be useful to identify genes responsible for phenotypic variation and adaptive divergence in nine-spined sticklebacks, as well as for constructing comparative gene maps of nine-spined and three-spined sticklebacks.
Nine-spined sticklebacks collected from the Baltic Sea (coastal; 60°12' N, 25°11' E), the 'Lake 1' (lake; 67°54' N, 20°50' E) and the Pyöreälampi (pond; 66°16' N, 29°26' E) were used in this study. The fish were sampled with seine nets or minnow traps in 2002 (Lake 1) and 2008 (Baltic Sea and Pyöreälampi). Total DNA was extracted from fin clips stored in 70-99% ethanol with a phenol-chloroform method  following proteinase K digestion.
Cross-species utility of three-spined stickleback SSR primers was tested for 158 SSR markers for physiologically important genes (gene-based SSRs) [Y. Shimada, T. Shikano and J. Merilä, unpublished] coupled with 101 markers derived from genomic libraries (genomic SSRs) and 87 markers derived from ESTs (EST-derived SSRs; Additional files 1 and 2) [72-76]. The genomic and EST-derived SSRs were classified according to the source information deposited in GenBank . The following factors potentially affecting cross-species amplification success were scored for each of the makers: SSR marker type (cf. gene-based, genomic and EST-derived SSRs), primer binding site (cf. exonic, intronic, intergenic and other combinations), average primer length, average and difference of GC content and melting temperature in forward and reverse primer pairs, as well as expected PCR product size. Primer binding sites were categorized into exonic, intronic and intergenic regions based on the Ensembl genebuild in the three-spined stickleback genome . Since information on untranslated regions was not available for a number of genes, we did not distinguish between coding and untranslated regions in the analyses. The primer parameters were calculated using BioEdit  under the actual PCR conditions (see below). The expected PCR product sizes were calculated based on the three-spined stickleback genome. The role of these factors was evaluated using generalized linear models as implemented in JMP 5 (SAS Inst. Inc.). In these tests, amplification success was treated as a binary dependent variable (successful amplification = 1, failed amplification = 0), SSR type and primer site as factors, and other parameters as covariates. Logit link function was used. For the successfully amplified loci, factors affecting incidence of polymorphism were evaluated using generalized linear models treating SSR type, SSR location (cf. exonic, intronic and intergenic regions) and SSR repeat motif (cf. di- and trinucleotide repeats) as factors. SSR location was categorized into exonic, intronic and intergenic regions based on the Ensembl genebuild in the three-spined stickleback genome. In this test, polymorphic locus was treated as a binary dependent variable (polymorphic = 1, monomorphic = 0) using logit link function.
Based on the literature on gene functions in teleosts, we selected 67 genes responsible for significant physiological - such as osmoregulation, thermal response, growth, disease and taste [Y. Shimada, T. Shikano and J. Merilä, unpublished] - and developmental functions [e.g. ] (Additional file 3). Genomic locations of these genes were identified in the three-spined stickleback genome following Shimada et al. [Y. Shimada, T. Shikano and J. Merilä, unpublished]. In brief, we searched the three-spined stickleback ESTs which correspond to target genes of this species or other teleosts in the GenBank database  and mapped them in the three-spined stickleback genome. The genomic range of respective genes was determined according to the Ensembl transcript and Genscan predictions (Additional file 3). Since the genomic region of the PITX1 was not available due to partially incomplete sequences of the three-spined stickleback genome, the sequence of this gene (GenBank: AY517634.1) was used.
SSRs were searched in the target genes and their flanking regions in the three-spined stickleback genome using Tandem repeats finder . In order to survey conserved regions for designing amplification and sequencing primers in nine-spined sticklebacks, the sequences of these genomic regions were subject to BLASTN searches against the currently available genome sequences of other teleosts, i.e. medaka, fugu, spotted green pufferfish and/or zebrafish . Conserved regions were determined by aligning the sequences of three-spined sticklebacks and those of other fish species detected by the BLASTN searches. Based on the location of SSRs and conserved regions, primer sequences for nine-spined sticklebacks were designed manually in one genomic region for each target gene, except for the GHRI and IGF-I, for which two and three regions were used, respectively (Additional file 3).
Two individuals of the Pyöreälampi were used for amplifying and sequencing the target genomic regions. One three-spined stickleback individual from the Baltic Sea (60°12' N, 25°11' E) was used as a positive control. Using a primer pair for respective target regions (Additional file 3), PCR amplifications were carried out in a 20 μl reaction volume consisting of 1× PCR buffer (Bioline), 1.5 mM MgCl2, 0.25 mM dNTP (Finnzymes), 0.15 U BIOTAQ DNA polymerase (Bioline), 5 pmol of each primer and approx. 40 ng of genomic DNA. The reactions were performed as follows: an initial degeneration step at 95°C for 3 min, followed by 30 s at 95°C, 30 s at 53-60°C and 60-120 s at 72°C for 35 cycles with a final extension at 72°C for 5 min (see Additional file 3 for optimal PCR conditions in each primer pair). Approximate size of the PCR amplicons was determined by electrophoresis on 1.5% agarose gel with a DNA ladder (GeneRuler™ DNA Ladder Mix, Fermentas). PCR products were purified using exonuclease I (New England Biolabs) and shrimp alkaline phosphatase (Roche) and directly sequenced in both forward and reverse directions with the same primers as those used in the PCRs. The sequencing reactions were performed using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) according to manufacture's instructions. Cycle sequencing products were purified by ethanol precipitation and analyzed on an ABI 3730xl DNA Analyzer (Applied Biosystems).
The sequences in forward and reverse directions of two individuals were aligned using CLUSTAL W  as implemented in MEGA 4  and edited by hand. For large PCR amplicons (≥1200 bp), the sequences in forward and reverse directions were separately aligned using two individuals (Additional file 4). As sequences were available only in one direction for four genomic regions even after retrials, sequences for these regions were aligned using two individuals (Additional file 4). The sequences were subject to BLASTN searches against the three-spined stickleback genome to ensure that they are mapped back to the correct locations in the genome. The homologous sequences in three-spined and nine-spined sticklebacks were aligned to compare SSR locations and motifs between them. This comparison was performed using SSRs with minimum repeat numbers of ten, five and four for one (mono-), two (di-) and three or longer (tri-, tetra-, penta- and hexanucleotide) repeat motifs, respectively. To further address SSR conservation in stickleback species, we also investigated if SSRs randomly derived from genomic libraries of Pungitius species are found at the homologous locations in the three-spined stickleback genome using publicly available SSR and flanking sequences of Pungitius pungitius (i.e. nine-spined stickleback) [; GenBank: AB473819-AB473831] and Pungitius sp. (Omono type) [; GenBank: AB300827-AB300851]. Out of the 38 sequences, six (GenBank: AB300830, AB300831, AB300841, AB300842, AB300844, AB300849) were excluded from the analyses because of low BLAST hit scores and alignment problems. To develop SSR markers for functionally important genes, primer sets were designed based on the sequences of nine-spined sticklebacks using WebSat . Primer sequences were deposited in GenBank under accession numbers GU553378-GU553434.
For the cross-species amplification test of three-spined stickleback primers, amplification success and polymorphism were determined using the following three-step procedure. Firstly, amplification was tested using four individuals from the Pyöreälampi and Baltic Sea (two individuals per population) with fluorescent labelled forward primers (FAM, HEX or TET) and GTTT-tailed reverse primers . As a positive control, one individual of the three-spined stickleback was used. PCRs were performed under optimal conditions for three-spined sticklebacks and conducted in a 10 μl reaction volume consisting of 1× PCR buffer (Bioline), 1.5 mM MgCl2, 0.2 mM dNTP (Finnzymes), 0.18 U BIOTAQ DNA polymerase (Bioline), 5 pmol of each primer and approx. 20 ng of template DNA. The reactions were performed as follows: an initial degeneration step at 95°C for 3 min, followed by 30 s at 95°C, 30 s at 53°C and 30 s at 72°C for 30 cycles with a final extension at 72°C for 5 min. Amplification success was determined by electrophoresis on 1.6% agarose gel. Secondly, for the loci that showed robust and specific amplification within the expected size range, polymorphism was investigated by genotyping 24 individuals from the Baltic Sea and the Lake 1 (12 individuals per population). For efficient screening, PCRs were carried out using the Qiagen Multiplex PCR Kit (Qiagen) in 10 μl reaction volumes containing 1× Qiagen Multiplex PCR Master Mix, 0.5× Q-Solution, 2 pmol of each primer and approx. 20 ng of template DNA. The reactions were performed by the following cycle: an initial activation step at 95°C for 15 min, followed by 30 s at 94°C, 90 s at 53°C and 60 s at 72°C for 30 cycles with a final extension at 60°C for 5 min. PCR products were visualized with a MegaBACE 1000 automated sequencer (Amersham Biosciences) and their sizes were determined with ET-ROX 550 size standard (Amersham Biosciences). Thirdly, genetic variability of the polymorphic loci identified with the 24 individuals was evaluated by genotyping a total of 24 individuals from each of the Baltic Sea, Lake 1 and Pyöreälampi populations using multiplex PCRs. For the SSR primers developed in nine-spined sticklebacks, polymorphism and genetic variability were evaluated following the procedures for the second and third steps. Since some of the SSR markers yielded clearer allele profiles at an annealing temperature of 55°C, this temperature was used for these loci instead of 53°C (see Table Table22 for an optimal annealing temperature for each primer pair). Alleles were scored using Fragment Profiler 1.2 (Amersham Biosciences) with visual inspection and manual corrections of alleles.
Locus and population specific gene diversities (HE)  were estimated using FSTAT 2.9.3 [85,86]. Within population and locus specific FIS were estimated for each population to detect possible deviations from Hardy-Weinberg equilibrium with 10 000 permutations using FSTAT 2.9.3. Sequential Bonferroni corrections  were applied to minimize type I errors. The presence of null alleles was tested using MICRO-CHECKER .
Factors affecting the levels of genetic variation were evaluated with general linear models where allele number or heterozygosity was treated as a dependent variable, SSR marker origin (cf. three-spined and nine-spined sticklebacks), SSR type, SSR location and SSR repeat motif (cf. mono-, di-, tri- and tetranucleotide repeats) as fixed factors and population as a random factor. These analyses were performed with JMP 5.
TS conceived the study, contributed to the gene selection and localization, prepared the molecular data, conducted the analyses and wrote the manuscript. JR selected the genes and conducted the molecular work. YS made contributions to the gene selection and localization. JM advised on the statistical analyses and contributed to writing the manuscript. All authors read and approved the final manuscript.
Genetic variability of three-spined stickleback SSR markers in three populations of nine-spined sticklebacks.
List of three-spined stickleback SSR markers that were not polymorphic or successfully amplified in nine-spined sticklebacks.
Target genes and their function and location in the theree-spined stickleback genome, and PCR and sequencing primers for nine-spined sticklebacks.
Homology of nine-spined stickleback sequences in the three-spined stickleback genome and comparative genomic location of SSRs in these species.
We thank people at the Oulanka Research Station, Abigel Gonda, Aki Hirvonen, John Loehr and Jarmo Saarikivi for help in obtaining samples. Thanks are also due to Craig Primmer for comments and checking the English. Our study was supported by the Academy of Finland and the Japan Society for the Promotion of Science.