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

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Mitonuclear genetic patterns of divergence in the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) along the Turkish seas

Biogeographical transition zones present good opportunities for studying the effect of the past ice ages on genetic structure of species because secondary contact zones of post-glacial lineages can be formed. In this study, we investigated the population genetic structure of the marbled rock crab, Pachygrapsus marmoratus along the coasts of Turkey. We genotyped 334 individuals from the Black Sea, the Turkish Straits System (TSS), the Aegean, and the Eastern Mediterranean basins. In order to reveal its evolutionary history and its population connectivity, we used mitochondrial CO1 region and five microsatellite loci. CO1 analyzes also included 610 additional samples from Genbank, which covered most of its distribution range. Both microsatellites and mtDNA showed decreased diversity in sampling sites of the TSS and the Black Sea as compared to those along the Aegean and the Levantine coasts. There is an especially strong geographical pattern in distributions of haplotypes in mtDNA, most probably as a result of genetic drift in the Black Sea and the Sea of Marmara (SoM). Microsatellite data analyses revealed two genetically distinct clusters of P. marmoratus (clusters C and M). While individuals belonging to cluster C are present in all the sampling locations, those belonging to cluster M are only detected along the Mediterranean coasts including the Aegean and the Levantine basins. These clusters shared similar haplotypes in the Mediterranean. Haplotypes of two sympatric clusters could be similar due to incomplete lineage sorting of ancestral polymorphisms. In order to retrieve the complex demographic history and to investigate evolutionary processes resulting in sympatric clusters in the Aegean Sea and the Levantine basin, mitochondrial markers with faster mutation rates than CO1 and/or SNP data will be useful.

An introgression breakthrough left by an anthropogenic contact between two ascidians

Human-driven translocations of species have diverse evolutionary consequences such as promoting hybridization between previously geographically isolated taxa. This is well illustrated by the solitary tunicate, Ciona robusta, native to the North East Pacific and introduced in the North East Atlantic. It is now co-occurring with its congener Ciona intestinalis in the English Channel, and C. roulei in the Mediterranean Sea. Despite their long allopatric divergence, first and second generation crosses showed a high hybridization success between the introduced and native taxa in the laboratory. However, previous genetic studies failed to provide evidence of recent hybridization between C. robusta and C. intestinalis in the wild. Using SNPs obtained from ddRAD-sequencing of 397 individuals from 26 populations, we further explored the genome-wide population structure of the native Ciona taxa. We first confirmed results documented in previous studies, notably (i) a chaotic genetic structure at regional scale, and (ii) a high genetic similarity between C. roulei and C. intestinalis, which is calling for further taxonomic investigation. More importantly, and unexpectedly, we also observed a genomic hotspot of long introgressed C. robusta tracts into C. intestinalis genomes at several locations of their contact zone. Both the genomic architecture of introgression, restricted to a 1.5 Mb region of chromosome 5, and its absence in allopatric populations suggest introgression is recent and occurred after the introduction of the non-native species. Overall, our study shows that anthropogenic hybridization can be effective in promoting introgression breakthroughs between species at a late stage of the speciation continuum.

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.

Genetic diversity and novel lineages in the cosmopolitan copepod Pleuromamma abdominalis in the Southeast Pacific

