Unexpected fine-scale population structure in a broadcast-spawning Antarctic marine mollusc

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

Heterogeneous microgeographic genetic structure of the common cockle (Cerastoderma edule) in the Northeast Atlantic Ocean: biogeographic barriers and environmental factors

Knowledge of genetic structure at the finest level is essential for the conservation of genetic resources. Despite no visible barriers limiting gene flow, significant genetic structure has been shown in marine species. The common cockle (Cerastoderma edule) is a bivalve of great commercial and ecological value inhabiting the Northeast Atlantic Ocean. Previous population genomics studies demonstrated significant structure both across the Northeast Atlantic, but also within small geographic areas, highlighting the need to investigate fine-scale structuring. Here, we analysed two geographic areas that could represent opposite models of structure for the species: (1) the SW British Isles region, highly fragmented due to biogeographic barriers, and (2) Galicia (NW Spain), a putative homogeneous region. A total of 9250 SNPs genotyped by 2b-RAD on 599 individuals from 22 natural beds were used for the analysis. The entire SNP dataset mostly confirmed previous observations related to genetic diversity and differentiation; however, neutral and divergent SNP outlier datasets enabled disentangling physical barriers from abiotic environmental factors structuring both regions. While Galicia showed a homogeneous structure, the SW British Isles region was split into four reliable genetic regions related to oceanographic features and abiotic factors, such as sea surface salinity and temperature. The information gathered supports specific management policies of cockle resources in SW British and Galician regions also considering their particular socio-economic characteristics; further, these new data will be added to those recently reported in the Northeast Atlantic to define sustainable management actions across the whole distribution range of the species.

Seascape genomics of coastal bottlenose dolphins along strong gradients of temperature and salinity

Heterogeneous seascapes and strong environmental gradients in coastal waters are expected to influence adaptive divergence, particularly in species with large population sizes where selection is expected to be highly efficient. However, these influences might also extend to species characterized by strong social structure, natal philopatry and small home ranges. We implemented a seascape genomic study to test this hypothesis in Indo-Pacific bottlenose dolphins (Tursiops aduncus) distributed along the environmentally heterogeneous coast of southern Australia. The datasets included oceanographic and environmental variables thought to be good predictors of local adaptation in dolphins and 8,081 filtered single nucleotide polymorphisms (SNPs) genotyped for individuals sampled from seven different bioregions. From a neutral perspective, population structure and connectivity of the dolphins were generally influenced by habitat type and social structuring. Genotype-environment association analysis identified 241 candidate adaptive loci and revealed that sea surface temperature and salinity gradients influenced adaptive divergence in these animals at both large- (1,000s km) and fine-scales (<100 km). Enrichment analysis and annotation of candidate genes revealed functions related to sodium-activated ion transport, kidney development, adipogenesis and thermogenesis. The findings of spatial adaptive divergence and inferences of putative physiological adaptations challenge previous suggestions that marine megafauna is most likely to be affected by environmental and climatic changes via indirect, trophic effects. Our work contributes to conservation management of coastal bottlenose dolphins subjected to anthropogenic disturbance and to efforts of clarifying how seascape heterogeneity influences adaptive diversity and evolution in small cetaceans.

Parallel shape divergence between ecotypes of the limpet Nacella concinna along the Antarctic Peninsula: a new model species for parallel evolution?

Parallel phenotypic divergence is the independent differentiation between phenotypes of the same lineage or species occupying ecologically similar environments in different populations. We tested in the Antarctic limpet Nacella concinna the extent of parallel morphological divergence in littoral and sublittoral ecotypes throughout its distribution range. These ecotypes differ in morphological, behavioural and physiological characteristics. We studied the lateral and dorsal outlines of shells and the genetic variation of the mitochondrial gene Cytochrome Oxidase subunit I from both ecotypes in 17 sample sites along more than 2,000 km. The genetic data indicate that both ecotypes belong to a single evolutionary lineage. The magnitude and direction of phenotypic variation differ between ecotypes across sample sites; completely parallel ecotype-pairs (i.e., they diverge in the same magnitude and in the same direction) were detected in 84.85% of lateral and 65.15% in dorsal view comparisons. Besides, specific traits (relative shell height, position of shell apex, and elliptical/pear-shape outline variation) showed high parallelism. We observed weak morphological covariation between the two shape shell views, indicating that distinct evolutionary forces and environmental pressures could be acting on this limpet shell shape. Our results demonstrate there is a strong parallel morphological divergence pattern in N. concinna along its distribution, making this Antarctic species a suitable model for the study of different evolutionary forces shaping the shell evolution of this limpet.

