INVITED REVIEW: Local adaptation and species segregation in two mussel (Mytilus edulis × Mytilus trossulus) hybrid zones

DNA adducts and polycyclic aromatic hydrocarbon (PAH) tissue levels in blue mussels (Mytilus spp.) from Nordic coastal sites

DNA adducts in gills and digestive gland, as well as polycyclic aromatic hydrocarbon (PAH) tissue levels were analysed in blue mussels (Mytilus spp.) from Nordic coastal areas (Iceland, Norway and Sweden) with diffuse or point sources of PAHs of various origins. Both DNA adduct and PAH tissue levels were generally low, indicating low PAH exposure to the mussels in the areas studied. DNA adducts were found to be higher in gills than in digestive gland of the mussels at all sites studied. Elevated DNA adduct levels in gills were found at 6 sites out of 18 compared to reference sites in respective coastal zones. Adduct levels ranged from 0.5 to 10 nmol adducts/mol normal nucleotides, being highest in mussels from Reykjavík harbour, Iceland (intertidal mussels), and from Fiskaatangen, Norway (subtidal mussels). Total PAH tissue levels in the mussels ranged between 40 and 11,670 ng/g dry wt., and were significantly correlated with DNA adduct levels (r(2)=0.73, p<0.001). PAH ratio values indicated that the PAHs were in most cases of pyrolytic origin, but with petrogenic input near harbours and an oil refinery.

Macoma balthicain the White and Barents Seas: properties of a widespread marine hybrid swarm (Mollusca: Bivalvia)

A main molecular subdivision in the circumpolar Macoma balthica complex has been described between Atlantic and Pacific taxa. In NE Europe, the clams of the White and Barents Seas, however, show deviant genetic structures. Using allozyme and mitochondrial DNA data, we explore the hypothesis that these deviations result from hybridization between an Atlantic (M. b. rubra) and an invading Pacific (M. b. balthica) lineage. A practically pure Atlantic Macoma extends from France north to the Varanger Peninsula (NE Norway), whereas populations farther east have genetic compositions intermediate between true Atlantic and true Pacific. Admixture estimates range from 32 to 90% Pacific contribution, with a notable deviation in a nearly pure Atlantic outpost in the Mezen Bay (NE White Sea). The pattern of variation is not one of a simple collinear mixing however. Different characters exhibit different degrees of introgression, and the relative introgression varies regionally. Yet, there are practically no interlocus genotypic disequilibria between the diverged loci, which brings out the White Sea-Barents Sea M. balthica as the best-documented marine animal hybrid swarms so far, arisen through amalgamation of genomes previously isolated since pre-Pleistocene times. On top of the main admixture pattern, strong geographical structuring is also seen in characters unrelated to the principal systematic distinction. The persistence of the regional patterns indicates restricted gene flow at the present time, despite the high dispersal potential of the species. The causes of this structuring could be in a complex history of colonization events and features of local hydrography enhancing isolation and divergence of populations.

Seasonal variability of free amino acids in two marine bivalves, Macoma balthica and Mytilus spp., in relation to environmental and physiological factors

The seasonal variability of the intracellular free amino acid (FAA) concentration was studied in 5 Macoma balthica populations and 7 Mytilus spp. populations along their European distribution. Because of the well known physiological role of FAA as organic osmolytes for salinity induced cell volume regulation in marine osmoconformers, FAA variations were compared in bivalve populations that were exposed to high vs. low intraannual salinity fluctuations. In general, seasonal FAA variations were more pronounced in M. balthica than in Mytilus spp. In both bivalve taxa from different locations in the Baltic Sea, highest FAA concentrations were found in autumn and winter and low FAA concentrations were measured in summer. Seasonal patterns were less pronounced in both taxa at locations with constant salinity conditions. In contrast to Baltic Sea populations, Atlantic and Mediterranean bivalves showed high FAA concentrations in summer and low values in winter, regardless of seasonal salinity fluctuations. Significant seasonal FAA variations at locations with constant salinity conditions showed that salinity appeared not to be the main factor in determining FAA concentrations. The seasonal patterns of the main FAA pool components, i.e. alanine, glycine and taurine, are discussed in the context of seasonal variations in environmental factors (salinity, temperature) and physiological state (glycogen content, reproductive stage).

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 x 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.

Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoan Celleporella hyalina