Across boundary currents, zooplankton are subject to strong oceanographic gradients and hence strong selective pressures. How such gradients interact with the speciation process of pelagic organisms is still poorly understood in the open ocean realm. Here we report on genetic diversity within the pelagic copepod Pleuromamma abdominalis in the poorly known Southeast Pacific region, with samples spanning an ocean gradient from coastal upwelling to the oligotrophic South Pacific Subtropical Gyre. We assessed variation in fragments of the mitochondrial (mt) genes cytochrome c oxidase subunit I (COI) and Cytochrome b as well as in the nuclear internal transcribed spacer (ITS) region and 28 S rRNA. Phylogenetic analyses revealed the presence of 8 divergent lineages occurring across the gradient with genetic distances in the range of 0.036 and 0.44 (mt genes), and GMYC species delimitation methods support their inference as distinct (undescribed) species. Genetic lineages occurring across the zonal gradient showed strong genetic structuring, with the presence of at least two new lineages within the coastal upwelling zone, revealing an unexpectedly high level of endemism within the Humboldt Current System. Multivariate analyses found strong correlation between genetic variation and surface chlorophyll-a and salinity, suggesting an important role for hydrographic gradients in maintaining genetic diversity. However, the presence of cryptic lineages within the upwelling zone cannot be easily accounted for by environmental heterogeneity and poses challenging questions for understanding the speciation process for oceanic zooplankton.

Single nucleotide polymorphisms reveal a genetic cline across the north-east Atlantic and enable powerful population assignment in the European lobster

Resolving stock structure is crucial for fisheries conservation to ensure that the spatial implementation of management is commensurate with that of biological population units. To address this in the economically important European lobster (Homarus gammarus), genetic structure was explored across the species' range using a small panel of single nucleotide polymorphisms (SNPs) previously isolated from restriction-site-associated DNA sequencing; these SNPs were selected to maximize differentiation at a range of both broad and fine scales. After quality control and filtering, 1,278 lobsters from 38 sampling sites were genotyped at 79 SNPs. The results revealed a pronounced phylogeographic break between the Atlantic and Mediterranean basins, while structure within the Mediterranean was also apparent, partitioned between lobsters from the central Mediterranean and the Aegean Sea. In addition, a genetic cline across the north-east Atlantic was revealed using both putatively neutral and outlier SNPs, but the precise driver(s) of this clinal pattern-isolation by distance, secondary contact, selection across an environmental gradient, or a combination of these factors-remains undetermined. Putatively neutral markers differentiated lobsters from Oosterschelde, an estuary on the Dutch coast, a finding likely explained by past bottlenecks and limited gene flow with adjacent North Sea populations. Building on the findings of our spatial genetic analysis, we were able to test the accuracy of assigning lobsters at various spatial scales, including to basin of origin (Atlantic or Mediterranean), region of origin and sampling location. The predictive model assembled using 79 SNPs correctly assigned 99.7% of lobsters not used to build the model to their basin of origin, but accuracy decreased to region of origin and again to sampling location. These results are of direct relevance to managers of lobster fisheries and hatcheries, and provide the basis for a genetic tool for tracing the origin of European lobsters in the food supply chain.

Genetic variation analysis of the cosmopolitan chaetognath Sagitta enflata in the northern South China Sea based on mitochondrial COI gene sequences

In this study, genetic diversity and population genetic structure of Sagitta enflata in the northern South China Sea were investigated by 623 bp fragment of mtDNA COI gene sequence. A total of 146 individuals were collected from nine stations and 92 different haplotypes were obtained. 485 variable sites (210 were parsimony informative and 275 were singleton variable sites), and no insertion or deletion was found. An analysis of molecular variance (AMOVA) and conventional population statistics (F_ST) revealed a low level of genetic differentiation among nine populations (F_ST = 0.14794, p < .05), indicating no geographical patterning among nine populations. The present results were able to provide a reference for the phylogenetic relationships and assessment of the genetic structure of S. enflata in the northern South China Sea.

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.

Parapatric genetic divergence among deep evolutionary lineages in the Mediterranean green crab, Carcinus aestuarii (Brachyura, Portunoidea, Carcinidae), accounts for a sharp phylogeographic break in the Eastern Mediterranean

BACKGROUND

