Structure and Genetic Variability of the Oceanic Whitetip Shark, Carcharhinus longimanus, Determined Using Mitochondrial DNA

Abstract

Information regarding population structure and genetic connectivity is an important contribution when establishing conservation strategies to manage threatened species. The oceanic whitetip shark, Carcharhinus longimanus, is a highly migratory, large-bodied, pelagic shark listed by the IUCN (International Union for Conservation of Nature) Red List as "vulnerable" throughout its range and “critically endangered” in the western north Atlantic. In 2014, the species was protected globally under Appendix II of CITES (Convention on International Trade in Endangered Species), limiting and regulating trade. This study used partial sequences of mitochondrial DNA (mtDNA) control region to determine the population genetic structure of oceanic whitetip sharks across the Atlantic and Indian Oceans. 724 base pairs were obtained from 215 individuals that identifed nine polymorphic sites and defined 12 distinct haplotypes. Total nucleotide diversity (π) was 0.0013 and haplotype diversity (h) was 0.5953. The Analysis of Molecular Variance (AMOVA) evidenced moderate levels of population structure (ɸ_ST = 0.1039) with restricted gene flow between the western and eastern Atlantic Ocean, and a strong relationship between the latter region and the Indian Ocean. Even though the oceanic whitetip is a highly migratory animal the results presented here show that their genetic variability is slightly below average of other pelagic sharks. Additionally, this study recommends that at least two populations in the Atlantic Ocean should be considered distinct (eastern and western Atlantic) and conservation efforts should be focused in areas with the greatest genetic diversity by environmental managers.

About the Authors:

Sâmia M. Camargo

Affiliations Laboratório de Biologia e Genética de Peixes, Instituto de Biociências de Botucatu, Universidade Estadual Paulista, (UNESP), Botucatu, São Paulo, Brasil, Laboratório de Genética Pesqueira e Conservação, Instituto do Mar, Universidade Federal de São Paulo (UNIFESP), Santos, São Paulo, Brasil

Rui Coelho

Affiliations Instituto Português do Mar e da Atmosfera (IPMA), Olhão, Portugal, Centro de Ciências do Mar (CCMAR), Universidade Algarve, Faro, Portugal

Demian Chapman

Affiliation: School of Marine and Atmospheric Science, Stony Brook University, Stony Brook, New York, United States of America

Lucy Howey-Jordan

Affiliation: Microwave Telemetry, Inc., Columbia, Maryland, United States of America

Edward J. Brooks

Affiliation: Shark Research and Conservation Program, Cape Eleuthera Institute, Eleuthera, The Bahamas

Daniel Fernando

Affiliations The Manta Trust, Catemwood House, Corscombe, Dorchester, United Kingdom, Department of Biology and Environmental Science, Linnaeus University, Lund, Sweden, Blue Resources, Colombo, Sri Lanka

Natalia J. Mendes

Affiliation: Laboratório de Biologia e Genética de Peixes, Instituto de Biociências de Botucatu, Universidade Estadual Paulista, (UNESP), Botucatu, São Paulo, Brasil

Fabio H. V. Hazin

Affiliation: Departamento de Pesca e Aquicultura, Universidade Federal Rural de Pernambuco, UFRPE, Recife, Pernambuco, Brasil

Claudio Oliveira

Affiliation: Laboratório de Biologia e Genética de Peixes, Instituto de Biociências de Botucatu, Universidade Estadual Paulista, (UNESP), Botucatu, São Paulo, Brasil

Miguel N. Santos

Affiliation: Instituto Português do Mar e da Atmosfera (IPMA), Olhão, Portugal

Fausto Foresti

Affiliation: Laboratório de Biologia e Genética de Peixes, Instituto de Biociências de Botucatu, Universidade Estadual Paulista, (UNESP), Botucatu, São Paulo, Brasil

Fernando F. Mendonça

* E-mail: [email protected]

Affiliation: Laboratório de Genética Pesqueira e Conservação, Instituto do Mar, Universidade Federal de São Paulo (UNIFESP), Santos, São Paulo, Brasil

Introduction

The effects of unsustainable fishing on populations of sharks and rays (elasmobranchs) has been well-documented globally [1], and studies have shown that over the last several decades many large migratory species commonly caught in large scale pelagic marine fisheries are rapidly declining [2, 3]. The International Union for Conservation of Nature (IUCN), which quantifies the conservation status of most taxa, currently lists over a thousand species of elasmobranchs [4, 5]. The highly migratory sharks are amongst the species with highest conservation concerns [6].

