First de novo whole genome sequencing and assembly of the pink-footed goose.

Annotated genomes can provide new perspectives on the biology of species. We present the first de novo whole genome sequencing for the pink-footed goose. In order to obtain a high-quality de novo assembly the strategy used was to combine one short insert paired-end library with two mate-pair libraries. The pink-footed goose genome was assembled de novo using three different assemblers and an assembly evaluation was subsequently performed in order to choose the best assembler. For our data, ALLPATHS-LG performed the best, since the assembly produced covers most of the genome, while introducing the fewest errors. A total of 26,134 genes were annotated, with bird species accounting for virtually all BLAST hits. We also estimated the substitution rate in the pink-footed goose, which can be of use in future demographic studies, by using a comparative approach with the genome of the chicken, the mallard and the swan goose. A substitution rate of 1.38×10-7 per nucleotide per generation was obtained when comparing the genomes of the two closely-related goose species (the pink-footed and the swan goose). Altogether, we provide a valuable tool for future genomic studies aiming at particular genes and regions of the pink-footed goose genome as well as other bird species.

Demographic inference from whole-genome and RAD sequencing data suggests alternating human impacts on goose populations since the last ice age.

We investigated how population changes and fluctuations in the pink-footed goose might have been affected by climatic and anthropogenic factors. First, genomic data confirmed the existence of two separate populations: western (Iceland) and eastern (Svalbard/Denmark). Second, demographic inference suggests that the species survived the last glacial period as a single ancestral population with a low population size (100-1,000 individuals) that split into the current populations at the end of the last glacial maximum with Iceland being the most plausible glacial refuge. While population changes during the last glaciation were clearly environmental, we hypothesize that more recent demographic changes are human-related: (1) the inferred population increase in the Neolithic is due to deforestation to establish new lands for agriculture, increasing available habitat for pink-footed geese, (2) the decline inferred during the Middle Ages is due to human persecution, and (3) improved protection explains the increasing demographic trends during the 20th century. Our results suggest both environmental (during glacial cycles) and anthropogenic effects (more recent) can be a threat to species survival.

Assessing pre- and post-zygotic barriers between North Atlantic eels (Anguilla anguilla and A. rostrata).

Elucidating barriers to gene flow is important for understanding the dynamics of speciation. Here we investigate pre- and post-zygotic mechanisms acting between the two hybridizing species of Atlantic eels: Anguilla anguilla and A. rostrata. Temporally varying hybridization was examined by analyzing 85 species-diagnostic single-nucleotide polymorphisms (SNPs; FST ⩾0.95) in eel larvae sampled in the spawning region in the Sargasso Sea in 2007 (N=92) and 2014 (N=460). We further investigated whether genotypes at these SNPs were nonrandomly distributed in post-F1 hybrids, indicating selection. Finally, we sequenced the mitochondrial ATP6 and nuclear ATP5c1 genes in 19 hybrids, identified using SNP and restriction site associated DNA (RAD) sequencing data, to test a previously proposed hypothesis of cytonuclear incompatibility leading to adenosine triphosphate (ATP) synthase dysfunction and selection against hybrids. No F1 hybrids but only later backcrosses were observed in the Sargasso Sea in 2007 and 2014. This suggests that interbreeding between the two species only occurs in some years, possibly controlled by environmental conditions at the spawning grounds, or that interbreeding has diminished through time as a result of a declining number of spawners. Moreover, potential selection was found at the nuclear and the cytonuclear levels. Nonetheless, one glass eel individual showed a mismatch, involving an American ATP6 haplotype and European ATP5c1 alleles. This contradicted the presence of cytonuclear incompatibility but may be explained by that (1) cytonuclear incompatibility is incomplete, (2) selection acts at a later life stage or (3) other genes are important for protein function. In total, the study demonstrates the utility of genomic data when examining pre- and post-zyotic barriers in natural hybrids.

Non-parallel divergence across freshwater and marine three-spined stickleback Gasterosteus aculeatus populations.

