Next Article in Journal
Telomere Length Changes during Critical Illness: A Prospective, Observational Study
Next Article in Special Issue
Population Genomics of an Anadromous Hilsa Shad Tenualosa ilisha Species across Its Diverse Migratory Habitats: Discrimination by Fine-Scale Local Adaptation
Previous Article in Journal
A Portal to Visualize Transcriptome Profiles in Mouse Models of Neurological Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Useful SNP Panel to Distinguish Two Cockle Species, Cerastoderma edule and C. glaucum, Co-Occurring in Some European Beds, and Their Putative Hybrids

1
Department of Zoology, Genetics and Physical Anthropology, ACUIGEN Group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002 Lugo, Spain
2
Centro de Investigacións Mariñas (CIMA), Consellería do Mar, Xunta de Galicia, 36620 Vilanova de Arousa, Spain
3
Department of Zoology, Genetics and Physical Anthropology, Mobile Genomes and Disease Group, CIMUS, Universidade de Santiago de Compostela, Campus of Santiago de Compostela, 15706 Santiago de Compostela, Spain
4
Department of Life Sciences, University of Alcalá, 28871 Alcalá de Henares, Spain
5
Research Centre for Experimental Marine Biology and Biotechnology (PIE), University of the Basque Country (UPV/EHU), 48620 Plentzia, Basque Country, Spain
6
Instituto de Acuicultura, Universidade de Santiago de Compostela, 15705 Santiago de Compostela, Spain
*
Author to whom correspondence should be addressed.
Present Address: Dipartimento di Scienze della Vita e Biotecnologia (SVeB), Università degli Studi di Ferrara, via Luigi Borsari, 46–44121 Ferrara, Italy.
Submission received: 29 August 2019 / Revised: 25 September 2019 / Accepted: 26 September 2019 / Published: 27 September 2019
(This article belongs to the Special Issue Genetic Diversity of Marine Populations)

Abstract

:
Cockles are highly appreciated mollusks and provide important services in coastal areas. The two European species, edible (Cerastoderma edule) and lagoon (Cerastoderma glaucum) cockles, are not easily distinguishable, especially when young. Interestingly, the species show different resistance to Marteilia cochillia, the parasite responsible for marteiliosis outbreaks, which is devastating cockle production in some areas. C. edule is severely affected by the parasite, while C. glaucum seems to be resistant, although underlying reasons are still unknown. Hybrids between both species might be interesting to introgress allelic variants responsible for tolerance, either naturally or through artificial selection, from lagoon into edible cockle. Here, we used 2b restriction site-associated DNA sequencing (2b–RAD) to identify single nucleotide polymorphisms (SNP) diagnostic for cockle discrimination (fixed for alternative allelic variants). Among the nine diagnostic SNPs selected, seven were validated using a SNaPshot assay in samples covering most of the distribution range of both species. The validated SNPs were used to check cockles that were suggested to be hybrids by a claimed diagnostic tool based on the internal transcribed spacers of the ribosomal RNA. Although these were shown to be false positives, we cannot rule out the fact that hybrids can occur and be viable. The SNP tool here developed will be valuable for their identification and management.

1. Introduction

