Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna.

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Anammox bacteria drive fixed nitrogen loss in hadal trench sediments.

Benthic N2 production by microbial denitrification and anammox is the largest sink for fixed nitrogen in the oceans. Most N2 production occurs on the continental shelves, where a high flux of reactive organic matter fuels the depletion of nitrate close to the sediment surface. By contrast, N2 production rates in abyssal sediments are low due to low inputs of reactive organics, and nitrogen transformations are dominated by aerobic nitrification and the release of nitrate to the bottom water. Here, we demonstrate that this trend is reversed in the deepest parts of the oceans, the hadal trenches, where focusing of reactive organic matter enhances benthic microbial activity. Thus, at ∼8-km depth in the Atacama Trench, underlying productive surface waters, nitrate is depleted within a few centimeters of the sediment surface, N2 production rates reach those reported from some continental margin sites, and fixed nitrogen loss is mainly conveyed by anammox bacteria. These bacteria are closely related to those known from shallow oxygen minimum zone waters, and comparison of activities measured in the laboratory and in situ suggest they are piezotolerant. Even the Kermadec Trench, underlying oligotrophic surface waters, exhibits substantial fixed N removal. Our results underline the role of hadal sediments as hot spots of deep-sea biological activity, revealing a fully functional benthic nitrogen cycle at high hydrostatic pressure and pointing to hadal sediments as a previously unexplored niche for anaerobic microbial ecology and diagenesis.

Stepping up to genome scan allows stock differentiation in the worldwide distributed blue shark Prionace glauca.

The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome-wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo-west Pacific samples (Southern population). Significant pairwise FST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies.

Subtle limits to connectivity revealed by outlier loci within two divergent metapopulations of the deep-sea hydrothermal gastropod Ifremeria nautilei.

Hydrothermal vents form archipelagos of ephemeral deep-sea habitats that raise interesting questions about the evolution and dynamics of the associated endemic fauna, constantly subject to extinction-recolonization processes. These metal-rich environments are coveted for the mineral resources they harbour, thus raising recent conservation concerns. The evolutionary fate and demographic resilience of hydrothermal species strongly depend on the degree of connectivity among and within their fragmented metapopulations. In the deep sea, however, assessing connectivity is difficult and usually requires indirect genetic approaches. Improved detection of fine-scale genetic connectivity is now possible based on genome-wide screening for genetic differentiation. Here, we explored population connectivity in the hydrothermal vent snail Ifremeria nautilei across its species range encompassing five distinct back-arc basins in the Southwest Pacific. The global analysis, based on 10,570 single nucleotide polymorphism (SNP) markers derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), depicted two semi-isolated and homogeneous genetic clusters. Demogenetic modeling suggests that these two groups began to diverge about 70,000 generations ago, but continue to exhibit weak and slightly asymmetrical gene flow. Furthermore, a careful analysis of outlier loci showed subtle limitations to connectivity between neighbouring basins within both groups. This finding indicates that migration is not strong enough to totally counterbalance drift or local selection, hence questioning the potential for demographic resilience at this latter geographical scale. These results illustrate the potential of large genomic data sets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

The Combined Effect of Haplodiplonty and Partial Clonality on Genotypic and Genetic Diversity in a Finite Mutating Population.

Partial clonality is known to affect the genetic composition and evolutionary trajectory of diplontic (single, free-living diploid stage) populations. However, many partially clonal eukaryotes exhibit life cycles in which somatic development occurs in both haploid and diploid individuals (haplodiplontic life cycles). Here, we studied how haplodiplontic life cycles and partial clonality structurally constrain, as immutable parameters, the reshuffling of genetic diversity and its dynamics in populations over generations. We assessed the distribution of common population genetic indices at different proportions of haploids, rates of clonality, mutation rates, and sampling efforts. Our results showed that haplodiplontic life cycles alone in finite populations affect effective population sizes and the ranges of distributions of population genetic indices. With nonoverlapping generations, haplodiplonty allowed the evolution of 2 temporal genetic pools that may diverge in sympatry due to genetic drift under full sexuality and clonality. Partial clonality in these life cycles acted as a homogenizing force between those 2 pools. Moreover, the combined effects of proportion of haploids, rate of clonality, and the relative strength of mutation versus genetic drift impacts the distributions of population genetics indices, rendering it difficult to transpose and use knowledge accumulated from diplontic or haplontic species. Finally, we conclude by providing recommendations for sampling and analyzing the population genetics of partially clonal haplodiplontic taxa.