This work investigated whether multiple freshwater populations of three-spined stickleback Gasterosteus aculeatus in different freshwater catchments in the Jutland Peninsula, Denmark, derived from the same marine populations show repeated adaptive responses. A total of 327 G. aculeatus collected at 13 sampling locations were screened for genetic variation using a combination of 70 genes putatively under selection and 26 neutral genes along with a marker linked to the ectodysplasin gene (eda), which is strongly correlated with plate armour morphs in the species. A highly significant genetic differentiation was found that was higher among different freshwater samples than between marine-freshwater samples. Tests for selection between marine and freshwater populations showed a very low degree of parallelism and no single nucleotide polymorphism was detected as outlier in all freshwater-marine pairwise comparisons, including the eda. This suggests that G. aculeatus is not necessarily the prime example of parallel local adaptation suggested in much of the literature and that important exceptions exist (i.e. the Jutland Peninsula). While marine populations in the results described here showed a high phenotype-genotype correlation at eda, a low association was found for most of the freshwater populations. The most extreme case was found in the freshwater Lake Hald where all low-plated phenotypes were either homozygotes for the allele supposed to be associated with completely plated morphs or heterozygotes, but none were homozygotes for the putative low-plated allele. Re-examination of data from seven G. aculeatus studies agrees in showing a high but partial association between phenotype-genotype at eda in G. aculeatus freshwater populations and that mismatches occur everywhere in the European regions studied (higher in some areas, i.e. Denmark). This is independent of the eda marker used.

Temporal changes in allele frequencies in a small marble trout Salmo marmoratus population threatened by extreme flood events.

The effect of extreme floods on the genetic composition of marble trout Salmo marmoratus living in Lipovscek, a tributary of the Soca River in Slovenia, which has been affected by multiple destructive flood events for centuries was investigated. By monitoring genetic variability during the period 2004-2011, apparent signatures of genetic erosion including a decline in observed and expected heterozygosities and allelic richness were observed. Contemporary effective population size was estimated between 11 and 55 individuals, which is congruent with census data. The data suggest asymmetric gene flow between the two sections of the river. The existence of substantial downstream migration (15-19%) was confirmed by paternity analysis. A small (1-3%) upstream migration was also suggested, which was confirmed by tagging data. Overall, low genetic diversity has not prevented the survival of the Lipovscek population, which might be a common feature of salmonid freshwater populations.

Signatures of natural selection between life cycle stages separated by metamorphosis in European eel.

BACKGROUND: Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels. RESULTS: We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle. CONCLUSIONS: The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Genome-wide single-generation signatures of local selection in the panmictic European eel.

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Speciation and demographic history of Atlantic eels (Anguilla anguilla and A. rostrata) revealed by mitogenome sequencing.

Processes leading to speciation in oceanic environments without obvious physical barriers remain poorly known. European and American eel (Anguilla anguilla and A. rostrata) spawn in partial sympatry in the Sargasso Sea. Larvae are advected by the Gulf Stream and other currents towards the European/North African and North American coasts, respectively. We analyzed 104 mitogenomes from the two species along with mitogenomes of other Anguilla and outgroup species. We estimated divergence time between the two species to identify major events involved in speciation. We also considered two previously stated hypotheses: one where the ancestral species was present in only one continent but was advected across the Atlantic by ocean current changes and another where population declines during Pleistocene glaciations led to increasing vicariance, facilitating speciation. Divergence time was estimated to ∼3.38 Mya, coinciding with the closure of the Panama Gateway that led to reinforcement of the Gulf Stream. This could have advected larvae towards European/North African coasts, in which case American eel would be expected to be the ancestral species. This scenario could, however, not be unequivocally confirmed by analyses of dN/dS, nucleotide diversity and effective population size estimates. Extended bayesian skyline plots showed fluctuations of effective population sizes and declines during glaciations, and thus also lending support to the importance of vicariance during speciation. There was evidence for positive selection at the ATP6 and possibly ND5 genes, indicating a role in speciation. The findings suggest an important role of ocean current changes in speciation of marine organisms.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms.

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Phylogenetic relationships and demographic histories of the Atherinidae in the Eastern Atlantic and Mediterranean Sea re-examined by Bayesian inference.

