Species status and population structure of mussels (Mollusca: Bivalvia: Mytilus spp.) in the Wadden Sea of Lower Saxony (Germany)

Three species of mussel (genus Mytilus) occur in Europe: M. edulis (Linnaeus 1758), M. galloprovincialis (Lamarck 1819) and M. trossulus (Gould, Boston Society of Natural History 3: 343?C348, 1850). Although these species are indigenous to the North Sea, the Mediterranean and the Baltic Sea, respectively, they form an extended patchy species complex along the coasts of Europe (??the Mytilus edulis complex??) and are able to hybridize where their distributions overlap. Recent studies examining the taxonomic status and genetic composition of Mytilus populations in the Netherlands and the British Isles have revealed introgressive hybridization processes within this species complex, with hints of an invasion of nonindigenous M. galloprovincialis into the North Sea. Furthermore, an extensive international mussel fishery industry in Europe (i.e., Great Britain, the Netherlands, Denmark, and Germany) is also in discussion for a possibly anthropogenically induced bioinvasion of nonindigenous Mytilus traits into the Wadden Sea area. Although it is assumed that the Wadden Sea of Germany comprises M. edulis only, this has never been confirmed in a molecular genetic study. To assess the situation for the Wadden Sea of Lower Saxony, we conducted the first molecular study of the Mytilus genus in the region. Taxonomic identification of 504 mussels from 13 intertidal mussel banks using the nDNA marker Me15/16 revealed a population composition of 99% M. edulis and 1% M. edulis X M. galloprovincialis hybrids. Hence, the Wadden Sea population is unaffected by range expansion of nonindigenous Mytilus traits. The genetic structure of the M. edulis populations was investigated using the phylogenetic and population genetics analyses of the mitochondrial DNA cytochrome-c-oxidase subunit I (COI) and the first variable domain of the control region (VD1), which were sequenced for >120 female individuals. These results showed a heterogeneous, panmictic population due to unrestricted gene flow. This can be attributed to extensive larval dispersal linked to the tidal circulation system in the back barrier basins of the Wadden Sea.     

Population genetic structure of mussels from the Baltic Sea

In a macrogeographic survey, the population genetic structure of mussels from various regions of the Baltic Sea, a large semi-enclosed brackish-water basin, was examined with reference toMytilus edulis andM. galloprovincialis samples from the North Sea, Irish coast and southern Portugal. Electrophoretically detectable variation was analysed at 6 polymorphic enzyme loci (Ap, Est-D, Lap-2, Odh, Pgi andPgm). Evidence was provided of a remarkably large amount of biochemical genetic differentiation among ecologically and morphologically divergent mussel populations in the Baltic. Patterns of allele frequencies in low-salinity populations from the area of the Baltic Proper were demonstrated to be widely homogeneous but contrast strongly with those of the western Baltic, the latter resembling populations from marine habitats of the North Sea. Associated with a pronounced salinity gradient, the spatial heterogeneity in gene-pool structure is indicated by steep clines of allele frequency changes in the area of the eastern Danish isles. The adaptive significance of the observed allozymic variation is suggested. From genetic distance estimates, the subdivision of population structure is discussed in relation to the significant amount of differentiation detected withinMytilus populations to date and to the evolutionary time required for the divergence of Baltic mussel populations. The allozymic data provide evidence for the genetic distinctiveness of mussels from the low-salinity areas of the Baltic. Their position at the specific or subspecific level of classification requires further consideration.     

Three species of Mytilus and their hybrids identified in a Scottish Loch: natives, relicts and invaders?

