Taking the chaos out of genetic patchiness: seascape genetics reveals ecological and oceanographic drivers of genetic patterns in three temperate reef species

Geographic genetic variation in the Coral Hawkfish, Cirrhitichthys oxycephalus (Cirrhitidae), in relation to biogeographic barriers across the Tropical Indo-Pacific

The Tropical Indo-Pacific (TIP) includes about two thirds of the world's tropical oceans and harbors an enormous number of marine species. The distributions of those species within the region is affected by habitat discontinuities and oceanographic features. As well as many smaller ones, the TIP contains seven large recognized biogeographic barriers that separate the Red Sea and Indian Ocean, the Indian from the Pacific Ocean, the central and eastern Pacific, the Hawaiian archipelago, the Marquesas and Easter Islands. We examined the genetic structuring of populations of Cirrhitichthys oxycephalus, a small cryptic species of reef fish, across its geographic range, which spans the longitudinal limits of the TIP. We assessed geographic variation in the mitochondrial cytb gene and the nuclear RAG1 gene, using 166 samples collected in 46 localities from the western to eastern edges of the TIP. Sequences from cytb show three well-structured groups that are separated by large genetic distances (1.58-2.96%): two in the Tropical Eastern Pacific (TEP), one at Clipperton Atoll another occupying the rest of that region and the third that ranges across the remainder of the TIP, from the central Pacific to the Red Sea and South Africa. These results indicate that the ~4,000 km wide Eastern Pacific Barrier between the central and eastern Pacific is an efficient barrier separating the two main groups. Further, the ~950 km of open ocean that isolates Clipperton Atoll from the rest of the TEP is also an effective barrier. Contrary to many other cases, various major and minor barriers from the Central Indo-Pacific to the Red Sea are not effective against dispersal by C. oxycephalus, although this species has not colonized the Hawiian islands and Easter Island. The nuclear gene partially supports the genetic structure evident in cytb, although all haplotypes are geographically mixed.

Oceanographical-driven dispersal and environmental variation explain genetic structure in an upwelling coastal ecosystem

The seascape comprises multiple environmental variables that interact with species biology to determine patterns of spatial genetic variation. The environment imposes spatially variable selective forces together with homogenizing and diverging drivers that facilitate or restrict dispersal, which is a complex, time-dependent process. Understanding how the seascape influences spatial patterns of genetic variation remains elusive, particularly in coastal upwelling systems. Here, we combine genome-wide SNP data, Lagrangian larval dispersal simulated over a hydrodynamic model, and ocean environmental information to quantify the relative contribution of ocean circulation and environmental heterogeneity as drivers of the spatial genetic structure of two congeneric intertidal limpets, Scurria scurra and S. araucana, along the central coast of Chile. We find that a genetic break observed in both limpet species coincides with a break in connectivity shown by the Lagrangian dispersal, suggesting that mean ocean circulation is an important seascape feature, in particular for S. scurra. For S. araucana, environmental variation appears as a better predictor of genetic structure than ocean circulation. Overall, our study shows broad patterns of seascape forcing on genetic diversity and contributes to our understanding of the complex ecological and evolutionary interactions along coastal upwelling systems.

Are Coastal Marine Larvae Dispersed Less Than Would Be Expected by Ocean Currents?

AbstractThe distance that offspring disperse from their parents affects how a species responds to habitat disturbance, climate change, and interspecific interactions. For many benthic species, this dispersal is via planktonic larvae, but the distance these larvae disperse is difficult to observe directly. Dispersal distance has usually been estimated indirectly by combining an observed quantity (e.g., the rate of spread of an invasive organism or genetic similarity between locations) with a model that links that quantity to the dispersal of larvae. The estimates of dispersal distance based on the speed of spread of invasive organisms have led many researchers to conclude that the larvae of most of these organisms disperse much less than would be expected if they were being passively transported by the expected ocean currents (Shanks et al.; Shanks). I argue that the discrepancy is instead caused by the choice of model linking dispersal distance to invasion speed. Their model neglected the impact of life history, population growth, and oceanographic parameters on invasion speed. When dispersal distance is estimated from a more complete model of invasion speed, it is found that larval dispersal distance is not much less than would be expected for larvae drifting in the observed ocean currents.

Ocean connectivity and habitat characteristics predict population genetic structure of seagrass in an extreme tropical setting

