Phylogeography of amphi-boreal fish: tracing the history of the Pacific herring Clupea pallasii in North-East European seas

Genetic structure of marine and lake forms of Pacific herring Clupea pallasii

The Pacific herring (Clupea pallasii) is one of the most important species in the commercial fisheries distributed in the North Pacific Ocean and the northeastern European seas. This teleost has marine and lake ecological forms a long its distribution in the Holarctic. However, the level of genetic differentiation between these two forms is not well known. In the present study, we used ddRAD-sequencing to genotype 54 specimens from twelve wild Pacific herring populations from the Kara Sea and the Russian part of the northwestern Pacific Ocean for unveiling the genetic structure of Pacific herring. We found that the Kara Sea population is significantly distinct from Pacific Ocean populations. It was demonstrated that lake populations of Pacific herring differ from one another as well as from marine specimens. Our results show that fresh and brackish water Pacific herring, which inhabit lakes, can be distinguished as a separate lake ecological form. Moreover, we demonstrate that each observed lake Pacific herring population has its own and unique genetic legacy.

Cycles of trans-Arctic dispersal and vicariance, and diversification of the amphi-boreal marine fauna

The amphi-boreal faunal element comprises closely related species and conspecific populations with vicarious distributions in the North Atlantic and North Pacific basins. It originated from an initial trans-Arctic dispersal in the Pliocene after the first opening of the Bering Strait, and subsequent inter-oceanic vicariance through the Pleistocene when the passage through the Arctic was severed by glaciations and low sea levels. Opportunities for further trans-Arctic dispersal have risen at times, however, and molecular data now expose more complex patterns of inter-oceanic affinities and dispersal histories. For a general view on the trans-Arctic dynamics and of the roles of potential dispersal-vicariance cycles in generating systematic diversity, we produced new phylogeographic data sets for amphi-boreal taxa in 21 genera of invertebrates and vertebrates, and combined them with similar published data sets of mitochondrial coding gene variation, adding up to 89 inter-oceanic comparisons involving molluscs, crustaceans, echinoderms, polychaetes, fishes and mammals. Only 39% of the cases correspond to a simple history of Pliocene divergence; in most taxonomical groups, the range of divergence estimates implies connections through the entire Pliocene-Pleistocene-Holocene time frame. Repeated inter-oceanic exchange was inferred for 23 taxa, and the latest connection was usually post-glacial. Such repeated invasions have usually led to secondary contacts and occasionally to widespread hybridization between the different invasion waves. Late- or post-glacial exchange was inferred in 46% of the taxa, stressing the importance of the relatively recent invasions to the current diversity in the North Atlantic. Individual taxa also showed complex idiosyncratic patterns and histories, and several instances of cryptic speciation were recognized. In contrast to a simple inter-oceanic vicariance scenario underlying amphi-boreal speciation, the data expose complex patterns of reinvasion and reticulation that complicate the interpretation of taxon boundaries in the region.

Genetic population structure of the pelagic mollusk Limacina helicina in the Kara Sea

Background

Pelagic pteropods Limacina helicina are widespread and can play an important role in the food webs and in biosedimentation in Arctic and Subarctic ecosystems. Previous publications have shown differences in the genetic structure of populations of L. helicina from populations found in the Pacific Ocean and Svalbard area. Currently, there are no data on the genetic structure of L. helicina populations in the seas of the Siberian Arctic. We assessed the genetic structure of L. helicina from the Kara Sea populations and compared them with samples from around Svalbard and the North Pacific.

Methods

We examined genetic differences in L. helicina from three different locations in the Kara Sea via analysis of a fragment of the mitochondrial gene COI. We also compared a subset of samples with L. helicina from previous studies to find connections between populations from the Atlantic and Pacific Oceans.

Results

65 individual L. helinica from the Kara Sea were sequenced to produce 19 different haplotypes. This is comparable with numbers of haplotypes found in Svalbard and Pacific samples (24 and 25, respectively). Haplotypes from different locations sampled around the Arctic and Subarctic were combined into two different groups: H1 and H2. The H2 includes sequences from the Kara Sea and Svalbard, was present only in the Atlantic sector of the Arctic. The other genetic group, H1, is widespread and found throughout all L. helicina populations. ϕ ST analyses also indicated significant genetic difference between the Atlantic and Pacific regions, but no differences between Svalbard and the Kara Sea.

Discussion

The obtained results support our hypothesis about genetic similarity of L. helicina populations from the Kara Sea and Svalbard: the majority of haplotypes belongs to the haplotype group H2, with the H1 group representing a minority of the haplotypes present. In contrast, in the Canadian Arctic and the Pacific Ocean only haplogroup H1 is found. The negative values of Fu’s Fs indicate directed selection or expansion of the population. The reason for this pattern could be an isolation of the Limacina helicina population during the Pleistocene glaciation and a subsequent rapid expansion of this species after the last glacial maximum.

The evolutionary history of the goby Elacatinus puncticulatus in the tropical eastern pacific: Effects of habitat discontinuities and local environmental variability