New insights into the population genetics of partially clonal organisms: When seagrass data meet theoretical expectations.

Seagrass meadows are among the most important coastal ecosystems in terms of both spatial extent and ecosystem services, but they are also declining worldwide. Understanding the drivers of seagrass meadow dynamics is essential for designing sound management, conservation and restoration strategies. However, poor knowledge of the effect of clonality on the population genetics of natural populations severely limits our understanding of the dynamics and connectivity of meadows. Recent modelling approaches have described the expected distributions of genotypic and genetic descriptors under increasing clonal rates, which may help us better understand and interpret population genetics data obtained for partial asexuals. Here, in the light of these recent theoretical developments, we revisited population genetics data for 165 meadows of four seagrass species. Contrasting shoot lifespan and rhizome turnover led to the prediction that the influence of asexual reproduction would increase along a gradient from Zostera noltii to Zostera marina, Cymodocea nodosa and Posidonia oceanica, with increasing departure from Hardy-Weinberg equilibrium (Fis ), mostly towards heterozygote excess, and decreasing genotypic richness (R). This meta-analysis provides a nested validation of this hypothesis at both the species and meadow scales through a significant relationship between Fis and R within each species. By empirically demonstrating the theoretical expectations derived from recent modelling approaches, this work calls for the use of Hardy-Weinberg equilibrium (Fis ) rather than only the strongly sampling-sensitive R to assess the importance of clonal reproduction (c), at least when the impact of selfing on Fis can be neglected. The results also emphasize the need to revise our appraisal of the extent of clonality and its influence on the dynamics, connectivity and evolutionary trajectory of partial asexuals in general, including in seagrass meadows, to develop the most accurate management strategies.

Climate-induced changes in the suitable habitat of cold-water corals and commercially important deep-sea fishes in the North Atlantic.

The deep sea plays a critical role in global climate regulation through uptake and storage of heat and carbon dioxide. However, this regulating service causes warming, acidification and deoxygenation of deep waters, leading to decreased food availability at the seafloor. These changes and their projections are likely to affect productivity, biodiversity and distributions of deep-sea fauna, thereby compromising key ecosystem services. Understanding how climate change can lead to shifts in deep-sea species distributions is critically important in developing management measures. We used environmental niche modelling along with the best available species occurrence data and environmental parameters to model habitat suitability for key cold-water coral and commercially important deep-sea fish species under present-day (1951-2000) environmental conditions and to project changes under severe, high emissions future (2081-2100) climate projections (RCP8.5 scenario) for the North Atlantic Ocean. Our models projected a decrease of 28%-100% in suitable habitat for cold-water corals and a shift in suitable habitat for deep-sea fishes of 2.0°-9.9° towards higher latitudes. The largest reductions in suitable habitat were projected for the scleractinian coral Lophelia pertusa and the octocoral Paragorgia arborea, with declines of at least 79% and 99% respectively. We projected the expansion of suitable habitat by 2100 only for the fishes Helicolenus dactylopterus and Sebastes mentella (20%-30%), mostly through northern latitudinal range expansion. Our results projected limited climate refugia locations in the North Atlantic by 2100 for scleractinian corals (30%-42% of present-day suitable habitat), even smaller refugia locations for the octocorals Acanella arbuscula and Acanthogorgia armata (6%-14%), and almost no refugia for P. arborea. Our results emphasize the need to understand how anticipated climate change will affect the distribution of deep-sea species including commercially important fishes and foundation species, and highlight the importance of identifying and preserving climate refugia for a range of area-based planning and management tools.

Rare sex or out of reach equilibrium? The dynamics of F IS in partially clonal organisms.