Sweepstake reproductive success and collective dispersal produce chaotic genetic patchiness in a broadcast spawner

A long-standing paradox of marine populations is chaotic genetic patchiness (CGP), temporally unstable patterns of genetic differentiation that occur below the geographic scale of effective dispersal. Several mechanisms are hypothesized to explain CGP including natural selection, spatiotemporal fluctuations in larval source populations, self-recruitment, and sweepstake reproduction. Discriminating among them is extremely difficult but is fundamental to understanding how marine organisms reproduce and disperse. Here, we report a notable example of CGP in the Antarctic limpet, an unusually tractable system where multiple confounding explanations can be discounted. Using population genomics, temporally replicated sampling, surface drifters, and forward genetic simulations, we show that CGP likely arises from an extreme sweepstake event together with collective larval dispersal, while selection appears to be unimportant. Our results illustrate the importance of neutral demographic forces in natural populations and have important implications for understanding the recruitment dynamics, population connectivity, local adaptation, and resilience of marine populations.

Divergence time and species delimitation of microbivalves in the Southern Ocean: the case of Kidderia species

The systematics of Subantarctic and Antarctic near-shore marine benthic invertebrates requires major revision and highlights the necessity to incorporate additional sources of information in the specimen identification chart in the Southern Ocean (SO). In this study, we aim to improve our understanding of the biodiversity of Kidderia (Dall 1876) through molecular and morphological comparisons of Antarctic and Subantarctic taxa. The microbivalves of the genus Kidderia are small brooding organisms that inhabit intertidal and shallow subtidal rocky ecosystems. This genus represents an interesting model to test the vicariance and dispersal hypothesis in the biogeography of the SO. However, the description of Kidderia species relies on a few morphological characters and biogeographic records that raise questions about the true diversity in the group. Here we will define the specimens collected with genetic tools, delimiting their respective boundaries across provinces of the SO, validating the presence of two species of Kidderia. Through the revision of taxonomic issues and species delimitation, it was possible to report that the Antarctic species is Kidderia subquadrata and the species recorded in the Subantarctic islands Diego Ramirez, South Georgia and the Kerguelen Archipelago is Kidderia minuta. The divergence time estimation suggests the origin and diversification of Kidderia lineages are related to historical vicariant processes probably associated with the separation of the continental landmasses close to the late Eocene.

Distinct fungal communities associated with different organs of the mangrove Sonneratia alba in the Malay Peninsula

Mangrove forests are key tropical marine ecosystems that are rich in fungi, but our understanding of fungal communities associated with mangrove trees and their various organs remains limited because much of the diversity lies within the microbiome. In this study, we investigated the fungal communities associated with the mangrove tree Sonneratia alba throughout Peninsular Malaysia and Singapore. At each sampling location, we collected leaves, fruits, pneumatophores and sediment samples and performed amplicon sequencing of the ribosomal internal transcribed spacer 1 to characterise the associated communities. Results show distinct fungal communities at each sampled location with further differentiation according to the plant part. We find a significant distance decay of similarity, particularly for sediment samples due to the greater variability of sediment environments relative to the more stable fungal habitats provided by living plant organs. We are able to assign taxonomy to the majority of sequences from leaves and fruits, but a much larger portion of the sequences recovered from pneumatophores and sediment samples could not be identified. This pattern underscores the limited mycological research performed in marine environments and demonstrates the need for a concerted research effort on multiple species to fully characterise the coastal microbiome and its role in the functioning of marine ecosystems.

Closely related octopus species show different spatial genetic structures in response to the Antarctic seascape