Three species of the mussel, Mytilus, occur in the North Atlantic region, M. edulis, M. galloprovincialis and M. trossulus, and hybrid zones are present where their distributions overlap. M. edulis is a native species in the UK. M. galloprovincialis originated in the Mediterranean and its distribution extends northwards along the Atlantic seaboard to Scotland. Baltic Sea mussels have a M. trossulus ancestry but are highly introgressed by M. edulis. In recent decades, farming of mussels on long-line rope culture systems has been introduced into Scotland. On farms in Loch Etive, a form of mussel with a fragile shell and a different shape to either M. edulis or M. galloprovincialis has been increasing in frequency over recent years. Samples of fragile shelled, normal strong shelled and intermediate mussel types were sampled from two farms in 2006 and compared with samples of M. edulis, M. galloprovincialis and M. trossulus from other sources where their species identity is well established. Abundance relative to depth, shell strength, condition index and shell morphology were analysed together with 5 allozyme loci and one nuclear DNA genetic marker (Me 15/16). The fragile shelled mussels, and many of those classed as intermediate, were identified as a mixture of M. trossulus and M. trossulus x M. edulis hybrids. This identification was strongly supported by both morphological and genetic data and is the first record of the presence of M. trossulus in UK waters. M. trossulus in Loch Etive are most likely to be a post-glacial relict population restricted to the low salinity area of the Loch that has recently increased in abundance due to commercial mussel growing activity. In addition, individual mussels of all three species and their hybrids were detected amongst Loch Etive mussels. This is the first genetic demonstration of all three species and their hybrids occurring together in one location in the Atlantic region and provides a unique opportunity to study the processes of speciation, divergence, and introgression in the genus Mytilus.     

Genetic Divergence Detected by ISSR Markers and Characterization of Microsatellite Regions in Mytilus Mussels

The wide distribution of microsatellites in mussels of the Mytilus edulis complex (Mytilidae) enables the analysis of inter-simple-sequence repeat (ISSR) markers. The aim of this investigation was to assess genetic differentiation in six sampling localities distributed along the European Atlantic coast to expose the potential of these markers in genetic studies requiring the detection of low polymorphism and as a source of sequences for developing microsatellite markers. We detected low genetic structuring within each member of the Mytilus edulis complex. Nei and Li distances and AMOVA clustered the individuals studied into two groups. On the basis of these results two sampling localities coming from the M. edulis × M. galloprovincialis mosaic hybrid zone in Western Europe were assigned to one species. On the other hand, mussels of a sampling locality in the Baltic Sea were not significantly different from a pure M. edulis locality supporting an extensive introgression of M. edulis in these individuals. Finally, 148 microsatellites were found in the sequences of 51 ISSR markers, and two polymorphic microsatellite markers were developed.     

Glacial history of the European marine mussels Mytilus,inferred from distribution of mitochondrial DNA lineages

Mussels of the genus Mytilus have been used to assess the circumglacial phylogeography of the intertidal zone. These mussels are representative components of the intertidal zone and have rapidly evolving mitochondrial DNA, suitable for high resolution phylogeographic analyses. In Europe, the three Mytilus species currently share mitochondrial haplotypes, owing to the cases of extensive genetic introgression. Genetic diversity of Mytilus edulis, Mytilus trossulus and Mytilus galloprovincialis was studied using a 900-bp long part of the most variable fragment of the control region from one of their two mitochondrial genomes. To this end, 985 specimens were sampled along the European coasts, at sites ranging from the Black Sea to the White Sea. The relevant DNA fragments were amplified, sequenced and analyzed. Contrary to the earlier findings, our coalescence and nested cladistics results show that only a single M. edulis glacial refugium existed in the Atlantic. Despite that, the species survived the glaciation retaining much of its diversity. Unsurprisingly, M. galloprovincialis survived in the Mediterranean Sea. In a relatively short time period, around the climatic optimum at 10 ky ago, the species underwent rapid expansion coupled with population differentiation. Following the expansion, further contemporary gene flow between populations was limited.     

Phylogeographic evidence for the postglacial colonization of the North and Baltic Sea coasts from inland glacial refugia by Triglochin maritima L.

We investigated the geographical distribution of genetic variation in 67 individuals of Triglochin maritima from 38 localities across Europe using AFLP markers. Analysis of genetic variation resulted in the recognition of two major genetic groups. Apart from few geographical outliers, these are distributed (1) along the Atlantic coasts of Portugal, Spain and France and (2) in the North Sea area, the Baltic Sea area, at central European inland localities, the northern Adriatic Sea coast and the Mediterranean coast of southwest France. Considering possible range shifts of T. maritima in reaction to Quaternary climatic changes as deduced from the present-day northern temperature limit of the species, Quaternary changes of coastline in the North Sea area and the very recent origin of the Baltic Sea, we conclude that the coastal populations of T. maritima in the North Sea and Baltic Sea areas originated from inland populations.     