The aim of our study is to examine the phylogenetic relationship, divergence times and demographic history of the five close-related Mediterranean and North-eastern Atlantic species/forms of Atherina using the full Bayesian framework for species tree estimation recently implemented in ∗BEAST. The inference is made possible by multilocus data using three mitochondrial genes (12S rRNA, 16S rRNA, control region) and one nuclear gene (rhodopsin) from multiple individuals per species available in GenBank. Bayesian phylogenetic analysis of the complete gene dataset produced a tree with strong support for the monophyly of each species, as well as high support for higher level nodes. An old origin of the Atherina group was suggested (19.2 MY), with deep split events within the Atherinidae predating the Messinian Salinity Crisis. Regional genetic substructuring was observed among populations of A. boyeri, with AMOVA and MultiDimensional Scaling suggesting the existence of five groupings (Atlantic/West Mediterranean, Adriatic, Greece, Black Sea and Tunis). The level of subdivision found might be consequence of the hydrographic isolation within the Mediterranean Sea. Bayesian inference of past demographic histories showed a clear signature of demographic expansion for the European coast populations of A. presbyter, possibly linked to post-glacial colonizations, but not for the Azores/Canary Islands, which is expected in isolated populations because of the impossibility of finding new habitats. Within the Mediterranean, signatures of recent demographic expansion were only found for the Adriatic population of A. boyeri, which could be associated with the relatively recent emergence of the Adriatic Sea.

Positive Darwinian selection in gamete recognition proteins of Strongylocentrotus sea urchins.

Gamete recognition proteins commonly experience positive Darwinian selection and evolve more rapidly than nonreproductive proteins, but the selective forces responsible for their adaptive diversification remain unclear. We examined the patterns of positive selection in the cognate interacting pair of proteins formed by sperm bindin and its egg receptor (EBR1) and in two regions of the sea urchin sperm receptor for egg jelly suREJ3 gene (exons 22 and 26) among four species of Strongylocentrotus sea urchins (S. purpuratus, S. droebachiensis, S. pallidus and S. franciscanus). The signatures of selection differed at each reproductive protein. A strong signal of positive selection was detected at bindin in all lineages even though the species compared had highly variable gamete traits and experience different intensities and forms of sexual selection and sexual conflict in nature. Weaker selection was observed at EBR1 but the small region studied precluded a clear understanding of the extent of sexual conflict between bindin and the EBR1 protein. At the suREJ3 locus, diversifying selection was observed in exon 22 but not exon 26, suggesting that these regions experience different selective pressures and evolutionary constraints. Positive selection was also detected within S. pallidus at suREJ-22 because of the presence of 12 amino acid replacement mutations segregating at frequencies >0.10. Our results suggest that sexual conflict may be the predominant evolutionary mechanism driving the rapid diversification of reproductive proteins between, and polymorphism within, strongylocentrotid sea urchins.

Inferring the demographic history of the Adriatic Flexopecten complex.

Phylogenetic analysis of DNA sequences can be a useful tool in taxonomical studies, which has important implications in terms of species traceability. The aim of our study is to solve the controversy regarding the taxonomical status of the Flexopecten group, which includes two putative sister species (F. glaber and F. proteus) that co-occur in the Adriatic Sea and are clearly distinguishable on the basis of morphometric characters. Our molecular analysis using three mitochondrial genes (COI, 12S rRNA and 16S rRNA) and one nuclear gene (Histone H3) suggest that F. glaber and F. proteus are the same species, as evidenced by both putative species appearing mixed in all genetic trees with no clustering according to species. Using a Bayesian approach, we inferred the demographic history of the Flexopecten group, which suggests that first F. glaber occurred in the Mediterranean, then F. glaber colonized the Adriatic sometime in the last 18,000 years, and finally the F. proteus morph appeared only recently in the Adriatic Sea. We propose F. proteus to be synonymized with F. glaber, which should have priority and be used in the future.

Genetic variability is unrelated to growth and parasite infestation in natural populations of the European eel (Anguilla anguilla).