Determining whether comparable processes drive genetic divergence among marine species is relevant to molecular ecologists and managers alike. Sympatric species with similar life histories might be expected to show comparable patterns of genetic differentiation and a consistent influence of environmental factors in shaping divergence. We used microsatellite loci to quantify genetic differentiation across the Scotia Arc in three species of closely related benthic octopods, Pareledone turqueti, P. charcoti, and Adelieledone polymorpha. The relative importance of environmental factors (latitude, longitude, depth, and temperature) in shaping genetic structure was investigated when significant spatial genetic structure was uncovered. Isolated populations of P. turqueti and A. polymorpha at these species' range margins were genetically different to samples close to mainland Antarctica; however, these species showed different genetic structures at a regional scale. Samples of P. turqueti from the Antarctic Peninsula, Elephant Island, and Signy Island were genetically different, and this divergence was associated primarily with sample collection depth. By contrast, weak or nonsignificant spatial genetic structure was evident across the Antarctic Peninsula, Elephant Island, and Signy Island region for A. polymorpha, and slight associations between population divergence and temperature or depth (and/or longitude) were detected. Pareledone charcoti has a limited geographic range, but exhibited no genetic differentiation between samples from a small region of the Scotia Arc (Elephant Island and the Antarctic Peninsula). Thus, closely related species with similar life history strategies can display contrasting patterns of genetic differentiation depending on spatial scale; moreover, depth may drive genetic divergence in Southern Ocean benthos.

Contrasting patterns of population connectivity between regions in a commercially important mollusc Haliotis rubra: integrating population genetics, genomics and marine LiDAR data

Estimating contemporary genetic structure and population connectivity in marine species is challenging, often compromised by genetic markers that lack adequate sensitivity, and unstructured sampling regimes. We show how these limitations can be overcome via the integration of modern genotyping methods and sampling designs guided by LiDAR and SONAR data sets. Here we explore patterns of gene flow and local genetic structure in a commercially harvested abalone species (Haliotis rubra) from southeastern Australia, where the viability of fishing stocks is believed to be dictated by recruitment from local sources. Using a panel of microsatellite and genomewide SNP markers, we compare allele frequencies across a replicated hierarchical sampling area guided by bathymetric LiDAR imagery. Results indicate high levels of gene flow and no significant genetic structure within or between benthic reef habitats across 1400 km of coastline. These findings differ to those reported for other regions of the fishery indicating that larval supply is likely to be spatially variable, with implications for management and long-term recovery from stock depletion. The study highlights the utility of suitably designed genetic markers and spatially informed sampling strategies for gaining insights into recruitment patterns in benthic marine species, assisting in conservation planning and sustainable management of fisheries.

Octopus vulgaris (Cuvier, 1797) in the Mediterranean Sea: Genetic Diversity and Population Structure

The common octopus, Octopus vulgaris Cuvier 1797, is a largely exploited cephalopod species in the Mediterranean Sea and the Atlantic Ocean, as well as along the coasts of Africa, Brazil and Japan, where its taxonomic identity is still debated. The assessment of its genetic structure is a pressing need to correctly manage the resource and to avoid overfishing and collapsing of local stocks. Here we analysed genetic variation and population structure of O. vulgaris using thirteen microsatellite loci in seven sampling localities from the Mediterranean Sea and one from the Atlantic Ocean. We also used a DNA barcoding approach by COI gene fragment to understand the phylogenetic relationships among the specimens here investigated and the ones whose sequences are available in literature. Our results reveal high levels of allelic richness and moderate heterozygosity in all samples investigated, and a pronounced differentiation of the Atlantic and Sicilian specimens. This latter aspect seems to support the isolation of the biota within the Strait of Messina. A certain degree of differentiation was detected among the other geographic samples within the Mediterranean Sea, which is more compatible with an island model than isolation by distance. The occurrence of null alleles affected more genetic diversity indices than population structure estimations. This study provides new insights about the genetic diversity and structure of O. vulgaris in the area of interest, which can be used as guidelines for a fisheries management perspective.

The Story of a Hitchhiker: Population Genetic Patterns in the Invasive Barnacle Balanus(Amphibalanus) improvisus Darwin 1854

