Phylogeography of the planktonic chaetognath Sagitta setosa reveals isolation in European seas

Highly restricted dispersal in habitat-forming seaweed may impede natural recovery of disturbed populations

Cystoseira sensu lato (Class Phaeophyceae, Order Fucales, Family Sargassaceae) forests play a central role in marine Mediterranean ecosystems. Over the last decades, Cystoseira s.l. suffered from a severe loss as a result of multiple anthropogenic stressors. In particular, Gongolaria barbata has faced multiple human-induced threats, and, despite its ecological importance in structuring rocky communities and hosting a large number of species, the natural recovery of G. barbata depleted populations is uncertain. Here, we used nine microsatellite loci specifically developed for G. barbata to assess the genetic diversity of this species and its genetic connectivity among fifteen sites located in the Ionian, the Adriatic and the Black Seas. In line with strong and significant heterozygosity deficiencies across loci, likely explained by Wahlund effect, high genetic structure was observed among the three seas (ENA corrected FST = 0.355, IC = [0.283, 0.440]), with an estimated dispersal distance per generation smaller than 600 m, both in the Adriatic and Black Sea. This strong genetic structure likely results from restricted gene flow driven by geographic distances and limited dispersal abilities, along with genetic drift within isolated populations. The presence of genetically disconnected populations at small spatial scales (< 10 km) has important implications for the identification of relevant conservation and management measures for G. barbata: each population should be considered as separated evolutionary units with dedicated conservation efforts.

Oceanic dispersal barriers in a holoplanktonic gastropod

Pteropods, a group of holoplanktonic gastropods, are regarded as bioindicators of the effects of ocean acidification on open ocean ecosystems, because their thin aragonitic shells are susceptible to dissolution. While there have been recent efforts to address their capacity for physiological acclimation, it is also important to gain predictive understanding of their ability to adapt to future ocean conditions. However, little is known about the levels of genetic variation and large-scale population structuring of pteropods, key characteristics enabling local adaptation. We examined the spatial distribution of genetic diversity in the mitochondrial cytochrome c oxidase I (COI) and nuclear 28S gene fragments, as well as shell shape variation, across a latitudinal transect in the Atlantic Ocean (35°N-36°S) for the pteropod Limacina bulimoides. We observed high levels of genetic variability (COI π = 0.034, 28S π = 0.0021) and strong spatial structuring (COI ΦST  = 0.230, 28S ΦST  = 0.255) across this transect. Based on the congruence of mitochondrial and nuclear differentiation, as well as differences in shell shape, we identified a primary dispersal barrier in the southern Atlantic subtropical gyre (15-18°S). This barrier is maintained despite the presence of expatriates, a gyral current system, and in the absence of any distinct oceanographic gradients in this region, suggesting that reproductive isolation between these populations must be strong. A secondary dispersal barrier supported only by 28S pairwise ΦST comparisons was identified in the equatorial upwelling region (between 15°N and 4°S), which is concordant with barriers observed in other zooplankton species. Both oceanic dispersal barriers were congruent with regions of low abundance reported for a similar basin-scale transect that was sampled 2 years later. Our finding supports the hypothesis that low abundance indicates areas of suboptimal habitat that result in barriers to gene flow in widely distributed zooplankton species. Such species may in fact consist of several populations or (sub)species that are adapted to local environmental conditions, limiting their potential for adaptive responses to ocean changes. Future analyses of genome-wide diversity in pteropods could provide further insight into the strength, formation and maintenance of oceanic dispersal barriers.

Genetic diversity and connectivity of Flaccisagitta enflata (Chaetognatha: Sagittidae) in the tropical Atlantic ocean (northeastern Brazil)

The phylogeography of the holoplanktonic chaetognath Flaccisagitta enflata was investigated in the Tropical Western Atlantic (TWA). Considering the cosmopolitan range of this species and the fact that its entire life cycle is planktonic, the central hypothesis of this study is that F. enflata exhibits connectivity due to its high dispersal capacity, forming a panmictic population among the study sites. The evaluated areas included neritic (Port of Recife–PR, and Tamandaré - TA) and oceanic (Fernando de Noronha Archipelago—FN, Rocas Atoll—RA, Guará seamount—GS and Saint Peter and Saint Paul’s Archipelago—SPSPA) locations of the Brazilian Blue Amazon. We used COI gene sequences as molecular marker. Partial sequences (425 bp) were obtained for 116 specimens and employed to reconstruct the phylogeny, build an haplotype network, evaluate gene flow through a migration model, and estimate diversity indices, population structuring and demographic history. High levels of haplotype diversity (mean: 0.98) and moderate to high levels of nucleotide diversity (mean: 0.023) were observed. The phylogeny and the haplotype network topologies showed some geographic clustering, indicating local structuring in GS and PR. This finding was supported by the AMOVA high global Φst (0.033, significant) and some pairwise Φst comparisons (7 out of 15 were significantly >0). Significant differences suggested lower levels of connectivity when GS population was compared to those of FN and SPSPA; as well as when TA was compared to FN. These results might be related to particularities of the oceanic dynamics which rules the TWA, sustaining such dissimilarities. Structuring was also observed between PR and all oceanic locations. We hypothesize that the topography of the port inlet, enclosured by a reef barrier, may constrain the water turnover ratio and thus migration rates of F. enflata in the TWA. Accordingly, Migrate-N yielded a four metapopulations model (PR ⇌ TA ⇌ SPSPA+FN ⇌ GS+RA) as the best (highest probability; ~0.90) to represent the structuring of F. enflata in the TWA. Therefore, the null hypothesis of one randomly mating population cannot be accepted. The demographic evaluation demonstrated that the neutral hypothesis of stable populations may not be rejected for most of the locations. This work is the start point to broaden the knowledge on the phylogeography and population genetic structure of a numerically dominant species in the Western Atlantic, with key role in the marine trophic web.

