Highly structured populations of copepods at risk to deep-sea mining: Integration of genomic data with demogenetic and biophysical modelling

Copepoda is the most abundant taxon in deep-sea hydrothermal vents, where hard substrate is available. Despite the increasing interest in seafloor massive sulphides exploitation, there have been no population genomic studies conducted on vent meiofauna, which are known to contribute over 50% to metazoan biodiversity at vents. To bridge this knowledge gap, restriction-site-associated DNA sequencing, specifically 2b-RADseq, was used to retrieve thousands of genome-wide single-nucleotide polymorphisms (SNPs) from abundant populations of the vent-obligate copepod Stygiopontius lauensis from the Lau Basin. SNPs were used to investigate population structure, demographic histories and genotype-environment associations at a basin scale. Genetic analyses also helped to evaluate the suitability of tailored larval dispersal models and the parameterization of life-history traits that better fit the population patterns observed in the genomic dataset for the target organism. Highly structured populations were observed on both spatial and temporal scales, with divergence of populations between the north, mid, and south of the basin estimated to have occurred after the creation of the major transform fault dividing the Australian and the Niuafo'ou tectonic plate (350 kya), with relatively recent secondary contact events (<20 kya). Larval dispersal models were able to predict the high levels of structure and the highly asymmetric northward low-level gene flow observed in the genomic data. These results differ from most studies conducted on megafauna in the region, elucidating the need to incorporate smaller size when considering site prospecting for deep-sea exploitation of seafloor massive sulphides, and the creation of area-based management tools to protect areas at risk of local extinction, should mining occur.

Contrasted phylogeographic patterns of hydrothermal vent gastropods along South West Pacific: Woodlark Basin, a possible contact zone and/or stepping-stone

Understanding drivers of biodiversity patterns is essential to evaluate the potential impact of deep-sea mining on ecosystems resilience. While the South West Pacific forms an independent biogeographic province for hydrothermal vent fauna, different degrees of connectivity among basins were previously reported for a variety of species depending on their ability to disperse. In this study, we compared phylogeographic patterns of several vent gastropods across South West Pacific back-arc basins and the newly-discovered La Scala site on the Woodlark Ridge by analysing their genetic divergence using a barcoding approach. We focused on six genera of vent gastropods widely distributed in the region: Lepetodrilus, Symmetromphalus, Lamellomphalus, Shinkailepas, Desbruyeresia and Provanna. A wide-range sampling was conducted at different vent fields across the Futuna Volcanic Arc, the Manus, Woodlark, North Fiji, and Lau Basins, during the CHUBACARC cruise in 2019. The Cox1-based genetic structure of geographic populations was examined for each taxon to delineate putative cryptic species and assess potential barriers or contact zones between basins. Results showed contrasted phylogeographic patterns among species, even between closely related species. While some species are widely distributed across basins (i.e. Shinkailepas tollmanni, Desbruyeresia melanioides and Lamellomphalus) without evidence of strong barriers to gene flow, others are restricted to one (i.e. Shinkailepas tufari complex of cryptic species, Desbruyeresia cancellata and D. costata). Other species showed intermediate patterns of isolation with different lineages separating the Manus Basin from the Lau/North Fiji Basins (i.e. Lepetodrilus schrolli, Provanna and Symmetromphalus spp.). Individuals from the Woodlark Basin were either endemic to this area (though possibly representing intermediate OTUs between the Manus Basin and the other eastern basins populations) or, coming into contact from these basins, highlighting the stepping-stone role of the Woodlark Basin in the dispersal of the South West Pacific vent fauna. Results are discussed according to the dispersal ability of species and the geological history of the South West Pacific.

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 harbor, 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. Demo-genetic 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 neighboring 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 datasets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark

With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.