Positive correlations between individual genetic heterozygosity and fitness-related traits (HFCs) have been observed in organisms as diverse as plants, marine bivalves, fish or mammals. HFCs are not universal and the strength and stability of HFCs seem to be variable across species, populations and ages. We analysed the relationship between individual genetic variability and two different estimators of fitness in natural samples of European eel, growth rate (using back-calculated length-at-age 1, 2 and 3) and parasite infestation by the swimbladder nematode Anguillicola crassus. Despite using a large data set of 22 expressed sequence tags-derived microsatellite loci and a large sample size of 346 individuals, no heterozygote advantage was observed in terms of growth rate or parasite load. The lack of association was evidenced by (i) nonsignificant global HFCs, (ii) a Multivariate General Linear Model showing no effect of heterozygosity on fitness components, (iii) single-locus analysis showing a lower number of significant tests than the expected false discovery rate, (iv) sign tests showing only a significant departure from expectations at one component, and, (v) a random distribution of significant single-locus HFCs that was not consistent across fitness components or sampling sites. This contrasts with the positive association observed in farmed eels in a previous study using allozymes, which can be explained by the nature of the markers used, with the allozyme study including many loci involved in metabolic energy pathways, while the expressed sequence tags-linked microsatellites might be located in genes or in the proximity of genes uncoupled with metabolism/growth.

Genetic composition of Atlantic and Mediterranean recruits of European eel Anguilla anguilla based on EST-linked microsatellite loci.

Anguilla anguilla glass eels arriving at two Mediterranean and two Atlantic sites were tested for differences in genetic composition between regions using a total of 23 microsatellite loci developed from an expressed sequence tag (EST) library. Hierarchical analysis of molecular variance indicated a non-significant difference between regions (Mediterranean v. Atlantic), which contrasted with the significant differences observed between samples within regions. The existence of a single spawning site for all A. anguilla individuals and extensive migration loop with great opportunity for mixing of individuals might explain the homogeneity in genetic composition found between regions. The observation of a (small-scale) pattern of genetic patchiness among intra-annual samples (arrival waves) within geographic regions does not conflict with the lack of (large-scale) geographic sub-structuring found between the Mediterranean and Atlantic regions, but most likely is a consequence of the strong dependence of A. anguilla on oceanic conditions in the Sargasso Sea that might result in a limited parental contribution to each spawning event. The comparison of Atlantic and Mediterranean A. anguilla glass eel recruits based on EST-linked microsatellite loci provides evidence supporting the hypothesis of panmixia A. anguilla across Europe.

Growth rate correlates to individual heterozygosity in the European eel, Anguilla anguilla L.

Heterozygosity-fitness correlations (HFCs) have been reported in populations of many species. We provide evidence for a positive correlation between genetic variability and growth rate at 12 allozyme loci in a catadromous marine fish species, the European eel (Anguilla anguilla L.). More heterozygous individuals show a significantly higher length and weight increase and an above average condition index in comparison with more homozygous individuals. To a lesser extent, six microsatellite loci show a similar pattern, with positive but not significant correlations between heterozygosity and growth rate. The HFCs observed could be explained by an effect of either direct allozyme over-dominance or associative overdominance. Selection affecting some of the allozyme loci would explain the greater strength of the HFCs found at allozymes in comparison with microsatellites and the lack of correlation between MLH at allozymes and MLH at microsatellites. Associative overdominance (where allozyme loci are merely acting as neutral markers of closely linked fitness loci) might provide an explanation for the HFCs if we consider that allozyme loci have a higher chance than microsatellites to be in linkage disequilibrium with fitness loci.

Population genetic structure in farm and feral American mink (Neovison vison) inferred from RAD sequencing-generated single nucleotide polymorphisms.

Feral American mink populations (), derived from mink farms, are widespread in Europe. In this study we investigated genetic diversity and genetic differentiation between feral and farm mink using a panel of genetic markers (194 SNP) generated from RAD sequencing data. Sampling included a total of 211 individuals from 14 populations, 4 feral and 10 from farms, the latter including a total of 7 color types (Brown, Black, Mahogany, Sapphire, White, Pearl, and Silver). Our study revealed similar low levels of genetic diversity in both farm and feral mink. Results are consistent with small effective population size as a consequence of line selection in the farms and founder effects of a few escapees from the farms in feral populations. Moderately high genetic differentiation was found between farm and feral animals, suggesting a scenario in which wild populations were founded from farm escapes a few decades ago. Currently, escapes and gene flow are probably limited. Genetic differentiation was higher among farm color types than among farms, consistent with line selection using few individuals to create the lines. Finally, no indications of inbreeding were found in either farm or feral samples, with significant negative values found in most farm samples, showing farms are successful in avoiding inbreeding.

