Skip to main content

Microsatellite cross-species amplification and utility in southern African elasmobranchs: A valuable resource for fisheries management and conservation



Similarly to the rest of the world, southern Africa’s diverse chondrichthyan fauna is currently experiencing high fishing pressures from direct and non-direct fisheries to satisfy market demands for shark products such as fins and meat. In this study, the development of microsatellite markers through cross-species amplification of primer sets previously developed for closely related species is reported as an alternative approach to de novo marker development. This included the design of four microsatellite multiplex assays and their cross-species utility in genetic diversity analysis of southern African elasmobranchs. As this study forms part of a larger project on the development of genetic resources for commercially important and endemic southern African species, Mustelus mustelus was used as a candidate species for testing these multiplex assays in down-stream applications.


Thirty five microsatellite primer sets previously developed for five elasmobranch species were selected from literature for testing cross-species amplification in 16 elasmobranch species occurring in southern Africa. Cross-species amplification success rates ranged from 28.6%-71.4%. From the successfully amplified microsatellites, 22 loci were selected and evaluated for levels of polymorphism, and four multiplex assays comprising of the 22 microsatellites were successfully constructed, optimised and characterised in a panel of 87 Mustelus mustelus individuals. A total of 125 alleles were observed across all loci, with the number of alleles ranging from 3–12 alleles. Cross-species amplification of the four optimised multiplex assays was further tested on 11 commercially important and endemic southern African elasmobranch species. Percentage of polymorphism ranged from 31.8%-95.5% in these species with polymorphic information content decreasing exponentially with evolutionary distance from the source species.


Cross-species amplification of the 35 microsatellites proved to be a time- and cost-effective approach to marker development in elasmobranchs and enabled the construction of four novel multiplex assays for characterising genetic diversity in a number of southern African elasmobranch species. This study successfully demonstrated the usefulness of these markers in down-stream applications such as genetic diversity assessment and species identification which could potentially aid in a more integrative, multidisciplinary approach to management and conservation of commercially important cosmopolitan and endemic elasmobranch species occurring in southern Africa.


The subclass Elasmobranchii (sharks, skates and rays) comprises a diverse group of over 1000 species, and is representative of one of the most ancient extant vertebrate lineages[1]. Recently, pressures from direct and non-direct fisheries have resulted in the depletion of elasmobranch populations globally[2]. Decline in wild populations of elasmobranchs is further compounded by their life history traits that are more similar to those of mammals (e.g. low fecundity, late maturity and long gestation periods) than those of teleost fishes[3, 4]. In comparison, elasmobranchs may not respond well to the high fishing pressures. This trend has been particularly pronounced for sharks due to unregulated harvesting to support an increase in demand for shark products (e.g. fins, meat, liver oil, skin and cartilage). A drastic reduction in population size (population bottleneck) can result in small populations experiencing the accumulating effects of inbreeding leading to severe loss of genetic diversity[5, 6]. These trends have previously been reported for species such as the basking shark (Cetorhinus maximus)[7] and the narrownose smoothhound shark (Mustelus schmitti)[8]. Assessing genetic diversity and population structure of wild populations is therefore important for sustainable long-term management of the global shark fishery industry.

Misidentification of shark species in fisheries operations is also a widespread concern[912], and molecular individual identification methods have been developed to alleviate this problem[4, 1317]. To integrate genetic knowledge with fisheries management, it is imperative for shark fisheries to report shark landings by species instead of lumping them into species- or family groups (e.g. houndsharks, carcharinids, hammerheads etc.). This stems from the difficulties involved with unambiguously identifying species within and across families[13, 18]e.g., carcharinids (Carcharhinus brachyurus, C. obscurus and C. plumbeus) and houndsharks (Mustelus mustelus, M. palumbes and Galeorhinus galeus) due to a high degree of conserved interspecific morphology[14]. Neglecting to report shark landings by species overlooks important differences in species susceptibility and population vulnerability to exploitation[15], and that in turn has important implications for species-specific conservation, management and trade monitoring programmes[19].

The general lack of molecular genetic markers (e.g. microsatellites) for many elasmobranch species impedes population and conservation genetic studies in that these markers can provide valuable information relating to population dynamics (spatial and temporal genetic variation) of individual species. Microsatellites are highly polymorphic due to their high mutation rate (between 10-3 and 10-4 mutations per gamete per generation) resulting in extensive length polymorphism[20, 21]. This makes microsatellite markers one of the most powerful molecular genetic tools with a remarkable array of applications ranging from genetic diversity[22, 23] and population structure inference[24, 25] to discerning genetic mating systems[26, 27] and the identification of species[2830].

Because the de novo development of microsatellites is challenging due to notoriously low rates of polymorphism in elasmobranchs[31], the development of microsatellite markers through cross-species amplification is the most effective alternative approach to de novo development of microsatellites and has recently also been reported in sharks[32]. Microsatellite cross-species amplification relies on the presence of conserved microsatellite flanking sequences[33], which in some organisms markedly demonstrate a high degree of conservation following millions of years of divergent evolution (e.g., 250 million years in sharks[28] and 470 million years in fish[34]). The success rate of microsatellite cross-species amplification has directly been correlated to the evolutionary distance between the source species and the target species[33, 35].

The underrepresentation of endemic taxa in many cross-species amplification studies is unfortunate as endemics should be of great interest for conservation of biodiversity on a regional scale. Southern Africa has one of the most diverse chondrichthyan faunas in the world, consisting of some 181 species in 44 families of which 34 species are endemic to southern Africa[36, 37]. Growing concerns regarding the sustainability of the southern African shark fishery, stemming from the local declines of cosmopolitan and endemic species, lead to stricter regulations being imposed so as to avert the collapse of natural populations[3840]. Accordingly, we report here the development of microsatellite markers through cross-species amplification of species-specific primers from closely related species. This included the design and optimisation of four microsatellite multiplex assays and their cross-species utility in genetic diversity analysis of 11 southern African elasmobranch species.