Geographical patterns of dominant bivalves and a polychaete in Europe: no metapopulations in the marine coastal zone?

The genetic diversity, differentiation and performance of some dominant invertebrates in the marine coastal zone of Europe are reviewed in order to discuss the use of the metapopulation concept in the marine coastal realm. A consistently high genetic diversity of the species studied (mussels of the Mytilus edulis complex, Baltic clams Macoma balthica and lugworms Arenicola marina), a low differentiation and an almost uniform ecophysiological performance (determined by growth, maximum length, level reserve constituents or stress resistance) all along the coast of Europe do not support the use of the metapopulation concept. Electronic Publication     

A first report on coexistence and hybridization of <Emphasis Type="Italic">Mytilus trossulus</Emphasis> and <Emphasis Type="Italic">M. edulis</Emphasis> mussels in Greenland

Semi-sessile Mytilus mussels are used as indicators of climate changes, but their geographic distribution is not sufficiently known in the Arctic. The aim of this study was to investigate the taxonomic status and genetic differentiation of Mytilus populations in a Northwest Greenlandic fjord at Maarmorilik, impacted by contaminations from a former mine. In this study, mussels were collected at three sites differing in exposure to environmental factors. A total of 54 polymorphic SNPs found in the Mytilus EST and DNA sequences analyzed were successfully applied to 256 individuals. The results provided the first evidence for the existence of M. trossulus in Greenland. The mussel from M. trossulus and M. edulis taxa are shown to coexist and hybridize in the fjord. The three studied sites were found to differ significantly in the distribution of taxa with a higher prevalence of M. trossulus in the inner fjord. The identified M. edulis × M. trossulus hybrids mostly had a hybrid index score of about 0.5, indicating a similar number of alleles characteristic for M. trossulus and M. edulis. There was a low number of backcrosses between ‘pure’ taxa and hybrids. This newly discovered hybrid zone between the two taxa is unique in comparison with the Canadian populations. As Mytilus mussels in Greenland hitherto have been regarded as the one taxon M. edulis, the results have importance for biogeography and future monitoring and environmental studies.     

Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?

While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°‐13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10–30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15–10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.     

Phylogeography of the common ragworm Hediste diversicolor (Polychaeta: Nereididae) reveals cryptic diversity and multiple colonization events across its distribution

Previous studies on the common ragworm Hediste diversicolor (Polychaeta: Nereididae) revealed a marked genetic fragmentation across its distribution and the occurrence of sibling taxa in the Baltic Sea. These results suggested that the phylogeographic patterns of H. diversicolor could reflect interactions between cryptic differentiation and multiple colonization events. This study aims to describe the large-scale genetic structuring of H. diversicolor and to trace the phylogeographic origins of the genetic types described in the Baltic Sea. Samples of H. diversicolor (2 <  n  < 28) were collected at 16 locations across the NE Atlantic coasts of Europe and Morocco and in the Mediterranean, Black and Caspian Seas and sequenced at two mitochondrial gene fragments (COI and cyt b , 345 and 290 bp, respectively). Bayesian analyses revealed deep phylogeographic splits yielding three main clades corresponding to populations (i) from the NE Atlantic coasts (from Germany to Morocco) and from part of the Western Mediterranean, (ii) from the Mediterranean Sea, and (iii) from the Black and Caspian Seas. These clades are further divided in well-supported subclades including populations from different regions of NE Atlantic and Mediterranean (i.e. Portugal/Morocco, Western Mediterranean, Adriatic Sea). The Baltic Sea comprises three sympatric lineages sharing a common evolutionary history with populations from NE Atlantic, Western Mediterranean and Black/Caspian Seas, respectively. Hence, the current patterns of genetic structuring of H. diversicolor appear as the result of allopatric isolation, multiple colonization events and possible adaptation to local environmental conditions.     

Recombination in Mitochondrial DNA of European Mussels <Emphasis Type="Italic">Mytilus</Emphasis>

