Development and characterization of thirteen novel microsatellite markers for use in Greenland sharks (Somniosus microcephalus), with cross-amplification in Pacific sleeper sharks (Somniosus pacificus)
BMC Research Notes volume 14, Article number: 28 (2021)
The objectives of this work are to isolate, develop, and characterize polymorphic microsatellite markers for use in Greenland sharks (Somniosus microcephalus).
Thirteen microsatellite loci were successfully amplified and yielded multi-locus genotypes for 36 S. microcephalus individuals from Grise Fjord (n = 16) and Svalbard (n = 20). Each locus yielded between 2 and 9 alleles and observed heterozygosity ranged from 0.11 to 0.70 when estimated across both sites. One locus and three loci deviated from HWE following Bonferroni correction, for individuals sampled from Grise Fjord and Svalbard, respectively. Cross-amplification was successful at every locus for five of the ten S. pacificus individuals.
Greenland sharks (Somniosus microcephalus) are long lived  and presumably late to mature  sharks capable of extensive migration  in the North Atlantic to Arctic marine environments. Previous molecular genetic work successfully differentiated Greenland sharks from other species in Somniosus  and revealed hybridization with Pacific sleeper sharks (Somniosus pacificus) [5, 6], however, knowledge of their population-level genetic variation is yet to be described.
Highly variable molecular genetic markers, such as microsatellites, can provide data to characterize population genetic structure, and have proven useful for describing spatial genetic variation in widely distributed elasmobranchs [7,8,9]. Here, we isolate, develop, and describe polymorphic microsatellite loci to provide molecular genetic markers for assessing the population-level genetic variation in the Greenland shark. Thirteen markers were identified as possible candidate loci for population analyses in Greenland shark samples collected from two locations in the Arctic. We also explored marker cross-amplification in Pacific sleeper sharks, providing evidence of utility in other Somniosus species.
Primer sets for 50 candidate microsatellite loci were constructed from reduced-representation genomic DNA sequence data from a single Somniosus microcephalus individual sampled from Resolute Bay (SRA accession: PRJNA655731). Genomic DNA from this individual was prepared using the QIAgen DNeasy Tissue Kit (Valencia, CA, USA) and shipped to Global Biologics (Columbia, MO) where a library was created with SPRI selection targeting 450 bp inserts. The library was pooled with other species libraries and sequenced on an Illumina HiSeq 2500. Approximately 4 million 250 bp PE reads were recovered; from these, tetrasat motifs with a minimum of five uninterrupted repeating motifs (e.g. GACAGACAGACAGACA) were targeted using STR finder in Galaxy  and PCR primer regions flanking each motif were identified using PRIMER3 . A total of 36 S. microcephalus sampled from two sites (Grise Fjord: n = 16; Svalbard: n = 20) were used to characterize and optimize the prospective microsatellite markers, and an additional ten S. pacificus individuals were used to test for cross-species amplification. For all individuals, genomic DNA was extracted and purified from fin tissue stored in either ethanol or RNAlater (Invitrogen, Carlsbad, CA, USA), using the Promega Wizard Genomic DNA Extraction Kit (Promega Corp, Madison, WI, USA), according to the manufacturer’s instructions.
PCR reactions occurred in 12.5 µL volumes containing 1.25 µL 10X PCR reaction buffer with 15 mM magnesium chloride (GenScript, Piscataway, NJ, USA), 0.25 µL 40 mM dNTP’s (APEX BioResearch Products), 0.25 µL 10 µM forward primer, 0.25 µL 10 µM reverse primer, 5 U Taq Polymerase (GenScript, Piscataway, NJ, USA), and 0.5 µL genomic DNA, with one locus (Smic2) requiring an additional 0.625 µL 15 mM magnesium chloride per reaction. Thermal-cycler conditions consisted of an initial denaturation at 94 °C for 2 min, followed by 30 cycles of 94 °C for 30 s (denaturation), 52–57 °C for 30 s (annealing) depending on marker identity (Table 1), 72 °C for 1 min (extension), and a final step at 72 °C for 1 min and 30 s. Three loci (Smic2, Smic24, Smic31) required half the time within each cycle (i.e., reduced from 30 to 15 s and 1 min to 30 s). PCR products were visualized on a 1.5% agarose gels for estimated size confirmation, then measured for precise fragment lengths using a Fragment Analyzer (Advanced Analytic Technologies, Inc., Ankeny, IA, USA) according to the manufacturer’s protocol. Fragment lengths were scored using PROsize 3.0 software (Advanced Analytic Technologies, Inc., Ankeny, IA, USA), and were binned across samples to account for variation between plate runs. All loci were analyzed for null alleles, stuttering, or large allele dropout using Micro-Checker 2.2.3 . The total number of alleles and observed and expected heterozygosity were calculated in GenAlEx 6.512b . Deviations from Hardy–Weinberg equilibrium (HWE) and detection of linkage disequilibrium among all loci pairs was determined in Genepop 4.2  for each site separately.
Results and discussion
Of the microsatellite loci screened in Somniosus microcephalus samples, 24 successfully amplified consistently and yielded indications of allelic polymorphism on agarose gels. Following initial agarose screening, 13 loci were further identified through fragment analysis as having clear fragment length peaks, and were then scored among 36 individuals. All other loci were deemed not suitable due to either monomorphic allele scoring, weak amplification, or the production of multiple (> 2) fragments following PCR and thermal-cycler optimizations (Additional file 1: Table S1).
Each locus produced two to seven alleles within each site, with observed and expected heterozygosity ranging from 0.063 to 0.750 and 0.117 to 0.770, respectively (Table 2). Four loci failed to conform to HWE following Bonferroni correction for multiple tests: Smic 24 at GF, and Smic1, Smic4, Smic16 at SV. Of the 13 loci, a single locus (Smic5) displayed signs of stuttering and possible scoring error, and homozygote excess in both sampling sites. Two additional loci (Smic4, Smic24) displayed homozygote excess in Grise Fjord samples as indicated by Microchecker.
Exact tests performed for each site resulted in deviations from HWE for individuals at one locus from Grise Fjord and three loci from Svalbard (Table 2). Linkage disequilibrium was not present at any locus pair in either site following Bonferroni correction. Cross-amplification at each of the 13 loci was successful in 5 of the 10 S. pacificus samples using the same cycling conditions. Four samples failed to amplify at 1 locus (2 samples at Smic15; 1 sample at Smic18; 1 sample at Smic1) and one sample at 2 loci (Smic1 and Smic18).
Due to the low sample sizes, it is not possible to make statements regarding the biological significance versus sampling artifacts for departures from HWE, homozygosity excess, and presence of null alleles. Nonetheless, these markers will be useful for exploring genetic variation and stock structure in Greenland sharks throughout their distribution. Furthermore, cross-amplification in the closely related Pacific sleeper shark (Somniosus pacificus) extends potential utility to other sleeper shark species.
Low sample size and/or a priori site-based population designation may be contributing to deviations from HWE, excess homozygosity, and presence of null alleles.
Availability of data and materials
The DNA sequences containing the microsatellite loci from which these primers were developed were submitted to Genbank (SRA accession: PRJNA655731). The datasets generated and/or analyzed during this study are available from the corresponding author upon request.
Polymerase chain reaction
- NA :
Number of alleles
- HO :
- HE :
- TA :
Nielson J, Hedeholm RB, Heinemeier J, Bushnell PG, Christiansen JS, Olsen J, et al. Eye lens radiocarbon reveals centuries of longevity in the Greenland shark (Somniosus microcephalus). Science. 2016;353(6300):702–4.
Yano K, Stevens JD, Compagno LJV. Distribution, reproduction and feeding of the Greenland shark Somniosus (Somniosus) microcephalus, with notes on two other sleeper sharks, Somniosus (Somniosus) pacificus and Somniosus (Somniosus) antarcticus. J Fish Biol. 2007;70:374–90.
Campana SE, Fisk AT, Klimley AP. Movements of Arctic and northwest Atlantic Greenland sharks (Somniosus microcephalus) monitored with archival satellite pop-up tags suggest long-range migrations. Deep-Sea Res. 2015;II(115):109–15.
Murray BW, Wang JY, Yang S, Stevens JD, Fisk A, Svavarsson J. Mitochondrial cytochrome b variation in sleeper sharks. Mar Biol. 2008;153:1015–22.
Hussey NE, Cosandey-Godin A, Walter RP, Hedges KJ, VanGerwen-Toyne M, Barkley AN, et al. Juvenile Greenland sharks Somniosus microcephalus (Bloch & Schneider, 1801) in the Canadian Arctic. Polar Biol. 2015;38(4):493–504.
Walter RP, Roy D, Hussey NE, Stelbrink B, Kovacs KM, Lydersen C, et al. Origins of the Greenland shark (Somniosus microcephalus): Impacts of ice-olation and introgression. Ecol Evolution. 2017;7:8113–25.
Karl SA, Castro ALF, Lopez JA, Charvet P, Burgess GH. Phylogeography and conservation of the bull shark (Carcharhinus leucas) inferred from mitochondrial and microsatellite DNA. Conserv Genet. 2011;12:371–82.
Vignaud TM, Maynard JA, Leblois R, Meekan MG, Vázquez-Juárez R, Ramírez-Macías D, et al. Genetic structure of populations of whale sharks among ocean basins and evidence for their historic rise and recent decline. Mol Ecol. 2014;23:2590–601.
Bester-van der Merwe AE, Bitalo D, Cuevas JM, Ovenden J, Hernández S, Da Silva S, McCord M, Roodt- Wilding R. Population genetics of Southern Hemisphere tope shark (Galeorhinus galeus): Intercontinental divergence and constrained gene flow at different geographical scales. PloS one. 2017;12(9):e0184481.
Afgan E, Baker D, Batut B, van den Beek M, Bouvier D, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018;46:537–44.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG. Primer3-new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. Microchecker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28:2537–9.
Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Heredity. 1995;86:248–9.
Aaron Fisk, Kit Kovaks, Christian Lydersen and Nigel Hussey for providing Greenland shark samples. Stacey McIntyre for assistance in molecular methods. Marcos Caprile and Anthony Nguyen for screening and comparing candidate sequence reads. Sharon Wildes and Cindy Tribuzio provided Pacific sleeper shark samples.
Funding for this work was provided to RPW by California State University, Fullerton startup funds.
Ethics approval and consent to participate
Consent for publication
All authors consent to the current version of the manuscript for submission.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Eleven additional microsatellite loci primer sets that were not further developed due to either lack of polymorphism or greater than two peaks during PCR amplification screening.
About this article
Cite this article
Swintek, M.A., Walter, R.P. Development and characterization of thirteen novel microsatellite markers for use in Greenland sharks (Somniosus microcephalus), with cross-amplification in Pacific sleeper sharks (Somniosus pacificus). BMC Res Notes 14, 28 (2021). https://doi.org/10.1186/s13104-021-05447-5