Results and discussion

Cross-species amplification

Development of microsatellite loci through cross-species amplification proved useful in establishing genetic markers for shark species that are commercially important and those (typically endemics) that are indirectly affected by fisheries’ operations. Amplification of the 35 microsatellites in 50 individuals from 16 different elasmobranch species (1–4 individuals per species) proved to be effective (Table 1). Cross-species amplification success rates or the percentage of microsatellites that amplified successfully ranged from 60.00%-71.40% in the Triakidae and Carcharhinidae families and from 28.57%-48.57% in the Scyliorhinidae, Sphyrnidae and Rajidae families (Figure 1). The higher success rates in the Triakidae and Carcharhinidae is expected as most microsatellites tested in this study were originally developed for species within the Triakidae family. Overall, the microsatellites showed less successful cross-species amplification to the taxa more divergent from the source species. Notably none of the individuals showed PCR amplification at any of the six Scyliorhinus canacula microsatellites. This may in part be attributed to S. canacula being more distantly related to the study species. The mean genetic distance between the taxa was 21.4 ± 1.7% (mean ± SD) (G. galeus as source species; Figure 2) and 18.7 ± 1.5% (M. canis as source species; Figure 3). Haploblepharus pictus could not be represented in the distance plot due to the lack of genetic information available in GenBank and Global Cestode Database: Elasmobranchs Specimens.

Table 1 Cross-species amplification of the 35 microsatellites among 16 elasmobranch species of southern Africa
Figure 1

Success rates of 35 microsatellite loci across five families of southern African elasmobranch species.

Figure 2

Cross-species amplification performance of Galeorhinus galeus microsatellites in 15 of the 16 elasmobranch species, and genetic divergence between G. galeus and the target species based on ND2 sequences.

Figure 3

Cross-species amplification performance of Mustelus canis microsatellites in 15 of the 16 elasmobranch species, and genetic divergence between M. canis and the target species based on ND2 sequences.

Results of cross-species amplification performance of G. galeus microsatellites exhibited a logarithmic regression function (Figure 2), albeit non-significant (slope within the 95% CI for no difference from zero), that may explain the general trend of negative correlation between cross-species amplification performance and genetic divergence seen across taxa[30, 33, 35, 41]. Cross-species amplification of microsatellite markers from source to target species is generally negatively correlated with evolutionary divergence[30, 4245]. A similar trend was not observed for the M. canis microsatellites (Figure 3). This can be due to different life history traits (i.e. mating system and generation time) and genome size (C value) between the source and target species, which have been previously found to have significant negative effects on cross-species amplification success[33]. However, apart from the source-target species evolutionary distance other factors, such as mutations in microsatellite flanking sequences, may affect the success rate of cross-species amplification. Since microsatellites are usually found in non-coding regions where the substitution rate is higher than in coding regions[46], these microsatellite flanking sequences which serve as regions for PCR primer design and binding sites are prone to mutations[35]. Mutations (indels) in these regions may therefore result in null alleles and in turn affect the patterns of cross-species amplification as demonstrated in birds[47] and salmonids[48].

Additionally, M. canis microsatellite loci were isolated from an enriched genomic library[49] whereas for G. galeus the microsatellites were developed by a high-throughput sequencing approach (Roche 454 pyrosequencing)[50]. Based on the observed data it is hypothesised that these different approaches may have influenced the cross-species performance possibly due to the different resolving power of each approach for capturing microsatellites distributed across different parts of the genome. Castoe et al.[51] argues that enrichment-based approaches commonly use a few specific repeated motifs, which are largely selected without prior knowledge of their abundance in the genome and therefore could introduce potential bias in genome representativeness. In contrast, microsatellite identification from randomly sequenced genomic regions (e.g. Roche 454 NGS) allows for an unbiased assessment of all types of microsatellite loci present in a genome[51].

Cross-species amplification of orthologous microsatellites, due to the presence of conserved microsatellite flanking sequences[33], can persist over millions of years following divergent evolution as previously shown in sharks (250 million years[28]) and in fish (470 million years[34]). This indicates lower mutation rates within microsatellite flanking regions in aquatic organisms[30, 34]. The reported cross-amplified microsatellite markers will thus contribute to establishing a molecular genetic marker repository for each of the Southern Africa elasmobranchs species included in this study. Future research efforts may be dedicated to generating microsatellite primers that have a high cross-species utility (e.g.[52]) as well as in silico mining of polymorphic microsatellite markers from expressed sequence tag data[53].

Multiplex assay characterisation

Twenty-two microsatellite loci that successfully cross-amplified across the study species and showed polymorphism in initial screening tests were used to develop four multiplex assays (MPS) comprised of at least five microsatellites each. These multiplex assays were characterised in a panel of 87 common smoothhound sharks (Mustelus mustelus) (Table 2). All 22 microsatellite loci were polymorphic across the multiplexes. With the exception of one locus, McaB22, all the microsatellite loci were found to deviate significantly from Hardy-Weinberg equilibrium most likely due to Wahlund effect as samples were pooled from diverse geographical locations for analysis. MicroChecker detected no significant genotyping errors but indicated that null alleles were present at two loci (Mh9 and Gg7). Slatkin’s exact test for neutrality indicated that two loci (McaB22 and Gg3) were candidates for being subjected to selection.