Extreme Mitogenomic Variation in Natural Populations of Chaetognaths

The extent of within-species genetic variation across the diversity of animal life is an underexplored problem in ecology and evolution. Although neutral genetic variation should scale positively with population size, mitochondrial diversity levels are believed to show little variation across animal species. Here, we report an unprecedented case of extreme mitochondrial diversity within natural populations of two morphospecies of chaetognaths (arrow worms). We determine that this diversity is composed of deep sympatric mitochondrial lineages, which are in some cases as divergent as human and platypus. Additionally, based on 54 complete mitogenomes, we observed mitochondrial gene order differences between several of these lineages. We examined nuclear divergence patterns (18S, 28S, and an intron) to determine the possible origin of these lineages, but did not find congruent patterns between mitochondrial and nuclear markers. We also show that extreme mitochondrial divergence in chaetognaths is not driven by positive selection. Hence, we propose that the extreme levels of mitochondrial variation could be the result of either a complex scenario of reproductive isolation, or a combination of large population size and accelerated mitochondrial mutation rate. These findings emphasize the importance of characterizing genome-wide levels of nuclear variation in these species and promote chaetognaths as a remarkable model to study mitochondrial evolution.

Unexpected mosaic distribution of two hybridizing sibling lineages in the teleplanically dispersing snail Stramonita haemastoma suggests unusual postglacial redistribution or cryptic invasion

Molecular approaches have proven efficient to identify cryptic lineages within single taxonomic entities. Sometimes these cryptic lineages maybe previously unreported or unknown invasive taxa. The genetic structure of the marine gastropod Stramonita haemastoma has been examined in the Western Mediterranean and North‐Eastern Atlantic populations with mtDNA COI sequences and three newly developed microsatellite markers. We identified two cryptic lineages, differentially fixed for alternative mtDNA COI haplogroups and significantly differentiated at microsatellite loci. The mosaic distribution of the two lineages is unusual for a warm‐temperate marine invertebrate with a teleplanic larval stage. The Atlantic lineage was unexpectedly observed as a patch enclosed in the north of the Western Mediterranean Sea between eastern Spain and the French Riviera, and the Mediterranean lineage was found in Macronesian Islands. Although cyto‐nuclear disequilibrium is globally maintained, asymmetric introgression occurs in the Spanish region where the two lineages co‐occur in a hybrid zone. A first interpretation of our results is mito‐nuclear discordance in a stable postglacial hybrid zone. Under this hypothesis, though, the location of genetic discontinuities would be unusual among planktonic dispersers. An alternative interpretation is that the Atlantic lineage, also found in Senegal and Venezuela, has been introduced by human activities in the Mediterranean area and is introgressing Mediterranean genes during its propagation, as theoretically expected. This second hypothesis would add an additional example to the growing list of cryptic marine invasions revealed by molecular studies.

Gene flow, population growth and a novel substitution rate estimate in a subtidal rock specialist, the black-faced blenny Tripterygion delaisi (Perciformes, Blennioidei, Tripterygiidae) from the Adriatic Sea

Population histories depend on the interplay between exogeneous and endogeneous factors. In marine species, phylogeographic and demographic patterns are often shaped by sea level fluctuations, water currents and dispersal ability. Using mitochondrial control region sequences (n = 120), we infer phylogeographic structure and historic population size changes of a common littoral fish species, the black-faced blenny Tripterygion delaisi (Perciformes, Blennioidei, Tripterygiidae) from the north-eastern Adriatic Sea. We find that Adriatic T. delaisi are differentiated from conspecific populations in the remaining Mediterranean, but display little phylogeographic structure within the Adriatic basin. The pattern is consistent with passive dispersal of planktonic larvae along cyclonic currents within the Adriatic Sea, but limited active dispersal of adults. Demographic reconstructions are consistent with recent population expansion, probably triggered by rising sea levels after the last glacial maximum (LGM). Placing the onset of population growth between the LGM and the warming of surface waters (18 000-13 000 years BP) and employing a novel expansion dating approach, we inferred a substitution rate of 2.61-3.61% per site per MY. Our study is one of only few existing investigations of the genetic structure of animals within the Adriatic basin and is the first to provide an estimate for mitochondrial control region substitution rates in blennioid fishes.

