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.

Transcriptomic features of Pecten maximus oocyte quality and maturation.

The king scallop Pecten maximus is a high valuable species of great interest in Europe for both fishery and aquaculture. Notably, there has been an increased investment to produce seed for enhancement programmes of wild scallop populations. However, hatchery production is a relatively new industry and it is still underdeveloped. Major hurdles are spawning control and gamete quality. In the present study, a total of 14 scallops were sampled in the bay of Brest (Brittany, France) to compare transcriptomic profiles of mature oocytes collected by spawning induction or by stripping. To reach such a goal, a microarray analysis was performed by using a custom 8x60K oligonucleotide microarray representing 45,488 unique scallop contigs. First we identified genes that were differentially expressed depending on oocyte quality, estimated as the potential to produce D-larvae. Secondly, we investigated the transcriptional features of both stripped and spawned oocytes. Genes coding for proteins involved in cytoskeletal dynamics, serine/threonine kinases signalling pathway, mRNA processing, response to DNA damage, apoptosis and cell-cycle appeared to be of crucial importance for both oocyte maturation and developmental competence. This study allowed us to dramatically increase the knowledge about transcriptional features of oocyte quality and maturation, as well as to propose for the first time putative molecular markers to solve a major bottleneck in scallop aquaculture.

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.

Characterization of the mantle transcriptome in bivalves: Pecten maximus, Mytilus edulis and Crassostrea gigas.

The calcareous shells secreted by bivalve molluscs display diverse and species specific structural compositions, which indicates possible divergent biomineralization processes. Thus, studying multiple mollusc species will provide a more comprehensive understanding of shell formation. Here, the transcriptomes of the mantle tissues responsible for shell deposition were characterized in three commercially relevant bivalve species. Using high-throughput sequencing and bioinformatics tools, de novo transcriptome assemblies of mantle tissues were generated for the mussel Mytilus edulis, the oyster Crassostrea gigas and the scallop Pecten maximus. These transcriptomes were annotated, and contigs with similarity to proteins known to have shell formation roles in other species were identified. Comparison of the shell formation specific proteins in the three bivalves indicates the possibility of species specific shell proteins.

Mining the transcriptomes of four commercially important shellfish species for single nucleotide polymorphisms within biomineralization genes.

Transcriptional profiling not only provides insights into patterns of gene expression, but also generates sequences that can be mined for molecular markers, which in turn can be used for population genetic studies. As part of a large-scale effort to better understand how commercially important European shellfish species may respond to ocean acidification, we therefore mined the transcriptomes of four species (the Pacific oyster Crassostrea gigas, the blue mussel Mytilus edulis, the great scallop Pecten maximus and the blunt gaper Mya truncata) for single nucleotide polymorphisms (SNPs). Illumina data for C. gigas, M. edulis and P. maximus and 454 data for M. truncata were interrogated using GATK and SWAP454 respectively to identify between 8267 and 47,159 high quality SNPs per species (total=121,053 SNPs residing within 34,716 different contigs). We then annotated the transcripts containing SNPs to reveal homology to diverse genes. Finally, as oceanic pH affects the ability of organisms to incorporate calcium carbonate, we honed in on genes implicated in the biomineralization process to identify a total of 1899 SNPs in 157 genes. These provide good candidates for biomarkers with which to study patterns of selection in natural or experimental populations.

Deciphering the molecular adaptation of the king scallop (Pecten maximus) to heat stress using transcriptomics and proteomics.

BACKGROUND: The capacity of marine species to survive chronic heat stress underpins their ability to survive warming oceans as a result of climate change. In this study RNA-Seq and 2-DE proteomics were employed to decipher the molecular response of the sub-tidal bivalve Pecten maximus, to elevated temperatures. RESULTS: Individuals were maintained at three different temperatures (15, 21 and 25 °C) for 56 days, representing control conditions, maximum environmental temperature and extreme warming, with individuals sampled at seven time points. The scallops thrived at 21 °C, but suffered a reduction in condition at 25 °C. RNA-Seq analyses produced 26,064 assembled contigs, of which 531 were differentially expressed, with putative annotation assigned to 177 transcripts. The proteomic approach identified 24 differentially expressed proteins, with nine identified by mass spectrometry. Network analysis of these results indicated a pivotal role for GAPDH and AP-1 signalling pathways. Data also suggested a remodelling of the cell structure, as revealed by the differential expression of genes involved in the cytoskeleton and cell membrane and a reduction in DNA repair. They also indicated the diversion of energetic metabolism towards the mobilization of lipid energy reserves to fuel the increased metabolic rate at the higher temperature. CONCLUSIONS: This work provides preliminary insights into the response of P. maximus to chronic heat stress and provides a basis for future studies examining the tipping points and energetic trade-offs of scallop culture in warming oceans.

Emerging patterns of genetic variation in the New Zealand endemic scallop Pecten novaezelandiae.