Understanding the ecological and evolutionary forces that determine the genetic structure and spread of invasive species is a key component of invasion biology. The bay barnacle, Balanus improvisus (= Amphibalanus improvisus), is one of the most successful aquatic invaders worldwide, and is characterised by broad environmental tolerance. Although the species can spread through natural larval dispersal, human-mediated transport through (primarily) shipping has almost certainly contributed to the current global distribution of this species. Despite its worldwide distribution, little is known about the phylogeography of this species. Here, we characterize the population genetic structure and model dispersal dynamics of the barnacle B. improvisus, and describe how human-mediated spreading via shipping as well as natural larval dispersal may have contributed to observed genetic variation. We used both mitochondrial DNA (cytochrome c oxidase subunit I: COI) and nuclear microsatellites to characterize the genetic structure in 14 populations of B. improvisus on a global and regional scale (Baltic Sea). Genetic diversity was high in most populations, and many haplotypes were shared among populations on a global scale, indicating that long-distance dispersal (presumably through shipping and other anthropogenic activities) has played an important role in shaping the population genetic structure of this cosmopolitan species. We could not clearly confirm prior claims that B. improvisus originates from the western margins of the Atlantic coasts; although there were indications that Argentina could be part of a native region. In addition to dispersal via shipping, we show that natural larval dispersal may play an important role for further colonisation following initial introduction.

Phylogeography of the Rock Shell Thais clavigera (Mollusca): Evidence for Long-Distance Dispersal in the Northwestern Pacific

The present-day genetic structure of a species reflects both historical demography and patterns of contemporary gene flow among populations. To precisely understand how these factors shape current population structure of the northwestern (NW) Pacific marine gastropod, Thais clavigera, we determined the partial nucleotide sequences of the mitochondrial COI gene for 602 individuals sampled from 29 localities spanning almost the whole distribution of T. clavigera in the NW Pacific Ocean (~3,700 km). Results from population genetic and demographic analyses (AMOVA, ΦST-statistics, haplotype networks, Tajima's D, Fu's FS, mismatch distribution, and Bayesian skyline plots) revealed a lack of genealogical branches or geographical clusters, and a high level of genetic (haplotype) diversity within each of studied population. Nevertheless, low but significant genetic structuring was detected among some geographical populations separated by the Changjiang River, suggesting the presence of geographical barriers to larval dispersal around this region. Several lines of evidence including significant negative Tajima's D and Fu's FS statistics values, the unimodally shaped mismatch distribution, and Bayesian skyline plots suggest a population expansion at marine isotope stage 11 (MIS 11; 400 ka), the longest and warmest interglacial interval during the Pleistocene epoch. The lack of genetic structure among the great majority of the NW Pacific T. clavigera populations may be attributable to high gene flow by current-driven long-distance dispersal of prolonged planktonic larval phase of this species.

Hierarchical population genetic structure in a direct developing antarctic marine invertebrate

Understanding the relationship between life-history variation and population structure in marine invertebrates is not straightforward. This is particularly true of polar species due to the difficulty of obtaining samples and a paucity of genomic resources from which to develop nuclear genetic markers. Such knowledge, however, is essential for understanding how different taxa may respond to climate change in the most rapidly warming regions of the planet. We therefore used over two hundred polymorphic Amplified Fragment Length Polymorphisms (AFLPs) to explore population connectivity at three hierachical spatial scales in the direct developing Antarctic topshell Margarella antarctica. To previously published data from five populations spanning a 1500 km transect along the length of the Western Antarctic Peninsula, we added new AFLP data for four populations separated by up to 6 km within Ryder Bay, Adelaide Island. Overall, we found a nonlinear isolation-by-distance pattern, suggestive of weaker population structure within Ryder Bay than is present over larger spatial scales. Nevertheless, significantly positive F_st values were obtained in all but two of ten pairwise population comparisons within the bay following Bonferroni correction for multiple tests. This is in contrast to a previous study of the broadcast spawner Nacella concinna that found no significant genetic differences among several of the same sites. By implication, the topshell's direct-developing lifestyle may constrain its ability to disperse even over relatively small geographic scales.