Oceanography promotes self-recruitment in a planktonic larval disperser

How long have we been mistaken? Multi-tools shedding light into the systematics of the widespread deep-water genus Madrepora Linnaeus, 1758 (Scleractinia)

Deep-water coral reefs are found worldwide and harbor biodiversity levels that are comparable to their shallow-water counterparts. However, the genetic diversity and population structure of deep-water species remain poorly explored, and historical taxonomical issues still need to be resolved. Here we used microsatellite markers as well as ultraconserved elements (UCE) and exons to shed light on the population structure, genetic diversity, and phylogenetic position of the genus Madrepora, which contains M. oculata, one of the most widespread scleractinian species. Population structure of 107 samples from three Southwestern Atlantic sedimentary basins revealed the occurrence of a cryptic species, herein named M. piresae sp. nov. (authored by Kitahara, Capel and Zilberberg), which can be found in sympatry with M. oculata. Phylogeny reconstructions based on 134 UCEs and exon regions corroborated the population genetic data, with the recovery of two well-supported groups, and reinforced the polyphyly of the family Oculinidae. In order to better accommodate the genus Madrepora, while reducing taxonomical confusion associated with the name Madreporidae, we propose the monogeneric family Bathyporidae fam. nov. (authored by Kitahara, Capel, Zilberberg and Cairns). Our findings advance the knowledge on the widespread deep-water genus Madrepora, resolve a long-standing question regarding the phylogenetic position of the genus, and highlight the need of a worldwide review of the genus.

Investigating population dynamics from parentage analysis in the highly endangered fan mussel Pinna nobilis

Understanding dispersal patterns is a major focus for conservation biology as it influences local survival and resilience in case of local disturbance, particularly for sessile species. Dispersal can be assessed through parentage analyses by estimating family structure and self-recruitment. This study documents the family structure of a pelagic spawner, Pinna nobilis, which is facing a major crisis that threatens its survival as most of its populations have been decimated by a parasite, Haplosporidium pinnae. In this context, we focused on a single population (Peyrefite, Banyuls-sur-mer, France) where 640 individuals were sampled in 2011, 2015, and 2018 and genotyped for 22 microsatellite markers. Genetic diversity was high and homogeneous among years, with mean allele numbers ranging between 13.6 and 14.8 and observed heterozygosities (H o) between 0.7121 and 0.7331. Low, but significant, genetic differentiations were found between 2011-2015 and 2015-2018. A parentage analysis described 11 clusters, including one prevailing, and revealed that 46.9% of individuals were involved in half-sib relationships, even between years, suggesting that source populations were recurrent year after year. There were few individuals resampled between years (30 in 2015 and 14 in 2018), indicating a rapid turnover. Considering the large number of half-sib relationships but the low number of relations per individual, we conclude that P. nobilis exhibit homogeneous reproductive success. Self-recruitment was not detected, making this population highly vulnerable as replenishment only relies on connectivity from neighboring populations. In the context of the pandemic caused by H. pinnae, these results will have to be considered when choosing a location to reintroduce individuals in potential future rescue plans.

Comparative population genetics of swimming crab host (Portunus pelagicus) and common symbiotic barnacle (Octolasmis angulata) in Vietnam

Background

By comparing spatial geographical structures of host populations with that of their symbionts light can be shed on their biological interactions, and the degree of congruence between host and symbiont phylogeographies should reflect their life histories and especially dispersal mechanisms.

Methods

Here, we analyzed the genetic diversity and structure of a host, the blue swimming crab, Portunus pelagicus, and its symbiotic pedunculate barnacle Octolasmis angulata from six location sites representing three geographic regions (north, central and south) along the Vietnam coastline. High levels of congruence in their phylogeographic patterns were expected as they both undergo planktonic larval stages.

Results

