Population genomics of seal lice provides insights into the postglacial history of northern European seals.

Genetic analyses of host-specific parasites can elucidate the evolutionary histories and biological features of their hosts. Here, we used population-genomic analyses of ectoparasitic seal lice (Echinophthirius horridus) to shed light on the postglacial history of seals in the Arctic Ocean and the Baltic Sea region. One key question was the enigmatic origin of relict landlocked ringed seal populations in lakes Saimaa and Ladoga in northern Europe. We found that that lice of four postglacially diverged subspecies of the ringed seal (Pusa hispida) and Baltic gray seal (Halichoerus grypus), like their hosts, form genetically differentiated entities. Using coalescent-based demographic inference, we show that the sequence of divergences of the louse populations is consistent with the geological history of lake formation. In addition, local effective population sizes of the lice are generally proportional to the census sizes of their respective seal host populations. Genome-based reconstructions of long-term effective population sizes revealed clear differences among louse populations associated with gray versus ringed seals, with apparent links to Pleistocene and Holocene climatic variation as well as to the isolation histories of ringed seal subspecies. Interestingly, our analyses also revealed ancient gene flow between the lice of Baltic gray and ringed seals, suggesting that the distributions of Baltic seals overlapped to a greater extent in the past than is the case today. Taken together, our results demonstrate how genomic information from specialized parasites with higher mutation and substitution rates than their hosts can potentially illuminate finer scale population genetic patterns than similar data from their hosts.

Phylogenomics and biogeography of sawflies and woodwasps (Hymenoptera, Symphyta).

Phylogenomic approaches have recently helped elucidate various insect relationships, but large-scale comprehensive analyses on relationships within sawflies and woodwasps are still lacking. Here, we infer the relationships and long-term biogeographic history of these hymenopteran groups using a large dataset of 354 UCE loci collected from 385 species that represent all major lineages. Early Hymenoptera started diversifying during the Early Triassic âˆ¼249 Ma and spread all over the ancient supercontinent Pangaea. We recovered Xyeloidea as a monophyletic sister group to other Hymenoptera and Pamphilioidea as sister to Unicalcarida. Within the diverse family Tenthredinidae, our taxonomically and geographically expanded taxon sampling highlights the non-monophyly of several traditionally defined subfamilies. In addition, the recent removal of Athalia and related genera from the Tenthredinidae into the separate family Athaliidae is supported. The deep historical biogeography of the group is characterised by independent dispersals and re-colonisations between the northern (Laurasia) and southern (Gondwana) palaeocontinents. The breakup of these landmasses led to ancient vicariance in several Gondwanan lineages, while interchange across the Northern Hemisphere has continued until the Recent. The little-studied African sawfly fauna is likewise a diverse mixture of groups with varying routes of colonization. Our results reveal interesting parallels in the evolution and biogeography of early hymenopterans and other ancient insect groups.

Loss of species and genetic diversity during colonization: Insights from acanthocephalan parasites in northern European seals.

Studies on host-parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple-enzyme restriction-site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population-genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra-host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies.

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections.

The Saimaa ringed seal (Pusa hispida saimensis) is endemic to Lake Saimaa in Finland. The subspecies is thought to have originated when parts of the ringed seal population of the Baltic region were trapped in lakes emerging due to postglacial bedrock rebound around 9000 years ago. During the 20th century, the population experienced a drastic human-induced bottleneck. Today encompassing a little over 400 seals with extremely low genetic diversity, it is classified as endangered. We sequenced sections of the mitochondrial control region from 60 up to 125-years-old museum specimens of the Saimaa ringed seal. The generated dataset was combined with publicly available sequences. We studied how genetic variation has changed through time in this subspecies and how it is phylogenetically related to other ringed seal populations from the Baltic Sea, Lake Ladoga, North America, Svalbard, and the White Sea. We observed temporal fluctuations in haplotype frequencies and loss of haplotypes accompanied by a recent reduction in female effective population size. In apparent contrast with the traditionally held view of the Baltic origin of the population, the Saimaa ringed seal mtDNA variation also shows affinities to North American ringed seals. Our results suggest that the Saimaa ringed seal has experienced recent genetic drift associated with small population size. The results further suggest that extant Baltic ringed seal is not representative of the ancestral population of the Saimaa ringed seal, which calls for re-evaluation of the deep history of this subspecies.