True cockles are members of the genus Cerastoderma, a group of bivalve mollusks widely distributed in the Mediterranean and Atlantic coasts of Europe [1]. Production of cockles, especially the species most valued as food, edible cockle (Cerastoderma edule (Linnaeus, 1758)), is especially relevant in Spain, the British Isles, and France, with production reaching up to 24,626 tons in 2017 in Europe [2]. European production of cockles is mostly exported to Galicia (NW Spain), where the species is marketed as fresh or canned food. In addition to food provisioning, cockles provide many ecosystem services [3], which are especially important for the development and support of many coastal areas. Edible cockle production has been severely affected by marteiliosis outbreaks caused by Martelia cochillia, especially in Galicia (NW Spain), where the parasite has devastated its production in the Southern estuaries (Rias). This parasite was first recorded in that area in 2012 [4], but its first description was reported in Catalonia (NE Spain) from cockles identified as C. edule [5], despite the typical cockle species occurring in the Mediterranean coasts is the lagoon cockle (Cerastoderma glaucum (Poiret, 1789)) [1]. Interestingly, in the Galician coast, only C. edule is affected by martelliosis, while C. glaucum seems to be resistant to this parasite [6]. The known geographic ranges of both Cerastoderma species [1] supports the fact that cockles occurring in the Mediterranean Sea are C. glaucum, and accordingly, the cockles infected with M. cochillia in Catalonia should be C. glaucum. This would suggest different susceptibility to this parasite between the Mediterranean (Catalonian) and the Atlantic (Galician) types of C. glaucum [6], thus adding new data to the convenience of splitting C. glaucum into different species or subspecies [7,8]. This information would refine the previous report suggesting that Catalonian and Galician populations would be included in the same lineage of the three identified for C. glaucum using restriction site-associated DNA Sequencing (RADseq) [9], which is widely distributed from the Western Mediterranean to the Baltic Sea [9]. Alternatively, the hypothesis of different virulence between the parasites M. cochillia occurring in Catalonia and Galicia should not be discarded. The reasons behind this differential resistance between C. edule and C. glaucum in Galicia are still unknown, and deserve special attention within the COCKLES EU project following a genomic approach in the most affected area, the Ria of Arousa (42°30′ N–08°56′ W, EAPA_458/2016). Furthermore, the tolerance of C. glaucum to M. cochillia makes putative hybrids especially interesting. Checking its viability and reproduction would be an interesting way to introgress allelic variants responsible for tolerance either naturally or through artificial selection in C. edule.
The two recognized cockle species, C. edule and C. glaucum, are not easy distinguishable at a glance, especially when they are young. However, there are some morphological characteristics that trained eyes can use for discrimination [10,11]. Hybrids have been claimed to be generated in captivity [12], but limited evidence that they were real hybrids was provided. Despite the two species showing differences in habitat preferences and lifestyle (see Tarnowska et al. [13] for a detailed description), they present overlapping geographical distributions in various sites of the European Atlantic coast [10,11], and overlapping spawning periods have been supported for these species [14]. Although no records of hybrids have been reported to date, their presence has been recently suggested (Seila Díaz, unpublished data) using a reported diagnostic tool based on the rDNA internal transcribed spacer (ITS) [15]. ITSs have been used to discriminate species quite frequently among parasites [16], but also for aquaculture species such as cockles [17] and mud crabs (Scilla sp.) [18]. ITSs represent a very useful and straightforward way to develop molecular tools for species discrimination because of the existence of conserved regions for universal primer design in the adjacent 18S, 5.8S, and 28S genes. However, the tandem constitution of rDNA genes and the multiple location and site polymorphism reported for Nucleolar Organizer Regions (NOR) reported in several species, including cockles [19], suggest caution and strict validation using consistent reference tools.
The development of cutting-edge genomic technologies and protocols in the last decade has opened the possibility to develop a wide range of genotyping by sequencing methods that enable simultaneously identifying and genotyping thousands of single nucleotide polymorphisms (SNP) at low cost [20]. This approach has also been applied to compare large portions of the genomes of individuals from closely related species to look for diagnostic markers that make possible suitable discrimination between taxa [21]. Furthermore, these SNP-based tools enable distinguishing F1, F2, and backcrosses with high confidence, thus allowing a refined evaluation of the hybridization degree of individuals and populations. Here, we describe the identification of a set of diagnostic SNPs for a reliable identification of C. edule, C. glaucum, and their putative hybrids in natural beds of the Atlantic Area, where geographical distribution of both species overlaps. Nine markers fixed for alternative allelic variants identified in “pure” individuals of the two species were validated in samples covering a wide distribution range and tuned up in a SNaPshot assay designed to provide a cheap and straightforward diagnostic tool. This SNP panel will be useful for management of extensive and intensive production areas where both species co-occur, and to address more sophisticated tools for studying the resistance of cockles to marteiliosis.

2. Materials and Methods

2.1. Samples for Genomic 2b Restriction Site-Associated DNA Sequencing (2b–RAD) Libraries

Thirty individuals of C. glaucum and 120 individuals of C. edule collected over 2017–2018 were genotyped using the 2b–RAD genotyping-by-sequencing (GBS) method [22] to identify diagnostic SNPs to distinguish the two species and their hybrids. Edible cockles were collected in four Atlantic European locations: Somme Bay, Miño, Campelo, and Ria Formosa (from north to south, 30 individuals collected in each sampling site) covering a wide distribution range (Table 1; Figure 1), while lagoon cockles were collected in Redondela (N = 30) (Figure 1).

2.2. DNA Extraction and Single Nucleotide Polymorphisms Calling and Genotyping

DNA was extracted from gill tissue using the E.Z.N.A. E-96 mollusk DNA kit (OMEGA Bio-tech, Norcross, GA, USA), following the manufacturer recommendations. SNP identification and selection, as well as genotyping and validation protocols, were similar to those described by Maroso et al. [21]. Briefly, genomic DNA was cut using the AlfI IIb restriction enzyme (RE), and then 2b–RAD libraries were constructed by joining adaptors to both fragment ends, followed by PCR amplification using primers targeting specific regions within adaptors. After PCR amplification, DNA from barcoded individuals was pooled for sequencing in an Illumina Next-seq machine. Ninety individuals were routinely multiplexed per run. After demultiplexing and filtering raw reads by quality (Phred value > 30) and the presence of the AlfI restriction site, Stacks 2.0 [23] was used to align all reads from both species to identify and characterize shared genomic loci (RAD–tags) between both species; this information was used to create a consensus catalog of loci and to call SNPs using the whole population data from both species using a de novo approach. A minimum of three identical reads was required, and up to three mismatches were allowed between reads to be considered as part of the same locus. The obtained SNP dataset was parsed to identify diagnostic SNPs (those with different alternative alleles fixed in the two analyzed species). This list was additionally filtered, and RAD–tags containing only one SNP and called in at least 15 C. glaucum and 30 C. edule specimens were retained. RAD–tags of this filtered dataset were then aligned against the C. edule draft genome (Tubío et al, unpublished) and only those aligning to a single genome site were retained. Sequences of the C. edule genome 150 bp up- and down-stream from the diagnostic SNP position were used to design external primers using the software Primer3 [24].