Mitochondrial DNA was long believed to be purely clonal and free from recombination. Major phylogenetic studies still depend on such assumptions. The peculiar genetic system of marine mussels Mytilus in which two divergent mitochondrial genomes exist provides a unique opportunity to study mtDNA recombination. Previous reports showed the existence of a few haplotypes having very strong recombination signal in the control region of mtDNA. Those recombinant variants have been found in a Baltic Sea population of Mytilus trossulus as well as in Mytilus galloprovincialis from the Black Sea. In both cases the mosaic genomes switched their transmission route and have been inherited paternally. In the present study rearranged mtDNA genomes found in all three European Mytilus species are described. The structure of their control region is a result of intra- and intermolecular recombination between mitochondrial genomes. Together with the phylogenetic reconstruction and geographic distribution, this suggests that two interlineage recombination events have occurred in the control region of mtDNA of European mussels Mytilus. Contrary to earlier observations, some of the mosaic genomes do not show any gender bias, which has important implications regarding the transmission and evolution of blue mussel mitochondrial genomes.     

The site of settlement indicates commensalism between bluemussel and its epibiont

Summary The site of settlement of barnacles (Balanus improvisus) attached on shells of bluemussels (Mytilus edulis) was mapped from a sample of mussels collected in the Baltic Sea. Most barnacles had settled near the siphonal apertures of the mussel. An experiment was made to measure the disadvantages and advantages that living in close association brings to barnacles and mussels. The barnacles on shells of living mussels were shown to grow significantly faster than those on empty mussel shells. Presence of barnacles had no effects on growth of mussels. The two-species association under study was demonstrated to be a case of commensalism.     

Bioinvasion threatens the genetic integrity of native diversity and a natural hybrid zone: smooth‐shelled blue mussels (Mytilus spp.) in the Strait of Magellan

Smooth‐shelled blue mussels of the Mytilus edulis species complex are widely distributed bivalve molluscs whose introductions threaten native marine biodiversity (non‐indigenous species – NIS). The aim of the present study was to identify the species and hybrids of Mytilus present in the Magellan Region (southern Chile). Results indicate that three mussel species of the Mytilus edulis complex are found in the region – M. edulis, M. chilensis (or the Southern Hemisphere lineage of Mytilus galloprovincialis), and M. galloprovincialis of Northern Hemisphere origin. For the first time, alleles of the introduced M. trossulus are reported from the Southern Hemisphere. In the Strait of Magellan the native Pacific blue mussel, Mytilus chilensis and the native Atlantic blue mussel, Mytilus edulis, meet and mix at a natural hybrid zone (about 125 km in length). This is the first record of a natural Mytilus hybrid zone in the Southern Hemisphere and is also the first record of the co‐occurrence of genes from all four Mytilus species in any one region. These results contribute to the knowledge of the biodiversity and delimitation of mussel species in southern South America, and highlight how introduced species may threaten the genetic integrity of native species through hybridization and introgression.     

Influence of genotype and geography on shell shape and morphometric trait variation among North Atlantic blue mussel (Mytilus spp.) populations

The influence of geography and genotype on shell shape (outline) and trait (morphometric) variation among North Atlantic blue mussels and their hybrids has been examined. Shape differences among reference taxa (Mytilus trossulus, Mytilus edulis and Mytilus galloprovincialis) were consistent with an association between taxon‐specific genes and shape genes. Newfoundland M. edulis × M. trossulus populations and northern Quebec M. trossulus populations exhibited an uncoupling of taxon‐specific genes from shape genes, whereas Nova Scotia M. trossulus populations and SW England M. edulis × M. galloprovincialis populations exhibited an association between taxon‐specific genes and shape genes. We found no evidence of a geographic effect (NE versus NW Atlantic) for shape variation, indicating that the genotype effect is stronger than any geographic effect at macrogeographic scales. Pronounced differences were observed in trait variability consistent with an association between taxon‐specific genes and trait genes in European populations, and trait divergence of New York M. edulis from all European mussels. Trait variability in mussels from Newfoundland, Nova Scotia and northern Quebec indicated an uncoupling of taxon genes from trait genes, whereas trait variability in SW England M. edulis × M. galloprovincialis populations was consistent with background genotype, indicating a strong association between taxon genes and trait genes. A pronounced macrogeographic split (NE versus NW Atlantic) regardless of taxonomy was observed, indicating that geography exerts a greater influence than genotype on trait variation at the macrogeographic scale. This is consistent with pronounced within‐taxon genetic divergence, indicative of different selection regimes or more likely of different evolutionary histories of mussels on either side of the North Atlantic. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 96 , 875–897.     