Based on the COI mtDNA markers, O. angulata populations showed higher genetic diversity in comparison with their host P. pelagicus (number of haplotype/individuals, haplotype and nucleotide diversity are 119/192, 0.991 ± 0.002 and 0.02; and 89/160, 0.913 ± 0.02 and 0.015, respectively). Pairwise Fst and AMOVA analyses showed a more pronounced population structure in the symbiotic barnacle than in its crab host. The DAPC analyses identified three genetic clusters. However, both haplotype networks and scatter plots supported connectivity of the host and the symbiotic barnacle throughout their distribution range, except for low subdivision of southern population. Isolation by distance were detected only for the symbiont O. angulata (R² = 0.332, P = 0.05), while dbMEM supported spatial structure of both partners, but only at MEM-1 (Obs. 0.2686, P < 0.01 and Obs. 0.2096, P < 0.01, respectively).

Genomics-informed models reveal extensive stretches of coastline under threat by an ecologically dominant invasive species

Explaining why some species are widespread, while others are not, is fundamental to biogeography, ecology, and evolutionary biology. A unique way to study evolutionary and ecological mechanisms that either limit species' spread or facilitate range expansions is to conduct research on species that have restricted distributions. Nonindigenous species, particularly those that are highly invasive but have not yet spread beyond the introduced site, represent ideal systems to study range size changes. Here, we used species distribution modeling and genomic data to study the restricted range of a highly invasive Australian marine species, the ascidian Pyura praeputialis This species is an aggressive space occupier in its introduced range (Chile), where it has fundamentally altered the coastal community. We found high genomic diversity in Chile, indicating high adaptive potential. In addition, genomic data clearly showed that a single region from Australia was the only donor of genotypes to the introduced range. We identified over 3,500 km of suitable habitat adjacent to its current introduced range that has so far not been occupied, and importantly species distribution models were only accurate when genomic data were considered. Our results suggest that a slight change in currents, or a change in shipping routes, may lead to an expansion of the species' introduced range that will encompass a vast portion of the South American coast. Our study shows how the use of population genomics and species distribution modeling in combination can unravel mechanisms shaping range sizes and forecast future range shifts of invasive species.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

Larval behaviour, dispersal and population connectivity in the deep sea

Ecosystem connectivity is an essential consideration for marine spatial planning of competing interests in the deep sea. Immobile, adult communities are connected through freely floating larvae, depending on new recruits for their health and to adapt to external pressures. We hypothesize that the vertical swimming ability of deep-sea larvae, before they permanently settle at the bottom, is one way larvae can control dispersal. We test this hypothesis with more than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3\times 10^{8}$$\end{document}3×108 simulated particles with a range of active swimming behaviours embedded within the currents of a high-resolution ocean model. Despite much stronger horizontal ocean currents, vertical swimming of simulated larvae can have an order of magnitude impact on dispersal. These strong relationships between larval dispersal, pathways, and active swimming demonstrate that lack of data on larval behaviour traits is a serious impediment to modelling deep-sea ecosystem connectivity; this uncertainty greatly limits our ability to develop ecologically coherent marine protected area networks.

Genome-wide SNP analyses reveal population structure of Portunus pelagicus along Vietnam coastline

The blue swimming crab (Portunus pelagicus Linnaeus, 1758) is one of the commercially exploited crab fishery resources in Vietnam. This is the first study to provide a broad survey of genetic diversity, population structure and migration patterns of P. pelagicus along the Vietnamese coastline. The crab samples were collected from northern, central and southern Vietnam. Here, we used a panel of single nucleotide polymorphisms (SNPs) generated from restriction site-associated DNA sequencing (RADseq). After removing 32 outlier loci, 306 putatively neutral SNPs from 96 individuals were used to assess fine-scale population structure of blue swimming crab. The mean observed heterozygosity (Ho) and expected heterozygosity (He) per locus was 0.196 and 0.223, respectively. Pairwise Fst and hierarchical AMOVA supported significant differentiation of central and northern from southern populations (P<0.01). Population structure analyses revealed that P. pelagicus in the south is a separate fisheries unit from the north and center. Contemporary migration patterns supported high migration between northern and central populations and restricted genetic exchange within the southern population. In contrast, historic gene flow provides strong evidence for single panmictic population. The results are useful for understanding current status of P. pelagicus in the wild under an environment changing due to natural and anthropogenic stresses, with implications for fisheries management.

