A broad transition zone between an inner Baltic hybrid swarm and a pure North Sea subspecies of Macoma balthica (Mollusca, Bivalvia)

The populations of the bivalve clam Macoma balthica in the low-salinity Northern Baltic Sea represent an admixture of two strongly diverged genomic origins, the Pacific Macoma balthica balthica (approx. 60% genomic contribution) and Atlantic Macoma balthica rubra (40%). Using allozyme and mtDNA characters, we describe the broad transition from this hybrid swarm to the pure M. b. rubra in the saline North Sea waters, spanning hundreds of kilometre distance. The zone is centred in the strong salinity gradient of the narrow Öresund strait and in the adjacent Western Baltic. Yet the multilocus clines show no simple and smoothly monotonic gradation: they involve local reversals and strong differences between neighbouring populations. The transitions in different characters are not strictly coincident, and the extent of introgression varies among loci. The Atlantic influence extends further into the Baltic in samples from the southern and eastern Baltic coasts than on the western coast, and further in deeper bottoms than at shallow (< 1 m) sites. This fits with the counterclockwise net circulation pattern and with a presumably weaker salinity barrier for invading Atlantic type larvae in saline deeper water, and corresponding facilitation of outwards drift of Baltic larvae in diluted surface waters. Genotypic disequilibria were strong particularly in the shallow-water samples of the steepest transition zone. This suggests larval mixing from different sources and limited interbreeding in that area, which makes a stark contrast to the evidence of thorough amalgamation of the distinct genomic origins in the inner Baltic hybrid swarm of equilibrium structure.     

Genetic, morphological, and feather isotope variation of migratory willow warblers show gradual divergence in a ring

The circular distribution of the willow warbler Phylloscopus trochilus around the Baltic Sea shares many features with the classic examples of ring species; however, the system is much younger. It has previously been shown that a secondary contact zone is located in central Scandinavia, where there are narrow clines for several morphological traits coincident with a migratory divide. Here we analyse multiple traits and genes from > 1700 males captured on breeding territories at 77 sites spread around the Baltic Sea to test the following hypothesis. If the secondary contact zone in Scandinavia is a result of divergence in two allopatric refuge populations during the last glaciation, we expect to find a similar secondary contact zone somewhere else around the circular distribution. Our results show that the trait clines were wider and displaced from each other along the eastern side of the Baltic Sea. Analyses of 12 microsatellite loci confirmed that the genome is very similar between the terminal forms ( F _ST = 0). Two AFLP-derived markers filtered out from a genomic scan instead appear to be maintained by selection. These markers exhibited steep clines at the secondary contact zone in Scandinavia, but as for the phenotypic traits, had vastly different cline centres east of the Baltic Sea. The trait clines along the ring distribution outside the Scandinavian secondary contact zone thus seem to have been shaped by independent action of selection or drift during the process of postglacial colonization.     

Diversity and trans-arctic invasion history of mitochondrial lineages in the North Atlantic Macoma balthica complex (Bivalvia: Tellinidae)

The history of repeated inter- or transoceanic invasions in bivalve mollusks of the circumpolar Macoma balthica complex was assessed from mtDNA COIII sequences. The data suggest that four independent trans-Arctic invasions, from the Pacific, gave rise to the current lineage diversity in the North Atlantic. Unlike in many other prominent North Atlantic littoral taxa, no evidence for (postinvasion) trans-Atlantic connections was found in the M. balthica complex. The earliest branch of the mtDNA tree is represented by the temperate-boreal North American populations (=Macoma petalum), separated from the M. balthica complex proper in the Early Pliocene at latest. The ensuing trans-Arctic invasions established the North European M. b. rubra, which now prevails on the North Sea and northeast Atlantic coasts, about two million years ago, and the currently northwest Atlantic M. balthica lineage in the Canadian Maritimes, in the Middle Pleistocene. The final reinvasion(s) introduced a lineage that now prevails in a number of North European marginal seas and is still hardly distinguishable from North Pacific mtDNA (M. b. balthica). We used coalescence simulation analyses to assess the age of the latest invasion from the Pacific to the northeast Atlantic. The results refute the hypothesis of recent, human-mediated reintroductions between northeast Pacific and the North European marginal seas in historical times. Yet they also poorly fit the alternative hypotheses of an early postglacial trans-Arctic invasion (< 11 thousand years ago), or an invasion during the previous Eemian interglacial (120 thousand years ago). Divergence time estimates rather fall in the Middle Weichselian before the Last Glacial Maximum, in conflict with the conventional thinking of trans-Arctic biogeographical connections; an early Holocene reinvasion may still be regarded as the most plausible scenario. Today, the most recently invaded Pacific mtDNA lineage is found admixed with the earlier established European Atlantic "rubra" lineage in the Baltic Sea and in Barents Sea populations east of the Varanger peninsula, and it is practically exclusive in the White and Pechora seas. Yet mtDNA does not always constitute an unequivocal taxonomic marker at individual level; the marginal populations represent hybrid swarms of the Atlantic and Pacific lineages in their nuclear genes.     