Recently, population genetic studies of Mediterranean marine species highlighted patterns of genetic divergence and phylogeographic breaks, due to the interplay between impacts of Pleistocene climate shifts and contemporary hydrographical barriers. These factors markedly shaped the distribution of marine organisms and their genetic makeup. The present study is part of an ongoing effort to understand the phylogeography and evolutionary history of the highly dispersive Mediterranean green crab, Carcinus aestuarii (Nardo, 1847), across the Mediterranean Sea. Recently, marked divergence between two highly separated haplogroups (genetic types I and II) of C. aestuarii was discerned across the Siculo-Tunisian Strait, suggesting an Early Pleistocene vicariant event. In order to better identify phylogeographic patterns in this species, a total of 263 individuals from 22 Mediterranean locations were analysed by comparing a 587 basepair region of the mitochondrial gene Cox1 (cytochrome oxidase subunit 1). The examined dataset is composed of both newly generated sequences (76) and previously investigated ones (187).

RESULTS

Our results unveiled the occurrence of a highly divergent haplogroup (genetic type III) in the most north-eastern part of the Mediterranean Sea. Divergence between the most distinct type III and the common ancestor of both types I and II corresponds to the Early Pleistocene and coincides with the historical episode of separation between types I and II. Our results also revealed strong genetic divergence among adjacent regions (separating the Aegean and Marmara seas from the remaining distribution zone) and confirmed a sharp phylogeographic break across the Eastern Mediterranean. The recorded parapatric genetic divergence, with the potential existence of a contact zone between both groups in the Ionian Sea and notable differences in the demographic history, suggest the likely impact of paleoclimatic events, as well as past and contemporary oceanographic processes, in shaping genetic variability of this species.

CONCLUSIONS

Our findings not only provide further evidence for the complex evolutionary history of the green crab in the Mediterranean Sea, but also stress the importance of investigating peripheral areas in the species' distribution zone in order to fully understand the distribution of genetic diversity and unravel hidden genetic units and local patterns of endemism.

Comparative phylogeography of the ocean planet

Understanding how geography, oceanography, and climate have ultimately shaped marine biodiversity requires aligning the distributions of genetic diversity across multiple taxa. Here, we examine phylogeographic partitions in the sea against a backdrop of biogeographic provinces defined by taxonomy, endemism, and species composition. The taxonomic identities used to define biogeographic provinces are routinely accompanied by diagnostic genetic differences between sister species, indicating interspecific concordance between biogeography and phylogeography. In cases where individual species are distributed across two or more biogeographic provinces, shifts in genotype frequencies often align with biogeographic boundaries, providing intraspecific concordance between biogeography and phylogeography. Here, we provide examples of comparative phylogeography from (i) tropical seas that host the highest marine biodiversity, (ii) temperate seas with high productivity but volatile coastlines, (iii) migratory marine fauna, and (iv) plankton that are the most abundant eukaryotes on earth. Tropical and temperate zones both show impacts of glacial cycles, the former primarily through changing sea levels, and the latter through coastal habitat disruption. The general concordance between biogeography and phylogeography indicates that the population-level genetic divergences observed between provinces are a starting point for macroevolutionary divergences between species. However, isolation between provinces does not account for all marine biodiversity; the remainder arises through alternative pathways, such as ecological speciation and parapatric (semiisolated) divergences within provinces and biodiversity hotspots.

Genetic structure of the crown-of-thorns seastar in the Pacific Ocean, with focus on Guam

Population outbreaks of the corallivorous crown-of-thorns seastar (COTS), Acanthaster 'planci' L., are among the most important biological disturbances of tropical coral reefs. Over the past 50 years, several devastating outbreaks have been documented around Guam, an island in the western Pacific Ocean. Previous analyses have shown that in the Pacific Ocean, COTS larval dispersal may be geographically restricted to certain regions. Here, we assess the genetic structure of Pacific COTS populations and compared samples from around Guam with a number of distant localities in the Pacific Ocean, and focused on determining the degree of genetic structure among populations previously considered to be isolated. Using microsatellites, we document substantial genetic structure between 14 localities from different geographical regions in the Pacific Ocean. Populations from the 14 locations sampled were found to be structured in three significantly differentiated groups: (1) all locations immediately around Guam, as well as Kingman Reef and Swains Island; (2) Japan, Philippines, GBR and Vanuatu; and (3) Johnston Atoll, which was significantly different from all other localities. The lack of genetic differentiation between Guam and extremely distant populations from Kingman Reef and Swains Island suggests potential long-distance dispersal of COTS in the Pacific.