Both historical and contemporary processes influence the genetic structure of species, but the relative roles of such processes are still difficult to access. Population genetic studies of species with recent evolutionary histories such as the New Zealand endemic scallop Pecten novaezelandiae (<1 Ma) permit testing of the effects of recent processes affecting gene flow and shaping genetic structure. In addition, studies encompassing the entire distributional range of species can provide insight into colonization processes. Analyses of genetic variation in P. novaezelandiae (952 individuals from 14 locations, genotyped at 10 microsatellite loci) revealed a weak but significant regional structure across the distributional range of the species, as well as latitudinal gradients of genetic diversity and differentiation: estimates of migration rates supported these patterns. Our results suggest that the observed genetic structure and latitudinal gradients reflect a stepping-stone model of colonization (north to south) and emerging divergence of populations as a result of ongoing limitations to gene flow and insufficient time to reach migration-drift equilibrium. The low levels of interpopulation and interregional genetic differentiation detected over hundreds of kilometres reflect the recent evolutionary history of P. novaezelandiae and stand in contrast to patterns reported for other evolutionary older species at the same spatial scale. The outcomes of this study contribute to a better understanding of evolutionary processes influencing the genetic variation of species and provide vital information on the genetic structure of P. novaezelandiae.

A novel satellite DNA isolated in Pecten jacobaeus shows high sequence similarity among molluscs.

The aim of this work is to investigate the sequence conservation and the evolution of repeated DNA in related species. Satellite DNA is a component of eukaryotic genomes and is made up of tandemly repeated sequences. These sequences are affected by high rates of mutation that lead to the occurrence of species-specific satellite DNAs, which are different in terms of both quantity and quality. In this work, a novel repetitive DNA family, named PjHhaI sat, is described in Pecten jacobaeus. The quantitative analyses revealed a different abundance of this element in the molluscan species investigated in agreement with the "library hypothesis" even if, in this case, at a high taxonomic level. In addition, the qualitative analysis demonstrated an astonishing sequence conservation not only among scallops but also in six other molluscan species belonging to three classes. These findings suggest that the PjHhaI sat may be considered as the most ancients of DNA described so far, which remained "frozen" during molluscan evolution. The widespread distribution of this sat DNA in molluscs as well as its long evolutionary preservation open up questions on the functional role of this element. A future challenge might be the identification of proteins or molecules which interact with the PjHhaI sat.

Deep sequencing of the mantle transcriptome of the great scallop Pecten maximus.

RNA-Seq transcriptome data were generated from mantle tissue of the great scallop, Pecten maximus. The consensus data were produced from a time course series of animals subjected to a 56-day thermal challenge at 3 different temperatures. A total of 26,064 contigs were assembled de novo, providing a useful resource for both the aquaculture community and researchers with an interest in mollusc shell production.

Deep transcriptome sequencing of Pecten maximus hemocytes: a genomic resource for bivalve immunology.

Pecten maximus, the king scallop, is a bivalve species with important commercial value for both fisheries and aquaculture, traditionally consumed in several European countries. Major problems in larval rearing, however, still limit hatchery-based seed production. High mortalities during early larval stages, likely related to bacterial pathogens, represent the most relevant bottleneck. To address this issue, understanding host defense mechanisms against microbes is extremely important. In this study next-generation RNA-sequencing was carried on scallop hemocytes. To enrich for immune-related transcripts, cDNA libraries from hemocytes challenged in vivo with inactivated-Vibrio anguillarum and in vitro with pathogen-associated molecular patterns, as well as unchallenged controls, were sequenced yielding 216,444,674 sequence reads. De novo assembly of the scallop hemocyte transcriptome consisted of 73,732 contigs (31% annotated). A total of 934 contigs encoded proteins with a known immune function, grouped into several functional categories. Particular attention was reserved to Toll-like receptors (TLRs), a family of pattern recognition receptors (PRRs) involved in non-self recognition. Through mining the scallop hemocyte transcriptome, at least four TLRs could be identified. The organization of canonical TLR domains demonstrated that single cysteine cluster and multiple cysteine cluster TLRs co-exist in this species. In addition, preliminary data concerning their mRNA level following bacterial challenge suggested that different members of this family could exhibit opposite responses to pathogenic stimuli. Finally, a global analysis of differential expression comparing gene-expression levels in in vitro and in vivo stimulated hemocytes against controls provided evidence on a large set of transcripts involved in the great scallop immune response.

A structural basis for substrate selectivity and stereoselectivity in octopine dehydrogenase from Pecten maximus.

Octopine dehydrogenase [N(2)-(D-1-carboxyethyl)-L-arginine:NAD(+) oxidoreductase] (OcDH) from the adductor muscle of the great scallop Pecten maximus catalyzes the reductive condensation of l-arginine and pyruvate to octopine during escape swimming. This enzyme, which is a prototype of opine dehydrogenases (OpDHs), oxidizes glycolytically born NADH to NAD(+), thus sustaining anaerobic ATP provision during short periods of strenuous muscular activity. In contrast to some other OpDHs, OcDH uses only l-arginine as the amino acid substrate. Here, we report the crystal structures of OcDH in complex with NADH and the binary complexes NADH/l-arginine and NADH/pyruvate, providing detailed information about the principles of substrate recognition, ligand binding and the reaction mechanism. OcDH binds its substrates through a combination of electrostatic forces and size selection, which guarantees that OcDH catalysis proceeds with substrate selectivity and stereoselectivity, giving rise to a second chiral center and exploiting a "molecular ruler" mechanism.

