Comparison between single and multi-locus approaches for specimen identification in Mytilus mussels

SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

Shedding light on variation in reproductive success through studies of population genetic structure in a Southeast Pacific Coast mussel

Phylogeography often focuses on the spatial dimension of genetic diversity, rarely including the temporal dynamics occurring interannually among local populations, which can provide insight into past variations in reproductive success. Currently, there is an intense aquaculture industry of Mytilus spp. on the Southeast Pacific Coast which depends entirely on the spat released by natural populations forming a relevant and sensitive social-ecological system. Temporal and spatial spat variability from natural mussel beds could be related to interannual reproductive dynamics with variable reproductive success and recruitment, which leave genetic signatures. Temporal and spatial genetic structure was evaluated in six natural beds in the Southeast Pacific (from 39°25'S to 43°07'S) on the most abundant and widespread Mytilus lineage detected, Mytilus cf. chilensis, in 4 consecutive years. Analyses included data from >180 individuals per year, with a total of 751 (mitochondrial COI) and 747 (nuclear H1) individuals, respectively. Overall, both markers showed high haplotype diversity and low spatial and temporal genetic differentiation. Likely, the high dispersal capacity of Mytilus cf. chilensis maintains population homogeneity and prevents diversity erosion. The slight differences in genetic variance of COI were better explained by differences among sites (space), and conversely, the H1 genetic variance was better explained by interannual (temporal) comparisons, which could explain temporal variability in spat availability. This study highlights the important insights achieved with the evaluation of both temporal and spatial population genetic structures in marine species with high reproductive output, which can condition the success and sustainability of the relevant social-ecological system.

Comparison of two commercial methods for smooth-shelled mussels (Mytilus spp.) species identification

Seafood international trade has increased the labeling requirements in standards and regulations to include product information that enable traders and consumers to make informed choices. The European Union (EU) Regulation No. 1379/2013 imposes the declaration of an official commercial designation and scientific names for all the fishery and aquaculture products to be offered for sale to the final consumers. DNA analyses are used to enforce this regulation and to test authenticity in processed foods. We compared the performance of two mono-locus approaches for species identification (SI) in 61 Mytilus mussels: the high-resolution melting analysis of the polyphenolic adhesive protein gene and the partial sequencing of the histone H1C gene. The H1C sequences were analyzed with five different methods. Both approaches show discrepancies in the identification of putative hybrids (0.0 < κ < 0.687 and 0.0 < MCC < 0.724). Excluding putative hybrids, methods show substantial to perfect agreement (0.772 < κ < 1.0 and 0.783 < MCC < 1.0). This study highlights the need to use standardized molecular tools, as well as to use multi-locus methods for SI of Mytilus mussels in testing laboratories.

Benefits and limitations of a new genome-based PCR-RFLP genotyping assay (GB-RFLP): A SNP-based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.

Species identification based on a semi-diagnostic marker: Evaluation of a simple conchological test for distinguishing blue mussels Mytilus edulis L. and M. trossulus Gould

Cryptic and hybridizing species may lack diagnostic taxonomic characters leaving researchers with semi-diagnostic ones. Identification based on such characters is probabilistic, the probability of correct identification depending on the species composition in a mixed population. Here we test the possibilities of applying a semi-diagnostic conchological character for distinguishing two cryptic species of blue mussels, Mytilus edulis and M. trossulus. These ecologically, stratigraphically and economically important molluscs co-occur and hybridize in many areas of the North Atlantic and the neighboring Arctic. Any cues for distinguishing them in sympatry without genotyping would save much research effort. Recently these species have been shown to statistically differ in the White Sea, where a simple character of the shell was used to distinguish two mussel morphotypes. In this paper, we analyzed the associations between morphotypes and species-specific genotypes based on an abundant material from the waters of the Kola Peninsula (White Sea, Barents Sea) and a more limited material from Norway, the Baltic Sea, Scotland and the Gulf of Maine. The performance of the "morphotype test" for species identification was formally evaluated using approaches from evidence-based medicine. Interspecific differences in the morphotype frequencies were ubiquitous and unidirectional, but their scale varied geographically (from 75% in the White Sea to 15% in the Baltic Sea). In addition, salinity-related variation of this character within M. edulis was revealed in the Arctic Barents Sea. For every studied region, we established relationships between the proportions of the morphotypes in the populations as well as between the proportions of the morphotypes in samples and the probabilities of mussels of different morphotypes being M. trossulus and M. edulis. We provide recommendations for the application of the morphotype test to mussels from unstudied contact zones and note that they may apply equally well to other taxa identified by semi-diagnostic traits.

Development and Validation of a Multi-Locus PCR-HRM Method for Species Identification in Mytilus Genus with Food Authenticity Purposes

DNA-based methods using informative markers such as single nucleotide polymorphism (SNPs) are suitable for reliable species identification (SI) needed to enforce compliance with seafood labelling regulations (EU No.1379/2013). We developed a panel of 10 highly informative SNPs to be genotyped by PCR-High resolution melting (HRM) for SI in the Mytilus genus through in silico and in vitro stages. Its fitness for purpose and concordance were assessed by an internal validation process and by the transference to a second laboratory. The method was applicable to identify M. chilensis, M. edulis, M. galloprovincialis and M. trossulus mussels, fresh, frozen and canned with brine, oil and scallop sauce, but not in preserves containing acetic acid (wine vinegar) and tomato sauce. False-positive and negative rates were zero. Sensitivity, expressed as limit of detection (LOD), ranged between 5 and 8 ng/μL. The method was robust against small variations in DNA quality, annealing time and temperature, primer concentration, reaction volume and HRM kit. Reference materials and 220 samples were tested in an inter-laboratory assay obtaining an “almost perfect agreement” (κ = 0.925, p < 0.001). In conclusion, the method was suitable for the intended use and to be applied in the seafood industry.

First description of a widespread Mytilus trossulus-derived bivalve transmissible cancer lineage in M. trossulus itself

Two lineages of bivalve transmissible neoplasia (BTN), BTN1 and BTN2, are known in blue mussels Mytilus. Both lineages derive from the Pacific mussel M. trossulus and are identified primarily by their unique genotypes of the nuclear gene EF1α. BTN1 is found in populations of M. trossulus from the Northeast Pacific, while BTN2 has been detected in populations of other Mytilus species worldwide but not in M. trossulus itself. Here we examined M. trossulus from the Sea of Japan (Northwest Pacific) for the presence of BTN. Using hemocytology and flow cytometry of the hemolymph, we confirmed the presence of disseminated neoplasia in our specimens. Cancerous mussels possessed the BTN2 EF1α genotype and two mitochondrial haplotypes with different recombinant control regions, similar to that of common BTN2 lineages. This is the first report of BTN2 in its original host species M. trossulus. A comparison of all available BTN and M. trossulus COI sequences suggests a common and recent origin of BTN2 diversity in populations of M. trossulus outside the Northeast Pacific, possibly in the Northwest Pacific.