Temporal Stability of Genetic Structure in a Mesopelagic Copepod

Although stochasticity in oceanographic conditions is known to be an important driver of temporal genetic change in many marine species, little is known about whether genetically distinct plankton populations can persist in open ocean habitats. A prior study demonstrated significant population genetic structure among oceanic gyres in the mesopelagic copepod Haloptilus longicornis in both the Atlantic and Pacific Oceans, and we hypothesized that populations within each gyre represent distinct gene pools that persist over time. We tested this expectation through basin-scale sampling across the Atlantic Ocean in 2010 and 2012. Using both mitochondrial (mtCOII) and microsatellite markers (7 loci), we show that the genetic composition of populations was stable across two years in both the northern and southern subtropical gyres. Genetic variation in this species was partitioned among ocean gyres (F_CT = 0.285, P < 0.0001 for mtCOII, F_CT = 0.013, P < 0.0001 for microsatellites), suggesting strong spatial population structure, but no significant partitioning was found among sampling years. This temporal persistence of population structure across a large geographic scale was coupled with chaotic genetic patchiness at smaller spatial scales, but the magnitude of genetic differentiation was an order of magnitude lower at these smaller scales. Our results demonstrate that genetically distinct plankton populations persist over time in highly-dispersive open ocean habitats, and this is the first study to rigorously test for temporal stability of large scale population structure in the plankton.

Multilocus evidence for globally distributed cryptic species and distinct populations across ocean gyres in a mesopelagic copepod

Zooplanktonic taxa have a greater number of distinct populations and species than might be predicted based on their large population sizes and open-ocean habitat, which lacks obvious physical barriers to dispersal and gene flow. To gain insight into the evolutionary mechanisms driving genetic diversification in zooplankton, we developed eight microsatellite markers to examine the population structure of an abundant, globally distributed mesopelagic copepod, Haloptilus longicornis, at 18 sample sites across the Atlantic and Pacific Oceans (n = 761). When comparing our microsatellite results with those of a prior study that used a mtDNA marker (mtCOII, n = 1059, 43 sample sites), we unexpectedly found evidence for the presence of a cryptic species pair. These species were globally distributed and apparently sympatric, and were separated by relatively weak genetic divergence (reciprocally monophyletic mtCOII lineages 1.6% divergent; microsatellite FST ranging from 0.28 to 0.88 across loci, P < 0.00001). Using both mtDNA and microsatellite data for the most common of the two species (n = 669 for microsatellites, n = 572 for mtDNA), we also found evidence for allopatric barriers to gene flow within species, with distinct populations separated by continental landmasses and equatorial waters in both the Atlantic and Pacific Ocean basins. Our study shows that oceanic barriers to gene flow can act as a mechanism promoting allopatric diversification in holoplanktonic taxa, despite the high potential dispersal abilities and pelagic habitat for these species.

Evidence of cryptic and pseudocryptic speciation in the Paracalanus parvus species complex (Crustacea, Copepoda, Calanoida)

Introduction

Many marine planktonic crustaceans such as copepods have been considered as widespread organisms. However, the growing evidence for cryptic and pseudo-cryptic speciation has emphasized the need of re-evaluating the status of copepod species complexes in molecular and morphological studies to get a clearer picture about pelagic marine species as evolutionary units and their distributions. This study analyses the molecular diversity of the ecologically important Paracalanus parvus species complex. Its seven currently recognized species are abundant and also often dominant in marine coastal regions worldwide from temperate to tropical oceans.

Results

COI and Cytochrome b sequences of 160 specimens of the Paracalanus parvus complex from all oceans were obtained. Furthermore, 42 COI sequences from GenBank were added for the genetic analyses. Thirteen distinct molecular operational taxonomic units (MOTU) and two single sequences were revealed with cladistic analyses (Maximum Likelihood, Bayesian Inference), of which seven were identical with results from species delimitation methods (barcode gaps, ABDG, GMYC, Rosenberg’s P(AB)). In total, 10 to 12 putative species were detected and could be placed in three categories: (1) temperate geographically isolated, (2) warm-temperate to tropical wider spread and (3) circumglobal warm-water species.

Conclusions

The present study provides evidence of cryptic or pseudocryptic speciation in the Paracalanus parvus complex. One major insight is that the species Paracalanus parvus s.s. is not panmictic, but may be restricted in its distribution to the northeastern Atlantic.