Macoma balthica in 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.     

Comparing clines on molecular and phenotypic traits in hybrid zones: a window on tension zone models

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.     

20-isosterols from the marine bivalve Macoma balthica

A fraction of 20-isosterols has been isolated from the Baltic Sea bivalve Macoma balthica. The main component of this fraction, (20S)-cholest-5-en-3 beta-ol, has been characterized by 300 MHz nuclear magnetic resonance spectroscopy. 20-Isosterols in M balthica probably originate from sea-bottom sediments, an important component of a diet of the animal studied.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.     

Molecular and morphological analysis of secondary contact zones of Cottus gobio in Fennoscandia: geographical discordance of character transitions

Northern Europe was postglacially colonized from different directions by distinct phylogeographical lineages of the bullhead Cottus gobio L. (Pisces: Scorpaeniformes). These lineages have then come into contact in coastal habitats of the currently brackish Baltic Sea and in the freshwaters north of it. We studied the patterns of intergradation in the contact zones in four morphometric and six molecular characters. In the north, intergradation between the western (W) and eastern (E) bullhead lineages is found both among rivers (west‐to‐east) and along individual rivers (south‐to‐north). The locations of the transition zones probably relate to the timing of the initial contact, subsequent Baltic shoreline displacement (i.e. emergence of the lower river reaches), and dispersal barriers caused by variations of coastal salinity. The transitions (clines) in different characters are, however, not geographically coincident. Mitochondrial DNA clines are generally found upstream and to the east of the other transitions, and GPI‐1 allozyme clines are mostly shifted downstream in the rivers, and west of the other transitions on the broader scale of the Baltic Sea. The location of the mtDNA clines may best reflect the initial contact between lineages, and the displacement of the other clines could result from dispersal being overall asymmetric (predominantly downstream) and sex‐biased (stronger in males). Alternatively, the non‐coincidence might reflect selection against deleterious cytonuclear character combinations. No clear evidence of reproductive incompatibility between the lineages was seen in local population structures; no remaining genetic correlations were observed locally among traits. In another transition area, a coastal transect in southern Finland, clinal patterns similar to those in the northern contact zone were recorded, but the population compositions could not be explained by simple in situ mixing of any of the putatively pure, invading refugial lineages. Probably, the bullhead stocks that initially came into contact in this southern study area already represented mixtures of the invading lineages. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 81 , 535–552.     

Evidence of a hybrid-zone in Atlantic cod (Gadus morhua) in the Baltic and the Danish Belt Sea revealed by individual admixture analysis

The study of hybrid zones is central to our understanding of the genetic basis of reproductive isolation and speciation, yet very little is known about the extent and significance of hybrid zones in marine fishes. We examined the population structure of cod in the transition area between the North Sea and the Baltic Sea employing nine microsatellite loci. Genetic differentiation between the North Sea sample and the rest increased along a transect to the Baltic proper, with a large increase in level of differentiation occurring in the Western Baltic area. Our objective was to determine whether this pattern was caused purely by varying degrees of mechanical mixing of North Sea and Baltic Sea cod or by interbreeding and formation of a hybrid swarm. Simulation studies revealed that traditional Hardy-Weinberg analysis did not have sufficient power for detection of a Wahlund effect. However, using a model-based clustering method for individual admixture analysis, we were able to demonstrate the existence of intermediate genotypes in all samples from the transition area. Accordingly, our data were explained best by a model of a hybrid swarm flanked by pure nonadmixed populations in the North Sea and the Baltic Sea proper. Significant correlation of gene identities across loci (gametic phase disequilibrium) was found only in a sample from the Western Baltic, suggesting this area as the centre of the apparent hybrid zone. A hybrid zone for cod in the ecotone between the high-saline North Sea and the low-saline Baltic Sea is discussed in relation to its possible origin and maintenance, and in relation to a classical study of haemoglobin variation in cod from the Baltic Sea/Danish Belt Sea, suggesting mixing of two divergent populations without interbreeding.     

Outlier Loci Detect Intraspecific Biodiversity amongst Spring and Autumn Spawning Herring across Local Scales

Herring, Clupea harengus, is one of the ecologically and commercially most important species in European northern seas, where two distinct ecotypes have been described based on spawning time; spring and autumn. To date, it is unknown if these spring and autumn spawning herring constitute genetically distinct units. We assessed levels of genetic divergence between spring and autumn spawning herring in the Baltic Sea using two types of DNA markers, microsatellites and Single Nucleotide Polymorphisms, and compared the results with data for autumn spawning North Sea herring. Temporally replicated analyses reveal clear genetic differences between ecotypes and hence support reproductive isolation. Loci showing non-neutral behaviour, so-called outlier loci, show convergence between autumn spawning herring from demographically disjoint populations, potentially reflecting selective processes associated with autumn spawning ecotypes. The abundance and exploitation of the two ecotypes have varied strongly over space and time in the Baltic Sea, where autumn spawners have faced strong depression for decades. The results therefore have practical implications by highlighting the need for specific management of these co-occurring ecotypes to meet requirements for sustainable exploitation and ensure optimal livelihood for coastal communities.     