Habitat discontinuities, temperature gradients, upwelling systems, and ocean currents, gyres and fronts, can affect distributions of species with narrow environmental tolerance or motility and influence the dispersal of pelagic larvae, with effects ranging from the isolation of adjacent populations to connections between them. The coast of the Tropical Eastern Pacific (TEP) is a highly dynamic environment, with various large gyres and upwelling systems, alternating currents and large rocky-habitat discontinuities, which may greatly influence the genetic connectivity of populations in different parts of the coast. Elacatinus puncticulatus is a cryptic, shallow-living goby that is distributed along the continental shore of virtually the entire TEP, which makes it a good model for testing the influence of these environmental characteristics in the molecular evolution of widespread species in this region. A multilocus phylogeny was used to evaluate the influence of habitat gaps, and oceanographic processes in the evolutionary history of E. puncticulatus throughout its geographical range in the TEP. Two well-supported allopatric clades (one with two allopatric subclades) were recovered, the geographic distribution of which does not correspond to any previously proposed major biogeographic provinces. These populations show strong genetic structure and substantial genetic distances between clades and sub-clades (cytb 0.8-7.3%), with divergence times between them ranging from 0.53 to 4.88 Mya, and recent population expansions dated at 170-130 Kya. The ancestral area of all populations appears to be the Gulf of Panama, while several isolation events have formed the phylogeographic patterns evident in this species. Local and regional oceanographic processes as well as habitat discontinuities have shaped the distribution patterns of the genetic lineages along the continental TEP. Large genetic distances, high genetic differentiation, and the results of species-tree and phylogenetic analyses indicate that E. puncticulatus comprises a complex of three allopatric species with an unusual geographic arrangement.

Glacial vicariance drives phylogeographic diversification in the amphi-boreal kelp Saccharina latissima

Glacial vicariance is regarded as one of the most prevalent drivers of phylogeographic structure and speciation among high-latitude organisms, but direct links between ice advances and range fragmentation have been more difficult to establish in marine than in terrestrial systems. Here we investigate the evolution of largely disjunct (and potentially reproductively isolated) phylogeographic lineages within the amphi-boreal kelp Saccharina latissima s. l. Using molecular data (COI, microsatellites) we confirm that S. latissima comprises also the NE Pacific S. cichorioides complex and is composed of divergent lineages with limited range overlap and genetic admixture. Only a few genetic hybrids were detected throughout a Canadian Arctic/NW Greenland contact zone. The degree of genetic differentiation and sympatric isolation of phylogroups suggest that S. latissima s. l. represents a complex of incipient species. Phylogroup distributions compared with paleo-environmental reconstructions of the cryosphere further suggest that diversification within S. latissima results from chronic glacial isolation in disjunct persistence areas intercalated with ephemeral interglacial poleward expansions and admixture at high-latitude (Arctic) contact zones. This study thus supports a role for glaciations not just in redistributing pre-existing marine lineages but also as a speciation pump across multi-glacial cycles for marine organisms otherwise exhibiting cosmopolite amphi-boreal distributions.

Displaced phylogeographic signals from Gyrodactylus arcuatus, a parasite of the three-spined stickleback Gasterosteus aculeatus, suggest freshwater glacial refugia in Europe

We examined the global mitochondrial phylogeography of Gyrodactylus arcuatus, a flatworm ectoparasite of three-spined stickleback Gasterosteus aculeatus. In accordance with the suggested high divergence rate of 13%/million years, the genetic variation of the parasite was high: haplotype diversity h=0.985 and nucleotide diversity π=0.0161. The differentiation among the parasite populations was substantial (Φst=0.759), with two main allopatric clades (here termed Euro and North) accounting for 54% of the total genetic variation. The diversity center of the Euro clade was in the Baltic Sea, while the North clade was spread across the Barents and White Seas. A single haplotype within the North clade was found in the western and eastern Pacific Ocean. Divergence of main clades was estimated to be circa 200 thousand years ago. Each main clade was further divided into six distinct subclades, estimated to have diverged in isolation since 135 thousand years ago. This second division corresponds approximately to the Eemian interglacial predating the last glacial maximum. A demographic expansion of the subclades is associated with colonisation of northern Europe since the last glacial maximum, circa 15-40 thousand years ago. The parasite phylogeny is most likely explained by sequential isolated bottlenecks and expansions in numerous allopatric refugia. The postglacial intermingling and high variation in the marine parasite populations, separately in the Baltic and Barents Seas, suggest low competition of divergent parasite matrilines, coupled with a large population size and high rate of dispersal of hosts. The genetic contribution of the assumed refugial fish populations maintaining the parasite during the last glacial maximum was not detected among the marine sticklebacks, which perhaps were infected after range expansion.

Latitudinal Gradient in Otolith Shape among Local Populations of Atlantic Herring (Clupea harengus L.) in Norway

Otolith shape analysis of Atlantic herring (Clupea harengus) in Norwegian waters shows significant differentiation among fjords and a latitudinal gradient along the coast where neighbouring populations are more similar to each other than to those sampled at larger distances. The otolith shape was obtained using quantitative shape analysis, the outlines were transformed with Wavelet and analysed with multivariate methods. The observed morphological differences are likely to reflect environmental differences but indicate low dispersal among the local herring populations. Otolith shape variation suggests also limited exchange between the local populations and their oceanic counterparts, which could be due to differences in spawning behaviour. Herring from the most northerly location (69°N) in Balsfjord, which is genetically more similar to Pacific herring (Clupea pallasii), differed in otolith shape from all the other populations. Our results suggest that the semi-enclosed systems, where the local populations live and breed, are efficient barriers for dispersal. Otolith shape can thus serve as a marker to identify the origin of herring along the coast of Norway.