Palaearctic willow-catkin sawflies: a revision of the amentorum species group of Euura (Hymenoptera, Tenthredinidae).

The Euura amentorum species group is Holarctic, and in Europe it is most species-rich in the North. Their larvae develop entirely within the female catkins of Salix species: some species bore in the central stalk, whereas others live outside this and feed mainly on the developing seeds. Eight Palaearctic species are treated here as valid, and a key to these is provided. Males of five species are known. Two new species are described from northern Europe: Euura pohjola sp. n. and E. ursaminor sp. n. First records of E. itelmena (Malaise, 1931) from the West Palaearctic are presented. We propose seven new synonymies: Pontopristia montana Lindqvist, 1961 (junior secondary homonym in Euura) with Euura freyja (Liston, Taeger & Blank, 2009); Pontopristia brevilabris Malaise, 1921, Amauronematus fennicus Lindqvist, 1944, Pontopristia boreoalpina Lindqvist, 1961, Pontopristia punctulata Lindqvist, 1961, and Amauronematus pyrenaeus Lacourt, 1995 with Euura microphyes (Frster, 1854); and Pteronidea holmgreni Lindqvist, 1968 with Nematus umbratus Thomson, 1871. Lectotypes are designated for: Amauronematus fennicus Lindqvist, 1944, Nematus amentorum Frster, 1854, Nematus suavis Ruthe, 1859, Pontopristia brevilabris Malaise, 1921, Pontopristia itelmena Malaise, 1931, Pontopristia kamtchatica Malaise, 1931, Pontopristia lapponica Malaise, 1921, Pontopristia latiserra Malaise, 1921, Pontopristia romani Malaise, 1921, and Pristiphora amentorum var. nigripleuris Enslin, 1916. Many new host plant associations are recorded.

A curated DNA barcode reference library for parasitoids of northern European cyclically outbreaking geometrid moths.

Large areas of forests are annually damaged or destroyed by outbreaking insect pests. Understanding the factors that trigger and terminate such population eruptions has become crucially important, as plants, plant-feeding insects, and their natural enemies may respond differentially to the ongoing changes in the global climate. In northernmost Europe, climate-driven range expansions of the geometrid moths Epirrita autumnata and Operophtera brumata have resulted in overlapping and increasingly severe outbreaks. Delayed density-dependent responses of parasitoids are a plausible explanation for the 10-year population cycles of these moth species, but the impact of parasitoids on geometrid outbreak dynamics is unclear due to a lack of knowledge on the host ranges and prevalences of parasitoids attacking the moths in nature. To overcome these problems, we reviewed the literature on parasitism in the focal geometrid species in their outbreak range and then constructed a DNA barcode reference library for all relevant parasitoid species based on reared specimens and sequences obtained from public databases. The combined recorded parasitoid community of E. autumnata and O. brumata consists of 32 hymenopteran species, all of which can be reliably identified based on their barcode sequences. The curated barcode library presented here opens up new opportunities for estimating the abundance and community composition of parasitoids across populations and ecosystems based on mass barcoding and metabarcoding approaches. Such information can be used for elucidating the role of parasitoids in moth population control, possibly also for devising methods for reducing the extent, intensity, and duration of outbreaks.

High levels of inbreeding with spatial and host-associated structure in lice of an endangered freshwater seal.