The oceanic whitetip shark, Carcharhinus longimanus, is a pelagic shark with worldwide distribution, found mainly in tropical and subtropical open waters [7, 8]. Historically grouped with the silky shark C. falciformis and the blue shark Prionace glauca, C. longimanus was once considered a highly abundant oceanic shark [7–11]. However, recent population assessments and anecdotal reports suggest that the oceanic whitetip shark is now only occasionally recorded in the Atlantic. A recent study by Coelho et al. [12] covering a wide area of the Atlantic in both hemispheres indicated that the oceanic whitetip shark bycatch in pelagic longline fisheries is less than 1% of the total elasmobranch catches. Furthermore, philopatric behavior, defined as the tendency of individuals to return or stay in their home areas, natal (birth) sites or adopted locales[13], was recently described in this species, highlighting the need for designated protective areas for oceanic whitetip sharks.

The IUCN Red List currently lists C. longimanus as "vulnerable" throughout its range and “critically endangered” in the western north Atlantic [14], and in 2014 was listed in Appendix II of CITES (Convention on International Trade in Endangered Species), regulating international trade and requirement of appropriate permits, including documentation validating that catch is legally sourced. A 2012 Pacific-wide stock assessment for oceanic whitetip sharks showed that the populations are currently overfished and stocks have declined to levels below maximum sustainable yield [15]. In the Atlantic and Indian Oceans, Ecological Risk Assessments for pelagic sharks, show that the oceanic whitetip shark has a relatively high vulnerability to pelagic fisheries, due to its low fecundity (average three pups per litter] and high susceptibility [16, 17]. Retaining oceanic whitetip sharks for commercial purposes has been prohibited in tuna RFMOs (Regional Fisheries Management Organizations) of all Oceans, including the ICCAT (International Commission for the Conservation of Atlantic Tunas) Atlantic convention area since 2010 (ICCAT Recommendation 10–07), the IOTC (Indian Ocean Tuna Commission) Indian Ocean convention area since 2013 (IOTC Resolution 13/06), the IATTC (Inter-American Tropical Tuna Commission) in the eastern Pacific since 2012 (Resolution C-11-10), and the WCPFC (Western and Central Pacific Fisheries Commission) in the western and central Pacific since 2013 (CMM 2011–04).

An aspect that is fundamental to sustainable use of marine resources and quantification of productivity of ocean ecosystems is the identification and maintenance of distinct stocks [18]. The conservation of the genetic variability is one of the basic objectives in regulatory programs assisting in the recovery of endangered species [19]. As such, the sequencing of the mitochondrial DNA control region (D-loop) is one of the most commonly applied method to population genetic studies of vertebrates, including sharks [20–26].

Population declines, such as that observed in the oceanic whitetip shark over the last 50 years, raises questions about the maintenance of genetic variability and the putative loss of evolutionary lineages. These issues are more relevant in species with population structure where the distribution of genetic characteristics is geographically limited. Therefore, this study aimed to describe the genetic variability and population structure of C. longimanus from the Atlantic and Indian Oceans by using control region of mitochondrial DNA with an intention to provide a framework for future assessments and tracking studies.

Materials and Methods

Sampling

Samples of C. longimanus were collected by fishery observers on Portuguese, Brazilian and American commercial pelagic longline vessels in several regions of the Atlantic and Indian Oceans. Tissue samples, including small muscle or fin fragments (<1cm³), were collected opportunistically during the normal fishing operations of the vessels, and were subsequently frozen in 95% ethanol. Tissue samples from the north Atlantic were collected in The Bahamas as part of an ongoing study of oceanic whitetip movement and life history. All samples were collected under the Cape Eleuthera Institute research permit (MAF/FIS/17 and MAF/FIS/34) and in accordance with Stony Brook University and Cape Eleuthera Institute regulations developed within the guidelines of the Association for the Study of Animal Behaviour and the Animal Behavior Society [27]. Sampling from the Portuguese fishery was conducted by the Portuguese Institute for the Ocean and Atmosphere (IPMA), within the scope of the European Data Collection Framework (PNAB/DCF). This sampling was conducted under an ICCAT permit for biological sampling of shark species in the Atlantic (ICCAT Recommendation 13/10) and an IOTC permit for biological sampling of oceanic whitetip sharks in the Indian Ocean (IOTC Resolution 13/06). No specific permissions were required for internationally transporting samples of C. longimanus before the species effectively became part of CITES in September 2014, all samples used in this study were collected before that date.

