Population genetics of a marine bivalve, Pinctada maxima, throughout the Indo-Australian Archipelago shows differentiation and decreased diversity at range limits

Genomic analysis of a reef-building coral, Acropora digitifera, reveals complex population structure and a migration network in the Nansei Islands, Japan

Understanding the structure and connectivity of coral populations is fundamental for developing marine conservation policies, especially in patchy environments such as archipelagos. The Nansei Islands, extending more than 1000 km in southwestern Japan, are characterized by high levels of biodiversity and endemism, supported by coral reefs, which make this region ideal for assessing genetic attributes of coral populations. In this study, we conducted population genomic analyses based on genome-wide, single-nucleotide polymorphisms (SNPs) of Acropora digitifera, a common species in the Nansei Islands. By merging newly obtained genome resequencing data with previously published data, we identified more than 4 million genome-wide SNPs in 303 colonies collected at 22 locations, with sequencing coverage ranging from 3.91× to 27.41×. While population structure analyses revealed genetic similarities between the southernmost and northernmost locations, separated by >1000 km, several subpopulations in intermediate locations suggested limited genetic admixture, indicating conflicting migration tendencies in the Nansei Islands. Although migration networks revealed a general tendency of northward migration along the Kuroshio Current, a substantial amount of southward migration was also detected, indicating important contributions of minor ocean currents to coral larval dispersal. Moreover, heterogeneity in the transition of effective population sizes among locations suggests different histories for individual subpopulations. The unexpected complexity of both past and present population dynamics in the Nansei Islands implies that heterogeneity of ocean currents and local environments, past and present, have influenced the population structure of this species, and similar unexpected population complexities may be expected for other marine species with similar reproductive modes.

Saving the sea cucumbers: Using population genomic tools to inform fishery and conservation management of the Fijian sandfish Holothuria (Metriatyla) scabra

The sea cucumber Holothuria (Metriatyla) scabra, known as sandfish, is a high-value tropical echinoderm central to the global bêche-de-mer (BDM) trade. This species has been heavily exploited across its natural range, with overharvesting and ineffective fishery management leaving stocks in the Pacific region heavily depleted. In Fiji, sandfish stocks have not recovered since a 1988 harvest ban, with surveys reporting declining populations and recruitment failure. Therefore, to inform fishery management policy for the wild sandfish resource and to guide hatchery-based restocking efforts, a high-resolution genomic audit of Fijian populations was carried out. A total of 6,896 selectively-neutral and 186 putatively-adaptive genome-wide SNPs (DArTseq) together with an independent oceanographic particle dispersal model were used to investigate genetic structure, diversity, signatures of selection, relatedness and connectivity in six wild populations. Three genetically distinct populations were identified with shallow but significant differentiation (average Fst = 0.034, p≤0.05), comprising (1) Lakeba island (Lau archipelago), (2) Macuata (Vanua Levu), and (3) individuals from Yasawa, Ra, Serua island and Kadavu comprising the final unit. Small reductions in allelic diversity were observed in marginal populations in eastern Fiji (overall mean A = 1.956 vs. Lau, A = 1.912 and Macuata, A = 1.939). Signatures of putative local adaptation were also discovered in individuals from Lakeba island, suggesting that they be managed as a discrete unit. An isolation-by-distance model of genetic structure for Fijian sandfish is apparent, with population fragmentation occurring towards the east. Hatchery-based production of juveniles is promising for stock replenishment, however great care is required during broodstock source population selection and juvenile releases into source areas only. The successful use of genomic data here has the potential to be applied to other sea cucumber species in Fiji, and other regions involved in the global BDM trade. While preliminary insights into the genetic structure and connectivity of sandfish in Fiji have been obtained, further local, regional and distribution-wide investigations are required to better inform conservation efforts, wild stock management and hatchery-based restocking interventions for this valuable invertebrate.

High gene flow in the silverlip pearl oyster Pinctada maxima between inshore and offshore sites near Eighty Mile Beach in Western Australia