Disjunct distribution of highly diverged mitochondrial lineage clade and population subdivision in a marine bivalve with pelagic larval dispersal

Mitochondrial DNA sequence data for 295 individuals of the marine bivalve Macoma balthica (L.) were collected from 10 sites across the European distribution, and from Alaska. The data were used to infer population subdivision history and estimate current levels of gene flow. Inferred historical biogeography was expected to be congruent with colonization of the Atlantic Ocean from the Pacific Ocean after the opening of the Bering Strait 3.5 Ma. In addition, the last glacial maximum, about 18000 years ago, was expected to have been responsible for most of the present-day distribution of molecular variation within Europe, because the area must have been recolonized after confinement to France and the south of the British Isles during the last glacial maximum. Current gene flow was hypothesized to be high, because the larvae of M. balthica spend 2-5 weeks drifting in the water column. The geographical distribution of one highly diverged haplotype clade was found to be disjunct and was encountered exclusively in samples from the Baltic Sea and Alaska. A molecular clock calibration for marine bivalve cytochrome-c-oxidase I dates this clade as having split off from the other haplotypes 9.8-39 Ma. Multiple colonizations of the Atlantic Ocean from the Pacific by M. balthica may explain the strong differences found between Baltic Sea and other European populations of this species. The sympatric occurrence of the highly diverged mitochondrial lineages in western parts of the Baltic Sea points to secondary admixture. With the use of coalescent analysis, population divergence times for French vs. other non-Baltic European populations ('Atlantic population assemblage') were estimated at a minimum of about 110000 years ago, well before the last glacial maximum 18000 years ago. Signatures of population divergence of M. balthica that appear to have originated during the Pleistocene have thus survived the last glacial maximum. Some of the populations within the Atlantic assemblage are currently isolated, while others appear to be connected by gene flow. Apparently, populations of this species can remain highly subdivided in spite of the potential for high gene flow, implying that their population and evolutionary dynamics can be independent.     

Benthic predator-prey interactions: evidence that adult Monoporeia affinis (Amphipoda) eat postlarval Macoma balthica (Bivalvia)

Predation by adults of the amphipod Monoporeia affinis on the plantigrade postlarval stage of the bivalve Macoma balthica was studied in the laboratory. We confirmed that M. affinis consumes small M. balthica. Amphipods offered 14C-labelled postlarvae took up the radioactive tracer, while those presented Rhodamine B-stained postlarvae acquired gut contents fluorescing strongly in orange, whereas control amphipods did not. Both labelling methods proved convenient to use in laboratory experiments, and are particularly useful when organisms lack structures that can be easily identified after being ingested, or when cross-over reactions may bias the results of immunoassays. The results reported here support the conclusion from earlier studies that predation by M. affinis on M. balthica can affect population dynamics of M. balthica and is likely to be an important structuring factor in the low-diversity benthic macrofauna community of the Baltic Sea.     

Seasonal variability in biomarkers in the bivalves Mytilus edulis and Macoma balthica from the northern Baltic Sea

Metallothionein level (MT), and acetylcholinesterase (AChE), catalase (CAT) and glutathione-S-transferase (GST) enzyme activities in the bivalves Mytilus edulis and Macoma balthica were investigated for seasonal variations from an inshore and an offshore site in the northern Baltic Sea. All the biomarkers showed variability, following mostly a similar pattern at both sites. Relationships between biomarkers and environmental factors and protein concentration and weight of target tissues were examined. In M. edulis, GST activity was related to Secchi depth, while in M. balthica a correlation with near-bottom oxygen saturation was observed. AChE activity correlated with the weight of the foot tissue of M. balthica. In both species, an integrated biomarker index indicated a stressed condition during the spring/early summer period. Strong seasonal variability in temperature and a concentrated period of food availability in spring-both governing the reproductive cycle of the bivalves-probably explains most of the observed natural variability in biomarkers in this sea area.     

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

BIOCHEMICAL PHENOTYPIC AND GENETIC STUDIES OF TWO INTRODUCED FIRE ANTS AND THEIR HYBRID (HYMENOPTERA: FORMICIDAE)

Two introduced fire ants, Solenopsis invicta and S. richteri, and their hybrid were studied using phenotypic markers from gas chromatographic analysis of hydrocarbons and venom alkaloids, as well as genetic markers from enzyme electrophoresis. Both methods show that extensive gene introgression is occurring over a distance of at least 120 km at the contact zone between the two forms in eastern Mississippi. Genetic analyses suggest that the hybrid population does not depart dramatically from panmixia. Also, recombinant genotypes predominate in the hybrid zone, indicating that F₁ hybrids are viable. Allele frequency clines through the hybrid zone are apparent for four polymorphic loci. Data sets generated by the chromatographic and electrophoretic methods are highly concordant in that they differentiate completely between the two forms and agree in designating colonies from the contact zone as hybrid or parental in a high proportion (90%) of cases. The two methods can serve as complementary tools for studying closely related but genetically distinct populations in this, and perhaps other, groups of insects.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

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