New view of population genetics of zooplankton: RAD-seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium-to-small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next-generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site-associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High-density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.

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.

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.

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.

Mitochondrial DNA reveals genetic structuring of Pinna nobilis across the Mediterranean Sea

Pinna nobilis is the largest endemic Mediterranean marine bivalve. During past centuries, various human activities have promoted the regression of its populations. As a consequence of stringent standards of protection, demographic expansions are currently reported in many sites. The aim of this study was to provide the first large broad-scale insight into the genetic variability of P. nobilis in the area that encompasses the western Mediterranean, Ionian Sea, and Adriatic Sea marine ecoregions. To accomplish this objective twenty-five populations from this area were surveyed using two mitochondrial DNA markers (COI and 16S). Our dataset was then merged with those obtained in other studies for the Aegean and Tunisian populations (eastern Mediterranean), and statistical analyses (Bayesian model-based clustering, median-joining network, AMOVA, mismatch distribution, Tajima's and Fu's neutrality tests and Bayesian skyline plots) were performed. The results revealed genetic divergence among three distinguishable areas: (1) western Mediterranean and Ionian Sea; (2) Adriatic Sea; and (3) Aegean Sea and Tunisian coastal areas. From a conservational point of view, populations from the three genetically divergent groups found may be considered as different management units.

Barriers to gene flow in the marine environment: insights from two common intertidal limpet species of the Atlantic and Mediterranean

Knowledge of the scale of dispersal and the mechanisms governing gene flow in marine environments remains fragmentary despite being essential for understanding evolution of marine biota and to design management plans. We use the limpets Patella ulyssiponensis and Patella rustica as models for identifying factors affecting gene flow in marine organisms across the North-East Atlantic and the Mediterranean Sea. A set of allozyme loci and a fragment of the mitochondrial gene cytochrome C oxidase subunit I were screened for genetic variation through starch gel electrophoresis and DNA sequencing, respectively. An approach combining clustering algorithms with clinal analyses was used to test for the existence of barriers to gene flow and estimate their geographic location and abruptness. Sharp breaks in the genetic composition of individuals were observed in the transitions between the Atlantic and the Mediterranean and across southern Italian shores. An additional break within the Atlantic cluster separates samples from the Alboran Sea and Atlantic African shores from those of the Iberian Atlantic shores. The geographic congruence of the genetic breaks detected in these two limpet species strongly supports the existence of transpecific barriers to gene flow in the Mediterranean Sea and Northeastern Atlantic. This leads to testable hypotheses regarding factors restricting gene flow across the study area.

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.

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.

Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern

We assessed the genetic structure of populations of the widely distributed sea cucumber Holothuria (Holothuria) mammata Grube, 1840, and investigated the effects of marine barriers to gene flow and historical processes. Several potential genetic breaks were considered, which would separate the Atlantic and Mediterranean basins, the isolated Macaronesian Islands from the other locations analysed, and the Western Mediterranean and Aegean Sea (Eastern Mediterranean). We analysed mitochondrial 16S and COI gene sequences from 177 individuals from four Atlantic locations and four Mediterranean locations. Haplotype diversity was high (H=0.9307 for 16S and 0.9203 for COI), and the haplotypes were closely related (π=0.0058 for 16S and 0.0071 for COI). The lowest genetic diversities were found in the Aegean Sea population. Our results showed that the COI gene was more variable and more useful for the detection of population structure than the 16S gene. The distribution of mtDNA haplotypes, the pairwise F(ST) values and the results of exact tests and amova revealed: (i) a significant genetic break between the population in the Aegean Sea and those in the other locations, as supported by both mitochondrial genes, and (ii) weak differentiation of the Canary and Azores Islands from the other populations; however, the populations from the Macaronesian Islands, Algarve and West Mediterranean could be considered to be a panmictic metapopulation. Isolation by distance was not identified in H. (H.) mammata. Historical events behind the observed findings, together with the current oceanographic patterns, were proposed and discussed as the main factors that determine the population structure and genetic signature of H. (H.) mammata.

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.