Genetic biodiversity in the Baltic Sea: species-specific patterns challenge management

Information on spatial and temporal patterns of genetic diversity is a prerequisite to understanding the demography of populations, and is fundamental to successful management and conservation of species. In the sea, it has been observed that oceanographic and other physical forces can constitute barriers to gene flow that may result in similar population genetic structures in different species. Such similarities among species would greatly simplify management of genetic biodiversity. Here, we tested for shared genetic patterns in a complex marine area, the Baltic Sea. We assessed spatial patterns of intraspecific genetic diversity and differentiation in seven ecologically important species of the Baltic ecosystem—Atlantic herring (Clupea harengus), northern pike (Esox lucius), European whitefish (Coregonus lavaretus), three-spined stickleback (Gasterosteus aculeatus), nine-spined stickleback (Pungitius pungitius), blue mussel (Mytilus spp.), and bladderwrack (Fucus vesiculosus). We used nuclear genetic data of putatively neutral microsatellite and SNP loci from samples collected from seven regions throughout the Baltic Sea, and reference samples from North Atlantic areas. Overall, patterns of genetic diversity and differentiation among sampling regions were unique for each species, although all six species with Atlantic samples indicated strong resistence to Atlantic-Baltic gene-flow. Major genetic barriers were not shared among species within the Baltic Sea; most species show genetic heterogeneity, but significant isolation by distance was only detected in pike and whitefish. These species-specific patterns of genetic structure preclude generalizations and emphasize the need to undertake genetic surveys for species separately, and to design management plans taking into consideration the specific structures of each species.     

Life on the margin: genetic isolation and diversity loss in a peripheral marine ecosystem, the Baltic Sea

Marginal populations are often isolated and under extreme selection pressures resulting in anomalous genetics. Consequently, ecosystems that are geographically and ecologically marginal might have a large share of genetically atypical populations, in need of particular concern in management of these ecosystems. To test this prediction, we analysed genetic data from 29 species inhabiting the low saline Baltic Sea, a geographically and ecologically marginal ecosystem. On average Baltic populations had lost genetic diversity compared to Atlantic populations: a pattern unrelated to dispersal capacity, generation time of species and taxonomic group of organism, but strongly related to type of genetic marker (mitochondrial DNA loci had lost c. 50% diversity, and nuclear loci 10%). Analyses of genetic isolation by geographic distance revealed clinal patterns of differentiation between Baltic and Atlantic regions. For a majority of species, clines were sigmoid with a sharp slope around the Baltic Sea entrance, indicating impeded gene flows between Baltic and Atlantic populations. Some species showed signs of allele frequencies being perturbed at the edge of their distribution inside the Baltic Sea. Despite the short geological history of the Baltic Sea (8000 years), populations inhabiting the Baltic have evolved substantially different from Atlantic populations, probably as a consequence of isolation and bottlenecks, as well as selection on adaptive traits. In addition, the Baltic Sea also acts a refuge for unique evolutionary lineages. This marginal ecosystem is thus vulnerable but also exceedingly valuable, housing unique genes, genotypes and populations that constitute an important genetic resource for management and conservation.     

Effect of salinity on filtration rates of mussels Mytilus edulis with special emphasis on dwarfed mussels from the low-saline Central Baltic Sea

The effect of salinity on the filtration rate of blue mussels, Mytilus edulis, from the brackish Great Belt (Denmark) and the low-saline Central Baltic Sea, respectively, was studied. First, we measured the effect of long-term (weeks) constant ambient salinities between 5 and 30 psu on the filtration rate of M. edulis collected in the Great Belt where the mean salinity is 17 psu. At salinities between 10 and 30 psu, the filtration rates did not vary much, but at 5 psu the filtration rates were significantly lower. Next, we studied dwarfed M. edulis (<25 mm shell length) from Central Baltic Sea (Askö, Sweden) where the mean salinity is 6.5 psu. The maximum filtration rate (F, ml min^?1 ind.^?1) as a function of shell length (L, mm) and dry weight of soft parts (W, mg) were found to be: F = 0.003L ^2.71 and F = 0.478W ^0.92, respectively, and these results indicate that the filtration rates of dwarfed Baltic Sea mussels are comparable to filtration rates of Great Belt mussels of similar size exposed to salinities >10 psu. When Baltic Sea mussels acclimatized to 20 psu in the laboratory were exposed to 6.5 psu this caused a drastic reduction in the filtration rate, but after about 2 days the previous high filtration rate was regained at 6.5 psu, and further, a similar pattern was observed when the 6.5 psu exposed mussels were finally re-exposed to 20 psu. The observed lack of Great Belt mussels to completely adjust to 5 psu, in contrast to the ease of Baltic Sea mussels to adjust back and forth between 6.5 and 20 psu, is remarkable and may perhaps be explained by different genotypes of Great Belt and Baltic Sea mussels.     

