Research note | Open | Published:
Development and characterization of 17 polymorphic microsatellite markers for the reef manta ray (Mobula alfredi)
BMC Research Notesvolume 12, Article number: 233 (2019)
Limited sample sizes are often a problem for species of conservation concern when using genetic tools to make population assessments. Lack of analytical power from small sample sizes can be compensated for by use of a large marker set. Here we report on development and characterization of 17 novel microsatellite markers for the reef manta ray (Mobula alfredi).
Loci were screened on 60 reef manta rays (M. alfredi) sampled from the east coast of Australia. The number of alleles per locus varied from 2 to 13 with observed heterozygosities ranging between 0.300 and 0.917. The development of these 17 additional markers increases the total number of microsatellite markers available for this species to 27.
The reef manta ray (Mobula alfredi) occurs in tropical and subtropical waters of the Indian and Pacific oceans . The species forms aggregations in coastal waters throughout its range, which provides for targeted ecotourism in many coastal regions . However, these aggregations may suffer depletion from competing by-catch and targeted fisheries, particularly for their highly-valued body parts that are used in ‘traditional medicines’ [3,4,5,6]. Like many large shark and ray species, M. alfredi is long-lived with conservative life history parameters (i.e. slow growth, late maturity, low fecundity) that make populations highly-susceptible to over-exploitation . As such, the global conservation concerns for this species have resulted in its classification as ‘Vulnerable to Extinction’ on the IUCN Red List , and has led to listing in Appendix II of the Convention on the International Trade in Endangered Species and Appendices I and II of the Convention on Migratory Species. Future research priorities identified for reef manta rays include investigation of population structure (regional and global) and connectivity, based on genetic analyses .
To date, the limited genetic work conducted on DNA from reef manta rays has been primarily in the context of broader phylogenetic studies [10, 11]. One recent phylogenetic study resulted in the collapse of the genus Manta into Mobula, and a reassignment of the position of M. alfredi within this genus . An earlier study used nuclear DNA obtained from M. alfredi sampled in Japanese waters to develop ten species specific microsatellite markers . However, preliminary studies on of M. alfredi from the east coast of Australia found that only eight of the previously developed markers were polymorphic, with the remaining markers either non-polymorphic or non-amplifiable. Collection of adequate sample sizes to overcome the limitation of few markers is often difficult, expensive and inappropriate for species of conservation concern. Therefore, we address this issue through the development of additional microsatellite markers for M. alfredi that may allow for robust estimation of population sizes without the need to sample large numbers of animals, and will facilitate future investigations into fine-scale structure (e.g. relatedness).
For studies focused on non-model species of conservation concern, microsatellite markers can be amplified from relatively small quantities of DNA and provide a relatively low cost per sample for genotyping (especially when working with pre-developed marker sets). While genotype-by-sequencing approaches (i.e. RADseq, DArTseq) are increasingly used for genomic applications to non-model species, the strict requirements surrounding DNA quality, quantity and contamination may limit application to historical, degraded or small amounts of sampled-DNA. Furthermore, microsatellite markers provide the potential for a standardised marker set that can facilitate comparisons of genetically distinct populations across broad temporal and geographical scales. Here we report the development of 17 polymorphic microsatellite markers for M. alfredi, to augment the 10 previously developed, and detail marker performance through the genotyping of 60 individual manta rays from the east coast of Australia. We anticipate that microsatellite markers developed here will inform ongoing regional studies on population connectivity, and may facilitate global population comparisons through the provision of an enlarged and published marker set.
Methods and results
Microsatellite primer sequences were developed de novo from DNA obtained from a Mobula alfredi voucher specimen from Lady Elliot Island Reef, Australia. A genomic library was developed with the TruSeq Nano DNA sample prep kit (Illumina, San Diego, CA) and sequenced using an Illumina HiSeq 2000 platform (San Diego, CA) that supported the acquisition of 2 × 125 base pair (bp) paired-end reads, following the protocols supplied by the manufacturer. Sequences between 150 and 400 bp were explored for microsatellite motifs using the software QDD v.3.1  following the microsatellite discovery protocol described by Vargas-Caro et al. . In brief, perfect microsatellite repeat motifs with > 10 repeating units were identified and loci with product sizes estimated at less than 100 bp were removed. Additional loci where flanking regions contained ‘runs’ (of a single base) or repeated sequence elements were excluded. Selected sequences consisting of repeat motifs and flanking sequences were blasted against NCBI GenBank using default parameters to exclude loci potentially located in coding regions. Primers with homology outside the target flanking region or that may amplify more than one locus were excluded by blasting back against the genomic library using GENEIOUS v 9.1.8 (http://www.geneious.com) . Of the remaining candidate primer pairs, we selected 48 to take through to development, giving priority to loci with tri or tetra nucleotide repeat motifs. For the 48 candidate microsatellite loci developed, a ‘CAG-tail’ was attached to the 5′ end of the forward primers allowing downstream use of fluorescent PCR product labelling. Loci were combined in multiplexed PCRs to avoid overlapping product sizes. The 5′ end of the reverse primer had a ‘GTTT-tail’ attached to ensure complete adenylation, thereby protecting PCR products from non-template nucleotide addition and facilitating accurate genotyping .
Genotyping was conducted using genomic DNA extracted from ethanol preserved tissue biopsies collected non-lethally (2015–2018) from M. alfredi in waters surrounding Lady Elliot Island (− 24.11°, 152.71°) and North Stradbroke Island (− 27.42°, 153.54°), Queensland, Australia. Possible sample duplication was avoided, as each manta ray’s unique ventral skin pattern was photographed at the time of biopsy collection . Loci were amplified in a 12 µl PCR reaction comprised of 1–2 µl of genomic DNA (10–15 ng), 0.5 µl primer stock, 6 µl of 2 × MyTaq mix (Bioline Australia), with remaining volume made up to 12 µl with Milli-Q H2O. PCR cycling conditions were a 95 °C for 60 s step, followed by 38 cycles of 95 °C for 60 s, 95 °C for 15 s, 55 °C for 15 s and 72 °C for 15 s, with a final extension step of 72 °C for 15 min. The diluted PCR products (1:10 dilution) were denatured with formamide and an internal size standard (Genescan 500-LIZ) was added. Fragment separation was performed with an ABI 3730 DNA Analyser. Allele scoring was undertaken using the GENEIOUS Microsatellite Plugin Version 1.4 (GENEIOUS v 9.1.8) . Loci were initially amplified in singleplex reaction against 6–12 individuals. Loci that successfully amplified were characterised against an additional 18–24 individuals to determine information content before selection for inclusion in a multiplex PCR. After removing loci that either did not amplify, produced an inconsistent product, or were mono-morphic, the 48 candidate loci were reduced down to a final set of 17 (Table 1).
A total of 60 M. alfredi from the Australian east coast population were genotyped with these 17 loci. To ensure there were no scoring errors or allelic dropout, all traces were checked using MICRO-CHECKER . Observed and expected heterozygosities were calculated using GENALEX v 6.5 . Locus MaQDD04 had a tri-nucleotide motif with the remaining markers exhibiting di-nucleotide repeats. The total number of alleles per locus ranged from 2 (MaQDD23) to 13 (MaQDD21 & MaQDD31). MaQDD04 displayed the highest observed heterozygosity (Ho = 0.88) with heterozygosity ranging widely for di-nucleotide loci (Ho = 0.35–0.80). Exact tests identified two loci (MaQDD23, p = 0.0212 & MaQDD39, p = 0.036) not in Hardy–Weinberg equilibrium, as calculated by GENEPOP v 4.2 . Linkage disequilibrium tests (GENEPOP v4.2) revealed significant linkage between only one pair of loci (MaQDD04 & MaQDD11, p = 0.00), after Bonferroni correction.
The need for high quality, standardized genetic markers resulted in the development of 17 novel microsatellite for reef manta rays. This marker set supplements microsatellite primers already available for the species , and enables research questions to be addressed that require the higher statistical power provided by a larger marker set. It is hoped that these additional markers also facilitate collaboration between studies of reef manta ray populations globally through the provision of named primers for which resultant genotypes and diversity metrics can be collaboratively compared and contrasted for reef manta ray populations globally. The combination of these novel markers with those already developed previously for reef manta rays  provides up to 27 markers for use in population structure and genetic effective population size investigations. The delineation of population structure and generation of population size estimates have been highlighted as research priorities to help guide conservation action for the reef manta ray .
At the time of primer development and manuscript preparation, the study did not have access to a sufficient number of samples from other locations for cross-population comparisons. Recent genotyping suggests marker applicability for M. alfredi sampled on Australia’s west coast (approximately 5600 km from east coast sampling locations) (author’s unpublished data). Due to sample availability and permit restrictions, this study did not test cross-species (within the family Mobulidae) primer amplification for the newly developed microsatellite loci.
Marshall AD, Compagno LJV, Bennett MB. Redescription of the genus Manta with resurrection of Manta alfredi (Krefft, 1868) (Chondrichthyes; Myliobatoidei; Mobulidae). Zootaxa. 2009;2301:1–28.
O’Malley MP, Lee-Brooks K, Medd HB. The global economic impact of manta ray watching tourism. PLoS ONE. 2013. https://doi.org/10.1371/journal.pone.0065051.
Couturier LIE, Marshall AD, Jaine FRA, Kashiwagi T, Pierce SJ, Townsend KA, et al. Biology, ecology and conservation of the Mobulidae. J Fish Biol. 2012. https://doi.org/10.1111/j.1095-8649.2012.03264.x.
Croll DA, Dewar H, Dulvy NK, Fernando D, Francis MP, Galván-Magaña F, et al. Vulnerabilities and fisheries impacts: the uncertain future of manta and devil rays. Aquat Conserv. 2016. https://doi.org/10.1002/aqc.2591.
O’Malley MP, Townsend KA, Hilton P, Heinrichs S, Stewart JD. Characterization of the trade in manta and devil ray gill plates in China and South-east Asia through trader surveys. Aquat Conserv. 2017. https://doi.org/10.1002/aqc.2670.
Steinke DI, Bernard AM, Horn RL, Hilton P, Hanner R, Shivji MS. DNA analysis of traded shark fins and mobulid gill plates reveals a high proportion of species of conservation concern. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-10123-5.
Dulvy NK, Pardo SA, Simpfendorfer CA, Carlson JK. Diagnosing the dangerous demography of manta rays using life history theory. Peer J. 2014. https://doi.org/10.7717/peerj.400.
Marshall AD, Kashiwagi T, Bennett MB, Deakos M, Stevens G, McGregor F, et al. Mobula alfredi (amended version of the 2011 assessment). IUCN Red List; 2018. http://dx.doi.org/10.2305/IUCN.UK. Accessed 13 Nov 2018.
Stewart JD, Jaine FRA, Armstrong AJ, Armstrong AO, Bennett MB, Burgess KB, et al. Research priorities to support effective manta and devil ray conservation. Front Mar Sci. 2018. https://doi.org/10.3389/fmars.2018.00314.
Kashiwagi T, Marshall AD, Bennett MB, Ovenden JR. The genetic signature of recent speciation in manta rays (Manta alfredi and M. birostris). Mol Phylogenet Evol. 2012. https://doi.org/10.1016/j.ympev.2012.03.020.
Walter RP, Kessel ST, Alhasan N, Fisk AT, Heath DD, Chekchak T, et al. First record of living Manta alfredi × Manta birostris hybrid. Mar Biodivers. 2014. https://doi.org/10.1007/s12526-013-0183-2.
White WT, Corrigan S, Yang L, Henderson AC, Bazinet AL, Swofford DL, et al. Phylogeny of the manta and devilrays (Chondrichthyes: mobulidae), with an updated taxonomic arrangement for the family. Zool J Linn Soc. 2018. https://doi.org/10.1093/zoolinnean/zlx018.
Kashiwagi T, Broderick D, Lance SL, Bennett MB, Ovenden JR. Development and characterization of ten microsatellite loci for the reef manta ray Manta alfredi. Conserv Genet Resour. 2012. https://doi.org/10.1007/s12686-012-9705-7.
Meglécz E, Pech N, Gilles A, Dubut V, Hingamp P, Trilles A, et al. QDD version 3.1: a user-friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate. Mol Ecol Resour. 2014. https://doi.org/10.1111/1755-0998.12271.
Vargas-Caro C, Bustamante C, Bennett MB, Ovenden JR. Towards sustainable fishery management for skates in South America: The genetic population structure of Zearaja chilensis and Dipturus trachyderma (Chondrichthyes, Rajiformes) in the south-east Pacific Ocean. PLoS ONE. 2017. https://doi.org/10.1371/journal.pone.0172255.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012. https://doi.org/10.1093/bioinformatics/bts199.
Brownstein MJ, Carpten JD, Smith JR. Modulation of Non-Templated Nucleotide Addition by Taq DNA Polymerase: Primer Modifications that Facilitate Genotyping. Biotechniques. 1996. https://doi.org/10.2144/96206st01.
Marshall AD, Dudgeon CL, Bennett MB. Size and structure of a photographically identified population of manta rays Manta alfredi in southern Mozambique. Mar Biol. 2011. https://doi.org/10.1007/s00227-011-1634-6.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004. https://doi.org/10.1111/j.1471-8286.2004.00684.x.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012. https://doi.org/10.1093/bioinformatics/bts460.
Raymond M, Rousset F. GENEPOP (Version 12): population genetics software for exact tests and ecumenicism. J Hered. 1995. https://doi.org/10.1093/oxfordjournals.jhered.a111573.
AJA carried out primer design filtering, genotyping, data analysis and manuscript writing; CB advised and assisted with primer design; CLD assisted with genotyping and was a major contributor in writing the manuscript; JRO assisted with primer design, genotyping, and was a major contributor in writing the manuscript; MBB provided funding, permits, ethical clearance, and feedback on the manuscript. All authors read and approved the final manuscript.
The authors wish to acknowledge the services of the Australian Equine Genetics Research Centre Genotyping & Sequencing Facility (AEGRC).
The authors declare that they have no competing interests.
Availability of data and materials
The datasets analysed during the current study are not at present publicly available due ongoing publication, but are available from the corresponding author upon reasonable request. Sequences have been uploaded to GenBank and can be found at their accession numbers listed in Table 1.
Consent for publication
Ethics approval and consent to participate
All research was conducted in accordance with the approval of the University of Queensland’s Animal Ethics Committee (Permit SBS/319/14/ARC/EA/LEIER), the Great Barrier Reef Marine Park Authority (Permit G16/37856.1), the Department of National Parks, Sports and Racing (Permit QS2015/CVL1440) and the Department of Agriculture and Fisheries (Permit 199045).
This work was supported by an Australian Research Council Linkage Grant (LP150100669). AA is supported by a University of Queensland Research Scholarship. The funders were not involved in the design, collection, analysis and interpretation of the data, nor in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.