Host-specialist parasites of endangered large vertebrates are in many cases more endangered than their hosts. In particular, low host population densities and reduced among-host transmission rates are expected to lead to inbreeding within parasite infrapopulations living on single host individuals. Furthermore, spatial population structures of directly-transmitted parasites should be concordant with those of their hosts. Using population genomic approaches, we investigated inbreeding and population structure in a host-specialist seal louse (Echinophthirius horridus) infesting the Saimaa ringed seal (Phoca hispida saimensis), which is endemic to Lake Saimaa in Finland, and is one of the most endangered pinnipeds in the world. We conducted genome resequencing of pairs of lice collected from 18 individual Saimaa ringed seals throughout the Lake Saimaa complex. Our analyses showed high genetic similarity and inbreeding between lice inhabiting the same individual seal host, indicating low among-host transmission rates. Across the lake, genetic differentiation among individual lice was correlated with their geographic distance, and assignment analyses revealed a marked break in the genetic variation of the lice in the middle of the lake, indicating substantial population structure. These findings indicate that movements of Saimaa ringed seals across the main breeding areas of the fragmented Lake Saimaa complex may in fact be more restricted than suggested by previous population-genetic analyses of the seals themselves.

Mitogenomics of the suborder Cottoidei (Teleostei: Perciformes): Improved assemblies, mitogenome features, phylogeny, and ecological implications.

We determined the mitogenome of Cyclopterus lumpus using a hybrid sequencing approach, and another four closely related species in the Liparidae based on available next-generation sequence data. We found that the mitogenome of C. lumpus was 17,266 bp in length, where the length and organisation were comparable to those reported for cottoids. However, we found a GC-homopolymer region in the intergenic space between tRNALeu2 and ND1 in liparids and cyclopterids. Phylogenetic reconstruction confirmed the monophyly of infraorders and firmly supported a sister-group relationship between Cyclopteridae and Liparidae. Purifying selection was the predominant force in the evolution of cottoid mitogenomes. There was significant evidence of relaxed selective pressures along the lineage of deep-sea fish, while selection was intensified in the freshwater lineage. Overall, our analysis provides a necessary expansion in the availability of mitogenomic sequences and sheds light on mitogenomic adaptation in Cottoidei fish inhabiting different aquatic environments.

Adipose gene expression profiles reveal insights into the adaptation of northern Eurasian semi-domestic reindeer (Rangifer tarandus).

Reindeer (Rangifer tarandus) are semi-domesticated animals adapted to the challenging conditions of northern Eurasia. Adipose tissues play a crucial role in northern animals by altering gene expression in their tissues to regulate energy homoeostasis and thermogenic activity. Here, we perform transcriptome profiling by RNA sequencing of adipose tissues from three different anatomical depots: metacarpal (bone marrow), perirenal, and prescapular fat in Finnish and Even reindeer (in Sakha) during spring and winter. A total of 16,212 genes are expressed in our data. Gene expression profiles in metacarpal tissue are distinct from perirenal and prescapular adipose tissues. Notably, metacarpal adipose tissue appears to have a significant role in the regulation of the energy metabolism of reindeer in spring when their nutritional condition is poor after winter. During spring, genes associated with the immune system are upregulated in the perirenal and prescapular adipose tissue. Blood and tissue parameters reflecting general physiological and metabolic status show less seasonal variation in Even reindeer than in Finnish reindeer. This study identifies candidate genes potentially involved in immune response, fat deposition, and energy metabolism and provides new information on the mechanisms by which reindeer adapt to harsh arctic conditions.

Sea ice reduction drives genetic differentiation among Barents Sea polar bears.