Limits to gene flow in a cosmopolitan marine planktonic diatom

The role of geographic isolation in marine microbial speciation is hotly debated because of the high dispersal potential and large population sizes of planktonic microorganisms and the apparent lack of strong dispersal barriers in the open sea. Here, we show that gene flow between distant populations of the globally distributed, bloom-forming diatom species Pseudo-nitzschia pungens (clade I) is limited and follows a strong isolation by distance pattern. Furthermore, phylogenetic analysis implies that under appropriate geographic and environmental circumstances, like the pronounced climatic changes in the Pleistocene, population structuring may lead to speciation and hence may play an important role in diversification of marine planktonic microorganisms. A better understanding of the factors that control population structuring is thus essential to reveal the role of allopatric speciation in marine microorganisms.

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.

Mito-nuclear discordance in the degree of population differentiation in a marine goby

An increasing number of phylogeographic studies on marine species shows discordant patterns in the degree of population differentiation between nuclear and mitochondrial markers. To understand better which factors have the potential to cause these patterns of discordance in marine organisms, a population genetic study was realized on the sand goby Pomatoschistus minutus (Pallas 1770; Gobiidae, Teleostei). Sand gobies from eight European locations were genotyped at eight microsatellite markers. Microsatellites confirmed the global phylogeographical pattern of P. minutus observed with mitochondrial DNA (mtDNA) markers and nuclear allozyme markers. Three groups consistent with the mitochondrial lineages were defined (the Mediterranean, Iberian and North Atlantic groups) and indications of a recent founder event in the northern Baltic Sea were found. Nevertheless, differences in the degree of population differentiation between the nuclear and mitochondrial markers were large (global F(ST)-values for microsatellites=0.0121; for allozymes=0.00831; for mtDNA=0.4293). Selection, sex-biased dispersal, homoplasy and a high effective population size are generally accepted as explanations for this mitonuclear discrepancy in the degree of population differentiation. In this study, selection on mtDNA and microsatellites, male-biased dispersal and homoplasy on microsatellite markers are unlikely to be a main cause for this discrepancy. The most likely reason for the discordant pattern is a recent demographical expansion of the sand goby, resulting in high effective population sizes slowing down the differentiation of nuclear DNA.

Complete mitochondrial genome sequences of the three pelagic chaetognaths Sagitta nagae, Sagitta decipiens and Sagitta enflata

The complete nucleotide sequences of the mitochondrial genomes were determined for the three pelagic chaetognaths, Sagitta nagae, Sagitta decipiens, and Sagitta enflata. The mitochondrial genomes of these species which were 11,459, 11,121, and 12,631bp in length, respectively, contained 14 genes (11 protein-coding genes, one transfer RNA gene, and two ribosomal RNA genes), and were found to have lost 23 genes that are present in the typical metazoan mitochondrial genome. The same mitochondrial genome contents have been reported from the benthic chaetognaths belonging to the family Spadellidae, Paraspadella gotoi and Spadella cephaloptera. Within the phylum Chaetognatha, Sagitta and Spadellidae are distantly related, suggesting that the gene loss occurred in the ancestral species of the phylum. The gene orders of the three Sagitta species are markedly different from those of the other non-Chaetognatha metazoans. In contrast to the region with frequent gene rearrangements, no gene rearrangements were observed in the gene cluster encoding COII-III, ND1-3, srRNA, and tRNA(met). Within this conserved gene cluster, gene rearrangements were not observed in the three Sagitta species or between the Sagitta and Spadellidae species. The gene order of this cluster was also assumed to be the ancestral state of the phylum.