BACKGROUND: Partially clonal organisms are very common in nature, yet the influence of partial asexuality on the temporal dynamics of genetic diversity remains poorly understood. Mathematical models accounting for clonality predict deviations only for extremely rare sex and only towards mean inbreeding coefficient [Formula: see text]. Yet in partially clonal species, both F IS < 0 and F IS > 0 are frequently observed also in populations where there is evidence for a significant amount of sexual reproduction. Here, we studied the joint effects of partial clonality, mutation and genetic drift with a state-and-time discrete Markov chain model to describe the dynamics of F IS over time under increasing rates of clonality. RESULTS: Results of the mathematical model and simulations show that partial clonality slows down the asymptotic convergence to F IS = 0. Thus, although clonality alone does not lead to departures from Hardy-Weinberg expectations once reached the final equilibrium state, both negative and positive F IS values can arise transiently even at intermediate rates of clonality. More importantly, such "transient" departures from Hardy Weinberg proportions may last long as clonality tunes up the temporal variation of F IS and reduces its rate of change over time, leading to a hyperbolic increase of the maximal time needed to reach the final mean [Formula: see text] value expected at equilibrium. CONCLUSION: Our results argue for a dynamical interpretation of F IS in clonal populations. Negative values cannot be interpreted as unequivocal evidence for extremely scarce sex but also as intermediate rates of clonality in finite populations. Complementary observations (e.g. frequency distribution of multiloci genotypes, population history) or time series data may help to discriminate between different possible conclusions on the extent of clonality when mean [Formula: see text] values deviating from zero and/or a large variation of F IS over loci are observed.

Species are hypotheses: avoid connectivity assessments based on pillars of sand.

Connectivity among populations determines the dynamics and evolution of populations, and its assessment is essential in ecology in general and in conservation biology in particular. The robust basis of any ecological study is the accurate delimitation of evolutionary units, such as populations, metapopulations and species. Yet a disconnect still persists between the work of taxonomists describing species as working hypotheses and the use of species delimitation by molecular ecologists interested in describing patterns of gene flow. This problem is particularly acute in the marine environment where the inventory of biodiversity is relatively delayed, while for the past two decades, molecular studies have shown a high prevalence of cryptic species. In this study, we illustrate, based on marine case studies, how the failure to recognize boundaries of evolutionary-relevant unit leads to heavily biased estimates of connectivity. We review the conceptual framework within which species delimitation can be formalized as falsifiable hypotheses and show how connectivity studies can feed integrative taxonomic work and vice versa. Finally, we suggest strategies for spatial, temporal and phylogenetic sampling to reduce the probability of inadequately delimiting evolutionary units when engaging in connectivity studies.

Disentangling the Influence of Mutation and Migration in Clonal Seagrasses Using the Genetic Diversity Spectrum for Microsatellites.

The recurrent lack of isolation by distance reported at regional scale in seagrass species was recently suggested to stem from stochastic events of large-scale dispersal. We explored the usefulness of phylogenetic information contained in microsatellite loci to test this hypothesis by using the Genetic Diversity Spectrum (GDS) on databases containing, respectively, 7 and 9 microsatellites genotypes for 1541 sampling units of Posidonia oceanica and 1647 of Cymodocea nodosa The simultaneous increase of microsatellite and geographic distances that emerges reveals a coherent pattern of isolation by distance in contrast to the chaotic pattern previously described using allele frequencies, in particular, for the long-lived P. oceanica These results suggest that the lack of isolation by distance, rather than the resulting from rare events of large-scale dispersal, reflects at least for some species a stronger influence of mutation over migration at the scale of the distribution range. The global distribution of genetic polymorphism may, therefore, result predominantly from ancient events of step-by-step (re)colonization followed by local recruitment and clonal growth, rather than contemporary gene flow. The analysis of GDS appears useful to unravel the evolutionary forces influencing the dynamics and evolution at distinct temporal and spatial scales by accounting for phylogenetic information borne by microsatellites, under an appropriate mutation model. This finding adds nuance to the generalization of the influence of large-scale dispersal on the dynamics of seagrasses.

Metagenome-assembled genomes of deep-sea sediments: changes in microbial functional potential lag behind redox transitions.