Species and hybrid identification of sturgeon caviar: a new molecular approach to detect illegal trade.

Overexploitation of wild populations due to the high economic value of caviar has driven sturgeons to near extinction. The high prices commanded by caviar on world markets have made it a magnet for illegal and fraudulent caviar trade, often involving low-value farmed caviar being sold as top-quality caviar. We present a new molecular approach for the identification of pure sturgeon species and hybrids that are among the most commercialized species in Europe and North America. Our test is based on the discovery of species-specific single nucleotide polymorphisms (SNPs) in the ribosomal protein S7, supplemented with the Vimentin gene and the mitochondrial D-loop. Test validations performed in 702 specimens of target and nontarget sturgeon species demonstrated a 100% identification success for Acipenser naccarii, A. fulvescens, A. stellatus, A. sinensis and A. transmontanus. In addition to species identification, our approach allows the identification of Bester and AL hybrids, two of the most economically important hybrids in the world, with 80% and 100% success, respectively. Moreover, the approach has the potential to identify many other existing sturgeon hybrids. The development of a standardized sturgeon identification tool will directly benefit trade law enforcement, providing the tools to monitor and regulate the legal trade of caviar and protect sturgeon stocks from illicit producers and traders, hence contributing to safeguarding this group of heavily threatened species.

A resource of genome-wide single-nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel.

Reduced representation genome sequencing such as restriction-site-associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single-nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome-wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome-wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19,703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long-term effective population size was estimated to range between 132,000 and 1,320,000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82,425 loci and 376,918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.

Detecting genome-wide gene transcription profiles associated with high pollution burden in the critically endangered European eel.

The European eel illustrates an example of a critically endangered fish species strongly affected by human stressors throughout its life cycle, in which pollution is considered to be one of the factors responsible for the decline of the stock. The objective of our study was to better understand the transcriptional response of European eels chronically exposed to pollutants in their natural environment. A total of 42 pre-migrating (silver) female eels from lowly, highly and extremely polluted environments in Belgium and, for comparative purposes, a lowly polluted habitat in Italy were measured for polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs) and brominated flame retardants (BFRs). Multipollutant level of bioaccumulation was linked to their genome-wide gene transcription using an eel-specific array of 14,913 annotated cDNAs. Shared responses to pollutant exposure were observed when comparing the highly polluted site in Belgium with the relatively clean sites in Belgium and Italy. First, an altered pattern of transcription of genes was associated with detoxification, with a novel European eel CYP3A gene and gluthatione S-transferase transcriptionally up-regulated. Second, an altered pattern of transcription of genes associated with the oxidative phosphorylation pathway, with the following genes involved in the generation of ATP being transcriptionally down-regulated in individuals from the highly polluted site: NADH dehydrogenase, succinate dehydrogenase, ubiquinol-cytochrome c reductase, cytochrome c oxidase and ATP synthase. Although we did not measure metabolism directly, seeing that the transcription level of many genes encoding enzymes involved in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated in the highly polluted site suggests that pollutants may have a significant effect on energy metabolism in these fish.

Isolation and characterization of expressed sequence tag-linked microsatellite loci for the European eel (Anguilla anguilla).

A European eel (Anguilla anguilla) expressed sequence tag database consisting of 795 contigs and 4008 singletons was screened for microsatellites sequences. Primers were designed to amplify 96 repeats, of which 86 gave good quality amplification products. Twenty-eight microsatellites were selected for further microsatellite genotyping. Only two loci were found to be monomorphic; out of the 26 polymorphic loci, number of alleles per locus ranged from two to 14, while the observed and expected heterozygosities ranged from 0.05 to 0.93, and from 0.05 to 0.95, respectively. All 28 primer sets tested revealed positive amplification in American eel (Anguilla rostrata).