2.3. Validation of Diagnostic Single Nucleotide Polymorphisms

Among the diagnostic loci retained, a small subset of 20 was selected for the final validation (see Results and Discussion). For this, external primers were checked in three C. glaucum and three C. edule specimens. Loci were amplified in 15 μL volume, including 1 μL template DNA (~30 ng), 1X PCR Gold Buffer (Applied Biosystems, Forest City, CA, USA), 1.5 mM of MgCl2, 100 µM of dNTP, 10 pmol of both forward and reverse PCR primers, and 0.5 U of Amplitaq GoldTM DNA polymerase (Applied Biosystems). Thermal cycling was conducted on a VerityTM 96-Well Thermal Cycler (Applied Biosystems) as follows: initial denaturation at 95 °C for 10 min, 35 cycles of denaturation at 94 °C for 45 s, annealing at 60 °C for 50 s, and extension at 72 °C for 50 s. There was a final extension step at 72 °C for 10 min. Only primer pairs showing a single neat band in agarose gels and the expected amplicon size in both species were selected for further steps. Selected RAD–tag amplicons were sequenced in three C. glaucum and three C. edule specimens following the ABI Prism BigDye™ Terminator v3.1 Cycle Sequencing Kit protocol in an ABI 3730xl DNA Genetic Analyzer (Applied Biosystems) in order to confirm the consensus sequences around the diagnostic SNPs, from which internal primers could be designed for genotyping. Annotation of these sequences was carried out using BLASTn with default parameters and e-value <1e–5 within the NCBI database. Selected SNPs were genotyped with the SNaPshot methodology (Applied Biosystems), which is based on the fluorescent detection of the SNP variants through a two-step reaction protocol. The first step involves the amplification by PCR of a region including the selected SNP (using the external primers), while the second step is a 1-base sequencing reaction from the adjacent SNP primer (i.e., internal primer). For SNP multiplexing, internal primers were designed with length differences by adding CAGT tails to the 5’ end of the primer sequence (see Table 2). Primer compatibility for multiplexing was checked using ThermoFisher web application [25] and one multiplex reaction following the methodology described by Vera et al. [26] was finally designed to analyze the nine selected SNPs (Table 2). For genotyping, 0.2 µM of each primer was used, both for PCR amplification with the external primers and for the SNaPshot reaction using internal primers. The annealing temperatures for PCR amplification and the SNaPshot reaction were 60 °C and 55 °C, respectively. The obtained SNaPshot assay was finally tested on 48 individuals morphologically identified of each species (total = 96), covering their main distribution ranges and including specimens previously genotyped by 2b–RAD GBS methodology (see Table 1). Finally, five different beds from Denmark (1) and Galicia (4) (total individuals = 46; see Table 1), where both species co-occur and suspected hybrids had been identified (Díaz et al. unpublished data) using the methodology described by Freire et al. [15], were analyzed with the developed SNP tool.

3. Results and Discussion

