Genetic structure of a commercially exploited bivalve, the great scallop Pecten maximus, along the European coasts

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus)

The genomic landscape of divergence—the distribution of differences among populations or species across the genome—is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12–22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

Genetic homogeneity of the critically endangered fan mussel, Pinna nobilis, throughout lagoons of the Gulf of Lion (North-Western Mediterranean Sea)

The fan mussel, Pinna nobilis, endemic to the Mediterranean Sea, is a critically endangered species facing mass mortality events in almost all of its populations, following the introduction of the parasite Haplosporidium pinnae. Such a unique pandemic in a marine organism, which spreads rapidly and with mortality rates reaching up to 100%, could lead to the potential extinction of the species. Only few regions, involving lagoon habitats, remain healthy throughout the entire Mediterranean Sea. This study describes the genetic structure of P. nobilis across the Gulf of Lion, including confined locations such as lagoons and ports. A total of 960 samples were collected among 16 sites distributed at 8 localities, and then genotyped using 22 microsatellite markers. Genetic diversity was high in all sites with mean allele numbers ranging between 10 and 14.6 and with observed heterozygosities (Ho) between 0.679 and 0.704. No genetic differentiation could be identified (F_ST ranging from 0.0018 to 0.0159) and the percentages of related individuals were low and similar among locations (from 1.6 to 6.5%). Consequently, all fan mussels, over the entire coastline surveyed, including those in the most geographically isolated areas, belong to a large genetically homogeneous population across the Gulf of Lion. Considering the ongoing mass mortality context, this result demonstrates that almost all of the genetic diversity of P. nobilis populations is still preserved even in isolated lagoons, which might represent a refuge habitat for the future of the species.

Evolution at two time frames: ancient structural variants involved in post-glacial divergence of the European plaice (Pleuronectes platessa)

Changing environmental conditions can lead to population diversification through differential selection on standing genetic variation. Structural variant (SV) polymorphisms provide examples of ancient alleles that in time become associated with novel environmental gradients. The European plaice (Pleuronectes platessa) is a marine flatfish showing large allele-frequency differences at two putative SVs associated with environmental variation. In this study, we explored the contribution of these SVs to population structure across the North East Atlantic. We compared genome-wide population structure using sets of RAD-sequencing SNPs with the spatial structure of the SVs. We found that in contrast to the rest of the genome, the SVs were only weakly associated with an isolation-by-distance pattern. Indeed, both SVs showed important variation in haplogroup frequencies, with the same haplogroup increasing both along the salinity gradient of the Baltic Sea, and found in high frequency in the northern-range margin of the Atlantic. Phylogenetic analyses suggested that the SV alleles are much older than the age of the Baltic Sea itself. These results suggest that the SVs are older than the age of the environmental gradients with which they currently co-vary. Altogether, our results suggest that the plaice SVs were shaped by evolutionary processes occurring at two time frames, firstly following their origin, ancient spread and maintenance in the ancestral populations, and secondly related to their current association with more recently formed environmental gradients such as those found in the North Sea-Baltic Sea transition zone.

Metal pollution as a potential threat to shell strength and survival in marine bivalves

Marine bivalve molluscs, such as scallops, mussels and oysters, are crucial components of coastal ecosystems, providing a range of ecosystem services, including a quarter of the world's seafood. Unfortunately, coastal marine areas often suffer from high levels of metals due to dumping and disturbance of contaminated material. We established that increased levels of metal pollution (zinc, copper and lead) in sediments near the Isle of Man, resulting from historical mining, strongly correlated with significant weakening of shell strength in king scallops, Pecten maximus. This weakness increased mortality during fishing and left individuals more exposed to predation. Comparative structural analysis revealed that shells from the contaminated area were thinner and exhibited a pronounced mineralisation disruption parallel to the shell surface within the foliated region of both the top and bottom valves. Our data suggest that these disruptions caused reduced fracture strength and hence increased mortality, even at subcritical contamination levels with respect to current international standards. This hitherto unreported effect is important since such non-apical responses rarely feed into environmental quality assessments, despite potentially significant implications for the survival of organisms exposed to contaminants. Hence our findings highlight the impact of metal pollution on shell mineralisation in bivalves and urge a reappraisal of currently accepted critical contamination levels.

The gene-rich genome of the scallop Pecten maximus

BACKGROUND

The king scallop, Pecten maximus, is distributed in shallow waters along the Atlantic coast of Europe. It forms the basis of a valuable commercial fishery and plays a key role in coastal ecosystems and food webs. Like other filter feeding bivalves it can accumulate potent phytotoxins, to which it has evolved some immunity. The molecular origins of this immunity are of interest to evolutionary biologists, pharmaceutical companies, and fisheries management.

FINDINGS

Here we report the genome assembly of this species, conducted as part of the Wellcome Sanger 25 Genomes Project. This genome was assembled from PacBio reads and scaffolded with 10X Chromium and Hi-C data. Its 3,983 scaffolds have an N50 of 44.8 Mb (longest scaffold 60.1 Mb), with 92% of the assembly sequence contained in 19 scaffolds, corresponding to the 19 chromosomes found in this species. The total assembly spans 918.3 Mb and is the best-scaffolded marine bivalve genome published to date, exhibiting 95.5% recovery of the metazoan BUSCO set. Gene annotation resulted in 67,741 gene models. Analysis of gene content revealed large numbers of gene duplicates, as previously seen in bivalves, with little gene loss, in comparison with the sequenced genomes of other marine bivalve species.

CONCLUSIONS

The genome assembly of P. maximus and its annotated gene set provide a high-quality platform for studies on such disparate topics as shell biomineralization, pigmentation, vision, and resistance to algal toxins. As a result of our findings we highlight the sodium channel gene Nav1, known to confer resistance to saxitoxin and tetrodotoxin, as a candidate for further studies investigating immunity to domoic acid.