Hadal sediments are hotspots of microbial activity in the deep sea and exhibit strong biogeochemical gradients. But although these gradients are widely assumed to exert selective forces on hadal microbial communities, the actual relationship between biogeochemistry, functional traits, and microbial community structure remains poorly understood. We tested whether the biogeochemical conditions in hadal sediments select for microbes based on their genomic capacity for respiration and carbohydrate utilization via a metagenomic analysis of over 153 samples from the Atacama Trench region (max. depth = 8085 m). The obtained 1357 non-redundant microbial genomes were affiliated with about one-third of all known microbial phyla, with more than half belonging to unknown genera. This indicated that the capability to withstand extreme hydrostatic pressure is a phylogenetically widespread trait and that hadal sediments are inhabited by diverse microbial lineages. Although community composition changed gradually over sediment depth, these changes were not driven by selection for respiratory or carbohydrate degradation capability in the oxic and nitrogenous zones, except in the case of anammox bacteria and nitrifying archaea. However, selection based on respiration and carbohydrate degradation capacity did structure the communities of the ferruginous zone, where aerobic and nitrogen respiring microbes declined exponentially (half-life = 125-419 years) and were replaced by subsurface communities. These results highlight a delayed response of microbial community composition to selective pressure imposed by redox zonation and indicated that gradual changes in microbial composition are shaped by the high-resilience and slow growth of microbes in the seafloor.

High connectivity across the fragmented chemosynthetic ecosystems of the deep Atlantic Equatorial Belt: efficient dispersal mechanisms or questionable endemism?

Chemosynthetic ecosystems are distributed worldwide in fragmented habitats harbouring seemingly highly specialized communities. Yet, shared taxa have been reported from highly distant chemosynthetic communities. These habitats are distributed in distinct biogeographical regions, one of these being the so-called Atlantic Equatorial Belt (AEB). Here, we combined genetic data (COI) from several taxa to assess the possible existence of cryptic or synonymous species and to detect the possible occurrence of contemporary gene flow among populations of chemosynthetic species located on both sides of the Atlantic. Several Evolutionary Significant Units (ESUs) of Alvinocarididae shrimp and Vesicomyidae bivalves were found to be shared across seeps of the AEB. Some were also common to hydrothermal vent communities of the Mid-Atlantic Ridge (MAR), encompassing taxa morphologically described as distinct species or even genera. The hypothesis of current or very recent large-scale gene flow among seeps and vents was supported by microsatellite analysis of the shrimp species Alvinocaris muricola/Alvinocaris markensis across the AEB and MAR. Two nonmutually exclusive hypotheses may explain these findings. The dispersion of larvae or adults following strong deep-sea currents, possibly combined with biochemical cues influencing the duration of larval development and timing of metamorphosis, may result in large-scale effective migration among distant spots scattered on the oceanic seafloor. Alternatively, these results may arise from the prevailing lack of knowledge on the ocean seabed, apart from emblematic ecosystems (chemosynthetic ecosystems, coral reefs or seamounts), where the widespread classification of endemism associated with many chemosynthetic taxa might hide wider distributions in overlooked parts of the deep sea.

Restriction to large-scale gene flow vs. regional panmixia among cold seep Escarpia spp. (Polychaeta, Siboglinidae).

The history of colonization and dispersal in fauna distributed among deep-sea chemosynthetic ecosystems remains enigmatic and poorly understood because of an inability to mark and track individuals. A combination of molecular, morphological and environmental data improves understanding of spatial and temporal scales at which panmixia, disruption of gene flow or even speciation may occur. Vestimentiferan tubeworms of the genus Escarpia are important components of deep -sea cold seep ecosystems, as they provide long-term habitat for many other taxa. Three species of Escarpia, Escarpia spicata [Gulf of California (GoC)], Escarpia laminata [Gulf of Mexico (GoM)] and Escarpia southwardae (West African Cold Seeps), have been described based on morphology, but are not discriminated through the use of mitochondrial markers (cytochrome oxidase subunit 1; large ribosomal subunit rDNA, 16S; cytochrome b). Here, we also sequenced the exon-primed intron-crossing Haemoglobin subunit B2 intron and genotyped 28 microsatellites to (i) determine the level of genetic differentiation, if any, among the three geographically separated entities and (ii) identify possible population structure at the regional scale within the GoM and West Africa. Results at the global scale support the occurrence of three genetically distinct groups. At the regional scale among eight sampling sites of E. laminata (n = 129) and among three sampling sites of E. southwardae (n = 80), no population structure was detected. These findings suggest that despite the patchiness and isolation of seep habitats, connectivity is high on regional scales.