Phylogeographic analysis of introns and mitochondrial DNA in the clam Ruditapes decussatus uncovers the effects of Pleistocene glaciations and endogenous barriers to gene flow

Studies on the phylogeography of species inhabiting the Mediterranean and the nearby coasts of the NE Atlantic Ocean (MEDAT) have found subdivision and/or phylogeographic structure in one or more of the Atlantic, western Mediterranean and eastern Mediterranean basins. This structure has been explained as the result of past population fragmentation caused by Pleistocene sea level changes and current patterns of marine circulation. However, the increasing use of nuclear markers has revealed that these two factors alone are not enough to explain the phylogeographic patterns, and an additional role has been suggested for endogenous barriers to gene flow or natural selection. In this article we examined the role of these factors in Ruditapes decussatus, a commercial clam species native to MEDAT. A genetic analysis of 11 populations was carried out by examining 6 introns with a PCR-RFLP technique. We found subdivision in three regions: Atlantic (ATL), western Mediterranean plus Tunisia (WMED), and Aegean and Adriatic seas (AEGAD). Two introns (Ech and Tbp) showed alleles that were restricted to AEGAD. Sequencing a subsample of individuals for these introns indicated that AEGAD-specific alleles were separate clades, thus revealing a phylogeographic brake at the WMED-AEGAD boundary. Sequencing of the mitochondrial COI locus confirmed this phylogeographic break. Dating of the AEGAD mitochondrial haplotypes and nuclear alleles with a Bayesian MCMC method revealed that they shared common ancestors in the Pleistocene. These results can be explained in the framework of Pleistocene sea level drops and patterns of gene flow in MEDAT. An additional observation was a lack of differentiation at COI between the ATL and WMED, in sharp contrast with 4 introns that showed clear genetic subdivision. Neutrality tests did not support the hypothesis of a selective sweep acting on mtDNA to explain the contrasting levels of differentiation between mitochondrial and nuclear markers across the ATL-WMED transition, and we argue that the difference between markers is best explained by the existence of an endogenous genetic barrier, rather than by a physical barrier to larval migration alone.

High cryptic diversity across the global range of the migratory planktonic copepods Pleuromamma piseki and P. gracilis

Although holoplankton are ocean drifters and exhibit high dispersal potential, a number of studies on single species are finding highly divergent genetic clades. These cryptic species complexes are important to discover and describe, as identification of common marine species is fundamental to understanding ecosystem dynamics. Here we investigate the global diversity within Pleuromamma piseki and P. gracilis, two dominant members of the migratory zooplankton assemblage in subtropical and tropical waters worldwide. Using DNA sequence data from the mitochondrial gene cytochrome c oxidase subunit II (mtCOII) from 522 specimens collected across the Pacific, Atlantic and Indian Oceans, we discover twelve well-resolved genetically distinct clades in this species complex (Bayesian posterior probabilities >0.7; 6.3–17% genetic divergence between clades). The morphologically described species P. piseki and P. gracilis did not form monophyletic groups, rather they were distributed throughout the phylogeny and sometimes co-occurred within well-resolved clades: this result suggests that morphological characters currently used for taxonomic identification of P. gracilis and P. piseki may be inaccurate as indicators of species’ boundaries. Cryptic clades within the species complex ranged from being common to rare, and from cosmopolitan to highly restricted in distribution across the global ocean. These novel lineages appear to be ecologically divergent, with distinct biogeographic distributions across varied pelagic habitats. We hypothesize that these mtDNA lineages are distinct species and suggest that resolving their systematic status is important, given the ecological significance of the genus Pleuromamma in subtropical-tropical waters worldwide.

Phylogeography of the sponge Suberites diversicolor in Indonesia: insights into the evolution of marine lake populations

The existence of multiple independently derived populations in landlocked marine lakes provides an opportunity for fundamental research into the role of isolation in population divergence and speciation in marine taxa. Marine lakes are landlocked water bodies that maintain a marine character through narrow submarine connections to the sea and could be regarded as the marine equivalents of terrestrial islands. The sponge Suberites diversicolor (Porifera: Demospongiae: Suberitidae) is typical of marine lake habitats in the Indo-Australian Archipelago. Four molecular markers (two mitochondrial and two nuclear) were employed to study genetic structure of populations within and between marine lakes in Indonesia and three coastal locations in Indonesia, Singapore and Australia. Within populations of S. diversicolor two strongly divergent lineages (A & B) (COI: p = 0.4% and ITS: p = 7.3%) were found, that may constitute cryptic species. Lineage A only occurred in Kakaban lake (East Kalimantan), while lineage B was present in all sampled populations. Within lineage B, we found low levels of genetic diversity in lakes, though there was spatial genetic population structuring. The Australian population is genetically differentiated from the Indonesian populations. Within Indonesia we did not record an East-West barrier, which has frequently been reported for other marine invertebrates. Kakaban lake is the largest and most isolated marine lake in Indonesia and contains the highest genetic diversity with genetic variants not observed elsewhere. Kakaban lake may be an area where multiple putative refugia populations have come into secondary contact, resulting in high levels of genetic diversity and a high number of endemic species.