Loss of Arctic sea ice owing to climate change is predicted to reduce both genetic diversity and gene flow in ice-dependent species, with potentially negative consequences for their long-term viability. Here, we tested for the population-genetic impacts of reduced sea ice cover on the polar bear (Ursus maritimus) sampled across two decades (1995-2016) from the Svalbard Archipelago, Norway, an area that is affected by rapid sea ice loss in the Arctic Barents Sea. We analysed genetic variation at 22 microsatellite loci for 626 polar bears from four sampling areas within the archipelago. Our results revealed a 3-10% loss of genetic diversity across the study period, accompanied by a near 200% increase in genetic differentiation across regions. These effects may best be explained by a decrease in gene flow caused by habitat fragmentation owing to the loss of sea ice coverage, resulting in increased inbreeding of local polar bears within the focal sampling areas in the Svalbard Archipelago. This study illustrates the importance of genetic monitoring for developing adaptive management strategies for polar bears and other ice-dependent species.

DNA barcoding reveals different cestode helminth species in northern European marine and freshwater ringed seals.

Three subspecies of the ringed seal (Pusa hispida) are found in northeastern Europe: P. h. botnica in the Baltic Sea, P. h saimensis in Lake Saimaa in Finland, and P. h. ladogensis in Lake Ladoga in Russia. We investigated the poorly-known cestode helminth communities of these closely related but ecologically divergent subspecies using COI barcode data. Our results show that, while cestodes from the Baltic Sea represent Schistocephalus solidus, all worms from the two lakes are identified as Ligula intestinalis, a species that has previously not been reported from seals. The observed shift in cestode communities appears to be driven by differential availability of intermediate fish host species in marine vs. freshwater environments. Both observed cestode species normally infect fish-eating birds, so further work is required to elucidate the health and conservation implications of cestode infections in European ringed seals, whether L. intestinalis occurs also in marine ringed seals, and whether the species is able to reproduce in seal hosts. In addition, a deep barcode divergence found within S. solidus suggests the presence of cryptic diversity under this species name.

Evolutionary transitions in broad tapeworms (Cestoda: Diphyllobothriidea) revealed by mitogenome and nuclear ribosomal operon phylogenetics.

Broad tapeworms (Diphyllobothriidea) are parasites whose adults are capable of infecting a wide range of freshwater, marine and terrestrial tetrapods including humans. Previous works examining the evolution of habitat and host use in this group have been hampered by the lack of a well-resolved phylogeny. In order to produce a robust phylogenetic framework for diphyllobothriideans, we sequenced the complete mitochondrial genome of 13 representatives, carefully chosen to cover the major clades, and two outgroup species representing the Spathebothriidea and Haplobothriidea. In addition, complementary data from the nuclear ribosomal operon was sequenced for 10 representative taxa. Mitogenomes and ssrDNA and lsrDNA were used towards elucidating the phylogenetic framework for the Diphyllobothriidea. The Cephalochlamydidae is confirmed as the earliest diverging diphyllobothriidean lineage, and Solenophoridae and Diphyllobothriidae are sister groups. We infer a probable freshwater origin of the diphyllobothriideans. The ancestral condition for life cycle complexity could not be unambiguously resolved. However, we infer exclusive use of a three-host life cycle following the origin of the Solenophoridae + Diphyllobothriidae. Regarding definitive host use, although we infer reptiles as the most likely ancestral condition, this result should be revisited with a more densely sampled phylogeny in future studies. Freshwater habitat is used by the early diverging lineages within the Solenophoridae + Diphyllobothriidae clade. For the latter, habitat use shifts between freshwater and marine environments, and definitive host use includes marine and terrestrial mammals and birds. We use mitochondrial genomes to distinguish Schistocephalus species occurring in different species of sticklebacks and demonstrate conspecificity of Ligula cf. intestinalis specimens collected from two Fennoscandian ringed seal subspecies.

Patterns of Microbiome Variation Among Infrapopulations of Permanent Bloodsucking Parasites.