An understanding of stock recruitment dynamics in fisheries is fundamental to successful management. Pinctada maxima is a bivalve mollusc widely distributed throughout the Indo-Pacific and is the main species targeted for cultured pearl and pearl shell production in Australia. Pearl production in Australia relies heavily on wild-caught individuals, the majority of which come from the Eighty Mile Beach region near Broome in Western Australia. In this study, we used a genotyping by sequencing approach to explore fine-scale patterns of genetic connectivity among inshore shallow and offshore deep populations of P. maxima near Eighty Mile Beach. Our results revealed high-levels of gene flow among inshore and offshore sites and no differences in genetic diversity between depths. Global estimates of genetic differentiation were low (F ST = 0.006) but significantly different from zero, and pairwise estimates of genetic differentiation among sites were significant in only 3% of comparisons. Moreover, Bayesian clustering detected no separation of inshore and offshore sample sites, and instead showed all samples to be admixed among sites, locations and depths. Despite an absence of any clear spatial clustering among sites, we identified a significant pattern of isolation by distance. In a dynamic environment like Eighty Mile Beach, genetic structure can change from year-to-year and successive dispersal and recruitment events over generations likely act to homogenize the population. Although we cannot rule out the null hypothesis of panmixia, our data indicate high levels of dispersal and connectivity among inshore and offshore fishing grounds.

Pleistocene expansion, anthropogenic pressure and ocean currents: Disentangling the past and ongoing evolutionary history of Patella aspera Röding, 1798 in the archipelago of Madeira

AIMS

Rising sea-level following the Last Glacial Maximum lead to fragmentation of coastal limpet populations between islands of the Archipelago of Madeira. This fragmentation is reinforced by recent heavy exploitation reducing effective population size on Madeira Island. We use the limpet P. aspera to understand how the role of processes at different time scales (i.e. changes in the sea level and overexploitation) can influence the genetic composition of an extant species, relating these processes to reproductive phenology and seasonal shifts in ocean currents.

LOCATION

Madeira Island, Porto Santo and Desertas (Archipelago of Madeira, NE Atlantic Ocean).

TAXON

The limpet Patella aspera.

METHODS

Twelve microsatellite genetic markers were used. A power analysis was used to evaluate the power of the microsatellite markers to detect a signal of population differentiation. Long-term past migrations were assessed using a Bayesian Markov Montecarlo approach in the software MIGRATE-n to estimate mutation-scaled migration rates (M = m/μ; m, probability of a lineage immigrating per generation; μ, mutation rate). Two scenarios were evaluated using an Approximate Bayesian Computation (ABC) in the software DIYABC 2.1 (i) Scenario 1: considered a population scenario from a reduced Ne at time t₃ to a higher Ne at time t₂; and (ii) Scenario 2 considering a reduction of Ne from a time t₃ to a time t₂.

RESULTS

Colonization of the archipelago by Portuguese settlers six centuries ago probably led to an important decrease in the genetic diversity of the species (N_e). Contemporary gene flow strongly support a pattern of high asymmetric connectivity explained by the reproductive phenology of the species and spatio-temporal seasonal changes in the ocean currents. Spatio-temporal reconstructions using Bayesian methods, including coalescent and Approximate Bayesian Computation (ABC) approaches, suggest changes in the migration patterns from highly symmetric to highly asymmetric connectivity with subtle population differentiation as consequence of post-glacial maximum sea level rise during the Holocene.

MAIN CONCLUSIONS

Our results suggest that anthropogenic activity could have had serious effects on the genetic diversity of heavily exploited littoral species since the end of the Pleistocene, probably accelerating in recent years.

Impacts of marine heatwaves on pearl oysters are alleviated following repeated exposure