High evolutionary potential of marine zooplankton

Open ocean zooplankton often have been viewed as slowly evolving species that have limited capacity to respond adaptively to changing ocean conditions. Hence, attention has focused on the ecological responses of zooplankton to current global change, including range shifts and changing phenology. Here, we argue that zooplankton also are well poised for evolutionary responses to global change. We present theoretical arguments that suggest plankton species may respond rapidly to selection on mildly beneficial mutations due to exceptionally large population size, and consider the circumstantial evidence that supports our inference that selection may be particularly important for these species. We also review all primary population genetic studies of open ocean zooplankton and show that genetic isolation can be achieved at the scale of gyre systems in open ocean habitats (100s to 1000s of km). Furthermore, population genetic structure often varies across planktonic taxa, and appears to be linked to the particular ecological requirements of the organism. In combination, these characteristics should facilitate adaptive evolution to distinct oceanographic habitats in the plankton. We conclude that marine zooplankton may be capable of rapid evolutionary as well as ecological responses to changing ocean conditions, and discuss the implications of this view. We further suggest two priority areas for future research to test our hypothesis of high evolutionary potential in open ocean zooplankton, which will require (1) assessing how pervasive selection is in driving population divergence and (2) rigorously quantifying the spatial and temporal scales of population differentiation in the open ocean. Recent attention has focused on the ecological responses of open ocean zooplankton to current global change, including range shifts and changing phenology. Here, we argue that marine zooplankton also are well poised for evolutionary responses to global change.

Population differentiation in the open sea: insights from the pelagic copepod Pleuromamma xiphias

Although a number of recent studies of marine holoplankton have reported significant genetic structure among populations, little is currently known about the biological and oceanographic processes that influence population connectivity in oceanic plankton. In order to examine how depth preferences influence dispersal in oceanic plankton, I characterized the genetic structure of a copepod with diel vertical migration (DVM) (Pleuromamma xiphias), throughout its global distribution, and compared these results to those expected given the interaction of this species' habitat depth with ocean circulation and bathymetry. Mitochondrial COI sequences from 651 individuals from 28 sites in the Indian, Pacific, and Atlantic Oceans revealed highly significant genetic differentiation both within and among ocean basins. Limited dispersal among distinct pelagic provinces seems to have played a major role in population differentiation in this species, with strong genetic breaks observed across known oceanographic fronts or current systems in all three ocean basins. The Indo-West Pacific (IWP) holds a highly distinct genetic population of this species that was sampled in both the western Pacific and eastern Indian Oceans. This suggests that the IWP does not act as a strong barrier to gene flow between basins, as expected, despite the relatively shallow water depth (<200 m) and vertically extensive (>400 m) diel migration of this species. A pattern of isolation by distance was observed in the Indian Ocean with genetic differentiation among samples down to spatial scales of ∼800 km, indicating that realized dispersal in P. xiphias occurs over much smaller spatial scales than in previously reported oceanic holoplankton. Given its highly regionalized population genetic structure, P. xiphias may have some capacity to adapt to local oceanographic conditions, and it should not be assumed that populations of this species in distinct pelagic biomes will respond in the same way to shared physical or climatic forcing.

High level of structural polymorphism driven by mobile elements in the Hox genomic region of the Chaetognath Spadella cephaloptera

Little is known about the relationships between genome polymorphism, mobile element dynamics, and population size among animal populations. The chaetognath species Spadella cephaloptera offers a unique perspective to examine this issue because they display a high level of genetic polymorphism at the population level. Here, we have investigated in detail the extent of nucleotide and structural polymorphism in a region harboring Hox1 and several coding genes and presumptive functional elements. Sequencing of several bacterial artificial chromosome inserts representative of this nuclear region uncovered a high level of structural heterogeneity, which is mainly caused by the polymorphic insertion of a diversity of genetic mobile elements. By anchoring this variation through individual genotyping, we demonstrated that sequence diversity could be attributed to the allelic pool of a single population, which was confirmed by detection of extensive recombination within the genomic region studied. The high average level of nucleotide heterozygosity provides clues of selection in both coding and noncoding domains. This pattern stresses how selective processes remarkably cope with intense sequence turnover due to substitutions, mobile element insertions, and recombination to preserve the integrity of functional landscape. These findings suggest that genome polymorphism could provide pivotal information for future functional annotation of genomes.

Lack of genetic structure in the jellyfish Pelagia noctiluca (Cnidaria: Scyphozoa: Semaeostomeae) across European seas

