Very fine-scale population genetic structure of sympatric asterinid sea stars with benthic and pelagic larvae: influence of mating system and dispersal potential

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.

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.

Reproduction in Urbanised Coastal Waters: Shallow-Water Sea Anemones (Entacmaea quadricolor and Stichodactyla haddoni) Maintain High Genetic Diversity and Panmixia

Sea anemones are sedentary marine animals that tend to disperse via planktonic larvae and are predicted to have high population connectivity in undisturbed habitats. We test whether two sea anemone species living in two different tidal zones of a highly disturbed marine environment can maintain high genetic connectivity. More than 1000 loci with single-nucleotide polymorphisms (SNPs) were obtained with double-digest RADseq for 81 Stichodactyla haddoni and 99 Entacmaea quadricolor individuals to test for population genetic structure. We find evidence that both species predominantly propagate via sexual reproduction, and asexual reproduction is limited. We observe panmixia that indicates the absence of effective dispersal barriers for these species living in a highly anthropogenically disturbed environment. This is positive news for both species that are also found in the aquarium trade. More fundamentally, our results suggest that inhabiting different parts of a shallow reef may not affect a species’ population connectivity nor favour asexual reproduction.

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.

Fine-scale hierarchical genetic structure and kinship analysis of the ascidian Pyura chilensis in the southeastern Pacific

Studying population structure and genetic diversity at fine spatial scales is key for a better understanding of demographic processes that influence population connectivity. This is particularly important in marine benthic organisms that rely on larval dispersal to maintain connectivity among populations. Here, we report the results of a genetic survey of the ascidian Pyura chilensis from three localities along the southeastern Pacific. This study follows up on a previous report that described a genetic break in this region among localities only 20 km apart. By implementing a hierarchical sampling design at four spatial levels and using ten polymorphic microsatellite markers, we test whether differences in fine-scale population structure explain the previously reported genetic break. We compared genetic spatial autocorrelations, as well as kinship and relatedness distributions within and among localities adjacent to the genetic break. We found no evidence of significant autocorrelation at the scale up to 50 m despite the low dispersal potential of P. chilensis that has been reported in the literature. We also found that the proportion of related individuals in close proximity (<1 km) was higher than the proportion of related individuals further apart. These results were consistent in the three localities. Our results suggest that the spatial distribution of related individuals can be nonrandom at small spatial scales and suggests that dispersal might be occasionally limited in this species or that larval cohorts can disperse in the plankton as clustered groups. Overall, this study sheds light on new aspects of the life of this ascidian as well as confirms the presence of a genetic break at 39°S latitude. Also, our data indicate there is not enough evidence to confirm that this genetic break can be explained by differences in fine-scale genetic patterns among localities.

Can transplanting enhance mobile marine invertebrates in ecologically engineered rock pools?

The field of eco-engineering has burgeoned in recent years in response to the proliferation of artificial structures. Adding water-retaining features to seawalls has been successful in increasing biodiversity relative to the surrounding structure. Artificial rock pools may not, however, completely mimic natural rock pools. Here, we compared natural colonisation, through dispersal and recruitment, of intertidal mobile species to water-retaining flowerpots on seawalls with that into rock pools. This represents the more usual 'passive' approach to eco-engineering where features are built to enhance biodiversity and are allowed to colonise naturally, as opposed to seeding or transplanting organisms to features. While flowerpots supported some mobile species not found on the seawall, other species common on natural shores did not recruit to flowerpots. Thus, in a second experiment we tested the effectiveness of an 'active' approach through transplanting mobile organisms to flowerpots to expedite the colonisation process. For the species examined, however, most individuals did not stay in the flowerpots for more than 24 h after being transplanted. Further understanding of the processes (e.g. dispersal distances, recruitment) influencing colonisation of eco-engineered habitats is needed to effectively inform management of marine infrastructure, particularly for projects targeted at restoration rather than enhancement.

Genome-wide comparisons reveal a clinal species pattern within a holobenthic octopod-the Australian Southern blue-ringed octopus, Hapalochlaena maculosa (Cephalopoda: Octopodidae)

The southern blue-ringed octopus, Hapalochlaena maculosa (Hoyle, 1883) lacks a planktonic dispersal phase, yet ranges across Australia's southern coastline. This species' brief and holobenthic life history suggests gene flow might be limited, leaving distant populations prone to strong genetic divergence. This study used 17,523 genome-wide SNP loci to investigate genetic structuring and local adaptation patterns of H. maculosa among eight sampling sites along its reported range. Within sites, interrelatedness was very high, consistent with the limited dispersal of this taxon. However, inbreeding coefficients were proportionally lower among sites where substructuring was not detected, suggesting H. maculosa might possess a mechanism for inbreeding avoidance. Genetic divergence was extremely high among all sites, with the greatest divergence observed between both ends of the distribution, Fremantle, WA, and Stanley, TAS. Genetic distances closely followed an isolation by geographic distance pattern. Outlier analyses revealed distinct selection signatures at all sites, with the strongest divergence reported between Fremantle and the other Western Australian sites. Phylogenetic reconstructions using the described sister taxon H. fasciata (Hoyle, 1886) further supported that the genetic divergence between distal H. maculosa sites in this study was equivalent to that of between established heterospecifics within this genus. However, it is advocated that taxonomic delineations within this species should be made with caution. These data indicate that H. maculosa forms a clinal species pattern across its geographic range, with gene flow present through allele sharing between adjacent populations. Morphological investigations are recommended for a robust resolution of the taxonomic identity and ecotype boundaries of this species.