Table 2 Characterisation of four multiplex assays for Mustelus mustelus based on 87 individuals from southern Africa

Multiplex assay cross-species amplification and efficiency in species identification

Cross-species amplification of the four multiplex assays was tested for 11 other southern African shark species (Table 3). The number of alleles observed was highest in G. galeus and M. palumbes, varying from 1 to 7, while the percentage of polymorphism (PP) for each marker ranged from 31.8%-95.5%. The polymorphic information content (PIC) decreased exponentially with evolutionary distance from the source species (Table 3) and the four multiplex assays showed the highest PIC in M. mustelus, M. palumbes and G. galeus.

Table 3 Multiplex transferability results of a total of 22 microsatellite loci showing the number of alleles per locus for the 11 elasmobranch species tested

The mean genetic diversity estimates for each species in terms of number of alleles (A N ), effective number of alleles (A E ), observed heterozygosity (H O ), expected heterozygosity (H E ) and PIC are shown in Figure 4. In group 1, the mean HE, AE and PIC varied from relatively low in C. brachyurus (mean HE = 0.230; mean AE = 1.4; mean PIC = 0.181) to relatively high in M. palumbes (mean HE = 0.653; mean AE = 3.3; mean PIC = 0.606). Group 2 exhibited similar patterns of genetic diversity that varied from moderate in S. zygaena (mean HE = 0.593; mean AE = 3.2; mean PIC = 0.554) to relatively high in P. pantherinum (mean HE = 0.662; mean AE = 3.4; mean PIC = 0.603). For group 3 with n = 4, the mean HE, AE and PIC ranged from relatively low in C. plumbeus (mean HE = 0.249; mean AE = 1.5; mean PIC = 0.193) to relatively high in C. obscurus (mean HE = 0.429; mean AE = 2.1; mean PIC = 0.367).

Figure 4

Mean genetic diversity estimates using 12 microsatellite loci shared between species: number of alleles ( A N ), effective number of alleles ( A E ), heterozygosity ( H E ) and polymorphic information content ( PIC ). Error bars represent standard error.

The genotypic distribution of the study species is depicted in Figure 5. Most of the study species could be differentiated on PC1 and PC2 of the PCoA plot as can be seen from individuals of each respective species clustering together. Individuals of the catshark species (H. edwardsii, P. pantherinum and H. pictus) however, were dispersed across quadrant 3 and 4. The PCoA also revealed that one of the M. mustelus individuals was misidentified as G. galeus. The identity of this particular individual was subsequently confirmed using the genetic identification method developed specifically for smoothhound sharks[16]. Briefly, this method involves using four primers (1 universal forward primer and 3 species-specific reverse primers) for the mitochondrial gene, nicotinamide adenine dehydrogenase subunit 2 (ND2), in a multiplex PCR reaction. The reverse primers amplify a fragment of different length for each species (M. asterias, 564 bp; M. mustelus, 392 bp; G. galeus, 671 bp) and can therefore be utilised for distinguishing species based on fragment size.

Figure 5

Principle coordinates analysis (PCoA) of study species based on 12 shared amplified microsatellite loci between species. Arrows depict misidentified/mislabelled individuals.

In this study, the potential use of microsatellite loci in species identification was successfully demonstrated using shared microsatellite loci between species. The polymorphic information of these microsatellite loci was characterised by low genetic variation as previously proposed for elasmobranchs[31]. The genotypic distribution of the study species could also be differentiated based on PCoA analysis. Markedly, the lack of differentiation between the catshark species (H. edwardsii and H. pictus) on the PCoA plot may be explained by the misidentification of the Haploblepharus species that is a common occurrence in the catshark family[54]. To further investigate whether the lack of differentiation detected with microsatellites was indeed due to misidentification, the cytochrome b (Cyt b) and ND2 genes were applied for species comparisons[54]. For both genes, sequence analyses revealed individuals with cryptic identification, suggesting that in the case of the catsharks, the microsatellites optimised in the current study were not successful in discriminating between the Haploblepharus species but could in the future aid in the identification of cryptic speciation within the catshark family.

In South Africa, the aforementioned misidentification issue is prominent in fishing operations (particularly in longline and trawl fisheries) where there is a high rate of incidental capture of non-target shark species[912]. This hinders the collection of reliable data on shark catch and trade on a species-specific basis making robust stock assessments and identification of overfished and potentially threatened species nearly impossible in most situations[4]. This was apparent in a study by Attwood et al.[12], which assessed bycatch in South Africa’s inshore trawl fishery based on observer records. In the aforementioned study, certain taxa were difficult to identify, and therefore taxonomic groups were lumped in species groups (e.g. Raja spp., Mustelus spp. and Scyliorhinidae), even though every attempt was made to analyse data at species level. Therefore, the molecular genetic makers developed by the current study may facilitate in obtaining species-specific catch data for stock assessment, characterising genetic diversity and delineating population genetic structure. This in turn will contribute to the implementation of future conservation and management plans on a species-specific level in southern Africa.


Cross-species amplification of available microsatellite loci to target species has proven to be more time- and cost-effective in comparison to the de novo development approach and permitted the cross-amplification of 22 markers across 12 elasmobranch species. Cross-species amplification of the four multiplex assays developed in the current study highlighted the usefulness of microsatellites for characterising genetic diversity and potentially also species identification of a number of commercially important and endemic elasmobranch species. The molecular genetic markers developed in this study and their usefulness in down-stream applications could therefore aid in a more integrative, multidisciplinary approach to conservation management of elasmobranchs in southern Africa.


Ethics statement