The genetic structure of the holopelagic scyphozoan Pelagia noctiluca was inferred based on the study of 144 adult medusae. The areas of study were five geographic regions in two European seas (Eastern Atlantic and Mediterranean Sea). A 655-bp sequence of mitochondrial cytochrome c oxidase subunit I (COI), and a 645-bp sequence of two nuclear internal transcribed spacers (ITS1 and ITS2) were analyzed. The protein coding COI gene showed a higher level of divergence than the combined nuclear ITS fragment (haplotype diversity 0.962 vs. 0.723, nucleotide diversity 1.16% vs. 0.31%). Phylogeographic analysis on COI gene revealed two clades, the larger consisting of specimens from all sampling sites, and the smaller mostly formed of specimens from the Mediterranean Sea. Haplotype diversity was very high throughout the sampled area, and within sample diversity was higher than diversity among geographical regions. No strongly supported genetically or geographically distinct groups of P. noctiluca were found. The results - long distance dispersal, insignificant F(ST) values, lack of isolation by distance - pointed toward an admixture among Mediterranean and East Atlantic populations.

Barcoding of arrow worms (Phylum Chaetognatha) from three oceans: genetic diversity and evolution within an enigmatic phylum

Arrow worms (Phylum Chaetognatha) are abundant planktonic organisms and important predators in many food webs; yet, the classification and evolutionary relationships among chaetognath species remain poorly understood. A seemingly simple body plan is underlain by subtle variation in morphological details, obscuring the affinities of species within the phylum. Many species achieve near global distributions, spanning the same latitudinal bands in all ocean basins, while others present disjunct ranges, in some cases with the same species apparently found at both poles. To better understand how these complex evolutionary and geographic variables are reflected in the species makeup of chaetognaths, we analyze DNA barcodes of the mitochondrial cytochrome oxidase c subunit I (COI) gene, from 52 specimens of 14 species of chaetognaths collected mainly from the Atlantic Ocean. Barcoding analysis was highly successful at discriminating described species of chaetognaths across the phylum, and revealed little geographical structure. This barcode analysis reveals hitherto unseen genetic variation among species of arrow worms, and provides insight into some species relationships of this enigmatic group.

Mitochondrial phylogeography of the European sprat (Sprattus sprattus L., Clupeidae) reveals isolated climatically vulnerable populations in the Mediterranean Sea and range expansion in the northeast Atlantic

We examined the genetic structure of the European sprat (Sprattus sprattus) by means of a 530-bp sequence of the mitochondrial control region from 210 fish originating from seven sampling localities of its distributional range. Phylogeographical analysis of 128 haplotypes showed a phylogenetic separation into two major clades with the Strait of Sicily acting as a barrier to gene flow between them. While no population differentiation was observed based on analysis of molecular variance and net nucleotide differences between samples of the Baltic Sea, the North Sea and the Bay of Biscay nor between the Black Sea and the Bosporus, a strong population differentiation between these samples and two samples from the Mediterranean Sea was found. Further, the biggest genetic distance was observed within the Mediterranean Sea between the populations of the Gulf of Lyon and the Adriatic Sea, indicating genetic isolation of these regions. Low genetic diversities and star-like haplotype networks of both Mediterranean Sea populations point towards recent demographic expansion scenarios after low population size, which is further supported by negative F(S) values and unimodal mismatch distributions with a low mean. Along the northeast Atlantic coast, a northwards range expansion of a large and stable population can be assumed. The history of a diverse but differentiated Black Sea population remains unknown due to uncertainties in the palaeo-oceanography of this sea. Our genetic data did not confirm the presently used classification into subspecies but are only preliminary in the absence of nuclear genetic analyses.

Cryptic ecological diversification of a planktonic estuarine copepod, Acartia tonsa

The recent discovery of cryptic species in marine holoplankton, organisms that 'drift' in oceanic currents throughout their life cycle, contrasts with their potential for long-distance passive dispersal and presumably high gene flow. These observations suggest that holoplankton species are adapting to surprisingly small-scale oceanographic features and imply either limited dispersal or strong selection gradients. Acartia tonsa is a widespread and numerically dominant estuarine copepod containing deep mitochondrial lineages within and among populations along the northwestern Atlantic coast. In this study, we intensively investigated A. tonsa populations in Chesapeake Bay with the goals of testing species status for the deep lineages and testing for their association with environmental features over space and time. Phylogenetic analyses of DNA sequences from mitochondrial cytochrome c oxidase I (mtCOI) and the nuclear ribosomal internal transcribed spacer (nITS) resolved two concordant monophyletic clades. Deep divergence between the two clades (13.7% uncorrected sequence divergence for mtCOI and 32.2% for nITS) and genealogical concordance within sympatric populations strongly suggest that the two clades represent reproductively isolated cryptic species. Based on restriction fragment length polymorphisms of mtCOI, representatives from the two clades were found consistently associated with contrasting salinity regimes (oligohaline vs. meso-polyhaline) with an overlap between 2 and 12 PSU in samples from 1995 to 2005. Finding these patterns in one of the best-known estuarine copepods reinforces the conclusion that marine biodiversity is underestimated, not only in terms of species numbers, but also with respect to niche partitioning and the potential importance of ecological divergence in marine holoplankton.

Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break?