Marine heatwaves (MHWs) have occurred with increasing duration, frequency and intensity in the past decade in the South China Sea, posing serious threats to marine ecosystems and fisheries. However, the impact of MHWs on marine bivalves - one of the most ecologically and economically important fauna in coastal ecosystems - remains largely unknown. Here, we investigated physiological responses of the pearl oyster, Pinctada maxima inhabiting a newly identified climate change hotspot (Beibu Gulf, South China Sea) to short-lasting and repeatedly-occurring MHWs scenarios. Following 3-day exposure to short-lasting MHWs scenarios with water temperature rapidly arising from 24 °C to 28 °C, 32 °C and 36 °C, respectively, mortality rates of pearl oysters increased, and especially they suffered 100% mortality at 36 °C. Activities of enzymes including acid phosphatase (ACP), alkaline phosphatase (AKP), glutathione (GSH) and level of malondialdehyde (MDA) increased significantly with increasing intensity and duration of MHWs, indicating thermal stress responses. When exposed to repeatedly-occurring MHWs scenarios, mortality rates of pearl oysters increased slightly, and thermal stress responses were alleviated, as exemplified by significant decreases in ACP, AKP, GSH and MDA activities compared with those during short-lasting MHWs scenarios, demonstrating the potential of P. maxima to acclimate rapidly to MHWs. These findings advance our understanding of how marine bivalves respond to MHWs scenarios varying in duration, frequency, and intensity.

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.

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.

Swept away: ocean currents and seascape features influence genetic structure across the 18,000 Km Indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

BACKGROUND

Genetic structure in many widely-distributed broadcast spawning marine invertebrates remains poorly understood, posing substantial challenges for their fishery management, conservation and aquaculture. Under the Core-Periphery Hypothesis (CPH), genetic diversity is expected to be highest at the centre of a species' distribution, progressively decreasing with increased differentiation towards outer range limits, as populations become increasingly isolated, fragmented and locally adapted. The unique life history characteristics of many marine invertebrates such as high dispersal rates, stochastic survival and variable recruitment are also likely to influence how populations are organised. To examine the microevolutionary forces influencing population structure, connectivity and adaptive variation in a highly-dispersive bivalve, populations of the black-lip pearl oyster Pinctada margaritifera were examined across its ~18,000 km Indo-Pacific distribution.

RESULTS

Analyses utilising 9,624 genome-wide SNPs and 580 oysters, discovered differing patterns of significant and substantial broad-scale genetic structure between the Indian and Pacific Ocean basins. Indian Ocean populations were markedly divergent (F st = 0.2534-0.4177, p < 0.001), compared to Pacific Ocean oysters, where basin-wide gene flow was much higher (F st = 0.0007-0.1090, p < 0.001). Partitioning of genetic diversity (hierarchical AMOVA) attributed 18.1% of variance between ocean basins, whereas greater proportions were resolved within samples and populations (45.8% and 35.7% respectively). Visualisation of population structure at selectively neutral loci resolved three and five discrete genetic clusters for the Indian and Pacific Oceans respectively. Evaluation of genetic structure at adaptive loci for Pacific populations (89 SNPs under directional selection; F st = 0.1012-0.4371, FDR = 0.05), revealed five clusters identical to those detected at neutral SNPs, suggesting environmental heterogeneity within the Pacific. Patterns of structure and connectivity were supported by Mantel tests of isolation by distance (IBD) and independent hydrodynamic particle dispersal simulations.

CONCLUSIONS

It is evident that genetic structure and connectivity across the natural range of P. margaritifera is highly complex, and produced by the interaction of ocean currents, IBD and seascape features at a broad scale, together with habitat geomorphology and local adaptation at regional levels. Overall population organisation is far more elaborate than generalised CPH predictions, however valuable insights for regional fishery management, and a greater understanding of range-wide genetic structure in a highly-dispersive marine invertebrate have been gained.

Unravelling the effects of gene flow and selection in highly connected populations of the silver-lip pearl oyster (Pinctada maxima)