Tracking genetic diversity in a large-scale oyster restoration program: effects of hatchery propagation and initial characterization of diversity on restored vs. wild reefs

The release of hatchery-propagated fish and shellfish is occurring on a global scale, but the genetic impacts of these practices are often not fully understood and rarely monitored. Slow recovery of depleted eastern oyster populations in the Chesapeake Bay, USA has prompted a hatchery-based restoration program focused in the Choptank River, Maryland consisting of the mass release of hatchery-produced juveniles from local, wild broodstock. To evaluate potential genetic effects of this program, we (1) examined changes in genetic diversity (allelic richness, heterozygosity) and the effective number of breeders (Nb) over the hatchery production cycle with microsatellite-based parentage of natural, mass- and controlled-spawned cohorts, and (2) compared genetic diversity and effective population size (Ne) of a restored reef to wild source populations. Mass-spawned cohorts showed high variance in reproductive contribution, particularly among males, leading to a 45% average reduction in Nb from spawning adult numbers and higher relatedness-lower magnitude reductions in heterozygosity and significant reductions in allelic richness were also observed. While controlled-spawns (single-male fertilizations of pooled eggs) reduced male variance, overall reproductive variance (Vk) remained high. Finally, oysters sampled from a restored reef displayed comparable Ne, genetic diversity, and relatedness to samples from wild populations, with no significant genetic differentiation among them. Overall, the hatchery-based results and initial field-based population genetic analyses suggest that despite reductions in diversity from parents to offspring owing to high Vk, enhancement with rotated, wild broodstock appears to have maintained genetic diversity in a restored reef population compared to proximal wild populations.

RAD sequencing sheds new light on the genetic structure and local adaptation of European scallops and resolves their demographic histories

Recent developments in genomics are advancing our understanding of the processes shaping population structure in wild organisms. In particular, reduced representation sequencing has facilitated the generation of dense genetic marker datasets that provide greater power for resolving population structure, investigating the role of selection and reconstructing demographic histories. We therefore used RAD sequencing to study the great scallop Pecten maximus and its sister species P. jacobeus along a latitudinal cline in Europe. Analysis of 219 samples genotyped at 82,439 single nucleotide polymorphisms clearly resolved an Atlantic and a Norwegian group within P. maximus as well as P. jacobeus, in support of previous studies. Fine-scale structure was also detected, including pronounced differences involving Mulroy Bay in Ireland, where scallops are commercially cultured. Furthermore, we identified a suite of 279 environmentally associated loci that resolved a contrasting phylogenetic pattern to the remaining neutral loci, consistent with ecologically mediated divergence. Finally, demographic inference provided support for the two P. maximus groups having diverged during the last glacial maximum and subsequently expanded, whereas P. jacobeus diverged around 95,000 generations ago and experienced less pronounced expansion. Our results provide an integrative perspective on the factors shaping genome-wide differentiation in a commercially important marine invertebrate.

Coupling molecular data and experimental crosses sheds light about species delineation: a case study with the genus Ciona

Molecular studies sometimes reveal evolutionary divergence within accepted species. Such findings can initiate taxonomic revision, as exemplified in the formerly recognized species Ciona intestinalis. While an increasing number of studies have examined the ecology, reproductive barriers and genetics of C. intestinalis and C. robusta, there are still much uncertainties regarding other species of this genus. Using experimental crosses and mitochondrial data, we investigated the evolutionary relationships among four native and introduced Ciona spp., found in sympatry in the Mediterranean Sea or English Channel. Outcome of 62 bi-parental reciprocal crosses between C. intestinalis, C. robusta, C. roulei and C. edwardsi showed that C. edwardsi is reproductively isolated from the other taxa, which is in agreement with its distinct location in the phylogenetic tree. Conversely, hybrids are easily obtained in both direction when crossing C. intestinalis and C. roulei, reinforcing the hypothesis of two genetically differentiated lineages but likely being from a same species. Altogether, this study sheds light on the evolutionary relationship in this complex genus. It also calls for further investigation notably based on genome-wide investigation to better describe the evolutionary history within the genus Ciona, a challenging task in a changing world where biological introductions are shuffling species distribution.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher Fst values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher F_st values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

Stock enhancement or sea ranching? Insights from monitoring the genetic diversity, relatedness and effective population size in a seeded great scallop population (Pecten maximus)

The mass release of hatchery-propagated stocks raises numerous questions concerning its efficiency in terms of local recruitment and effect on the genetic diversity of wild populations. A seeding program, consisting of mass release of hatchery-produced juveniles in the local naturally occurring population of great scallops (Pecten maximus L.), was initiated in the early 1980s in the Bay of Brest (France). The present study aims at evaluating whether this seeding program leads to actual population enhancement, with detectable effects on genetic diversity and effective population size, or consists of sea ranching with limited genetic consequences on the wild stock. To address this question, microsatellite-based genetic monitoring of three hatchery-born and naturally recruited populations was conducted over a 5-year period. Results showed a limited reduction in allelic richness but a strong alteration of allelic frequencies in hatchery populations, while genetic diversity appeared very stable over time in the wild populations. A temporal increase in relatedness was observed in both cultured stock and wild populations. Effective population size (Ne) estimates were low and variable in the wild population. Moreover, the application of the Ryman-Laikre model suggested a high contribution of hatchery-born scallops to the reproductive output of the wild population. Overall, the data suggest that the main objective of the seeding program, which is stock enhancement, is fulfilled. Moreover, gene flow from surrounding populations and/or the reproductive input of undetected sub-populations within the bay may buffer the Ryman-Laikre effect and ensure the retention of the local genetic variability.