The geological history of the Mediterranean Sea, its hydrography and connection with the Atlantic Ocean have been well documented. Despite a wealth of historical and oceanographic data, the Atlantic-Mediterranean transition remains controversial at the biological level as there are discordant results regarding the biogeographical separation between the Atlantic and Mediterranean biota. The opening of the Strait of Gibraltar at the end of the Messinian Salinity Crisis (some 5.33 million years ago), removed the land barrier that impeded the marine biota allowing it to disperse freely into the Mediterranean Sea. However, present day genetic patterns suggest a limitation to gene flow for some marine species, preventing population admixture. In the last few years, a large number of studies have challenged the hypothesis of the Strait of Gibraltar representing a phylogeographical break. A review of more than 70 papers reveals no obvious relationship between either dispersal ability or life history, and observed patterns of partial or complete genetic isolation between Atlantic and Mediterranean populations. We re-analysed a selection of this large body of data (20 studies in total) in order to provide a homogeneous and coherent view on the generality of the phylogeographical patterns and the presence of a phylogeographical barrier. This offered the opportunity to summarize the state of the art on this matter and reach some general conclusions on the evolutionary history across the Atlantic-Mediterranean range. Geographically, some species in the transition zone showed step changes of allele frequencies associated with the Almeria-Oran Front rather than with the Strait of Gibraltar itself. A major part of the data describe evolutionary events well within the time frame of the Quaternary age as very few taxa pre-date closure of the Tethys Sea. Results point to a combined signature of vicariance, palaeoclimate fluctuation and life-history traits on the Atlantic-Mediterranean phylogeographical patterns. Principal component analysis failed to show any particular association between biological traits and genetic variables. It would argue that organismal determinism may play a far less significant role than marine biogeographers have generally believed.

Testing hypotheses of population structuring in the Northeast Atlantic Ocean and Mediterranean Sea using the common cuttlefish Sepia officinalis

Population structuring in species inhabiting marine environments such as the Northeast Atlantic Ocean (NEA) and Mediterranean Sea (MS) has usually been explained based on past and present physical barriers to gene flow and isolation by distance (IBD). Here, we examined the relative importance of these factors on population structuring of the common cuttlefish Sepia officinalis by using methods of phylogenetic inference and hypothesis testing coupled with coalescent and classical population genetic parameter estimation. Individuals from 10 Atlantic and 15 Mediterranean sites were sequenced for 659 bp of the mitochondrial COI gene (259 sequences). IBD seems to be the main factor driving present and past genetic structuring of Sepia populations across the NEA-MS, both at large and small geographical scales. Such an evolutionary process agrees well with some of the biological features characterizing this cuttlefish species (short migrations, nektobenthic habit, benthic eggs hatching directly to benthic juveniles). Despite the many barriers to migration/gene flow suggested in the NEA-MS region, genetic population fragmentation due to past isolation of water masses (Pleistocene; 0.56 million years ago) and/or present-day oceanographic currents was only detected between the Aegean-Ionian and western Mediterranean Seas. Restricted gene flow associated with the Almería-Oran hydrographic front was also suggested between southern and eastern Spanish populations. Distinct population boundaries could not be clearly determined, except for the Aegean-Ionian stock. Two Atlantic and five Mediterranean samples showed evidence of current decline in genetic diversity, which may indicate over-exploitation of Sepia in both marine regions.

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.

Spatial and temporal genetic structure of the planktonic Sagitta setosa (Chaetognatha) in European seas as revealed by mitochondrial and nuclear DNA markers

Little is known about the spatial and temporal scales at which planktonic organisms are genetically structured. A previous study of mitochondrial DNA (mtDNA) in the holoplanktonic chaetognath Sagitta setosa revealed strong phylogeographic structuring suggesting that Northeast (NE) Atlantic, Mediterranean and Black Sea populations are genetically disjunct. The present study used a higher sampling intensity and a combination of mitochondrial and four microsatellite markers to reveal population structuring between and within basins. Between basins, both marker sets indicated significant differentiation confirming earlier results that gene flow is probably absent between the respective S. setosa populations. At the within-basin scale, we found no evidence of spatial or temporal structuring within the NE Atlantic. In the Mediterranean basin, both marker sets indicated significant structuring, but only the mtDNA data indicated a sharp genetic division between Adriatic and all other Mediterranean populations. Data were inconclusive about population structuring in the Black Sea. The levels of differentiation indicated by the two marker sets differed substantially, with far less pronounced structure detected by microsatellite than mtDNA data. This study also uncovered the presence of highly divergent mitochondrial lineages that were discordant with morphology, geography and nuclear DNA. We thus propose the hypothesis that highly divergent mitochondrial lineages may be present within interbreeding S. setosa populations.