Biogeography revisited with network theory: retracing the history of hydrothermal vent communities.

Defining biogeographic provinces to understand the history and evolution of communities associated with a given kind of ecosystem is challenging and usually requires a priori assumptions to be made. We applied network theory, a holistic and exploratory method, to the most complete database of faunal distribution available on oceanic hydrothermal vents, environments which support fragmented and unstable ecosystems, to infer the processes driving their worldwide biogeography. Besides the identification of robust provinces, the network topology allowed us to identify preferential pathways that had hitherto been overlooked. These pathways are consistent with the previously proposed hypothesis of a role of plate tectonics in the biogeographical history of hydrothermal vent communities. A possible ancestral position of the Western Pacific is also suggested for the first time. Finally, this work provides an innovative example of the potential of network tools to unravel the biogeographic history of faunal assemblages and to supply comprehensive information for the conservation and management of biodiversity.

Entangled effects of allelic and clonal (genotypic) richness in the resistance and resilience of experimental populations of the seagrass Zostera noltii to diatom invasion.

BACKGROUND: The relationship between species diversity and components of ecosystem stability has been extensively studied, whilst the influence of the genetic component of biodiversity remains poorly understood. Here we manipulated both genotypic and allelic richness of the seagrass Zostera noltii, in order to explore their respective influences on the resistance of the experimental population to stress. Thus far intra-specific diversity was seldom taken into account in management plans, and restoration actions showed very low success. Information is therefore needed to understand the factors affecting resistance and resilience of populations. RESULTS: Our results show a positive influence of both allelic and genotypic richness on the resistance of meadows to environmental perturbations. They also show that at the low genotypic (i.e. clonal) richness levels used in prior experimental approaches, the effects of genotypic and allelic richness could not be disentangled and allelic richness was a likely hidden treatment explaining at least part of the effects hitherto attributed to genotypic richness. CONCLUSIONS: Altogether, these results emphasize the need to acknowledge and take into account the interdependency of both genotypic and allelic richness in experimental designs attempting to estimate their importance alone or in combination. A positive influence of allelic richness on resistance to perturbations, and of allelic richness combined with genotypic richness on the recovery (resilience) of the experimental populations is supported by differential mortality. These results, on the key species structuring of one of the most threatened coastal ecosystem worldwide, seagrass meadows, support the need to better take into account the distinct compartments of clonal and genetic diversity in management strategies, and in possible restoration plans in the future.

What lies underneath: conserving the oceans' genetic resources.

The marine realm represents 70% of the surface of the biosphere and contains a rich variety of organisms, including more than 34 of the 36 living phyla, some of which are only found in the oceans. The number of marine species used by humans is growing at unprecedented rates, including the rapid domestication of marine species for aquaculture and the discovery of natural products and genes of medical and biotechnological interest in marine biota. The rapid growth in the human appropriation of marine genetic resources (MGRs), with over 18,000 natural products and 4,900 patents associated with genes of marine organisms, with the latter growing at 12% per year, demonstrates that the use of MGRs is no longer a vision but a growing source of biotechnological and business opportunities. The diversification of the use of marine living resources by humans calls for an urgent revision of the goals and policies of marine protected areas, to include the protection of MGRs and address emerging issues like biopiracy or benefit sharing. Specific challenges are the protection of these valuable resources in international waters, where no universally accepted legal framework exists to protect and regulate the exploitation of MGRs, and the unresolved issues on patenting components of marine life. Implementing steps toward the protection of MGRs is essential to ensure their sustainable use and to support the flow of future findings of medical and biotechnological interest.

Distribution and genomic variation of ammonia-oxidizing archaea in abyssal and hadal surface sediments.