All specimens of both species used for identifying diagnostic SNPs were retained for further analysis, except one C. glaucum sample, which was eliminated before filtering due to the low number of raw reads. On average, 2,768,161 reads per specimen were retained after quality filtering. Stacks analysis identified a total of 13,715 SNPs in RAD–tags, including reads from both species, 2057 diagnostic and 11,658 restricted to one or other of the two species (Tables S1 and S2, respectively). Among the latter, 1658 were private of C. edule (1190 with a minimum allele frequency (MAF) ≥ 0.05; 71.8%) and 10,000 of C. glaucum (3028 with MAF ≥ 0.05; 30.3%). These figures, despite being representative of a particular fraction of the genome (RAD–tags including diagnostic SNPs), strongly suggest a much higher genetic diversity in C. glaucum than in C. edule, especially regarding the number of SNPs per physical length; furthermore, much of the difference was due to the higher proportion of rare alleles in C. glaucum than in C. edule (MAF < 0.05: 69.7% vs. 28.2%, respectively), suggestive of higher effective population size and/or no recent bottlenecks. Moreover, the most representative C. glaucum beds were located within the Southern glacial refugia of the species (i.e., Iberian Peninsula), where higher genetic diversity is expected [13]. As usually reported in mollusk species [27,28], a global significant heterozygote deficit was observed in both species (FIS > 0; P = 0), but was significantly higher in C. glaucum than in C. edule (FIS 0.252 vs. 0.129).
The 2057 diagnostic loci were identified in 1534 different RAD–tags, meaning that more than one diagnostic SNP was detected in many RAD–tags; nearly 50% contained two diagnostic loci and six RAD–tags up to three diagnostic SNPs (Table S1). Only 157 RAD–tags contained a single diagnostic SNP and no other polymorphism (i.e., private SNPs). Since these 157 RAD–tags showed no genetic variation within species, they were considered the most suitable ones for designing primers, taking the C. edule genome as a reference. After eliminating RAD–tags shared by a low number of individuals of the two species, 60 were retained, and among these, 46 matched to unique positions in the C. edule genome (Tubío et al., unpublished data). A final selection of 20 RAD–tags was randomly made for the final validation step using the SNaPshot protocol, but the remaining 26 RAD–tags represent a useful repository for increasing the number of diagnostic SNPs if necessary.
Twenty primer pairs were designed with Primer 3 using the C. edule genome as a reference and their performance for PCR amplification was checked on agarose gels in three C. edule individuals from Ría de Noia (SNO) and three C. glaucum from Redondela (SRG). Nine loci showed a single band of expected size, thus being considered suitable for the last step of validation using the SNaPshot methodology (Cerast_31, Cerast_450, Cerast_456, Cerast_586, Cerast_1173A, Cerast_1255, Cerast_1316, Cerast_1400, and Cerast_2530, Table 2). The sequences of these nine RAD–tags confirmed that those markers showed a single diagnostic difference between both species (GenBank Accession Numbers: MN178488–MN178496). Four out of the nine sequences matched against annotated genes, but these markers were not associated with either mtDNA or ribosomal genes (Table S1). Internal primers were then designed adjacent to the SNP position for these nine loci to achieve the best compatibility for multiplexing the SNaPshot mini-sequencing reaction (Table 2). The obtained genotypes in the individuals analyzed with both techniques (i.e., 2b–RAD GBS and SNaPshot) were identical (Table S3), confirming the correspondence between them and giving additional support to the designed molecular tool. When individuals from distant beds not included in the initial RADseq screening were genotyped, diagnostic differences remained for all loci except for Cerast_586 and Cerast_1173A, which showed a low polymorphism for C. glaucum in Mediterranean samples (Table S3). Previous molecular tools have been reported for the identification of cockle species, but they were focused on maternally inherited mtDNA [29], on dominant random amplified polymorphic DNA (RAPDs) patterns [30] and on the ITSs of the tandemly arranged rDNA gene family [15,31], all of them displaying different limitations for the identification of hybrids. The seven SNPs here reported as diagnostic and developed by Next Generation Sequencing (NGS) and Genotyping-By-Sequencing (GBS) methodologies allow for a more refined analysis of hybridization between both species if it exists, detecting with fully confidence F1 hybrids and with reasonable accuracy their backcrosses (BC) and F2 hybrids (P > 0.95).
The 10 suspected hybrids detected (Díaz et al., unpublished data) using the methodology ITS region [15] were analyzed with the tool here described, and all of them showed pure multilocus genotypes either for C. edule (three individuals) or C. glaucum (seven individuals) (Table S4), strongly suggesting that they are not hybrids. In fact, the most intense ITS band shown on agarose gels by these individuals corresponded with the assigned species with the SNP tool (data not shown). Although we cannot fully discard a recurrent backcross introgression from one species into another occurring several generations ago, a more plausible explanation for this observation is the existence of an intraspecific rDNA polymorphism in cockles, as previously reported in C. glaucum [19] and in many other species [32,33,34,35]. As outlined above, the probability of being a BC or F2 hybrid with pure species genotype would be < 0.05 for seven diagnostic loci, and this would be much lower for the multiple individuals detected showing that pattern. Despite the fact that this observation does not fully invalidate the ITS tool by Freire et al. [15], methodological modifications (e.g., changes in annealing temperature) should be addressed to avoid the appearance of unspecific bands in the gels.
In summary, our data support the usefulness of the seven diagnostic SNPs reported from a large sample of individuals covering the main distribution range for discriminating edible and lagoon cockles; furthermore, a larger SNP diagnostic repository is available in case more markers were needed depending on the goals of the study. In addition, our study highlights the usefulness of NGS strategies to identify diagnostic differences in the genomes of closely related species for the analysis of hybridization in the wild, useful for management of genetic resources in productive areas or in those areas where cockles provide ecosystem services.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/10/10/760/s1, Table S1: Diagnostic single nucleotide polymorphisms (SNPs) for both cockle species analyzed in the present study. SNP position in the RADtag sequence (SNP position), SNP variants (Allele 1 and Allele 2), number of analyzed individuals (N), marker name (see main text and Table 2) and annotation against NCBI database are shown, Table S2: Private alleles detected for each cockle species. SNPs variants (Allele 1 and Allele 2), average number of reads per marker (Av. Nb reads), number of individuals analyzed (N), observed and expected heterozygosities (Ho and He, respectively), and inbreeding coefficients (Fis) are shown, Table S3: Genotypes obtained for the 96 individuals (48 C. edule and 48 C. glaucum) tested with the SNaPshot technique. Expected genotypes for each species with 2b–RAD–GBS technique and those previously found with this technique in control individuals are also shown. Morphological and molecular species identification is included for each individual. Bed codes are shown on Table 1, Table S4: Genotypes obtained with the SNaPshot technique for the 46 individuals from beds where the presence of hybrids has been suggested. For each individual, their identification using the methodology described by Freire et al. (2011) using ITS region (ITS) and the SNaPshot methodology described in the present study (SNaPshot) are shown. The 10 suspected hybrids and the multilocus genotypes for each species are shown in yellow and green, respectively. Bed codes are shown on Table 1.