Many marine organisms often display weak levels of population genetic structuring as a result of both environmental characteristics (e.g., ocean currents) and life history traits (e.g., widely dispersed planktonic larval stages) maintaining high levels of gene flow. This can lead to the assumption that these organisms can be managed as a single stock based on high levels of population connectivity. However, this neglects to account for other micro-evolutionary forces such as selection, which also shape these populations. This study utilizes 1130 genome-wide SNP loci to unravel the effects of gene flow and selection shaping three highly connected populations of the silver-lip pearl oyster (Pinctada maxima) in the ecologically and economically important Indo-Pacific region (Aru, Bali, and West Papua). Twenty-two loci under directional selection were identified amongst the populations, providing further supporting evidence of strong local adaptation (i.e., G×E effects) among populations in this region. Global Fst values for directional outliers (0.348) were up to eight times greater than for neutral markers (0.043). Pairwise Fst comparisons between Aru and Bali revealed the largest directional differences (0.488), while Bali and West Papua had the least (0.062). Unrooted neighbour-joining (NJ) distance trees and genetic diversity indices of directional outliers revealed that individuals from Bali and West Papua had reduced allelic variation (MAFavg=0.144, Ho=0.238 and MAFavg=0.232, Ho=0.369, respectively) compared to Aru (MAFavg=0.292, Ho=0.412). This indicates that directional selection is most likely acting upon genetic variation within the Bali and West Papua populations. NJ distance trees, discriminant analysis of principal components, and Fst analyses of directional outliers revealed two divergent groups ("Bali/West Papua"; "Aru") that had previously gone unrecognized. This study not only illustrates that relatively strong local adaptive forces are occurring despite high gene flow, but identifies the populations that are most likely experiencing selection. Additionally, this study highlights the need to understand all micro-evolutionary forces acting on populations when resolving stock structure.

netview p: a network visualization tool to unravel complex population structure using genome-wide SNPs

Network-based approaches are emerging as valuable tools for the analysis of complex genetic structure in wild and captive populations. netview p combines data quality control with the construction of population networks through mutual k-nearest neighbours thresholds applied to genome-wide SNPs. The program is cross-platform compatible, open-source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs. The pipeline was used for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive silver-lipped pearl oysters (n = 415, SNPs = 1107), wild populations of the European hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and grey wolves from North America (n = 239, SNPs = 78 255). The population networks effectively visualize large- and fine-scale genetic structure within and between populations, including family-level structure and relationships. netview p comprises a network-based addition to other population analysis tools and provides user-friendly access to a complex network analysis pipeline through implementation in python.

Genome-wide SNP validation and mantle tissue transcriptome analysis in the silver-lipped pearl oyster, Pinctada maxima

Pearl oysters are not only farmed for their gemstone quality pearls worldwide, but they are also becoming important model organisms for investigating genetic mechanisms of biomineralisation. Despite their economic and scientific significance, limited genomic resources are available for this important group of bivalves, hampering investigations into identifying genes that regulate important pearl quality traits and unique biological characteristics (i.e. biomineralisation). The silver-lipped pearl oyster, Pinctada maxima, is one species where there is interest in understanding genes that regulate commercially important pearl traits, but presently, there is a dearth of genomic information. The objective of this study was to develop and validate a large number of type I genome-wide single nucleotide polymorphisms (SNPs) for P. maxima suitable for high-throughput genotyping. In addition, sequence annotations and Gene Ontology terms were assigned to a large mantle tissue 454 expressed sequence tag assembly (96,794 contigs) and information on known bivalve biomineralisation genes was incorporated into SNP discovery. The SNP discovery effort resulted in the de novo identification of 172,625 SNPs, of which 9,108 were identified as high value [minor allele frequency (MAF)≥ 0.15, read depth  ≥ 8]. Validation of 2,782 of these SNPs using Illumina iSelect Infinium genotyping technology returned some of the highest assay conversion (86.6 %) and validation (59.9 %; mean MAF 0.28) rates observed in aquaculture species to date. Genomic resources presented here will be pivotal to future research investigating the biological mechanisms behind biomineralisation and will form a strong foundation for genetic selective breeding programs in the P. maxima pearling industry.