A total of 215 specimens were sampled for this study. Specifically, 206 were sampled from various regions of the Atlantic Ocean, including 28 from the Northwest Atlantic, 51 from the Western Equatorial Atlantic, 54 from the Northeast Tropical Atlantic, 50 from the Meso-Atlantic, 17 from the Southeast Atlantic and six from the Southwest Atlantic. In addition, nine samples were collected from the Indian Ocean (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Geographic distribution of samples of Carcharhinus longimanus with the network haplotypes analyzed and compiled from the sequences of the mitochondrial DNA control region.

https://doi.org/10.1371/journal.pone.0155623.g001

DNA extraction and sequencing

Genomic DNA was extracted using the NucleoSpin® Tissue XS Kit (Macherey & Nagel, Düren, Germany). Partial sequences of the control region of the mitochondrial DNA were attained according to Mendonça et al. [28]. Individual reactions were performed with approximately 30 ng DNA template, 3.2 pmol primer, 1 μl terminator mix, and 5 μl Better Buffer (The Gel Co.) in a total volume of 15 μl. PCR sequencing profiles consisted of an initial denaturation step of 4 min at 96°C, followed by 30 cycles of 30s at 96°C, 15s at 50°C, and 4 min at 60°C. Cycle sequencing was performed with BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems). Sequencing was completed on an automated sequencer, ABI 3130, manufactured by Applied Biosystems. The consensus sequences were assembled and edited using Geneious 4.8.5 [29] and aligned with Muscle algorithm implemented within Geneious 4.8.5. The generated haplotypes were deposited in GenBank under accession numbers KT160318 to KT160329.

Population genetics analyses

The relative nucleotide composition, number of polymorphic sites, haplotype diversity (h), nucleotide diversity (π), and number of pairwise nucleotide differences among populations were calculated using ARLEQUIN 3.5.1.3 software [30]. To estimate levels of genetic divergence among populations of C. longimanus index Φ_ST was calculated using the Analysis of Molecular Variance (AMOVA) [31] under the nucleotide evolution model of Tamura-Nei [32]. The Φ_ST estimates were tested with 1.000 non-parametric bootstrap pseudo-replicates using the ARLEQUIN 3.5.1.3 software, and the p-values were adjusted for simultaneous pairwise comparisons using the sequential Bonferroni correction [33]. The AMOVA test considered the possibilities of the maximum and minimum number of populations using hypothetic groups of population structure. Groupings that were most and least inclusive were tested, first by geographic proximity of the collection points, considering maritime currents; and secondly, by using the geographical distribution and correlation between the North and South, and East and West Atlantic. A minimum-spanning haplotype network was estimated using the Network 4.611 program [34]. The mismatch analysis in ARLEQUIN 3.5.1.3 included a raggedness index to determine goodness-of-fit to a unimodal distribution [35].

Results

Among the 215 specimens of C. longimanus analyzed, we obtained a total matrix with 724 bp with 9 polymorphic sites, yielding 12 haplotypes (Table 1).

[Figure omitted. See PDF.]

Table 1. Polymorphisms found in the twelve haplotypes of Carcharhinus longimanus along several regions of the Atlantic and Indian Oceans.

https://doi.org/10.1371/journal.pone.0155623.t001