Ammonia-oxidizing archaea of the phylum Thaumarchaeota play a central role in the biogeochemical cycling of nitrogen in benthic sediments, at the interface between pelagic and subsurface ecosystems. However, our understanding of their niche separation and of the processes controlling their population structure in hadal and abyssal surface sediments is still limited. Here, we reconstructed 47 AOA metagenome-assembled genomes (MAGs) from surface sediments of the Atacama and Kermadec trench systems. They formed deep-sea-specific groups within the family Nitrosopumilaceae and were assigned to six amoA gene-based clades. MAGs from different clades had distinct distribution patterns along oxygen-ammonium counter gradients in surface sediments. At the species level, MAGs thus seemed to form different ecotypes and follow deterministic niche-based distributions. In contrast, intraspecific population structure, defined by patterns of Single Nucleotide Variants (SNV), seemed to reflect more complex contributions of both deterministic and stochastic processes. Firstly, the bathymetric range had a strong effect on population structure, with distinct populations in abyssal plains and hadal trenches. Then, hadal populations were clearly separated by trench system, suggesting a strong isolation-by-topography effect, whereas abyssal populations were rather controlled by sediment depth or geographic distances, depending on the clade considered. Interestingly, genetic variability between samples was lowest in sediment layers where the mean MAG coverage was highest, highlighting the importance of selective pressure linked with each AOA clade's ecological niche. Overall, our results show that deep-sea AOA genome distributions seem to follow both deterministic and stochastic processes, depending on the genomic variability scale considered.

Not out of the Mediterranean: Atlantic populations of the gorgonian Paramuricea clavata are a separate sister species under further lineage diversification.

The accurate delimitation of species boundaries in nonbilaterian marine taxa is notoriously difficult, with consequences for many studies in ecology and evolution. Anthozoans are a diverse group of key structural organisms worldwide, but the lack of reliable morphological characters and informative genetic markers hampers our ability to understand species diversification. We investigated population differentiation and species limits in Atlantic (Iberian Peninsula) and Mediterranean lineages of the octocoral genus Paramuricea previously identified as P. clavata. We used a diverse set of molecular markers (microsatellites, RNA-seq derived single-copy orthologues [SCO] and mt-mutS [mitochondrial barcode]) at 49 locations. Clear segregation of Atlantic and Mediterranean lineages was found with all markers. Species-tree estimations based on SCO strongly supported these two clades as distinct, recently diverged sister species with incomplete lineage sorting, P. cf. grayi and P. clavata, respectively. Furthermore, a second putative (or ongoing) speciation event was detected in the Atlantic between two P. cf. grayi color morphotypes (yellow and purple) using SCO and supported by microsatellites. While segregating P. cf. grayi lineages showed considerable geographic structure, dominating circalittoral communities in southern (yellow) and western (purple) Portugal, their occurrence in sympatry at some localities suggests a degree of reproductive isolation. Overall, our results show that previous molecular and morphological studies have underestimated species diversity in Paramuricea occurring in the Iberian Peninsula, which has important implications for conservation planning. Finally, our findings validate the usefulness of phylotranscriptomics for resolving evolutionary relationships in octocorals.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems.

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

Spatial coalescent connectivity through multi-generation dispersal modelling predicts gene flow across marine phyla.

Gene flow governs the contemporary spatial structure and dynamic of populations as well as their long-term evolution. For species that disperse using atmospheric or oceanic flows, biophysical models allow predicting the migratory component of gene flow, which facilitates the interpretation of broad-scale spatial structure inferred from observed allele frequencies among populations. However, frequent mismatches between dispersal estimates and observed genetic diversity prevent an operational synthesis for eco-evolutionary projections. Here we use an extensive compilation of 58 population genetic studies of 47 phylogenetically divergent marine sedentary species over the Mediterranean basin to assess how genetic differentiation is predicted by Isolation-By-Distance, single-generation dispersal and multi-generation dispersal models. Unlike previous approaches, the latter unveil explicit parents-to-offspring links (filial connectivity) and implicit links among siblings from a common ancestor (coalescent connectivity). We find that almost 70 % of observed variance in genetic differentiation is explained by coalescent connectivity over multiple generations, significantly outperforming other models. Our results offer great promises to untangle the eco-evolutionary forces that shape sedentary population structure and to anticipate climate-driven redistributions, altogether improving spatial conservation planning.