While interspecific variation in microbiome composition can often be readily explained by factors such as host species identity, there is still limited knowledge of how microbiomes vary at scales lower than the species level (e.g., between individuals or populations). Here, we evaluated variation in microbiome composition of individual parasites among infrapopulations (i.e., populations of parasites of the same species living on a single host individual). To address this question, we used genome-resolved and shotgun metagenomic data of 17 infrapopulations (balanced design) of the permanent, bloodsucking seal louse Echinophthirius horridus sampled from individual Saimaa ringed seals Pusa hispida saimensis. Both genome-resolved and read-based metagenomic classification approaches consistently show that parasite infrapopulation identity is a significant factor that explains both qualitative and quantitative patterns of microbiome variation at the intraspecific level. This study contributes to the general understanding of the factors driving patterns of intraspecific variation in microbiome composition, especially of bloodsucking parasites, and has implications for understanding how well-known processes occurring at higher taxonomic levels, such as phylosymbiosis, might arise in these systems.

Microbiomes of willow-galling sawflies: effects of host plant, gall type, and phylogeny on community structure and function.

While free-living herbivorous insects are thought to harbor microbial communities composed of transient bacteria derived from their diet, recent studies indicate that insects that induce galls on plants may be involved in more intimate host-microbe relationships. We used 16S rDNA metabarcoding to survey larval microbiomes of 20 nematine sawfly species that induce bud or leaf galls on 13 Salix species. The 391 amplicon sequence variants (ASVs) detected represented 69 bacterial genera in six phyla. Multi-variate statistical analyses showed that the structure of larval microbiomes is influenced by willow host species as well as by gall type. Nevertheless, a "core" microbiome composed of 58 ASVs is shared widely across the focal galler species. Within the core community, the presence of many abundant, related ASVs representing multiple distantly related bacterial taxa is reflected as a statistically significant effect of bacterial phylogeny on galler-microbe associations. Members of the core community have a variety of inferred functions, including degradation of phenolic compounds, nutrient supplementation, and production of plant hormones. Hence, our results support suggestions of intimate and diverse interactions between galling insects and microbes and add to a growing body of evidence that microbes may play a role in the induction of insect galls on plants.

Genomes of the willow-galling sawflies Euura lappo and Eupontania aestiva (Hymenoptera: Tenthredinidae): a resource for research on ecological speciation, adaptation, and gall induction.

Hymenoptera is a hyperdiverse insect order represented by over 153,000 different species. As many hymenopteran species perform various crucial roles for our environments, such as pollination, herbivory, and parasitism, they are of high economic and ecological importance. There are 99 hymenopteran genomes in the NCBI database, yet only five are representative of the paraphyletic suborder Symphyta (sawflies, woodwasps, and horntails), while the rest represent the suborder Apocrita (bees, wasps, and ants). Here, using a combination of 10X Genomics linked-read sequencing, Oxford Nanopore long-read technology, and Illumina short-read data, we assembled the genomes of two willow-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae: Euurina): the bud-galling species Euura lappo and the leaf-galling species Eupontania aestiva. The final assembly for E. lappo is 259.85 Mbp in size, with a contig N50 of 209.0 kbp and a BUSCO score of 93.5%. The E. aestiva genome is 222.23 Mbp in size, with a contig N50 of 49.7 kbp and a 90.2% complete BUSCO score. De novo annotation of repetitive elements showed that 27.45% of the genome was composed of repetitive elements in E. lappo and 16.89% in E. aestiva, which is a marked increase compared to previously published hymenopteran genomes. The genomes presented here provide a resource for inferring phylogenetic relationships among basal hymenopterans, comparative studies on host-related genomic adaptation in plant-feeding insects, and research on the mechanisms of plant manipulation by gall-inducing insects.

Genetic modification of the flavonoid pathway alters growth and reveals flexible responses to enhanced UVB - Role of foliar condensed tannins.