Hidden diversity and host specificity in cycliophorans: a phylogeographic analysis along the North Atlantic and Mediterranean Sea

In order to elucidate the evolutionary history and the population structure of the members of the phylum Cycliophora, which live commensally on three species of lobsters, we studied sequence variation in the mitochondrial gene cyctochrome c oxidase subunit I. Overall 242 sequences from 16 locations on both coasts of the North Atlantic, including the North Sea and the Mediterranean, were analysed, revealing 28 haplotypes, with a maximum sequence divergence of 16.6%. Total genetic diversity was high (h = 0.8322, pi = 0.0898), as it was for the commensals on Homarus americanus (17 haplotypes, h = 0.7506, pi = 0.0504). However, it was low for commensals on Nephrops norvegicus (6 haplotypes, h = 0.3899, pi = 0.0035), and intermediate for cycliophorans on Homarus gammarus (5 haplotypes, h = 0.3020, pi = 0.0140). Although two of the host lobsters co-inhabit the coastal waters of Europe, a strong genetic structure (78.45% of the observed genetic variation) was detected among populations on all host species, indicating the existence of a reproductively isolated species on each lobster. In addition, genetic structure over long distances exists among populations on each host species. Such patterns can be explained by the limited dispersal ability of the cycliophoran chordoid larva. Demographic and phylogenetic analyses suggest old and possibly cryptic populations present on H. americanus and H. gammarus, while the latter may have experienced recent bottlenecks, perhaps during Pleistocene glaciations. Populations on N. norvegicus appear to be of recent origin.

Chloroplast DNA phylogeography of European ashes,Fraxinussp. (Oleaceae): roles of hybridization and life history traits

We investigated range-wide phylogeographic variation in three European ash species (Fraxinus sp., Oleaceae). Chloroplast DNA (cpDNA) microsatellites were typed in the thermophilous Fraxinus angustifolia and Fraxinus ornus and the observed haplotypes and the geographic distribution of diversity were compared to cpDNA data previously obtained in the more cold-tolerant Fraxinus excelsior. We found wide-ranging haplotype sharing between the phylogenetically close F. angustifolia and F. excelsior, suggesting hybridization (i) in common glacial refuges in the Iberian Peninsula, northern Italy, the eastern and/or Dinaric Alps and the Balkan Peninsula, and/or (ii) during postglacial recolonization. The data allowed us to propose additional glacial refuges for F. angustifolia in southern Italy and in Turkey, and populations from the latter region were particularly polymorphic. There was evidence for refuge areas in Italy, the Balkan Peninsula and Turkey for F. ornus, which did not share any single chloroplast haplotype with the other species. In both F. angustifolia and F. ornus, cpDNA diversity (h(S) = 0.027 and h(S) = 0.009, respectively) was lower and fixation levels (G(ST) = 0.964 and G(ST) = 0.983, respectively) higher than in sympatric F. excelsior (h(S) = 0.096, G(ST) = 0.870). These diversity patterns could be due to temperature tolerance or the demographic history.

Mitochondrial DNA reveals a mosaic pattern of phylogeographical structure in Atlantic and Mediterranean populations of anchovy (Engraulis encrasicolus)

This study extends the geographic coverage of a previous study of mitochondrial DNA restriction fragment length polymorphism in European anchovy. Both studies together include 24 samples representing 17 localities extending from the Black Sea, through the Mediterranean Sea to the eastern Atlantic as far south as Dakar, Senegal. Eighty-eight haplotypes define two clades (A and B) separated by 3.2% sequence divergence. Clade A has a star-like genealogy indicative of a recent population expansion. Clade B has a more complex genealogy, consisting of several haplotypes at intermediate frequencies. The distributions of these clades consist of a mosaic with abrupt changes between some areas and gradients between other areas. Clade A predominates the Black and Aegean seas, but is present throughout the Mediterranean. Unexpectedly, new data show that clade A is also at a high frequency in the Atlantic, from Portugal to at least Senegal. Overall, the level of genetic differentiation among populations is high (F(ST)=0.148, p<0.0001), with the greatest differences between basins. AMOVA reveals four main geographical groups: Atlantic, central Mediterranean, Aegean Sea, and Black Sea. Mismatch distribution clearly indicates historical bottleneck and population expansion for clade A, while for clade B such evidence is equivocal. This difference may reflect a range expansion for both clades, but with higher gene flow (Nm values) between demes for clade A. Both contemporary and historical processes are important in shaping the complex genetic population structure of European anchovy.