Careful with understudied phyla: the case of chaetognath

Background

A recent study by Barthélémy et al. described a set of ribosomal protein (RP) genes extracted from a collection of expressed sequence tags (ESTs) of the chaetognath (arrow worm) Spadella cephaloptera. Three main conclusions were drawn in this paper. First, the authors stated that RP genes present paralogous copies, which have arisen through allopolyploidization. Second, they reported two alternate nucleotide stretches conserved within the 5' untranslated regions (UTR) of multiple ribosomal cDNAs and they suggested that these motifs are involved in the differential transcriptional regulation of paralogous RP genes. Third, they claimed that the phylogenetic position of chaetognaths could not be accurately inferred from a RP dataset because of the persistence of two problems: a long branch attraction (LBA) artefact and a compositional bias.

Results

We reconsider here the results described in Barthélémy et al. and question the evidence on which they are based. We find that their evidence for paralogous copies relies on faulty PCR experiments since they attempted to amplify DNA fragments absent from the genomic template. Our PCR experiments proved that the conserved motifs in 5'UTRs that they targeted in their amplifications are added post-transcriptionally by a trans-splicing mechanism. Then, we showed that the lack of phylogenetic resolution observed by these authors is due to limited taxon sampling and not to LBA or to compositional bias. A ribosomal protein dataset thus fully supports the position of chaetognaths as sister group of all other protostomes. This reinterpretation demonstrates that the statements of Barthélémy et al. should be taken with caution because they rely on inaccurate evidence.

Conclusion

The genomic study of an unconventional model organism is a meaningful approach to understand the evolution of animals. However, the previous study came to incorrect conclusions on the basis of experiments that omitted validation procedures.

Chaetognath transcriptome reveals ancestral and unique features among bilaterians

BACKGROUND

The chaetognaths (arrow worms) have puzzled zoologists for years because of their astonishing morphological and developmental characteristics. Despite their deuterostome-like development, phylogenomic studies recently positioned the chaetognath phylum in protostomes, most likely in an early branching. This key phylogenetic position and the peculiar characteristics of chaetognaths prompted further investigation of their genomic features.

RESULTS

Transcriptomic and genomic data were collected from the chaetognath Spadella cephaloptera through the sequencing of expressed sequence tags and genomic bacterial artificial chromosome clones. Transcript comparisons at various taxonomic scales emphasized the conservation of a core gene set and phylogenomic analysis confirmed the basal position of chaetognaths among protostomes. A detailed survey of transcript diversity and individual genotyping revealed a past genome duplication event in the chaetognath lineage, which was, surprisingly, followed by a high retention rate of duplicated genes. Moreover, striking genetic heterogeneity was detected within the sampled population at the nuclear and mitochondrial levels but cannot be explained by cryptic speciation. Finally, we found evidence for trans-splicing maturation of transcripts through splice-leader addition in the chaetognath phylum and we further report that this processing is associated with operonic transcription.

CONCLUSION

These findings reveal both shared ancestral and unique derived characteristics of the chaetognath genome, which suggests that this genome is likely the product of a very original evolutionary history. These features promote chaetognaths as a pivotal model for comparative genomics, which could provide new clues for the investigation of the evolution of animal genomes.

Transcriptome Analysis of ESTs from a Chaetognath Reveals a Deep-Branching Clade of Retrovirus-Like Retrotransposons