Do larval types affect genetic connectivity at sea? Testing hypothesis in two sibling marine gastropods with contrasting larval development

In marine environments, connectivity among populations of benthic invertebrates is provided primarily by dispersion of larvae, with the duration of pelagic larval phase (PLD) supposed to represent one of the major factor affecting connectivity. In marine gastropods, PLD is linked to specific larval development types, which may be entirely intracapsular (thus lacking a pelagic dispersal), or include a short pelagic lecithotrophic or a long planktotrophic phase. In the present study, we investigated two sibling species of the cosmopolitan neogastropod genus Columbella (commonly known as dove shells): Columbella adansoni Menke, 1853, from the Macaronesian Atlantic archipelagos, with planktotrophic development, and Columbella rustica Linnaeus, 1758, from the Mediterranean Sea, with intracapsular development. We expected to find differences between these two sister species, in terms of phylogeographic structure, levels of genetic diversification and spatial distribution of genetic diversity, if PLD was actually a relevant factor affecting connectivity. By analysing the sequence variation at the cytochrome c oxidase subunit I (COI) in 167 specimens of the two species, collected over a comparable geographic range, we found that Columbella adansoni, the species with planktotrophic development, and thus longer PLD, showed no phylogeographic structure, lower levels of genetic diversity, interpopulational variance lower than the intrapopulational one and no spatial structure in the distribution of the genetic diversity; Columbella rustica, the species with intracapsular development, thus with evidently lower dispersal abilities, showed a clear phylogeographic structure, higher levels of genetic diversity, high interpopulational and low intrapopulational variance, and a clear signature of global spatial structure in the distribution of the genetic diversity. Thus, in this study, two sibling species differing almost only in their larval ecology (and PLD), when compared for their genetic variation showed patterns supporting the hypothesis that PLD is a major factor affecting genetic connectivity. Therefore, it seems reasonable to expect that the ecological attributes of the marine communities - also in terms of the variation in larval ecology of the species involved - are taken into the due consideration in conservation actions, like the design of marine protected areas networks.

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.

Low genetic diversity and recent demographic expansion in the red starfish Echinaster sepositus (Retzius 1816)

Understanding the phylogeography and genetic structure of populations and the processes responsible of patterns therein is crucial for evaluating the vulnerability of marine species and developing management strategies. In this study, we explore how past climatic events and ongoing oceanographic and demographic processes have shaped the genetic structure and diversity of the Atlanto-Mediterranean red starfish Echinaster sepositus. The species is relatively abundant in some areas of the Mediterranean Sea, but some populations have dramatically decreased over recent years due to direct extraction for ornamental aquariums and souvenir industries. Analyses across most of the distribution range of the species based on the mitochondrial cytochrome c oxidase subunit I gene and eight microsatellite loci revealed very low intraspecific genetic diversity. The species showed a weak genetic structure within marine basins despite the a priori low dispersal potential of its lecithotrophic larva. Our results also revealed a very recent demographic expansion across the distribution range of the species. The genetic data presented here indicate that the species might be highly vulnerable, due to its low intraspecific genetic diversity.

Reduced genetic diversity and increased reproductive isolation follow population-level loss of larval dispersal in a marine gastropod

Population-level consequences of dispersal ability remain poorly understood, especially for marine animals in which dispersal is typically considered a species-level trait governed by oceanographic transport of microscopic larvae. Transitions from dispersive (planktotrophic) to nondispersive, aplanktonic larvae are predicted to reduce connectivity, genetic diversity within populations, and the spatial scale at which reproductive isolation evolves. However, larval dimorphism within a species is rare, precluding population-level tests. We show the sea slug Costasiella ocellifera expresses both larval morphs in Florida and the Caribbean, regions with divergent mitochondrial lineages. Planktotrophy predominated at 11 sites, 10 of which formed a highly connected and genetically diverse Caribbean metapopulation. Four populations expressed mainly aplanktonic development and had markedly reduced connectivity, and lower genetic diversity at one mitochondrial and six nuclear loci. Aplanktonic dams showed partial postzygotic isolation in most interpopulation crosses, regardless of genetic or geographic distance to the sire's source, suggesting that outbreeding depression affects fragmented populations. Dams from genetically isolated and neighboring populations also exhibited premating isolation, consistent with reinforcement contingent on historical interaction. By increasing self-recruitment and genetic drift, the loss of dispersal may thus initiate a feedback loop resulting in the evolution of reproductive isolation over small spatial scales in the sea.