Accumulation of certain phenolics is a well-known response of plants to enhanced UVB radiation (280-315 nm), but few experiments have compared the relative importance of different phenolic groups for UVB resilience. To study how an altered phenolic profile affects the responses and resilience of silver birch (Betula pendula) to enhanced UVB, we used RNA interference (RNAi) targeting dihydroflavonol reductase (DFR), anthocyanidin synthase (ANS), or anthocyanidin reductase (ANR) to change the accumulation of phenolics. The unmodified control line and RNAi-modified plants were grown for 51 days under ambient or +32% enhanced UVB dose in a greenhouse. RNAi greatly affected phenolic profile and plant growth. There were no interactive effects of RNAi and UVB on growth or photosynthesis, which indicates that the RNAi and unmodified control plants were equally resilient. UVB enhancement led to an accumulation of foliar flavonoids and condensed tannins, and an increase in the density of stem glands and glandular trichomes on upper leaf surfaces in both the control and RNAi-modified plants. Our results do not indicate a photoprotective role for condensed tannins. However, decreased growth of high-flavonoid low-tannin DFRi and ANRi plants implies that the balance of flavonoids and condensed tannins might be important for normal plant growth.

Determinants of genetic variation across eco-evolutionary scales in pinnipeds.

The effective size of a population (Ne), which determines its level of neutral variability, is a key evolutionary parameter. Ne can substantially depart from census sizes of present-day breeding populations (NC) as a result of past demographic changes, variation in life-history traits and selection at linked sites. Using genome-wide data we estimated the long-term coalescent Ne for 17 pinniped species represented by 36 population samples (total n = 458 individuals). Ne estimates ranged from 8,936 to 91,178, were highly consistent within (sub)species and showed a strong positive correlation with NC ([Formula: see text] = 0.59; P = 0.0002). Ne/NC ratios were low (mean, 0.31; median, 0.13) and co-varied strongly with demographic history and, to a lesser degree, with species' ecological and life-history variables such as breeding habitat. Residual variation in Ne/NC, after controlling for past demographic fluctuations, contained information about recent population size changes during the Anthropocene. Specifically, species of conservation concern typically had positive residuals indicative of a smaller contemporary NC than would be expected from their long-term Ne. This study highlights the value of comparative population genomic analyses for gauging the evolutionary processes governing genetic variation in natural populations, and provides a framework for identifying populations deserving closer conservation attention.

Genome- and transcriptome-derived microsatellite loci in lumpfish Cyclopterus lumpus: molecular tools for aquaculture, conservation and fisheries management.

The lumpfish Cyclopterus lumpus is commercially exploited in numerous areas of its range in the North Atlantic Ocean, and is important in salmonid aquaculture as a biological agent for controlling sea lice. Despite the economic importance, few genetic resources for downstream applications, such as linkage mapping, parentage analysis, marker-assisted selection (MAS), quantitative trait loci (QTL) analysis, and assessing adaptive genetic diversity are currently available for the species. Here, we identify both genome- and transcriptome-derived microsatellites loci from C. lumpus to facilitate such applications. Across 2,346 genomic contigs, we detected a total of 3,067 microsatellite loci, of which 723 were the most suitable ones for primer design. From 116,555 transcriptomic unigenes, we identified a total of 231,556 microsatellite loci, which may indicate a high coverage of the available STRs. Out of these, primer pairs could only be designed for 6,203 loci. Dinucleotide repeats accounted for 89 percent and 52 percent of the genome- and transcriptome-derived microsatellites, respectively. The genetic composition of the dominant repeat motif types showed differences from other investigated fish species. In the genome-derived microsatellites AC/GT (67.8 percent), followed by AG/CT (15.1 percent) and AT/AT (5.6 percent) were the major motifs. Transcriptome-derived microsatellites showed also most dominantly the AC/GT repeat motif (33 percent), followed by A/T (26.6 percent) and AG/CT (11 percent). Functional annotation of microsatellite-containing transcriptomic sequences showed that the majority of the expressed sequence tags encode proteins involved in cellular and metabolic processes, binding activity and catalytic reactions. Importantly, STRs linked to genes involved in immune system process, growth, locomotion and reproduction were discovered in the present study. The extensive genomic marker information reported here will facilitate molecular ecology studies, conservation initiatives and will benefit many aspects of the breeding programmes of C. lumpus.