Population expansion and genetic structure in Carcharhinus brevipinna in the southern Indo-Pacific

BACKGROUND

Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa.

METHODS AND FINDINGS

Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and F ST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717-0.03508, p values ≤ 0.0013; pairwise F ST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia's east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).

CONCLUSIONS

The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

Longitudinal differentiation among pelagic populations in a planktic foraminifer

Evolutionary processes in marine plankton have been assumed to be dependent on the oceanic circulation system, which transports plankton between populations in marine surface waters. Gene flow facilitated by oceanic currents along longitudinal gradients may efficiently impede genetic differentiation of pelagic populations in the absence of confounding marine environmental effects. However, how responsible oceanic currents are for the geographic distribution and dispersal of plankton is poorly understood. We examined the phylogeography of the planktic foraminifer Pulleniatina obliquiloculata in the Indo-Pacific Warm Pool (IPWP) by using partial small subunit ribosomal DNA (SSU rDNA) sequences. We found longitudinal clines in the frequencies of three distinct genetic types in the IPWP area. These frequencies were correlated with environmental factors that are characteristic of three water masses in the IPWP. Noteworthy, populations inhabiting longitudinally distant water masses at the Pacific and Indian sides of the IPWP were genetically different, despite transportation of individuals via oceanic currents. These results demonstrate that populations of pelagic plankton have diverged genetically among different water masses within a single climate zone. Changes of the oceanic circulation system could have impacted the geographic patterns of dispersal and divergence of pelagic plankton.

In silico whole-genome EST analysis reveals 2322 novel microsatellites for the silver-lipped pearl oyster, Pinctada maxima

Molecular stock improvement techniques such as marker assisted selection have great potential in accelerating selective breeding programmes for animal production industries. However, the discovery and application of trait/marker associations usually requires a large number of genome-wide polymorphic loci. Here, we present 2322 unique microsatellites for the silver-lipped pearl oyster, Pinctada maxima, a species of aquaculture importance throughout the Indo-Australian Archipelago for production of the highly valued South Sea pearl. More than 1.2 million Roche 454 expressed sequence tag (EST) reads were screened for microsatellite repeat motifs. A total of 12,604 sequences contained either a di, tri, tetra, penta or hexa microsatellite repeat motif (n ≥ 6), with 6435 of these sequences having sufficient flanking regions for primer development. All identified microsatellites with designed primers were condensed into 2322 unique clusters (i.e., unique loci) of which 360 were shown to be polymorphic based on multiple sequence reads with different repeat motifs. Genotyping of five microsatellite loci demonstrated that in silico evaluation of polymorphism levels was a very useful method for identification of polymorphic loci, with the variation uncovered being a lower bound. Gene Ontology annotations of sequences containing microsatellites suggest that most are derived from a diverse array of unique genes. This EST derived microsatellite database will be a valuable resource for future studies in genetic map construction, diversity analysis, quantitative trait loci analysis, association mapping and marker assisted selection, not only for P. maxima, but also closely related species within the genus Pinctada.

DNA barcoding of marine metazoa

More than 230,000 known species representing 31 metazoan phyla populate the world's oceans. Perhaps another 1,000,000 or more species remain to be discovered. There is reason for concern that species extinctions may out-pace discovery, especially in diverse and endangered marine habitats such as coral reefs. DNA barcodes (i.e., short DNA sequences for species recognition and discrimination) are useful tools to accelerate species-level analysis of marine biodiversity and to facilitate conservation efforts. This review focuses on the usual barcode region for metazoans: a approximately 648 base-pair region of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Barcodes have also been used for population genetic and phylogeographic analysis, identification of prey in gut contents, detection of invasive species, forensics, and seafood safety. More controversially, barcodes have been used to delimit species boundaries, reveal cryptic species, and discover new species. Emerging frontiers are the use of barcodes for rapid and increasingly automated biodiversity assessment by high-throughput sequencing, including environmental barcoding and the use of barcodes to detect species for which formal identification or scientific naming may never be possible.