The marine bryozoan Celleporella hyalina is a species complex composed of many highly divergent and mostly allopatric genetic lineages that are reproductively isolated but share a remarkably similar morphology. One such lineage commonly encrusts macroalgae throughout the NE Atlantic coast. To explore the processes leading to geographical diversification, reproductive isolation and speciation in this taxon, we (i) investigated NE Atlantic C. hyalina mitochondrial DNA phylogeography, and (ii) used breeding trials between geographical isolates to ascertain reproductive isolation. We find that haplotype diversity is geographically variable and there is a strong population structure, with significant isolation by distance. NE Atlantic C. hyalina is structured into two main parapatric lineages that appear to have had independent Pleistocene histories. Range expansions have resulted in two contact zones in Spain and W Ireland. Lineage 1 is found from Ireland to Spain and has low haplotype diversity, with closely related haplotypes, suggesting a recent population expansion into the Irish Sea, S Ireland, S England and Spain. Lineage 2 is found from Iceland to Spain and has high haplotype diversity. Complete reproductive isolation was found between some geographical isolates representing both lineages, whereas it was incomplete or asymmetric between others, suggesting these latter phylogeographical groups probably represent incipient species. The phylogeographical distribution of NE Atlantic C. hyalina does not fall easily into a pattern of southern refugia, and we discuss likely differences between terrestrial and marine system responses to Pleistocene glacial cycles.

Development of mussel mRNA profiling: Can gene expression trends reveal coastal water pollution?

Marine bivalves of the genus Mytilus are intertidal filter-feeders commonly used as biosensors of coastal pollution. Mussels adjust their functions to ordinary environmental changes, e.g. temperature fluctuations and emersion-related hypoxia, and react to various contaminants, accumulated from the surrounding water and defining a potential health risk for sea-food consumers. Despite the increasing use of mussels in environmental monitoring, their genome and gene functions are largely unexplored. Hence, we started the systematic identification of expressed sequence tags and prepared a cDNA microarray of Mytilus galloprovincialis including 1714 mussel probes (76% singletons, approximately 50% putatively identified transcripts) plus unrelated controls. To assess the potential use of the gene set represented in MytArray 1.0, we tested different tissues and groups of mussels. The resulting data highlighted the transcriptional specificity of the mussel tissues. Further testing of the most responsive digestive gland allowed correct classification of mussels treated with mixtures of heavy metals or organic contaminants (expression changes of specific genes discriminated the two pollutant cocktails). Similar analyses made a distinction possible between mussels living in the Venice lagoon (Italy) at the petrochemical district and mussels close to the open sea. The suggestive presence of gene markers tracing organic contaminants more than heavy metals in mussels from the industrial district is consistent with reported trends of chemical contamination. Further study is necessary in order to understand how much gene expression profiles can disclose the signatures of pollutants in mussel cells and tissues. Nevertheless, the gene expression patterns described in this paper support a wider characterization of the mussel transcriptome and point to the development of novel environmental metrics.

Doubly uniparental inheritance is associated with high polymorphism for rearranged and recombinant control region haplotypes in Baltic Mytilus trossulus

Many bivalve species, including mussels of the genus Mytilus, are unusual in having two mtDNA genomes, one inherited maternally (the F genome) and the other inherited paternally (the M genome). The sequence differences between the genomes are usually great, indicating ancient divergence predating speciation events. However, in Mytilus trossulus from the Baltic, both genomes are similar to the F genome from the closely related M. edulis. This study analyzed the mtDNA control region structure in male and female Baltic M. trossulus mussels. We show that a great diversity of structural rearrangements is present in both sexes. Sperm samples are dominated by recombinant haplotypes with M. edulis M-like control region segments, some having large duplications. By contrast, the rearranged haplotypes that dominate in eggs lack segments from this M genome. The rearrangements can be explained by a combination of tandem duplication, deletion, and intermolecular recombination. An evolutionary pathway leading to the recombinant haplotypes is suggested. The data are also considered in relation to the hypothesis that the M. edulis M-like control region sequence is necessary to confer the paternal role on genomes that are otherwise F-like.

A genomic scan for divergent selection in a secondary contact zone between Atlantic and Gulf of Mexico oysters,Crassostrea virginica