In the overall analyses, we found low relative levels of haplotype diversity h = 0.5953 and nucleotide diversity π = 0.0013, with the greatest diversities found in the Southwest Atlantic (h = 0.8000, π = 0.00184) and in the Meso-Atlantic (h = 0.7571, π = 0.00192. Haplotype numbers three (H3 = 131 individuals) and one (H1 = 32 individuals) were the most common and were found in all sampled regions, representing 75.8% of the analyzed oceanic whitetip shark specimens. Numbers of polymorphic sites, haplotypes, nucleotide diversity and haplotypic diversity for different samples are shown in Tables 2 and 3.

[Figure omitted. See PDF.]

Table 2. Geographical distribution of haplotypes of Carcharhinus longimanus delimited by the different areas of the Atlantic and Indian Oceans.

https://doi.org/10.1371/journal.pone.0155623.t002

[Figure omitted. See PDF.]

Table 3. Population statistics of Carcharhinus longimanus—n, number of individuals; P, polymorphic sites; Nh, number of haplotypes; h, haplotype diversity; π, nucleotide diversity.

https://doi.org/10.1371/journal.pone.0155623.t003

The AMOVA test was adopted following hypothetic groups of population structure. While examining different scenarios of structuring, a greater difference (ɸ_ST = 0.1039, P <0.001) was observed when populations were grouped in two geographic regions: Eastern (Northeast Atlantic, Southeast Atlantic and Meso-Atlantic) and Western (Northwest Atlantic, Western Equatorial Atlantic and Southwest Atlantic) Atlantic Ocean. In this case, the indices related to genetic variability were h = 0.716 and π = 0.00165 in West Atlantic, and h = 0.471 and π = 0.00104 in the East Atlantic. Due to the low sampling rate in the the Indian Ocean, the AMOVA did not detect a statistically significant structure when left as a distinct group. However, in the F_ST pairwise analysis an absence of structure was observed between the Indian Ocean and the East Atlantic groups (Table 4).

[Figure omitted. See PDF.]

Table 4. Differentiation index (F_ST) between pairs of sampled regions in the Atlantic and Indian Oceans.

Pairwise F_ST numbers are below the diagonal line, and the significance of the p-values are represented above the diagonal (+ represents a statistically significant difference with p-value < 0.01).

https://doi.org/10.1371/journal.pone.0155623.t004

A statistical parsimony haplotype network was constructed using the Network 4.611 program. More than 60% analyzed sequences shared a single haplotype (H3), which was found in all sampled regions (Fig 1).

Discussion

Population genetics of the oceanic whitetip shark

The results show low levels of genetic diversity for the oceanic whitetip shark in the Atlantic and also in the Indian Ocean, although the latter has been evaluated with a small number of individuals and in a relatively small area. Population genetic studies with marine vertebrates, such as bony fishes, cetaceans and other elasmobranchs (including pelagic sharks) also tend to report low values of nucleotide and haplotypic diversity in the Atlantic Ocean [36–39]. Although low levels of nucleotide diversity are not considered standard for elasmobranchs, similar results to those obtained for oceanic whitetip sharks have been found in other highly migratory pelagic species including basking shark (Cetorhinus maximus), crocodile shark (Pseudocarcharias kamoharai), the whale shark (Rhincodon typus), pelagic thresher shark (Alopias pelagicus) and the white shark (Carcharodon carcharias) [25, 26, 38, 40–44]. The combination of low haplotypic diversity with a single haplotype shared by most of the individuals is also found in the nurse shark (Ginglymostoma cirratum) [39], white shark (Carcharodon carcharias) [43] and sicklefin lemon shark (Negaprion acutidens) [44].

Considering two distinct groups for the oceanic whitetip shark in the Atlantic Ocean and based on the hypothetic population groupings and genetic structuring simulations, restrictions to gene flow with moderate genetic divergence were observed, with one population occurring in the West and another in the East Atlantic. Population structure between the eastern and western portions of an oceanic basin was observed by Cardeñosa et al [41] for the thresher shark with strong differentiation in the Pacific. In the latest study on population genetic structure of the white shark, presented by O'Leary et al. [43], a similar structure was observed to the one in our study, being the Φ_ST = 0.10 with sequences of mtDNA control region and F_ST = 0.1057 with microsatellite marker, between the northwest Atlantic and southern Africa.

It would not be surprising to observe panmitic populations given the migration potential of this species. However, ocean currents and their interactions produce the oceanic barriers, organismal dispersal potential, behavior and organismal environmental tolerances which all contribute to patterns of diversification in the oceans [45, 46]. Scientific findings indicate that the western Atlantic is an area of both origin and accumulation of biodiversity [47]. This is supported by population growth rates of C. longimanus as rates of haplotype and nucleotide diversity are approximately 35% higher in western Atlantic. The haplotype network in star contraction generated from the haplotype H3 suggests an expansion event, possibly from the individuals of the eastern Atlantic, where the highest rates were observed in diversity with higher frequency on the H3.

Currently, very little is known about the movements and habitat use of oceanic whitetips in the Atlantic. The most comprehensive movement study of this species used Pop-up satellite archival tags applied to 11 mature females and one mature male,.near Cat Island in the central Bahamas. In that study, the maximum individual displacement from the tagging site dispersed 290–1940 km after 30–245 days at liberty, with individuals moving to several different destinations. The tag affixed to the male shark never reported. All of the tagged oceanic whitetips remained within 500 km of their tagging location for the first 30 days of their tracks and none of them moved away from the northwestern Atlantic [27].

The shortest route between the western and eastern Atlantic is between Brazil and Guinea-Bissau, requiring an oceanic crossing measuring ~2400 km. Although this distance does not seem to be a barrier to potential oceanic whitetip migration, such a trans-Atlantic route does not seem to be an aspect of the female behavior of the species. In those regions the oceanic whitetip shark populations in the western and east Atlantic are genetically differentiated from one another at mitochondrial loci, at least by female lineages. Philopatry is one of the factors that could also influence the reduction of gene flow across the Atlantic, by the migration of females to their original birth place [27]. Philopatry is well documented in many other sharks reviewed in Chapman et al. [48] and trans-Atlantic structure may have developed in oceanic whitetips because females remain within or return to give birth on one side of the basin or the other. A survey of biparentally inherited genetic markers is needed to determine whether or not there is male-mediated gene flow across the Atlantic.

The haplotype network denotes genetic differences among individuals without showing correlations between the different branches and their geographic distributions, with a characteristic pattern of an unstructured population. Although dispersion of fishes between the Indian and Atlantic Oceans are rarely observed, mtDNA data show that some species of tropical and sub-tropical fish found in the Indian Ocean may cross the barrier of the Benguela current [49]. The water transport between the Indian and Atlantic Oceans is summarized in [50]. A strong genetic connectivity between the Indian Ocean group and the Eastern Atlantic groups through the pairwise F_ST, despite low statistical significance, highlight the dispersal capability of C. longmanus. Considering the movement of individuals between Indian and eastern Atlantic, the absence of significant genetic structure can indicate the existence of only one genetic stock of oceanic whitetip sharks around the African continent, both in the eastern Atlantic and western Indian Ocean. The warm Agulhas current passes the Cape of Good Hope and is incorporated into the cool Benguela current. The Agulhas current is one of the strongest currents in the world, and it may help the transport of oceanic whitetip sharks from the southwest Indian Ocean to the Southeast Atlantic. However, the question on whether the flow of these sharks can also occur in the reverse order remains unknown.

More genetic studies and their importance in the evolutionary history of elasmobranchs with oceanic distribution, as well as for the management of commercially exploited shark species, are needed. Phylogeographic studies presenting extensive sampling and appropriate geographic scales including habitats with different ecological characteristics (such as oceanic and continental, tropical or subtropical) have a greater chance of detecting evolutionarily significant units, which would facilitate the development and application of appropriate conservation measures [51].

Conservation of the oceanic whitetip shark

Though shark populations are often impacted by overfishing, pelagic species, including C. longimanus, ubiquitous in all tropical oceans, with its susceptibility to incidental longline capture, make it prone to population declines, prompting the IUCN to designate it “vulnerable”. Considering the important relationship between the sustainability of a natural population and their genetic diversity, there is an urgent need to uncover information about global levels of genetic diversity and geographical distribution as well as patterns of gene flow and the detection of significant evolutionary units.

Here we present evidence that two populations of oceanic whitetip sharks exist in the Atlantic, specifically one in the western and another in the eastern regions of the basin. Despite their highly migratory nature, barriers to the gene flow of oceanic whitetips in the Atlantic result in two genetically distinct and demographically independent populations. This structure should be incorporated into assessments and monitoring of this species. The factors that restrict gene flow in the Atlantic may be absent between the east Atlantic and parts of the Indian Ocean; as there appears to be connectivity between those two regions. Isolated populations, especially if reduced by fishing or other sources of mortality, are vulnerable to inbreeding as reproducing adults have an elevated probability of mating with relatives [52]. Nonrandom matings generate deviations in allele frequencies with heterozygous deficiencies. Therefore, resource managers, tasked with the conservation of this species, should recognize that intrinsic reproductive barriers need to be considered concurrently with immigration rates, for sustaining at-risk populations.In addition to the low levels of genetic variability found for whitetip sharks throughout the study area, an important difference was observed between the two populations, with the haplotype diversity of the eastern population 34.2% less than the western population. This difference was even more striking for the nucleotide diversity with 36.9% less in the eastern population. These low genetic variability rates found may represent a dramatic risk to the adaptive potential of the species leading to a weaker ability to respond to environmental changes, and consequently, could promote extinction of some lineages in the future. Describing differences in indices of genetic diversity allows for managers to prioritize conservation efforts between populations with the highest diversity index, maintaining evolutionary potential, or populations with the lowest genetic variability, in an effort to prevent further inbreeding and reduced reproductive success. In order to prevent further population declines, we suggest that global genetic diversity should be preserved for all haplotypes through multinational cooperation, and particularly for the two distinct populations of oceanic whitetips identified in the Atlantic.

Acknowledgments

Sampling from the Portuguese fishery in the Atlantic and Indian Oceans was conducted by the Portuguese Institute for the Ocean and Atmosphere (IPMA), Portugal, within the scope of the European Data Collection Framework (PNAB/DCF). The authors thank all collaborators who contributed with samples for this study, particularly Annabelle Brooks, Sean Williams, Debra Abercrombie, Lance Jordan and Mark Bond. Funding for north Atlantic sampling was provided by grants from The Moore Bahamas Foundation and Save our Seas Foundation. Thanks are also due to the technicians from IPMA for collecting and carefully handling samples during the Portuguese pelagic longline fishery observer program. R. Coelho is supported by an Investigador-FCT contract from the Portuguese Foundation for Science and Technology (FCT, Fundação para a Ciência e Tecnologia) supported by the EU European Social Fund and the Programa Operacional Potencial Humano (Ref: IF/00253/2014). The Microwave Telemetry, Inc. provided support in the form of salaries for author L Howey-Jordan. This work was funded by FAPESP—Fundação de Amparo à Pesquisa do Estado de São Paulo (Ref: BIOTA 2011/23787-0) and coordinated by FFM. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author Contributions

Conceived and designed the experiments: SMC RC CO FF FFM. Performed the experiments: SMC NJM. Analyzed the data: SMC RC FFM. Contributed reagents/materials/analysis tools: SMC RC DC EJB DF FHVH CO MNS FF FFM. Wrote the paper: SMC RC DC LHJ EJB NJM FHVH CO MNS FF FFM.

Citation: Camargo SM, Coelho R, Chapman D, Howey-Jordan L, Brooks EJ, Fernando D, et al. (2016) Structure and Genetic Variability of the Oceanic Whitetip Shark, Carcharhinus longimanus, Determined Using Mitochondrial DNA. PLoS ONE 11(5): e0155623. https://doi.org/10.1371/journal.pone.0155623

References

1. 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 Journal of Marine Science: Journal du Conseil. 2000;57(3):476–94.

2. Baum JK, Myers RA, Kehler DG, Worm B, Harley SJ, Doherty PA. Collapse and conservation of shark populations in the Northwest Atlantic. Science. 2003;299(5605):389–92. pmid:12532016

3. Burgess GH, Beerkircher LR, Cailliet GM, Carlson JK, Cortés E, Goldman KJ, et al. Is the collapse of shark populations in the Northwest Atlantic Ocean and Gulf of Mexico real? Fisheries. 2005;30(10):19–26.

4. Camhi M, Fowler S, Musick J, Bräutigam A, S. F. Sharks and their relatives. 1998. 20 p.

5. IUCN. The IUCN Red List of Threatened Species 2014 [cited Version 2014.3].

6. Dulvy NK, Baum JK, Clarke S, Compagno LJ, Cortes E, Domingo A, et al. You can swim but you can't hide: the global status and conservation of oceanic pelagic sharks and rays. Aquatic Conservation: Marine and Freshwater Ecosystems. 2008;18(5):459.

7. Backus RH, Springer S, Arnold EL. A contribution to the natural history of the white-tip shark, Pterolamiops longimanus (Poey). Deep Sea Research (1953). 1956;3(3):178–88.

8. Bonfil R. Overview of world elasmobranch fisheries: Food & Agriculture Org.; 1994.

9. Castro JI, Woodley CM, Brudek RL. A preliminary evaluation of the status of shark species: Food & Agriculture Org.; 1999.

10. Compagno LJV. Species Catalogue. Sharks of the World. A, Annotated and Illustrated Catalogue of Shark Species Known to Date. Part 2. Carcharhiniformes. 1984. 251–655 p.

11. Hoenig J, Gruber S, Pratt H, Taniuchi T. Elasmobranchs as living resources: advances in the biology, ecology, systematics, and the status of the fisheries. Elasmobranchs as living resources: advances in the biology, ecology, systematics, and the status of the fisheries. 1990.

12. Coelho R, Fernandez-Carvalho J, Lino PG, Santos MN. An overview of the hooking mortality of elasmobranchs caught in a swordfish pelagic longline fishery in the Atlantic Ocean. Aquatic Living Resources. 2012;25(04):311–9.

13. Mayr E. Animal species and evolution. Animal species and their evolution. 1963.

14. Baum J, Medina E, M MJS. Carcharhinus longimanus. IUCN Red List of Threatened Species. 2006; IUCN 2010(Version 2010.1).

15. Rice J, Harley S. Stock assessment of oceanic whitetip sharks in the western and central Pacific Ocean. Scientific Committee Eighth Regular Session of WCPFC, Busan Republic of Korea. 2012.

16. Cortés E, Domingo A, Miller P, Forselledo R, Mas F, Arocha F, et al. K. Yokawa. 2012. Expanded ecological risk assessment of pelagic sharks caught in Atlantic pelagic longline fisheries. Report of the Standing Committed on Research and Statistics (SCRS), ICCAT, Madrid, Spain, October 1–5, 2012. SCRS/2012/167.

17. Murua H, Cohelo R, Neves M, Arrizabalaga H, Yokawa K, Romanov E, et al. Preliminary Ecological Risk Assessment (ERA) for shark species caught in fisheries managed by the Indian Ocean Tuna Commission (IOTC). IOTC Techical Paper. 2012;1:1–26.

18. Carvalho GR, Hauser L. Molecular genetics and the stock concept in fisheries. Reviews in Fish Biology Fisheries 1994;4:326–50.

19. Lindenmayer D, Lacy R. A simulation study of the impacts of population subdivision on the mountain brushtail possum Trichosurus caninus Ogilby (Phalangeridae: Marsupialia) in south-eastern Australia. I. Demographic stability and population persistence. Biological Conservation. 1995;73(2):119–29.

20. Clarke CR, Karl SA, Horn RL, Bernard AM, Lea JS, Hazin FH, et al. Global mitochondrial DNA phylogeography and population structure of the silky shark, Carcharhinus falciformis. Marine Biology. 2015;162(5):945–55.

21. Mendonca FF, Oliveira C, Gadig OBF, Foresti F. Diversity and genetic population structure of the Brazilian sharpnose shark Rhizoprionodon lalandii. Aquatic Conservation-Marine and Freshwater Ecosystems. 2013;23(6):850–7.

22. Karl SA, Castro ALF, Lopez JA, Charvet P, Burgess GH. Phylogeography and conservation of the bull shark (Carcharhinus leucas) inferred from mitochondrial and microsatellite DNA. Conservation Genetics. 2011;12(2):371–82.

23. Portnoy DS, McDowell JR, Heist EJ, Musick JA, Graves JE. World phylogeography and male‐mediated gene flow in the sandbar shark, Carcharhinus plumbeus. Molecular Ecology. 2010;19(10):1994–2010. pmid:20406387

24. Ovenden JR, Kashiwagi T, Broderick D, Giles J, Salini J. The extent of population genetic subdivision differs among four co-distributed shark species in the Indo-Australian archipelago. BMC Evolutionary Biology. 2009;9.

25. Hoelzel AR, Shivji MS, Magnussen J, Francis MP. Low worldwide genetic diversity in the basking shark (Cetorhinus maximus). Biol Lett. 2006;2(4):639–42. pmid:17148309

26. Ferrette BLD, Mendonca FF, Coelho R, de Oliveira PGV, Hazin FHV, Romanov EV, et al. High Connectivity of the Crocodile Shark between the Atlantic and Southwest Indian Oceans: Highlights for Conservation. Plos One. 2015;10(2).

27. Howey-Jordan LA, Brooks EJ, Abercrombie DL, Jordan LK, Brooks A, Williams S, et al. Complex movements, philopatry and expanded depth range of a severely threatened pelagic shark, the oceanic whitetip (Carcharhinus longimanus) in the western North Atlantic. PLoS One. 2013;8(2):e56588. pmid:23437180

28. Mendonça FF, Oliveira C, Gadig OBF, Foresti F. Populations analysis of the Brazilian Sharpnose Shark Rhizoprionodon lalandii (Chondrichthyes: Carcharhinidae) on the São Paulo coast, Southern Brazil: inferences from mt DNA sequences. Neotropical Ichthyology. 2009;7(2):213–6.

29. 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;28(12):1647–9. pmid:22543367

30. Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular ecology resources. 2010;10(3):564–7. pmid:21565059

31. Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131(2):479–91. pmid:1644282

32. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular biology and evolution. 1993;10(3):512–26. pmid:8336541

33. Rice WR. Analyzing tables of statistical tests. Evolution. 1989:223–5.

34. Bandelt H-J, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Molecular biology and evolution. 1999;16(1):37–48. pmid:10331250

35. Harpending H. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Human biology. 1994:591–600. pmid:8088750

36. Bradman H, Grewe P, Appleton B. Direct comparison of mitochondrial markers for the analysis of swordfish population structure. Fisheries Research. 2011;109(1):95–9.

37. Mesnick SL, Taylor BL, Archer FI, Martien KK, Trevino SE, Hancock-Hanser BL, et al. Sperm whale population structure in the eastern and central North Pacific inferred by the use of single‐nucleotide polymorphisms, microsatellites and mitochondrial DNA. Molecular Ecology Resources. 2011;11(s1):278–98.

38. Phillips NM, Chaplin JA, Morgan DL, Peverell SC. Population genetic structure and genetic diversity of three critically endangered Pristis sawfishes in Australian waters. Marine Biology. 2011;158(4):903–15.

39. Karl SA, Castro AL, Garla RC. Population genetics of the nurse shark (Ginglymostoma cirratum) in the western Atlantic. Marine Biology. 2012;159(3):489–98.

40. Schmidt JV, Schmidt CL, Ozer F, Ernst RE, Feldheim KA, Ashley MV, et al. Low genetic differentiation across three major ocean populations of the whale shark, Rhincodon typus. PLoS One. 2009;4(4):e4988. pmid:19352489

41. Cardeñosa D, Hyde J, Caballero S. Genetic Diversity and Population Structure of the Pelagic Thresher Shark (Alopias pelagicus) in the Pacific Ocean: Evidence for Two Evolutionarily Significant Units. 2014.

42. Castro A, Stewart B, Wilson S, Hueter R, Meekan M, Motta P, et al. Population genetic structure of Earth's largest fish, the whale shark (Rhincodon typus). Molecular Ecology. 2007;16(24):5183–92. pmid:18092992

43. O'Leary SJ, Feldheim KA, Fields AT, Natanson LJ, Wintner S, Hussey N, et al. Genetic Diversity of White Sharks, Carcharodon carcharias, in the Northwest Atlantic and Southern Africa. J Hered. 2015.

44. Schultz JK, Feldheim KA, Gruber SH, Ashley MV, McGovern TM, Bowen BW. Global phylogeography and seascape genetics of the lemon sharks (genus Negaprion). Mol Ecol. 2008;17(24):5336–48. pmid:19121001

45. Joyeux JC, Floeter SR, Ferreira CEL, Gasparini JL. Biogeography of tropical reef fishes: the South Atlantic puzzle. Journal of Biogeography. 2001;28(7):831–41.

46. Hellberg ME. Gene flow and isolation among populations of marine animals. 2009.

47. Rocha LA, Bowen BW. Speciation in coral‐reef fishes. Journal of Fish Biology. 2008;72(5):1101–21.

48. Chapman DD, Feldheim KA, Papastamatiou YP, Hueter RE. There and Back Again: A Review of Residency and Return Migrations in Sharks, with Implications for Population Structure and Management. Annual review of marine science. 2015;7:547–70. pmid:25251267

49. Rocha LA, Robertson DR, Rocha CR, Van Tassell JL, Craig MT, Bowen BW. Recent invasion of the tropical Atlantic by an Indo-Pacific coral reef fish. Mol Ecol. 2005;14(13):3921–8. pmid:16262848

50. Tomczak M, Godfrey JS. Regional oceanography: an introduction: Elsevier; 2013.

51. Rocha L, Craig M, Bowen B. Phylogeography and the conservation of coral reef fishes. Coral Reefs. 2007;26(3):501–12.

52. Wright S. Evolution and the genetics of populations, volume 3: experimental results and evolutionary deductions: University of Chicago press; 1984.