Putative reaction mechanism of heterologously expressed octopine dehydrogenase from the great scallop, Pecten maximus (L).

cDNA for octopine dehydrogenase (ODH) from the adductor muscle of the great scallop, Pecten maximus, was cloned using 5'- and 3'-RACE. The cDNA comprises an ORF of 1197 nucleotides and the deduced amino acid sequence encodes a protein of 399 amino acids. ODH was heterologously expressed in Escherichia coli with a C-terminal penta His-tag. ODH-5His was purified to homogeneity using metal-chelate affinity chromatography and Sephadex G-100 gel filtration. Recombinant ODH had kinetic properties similar to those of wild-type ODH isolated from the scallop's adductor muscle. Site-directed mutagenesis was used to elucidate the involvement of several amino acid residues for the reaction catalyzed by ODH. Cys148, which is conserved in all opine dehydrogenases known to date, was converted to serine or alanine, showing that this residue is not intrinsically important for catalysis. His212, Arg324 and Asp329, which are also conserved in all known opine dehydrogenase sequences, were subjected to site-directed mutagenesis. Modification of these residues revealed their importance for the catalytic activity of the enzyme. Conversion of each of these residues to alanine resulted in strong increases in K(m) and decreases in k(cat) values for pyruvate and L-arginine, but had little effect on the K(m) and k(cat) values for NADH. Assuming a similar structure for ODH compared with the only available structure of a bacterial opine dehydrogenase, these three amino acids may function as a catalytic triad in ODH similar to that found in lactate dehydrogenase or malate dehydrogenase. The carboxyl group of pyruvate is then stabilized by Arg324. In addition to orienting the substrate, His212 will act as an acid-base catalyst by donating a proton to the carbonyl group of pyruvate. The acidity of this histidine is further increased by the proximity of Asp329.

Karyotype and chromosomal location of 18S-28S and 5S ribosomal DNA in the scallops Pecten maximus and Mimachlamys varia (Bivalvia: Pectinidae).

This work describes the karyotype and chromosomal location of the ribosomal DNA (rDNA) of Pecten maximus and Mimachlamys varia, two commercial scallop species from Europe. According to the chromosome centromeric index values found, the karyotype of P. maximus is composed of 1 metacentric, 2 metacentric-submetacentric, 1 telocentric-subtelocentric and 15 telocentric pairs, and that of M. varia of 4 metacentric, 2 subtelocentric-submetacentric, 9 subtelocentric, 3 subtelocentric-telocentric and 1 telocentric-subtelocentric pairs. In P. maximus, 18S-28S rDNA was located by FISH on a metacentric-submetacentric pair, and in M. varia on a subtelocentric-submetacentric pair using both silver staining and FISH. PCR amplification of the 5S rDNA unit yielded a single product of about 460 bp (P. maximus) and 450 bp (M. varia), that used as probe revealed a 5S rDNA site on a telocentric pair in P. maximus and a subtelocentric pair in M. varia. Two-color FISH or sequential silver staining of 5S rDNA-FISH-metaphases corroborated that the two gene families are located on different chromosomes in both species. A comparative analysis of the data allowed the inference of karyotypic relationships within scallops.

Identification of Mytilus spp. and Pecten maximus in Irish waters by standard PCR of the 18S rDNA gene and multiplex PCR of the 16S rDNA gene.

Two molecular protocols for the identification of mussel and scallop have been developed using specific primers targeting the mitochondrial 16S ribosomal DNA gene and the nuclear 18S ribosomal DNA gene. Primers for the mitochondrial 16S ribosomal DNA gene in multiplex polymerase chain reaction (PCR) protocols yielded diagnostic DNA fragments for the mussels Mytilus edulis, Mytilus galloprovincialis, and the hybrid Mytilus edulis/galloprovincialis (335 bp), the king scallop Pecten maximus (382 bp) and the black scallop Mimachlamys varia (398 bp). DNA from the queen scallop Aequipecten opercularis showed no consistent PCR amplification of the 16S rDNA gene. Primers for the nuclear 18S rDNA gene in standard PCR protocols yielded similar-sized, diagnostic DNA fragments (approx. 190 bp) for the mussels Mytilus edulis, Mytilus galloprovincialis, and the hybrid Mytilus edulis/galloprovincialis, the king scallop Pecten maximus, the black scallop Mimachlamys varia, and the queen scallop Aequipecten opercularis. Both protocols have been tested with Mytilus spp., P. maximus, and 6 other bivalve species from a wide range of locations in Irish and European waters. Cross reaction of the specific primers with DNA template from any of the 6 other bivalve species was not observed. Rapid DNA extraction using FTA Card technology and the16S rDNA primers allowed for the detection of at least 10 mussel larvae in a subsample of natural plankton.