The degree of population structure within species often varies considerably among loci. This makes it difficult to determine whether observed variance reflects neutral-drift stochasticity or locus-specific selection at one or more loci. This uncertainty is exacerbated when evolutionary equilibrium cannot be assumed and/or admixture potentially inflates genomic variance. Thus, the value of a 'genome scan', where locus-specific summary statistics are compared with a simulated neutral distribution among loci, may be limited in secondary contact zones if the null distribution is sensitive to the timing of secondary contact. Of particular interest here is the wide variance previously observed in locus-specific patterns of population structure between Atlantic and Gulf of Mexico populations of eastern oyster, Crassostrea virginica. To test the robustness of an equilibrium null model, we compared among-locus distributions of F(ST) simulated under migration-drift equilibrium and several nonequilibrium secondary contact histories. We then tested for evidence of divergent selection between two oyster populations on either side of a secondary contact zone using 215 amplified fragment length polymorphism (AFLP) loci. Constant-migration equilibrium and nonequilibrium secondary-contact simulations produced equivalent distributions of F(ST) when anchored by the global mean F(ST) observed in oysters, 0.0917. The 99th quantile of simulated neutral F(ST) encompassed most of the variation among oyster loci. Three AFLP loci exhibited F(ST) values higher than this threshold. Although no locus was significant after correcting for multiple tests, our results show in geographically clinal organisms: AFLPs can efficiently characterize the genomic distribution of F(ST); equilibrium models can be used to evaluate outliers; these procedures help focus research on genomic regions of interest.

Geographic Variation and Positive Selection on M7 Lysin, an Acrosomal Sperm Protein in Mussels (Mytilus spp.)

Successful fertilization in free-spawning marine organisms depends on the interactions between genes expressed on the surfaces of eggs and sperm. Positive selection frequently characterizes the molecular evolution of such genes, raising the possibility that some common deterministic process drives the evolution of gamete recognition genes and may even be important for understanding the evolution of prezygotic isolation and speciation in the marine realm. One hypothesis is that gamete recognition genes are subject to selection for prezygotic isolation, namely, reinforcement. In a previous study, positive selection on the gene coding for the acrosomal sperm protein M7 lysin was demonstrated among allopatric populations of mussels in the Mytilus edulis species group (M. edulis, Mytilus galloprovincialis, and Mytilus trossulus). Here, we expand sampling to include M7 lysin haplotypes from populations where mussel species are sympatric and hybridize to determine whether there is a pattern of reproductive character displacement (RCD), which would be consistent with reinforcement driving selection on this gene. We do not detect a strong pattern of RCD; neither are there unique haplotypes in sympatry nor is there consistently greater population structure in comparisons involving sympatric populations. One distinct group of haplotypes, however, is strongly affected by natural selection, and this group of haplotypes is found within M. galloprovincialis populations throughout the Northern Hemisphere concurrent with haplotypes common to M. galloprovincialis and M. edulis. We suggest that balancing selection, perhaps resulting from sexual conflicts between sperm and eggs, maintains old allelic diversity within M. galloprovincialis.

INVITED REVIEW: 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.

Genetic Clines in the Bay of Biscay Provide Estimates of Migration for Sardina pilchardus

Nine allozymic loci in 1,635 individuals of Sardina pilchardus obtained at 33 sites ranging from the North to the South limits of the Bay of Biscay were analyzed to provide a description of the genetic structure of the sardine population. Individual body size and age were also recorded and analyzed. In the study population, weak but significant genetic differences were found, and a cline was observed between multilocus heterozygosity and longitude. The cline was predominantly driven by allelic frequencies of two loci, PGM-1* and PEP-lt*, and using a cline model, we estimated a migration rate of 103.1 km/gen (dispersal distance per generation). In addition, we observed that the cline was linked to biological data such as mean length and mean age of the fish. Two hypotheses may explain this cline: mixing of two different populations in the Bay of Biscay or a selective process. The weak genetic differences, the important dispersal distance per generation, and the link between genetic and biological data suggest that selection is likely to be the primary factor that maintains the cline.

Genetic structure and environmental heterogeneity in the European hake (Merluccius merluccius)

This study aimed at assessing the genetic structure and the state of the stocks of the European hake (Merluccius merluccius). To this end, 15 samples were taken from the whole range of the species and analysed using allozymes. Since 11 samples were taken from the poorly studied Mediterranean Sea, the results obtained provided a complete picture of the hake's genetic structure and an initial insight into its relationships with environmental features. Atlantic and Mediterranean hake populations are separated by the Almeria-Oran front. This area has been proved to be the boundary between Atlantic and Mediterranean stocks of many marine organisms, but some doubt exists concerning the efficaciousness of the local gyres as barriers to the gene flow. Our data have evidenced a latitudinal cline at loci Gapdh and Gpi-2 within the Mediterranean Sea, with a further steep change across the Almeria-Oran front. The genetic pattern showed a strong correlation with the values of the salinity both at the surface and at -320 m and of the salinity + temperature at the surface, suggesting a role for these parameters in maintaining the genetic differentiation among the two population groups through selective processes. Finally, the levels of genetic variability were found to be slightly lower in the depleted Atlantic stock than in the Mediterranean one.