Figures
Abstract
The population dynamics of shark species are generally poorly described because highly mobile marine life is challenging to investigate. Here we investigate the genetic population structure of the blacktip reef shark (Carcharhinus melanopterus) in French Polynesia. Five demes were sampled from five islands with different inter-island distances (50–1500 km). Whether dispersal occurs between islands frequently enough to prevent moderate genetic structure is unknown. We used 11 microsatellites loci from 165 individuals and a strong genetic structure was found among demes with both F-statistics and Bayesian approaches. This differentiation is correlated with the geographic distance between islands. It is likely that the genetic structure seen is the result of all or some combination of the following: low gene flow, time since divergence, small effective population sizes, and the standard issues with the extent to which mutation models actually fit reality. We suggest low levels of gene flow as at least a partial explanation of the level of genetic structure seen among the sampled blacktip demes. This explanation is consistent with the ecological traits of blacktip reef sharks, and that the suitable habitat for blacktips in French Polynesia is highly fragmented. Evidence for spatial genetic structure of the blacktip demes we studied highlights that similar species may have populations with as yet undetected or underestimated structure. Shark biology and the market for their fins make them highly vulnerable and many species are in rapid decline. Our results add weight to the case that total bans on shark fishing are a better conservation approach for sharks than marine protected area networks.
Citation: Vignaud T, Clua E, Mourier J, Maynard J, Planes S (2013) Microsatellite Analyses of Blacktip Reef Sharks (Carcharhinus melanopterus) in a Fragmented Environment Show Structured Clusters. PLoS ONE 8(4): e61067. https://doi.org/10.1371/journal.pone.0061067
Editor: Paul Sunnucks, Monash University, Australia
Received: July 21, 2011; Accepted: March 8, 2013; Published: April 9, 2013
Copyright: © 2013 Vignaud et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The DIREN Polynesia provided funding for the fieldwork in Maria. The Service de la Pêche (Polynesia) provided funding for the fieldwork in Fakahina. The Ministère de l'Ecologie du Développement Durable et de l'Energie, The FondsPacifique, CRISP and IMODEL provided funding or assisted with fieldwork, logistics and analysis. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Sharks have been considered keystone species driving the evolution of marine ecosystems [1], [2]. It is a concern therefore that their biological characteristics (e.g., low fecundity, late maturity, long gestation periods) mean shark populations are highly vulnerable and sensitive to overfishing and habitat degradation. Declining populations of sharks and other top predators from marine ecosystems is likely to have cascading effects on trophic structure and ecosystem dynamics [3]. Shark populations around the world have already suffered drastic declines, ranging from 50 to more than a 90% reduction in stocks, depending on species and area, and resulting in 34% of the oceanic species facing extinction [4]. Estimates indicate that even in recent years several dozen million sharks have been killed every year [5]. These alarming figures and media attention have raised awareness of the practice of shark finning and protecting sharks has become a worldwide conservation priority.
Effective conservation planning for any organism can partly rely on defining the spatial extent of populations, the levels of exchange (connectivity) among them and their distribution in space [6], [7]. Sharks are usually described as highly mobile species that often have large home ranges. As a consequence conventional marine protected areas (MPAs) and MPA networks are not likely to be effective in protecting shark populations [8] since home ranges are often larger than typical MPA sizes. Correctly estimating the connectivity between demes or populations is difficult and creates issues for and raises the uncertainty of the population models used to plan MPA networks [9]. Shark ecological traits generally differ significantly from those of most marine teleosts. For instance, sharks do not have a larval phase and for most species dispersal occurs only during adulthood [10]. Also, effective population sizes are much smaller for sharks relative to those of small reef fish species. These characteristics have important implications for the dynamics and evolution of shark populations.
Several studies that have found significant genetic heterogeneity among populations of shark species have included samples between oceans (e.g. [11]–[13]) leading to the suggestion of genetic homogeneity within ocean basins except in rare occasions (e.g. [14], [15]). When genetic heterogeneity has been found on smaller scales for some shark populations, evidence for explanations like a bottleneck and/or recent colonization has been presented, or the subject species is slow-moving (e.g. [14]–[22]). Here we assess the level of genetic structure among demes of blacktip reef sharks in the island chains of French Polynesia in the south-central Pacific Ocean. Given the literature, it is hard to hypothesize about blacktip reef shark dispersal and gene exchange. They are small coastal sharks commonly residing within coastal habitats that can be separated from other suitable habitat by hundreds of km, which could suggest limited dispersal. However, they have colonized remote islands, are fast active-swimming species, and a recent study has shown that blacktip reef shark females in French Polynesia will cross open water between islands as much as 50 km apart [23]. Whether dispersal occurs in the fragmented environment of French Polynesia frequently enough to prevent moderate genetic structure is unknown.
French Polynesia is a prime example of a fragmented marine environment for a shark species known to inhabit coastlines; there are 130 islands in French Polynesia spread across 2,500,000 km2 of ocean [24]. Coral reef fish communities and blacktip reef sharks are concentrated and aggregated around islands separated by large bodies of open deep water. The open water separating the islands may act as at least a partial barrier to dispersal; the waters between the islands can be >1000 m deep and some islands are separated by >100 km. Blacktip reef sharks are viviparous and produce a few (2–5) pups generally once a year. The gestation period is about 10–11 months, and it is likely that sexual maturity is not reached before 4–5 years old [25], [26]. They are easy to catch in shallow waters either in the lagoon or outside the barrier reef. Shark populations have been over-fished in many coastal areas [6], [27]–[30] but fishing pressure on sharks is still relatively low in French Polynesia. Human population numbers are low, the local culture favors sharks and laws have recently been passed to protect sharks. Shark population dynamics are thus likely to be as close to natural as possible. Some sharks are killed in remote islands to prevent damage to fishing gear and traps but direct human pressure on blacktip reef sharks in French Polynesia is low.
We tested for the presence of genetic structure among 5 shark demes from the reef and coastal habitat surrounding 5 different islands with varying inter-island distances (50 to 1500 km). We use the word ‘deme’ here to define a group of individuals that have been sampled around the same island and because these may not be populations per se.
Methods
Study area and sampling
Locations of sampling sites in French Polynesia are detailed in Figure 1. A total of 165 individuals were sampled from 5 locations: Moorea (high island, N = 38), Tetiaroa (small semi-closed atoll, N = 12), Rangiroa (large open atoll, N = 35), Fakahina (small closed atoll, N = 41) and Maria (small closed atoll, N = 39). These locations ensure a variety of distances between sampling sites. In each location, sharks were fished with barbless circled hooks and metal lines from the shore or from a boat. The animals were sized, sexed, and a small tissue sample was taken from either a dorsal or anal fin. Each shark was then released in good condition. Fin clip samples were conserved in 90% Ethanol.
Distances (in km) between the sample locations are shown in the matrix top right.
Laboratory procedures
DNA samples were extracted following the PURAGENE™ DNA Purification Kit procedure for 5–10 mg fresh tissue. PCR was conducted using 4 multiplexes for 11 microsatellite loci, and following QIAGEN® Multiplexes PCR Kit procedure. The loci characteristics and thermal cycling conditions used are described in Table S1 and were taken and adapted from [31]–[33] and optimized for the target species and multiplex PCR. PCR products were analyzed using a capillary electrophoresis sequencing Beckman CEQ™8000, and the fragment sizes were manually read using the program CEQ™8000.
Data analyses
A first set of analyses was computed to test the quality of the 11 markers. Microcheckerv2.2.3 [34] was used to check potential genotyping errors and null allele(s). GenAlEx [35] was used to test for the Hardy-Weinberg equilibrium for each locus in the entire dataset, and for each locus in each deme. GenAlEx was also used to compute general loci information such as allelic richness and private alleles.
We used one genotypic clustering approach and two allele-based statistics approaches (FST and RST) [36]–[39] to measure the genetic structure between the 5 demes. The clustering approach used is the one implemented in Adegenet [40] for R [41], based on a Bayesian algorithm for the Discriminant Analysis of Principal Component (DAPC) [42]. The optimal number of principal components chosen for the analysis was of 21, following indications of alpha scores [42], [43].
Pairwise FST and RST values were calculated using the AMOVA method implemented in GenAlEx, in order to measure the genetic differences between demes. These values were also calculated for each locus to help understand potential discrepancies between the global pairwise FST and RST values. Additionally, an AMOVA was performed using GenAlEx to determine the global level of differentiation among clusters.
We used ARES [44] to estimate the allelic richness for each deme and to extrapolate it to the higher number of individuals among samples, for comparison among demes (Table S2). We used the program Bottleneck [45] to check for bottlenecks and changes in effective population size. To inform this analysis we checked heterozygote excesses, using the 3 different models: the infinite alleles model (IAM), the stepwise-mutation model (SMM) and the two-phase microsatellite evolution model (TPM) with 70% SMM and 30% IAM.
Because of the low number of demes, the correlation between geographic distance and genetic differences was examined by jackknifing over populations, using Isolation By Distance [46]. This approach calculates the slope and intercept of the relationship between genetic distances and geographic distances using a Reduced Major Axis regression. We used the log of FST values and the log of geographic distances, as suggested in [47] when using distances in two dimensions.
Considering the stepwise mutation model, we used the program SPAGeDI [48] to test if stepwise-like mutations contribute to genetic differentiation, using the RST/pRST test described in [49]. Loci with 3 or less alleles (Cli103, Cli2 and LS20) were excluded from the analysis. Finally, a sex-bias dispersal test was performed using GenAlEx following the method described in [49], [50], and deleting all juveniles and small sub-adults (N = 35) from the dataset.
Results
Heterozygosity ranged from 0.077 (Cli2) to 0.865 (LS53) with a mean of 0.492 (Table 1). Overall, none of the loci showed significant departures from the Hardy-Weinberg equilibrium (all p-values>0.05). However, when looking at the patterns from each deme, 3 loci diverged from the Hardy-Weinberg equilibrium in the deme from Moorea (Cli107, Cli111, LS54) and 2 in the deme from Rangiroa (Cli108 and LS32). Cli2 is monomorphic in the deme from Fakahina. Within all samples (165 individuals) and the 11 microsatellites we found 79 alleles including 14 private alleles distributed in 10 of the 11 loci. Allele frequencies for each deme are ordered by size in Table S3, which has a histogram describing differences in allele frequency between demes. Differences in allele frequencies are used here to estimate genetic differences between demes. These raw frequencies, for each deme and each locus, are very different in some of our cases. As an example, the allele 220 in the locus LS20 is present 5.7% of the time in the Rangiroa deme but 42.4% of the time in the Fakahina deme.
Genetic structure
The DAPC results for the 21 principal components found show a clear genetic structure for all demes.
Overall, individuals are associated with their original deme. The membership probability is particularly strong for Fakahina as a majority of individuals are 100% associated with the Fakahina deme. A large number of unsampled demes are likely reinforcing the structure of isolated demes. However, individuals among the 4 other demes also show strong membership probability to their original deme, but with more individuals having mixed signals than it was the case with Fakahina. There is very little genetic signal from Fakahina in other demes. In contrast, Maria has mixed individuals with signal mostly from Moorea and Rangiroa. Some outliers are observed in all clusters, suggesting potential recent migrations. These results are depicted in Figure 2, including the sex and size of a few outliers. The outliers are mostly male around 115 cm in total length, roughly the size at which they become mature. The lowest Bayesian Information Criterion (BIC) values were between 4 and 6, although there is no clear elbow in BIC value curve.
The optimal number of principal components found for the analysis was 21 based on the trend in alpha scores (a). The BIC (Bayesian Information Criterion) values are shown in relation to the number of genetic clusters in (b). Each vertical bar represents an individual in the DAPC diagram (compoplot) shown as (c), and each color represents the probability of belonging to one of the genetic clusters. Some outliers have been noted on the top of the figure: M = male, F = female, followed by the total length of each shark in cm.
Pairwise FST results show high values for the Fakahina deme (vs. each of the other demes; FST = 0.097, 0.098, 0.142 and 0.148) and for the Maria deme (vs. each of the other demes; FST = 0.080, 0.095, 0.148 and 0.106, Table 2). The FST values associated with Fakahina are higher than those associated with Maria, even though Maria is further from the other islands. The Moorea, Tetiaroa and Rangiroa demes are less differentiated among each other (FST = 0.038, 0.025 and 0.043). FST results are highly significant for all deme pairs (p = 0.01). Fakahina RST values are of 0.301 versus Moorea, 0.437 versus Tetiaroa, 0.111 versus Rangiroa and 0.081 versus Maria (Table 2). Maria RST values, except for versus Fakahina, are non-significant and small (0.009, 0, 0.003). Rangiroa RST values versus Moorea and Tetiaroa are small (0.0041 and 0.126) but significant (p = 0.017 and 0.008). The Moorea versus Tetiaroa RST value is 0.054, and is non-significant (p = 0.063).
AMOVA for 3 geographic clusters (Fakahina vs. Maria vs. others) indicated that 8% of the variation was caused by variation among geographic clusters, 38% by the variation within individuals and 54% by the variation among individuals. When an AMOVA was performed for either 4 or 5 groups, the ‘among populations’ value was 7%. The Isolation by Distance shows that genetic and geographic differences are correlated (R2 = 0.739); the slope after jackknifing over the 5 demes (0.3654) has a large standard error (0.275, Figure 3).
Jackknifing slope over the 5 demes is 0.3654±0.275 (SE). The data points denote two-location comparisons: Mo = Moorea, R = Rangiroa, F = Fakahina, T = Tetiaroa and Ma = Maria (R2 = 0.739).
General population information
As a proxy of effective population size, we estimated values of allelic richness for each deme. To compare these values we had to restrict (Fakahina)/extrapolate (Moorea, Tetiaroa, Rangiroa, Maria) the number of samples used for the analysis to 40. The estimated values of allelic richness ranged from 68 for Rangiroa to 44 for Fakahina and Maria, 57 for Moorea and 53 for Tetiaroa. Despite these differences in allelic richness, no evidence of a bottleneck was found, using the IAM, SMM and TPM (all>0.05). The microsatellite allele frequency distribution was L-shaped for each locality, indicating that these were all at mutation-drift equilibrium for the loci used. The RST/pRST results on the 8 loci with more than 3 alleles are significant for Moorea-Fakahina (P(one sided test, H1: obs<exp) = 0.984) and for Tetiaroa-Fakahina (P(one sided test, H1: obs<exp) = 0.977) (Table 3). There are large differences though when looking at the results for each locus, and the high variance in the results for Tetiaroa may be due to a low sample size.
Discussion
In combination, the FST and RST estimators can provide valuable information about the genetic structure of populations with some provisos [51]. These measures have different qualities and work under different assumptions, none of which are fully known or constant in all cases. One of the reasons for this involves the highly versatile behavior of microsatellites. Microsatellites are not strictly in the Infinite Allele Model (IAM) or the Stepwise Mutation Model (SMM), and many have varying mutation rates. Other reasons to carefully interpret FST and RST results include the presence of evolutionary forces like genetic drift [52] or ghost populations [53].
Even given the above, in the case presented here pairwise FST values are high, suggesting strong genetic differences among demes. This is especially true for the Fakahina and Maria demes, which are further from Moorea, Tetiaroa and Rangiroa than any of these are from each other. The Isolation by Distance analysis shows that genetic differentiation (observed from FST values) mostly increases with the distance separating demes. With only 5 sampled demes, it is hard to make general conclusions based on the Isolation by Distance results. Pairwise FST values are slightly higher though for Fakahina than for Maria, even though Maria is furthest the other islands sampled. This could be due to differences between these islands in the degree of isolation. The island closest to Fakahina is more than 200 km away while 4 islands in the Tuamotu chain are within 100 km of Maria. The proximity of islands in the Tuamotu chain to Maria is likely promoting some gene flow. The effective population size estimates of Fakahina and Maria are exactly the same based on allelic richness and so probably do not affect the comparison of their FST values.
FST values are known to underestimate genetic differentiation in the case of strong genetic structure [53]. It is thus possible to tentatively conclude that there is poor gene flow between Fakahina and Maria, and between the other sampled demes and Fakahina or Maria, and a moderated gene flow in the cluster Moorea-Tetiaroa-Rangiroa. The Bayesian approach results visually strongly suggest genetic differentiation, adding further support to this conclusion (see Figure 2). These results probably indicate a limited gene flow between the investigated blacktip shark demes. The Bayesian approach can result in the detection of a number of “populations”, usually written as ‘K’, [54]. We do not present a result for K here because it does not capture the more complex ecological reality of our case. Isolation by Distance and the presence of numerous unsampled demes may affect clustering approaches. These effects can create groups that do not necessarily exist if the sampled groups represent only a few points along a potentially graded genetic continuum. Even so, high levels of gene flow rapidly overwhelm the influence of other factors, leading to less structured demes than we see here. This does not mean that any of the demes, even Fakahina and Maria, are completely isolated as they are probably mixed with closer unsampled demes. The level of genetic structure seen in the study demes may also be influenced by small population sizes associated with genetic drift.
The Fakahina deme appears more differentiated when using the RST, potentially indicating that the gene flow is not sufficient to counter the effect of a higher level of genetic divergence. The only significant RST/pRST results are associated with Fakahina. Variance is high for the results among loci but this result is consistent with the particular geographic isolation of Fakahina. Migration rates are unlikely to be sufficient at Fakahina to keep allele sizes from diverging from the other demes. When these new microsatellite alleles appear in a deme, it means that gene flow has been especially low for many generations.
The patterns are not consistent though – there are never more than 4 loci significant for the same pair of demes. As with the FST results this suggests structuration by genetic drift with some level of gene flow. Surprisingly, RST values associated with Maria (except for Maria-Fakahina) are small and non-significant. The small RST values associated with Maria deme probably come from loci that do not fit the SMM behind the RST calculation. We suggest the RST values between demes are misleading in comparison to the results of FST value comparisons. Our case for this is that the FST values, the results of the Bayesian approach and the allele distribution (see Table S3) all suggest that Maria deme is strongly differentiated from any other deme.
Tests for sex-biased dispersal seem to hint at a pattern of male-based dispersal (non-significant result) but the sample sizes are likely too low to detect a pattern. Most of the outliers found in the clustering analysis are indeed males, supporting the suggestion that male-based dispersal is a life-history trait of these populations. This is aligned with the findings of other similar studies, like that of Carcharhinus limbatus in [55], and reviewed in [56]. Females also have the capacity to move from one island to another and have been shown to cross open waters in French Polynesia to give birth in the waters near islands where they may have been born [23]. Consistently giving birth in environments known to be suitable and safe could be part of a powerful evolutionary process. This could explain the lower detection of female dispersal in this study as the range of females ensures they can come back to a nursery area regularly (every year in some case) to give birth. It is unlikely though that philopatric behavior would always occur, or would have prevented the species from colonizing new areas.
Understanding the complex process driving genetic structure can be tricky, and results must be treated with caution (reviews and examples in [52], [57]–[59]). It is likely that the genetic structure seen for the blacktip demes is the result of all or some combination of the following: low gene flow, time since divergence, small effective population sizes, and the standard issues with the extent to which mutation models actually fit reality. We suggest low levels of gene flow as at least a partial explanation of the level of genetic structure seen among the sampled blacktip demes. This explanation is consistent with the ecological traits of blacktip reef sharks, and that the suitable habitat in French Polynesia is highly fragmented.
Evidence for genetic structuration of the blacktip demes we studied highlights that similar species may have populations with as yet undetected or underestimated structure. The diversity of genes tested and techniques used in the literature make the possible comparisons with other shark studies irrelevant. It is noteworthy though that genetic subdivision has been found when investigating active swimming coastal shark species. These include Carcharhinus sorrah in the Indo-Australian archipelago [14], and Carcharhinus limbatus in the northwestern Atlantic Ocean [15]. Once evidence for structure is found, logical follow-on research involves helping understand the causes. Sequencing more genes and analyzing the overall genetics with coalescent models represents planned future work with blacktips in French Polynesia. Such research may increase our understanding of the relative influence of levels of gene exchange, genetic drift and divergence time on the levels of genetic structure we document here.
Finding moderate levels of genetic structure in the blacktip reef shark demes we studied has implications for management and conservation planning. Both the adults (as reproducers and potential dispersers) and juveniles/sub-adults (as future adults) need to be protected to ensure shark populations are self-sustaining. Range sizes for sharks are difficult to assess but are highly likely to be larger than the proposed areas management agencies and conservationists plan to have in a marine protected area (MPA) network. Further, shark populations are naturally small, at least in fragmented environments. Protecting large areas will not result in the population and biomass increases seen in teleost fish protected in MPAs. For these reasons it is important to consider shark populations in conservation planning but MPAs are not likely to be the best solution for protecting sharks. The management ‘best practice’ approach to protecting most shark populations is a total ban of shark fishing everywhere this is possible. The Bahamas, French Polynesia, and Palau are examples of locations where such a ban has already been implemented. These nation-wide regulations with other recently implemented shark sanctuaries represent steps in the right direction. Other management actions like MPAs remain critical and can have positive effects on shark populations but are insufficient to avoid the severe declines caused by shark fishing.
Supporting Information
Table S1.
Locus description with sequence, repeat type and hybridization temperatures. All ‘Cli’ are taken from [31]. LS20, LS32 and LS53 are taken from [32]. LS54 and LS75 are taken from [33].
https://doi.org/10.1371/journal.pone.0061067.s001
(TIF)
Table S2.
Extrapoled allelic richness and means. Tetiaroa was excluded because of a low number of samples (12).
https://doi.org/10.1371/journal.pone.0061067.s002
(TIF)
Table S3.
Allele frequencies for each deme ordered by size (left) and presented as histograms.
https://doi.org/10.1371/journal.pone.0061067.s003
(TIF)
Acknowledgments
We thank the many local community members that provided field assistance at the study sites, and the NO Braveheart crew for assisting in the Actéons. Pablo Saenz and Cecile Fauvelot assisted with genetic analysis and Elisabeth Rochel assisted with microsatellite screening. Versions of the ms benefited from comments from Pablo Saenz-agudelo, Jonathan Werry, Judith Bakker, the editor Paul Sunnucks and anonymous reviewers.
Author Contributions
Conceived and designed the experiments: TV SP. Performed the experiments: TV J. Mourier EC SP. Analyzed the data: TV. Contributed reagents/materials/analysis tools: TV SP. Wrote the paper: TV J. Maynard.
References
- 1. Stevens JD, Bonfil R, Dulvy NK, Walker PA (2000) The effects of fishing on sharks, rays, and chimaeras (chondrichthyans), and the implications for marine ecosystems. ICES J Mar Sci 57: 476–494.
- 2. Bascompte J, Melian CJ, Sala E (2005) Interaction strength combinations and the overfishing of a marine food web. Proc Natl Acad Sci USA 102: 5443–5447.
- 3. Myers RA, Baum JK, Shepherd T, Powers SP, Peterson CH (2007) Cascading effects of the loss of apex predatory sharks from a coastal ocean. Science 315: 1846.
- 4.
Camhi MD, Valenti SV, Fordham SV, Fowler SL, Gibson C (2009) The conservation status of pelagic sharks and rays. Report of the IUCN Shark Specialist Group. Pelagic Shark Red List Workshop. IUCN Species Survival Commission Shark Specialist Group. Newbury, UK. x + 78p
- 5. Clarke S, McAllister MK, Milner-Gulland EJ, Kirkwood GP, Michielsens CGJ, et al. (2006) Global estimates of shark catches using trade records from commercial markets. Ecology Letters 9: 1115–1126.
- 6. McCook LJ, Almany GR, Berumen ML, Day JC, Green A, et al. (2009) The challenge of incorporating connectivity science into coral reef management now: principles and practice. Coral Reef 28: 353–366.
- 7. Botsford LW, White JW, Coffroth MA, Paris CB, Planes S, et al. (2009) Connectivity and resilience of coral reef metapopulations in MPAs: matching empirical efforts to predictive needs. Coral Reef 28: 327–337.
- 8. Robbins WD, Hisano M, Connolly SR, Choat JH (2006) Ongoing Collapse of Coral-Reef Shark Populations. Current Biology 16: 2314–2319.
- 9. Planes S, Jones G, Thorrold S (2009) Larval dispersal connects fish populations in a network of marine protected areas. Proceedings of the National Academy of Sciences 106 (14) 5693–5697.
- 10.
Musick JA, Harbin MM, Compagno LJV (2004) Historical zoogeography of the Selachii. In: Biology of Sharks and Their Relatives (Carrier JC, Musick JA, Heithaus MR). CRC press, Boca Raton, FL, 33–78
- 11. Schultz JK, Feldheim KA, Gruber SH, Ashley MV, McGovern TM, et al. (2008) Global phylogeography and seascape genetics of the lemon sharks (genus Negaprion). Mol Ecol 17 (24) 5336–48.
- 12. Benavides MT, Horn RL, Feldheim KA, Shivji MS, Clarke SC, et al. (2011) Global phylogeography of the dusky shark Carcharhinus obscurus: implications for fisheries management and monitoring the shark fin trade. Endangered Species Research 14 (1) 13–22.
- 13. Daly-engel TS, Seraphin KD, Holland KN, Coffey JP, Nance HA, et al. (2012) Global Phylogeography with Mixed-Marker Analysis Reveals Male-Mediated Dispersal in the Endangered Scalloped Hammerhead Shark (Sphyrna lewini). PLoS ONE 7 (1) e29986.
- 14. Ovenden JR, Kashiwagi T, Broderick D, Giles J, Salini J (2009) The extent of population genetic subdivision differs among four co-distributed shark species in the Indo-Australian archipelago. BMC evolutionary biology 9: 40.
- 15. Keeney DB, Heupel MR, Hueter RE, Heist EJ (2005) Microsatellite and mitochondrial DNA analyses of the genetic structure of blacktip shark (Carcharhinus limbatus) nurseries in the northwestern Atlantic, Gulf of Mexico, and Caribbean Sea. Mol Ecol 14: 1911–1923.
- 16. Heist EJ, Musick JA, Graves JE (1996) Mitochondrial DNA Diversity and Divergence Among Sharpnose Sharks, Rhizoprionodon terraenovae, From the Gulf of Mexico and Mid-Atlantic Bight. Fish Bull 94: 664–668.
- 17. Ahonen H, Harcourt RG, Stow A (2009) Nuclear and mitochondrial DNA reveals isolation of imperilled grey nurse shark populations (Carcharius taurus). Mol Ecol 18: 4409–4421.
- 18. Zenger KR, Stow AJ, Peddemors V, Briscoe DA, Harcourt RG (2006) Wide-spread utility of highly informative AFLP molecular markers across divergent shark species. Journal of Heredity 97: 607–611.
- 19. Feldheim KA, Stow AJ, Ahonen H, Chapman DD, Shivji M, et al. (2007) Polymorphic microsatellite markers for studies of the conservation and reproductive genetics of imperilled sand tiger sharks (Carcharias taurus). Mol Ecol Resour 7: 1366–1368.
- 20. Corrigan S, Huveneers C, Schwartz TS, Harcourt RG, Beheregaray LB (2008) Genetic and reproductive evidence for two species of ornate wobbegong shark on the Australian East Coast. J Fish Biol 73: 1662–1675.
- 21. Lewallen EA, Anderson TW, Bohonak AJ (2007) Genetic structure of leopard shark (Triakissemifasciata) popula- tions in California waters. Mar Biol 152: 599–60.
- 22. Gaida IH (1997) Population structure of the Pacific angel shark, Squatinacalifornica (Squatiniformes, Squatinidae), around the California Channel Islands. Copeia 1997: 738–744.
- 23. Mourier J, Planes S (2013) Direct genetic evidence for reproductive philopatry and associated fine-scale migrations in female blacktip reef sharks (Carcharhinus melanopterus) in French Polynesia. Mol Ecol 22 (1) 201–14.
- 24.
Kingfisher Geography Encyclopedia. ISBN 1-85613-582-9. Page 546
- 25.
Compagno LJV (1984) FAO Species Catalogue. Vol 4: Sharks of the World. An annotated and illustrated catalogue of shark species known to date. Part 2: Carcharhiniformes. Food and Agriculture Organization of the United Nations, Rome.
- 26. Papastamatiou YP, Lowe CG, Caselle JE, Friedlander AM (2009) Scale-dependent effects of habitat on movements and path structure of reef sharks at a predator-dominated atoll. Ecology 90(4): 996–1008.
- 27. Baum JK, Myers RA, Kehler DG, Worm B, Harley SJ, et al. (2003) Collapse and Conservation of Shark Populations in the Northwest Atlantic. Science 299 (5605) 389.
- 28.
Bonfil R (1994) Overview of world elasmobranch fisheries. FAO Fisheries Technical Paper, 341, 119pp
- 29.
Watts S, Wu V (2005) At rock bottom: the declining sharks of the eastern tropical pacific. WildAid report.
- 30. Heupel MR, Williams AJ, Welch DJ, Ballagh A, Mapstone BD, et al. (2009) Effects of fishing on tropical reef associated shark populations on the Great Barrier Reef. Fisheries Research 95: 350–361.
- 31. Keeney DB, Heist EJ (2003) Characterization of microsatellite loci isolated from the blacktip shark and their utility in requiem and hammerhead sharks. Mol Ecol Notes 3: 501–504.
- 32. Feldheim KA, Ashley MV, Gruber SH (2001) Population genetic structure of the lemon shark (Negaprion brevirostris) in the western Atlantic: DNA microsatellite variation. Mol Ecol 10: 295–303.
- 33. Feldheim KA, Gruber SH, Ashley MV (2002) The Breeding Biology of lemon sharks at a Tropical Nursery Lagoon. Proceedings of the Royal Society (London) B 269 (1501) 1655–1661.
- 34. van Oosterhoot C, Hutchinson W, Wills D, Shipley P (2004) MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538.
- 35. Peakall R, Smouse PE (2006) GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6: 288–295.
- 36. Weir BS, Cockerham CC (1984) Estimating F-stastistics for the analysis of population structure. Evolution 38: 1358–1370.
- 37. Michalakis Y, Excoffier L (1996) A genetic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics 142: 1061–1064.
- 38. Waples RS, Gaggiotti O (2006) What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol Ecol 15: 1419–1439.
- 39. Guillot G, Leblois R, Coulon A, Frantz A (2009) Statistical methods in spatial genetics. Mol Ecol 18: 4734–4756.
- 40. Jombart T (2008) adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24: 1403–1405.
- 41.
R Development Core Team (2011) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0
- 42. Jombart T, Devillard S, Balloux F (2010) Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genetics 11: 94.
- 43.
Jombart T (2011) A tutorial for Discriminant Analysis of Principal Components (DAPC) using adegenet 1.3-0
- 44. van Loon E, Cleary DFR, Fauvelot C (2007) ARES: software to compare allelic richness between uneven samples. Mol Ecol Notes 7: 579–582.
- 45. Cornuet JM, Luikart G (1997) Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 144: 2001–2014.
- 46. Jensen JL, Bohonak AJ, Kelley ST (2005) Isolation by distance, web service. BMC Genetics 6: 13.v.3.16.
- 47. Slatkin M (1993) Isolation by distance in equilibrium and non-equilibrium populations. Evolution 47: 264–279.
- 48. Hardy OJ, Vekemans X (2002) SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes 2: 618–620.
- 49. Hardy OJ, Charbonnel N, Freville H, Heuertz M (2003) Microsatellite allele sizes: a simple test to assess their significance on genetic differentiation. Genetics 163 (4) 1467–1482.
- 50. Mossman CA, Waser PM (1999) Genetic detection of sex-biased dispersal. Mol Ecol 8: 1063–1067.
- 51. Holsinger KE, Weir BS (2009) Genetics in geographically structured populations: defining, estimating and interpreting FST. Nature Reviews Genetics 10: 639–650.
- 52. Marko PB, Hart MW (2011) The complex analytical landscape of gene flow inference. Trends Ecol Evol 26 (9) 448–456.
- 53. Slatkin M (2005) Seeing ghosts: the effect of unsampled populations on migration rates estimated for sampled populations. Mol Ecol 14: 67–73.
- 54. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620.
- 55. Keeney DB, Heist EJ (2006) Worldwide phylogeography of the blacktip shark (Carcharhinus limbatus) inferred from mitochondrial DNA reveals isolation of western Atlantic populations coupled with recent Pacific dispersal. Mol Ecol 15: 3669–3679.
- 56. Hueter RE, Heupel MR, Heist EJ, Keeney DB (2005) Evidence of philopatry in sharks and implications for the management of shark fisheries. Journal of Northwest Atlantic Fishery Science 35: 239–247.
- 57. Lowe WH, Allendorf FW (2010) What can genetics tell us about population connectivity? Mol Ecol 19: 3038–3051.
- 58. Sunnucks P (2011) Towards modelling persistence of woodland birds: the role of genetics. Emu 111: 19–39.
- 59. Balloux F, Lugon-Moulin N (2002) The estimation of population differentiation with microsatellite markers. Mol Ecol 11: 155–165.