Author Contributions

Conceptualization, M.V. and P.M.; methodology, F.M. and C.P.d.G.; software, F.M. and C.P.d.G.; validation, F.M., C.P.d.G. and M.V.; formal analysis, F.M.; investigation, F.M.; resources, D.I., A.C., S.D., A.V., P.M.; data curation, F.M., C.P.d.G. and M.V.; writing—original draft preparation, F.M. and M.V.; writing—review and editing, C.P.d.G., D.I., A.C., S.D., A.V. and P.M.; visualization, M.V. and P.M.; supervision, M.V. and P.M.; project administration, M.V. and P.M.; funding acquisition, P.M.

Funding

This research funded by INTERREG EUROPEAN PROGRAMME, grant number EAPA_458/2016 (COCKLES project).

Acknowledgments

Authors wish to thank L. Insua, S. Sánchez-Darriba and S. Gómez for their technical support. Authors are also indebted with COCKLES, VIVALDI and Scuba Cancers partners, together with M.L. Conde-Varela, M.T. Fernández-Núñez, M. García-Graña, P. Luttikhuizen, S. Pereira, A. Simón, and L. Solís for providing many of the samples analyzed.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Malham, S.K.; Hutchinson, T.H.; Longshaw, M. A review of the biology of European cockles (Cerastoderma spp.). J. Mar. Biol. Assoc. UK 2012, 92, 1563–1577. [Google Scholar] [CrossRef]
  2. Food and Agriculture Organization of the United Nations (FAO). Statistics. Available online: http://www.fao.org/fishery/statistics/global-production/en (accessed on 20 August 2019).
  3. Donadi, S.; van der Heide, T.; van der Zee, E.M.; Eklöf, J.S.; de Koppel, J.V.; Weerman, E.J.; Piersma, T.; Olff, H.; Eriksson, B.K. Cross-habitat interactions among bivalve species control community structure on intertidal flats. Ecology 2013, 94, 489–498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Villalba, A.; Iglesias, D.; Ramilo, A.; Darriba, S.; Parada, J.M.; No, E.; Abollo, E.; Molares, J.; Carballal, M.J. Cockle Cerastoderma edule fishery collapse in the Ría de Arousa (Galicia, NW Spain) associated with the protistan parasite Marteilia cochillia. Dis. Aquat. Org. 2014, 109, 55–80. [Google Scholar] [CrossRef] [PubMed]
  5. Carrasco, N.; Hine, P.M.; Durfort, M.; Andree, K.B.; Malchus, N.; Lacuesta, B.; González, M.; Roque, A.; Rodgers, C.; Furones, M.D. Marteilia cochillia sp. nov., a new Marteilia species affecting the edible cockle Cerastoderma edule in European Waters. Aquaculture 2013, 412–413, 223–230. [Google Scholar] [CrossRef]
  6. Carballal, M.J.; Iglesias, D.; Darriba, S.; Cao, A.; Mariño, J.C.; Ramilo, A.; No, E.; Villalba, A. Parasites, pathological conditions and resistance to Marteilia cochillia in lagoon cockle Cerastoderma glaucum from Galicia (NW Spain). Dis. Aquat. Org. 2016, 122, 137–152. [Google Scholar] [CrossRef] [PubMed]
  7. Brock, V. Genetic relations between the bivalves Cardium (Cerastoderma) edule, Cardium lamarcki and Cardium glaucum, studied by means of crossed immunoelectrophoresis. Mar. Biol. 1987, 93, 493–498. [Google Scholar] [CrossRef]
  8. Ladhar-Chaabouni, R.; Hamza-Chaffai, A.; Hardivillier, Y.; Chénais, B.; Denis, F. A pilot study of genetic differentiation between two phenotypes of a Mediterranean population of the bivalve Cerastoderma glaucum and genetic discrimination with other Cerastoderma glaucum and Cerastoderma edule populations outside the Mediterranean. Mar. Ecol. 2010, 31, 355–363. [Google Scholar] [CrossRef]
  9. Sromek, L.; Forcioli, D.; Lasota, R.; Furla, P.; Wolowicz, M. Next-generation phylogeography of the cockle Cerastoderma glaucum: Highly heterogeneous genetic differentiation in a lagoon species. Ecol. Evol. 2019, 9, 4667–4682. [Google Scholar] [CrossRef]
  10. Brock, V. Morphological and biochemical criteria for the separation of Cardium glaucum (Bruguière) from Cardium edule (L.). Ophelia 1978, 17, 207–214. [Google Scholar] [CrossRef]
  11. Parada, J.M. Validation of lateral visibility of the ligament as a characteristic for fast discrimination between juveniles of Cerastoderma edule and C. glaucum (Mollusca, Bivalvia). Nova Acta Cient. Compost. 2018, 25, 1–8. [Google Scholar]
  12. Kingston, P.F. Interspecific hybridization in Cardium. Nature 1973, 243, 360. [Google Scholar] [CrossRef]
  13. Tarnowska, K.; Krakau, M.; Jacobsen, S.; Wolowick, M.; Féral, J.P.; Chenuil, A. Comparative phylogeography of two sister (congeneric) species of cardiid bivalve: Strong influence of habitat, life history and post-glacial history. Estuar. Coast. Shelf Sci. 2012, 107, 150–158. [Google Scholar] [CrossRef]
  14. Yankson, K. Reproductive cycles of Cerastoderma glaucum (Bruguiére) and C. edule (L.) with special reference to the effects of the 1981-82 severe winter. J. Molluscan Stud. 1986, 52, 6–14. [Google Scholar] [CrossRef]
  15. Freire, R.; Arias, A.; Méndez, J.; Insua, A. Identification of European comercial cockles (Cerastoderma edule and C. glaucum) by species-specific PCR amplification of the ribosomal DNA ITS region. Eur. Food Res. Technol. 2011, 232, 83–86. [Google Scholar] [CrossRef]
  16. Toz, S.O.; Tetik, A.; Ertabaklar, H.; Zeyrek, F.Y.; Culha, G.; Ozbel, Y. Optimisation of diagnostic real time PCR assays for ITS-1 and ITS-2 regions of Leishmania Genus. Trop. Med. Int. Health 2009, 14, 181–182. [Google Scholar]
  17. Freire, R.; Arias, A.; Mendez, J.; Insua, A. Sequence variation of the internal transcribed spacer (ITS) region of ribosomal DNA in Cerastoderma species (Bivalvia: Cardiidae). J. Molluscan Stud. 2010, 76, 77–86. [Google Scholar] [CrossRef]
  18. Vince Cruz-Abeledo, C.C.; Ting, K.; Ablan-Lagman, M.C. Identification of a species diagnostic character for instar and juvenile mud crabs (Genus Scylla). Aquaculture 2018, 491, 205–209. [Google Scholar] [CrossRef]
  19. Thiriot-Quievreux, C.; Wolowicz, M. Karyotypes of Cerastoderma glucum (Bivalvia) from Baltic and Mediterranean populations. Hydrobiologia 1996, 324, 149–155. [Google Scholar] [CrossRef]
  20. Robledo, D.; Palaiokostas, C.; Bargelloni, L.; Martinez, P.; Houston, R. Applications of genotyping by sequencing in aquaculture breeding and genetics. Rev. Aquacult. 2018, 10, 670–682. [Google Scholar] [CrossRef]
  21. Maroso, F.; Casanova, A.; do Prado, F.D.; Bouza, C.; Pardo, B.G.; Blanco, A.; Hermida, M.; Fernandez, C.; Vera, M.; Martinez, P. Species identification of two closely exploited flatfish, turbot (Scophthalmus maximus) and brill (Scophthalmus rhombus), using a ddRADseq genomic approach. Aquat. Conserv. Mar. Freshw. Ecosys. 2018, 28, 1253–1260. [Google Scholar] [CrossRef]
  22. Wang, S.; Meyer, E.; McKay, J.K.; Matz, M.V. 2b-RAD: A simple and flexible method for genome-wide genotyping. Nat. Methods 2012, 9, 808–810. [Google Scholar] [CrossRef] [PubMed]
  23. Catchen, J.; Hohenlohe, P.A.; Bassham, S.; Amores, A.; Cresko, W.A. Stacks: An analysis tool set for population genomics. Mol. Ecol. 2013, 22, 3124–3140. [Google Scholar] [CrossRef] [PubMed]
  24. Rozen, S.; Skaletsky, H. Primer3 on the WWW for general users and for biologist programmers. In Bioinformatics Methods and Protocols: Methods in Molecular Biology; Krawetz, S., Misener, S.E., Eds.; Humana Press: Totowa, NJ, USA, 2000; pp. 365–386. [Google Scholar]
  25. ThermoFisher. Multiple Primer Analyzer. Available online: https://www.thermofisher.com/es/es/home/brands/thermo-scientific/molecular-biology/molecular-biology-learning-center/molecular-biology-resource-library/thermo-scientific-web-tools/multiple-primer-analyzer.html (accessed on 3 June 2019).
  26. Vera, M.; Álvarez-Dios, J.A.; Millán, A.; Pardo, B.G.; Bouza, C.; Hermida, M.; Fernandez, C.; de la Herran, R.; Molina-Luzon, M.J.; Martínez, P. Validation of single nucleotide polymorphism (SNP) markers from an immune Expressed Sequence Tag (EST) turbot, Scophthalmus maximus, database. Aquaculture 2011, 313, 31–41. [Google Scholar] [CrossRef]
  27. Gardeström, J.; Pereyra, R.T.; Andre, C. Characterization of six microsatellite loci in the Baltic blue mussel Mytilus trossulus and cross-species amplification in North Sea Mytilus edulis. Conserv. Genet. 2008, 9, 1003–1005. [Google Scholar] [CrossRef]
  28. Diz, A.P.; Presa, P. The genetic diversity pattern of Mytilus galloprovincialis in Galician Rias (NW Iberian estuaries). Aquaculture 2009, 287, 278–285. [Google Scholar] [CrossRef]
  29. Brock, V.; Christiansen, G. Evolution of Cardium (Cerastoderma) edule, C. lamarcki and C. glaucum: Studies of DNA-variation. Mar. Biol. 1989, 102, 505–511. [Google Scholar] [CrossRef]
  30. André, C.; Lindegarth, M.; Jonsson, P.R.; Sundberg, P. Species identification of bivalve larvae using random amplified polymorphic DNA (RAPD): Differentiation berween Cerastoderma edule and C. lamarcki. J. Mar. Biol. Assoc. UK 1999, 79, 563–565. [Google Scholar] [CrossRef]
  31. Freire, R.; Insua, A.; Mendez, J. Cerastoderma glaucum 5S ribosomal DNA: Characterization of the repeat unit, divergence with respect to Cerastoderma edule, and PCR-RFLPs for the identification of both cockles. Genome 2005, 48, 427–442. [Google Scholar] [CrossRef] [PubMed]
  32. Castro, J.; Rodriguez, S.; Pardo, B.G.; Sanchez, L.; Martinez, P. Population analysis of an unusual NOR-site polymorphism in brown trout (Salmo trutta L.). Heredity 2001, 86, 291–302. [Google Scholar] [CrossRef]
  33. Singh, M.; Kumar, R.; Nagpure, N.S.; Kushwaha, B.; Mani, I.; Lakra, W.S. Extensive NOR site polymorphism in geographically isolated populations of Golden mahseer, Tor putitora. Genome 2009, 52, 783–789. [Google Scholar] [CrossRef]
  34. Alquezar, D.E.; Hemmerter, S.; Cooper, R.D.; Beebe, N.W. Incomplete concerted evolution and reproductive isolation at the rDNA locus uncovers nine cryptic species within Anopheles longirostris from Papua New Guinea. BMC Evol. Biol. 2010, 10, 392. [Google Scholar] [CrossRef] [PubMed]
  35. Lunerova, J.; Renny-Byfield, S.; Matyasek, R.; Leitch, A.; Kovarik, A. Concerted evolution rapidly eliminates sequence variation in rDNA coding regions but not in intergenic spacers in Nicotiana tabacum allotetraploid. Plant Syst. Evol. 2017, 303, 1043–1060. [Google Scholar] [CrossRef]