The collection of specimens from various shark species used in this study complied with the Convention on Biological Diversity ( and the Convention on the Trade in Endangered Species of Wild Fauna and Flora ( All permits to collect finclip or muscle tissue for research purposes were granted by the Department of Agriculture, Forestry and Fisheries (Republic of South Africa).

Study species and DNA extraction

Sixteen elasmobranch species occurring in southern African waters belonging to five families within two orders were selected for cross-species amplification (Additional file1). Where possible, specimens were collected from at least two sampling locations to better capture allelic diversity present within populations of each respective species. However, due to opportunistic sampling for a majority of the study species, samples were obtained from only a single location (Additional file1). Muscle tissue or finclips were preserved in 99% ethanol and stored at room temperature until further use. Total genomic DNA was isolated using the standard cetyltrimethylammonium bromide (CTAB) method of Saghai-Maroof et al.[55]. The extracted DNA was quantified using a NanoDrop ND-1000 spectrophotometer v.3.0.1 (NanoDrop®). For testing cross-species amplification, each DNA sample was adjusted to a working concentration of 50 ng/μl and stored at -20°C prior to polymerase chain reaction (PCR) analysis.

Microsatellite primer transfer

A total of 35 microsatellite markers previously developed in five elasmobranch species (Raja pulchra[44], Scyliorhinus canacula[45], M. canis[49], G. galeus[50] and M. henlei[56, 57]) were selected for testing cross-species amplification. Primer sequences and annealing temperature (TA) of each primer set optimised for each respective source species are shown in Table S2 (Additional file2). Polymerase chain reaction conditions optimised for the majority of the source species [G. galeus, M. henlei, R. pulchra and S. canacula] were applied for cross-species microsatellite examinations except for the M. canis ([56] PCR protocol) and some R. pulchra (modified[57] PCR protocol) primer sets. Polymerase chain reactions for all individuals were executed in a GeneAmp® PCR System 2700.

The PCR amplicons were visualised on a 2% agarose gel stained with ethidium bromide together with negative controls and Promega 100 bp molecular size ladder for preliminary size determination. Success or failure of PCR amplification in cross-species trials was determined simply on the basis of whether band intensity was sufficient to score alleles. In most instances, less stringent PCR conditions were not employed in the cross-species assays so as to minimise the risk of amplification of non-orthologous loci in the target species. The number of markers that showed amplification success at all or a percentage of individuals in the target species (“+/++” in Table 1) were counted as an index to measure the cross-species microsatellite amplification performance.

Multiplex design and optimisation

As this study forms part of a larger project on the development of genetic resources for commercially important and endemic species of southern Africa, Mustelus mustelus was used as a candidate species for testing of the four novel multiplex assays. Levels of polymorphism were initially assessed at all the successfully cross-amplified microsatellite loci in a panel of eight individuals of M. mustelus. Amplicons were subjected to electrophoresis for two hours at 150 volts on a 12% polyacrylamide gel to detect size variants. Microsatellites were considered to be polymorphic when two bands were distinguishable in a single individual (i.e. heterozygote) and/or clear size differences were detected between different individuals.

Twenty-two polymorphic microsatellite loci were selected, and primers fluorescently labelled and optimised in four multiplex assays (5–6 loci per MPS) using a strategy outlined by Guichoux et al.[58] with one of the following dyes: FAM, VIC, PET, or NED. The use of different dyes was to facilitate co-amplification of multiple microsatellite markers in a single reaction for cost- and time-efficient genotyping (Multiplex PCR).

After optimisation of the newly designed MPS (MPS1, MPS2, MPS3 and MPS4), a panel of 87 M. mustelus individuals from across the distribution range in southern Africa was genotyped for marker characterisation purposes. The multiplex assays were then tested on 11 additional species to show their overall application in genetic diversity and population structure analysis. Finally, a total of 12 microsatellite loci that were successfully genotyped across the study species (Mh1, Mh25, Mca25, McaB39, McaB5, McaB22, Mca33, Gg15, Gg17, Gg18, Gg22 and Gg23) were selected to demonstrate the potential use of microsatellite loci in species identification.

The percentage of polymorphism (PP) was calculated using the formula:

PP = N P N T X 100

where N P is the total number of polymorphic loci and N T is the total number of loci multiplied by 100.

For the multiplex reaction, the Qiagen Multiplex PCR kit was used and PCR conducted according to the manufacturer’s instructions except for varying TA, 59°C for MPS1, MPS3 and MPS4; and 56°C for MPS2. For subsequent analysis on an ABI 3730XL DNA Analyzer, PCR products were diluted in distilled water and fragment analysis performed together with the LIZ600 internal size standard. Individual genotypes were scored based on fragment size via Peak Scanner® software v.1 (Life Technologies). AutoBin v.0.9 macro for Excel (; see[58]) was used to detect discreet size variants where allele binning of genotype data obtained from Peak Scanner® software v.1 was based on raw size.

Genetic diversity analysis

MicroChecker v.2.2.3[59] was used to evaluate the presence of genotypic errors caused by allele dropout, stuttering and null alleles. Null allele frequencies (Fr NULL ) were calculated using the Brookfield 1 estimator implemented in this program. Locus-specific fixation index (F IS ) and over all loci was estimated to measure departure from Hardy-Weinberg equilibrium using the exact probability test (20 batches, Dememorization; 10000 and 5000 iterations) using Genepop v.4.0[60]. Linkage disequilibrium between all pairs of loci was calculated using an exact test implemented also in Genepop. Slatkin’s exact test (1000 permutations) for neutrality, based on Ewens-Watterson sampling theory[61] was used to detect loci under selection as implemented in Arlequin v.[62]. The number (A N ) of alleles at each microsatellite locus, as well as the effective number of alleles A E :

A E = 1 / i = 1 n p i 2

where p i is the frequency of the ith allele and n is the number of alleles was calculated using the GenAlEx v.6.5 program[63]. The proportion of individual samples that were heterozygous [direct count heterozygosity (H O ) and expected under Hardy-Weinberg equilibrium (H E )] was calculated using MsatTools[64]. MsatTools was also used to calculate the polymorphic information content (PIC) of each marker according to the following equation in[65]:

PIC = 1 - i = 1 n p i 2 - 2 i = 1 n - 1 j = i + 1 n p i 2 p j 2 where p i and p j are the frequency of the ith and jth allele respectively and n is the number of alleles.

Direct comparison of genetic diversity estimates (HE,A E and PIC) across the 11 species was not plausible due to the different sample sizes that were used. Species were therefore grouped into three groups according to sample size: (1) M. mustelus, M. palumbes, G. galeus, C. brachyurus and H. pictus (n = 8); (2) P. pantherinum, S. lewini and S. zygaena (n = 5) and (3) C. obscurus, C. limbatus, C. plumbeus and H. edwardsii (n = 4). The potential use of microsatellite data for species-assignment was assessed through principle coordinate analysis (PCoA) in GenAlEx v.6.5 using genetic distances between individuals.

To evaluate cross-species amplification performance, DNA sequences derived from the mitochondrial ND2 gene (1044 bp) of each species were downloaded from GenBank and Global Cestode Database: Elasmobranchs Specimens ( (Additional file3). The genetic distance of the study taxa was estimated using the Kimura 2-parameter model with the rate variation among sites modelled with a gamma distribution (shape parameter = 5) implemented in MEGA v.5[66].



Confidence interval

PP :

Percentage of polymorphism

P :



Polymorphic information content

A N :

Number of alleles

A E :

Effective number of alleles

H O :

Observed heterozygosity

H E :

Expected heterozygosity

F IS :

Inbreeding coefficient


Null allele frequency

P E-W :

Ewans-Watterson probability


Kimura 2-parameter model

ND2 :

Nicotinamide adenine dehydrogenase subunit 2


Mustelus mustelus


Mustelus palumbes


Galeorhinus galeus


Scylliogaleus quecketti


Carcharhinus brachyurus


Carcharhinus limbatus


Carcharhinus obscurus


Carcharhinus plumbeus


Haploblepharus pictus


Haploblepharus edwardsii


Poroderma africanum


Poroderma pantherinum


Sphyrna lewini


Sphyrna zygaena


Raja straeleni


Raja alba.


  1. 1.

    Corrigan S, Beheregaray LB: A recent shark radiation: molecular phylogeny, biogeography and speciation of wobbegong sharks (family: Orectolobidae). Mol Phyl Evol. 2009, 52: 205-216. 10.1016/j.ympev.2009.03.007.

    CAS  Article  Google Scholar 

  2. 2.

    Worm B, Davis B, Kettemer L, Ward-Paige CA, Chapman D, Heithaus MR, Kessel ST, Gruber SH: Global catches, exploitation rates, and rebuilding options for sharks. Mar Policy. 2013, 40: 194-204.

    Article  Google Scholar 

  3. 3.

    Stevens JD, Bonfil R, Dulvy NK, Walker PA: The effects of fishing on sharks, rays, and chimaeras (chondrichthyans), and the implications for marine ecosystems. ICES J Mar Sci. 2000, 57: 476-494. 10.1006/jmsc.2000.0724.

    Article  Google Scholar 

  4. 4.

    Shivji M, Clarke S, Pank M, Natanson L, Kohler N, Stanhope M: Genetic identification of pelagic shark body parts for conservation and trade monitoring. Conserv Biol. 2002, 16: 1036-1047. 10.1046/j.1523-1739.2002.01188.x.

    Article  Google Scholar 

  5. 5.

    Nei M, Maruyama T, Chakraborty R: Bottleneck effect and genetic variability in populations. Evolution. 1975, 29: 1-10. 10.2307/2407137.

    Article  Google Scholar 

  6. 6.

    Glenn TC, Stephan W, Braun MJ: Effects of a population bottleneck on whooping crane mitochondrial DNA variation. Conserv Biol. 1999, 13: 1097-1107. 10.1046/j.1523-1739.1999.97527.x.

    Article  Google Scholar 

  7. 7.

    Hoelzel AR, Shiviyi MS, Magnussen J, Francis MP: Low worldwide genetic diversity in the basking shark. Biol Lett. 2006, 2: 630-642.

    Google Scholar 

  8. 8.

    Pereyra S, García G, Miller P, Oviedo S, Domingo A: Low genetic diversity and population structure of the narrownose shark (Mustelus schmitti). Fish Res. 2010, 106: 468-473. 10.1016/j.fishres.2010.09.022.

    Article  Google Scholar 

  9. 9.

    Myers RA, Worm B: Rapid worldwide depletion of predatory fish communities. Nature. 2003, 423: 280-283. 10.1038/nature01610.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Barker MJ, Schluessel V: Managing global shark fisheries: suggestions for prioritising management strategies. Aquat Conserv Mar and Freshwat Ecosyst. 2005, 15: 325-347. 10.1002/aqc.660.

    Article  Google Scholar 

  11. 11.

    Petersen SL, Honig MB, Ryan PG, Underhill LG, Compagno LJV: Pelagic shark bycatch in the pelagic longline fishery off Southern Africa. Understanding and Mitigating Vulnerable Bycatch in Southern African Trawl and Longline Fisheries. Edited by: Petersen SL, Nel DC, Ryan PG, Underhill LG. 2008, WWF South Africa Report Series - 2008/Marine/002

    Google Scholar 

  12. 12.

    Attwood CG, Peterson SL, Kerwath SE: Bycatch in South Africa's inshore trawl fishery as determined from observer records. ICES J Mar Sci. 2011, 68: 2163-2174. 10.1093/icesjms/fsr162.

    Article  Google Scholar 

  13. 13.

    Pank M, Shivji MS, Stanhope M, Natanson L, Kohler N: Rapid and simultaneous identification of body parts from the morphologically similar sharks Carcharhinus obscurus and Carcharhinus plumbeus (Carcharhinidae) using multiplex PCR. Mar Biotechnol. 2001, 3: 231-240. 10.1007/s101260000071.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    Domingues RR, de Amorim AF, Hilsdorf AWS: Genetic identification of Carcharhinus sharks from the southwest Atlantic Ocean (Chondrichthyes: Carcharhiniformes). J Appl Ichthyol. 2013, 29: 738-742. 10.1111/jai.12154.

    CAS  Article  Google Scholar 

  15. 15.

    Abercrombie DL, Clarke SC, Shivji MS: Global-scale genetic identification of hammerhead sharks: Application to assessment of the international fin trade and law enforcement. Conserv Genet. 2005, 6: 775-788. 10.1007/s10592-005-9036-2.

    CAS  Article  Google Scholar 

  16. 16.

    Farrell ED, Clarke MW, Mariani S: A simple genetic identification method for northeast Atlantic smoothhound sharks (Mustelus spp.). ICES J Mar Sci. 2009, 66: 561-565. 10.1093/icesjms/fsn218.

    Article  Google Scholar 

  17. 17.

    Naylor GJP, Caira JN, Jensen K, Rosana KAM, White WT, Last PR: A DNA sequence based approach to the identification of shark and ray species and its implications for global elasmobranch diversity and parasitology. Bull Am Mus of Nat Hist. 2012, 367: 1-262.

    Article  Google Scholar 

  18. 18.

    Da Silva C, Bürgener M: South Africa’s demersal shark meat harvest. Traffic Bulletin. 2007, 21: 55-65.

    Google Scholar 

  19. 19.

    Clarke SC: Understanding pressures on fishery resources through trade statistics: A pilot study of four products in the Chinese dried seafood market. Fish Fish. 2004, 5: 53-74. 10.1111/j.1467-2960.2004.00137.x.

    Article  Google Scholar 

  20. 20.

    Tautz D: Hypervariability of simple sequences as a general source for polymorphic DNA markers. Nucleic Acids Res. 1989, 17: 6463-6471. 10.1093/nar/17.16.6463.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  21. 21.

    Weber JL, Wong C: Mutation of human short tandem repeats. Hum Mol Genet. 1993, 2: 1123-1128. 10.1093/hmg/2.8.1123.

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Edwards CJ, Soulsbury CD, Statham MJ, Ho SYW, Wall D, Dolf D, Iossa G, Baker PJ, Harris S, Sacks BN, Bradley DG: Temporal genetic variation of the red fox, Vulpes vulpes, across western Europe and the British Isles. Quat Sci Rev. 2012, 57: 95-104.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Ozerov MY, Veselov AE, Lumme J, Primmer CR: Temporal variation of genetic composition in Atlantic salmon populations from the Western White Sea Basin: influence of anthropogenic factors?. BMC Genet. 2013, 14: 88-

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Urian KW, Hofmann S, Wells RS, Read AJ: Fine-scale population structure of bottlenose dolphins (Tursiops truncatus) in Tampa Bay, Florida. Mar Mamm Sci. 2009, 25: 619-638. 10.1111/j.1748-7692.2009.00284.x.

    Article  Google Scholar 

  25. 25.

    Mourier J, Mills SC, Planes S: Population structure, spatial distribution and life‒history traits of blacktip reef sharks Carcharhinus melanopterus. J Fish Biol. 2013, 82: 979-993. 10.1111/jfb.12039.

    PubMed  CAS  Article  Google Scholar 

  26. 26.

    Mobley KB, Amundsen T, Forsgren E, Svensson PA, Jones AG: Multiple mating and low incidence of cuckoldry for nest-holding males in the two-spotted goby, Gobiusculus flavescens. BMC Evol Biol. 2009, 9: 6-10.1186/1471-2148-9-6.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Boomer JJ, Harcourt RG, Francis MP, Walker TI, Braccini JM, Stow AJ: Frequency of multiple paternity in gummy shark, Mustelus antarcticus, and rig, Mustelus lenticulatus, and the implications of mate encounter rate, postcopulatory influences, and reproductive mode. J Hered. 2013, 104: 371-379. 10.1093/jhered/est010.

    PubMed  Article  Google Scholar 

  28. 28.

    Martin AP, Pardini AT, Noble LR, Jones CS: Conservation of a dinucleotide simple sequence repeat locus in sharks. Mol Phyl Evol. 2002, 23: 205-213. 10.1016/S1055-7903(02)00001-5.

    CAS  Article  Google Scholar 

  29. 29.

    Costa M, Fernandes C, Rodrigues M, Santos-Reis M, Bruford MW: A panel of microsatellite markers for genetic studies of European polecats (Mustela putorius) and ferrets (Mustela furo). Eur J Wildl Res. 2012, 58: 629-633. 10.1007/s10344-012-0627-1.

    Article  Google Scholar 

  30. 30.

    Sekino M, Hara M: Individual assignment tests proved genetic boundaries in a species complex of Pacific abalone (genus Haliotis). Conserv genet. 2007, 8: 823-841. 10.1007/s10592-006-9229-3.

    CAS  Article  Google Scholar 

  31. 31.

    Dudgeon CL, Blower DC, Broderick D, Giles JL, Holmes BJ, Kashiwagi T, Krück NC, Morgan JAT, Tillett BJ, Ovenden JR: A review of the application of molecular genetics for fisheries management and conservation of sharks and rays. J Fish Biol. 2012, 80: 1789-1843. 10.1111/j.1095-8649.2012.03265.x.

    PubMed  CAS  Article  Google Scholar 

  32. 32.

    Boomer JJ, Stow AJ: Rapid isolation of the first set of polymorphic microsatellite loci from the Australian gummy shark, Mustelus antarticus and their utility across divergent shark taxa. Conserv Genet Resour. 2010, 2: 393-395. 10.1007/s12686-010-9274-6.

    Article  Google Scholar 

  33. 33.

    Barbara T, Palma-Silva C, Paggi GM, Bered F, Fay MF, Lexer C: Cross-species transfer of nuclear microsatellite markers: potential and limitations. Mol Ecol. 2007, 16: 3759-3767. 10.1111/j.1365-294X.2007.03439.x.

    PubMed  Article  Google Scholar 

  34. 34.

    Rico C, Rico I, Hewitt G: 470 million years of conservation of microsatellite loci among fish species. Proc R Soc Lond B. 1996, 263: 549-557. 10.1098/rspb.1996.0083.

    CAS  Article  Google Scholar 

  35. 35.

    Primmer CR, Painter JN, Koskinen MT, Palo JU, Merilä J: Factors affecting avian cross-species microsatellite amplification. J Avian Biol. 2005, 36: 348-360. 10.1111/j.0908-8857.2005.03465.x.

    Article  Google Scholar 

  36. 36.

    Compagno LJV: An overview of chondrichthyans systematics and biodiversity in Southern Africa. Trans R Soc S Afr. 2000, 54: 75-120.

    Article  Google Scholar 

  37. 37.

    Department of Agriculture, Forestry and Fisheries (DAFF): South Africa’s National Plan of Action for the Conservation and Management of Sharks. []

  38. 38.

    Navigating Global Shark Conservation: Current Measures and Gaps. []

  39. 39.

    Namibia's National Plan of Action to prevent, deter and eliminate illegal, unreported and unregulated fishing. []

  40. 40.

    Pierce SJ, Trerup M, Williams C, Tilley A, Marshall AD, Raba N: Shark fishing in Mozambique: A preliminary assessment of artisanal fisheries. 2008, Maputo: Eyes on the Horizon, 1-28.

    Google Scholar 

  41. 41.

    Karaiskou N, Buggiotti L, Leder E, Primmer CR: High degree of transferability of 86 newly developed zebra finch EST-linked microsatellite markers in 8 bird species. J Hered. 2008, 99: 688-693. 10.1093/jhered/esn052.

    PubMed  CAS  Article  Google Scholar 

  42. 42.

    FitzSimmons NN, Moritz C, Moore S: Conservation and dynamics of microsatellite loci over 300 million years of marine turtle evolution. Mol Biol Evol. 1995, 12: 432-440.

    PubMed  CAS  Google Scholar 

  43. 43.

    Primmer CR, Møller AP, Ellegren H: A wide-range survey of cross-species microsatellite amplification in birds. Mol Ecol. 1996, 5: 365-378. 10.1111/j.1365-294X.1996.tb00327.x.

    PubMed  CAS  Article  Google Scholar 

  44. 44.

    Kang J, Park J, Jo H: Rapid development of microsatellite markers with 454 pyrosequencing in a vulnerable fish, the mottled skate, Raja pulchra. Int J Mol Sci. 2012, 13: 7199-7211. 10.3390/ijms13067199.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  45. 45.

    Griffiths AM, Casane D, McHugh M, Wearmouth VJ, Sims DW, Genner MJ: Characterisation of polymorphic microsatellite loci in the small-spotted catshark (Scyliorhinus canicula L.). Conserv Genet Resour. 2011, 3: 705-709. 10.1007/s12686-011-9438-z.

    Article  Google Scholar 

  46. 46.

    Zane L, Bargelloni L, Patarnello T: Strategies for microsatellite isolation: a review. Mol Ecol. 2002, 11: 1-16. 10.1046/j.0962-1083.2001.01418.x.

    PubMed  CAS  Article  Google Scholar 

  47. 47.

    Primmer CR, Ellegren H: Patterns of molecular evolution in avian microsatellites. Mol Biol Evol. 1998, 15: 997-1008. 10.1093/oxfordjournals.molbev.a026015.

    PubMed  CAS  Article  Google Scholar 

  48. 48.

    Angers B, Bernatchez L: Complex evolution of a salmonid microsatellite locus and its consequences in inferring allelic divergence from size information. Mol Biol Evol. 1997, 14: 230-238. 10.1093/oxfordjournals.molbev.a025759.

    PubMed  CAS  Article  Google Scholar 

  49. 49.

    Giresi M, Renshaw MA, Portnoy DS, Gold JR: Isolation and characterization of microsatellite markers for the dusky smoothhound shark, Mustelus canis. Conserv Genet Resour. 2012, 4: 101-104. 10.1007/s12686-011-9484-6.

    Article  Google Scholar 

  50. 50.

    Chabot C, Nigenda S: Characterization of 13 microsatellite loci for the tope shark, Galeorhinus galeus, discovered with next generation sequencing and their utility for eastern Pacific smooth-hound sharks (Mustelus). Conserv Genet Resour. 2011, 3: 553-555. 10.1007/s12686-011-9402-y.

    Article  Google Scholar 

  51. 51.

    Castoe TA, Poole AW, Gu W, Jason De Koning AP, Daza JM, Smith EN, Pollock DD: Rapid identification of thousands of copperhead snake (Agkistrodon contortrix) microsatellite loci from modest amounts of 454 shotgun genome sequence. Mol Ecol Res. 2010, 10: 341-347. 10.1111/j.1755-0998.2009.02750.x.

    CAS  Article  Google Scholar 

  52. 52.

    Dawson DA, Horsburgh GJ, Küpper C, Stewart IR, Ball AD, Durrant KL, Burke T: New methods to identify conserved microsatellite loci and develop primer sets of high cross-species utility–as demonstrated for birds. Mol Ecol Res. 2010, 10: 475-494. 10.1111/j.1755-0998.2009.02775.x.

    CAS  Article  Google Scholar 

  53. 53.

    Hoffman JI, Nichols HJ: A novel approach for mining polymorphic microsatellite markers in silico. PLoS One. 2011, 6: e23283-10.1371/journal.pone.0023283.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  54. 54.

    Human BA: A taxonomic revision of the catshark genus Haploblepharus Garman 1913 (Chondrichthyes: Carcharhiniformes: Scyliorhinidae). Zootaxa. 2007, 1451: 1-40.

    Google Scholar 

  55. 55.

    Saghai-Maroof MA, Solima KM, Jorgenson RA, Allard RW: Ribosomal DNA spacerlength polymorphisms in barley: Mendelian inheritance, chromosomal location, and population dynamics. Proc Natl Acad Sci U S A. 1984, 81: 8014-8018. 10.1073/pnas.81.24.8014.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  56. 56.

    Chabot CL: Characterization of 11 microsatellite loci for the brown smooth-hound shark, Mustelus henlei, discovered with next-generation sequencing. Conserv Genet Resour. 2012, 4: 23-25. 10.1007/s12686-011-9464-x.

    Article  Google Scholar 

  57. 57.

    Byrne RJ, Avis JC: Genetic mating system of the brown smoothhound shark (Mustelus henlei), including a literature review of multiple paternity in other elasmobranch species. Mar Biol. 2012, 159: 749-756. 10.1007/s00227-011-1851-z.

    Article  Google Scholar 

  58. 58.

    Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, Lepoittevin C, Malaus T, Revardel E, Salin F, Petit RJ: Current trends in microsatellite genotyping. Mol Ecol Resour. 2011, 11: 591-611. 10.1111/j.1755-0998.2011.03014.x.

    PubMed  CAS  Article  Google Scholar 

  59. 59.

    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-538. 10.1111/j.1471-8286.2004.00684.x.

    CAS  Article  Google Scholar 

  60. 60.

    Rousset F: GENEPOP’007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.

    PubMed  Article  Google Scholar 

  61. 61.

    Slatkin M: An exact test for neutrality based on the Ewens sampling distribution. Genet Res. 1994, 64: 71-74. 10.1017/S0016672300032560.

    PubMed  CAS  Article  Google Scholar 

  62. 62.

    Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.

    CAS  PubMed Central  Google Scholar 

  63. 63.

    Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in Excel: population genetic software for teaching and research – an update. Bioinformatics. 2012, 28: 2537-2539. 10.1093/bioinformatics/bts460.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  64. 64.

    Park S: The Excel microsatellite toolkit. []

  65. 65.

    Botstein D, White RL, Skolnick M, Davis RW: Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980, 32: 314-331.

    PubMed  CAS  PubMed Central  Google Scholar 

  66. 66.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

Download references


The authors gratefully acknowledge and thank the KwaZulu-Natal Sharks Board, Oceans Research, South African Department of Agriculture, Forestry and Fisheries (DAFF), South African Shark Conservancy and White Shark Africa for providing valuable samples. The authors also thank Michelle-Angelique Hallendorff for editorial revision and Clint Rhode and Luca Mirimin for valuable discussions. This work was funded by the National Research Foundation, South Africa.

Author information



Corresponding author

Correspondence to Aletta E Bester-van der Merwe.

Additional information

Competing interests

The authors have no competing interests to declare.

Authors’ contributions

SNM performed cross-species amplification, microsatellite genotyping, genetic data analyses, organised the samples and drafted the manuscript. CR participated in cross-species amplification and microsatellite genotyping. RR-W participated in experimental design and coordination and contributed to manuscript preparation. AEB-vdM conceived the study, provided funds, participated in its design and coordination and contributed to manuscript preparation. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1:The 16 elasmobranch species of southern Africa selected for cross-species amplification, including family, species, distribution and sampling locations.(DOC 254 KB)

The 35 putative microsatellite markers developed from five closely related species for cross-species amplification in the study taxa, including the primers sequence, microsatellite repeat motif, annealing temperature (T

Additional file 2: A ) and GenBank accession numbers.(DOC 80 KB)


Additional file 3: ND2 sequence information of the study taxa used to estimate the genetic distance to evaluate cross-species performance, including ID Verified, availability of images (yes or no) which are available in the on-line host specimen database ( and GenBank accession numbers. (DOC 46 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Maduna, S.N., Rossouw, C., Roodt-Wilding, R. et al. Microsatellite cross-species amplification and utility in southern African elasmobranchs: A valuable resource for fisheries management and conservation. BMC Res Notes 7, 352 (2014).

Download citation


  • Cross-species amplification
  • Microsatellites
  • Multiplex assays
  • Genetic diversity
  • Species identification
  • Conservation management