Understanding patterns of gene flow and processes driving genetic differentiation is important for a broad range of conservation practices. In marine organisms, genetic differentiation among populations is influenced by a range of spatial, oceanographic, and environmental factors that are attributed to the seascape. The relative influences of these factors may vary in different locations and can be measured using seascape genetic approaches. Here, we applied a seascape genetic approach to populations of the seagrass, Thalassia hemprichii, at a fine spatial scale (~80 km) in the Kimberley coast, western Australia, a complex seascape with strong, multidirectional currents greatly influenced by extreme tidal ranges (up to 11 m, the world's largest tropical tides). We incorporated genetic data from a panel of 16 microsatellite markers, overwater distance, oceanographic data derived from predicted passive dispersal on a 2 km-resolution hydrodynamic model, and habitat characteristics from each meadow sampled. We detected significant spatial genetic structure and asymmetric gene flow, in which meadows 12-14 km apart were less connected than ones 30-50 km apart. This pattern was explained by oceanographic connectivity and differences in habitat characteristics, suggesting a combined scenario of dispersal limitation and facilitation by ocean current with local adaptation. Our findings add to the growing evidence for the key role of seascape attributes in driving spatial patterns of gene flow. Despite the potential for long-distance dispersal, there was significant genetic structuring over small spatial scales implicating dispersal and recruitment bottlenecks and highlighting the importance of implementing local-scale conservation and management measures.

Dispersive currents explain patterns of population connectivity in an ecologically and economically important fish

How to identify the drivers of population connectivity remains a fundamental question in ecology and evolution. Answering this question can be challenging in aquatic environments where dynamic lake and ocean currents coupled with high levels of dispersal and gene flow can decrease the utility of modern population genetic tools. To address this challenge, we used RAD-Seq to genotype 959 yellow perch (Perca flavescens), a species with an ~40-day pelagic larval duration (PLD), collected from 20 sites circumscribing Lake Michigan. We also developed a novel, integrative approach that couples detailed biophysical models with eco-genetic agent-based models to generate "predictive" values of genetic differentiation. By comparing predictive and empirical values of genetic differentiation, we estimated the relative contributions for known drivers of population connectivity (e.g., currents, behavior, PLD). For the main basin populations (i.e., the largest contiguous portion of the lake), we found that high gene flow led to low overall levels of genetic differentiation among populations (F _ST = 0.003). By far the best predictors of genetic differentiation were connectivity matrices that were derived from periods of time when there were strong and highly dispersive currents. Thus, these highly dispersive currents are driving the patterns of population connectivity in the main basin. We also found that populations from the northern and southern main basin are slightly divergent from one another, while those from Green Bay and the main basin are highly divergent (F _ST = 0.11). By integrating biophysical and eco-genetic models with genome-wide data, we illustrate that the drivers of population connectivity can be identified in high gene flow systems.

Population structure and genome-wide evolutionary signatures reveal putative climate-driven habitat change and local adaptation in the large yellow croaker

The large yellow croaker (Larimichthys crocea) is one of the most economically valuable marine fish in China and is a notable species in ecological studies owing to a serious collapse of wild germplasm in the past few decades. The stock division and species distribution, which have important implications for ecological protection, germplasm recovery, and fishery resource management, have been debated since the 1960s. However, it is still uncertain even how many stocks exist in this species. To address this, we evaluated the fine-scale genetic structure of large yellow croaker populations distributed along the eastern and southern Chinese coastline based on 7.64 million SNP markers. Compared with the widely accepted stock boundaries proposed in the 1960s, our results revealed that a climate-driven habitat change probably occurred between the Naozhou (Nanhai) Stock and the Ming-Yuedong (Mindong) Stock. The boundary between these two stocks might have shifted northwards from the Pearl River Estuary to the northern area of the Taiwan Strait, accompanied by highly asymmetric introgression. In addition, we found divergent landscapes of natural selection between the stocks inhabiting northern and southern areas. The northern population exhibited highly agminated signatures of strong natural selection in genes related to developmental processes, whereas moderate and interspersed selective signatures were detected in many immune-related genes in the southern populations. These findings establish the stock status and genome-wide evolutionary landscapes of large yellow croaker, providing a basis for conservation, fisheries management and further evolutionary biology studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00165-2.

Historic and contemporary biogeographic perspectives on range-wide spatial genetic structure in a widespread seagrass

Historical and contemporary processes drive spatial patterns of genetic diversity. These include climate-driven range shifts and gene flow mediated by biogeographical influences on dispersal. Assessments that integrate these drivers are uncommon, but critical for testing biogeographic hypotheses. Here, we characterize intraspecific genetic diversity and spatial structure across the entire distribution of a temperate seagrass to test marine biogeographic concepts for southern Australia. Predictive modeling was used to contrast the current Posidonia australis distribution to its historical distribution during the Last Glacial Maximum (LGM). Spatial genetic structure was estimated for 44 sampled meadows from across the geographical range of the species using nine microsatellite loci. Historical and contemporary distributions were similar, with the exception of the Bass Strait. Genetic clustering was consistent with the three currently recognized biogeographic provinces and largely consistent with the finer-scale IMCRA bioregions. Discrepancies were found within the Flindersian province and southwest IMCRA bioregion, while two regions of admixture coincided with transitional IMCRA bioregions. Clonal diversity was highly variable but positively associated with latitude. Genetic differentiation among meadows was significantly associated with oceanographic distance. Our approach suggests how shared seascape drivers have influenced the capacity of P. australis to effectively track sea level changes associated with natural climate cycles over millennia, and in particular, the recolonization of meadows across the Continental Shelf following the LGM. Genetic structure associated with IMCRA bioregions reflects the presence of stable biogeographic barriers, such as oceanic upwellings. This study highlights the importance of biogeography to infer the role of historical drivers in shaping extant diversity and structure.

Sardines at a junction: Seascape genomics reveals ecological and oceanographic drivers of variation in the NW Mediterranean Sea

By evaluating genetic variation across the entire genome, one can address existing questions in a novel way while raising new ones. The latter includes how different local environments influence adaptive and neutral genomic variation within and among populations, providing insights into local adaptation of natural populations and their responses to global change. Here, under a seascape genomic approach, ddRAD data of 4609 single nucleotide polymorphisms (SNPs) from 398 sardines (Sardina pilchardus) collected in 11 Mediterranean and one Atlantic site were generated. These were used along with oceanographic and ecological information to detect signals of adaptive divergence with gene flow across environmental gradients. The studied sardines constitute two clusters (F_ST  = 0.07), a pattern attributed to outlier loci, highlighting putative local adaptation. The trend in the number of days with sea surface temperature above 19°C, a critical threshold for successful sardine spawning, was crucial at all levels of population structuring with implications on the species' key biological processes. Outliers link candidate SNPs to the region's environmental heterogeneity. Our findings provide evidence for a dynamic equilibrium in which population structure is maintained by physical and ecological factors under the opposing influences of migration and selection. This dynamic in a natural system warrants continuous monitoring under a seascape genomic approach that might benefit from a temporal and more detailed spatial dimension. Our results may contribute to complementary studies aimed at providing deeper insights into the mechanistic processes underlying population structuring. Those are key to understanding and predicting future changes and responses of this highly exploited species in the face of climate change.

Species' attributes predict the relative magnitude of ecological and genetic recovery following mass mortality

Theoretically, species' characteristics should allow estimation of dispersal potential and, in turn, explain levels of population genetic differentiation. However, a mismatch between traits and genetic patterns is often reported for marine species, and interpreted as evidence that life-history traits do not influence dispersal. Here, we couple ecological and genomic methods to test the hypothesis that species with attributes favouring greater dispersal potential-e.g., longer pelagic duration, higher fecundity and larger population size-have greater realized dispersal overall. We used a natural experiment created by a large-scale and multispecies mortality event which created a "clean slate" on which to study recruitment dynamics, thus simplifying a usually complex problem. We surveyed four species of differing dispersal potential to quantify the abundance and distribution of recruits and to genetically assign these recruits to probable parental sources. Species with higher dispersal potential recolonized a broader extent of the impacted range, did so more quickly and recovered more genetic diversity than species with lower dispersal potential. Moreover, populations of taxa with higher dispersal potential exhibited more immigration (71%-92% of recruits) than taxa with lower dispersal potential (17%-44% of recruits). By linking ecological with genomic perspectives, we demonstrate that a suite of interacting life-history and demographic attributes do influence species' realized dispersal and genetic neighbourhoods. To better understand species' resilience and recovery in this time of global change, integrative eco-evolutionary approaches are needed to more rigorously evaluate the effect of dispersal-linked attributes on realized dispersal and population genetic differentiation.

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.

Biophysical models of dispersal contribute to seascape genetic analyses

Dispersal is generally difficult to directly observe. Instead, dispersal is often inferred from genetic markers and biophysical modelling where a correspondence indicates that dispersal routes and barriers explain a significant part of population genetic differentiation. Biophysical models are used for wind-driven dispersal in terrestrial environments and for propagules drifting with ocean currents in the sea. In the ocean, such seascape genetic or seascape genomic studies provide promising tools in applied sciences, as actions within management and conservation rely on an understanding of population structure, genetic diversity and presence of local adaptations, all dependent on dispersal within the metapopulation. Here, we surveyed 87 studies that combine population genetics and biophysical models of dispersal. Our aim was to understand if biophysical dispersal models can generally explain genetic differentiation. Our analysis shows that genetic differentiation and lack of genetic differentiation can often be explained by dispersal, but the realism of the biophysical model, as well as local geomorphology and species biology also play a role. The review supports the use of a combination of both methods, and we discuss our findings in terms of recommendations for future studies and pinpoint areas where further development is necessary, particularly on how to compare both approaches. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double‐digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long‐term (~30‐year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation‐by‐distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

High-throughput SNPs dataset reveal restricted population connectivity of marine gastropod within the narrow distribution range of peripheral oceanic islands

Molecular studies based on the high resolution genetic markers help us to grasp the factor shaping the genetic structure of marine organisms. Ecological factors linking to life history traits have often explained the process of genetic structuring in open and connectable oceanic environments. Besides, population genetic divergence can be affected by fragmented habitat, oceanic current, and past geographical events. In the present study, we demonstrated the genetic differentiation of marine gastropod Monodonta sp. within a narrow range of peripheral oceanic islands, the Ogasawara Islands. Genetic analyses were performed not only with a mitochondrial DNA marker but also with a high-throughput SNPs dataset obtained by ddRAD-seq. The results of the mtDNA analyses did not show genetic divergence among populations, while the SNPs dataset detected population genetic differentiation. Population demographic analyses and gene flow estimation suggested that the genetic structure was formed by sea level fluctuation associated with the past climatic change and regulated by temporal oceanographic conditions. These findings provide important insights into population genetic patterns in open and connectable environments.

Comparing the genetic diversity and population structure of sister marine snails having contrasting habitat specificity

BACKGROUND

To grasp the processes of spatial genetic structuring in open and connectable marine environments is the principal study goal in molecular biological studies. Comparative seascape genetics using multiple species are a powerful approach to understand the physical geographic and oceanographic effects on genetic variation. Besides, species-specific ecological traits such as dispersal abilities and habitat specificity are important factors for spatial genetic structuring.

METHODS AND RESULTS

We focused on the sister marine snail species Tegula kusairo and T. xanthostigma around the Japanese mainland, which have contrasting habitat specificities for wave strength. Tegula kusairo only inhabits sheltered coastal environments, while T. xanthostigma is found mainly on wave-exposed rocky shores facing the open sea. We estimated their genetic diversity indices and levels of population differentiation based on mtDNA. We found that the genetic diversity of T. kusairo was lower than that of T. xanthostigma, while their level of population genetic differentiation was higher than that of T. xanthostigma. Namely, the species specific to weak wave environments had a higher level of population genetic differentiation than the species specific to strong wave action.

CONCLUSION

Ecological traits linked not only to dispersal abilities but also to habitat specificity can influence genetic variation in a pair of closely related sister species distributed in the same seascape.

A paradigm for the recognition of cryptic trematode species in tropical Indo-west Pacific fishes: the problematic genus Preptetos (Trematoda: Lepocreadiidae)

Molecular data have transformed approaches to trematode taxonomy by providing objective evidence for the delineation of species. However, although the data are objective, the interpretation of these data regarding species boundaries is subjective, especially when different markers conflict. Conserved markers can lead to an underestimation of richness and those used for finer species delineation have the capacity to inflate species recognition, perhaps unrealistically. Here we examine molecular and morphological evidence for species recognition in an especially confusing system, the lepocreadiid genus Preptetos Pritchard, 1960 in acanthuriform fishes of the tropical Indo-west Pacific. We consider species boundaries within this genus based on combined data (ITS2 and 28S rDNA; cox1 mtDNA and morphometrics) for substantial new collections. Delineation of species using only morphological data suggest fewer species than analysis of the sequence data; the latter suggests the presence of potential cryptic species and analysis of different markers suggests the presence of differing numbers of species. We conclude that an integrative interpretation creates the most satisfying taxonomic hypothesis. In the light of the new data, we have chosen and propose a model of trematode species recognition that demands reciprocal monophyly in the most discriminating available molecular marker plus distinction in morphology or host distribution. By invoking these criteria, we distinguish eight species in our new tropical Indo-west Pacific collections. Six of these are new (Preptetos allocaballeroi n. sp., Preptetos paracaballeroi n. sp., Preptetos pearsoni n. sp., Preptetos prudhoei n. sp., Preptetos quandamooka n. sp. and Preptetos zebravaranus n. sp.) and we continue to recognise Preptetos cannoni Barker, Bray & Cribb, 1993 and Preptetos laguncula Bray and Cribb, 1996. Notably; two of the new species, P. allocaballeroi n. sp. and P. paracaballeroi n. sp., are morphologically cryptic relative to each other. Our criteria lead us to recognise, as species, populations with unvarying morphology and similar host relationships but which may have a complex population structure over their range. In our view, this paradigm has the capacity to render tractable the interpretation of the species status of the huge trematode fauna of the tropical Indo-west Pacific.

Habitat-linked genetic structure for white-crowned sparrow (Zonotrichia leucophrys): Local factors shape population genetic structure

Ecological, environmental, and geographic factors all influence genetic structure. Species with broad distributions are ideal systems because they cover a range of ecological and environmental conditions allowing us to test which components predict genetic structure. This study presents a novel, broad geographic approach using molecular markers, morphology, and habitat modeling to investigate rangewide and local barriers causing contemporary genetic differentiation within the geographical range of three white-crowned sparrow (Zonotrichia leucophrys) subspecies: Z. l. gambelii, Z. l. oriantha, and Z. l. pugetensis. Three types of genetic markers showed geographic distance between sampling sites, elevation, and ecosystem type are key factors contributing to population genetic structure. Microsatellite markers revealed white-crowned sparrows do not group by subspecies, but instead indicated four groupings at a rangewide scale and two groupings based on coniferous and deciduous ecosystems at a local scale. Our analyses of morphological variation also revealed habitat differences; sparrows from deciduous ecosystems are larger than individuals from coniferous ecosystems based on principal component analyses. Habitat modeling showed isolation by distance was prevalent in describing genetic structure, but isolation by resistance also had a small but significant influence. Not only do these findings have implications concerning the accuracy of subspecies delineations, they also highlight the critical role of local factors such as habitat in shaping contemporary population genetic structure of species with high dispersal ability.

Effects of natural and seminatural elements on the composition and dispersion of carabid beetles inhabiting an agroecosystem in Northern Italy

The natural and seminatural components of agricultural landscapes play a key role in maintaining a high level of biodiversity. Being the Po Valley one of the most human-dominated and intensively cultivated landscapes in Europe, we investigated the effect of no-crop habitats on carabid richness and composition and evaluated the role of tree row as corridor for forest carabid dispersion. Carabids were sampled with 70 pitfall traps arranged in 35 sampling plots along three parallel transects (80, 100, and 140 m long) and encompassing five different habitats: tree row, tree row edge, grassland, forest edge, and forest. We found 5,615 individuals belonging to 55 species. Despite the similarity in species richness, all the habitats investigated showed a peculiar and distinct species assemblage. The main distinction was between the "open habitat" cluster composed of grassland and tree row edge and the "forest" cluster composed of forest, tree row, and forest edge. We found that forest species are able to penetrate the grassland matrix up to 30 m from the forest edge and that a distance of no more than 60 m between tree row and forest can allow the passage of up to 50% of the forest species. Beyond this distance, the grassland matrix becomes a barrier, preventing them from reaching other suitable habitats. Our findings confirm the importance of maintaining different types of natural habitats to significantly increase biodiversity in an intensively cultivated agroecosystem and demonstrated the role of linear elements as a corridor and "stepping stones" for many forest species.

Temporal and spatial variation in population structure among brooding sea stars in the genus Leptasterias

Temporal genetic studies of low-dispersing organisms are rare. Marine invertebrates lacking a planktonic larval stage are expected to have lower dispersal, low gene flow, and a higher potential for local adaptation than organisms with planktonic dispersal. Leptasterias is a genus of brooding sea stars containing several cryptic species complexes. Population genetic methods were used to resolve patterns of fine-scale population structure in central California Leptasterias species using three loci from nuclear and mitochondrial genomes. Historic samples (collected between 1897 and 1998) were compared to contemporary samples (collected between 2008 and 2014) to delineate changes in species distributions in space and time. Phylogenetic analysis of contemporary samples confirmed the presence of a bay-localized clade and revealed the presence of an additional bay-localized and previously undescribed clade of Leptasterias. Analysis of contemporary and historic samples indicates two clades are experiencing a constriction in their southern range limit and suggests a decrease in clade-specific abundance at sites at which they were once prevalent. Historic sampling revealed a dramatically different distribution of diversity along the California coastline compared to contemporary sampling and illustrates the importance of temporal genetic sampling in phylogeographic studies. These samples were collected prior to significant impacts of Sea Star Wasting Disease (SSWD) and represent an in-depth analysis of genetic structure over 117 years prior to the SSWD-associated mass die-off of Leptasterias.

Evidence that spillover from Marine Protected Areas benefits the spiny lobster (Panulirus interruptus) fishery in southern California

Marine Protected Areas (MPAs) are designed to enhance biodiversity and ecosystem services. Some MPAs are also established to benefit fisheries through increased egg and larval production, or the spillover of mobile juveniles and adults. Whether spillover influences fishery landings depend on the population status and movement patterns of target species both inside and outside of MPAs, as well as the status of the fishery and behavior of the fleet. We tested whether an increase in the lobster population inside two newly established MPAs influenced local catch, fishing effort, and catch-per-unit-effort (CPUE) within the sustainable California spiny lobster fishery. We found greater build-up of lobsters within MPAs relative to unprotected areas, and greater increases in fishing effort and total lobster catch, but not CPUE, in fishing zones containing MPAs vs. those without MPAs. Our results show that a 35% reduction in fishing area resulting from MPA designation was compensated for by a 225% increase in total catch after 6-years, thus indicating at a local scale that the trade-off of fishing ground for no-fishing zones benefitted the fishery.

Extreme seascape drives local recruitment and genetic divergence in brooding and spawning corals in remote north-west Australia

Management strategies designed to conserve coral reefs threatened by climate change need to incorporate knowledge of the spatial distribution of inter- and intra-specific genetic diversity. We characterized patterns of genetic diversity and connectivity using single nucleotide polymorphisms (SNPs) in two reef-building corals to explore the eco-evolutionary processes that sustain populations in north-west Australia. Our sampling focused on the unique reefs of the Kimberley; we collected the broadcast spawning coral Acropora aspera (n = 534) and the brooding coral Isopora brueggemanni (n = 612) across inter-archipelago (tens to hundreds of kilometres), inter-reef (kilometres to tens of kilometres) and within-reef (tens of metres to a few kilometres) scales. Initial analysis of A. aspera identified four highly divergent lineages that were co-occurring but morphologically similar. Subsequent population analyses focused on the most abundant and widespread lineage, Acropora asp-c. Although the overall level of geographic subdivision was greater in the brooder than in the spawner, fundamental similarities in patterns of genetic structure were evident. Most notably, limits to gene flow were observed at scales <35 kilometres. Further, we observed four discrete clusters and a semi-permeable barrier to dispersal that were geographically consistent between species. Finally, sites experiencing bigger tides were more connected to the metapopulation and had greater gene diversity than those experiencing smaller tides. Our data indicate that the inshore reefs of the Kimberley are genetically isolated from neighbouring oceanic bioregions, but occasional dispersal between inshore archipelagos is important for the redistribution of evolutionarily important genetic diversity. Additionally, these results suggest that networks of marine reserves that effectively protect reefs from local pressures should be spaced within a few tens of kilometres to conserve the existing patterns of demographic and genetic connectivity.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Travel with your kin ship! Insights from genetic sibship among settlers of a coral damselfish

Coral reef fish larvae are tiny, exceedingly numerous, and hard to track. They are also highly capable, equipped with swimming and sensory abilities that may influence their dispersal trajectories. Despite the importance of larval input to the dynamics of a population, we remain reliant on indirect insights to the processes influencing larval behavior and transport. Here, we used genetic data (300 independent single nucleotide polymorphisms) derived from a light trap sample of a single recruitment event of Dascyllus abudafur in the Red Sea (N = 168 settlers). We analyzed the genetic composition of the larvae and assessed whether kinship among these was significantly different from random as evidence for cohesive dispersal during the larval phase. We used Monte Carlo simulations of similar-sized recruitment cohorts to compare the expected kinship composition relative to our empirical data. The high number of siblings within the empirical cohort strongly suggests cohesive dispersal among larvae. This work highlights the utility of kinship analysis as a means of inferring dynamics during the pelagic larval phase.

A review of subtidal kelp forests in Ireland: From first descriptions to new habitat monitoring techniques

AIM

Kelp forests worldwide are important marine ecosystems that foster high primary to secondary productivity and multiple ecosystem services. These ecosystems are increasingly under threat from extreme storms, changing ocean temperatures, harvesting, and greater herbivore pressure at regional and global scales, necessitating urgent documentation of their historical to present-day distributions. Species range shifts to higher latitudes have already been documented in some species that dominate subtidal habitats within Europe. Very little is known about kelp forest ecosystems in Ireland, where rocky coastlines are dominated by Laminaria hyperborea. In order to rectify this substantial knowledge gap, we compiled historical records from an array of sources to present historical distribution, kelp and kelp forest recording effort over time, and present rational for the monitoring of kelp habitats to better understand ecosystem resilience.

LOCATION

Ireland (Northern Ireland and Éire).

METHODS

Herbaria, literature from the Linnaean society dating back to late 1700s, journal articles, government reports, and online databases were scoured for information on L. hyperborea. Information about kelp ecosystems was solicited from dive clubs and citizen science groups that are active along Ireland's coastlines.

RESULTS

Data were used to create distribution maps and analyze methodology and technology used to record L. hyperborea presence and kelp ecosystems within Ireland. We discuss the recent surge in studies on Irish kelp ecosystems, fauna associated with kelp ecosystems that may be used as indicators of ecosystem health and suggest methodologies for continued monitoring.

MAIN CONCLUSIONS

While there has been a steady increase in recording effort of the dominant subtidal kelp forest species, L. hyperborea, only recently have studies begun to address other important eco-evolutionary processes at work in kelp forests including connectivity among kelp populations in Ireland. Further monitoring, using suggested methodologies, is required to better understand the resilience of kelp ecosystems in Ireland.

Asymmetrical gene flow in five co-distributed syngnathids explained by ocean currents and rafting propensity

Ocean circulation driving macro-algal rafting is believed to serve as an important mode of dispersal for many marine organisms, leading to predictions on population-level genetic connectivity and the directionality of effective dispersal. Here, we use genome-wide single nucleotide polymorphism data to investigate whether gene flow directionality in two seahorses (Hippocampus) and three pipefishes (Syngnathus) follows the predominant ocean circulation patterns in the Gulf of Mexico and northwestern Atlantic. In addition, we explore whether gene flow magnitudes are predicted by traits related to active dispersal ability and habitat preference. We inferred demographic histories of these co-distributed syngnathid species, and coalescent model-based estimates indicate that gene flow directionality is in agreement with ocean circulation data that predicts eastward and northward macro-algal transport. However, the magnitude to which ocean currents influence this pattern appears strongly dependent on the species-specific traits related to rafting propensity and habitat preferences. Higher levels of gene flow and stronger directionality are observed in Hippocampus erectus, Syngnathus floridae and Syngnathus louisianae, which closely associated with the pelagic macro-algae Sargassum spp., compared to Hippocampus zosterae and the Syngnathus scovelli/Syngnathus fuscus sister-species pair, which prefer near shore habitats and are weakly associated with pelagic Sargassum. This study highlights how the combination of population genomic inference together with ocean circulation data can help explain patterns of population structure and diversity in marine ecosystems.

Time matters: genetic composition and evaluation of effective population size in temperate coastal fish species

Background

Extensive knowledge on the genetic characterization of marine organisms has been assembled, mainly concerning the spatial distribution and structuring of populations. Temporal monitoring assesses not only the stability in genetic composition but also its trajectory over time, providing critical information for the accurate forecast of changes in genetic diversity of marine populations, particularly important for both fisheries and endangered species management. We assessed fluctuations in genetic composition among different sampling periods in the western Portuguese shore in three fish species.

Methods

White seabream Diplodus sargus, sand smelt Atherina presbyter and shanny Lipophrys pholis were chosen, because of their genetic patterns in distinct ecological environments, insight into historical and contemporary factors influencing population effective size (N_e), and degree of commercial exploitation. Samples were obtained near Lisbon between 2003 and 2014 and screened for genetic variation with mitochondrial and nuclear markers. Analyses included genealogies, genetic diversities, temporal structures and contemporary N_e.

Results

For mtDNA no temporal structure was detected, while for nDNA significant differences were recorded between some sampling periods for the shanny and the sand smelt. Haplotype networks revealed deep genealogies, with various levels of diversification. The shanny revealed a smaller N_e/generation when compared to the other species, which, in turn, revealed no evidence of genetic drift for most study periods. These results highlight the fact that temporal variations in genetic pool composition should be considered when evaluating the population structure of fish species with long distance dispersal, which are more vulnerable to recruitment fluctuations.

Seascape genetics of the stalked kelp Pterygophora californica and comparative population genetics in the Santa Barbara Channel

We conducted a population genetic analysis of the stalked kelp, Pterygophora californica, in the Santa Barbara Channel, California, USA. The results were compared with previous work on the genetic differentiation of giant kelp, Macrocystis pyrifera, in the same region. These two sympatric kelps not only share many life history and dispersal characteristics but also differ in that dislodged P. californica does not produce floating rafts with buoyant fertile sporophytes, commonly observed for M. pyrifera. We used a comparative population genetic approach with these two species to test the hypothesis that the ability to produce floating rafts increases the genetic connectivity among kelp patches in the Santa Barbara Channel. We quantified the association of habitat continuity and oceanographic distance with the genetic differentiation observed in stalked kelp, like previously conducted for giant kelp. We compared both overall (across all patches) and pairwise (between patches) genetic differentiation. We found that oceanographic transit time, habitat continuity, and geographic distance were all associated with genetic connectivity in P. californica, supporting similar previous findings for M. pyrifera. Controlling for differences in heterozygosity between kelp species using Jost's DEST , we showed that global differentiation and pairwise differentiation were similar among patches between the two kelp species, indicating that they have similar dispersal capabilities despite their differences in rafting ability. These results suggest that rafting sporophytes do not play a significant role in effective dispersal of M. pyrifera at ecologically relevant spatial and temporal scales.

High Connectivity of the White Seabream (Diplodus sargus, L. 1758) in the Aegean Sea, Eastern Mediterranean Basin

Population dynamics in the marine realm can shape species' spatial structure and genetic variability between given geographical areas. Connectivity is an important factor of species' population structure. In this study, we examined the genetic diversity and structure of white seabream (Diplodus sargus, L. 1758) in the eastern Mediterranean basin, using a panel of four microsatellite markers. Recorded low F_ST values within the study area indicate little evidence of genetic differentiation among populations. Results suggest high gene flow which may imply near-panmixia between populations, indicating the possibility of a probable movement of adult migrants, or strong passive drift at sea in early life stages of the species. To this extent, bibliographically speaking, different species within the Sparidae family favor altered population dynamics patterns with respect to local populations and genetic divergence, in the context of the molecular marker used.

Asymmetric oceanographic processes mediate connectivity and population genetic structure, as revealed by RADseq, in a highly dispersive marine invertebrate (Parastichopus californicus)

Marine populations are typically characterized by weak genetic differentiation due to the potential for long-distance dispersal favouring high levels of gene flow. However, strong directional advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (F_ST  = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R²  = .49) than geographic distance (IBD) (R²  = .18), and directional processes played an important role in shaping fine-scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.

Isolation by environment in the highly mobile olive ridley turtle (Lepidochelys olivacea) in the eastern Pacific

Spatial and temporal scales at which processes modulate genetic diversity over the landscape are usually overlooked, impacting the design of conservation management practices for widely distributed species. We examine processes shaping population divergence in highly mobile species by re-assessing the case of panmixia in the iconic olive ridley turtle from the eastern Pacific. We implemented a biophysical model of connectivity and a seascape genetic analysis based on nuclear DNA variation of 634 samples collected from 27 nesting areas. Two genetically distinct populations largely isolated during reproductive migrations and mating were detected, each composed of multiple nesting sites linked by high connectivity. This pattern was strongly associated with a steep environmental gradient and also influenced by ocean currents. These findings relate to meso-scale features of a dynamic oceanographic interface in the eastern tropical Pacific (ETP) region, a scenario that possibly provides different cost-benefit solutions and selective pressures for sea turtles during both the mating and migration periods. We reject panmixia and propose a new paradigm for olive ridley turtles where reproductive isolation due to assortative mating is linked to its environment. Our study demonstrates the relevance of integrative approaches for assessing the role of environmental gradients and oceanographic currents as drivers of genetic differentiation in widely distributed marine species. This is relevant for the conservation management of species of highly mobile behaviour, and assists the planning and development of large-scale conservation strategies for the threatened olive ridley turtles in the ETP.

Combining Bayesian genetic clustering and ecological niche modeling: Insights into wolf intraspecific genetic structure

The distribution of intraspecific genetic variation and how it relates to environmental factors is of increasing interest to researchers in macroecology and biogeography. Recent studies investigated the relationships between the environment and patterns of intraspecific genetic variation across species ranges but only few rigorously tested the relation between genetic groups and their ecological niches. We quantified the relationship of genetic differentiation (F ST) and the overlap of ecological niches (as measured by n-dimensional hypervolumes) among genetic groups resulting from spatial Bayesian genetic clustering in the wolf (Canis lupus) in the Italian peninsula. Within the Italian wolf population, four genetic clusters were detected, and these clusters showed different ecological niches. Moreover, different wolf clusters were significantly related to differences in land cover and human disturbance features. Such differences in the ecological niches of genetic clusters should be interpreted in light of neutral processes that hinder movement, dispersal, and gene flow among the genetic clusters, in order to not prematurely assume any selective or adaptive processes. In the present study, we found that both the plasticity of wolves-a habitat generalist-to cope with different environmental conditions and the occurrence of barriers that limit gene flow lead to the formation of genetic intraspecific genetic clusters and their distinct ecological niches.

Geographic isolation and larval dispersal shape seascape genetic patterns differently according to spatial scale

Genetic variation, as a basis of evolutionary change, allows species to adapt and persist in different climates and environments. Yet, a comprehensive assessment of the drivers of genetic variation at different spatial scales is still missing in marine ecosystems. Here, we investigated the influence of environment, geographic isolation, and larval dispersal on the variation in allele frequencies, using an extensive spatial sampling (47 locations) of the striped red mullet (Mullus surmuletus) in the Mediterranean Sea. Univariate multiple regressions were used to test the influence of environment (salinity and temperature), geographic isolation, and larval dispersal on single nucleotide polymorphism (SNP) allele frequencies. We used Moran's eigenvector maps (db-MEMs) and asymmetric eigenvector maps (AEMs) to decompose geographic and dispersal distances in predictors representing different spatial scales. We found that salinity and temperature had only a weak effect on the variation in allele frequencies. Our results revealed the predominance of geographic isolation to explain variation in allele frequencies at large spatial scale (>1,000 km), while larval dispersal was the major predictor at smaller spatial scale (<1,000 km). Our findings stress the importance of including spatial scales to understand the drivers of spatial genetic variation. We suggest that larval dispersal allows to maintain gene flows at small to intermediate scale, while at broad scale, genetic variation may be mostly shaped by adult mobility, demographic history, or multigenerational stepping-stone dispersal. These findings bring out important spatial scale considerations to account for in the design of a protected area network that would efficiently enhance protection and persistence capacity of marine species.

Regional variation in drivers of connectivity for two frog species (Rana pretiosa and R. luteiventris) from the U.S. Pacific Northwest

Comparative landscape genetics has uncovered high levels of variability in which landscape factors affect connectivity among species and regions. However, the relative importance of species traits versus environmental variation for predicting landscape patterns of connectivity is unresolved. We provide evidence from a landscape genetics study of two sister taxa of frogs, the Oregon spotted frog (Rana pretiosa) and the Columbia spotted frog (Rana luteiventris) in Oregon and Idaho, USA. Rana pretiosa is relatively more dependent on moisture for dispersal than R. luteiventris, so if species traits influence connectivity, we predicted that connectivity among R. pretiosa populations would be more positively associated with moisture than R. luteiventris. However, if environmental differences are important drivers of gene flow, we predicted that connectivity would be more positively related to moisture in arid regions. We tested these predictions using eight microsatellite loci and gravity models in two R. pretiosa regions and four R. luteiventris regions (n = 1,168 frogs). In R. pretiosa, but not R. luteiventris, connectivity was positively related to mean annual precipitation, supporting our first prediction. In contrast, connectivity was not more positively related to moisture in more arid regions. Various temperature metrics were important predictors for both species and in all regions, but the directionality of their effects varied. Therefore, the pattern of variation in drivers of connectivity was consistent with predictions based on species traits rather than on environmental variation.

Population genetic signatures of a climate change driven marine range extension

Shifts in species distribution, or 'range shifts', are one of the most commonly documented responses to ocean warming, with important consequences for the function and structure of ecosystems, and for socio-economic activities. Understanding the genetic signatures of range shifts can help build our knowledge of the capacity of species to establish and persist in colonised areas. Here, seven microsatellite loci were used to examine the population connectivity, genetic structure and diversity of Octopus tetricus, which has extended its distribution several hundred kilometres polewards associated with the southwards extension of the warm East Australian Current along south-eastern Australia. The historical distribution and the range extension zones had significant genetic differences but levels of genetic diversity were comparable. The population in the range extension zone was sub-structured, contained relatively high levels of self-recruitment and was sourced by migrants from along the entire geographic distribution. Genetic bottlenecks and changes in population size were detected throughout the range extension axis. Persistent gene flow from throughout the historical zone and moderate genetic diversity may buffer the genetic bottlenecks and favour the range extension of O. tetricus. These characteristics may aid adaptation, establishment, and long-term persistence of the population in the range extension zone.

Seascape genetics and biophysical connectivity modelling support conservation of the seagrass Zostera marina in the Skagerrak-Kattegat region of the eastern North Sea

Maintaining and enabling evolutionary processes within meta-populations are critical to resistance, resilience and adaptive potential. Knowledge about which populations act as sources or sinks, and the direction of gene flow, can help to focus conservation efforts more effectively and forecast how populations might respond to future anthropogenic and environmental pressures. As a foundation species and habitat provider, Zostera marina (eelgrass) is of critical importance to ecosystem functions including fisheries. Here, we estimate connectivity of Z. marina in the Skagerrak-Kattegat region of the North Sea based on genetic and biophysical modelling. Genetic diversity, population structure and migration were analysed at 23 locations using 20 microsatellite loci and a suite of analytical approaches. Oceanographic connectivity was analysed using Lagrangian dispersal simulations based on contemporary and historical distribution data dating back to the late 19th century. Population clusters, barriers and networks of connectivity were found to be very similar based on either genetic or oceanographic analyses. A single-generation model of dispersal was not realistic, whereas multigeneration models that integrate stepping-stone dispersal and extant and historic distribution data were able to capture and model genetic connectivity patterns well. Passive rafting of flowering shoots along oceanographic currents is the main driver of gene flow at this spatial-temporal scale, and extant genetic connectivity strongly reflects the "ghost of dispersal past" sensu Benzie, 1999. The identification of distinct clusters, connectivity hotspots and areas where connectivity has become limited over the last century is critical information for spatial management, conservation and restoration of eelgrass.

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment

BACKGROUND

Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist.

RESULTS

Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region.

CONCLUSIONS

Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Independent estimates of marine population connectivity are more concordant when accounting for uncertainties in larval origins

Marine larval dispersal is a complex biophysical process that depends on the effects of species biology and oceanography, leading to logistical difficulties in estimating connectivity among populations of marine animals with biphasic life cycles. To address this challenge, the application of multiple methodological approaches has been advocated, in order to increase confidence in estimates of population connectivity. However, studies seldom account for sources of uncertainty associated with each method, which undermines a direct comparative approach. In the present study we explicitly account for the statistical uncertainty in observed connectivity matrices derived from elemental chemistry of larval mussel shells, and compare these to predictions from a biophysical model of dispersal. To do this we manipulate the observed connectivity matrix by applying different confidence levels to the assignment of recruits to source populations, while concurrently modelling the intrinsic misclassification rate of larvae to known sources. We demonstrate that the correlation between the observed and modelled matrices increases as the number of observed recruits classified as unknowns approximates the observed larval misclassification rate. Using this approach, we show that unprecedented levels of concordance in connectivity estimates (r = 0.96) can be achieved, and at spatial scales (20-40 km) that are ecologically relevant.

Population genetic structure in Sabatieria (Nematoda) reveals intermediary gene flow and admixture between distant cold seeps from the Mediterranean Sea

BACKGROUND

There is a general lack of information on the dispersal and genetic structuring for populations of small-sized deep-water taxa, including free-living nematodes which inhabit and dominate the seafloor sediments. This is also true for unique and scattered deep-sea habitats such as cold seeps. Given the limited dispersal capacity of marine nematodes, genetic differentiation between such geographically isolated habitat patches is expected to be high. Against this background, we examined genetic variation in both mitochondrial (COI) and nuclear (18S and 28S ribosomal) DNA markers of 333 individuals of the genus Sabatieria, abundantly present in reduced cold-seep sediments. Samples originated from four Eastern Mediterranean cold seeps, separated by hundreds of kilometers, and one seep in the Southeast Atlantic.

RESULTS

Individuals from the Mediterranean and Atlantic were divided into two separate but closely-related species clades. Within the Eastern Mediterranean, all specimens belonged to a single species, but with a strong population genetic structure (ΦST = 0.149). The haplotype network of COI contained 19 haplotypes with the most abundant haplotype (52% of the specimens) shared between all four seeps. The number of private haplotypes was high (15), but the number of mutations between haplotypes was low (1-8). These results indicate intermediary gene flow among the Mediterranean Sabatieria populations with no evidence of long-term barriers to gene flow.

CONCLUSIONS

The presence of shared haplotypes and multiple admixture events indicate that Sabatieria populations from disjunct cold seeps are not completely isolated, with gene flow most likely facilitated through water current transportation of individuals and/or eggs. Genetic structure and molecular diversity indices are comparable to those of epiphytic shallow-water marine nematodes, while no evidence of sympatric cryptic species was found for the cold-seep Sabatieria.

Distinct genetic differentiation and species diversification within two marine nematodes with different habitat preference in Antarctic sediments

Background

Dispersal ability, population genetic structure and species divergence in marine nematodes are still poorly understood, especially in remote areas such as the Southern Ocean. We investigated genetic differentiation of species and populations of the free-living endobenthic nematode genera Sabatieria and Desmodora using nuclear 18S rDNA, internal transcribed spacer (ITS) rDNA, and mitochondrial cytochrome oxidase I (COI) gene sequences. Specimens were collected at continental shelf depths (200–500 m) near the Antarctic Peninsula, Scotia Arc and eastern side of the Weddell Sea. The two nematode genera co-occurred at all sampled locations, but with different vertical distribution in the sediment. A combination of phylogenetic (GMYC, Bayesian Inference, Maximum Likelihood) and population genetic (AMOVA) analyses were used for species delimitation and assessment of gene flow between sampling locations.

Results

Sequence analyses resulted in the delimitation of four divergent species lineages in Sabatieria, two of which could not be discriminated morphologically and most likely constitute cryptic species. Two species were recognised in Desmodora, one of which showed large intraspecific morphological variation. Both genera comprised species that were restricted to one side of the Weddell Sea and species that were widely spread across it. Population genetic structuring was highly significant and more pronounced in the deeper sediment-dwelling Sabatieria species, which are generally less prone to resuspension and passive dispersal in the water column than surface Desmodora species.

Conclusions

Our results indicate that gene flow is restricted at large geographic distance in the Southern Ocean, which casts doubt on the efficiency of the Weddell gyre and Antarctic Circumpolar Current in facilitating circum-Antarctic nematode species distributions. We also show that genetic structuring and cryptic speciation can be very different in nematode species isolated from the same geographic area, but with different habitat preferences (surface versus deeper sediment layers).

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0968-1) contains supplementary material, which is available to authorized users.

Considering reefscape configuration and composition in biophysical models advance seascape genetics

Previous seascape genetics studies have emphasized the role of ocean currents and geographic distances to explain the genetic structure of marine species, but the role of benthic habitat has been more rarely considered. Here, we compared the population genetic structure observed in West Pacific giant clam populations against model simulations that accounted habitat composition and configuration, geographical distance, and oceanic currents. Dispersal determined by geographical distance provided a modelled genetic structure in better agreement with the observations than dispersal by oceanic currents, possibly due to insufficient spatial resolution of available oceanographic and coastal circulation models. Considering both habitat composition and configuration significantly improved the match between simulated and observed genetic structures. This study emphasizes the importance of a reefscape genetics approach to population ecology, evolution and conservation in the sea.

Life-history predicts past and present population connectivity in two sympatric sea stars

Life‐history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic‐developing sea star, Parvulastra exigua, has lower levels of within‐ and among‐population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic‐developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic‐developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post‐Pleistocene introgression.