Seascape genetics of the spiny lobster Panulirus homarus in the Western Indian Ocean: Understanding how oceanographic features shape the genetic structure of species with high larval dispersal potential

This study examines the fine-scale population genetic structure and phylogeography of the spiny lobster Panulirus homarus in the Western Indian Ocean. A seascape genetics approach was used to relate the observed genetic structure based on 21 microsatellite loci to ocean circulation patterns, and to determine the influence of latitude, sea surface temperature (SST), and ocean turbidity (KD490) on population-level processes. At a geospatial level, the genetic clusters recovered corresponded to three putative subspecies, P. h. rubellus from the SW Indian Ocean, P. h. megasculptus from the NW Indian Ocean, and P. h. homarus from the tropical region in-between. Virtual passive Lagrangian particles advected using satellite-derived ocean surface currents were used to simulate larval dispersal. In the SW Indian Ocean, the dispersion of particles tracked over a 4-month period provided insight into a steep genetic gradient observed at the Delagoa Bight, which separates P. h. rubellus and P. h. homarus. South of the contact zone, particles were advected southwestwards by prevailing boundary currents or were retained in nearshore eddies close to release locations. Some particles released in southeast Madagascar dispersed across the Mozambique Channel and reached the African shelf. Dispersal was characterized by high seasonal and inter-annual variability, and a large proportion of particles were dispersed far offshore and presumably lost. In the NW Indian Ocean, particles were retained within the Arabian Sea. Larval retention and self-recruitment in the Arabian Sea could explain the recent genetic divergence between P. h. megasculptus and P. h. homarus. Geographic distance and minimum SST were significantly associated with genetic differentiation in multivariate analysis, suggesting that larval tolerance to SST plays a role in shaping the population structure of P. homarus.

High level of genetic connectivity in a deep-water reef fish, Caulolatilus microps

The life-history characteristics of blueline tilefish make them particularly vulnerable to overfishing. Although North Carolina, U.S.A. was previously considered to be the northern extent of the range of C. microps, concentrations have recently been discovered in Virginia and Maryland, with reports as far north as Montauk, New York. Attempts to manage the fishery have been hampered by a lack of information about whether the U.S. East Coast includes multiple stocks. To assess the appropriateness of alternate management options, we used 25 variable microsatellite loci and sequencing of the mitochondrial (mt)DNA control region (CR) to evaluate the genetic structure of 490 C. microps sampled from across the U.S. East Coast range. Pairwise comparisons of genetic differentiation among collection locations based on both nuclear microsatellite and mtCR sequence data were all low and non significant. No significant autocorrelation was observed across multiple distance classes, consistent with widespread dispersal., Although the assumed sedentary nature of adult C. microps suggests population structuring, the genetic data were consistent with the presence of sufficient gene flow to prevent the accumulation of significant genetic differences and supports management of C. microps as a single stock along the U.S. East Coast.

Consistency in the supply of larval fishes among coral reefs in French Polynesia

For marine fishes with a bipartite life cycle, pelagic larval dispersal can shape the distribution, connectivity, composition and resilience of adult populations. Numerous studies of larval dispersal, and associated settlement and recruitment processes, have examined the relationship between population connectivity and oceanographic features. However, relatively little is known about spatial and temporal variation in the abundance of larvae settling among different reefs and the extent to which the species assemblage of larvae settling at one location is reflective of the assemblage in neighbouring areas. Here, using crest nets, which provide a non-selective measure of the total abundance and assemblage of larvae settling to a reef (i.e. larval supply), we collected larval coral reef fishes at five locations surrounding two spatially disparate French Polynesian islands: Moorea and Nengo-Nengo. Overall, larval settlement patterns were correlated with the lunar cycle, with larval abundance peaking during the new moon. Although there were some spatial differences in larval supply among the five monitored sites, settlement patterns were largely consistent, even at the species level, irrespective of factors such as coastline orientation or distance between sites. This study provides further insights into the mechanisms driving patterns of dispersal and settlement of larval fishes over large spatial scales.

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.