Influence of the larval phase on connectivity: strong differences in the genetic structure of brooders and broadcasters in the Ophioderma longicauda species complex

Closely related species with divergent life history traits are excellent models to infer the role of such traits in genetic diversity and connectivity. Ophioderma longicauda is a brittle star species complex composed of different genetic clusters, including brooders and broadcasters. These species diverged very recently and some of them are sympatric and ecologically syntopic, making them particularly suitable to study the consequences of their trait differences. At the scale of the geographic distribution of the broadcasters (Mediterranean Sea and northeastern Atlantic), we sequenced the mitochondrial marker COI and genotyped an intron (i51) for 788 individuals. In addition, we sequenced 10 nuclear loci newly developed from transcriptome sequences, for six sympatric populations of brooders and broadcasters from Greece. At the large scale, we found a high genetic structure within the brooders (COI: 0.07 < F(ST) < 0.65) and no polymorphism at the nuclear locus i51. In contrast, the broadcasters displayed lower genetic structure (0 < F(ST) < 0.14) and were polymorphic at locus i51. At the regional scale, the multilocus analysis confirmed the contrasting genetic structure between species, with no structure in the broadcasters (global F(ST) < 0.001) and strong structure in the brooders (global F(ST) = 0.49), and revealed a higher genetic diversity in broadcasters. Our study showed that the lecithotrophic larval stage allows on average a 50-fold increase in migration rates, a 280-fold increase in effective size and a threefold to fourfold increase in genetic diversity. Our work, investigating complementary genetic markers on sympatric and syntopic taxa, highlights the strong impact of the larval phase on connectivity and genetic diversity.

Temporal genetic structure in a poecilogonous polychaete: the interplay of developmental mode and environmental stochasticity

Background

Temporal variation in the genetic structure of populations can be caused by multiple factors, including natural selection, stochastic environmental variation, migration, or genetic drift. In benthic marine species, the developmental mode of larvae may indicate a possibility for temporal genetic variation: species with dispersive planktonic larvae are expected to be more likely to show temporal genetic variation than species with benthic or brooded non-dispersive larvae, due to differences in larval mortality and dispersal ability. We examined temporal genetic structure in populations of Pygospio elegans, a poecilogonous polychaete with within-species variation in developmental mode. P. elegans produces either planktonic, benthic, or intermediate larvae, varying both among and within populations, providing a within-species test of the generality of a relationship between temporal genetic variation and larval developmental mode.

Results

In contrast to our expectations, our microsatellite analyses of P. elegans revealed temporal genetic stability in the UK population with planktonic larvae, whereas there was variation indicative of drift in temporal samples of the populations from the Baltic Sea, which have predominantly benthic and intermediate larvae. We also detected temporal variation in relatedness within these populations. A large temporal shift in genetic structure was detected in a population from the Netherlands, having multiple developmental modes. This shift could have been caused by local extiction due to extreme environmental conditions and (re)colonization by planktonic larvae from neighboring populations.

Conclusions

In our study of P. elegans, temporal genetic variation appears to be due to not only larval developmental mode, but also the stochastic environment of adults. Large temporal genetic shifts may be more likely in marine intertidal habitats (e.g. North Sea and Wadden Sea) which are more prone to environmental stochasticity than the sub-tidal Baltic habitats. Sub-tidal and/or brackish (less saline) habitats may support smaller P. elegans populations and these may be more susceptible to the effects of random genetic drift. Moreover, higher frequencies of asexual reproduction and the benthic larval developmental mode in these populations leads to higher relatedness and contributes to drift. Our results indicate that a general relationship between larval developmental mode and temporal genetic variation may not exist.

Shallow gene pools in the high intertidal: extreme loss of genetic diversity in viviparous sea stars (Parvulastra)

We document an extreme example of reproductive trait evolution that affects population genetic structure in sister species of Parvulastra cushion stars from Australia. Self-fertilization by hermaphroditic adults and brood protection of benthic larvae causes strong inbreeding and range-wide genetic poverty. Most samples were fixed for a single allele at nearly all nuclear loci; heterozygotes were extremely rare (0.18%); mitochondrial DNA sequences were more variable, but few populations shared haplotypes in common. Isolation-with-migration models suggest that these patterns are caused by population bottlenecks (relative to ancestral population size) and low gene flow. Loss of genetic diversity and low potential for dispersal between high-intertidal habitats may have dire consequences for extinction risk and potential for future adaptive evolution in response to climate and other selective agents.