Chara baltica Bruzelius 1824 and Chara intermedia A. Braun 1859—Distinct species or habitat specific modifications?

Species separation in charophytes is primarily based on vegetative morphology, but the taxonomic delineation at species level worldwide is inconsistent. Here we compare ecophysiological characteristics of selected Chara species with vegetative morphological traits and genetic data. Four populations of the Chara baltica-Chara intermedia species cluster from locations along a north-south gradient through Europe were investigated. Physiological differences indicate habitat-specific adaptations. The Baltic Sea population from Hiddensee was found to be separated from a freshwater river population close to Munich with respect to light and salinity acclimation capabilities. However, a population from a brackish Mediterranean pool near Montpellier and one from the continental brackish water lake “Salziger See” near Halle, Germany, show intermediate ecophysiological characteristics, suggesting a continuum. Genetic analyses using AFLP (amplified fragment length polymorphism) match the results of the physiological analyses. The individuals of the Mediterranean and the Salziger See neither cluster to the typical C. baltica from the Baltic Sea nor to the C. intermedia from the freshwater habitat, but instead take an intermediate position. These results are stable against the background of a larger AFLP dataset on charophytes. Morphological analysis revealed no distinct groups and we therefore conclude that the C. baltica-C. intermedia cluster forms a physiological, morphological and genetic continuum.     

Différenciation allozymique multilocus des populations de moule Mytilus galloprovincialis Lmk. des côtes marocainesMultilocus allozyme differentiation of mussel populations Mytilus galloprovincialis Lmk. from Moroccan coasts.

Along the Moroccan coasts, the systematic status of Mytilus populations have been for a long time uncertain and confused, due to the use of unreliable morphometric criteria. In the present study, allozyme markers reveal the exclusive existence of M. galloprovincialis on Mediterranean and Atlantic coasts. Nei’s genetic distances are low and reflect a high gene flow between Atlantic and Mediterranean populations. However, a significant multilocus discontinuity revealed by F-statistics separate southern Atlantic populations from Mediterranean and north Atlantic ones and could be explained by a gene flow breaking because of a larval dispersal decrease, due to a sea surface current direction change from Cap Ghir towards the Canaries archipelago, and probably by differential selection effects in these two geographic areas.     

Optimal foraging versus shared doom effects: interactive influence of mussel size and epibiosis on predator preference

In the western Baltic Sea, the highly competitive blue mussel Mytilus edulis tends to monopolize shallow water hard substrata. In many habitats, mussel dominance is mainly controlled by the generalist predator Carcinus maenas. These predator-prey interactions seem to be affected by mussel size (relative to crab size) and mussel epibionts.There is a clear relationship between prey size and predator size as suggested by the optimal foraging theory: Each crab size class preferentially preys on a certain mussel size class. Preferred prey size increases with crab size.Epibionts on Mytilus, however, influence this simple pattern of feeding preferences by crabs. When offered similarly sized mussels, crabs prefer Balanus-fouled mussels over clean mussels. There is, however, a hierarchy of factors: the influence of attractive epibiotic barnacles is weaker than the factor ‘mussel size’. Testing small mussels against large mussels, presence or absence of epibiotic barnacles does not significantly alter preferences caused by mussel size. Balanus enhanced crab predation on mussels in two ways: Additional food gain and, probably more important, improvement in handling of the prey. The latter effect is illustrated by the fact that artificial barnacle mimics increased crab predation on mussels to the same extent as do live barnacles.We conclude that crab predation preferences follows the optimal foraging model when prey belong to different size classes, whereas within size classes crab preferences is controlled by epibionts.