Chaetognaths constitute a small marine phylum exhibiting several characteristic which are highly unusual in animal genomes, including two classes of both rRNA and protein ribosomal genes. As in this phylum presence of retrovirus-like elements has never been documented, analysis of a published expressed sequence tag (EST) collection of the chaetognath Spadella cephaloptera has been made. Twelve sequences representing transcript sections of reverse transcriptase domain of active retrotransposons were isolated from~11,000 ESTs. Five of them are originated from Gypsy retrovirus-like elements, whereas the other are transcripts from a Bel-Pao LTR-retrotransposon, a Penelope-like element and LINE retrotransposons. Moreover, a part of a putative integrase has also been found. Phylogenetic analyses suggest a deep-branching clade of the retrovirus-like elements, which is in agreement with the probably Cambrian origin of the phylum. Moreover, retrotransposons have not been found in telomeric-like transcripts which are probably constituted by both vertebrate and arthropod canonical repeats.

Genetic structuring of the temperate gorgonian coral (Corallium rubrum) across the western Mediterranean Sea revealed by microsatellites and nuclear sequences

In the past decades, anthropogenic disturbance has increased in marine costal habitats, leading to dramatic shifts in population size structure of various marine species. In the temperate region, the gorgonian coral (Corallium rubrum) is one of the major disturbed species, mostly owing to the exploitation of its red skeleton for jewellery purposes. Red coral is a gonocoric species inhabiting subtidal rocky habitats in the Mediterranean and neighbouring Atlantic coasts. In order to investigate the spatial genetic structuring of C. rubrum, five microsatellite markers and the ribosomal internal transcribed spacer 1 (ITS-1) sequences were analysed in 11 samples from the northwestern Mediterranean Sea. Phylogenetic reconstructions obtained from ITS-1 sequences analysis showed that samples from Minorca and Sardinia islands were the most divergent while the three samples from the Tuscan archipelago (Argentario, Giannutri and Elba) appeared genetically homogeneous. Both markers revealed a strong structuring over large spatial scales (though to a different extent) with no isolation by distance pattern. Microsatellite-based F(ST) estimates appeared much smaller than ITS-based estimates and not significantly correlated, likely due to microsatellite allele size homoplasy typical of these highly polymorphic loci. Our study shows that the absence of clear patterns of genetic structuring over large spatial scales together with strong genetic structuring should be interpreted with caution because such patterns may hide underlying small-scale genetic structuring. Our results further confirm that effective larval dispersal in red coral is highly restricted in the northwestern Mediterranean Sea, suggesting that an increase of anthropogenic disturbance could aggravate the disappearance of red coral, not only along the Mediterranean coasts but also, and with more intensity, in the main Mediterranean islands.

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.

A molecular assessment of northeast Pacific Alaria species (Laminariales, Phaeophyceae) with reference to the utility of DNA barcoding

Despite their relatively complex morphologies, species in the genus Alaria Greville are notoriously difficult to identify with certainty. Morphological characters, often influenced by environmental factors, make individuals in similar habitats artificially appear related. Species identification would, therefore, benefit greatly from the application of molecular tools. We applied DNA barcoding, using the 5' end of the cytochrome c oxidase I (coxI-5') gene from the mitochondrial genome, to define species limits and relationships in northeast Pacific populations of Alaria. This emerging technique is being employed to catalogue species diversity worldwide, particularly among animals, and it has been shown to be sensitive enough to discriminate between closely related species. However, the utility of this marker for identifying or categorizing the majority of life remains unclear. We compared the resolution obtained with this marker to two other molecular systems commonly used in algal research: the nuclear internal transcribed spacer (ITS) of the ribosomal cistron, and the plastid Rubisco operon spacer (rbcSp). In agreement with previous results, Alaria fistulosa Postels & Ruprecht, with its distinct morphological, ecological and molecular features, stands apart from the other species in the genus and we establish Druehlia gen. nov. to accommodate it. For the remaining isolates, distinct mitochondrial haplotypes resolved with the barcode data indicate a period of genetic isolation for at least three incipient species in the northeast Pacific, whereas unexpected levels and patterns of ITS variation, as well as the extreme morphological plasticity found among these isolates, have most probably resulted from a recent collapse in species barriers. The cloning of ITS amplicons revealed multiple ITS copies in several individuals, further supporting this hypothesis.