Figure 1. Geographical situation of the natural beds analyzed in the present studies. Codes are shown on Table 1.
Figure 1. Geographical situation of the natural beds analyzed in the present studies. Codes are shown on Table 1.
Genes 10 00760 g001
Table 1. Cockle beds used in the present study. Individuals used for genomic libraries for 2b restriction site-associated DNA sequencing genotyping by sequencing (2b–RAD GBS) and those analyzed with the SNaPshot molecular tool (SNaPshot) developed in the present study are shown.
Table 1. Cockle beds used in the present study. Individuals used for genomic libraries for 2b restriction site-associated DNA sequencing genotyping by sequencing (2b–RAD GBS) and those analyzed with the SNaPshot molecular tool (SNaPshot) developed in the present study are shown.
Bed (Drainage)CodeGeographical CoordinatesNSpecimensTechnique
Cerastoderma eduleCerastoderma glaucumSuspected Hybrids2bRAD GBSSNaPshot
Texel—Netherlands (Atlantic)NTX53°00.22′ N,
4°46.23′ E
88 8
Somme Bay—France (Atlantic)FBS50°12.08′ N,
1°37.62′ W
3030 308
Burry—Wales (Atlantic)WBY51°38.55′ N,
4°09.98′ W
88 8
Dundalk Bay—Ireland (Atlantic)IDA53°53.05′ N,
6°20.48′ W
88 8
Miño–Galicia Spain (Atlantic)SMI43°21.69′ N,
8°12.35′ W
3030 30
Ria de Noia—Spain (Atlantic)SNO42°47.42′ N,
8°55.36′ W
88 8
Campelo—Spain (Atlantic)SCA42°25.25′ N,
8°41.09′ W
3030 30
Redondela—Spain (Atlantic)SRG42°17.63′ N,
8°37.24′ W
30 30 3019
Tejo Mouth–Portugal (Atlantic)PTE38°46.00′ N,
9°02.00′ W
88 8
Ria Formosa–Portugal (Atlantic)PRF36°59.86′ N,
7°49.81′ W
3030 30
Delta del Ebro—Spain (Mediterranean)EDCE40°41.75′ N,
0°45.12′ E
29 29 29
Nykobing Mors—Denmark (Atlantic) *DNCE56°46.91′ N,
8°52.27′ W
11623 11
Espasante—Spain (Atlantic) *EECE43°43.44′ N,
7°48.26′ W
98 1 9
Carril—Spain (Atlantic) *ECCE42°36.75′ N,
8°46.50′ W
10631 10
Combarro—Spain (Atlantic) *EACE42°25.92′ N,
8°42.30′ W
74 3 7
Lourizán—Spain (Atlantic) *ELCE42°24.53′ N,
8°40.66′ W
97 2 9
* Locations where suspected hybrids have been detected.
Table 2. External and internal primers (IP) for the two multiplex reactions developed in the present study. Lowercase letters in internal primers indicate the GACT nucleotide tails. Single nucleotide polymorphisms (SNPs) variants (with their position within their submitted GenBank sequence between parentheses), sense of the internal primer (IP sense; F: Forward, R: Reverse) and GenBank accession numbers for the sequences used to design the markers (length in bp of the sequence between parentheses) are shown.
Table 2. External and internal primers (IP) for the two multiplex reactions developed in the present study. Lowercase letters in internal primers indicate the GACT nucleotide tails. Single nucleotide polymorphisms (SNPs) variants (with their position within their submitted GenBank sequence between parentheses), sense of the internal primer (IP sense; F: Forward, R: Reverse) and GenBank accession numbers for the sequences used to design the markers (length in bp of the sequence between parentheses) are shown.
MarkerSNPs Variants C.edule/C. glaucumExternal Primers (5′-3′)Tailed Internal Primer (5′-3′)IP SenseGenBank Accession Number
Cerast_1173AA/G
(129)
F: GGGACGGCACTTTTCACAAT
R: TGGTGCTGTGGATGAATCGA
gactTGGTGGGCACTTGGATGCRMN178492 (240)
Cerast_1316T/C
(55)
F: TAGACAAAACAGGCCTACGC
R: TCGTGATCTGCCAAAGGTTT
gactgactgacGCAACAGATTGCCAGCTGTRMN178494 (171)
Cerast_1400G/A
(99)
F: AGCACGGTTGTTGATTGGAC
R: TTTAAGCCAGGGTCCTCGG
gactgactgactACTCTTCTTCATGGTTGAAAAGTCFMN178495 (218)
Cerast_456C/T
(90)
F: CAGCTTGGCATAACGTCACC
R: TTATGCCTTGCGAATGTCCG
gactgactgactgactgactgaCAGAAGGATGCGGCATTGTFMN178490 (190)
Cerast_31A/G
(141)
F: CAGACCAGGCAAACACATCA
R: CAGCTCGCATTTGTTCCCTT
gactgactgactgactgactgacATGATTAAGCAAGCTACTGCTAGFMN178488 (241)
Cerast_1255A/T
(115)
F: AATCGTTCATCATGTCCCGC
R: AGTTGGTGTTCACAATTCCCC
gactgactgactgactgactgactgactgaATTGATTCGCAGTGTTTTGCTRMN178493 (192)
Cerast_586C/T
(154)
F: TGAATCTGTCCGCATCCTGA
R: GATAAAACCTAACAGGTGGGC
gactgactgactgactgactgactgactgactgactgacATGGGCATGCGCAAAGGFMN178491 (213)
Cerast_450A/T
(165)
F: TTCACTCCACAACGAATCCA
R: CCGGTACCCCAAACAATATAACA
gactgactgactgactgactgactgactgactgactgactCATTATAAATTCCTAGCGAGCAGAFMN178489 (250)
Cerast_2530T/A
(126)
F: TGTGATTTGTGTGGTGCTGT
R: GTGTCTTCTAACATGGCATGCT
gactgactgactgactgactgactgactgactgactgactgactgactTGAATTTTGGCATGTTTTTGCTCTAGRMN178496 (243)

Share and Cite

MDPI and ACS Style

Maroso, F.; Pérez de Gracia, C.; Iglesias, D.; Cao, A.; Díaz, S.; Villalba, A.; Vera, M.; Martínez, P. A Useful SNP Panel to Distinguish Two Cockle Species, Cerastoderma edule and C. glaucum, Co-Occurring in Some European Beds, and Their Putative Hybrids. Genes 2019, 10, 760. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10100760

AMA Style

Maroso F, Pérez de Gracia C, Iglesias D, Cao A, Díaz S, Villalba A, Vera M, Martínez P. A Useful SNP Panel to Distinguish Two Cockle Species, Cerastoderma edule and C. glaucum, Co-Occurring in Some European Beds, and Their Putative Hybrids. Genes. 2019; 10(10):760. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10100760

Chicago/Turabian Style

Maroso, Francesco, Celia Pérez de Gracia, David Iglesias, Asunción Cao, Seila Díaz, Antonio Villalba, Manuel Vera, and Paulino Martínez. 2019. "A Useful SNP Panel to Distinguish Two Cockle Species, Cerastoderma edule and C. glaucum, Co-Occurring in Some European Beds, and Their Putative Hybrids" Genes 10, no. 10: 760. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10100760

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop