# From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator

Floriaan Devloo-delva, Christopher P Burridge, Peter M Kyne, Juerg M Brunnschweiler, Demian D Chapman, Patricia Charvet, Xiao Chen, Geremy

Cliff, Ryan Daly, J. Marcus Drymon, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Floriaan Devloo-delva, Christopher P Burridge, Peter M Kyne, Juerg M Brunnschweiler, Demian D Chapman, et al. From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator. Ecology and Evolution, 2023, 13 (2), 10.1002/ece3.9837. hal-04051570

HAL Id: hal-04051570

https://hal.univ-reunion.fr/hal-04051570

Submitted on 30 Mar 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

RESEARCH ARTICLE

Ecology and Evolution

Open Access

WILEY

# From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator

Floriaan Devloo-Delva<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID

Christopher P. Burridge<fnref n="3" /> ID

Peter M. Kyne<fnref n="4" /> ID

Juerg M. Brunnschweiler<fnref n="5" /> ID

Demian D. Chapman<fnref n="6" /> ID

Patricia Charvet<fnref n="7" /> ID

Xiao Chen<fnref n="8" /> |

Ryan Daly<fnref n="10" />,<fnref n="11" /> ID

J. Marcus Drymon<fnref n="12" />,<fnref n="13" /> ID

Mario Espinoza<fnref n="14" /> ID

Daniel Fernando<fnref n="15" /> ID

Laura Garcia Barcia<fnref n="6" /> ID

Kerstin Glaus<fnref n="16" /> ID

Blanca I. González-Garza<fnref n="17" /> |

Michael I. Grant<fnref n="18" /> ID

Rasanthi M. Gunasekera<fnref n="1" /> ID

Sebastian Hernandez<fnref n="19" />,<fnref n="20" /> ID

Susumu Hyodo<fnref n="21" /> ID

Rima W. Jabado<fnref n="18" />,<fnref n="22" /> ID

Sébastien Jaquemet²³ ID

Grant Johnson²⁴

James T. Ketchum<fnref n="17" />,²⁵ ID

Hélène Magalon²³ ID

James R. Marthick²⁶

Frederik H. Mollen²⁷ ID

Stefano Mona²⁸,²⁹ ID

Gavin J. P. Naylor³⁰ ID

John E. G. Nevill³¹

Nicole M. Phillips³² ID

Richard D. Pillans³³ ID

Bautisse D. Postaire²³ ID

Amy F. Smoothey³⁴ ID

Katsunori Tachihara³⁵ |

Bree J. Tillet³⁶ ID

Jorge A. Valerio-Vargas<fnref n="14" /> ID

Pierre Feutry<fnref n="1" /> ID

<footnote n="1">¹Oceans and Atmosphere, CSIRO, Hobart, Tasmania, Australia</footnote>

<footnote n="2">²Quantitative Marine Science, Institute for Marine and Antarctic Studies, University of Tasmania, Hobart, Tasmania, Australia</footnote>

<footnote n="3">³Discipline of Biological Sciences, School of Natural Sciences, University of Tasmania, Hobart, Tasmania, Australia</footnote>

<footnote n="4">⁴Research Institute for the Environment and Livelihoods, Charles Darwin University, Darwin, Northern Territory, Australia</footnote>

<footnote n="5">⁵Independent Researcher, Zurich, Switzerland</footnote>

<footnote n="6">⁶Department of Biological Sciences, Florida International University, North Miami, Florida, USA</footnote>

<footnote n="7">⁷Programa de Pós-graduação em Sistemática, Uso e Conservação da Biodiversidade, Universidade Federal do Ceará (PPGSis - UFC), Fortaleza, Brazil</footnote>

<footnote n="8">⁸College of Veterinary Medicine, South China Agricultural University, Guangzhou, China</footnote>

<footnote n="9">⁹KwaZulu-Natal Sharks Board, Umhlanga 4320, South Africa and School of Life Sciences, University of KwaZulu-Natal, Durban, South Africa</footnote>

<footnote n="10">¹⁰Oceanographic Research Institute, South African Association for Marine Biological Research, Point, Durban, South Africa</footnote>

<footnote n="11">¹¹South African Institute for Aquatic Biodiversity, Mkhanda, South Africa</footnote>

<footnote n="12">¹²Coastal Research and Extension Center, Mississippi State University, Biloxi, Mississippi, USA</footnote>

<footnote n="13">¹³Mississippi-Alabama Sea Grant Consortium, Ocean Springs, Mississippi, USA</footnote>

<footnote n="14">¹⁴Centro de Investigación en Ciencias del Mar y Limnología & Escuela de Biología, Universidad de Costa Rica, San Pedro de Montes de Oca, San José, Costa Rica</footnote>

<footnote n="15">¹⁵Blue Resources Trust, Colombo, Sri Lanka</footnote>

<footnote n="16">¹⁶Faculty of Science, Technology and Environment, School of Marine Studies, The University of the South Pacific, Suva, Fiji</footnote>

<footnote n="17">¹⁷Pelagios-Kakunja, La Paz, Mexico</footnote>

<footnote n="18">¹⁸College of Science and Engineering, Centre for Sustainable Tropical Fisheries and Aquaculture, James Cook University, Townsville, Queensland, Australia</footnote>

<footnote n="19">¹⁹Biomolecular Laboratory, Center for International Programs, Universidad VERITAS, San José, Costa Rica</footnote>

<footnote n="20">²⁰Sala de Colecciones, Facultad de Ciencias del Mar, Universidad Católica del Norte, Coquimbo, Chile</footnote>

<footnote n="21">²¹Laboratory of Physiology, Atmosphere and Ocean Research Institute, University of Tokyo, Kashiwa, Chiba, Japan</footnote>

<footnote n="22">²²Elasmo Project, Dubai, United Arab Emirates</footnote>

<footnote>This is an open access article under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>

<footnote>© 2023 The Authors. *Ecology and Evolution* published by John Wiley & Sons Ltd.</footnote>
**Correspondence**

Floriaan Devloo-Delva, Oceans and Atmosphere, CSIRO, Hobart, TAS, Australia.

Email: floriaan.devloo-delva@csiro.au

**Funding information**

Sea World Research and Rescue
Foundation, Grant/Award Number:
No27; National Environmental Science
Program, Marine Biodiversity Hub, Grant/
Award Number: A1; CSIRO, Grant/Award
Number: Ord River Research Off set

**Abstract**

The Bull Shark (*Carcharhinus leucas*) faces varying levels of exploitation around the world due to its coastal distribution. Information regarding population connectivity is crucial to evaluate its conservation status and local fishing impacts. In this study, we sampled 922 putative Bull Sharks from 19 locations in the first global assessment of population structure of this cosmopolitan species. Using a recently developed DNA-capture approach (DArTcap), samples were genotyped for 3400 nuclear markers. Additionally, full mitochondrial genomes of 384 Indo-Pacific samples were sequenced. Reproductive isolation was found between and across ocean basins (eastern Pacific, western Atlantic, eastern Atlantic, Indo-West Pacific) with distinct island populations in Japan and Fiji. Bull Sharks appear to maintain gene flow using shallow coastal waters as dispersal corridors, whereas large oceanic distances and historical land-bridges act as barriers. Females tend to return to the same area for reproduction, making them more susceptible to local threats and an important focus for management actions. Given these behaviors, the exploitation of Bull Sharks from insular populations, such as Japan and Fiji, may instigate local decline that cannot readily be replenished by immigration, which can in turn affect ecosystem dynamics and functions. These data also supported the development of a genetic panel to ascertain the population of origin, which will be useful in monitoring the trade of fisheries products and assessing population-level impacts of this harvest.

**KEYWORDS**

close-kin, DArTseq, DNA forensics, genetic connectivity, mitogenome, provenance

**TAXONOMY CLASSIFICATION**

Genomics, Population genetics

# 1 | INTRODUCTION

Understanding the population structure of a species and the barriers that disrupt dispersal is important to accurately assess the global

conservation status and manage the risk of local extinction. This is especially true for species of commercial importance (Begg et al., 1999) or conservation concern (Moritz, 1994), which are impacted disproportionately by anthropogenic or environmental pressures. Dispersal

<footnote n="23"><sup>23</sup>UMR ENTROPIE (Université de La Réunion, Université de Nouvelle-Calédonie, IRD, CNRS, IFREMER), Faculté des Sciences et Technologies, Université de La Réunion, Cedex 09, La Réunion, France</footnote>

<footnote n="24"><sup>24</sup>Department of Industry, Tourism and Trade, Aquatic Resource Research Unit, Darwin, Northern Territory, Australia</footnote>

<footnote n="25"><sup>25</sup>MigraMar, Olema, California, USA</footnote>

<footnote n="26"><sup>26</sup>Menzies Institute for Medical Research, University of Tasmania, Hobart, Tasmania, Australia</footnote>

<footnote n="27"><sup>27</sup>Elasmobranch Research, Bonheiden, Belgium</footnote>

<footnote n="28"><sup>28</sup>Institut de Systématique, Evolution, Biodiversité, ISYEB (UMR 7205), Muséum National d'Histoire Naturelle, CNRS, Sorbonne Université, EPHE, Université des Antilles, Paris, France</footnote>

<footnote n="29"><sup>29</sup>EPHE, PSL Research University, Paris, France</footnote>

<footnote n="30"><sup>30</sup>Florida Museum of Natural History, University of Florida, Gainesville, Florida, USA</footnote>

<footnote n="31"><sup>31</sup>Environment Seychelles, Victoria, Seychelles</footnote>

<footnote n="32"><sup>32</sup>School of Biological, Environmental and Earth Sciences, The University of Southern Mississippi, Hattiesburg, Mississippi, USA</footnote>

<footnote n="33"><sup>33</sup>Oceans and Atmosphere, CSIRO, Dutton Park, Queensland, Australia</footnote>

<footnote n="34"><sup>34</sup>NSW Department of Primary Industries, Fisheries Research, Sydney Institute of Marine Science, Mosman, New South Wales, Australia</footnote>

<footnote n="35"><sup>35</sup>Laboratory of Fisheries Biology and Coral Reef Studies, Faculty of Science, University of Ryukyus, Nishihara, Okinawa, Japan</footnote>

<footnote n="36"><sup>36</sup>Translational Research Institute, University of Queensland Diamantina Institute, Brisbane, Queensland, Australia</footnote>
can promote genetic connectivity across patches of suitable habitat (Ronce, 2007), but physical barriers and behaviors can ultimately limit gene flow, even when dispersal potential is high, and result in demographically independent populations (Waples & Gaggiotti, 2006). When genetically isolated populations are reduced in size due to unsustainable harvest, there is a risk of inbreeding depression and the loss of genetic diversity, without the chance to be 'rescued' by individuals dispersing from adjacent populations, increasing the likelihood of population extinctions (Frankham et al., 2017).

In marine taxa, gene flow tends to be restricted by environmental or biogeographic barriers, movement ecologies, and habitat preferences (Bowen et al., 2016; Dudgeon et al., 2012; Hirschfeld et al., 2021; Rocha et al., 2007). In large-bodied coastal species with global distributions, such as marine turtles, cetaeans, and many elasmobranchs (sharks and rays), large-scale marine biogeographic barriers shape population genetic structure (Dutton et al., 1999; Fontaine et al., 2007). For example, the Scalloped Hammerhead (Sphyrna lewini) shows genetic connectivity along the continental margins, yet limited gene flow across the East Pacific barrier, the Mid-Atlantic barrier, and the Isthmus of Panama (Daly-Engel et al., 2012; Green, Appleyard, et al., 2022). The permeability of environmental barriers, such as the Indo-Australian Archipelago, changes across time and space, and consequently determines the observed distribution of genetic variation (Cowman & Bellwood, 2013). The effect of these ocean barriers on the spatial structuring of populations is essential knowledge for management, given that threats, such as overfishing and habitat modification, should be assessed and managed at biologically relevant spatial scales. For vagile marine taxa this often requires cooperative strategies between nations (e.g., International Whaling Commission, Indian Ocean Tuna Commission, and Western and Central Pacific Fisheries Commission).

Delineation of population structure in species with high mobilities and large population sizes have recently been improved with the use of genomic data (Layton et al., 2020; Luikart et al., 2019; Oleksiak & Rajora, 2020; Ovenden et al., 2018). Complexity-reduction genome-scan methods, such as Diversity Arrays Technology sequencing (DArTseq, Jaccoud et al., 2001) and targeted DNA-capture approaches, including Rapture and DArTcap (Ali et al., 2016; Feutry et al., 2020), have been widely used to assess genetic diversity and reproductive connectivity in natural populations (e.g., Green et al., 2019; Komoroske et al., 2019). Recently, these methods have also been applied to DNA forensics or traceability studies in nonmodel species, with the objective to identify species, sex, provenance, and close-kin relationships (e.g., Arenas et al., 2017; Feutry et al., 2017; Nielsen et al., 2012; Stovall et al., 2018). Moreover, in taxa with slow mitochondrial DNA (mtDNA) mutation rates, such as elasmobranchs and marine turtles (Avise et al., 1992; Martin et al., 1992), the sequencing of full mitochondrial genomes (mitogenomes) instead of single genes has improved the fine-scale resolution of matrilineal population structure (Feutry et al., 2014).

The Bull Shark (Carcharhinus leucas) is a cosmopolitan species that occupies tropical, subtropical, and temperate coastal waters and has an important ecological role in freshwater, estuarine, and marine environments (Matich et al., 2011; Smoothey et al., 2019; Trystram et al., 2017). This species experiences variable degrees of exploitation within its range and is assessed as Vulnerable on the IUCN Red List of Threatened Species (Rigby et al., 2021). The main threats affecting Bull Sharks are small- and large-scale fisheries for meat and fins (Glaus et al., 2015; Holmes et al., 2009), and shark control programs which directly target this species (Blaison et al., 2015; Dudley & Simpfendorfer, 2006; Niella et al., 2021). Because of its global catch, the Bull Shark is also frequently found in the international shark fin and meat trade, raising questions regarding the origin and trade routes of the fished products (Cardeñosa et al., 2022; Clarke et al., 2006; Fields et al., 2018). Unsustainable exploitation may result in population declines and negative ecosystem impacts (Ferretti et al., 2010; MacNeil et al., 2020).

The Bull Shark has shown capacity for long-distance coastal movement (Brunnschweiler et al., 2010; Daly et al., 2014; Espinoza et al., 2016, 2021; Heupel et al., 2015; Lea et al., 2015) and genetic connectivity along continental shelves (Glaus et al., 2020; Pirog et al., 2019; Testerman, 2014). However, females exhibit reproductive philopatry to estuarine habitats at small spatial scales (~100km; Karl et al., 2011; Sandoval Laurrabaquio-Alvarado et al., 2021; Tillett et al., 2012). Genetic connectivity of Bull Sharks has been studied within specific regions, such as the western Atlantic (Karl et al., 2011; Sandoval Laurrabaquio-Alvarado et al., 2019, 2021) and the Indo-West Pacific (Deng et al., 2019; Glaus et al., 2020; Kitamura et al., 1996; Pirog et al., 2019; Tillett et al., 2012). Yet, no studies have investigated genetic connectivity across its global distribution, which can provide crucial data to assess the population-level impacts of threats and identify the origin of fisheries products.

This study aims to apply genomic techniques (full mitogenomes and nuclear Single Nucleotide Polymorphisms or SNPs) to investigate the population genetic structure of Bull Sharks at global and local scales. We hypothesize that gene flow is limited by large-scale biogeographic barriers, but that connectivity occurs along continuous coastlines. We further evaluate the forensic power of genomic data to assign sample provenance and establish a diagnostic SNP panel to aid monitoring the origin and trade of fisheries products and assess population-level impacts of this trade. To achieve these goals, we analyze a global set of samples with three different genomic sequencing approaches (DArTseq, DArTcap, and mitogenomes), which allows us to assess population structure and gene flow across putative barriers for dispersal (such as open-ocean expanses, strong temperature gradients, and historical land-bridges). At a smaller spatial scale, contemporary reproductive connectivity is estimated by examining the spatial distribution of closely related individuals, such as full siblings and cross-cohort half siblings.
2 | MATERIALS AND METHODS

2.1 | Sample collection and DNA extraction

A total of 922 putative Bull Shark samples (muscle or fin clip) were
collected between 1980 and 2019 from 19 different countries or
water bodies (termed 'sampling locations') around the globe in all
major ocean basin regions: the eastern Pacific (E-PAC), western
Atlantic (W-ATL), eastern Atlantic (E-ATL), and the Indo-West Pacific
(IWP; Figure 1). Samples from Brazil, eastern Indian Ocean (E-IO),
Fiji, and numerous Australian samples were sourced from previ-
ously published genetic studies (n = 175; Glaus et al., 2020; Karl
et al., 2011; Pirog et al., 2019; Tillett et al., 2012) while all other
samples were novel. The samples from Japan included 10 individuals
from the aquarium (Okinawa, locally sourced animals) with known
pedigree and 38 wild-caught samples from the Urauchi River. Each
sampling location had an approximately equal sex ratio, and total
length (TL) ranged from 26.4 (in utero) to 406.0 cm (Supporting
Information section 2). Sharks smaller than 150cm TL (68% of all
individuals) were considered juveniles with limited dispersal capac-
ity (Heupel et al., 2015; Pillans & Franklin, 2004; Simpfendorfer &
Milward, 1993). DNA was extracted using the Qiagen Blood and
Tissue kit following the standard protocol (Qiagen Inc., Valencia,

California, USA). DNA quality and quantity was assessed on a 1%
agarose gel, stained with SYBR safe (Invitrogen, USA), and with the
NanoDrop ND-8000 spectrophotometer (Thermo Fisher Scientific,
Waltham, Massachusetts, USA).

2.2 | SNP bait design and DArTcap genotyping

Initially, 188 samples with a minimum of 10 samples per location
were genotyped following the DArTseq approach (according to
Feutry et al., 2017). De novo SNP calling was performed with a
proprietary software (DArTsoft14). The resulting markers were
filtered using the 'filter_rad' function, as implemented in radiator
v1.1.5 (Gosselin et al., 2020). Filtering thresholds were chosen
based on the empirical distribution of each filtering parameter
(see Supporting Information sections 4–6). Briefly, the data were
filtered for: (1) low DArT reproducibility based on replicated li-
braries (technical replicates), (2) monomorphic markers, (3) high
degree of missing data, high levels of inferred heterozygosity,
and low sequencing coverage per individual, (4) low minor allele
count, (5) low and high SNP sequencing coverage to avoid unreli-
able SNP calls or paralogous sequences, (6) high missing data per
SNP, (7) too many SNPs per sequence due to suboptimal sequence

<figure><img src="image_4.png" /><figcaption>FIGURE 1 Map indicating the *Carcharhinus leucas* sampling locations with red circles and the known species range distribution in yellow. The sample sizes for the SNP data are underlined with the number of samples before/after data filtering. The sample sizes for the mitogenome data are in *italics* with the number of samples before/after data filtering. Putative barriers for gene flow are indicated by dashed lines. GOC, Gulf of California; COR, Costa Rica; BRZ, Brazil; CAR, Caribbean Sea; GOM, Gulf of Mexico; WNA, Western North Atlantic; SIL, Sierra Leone; SAF, South Africa; MOZ, Mozambique; RUN, Réunion Island; SEY, Seychelles; ARP, Arabian Peninsula; SRL, Sri Lanka; TAI, Thailand; IND, Indonesia; PNG, Papua New Guinea; AUS, Australia; JAP, Japan; and FIJ, Fiji. Australian sampling locations were presented as an additional inset: FZR, Fitzroy River; VIR, Victoria River; DAR, Daly River; ADR, Adelaide River; DWC, Darwin Coastal; SAR, South Alligator River; EAR, East Alligator River; BMB, Blue Mud Bay; ROR, Roper River; TOR, Towns River; WER, Wenlock River; TRI, Trinity Inlet; CLR, Clarence River; SYH, Sydney Harbor; and UNK, Australian fisheries samples from unknown origin.</figcaption></figure>
clustering, (8) short-distance linkage by keeping one SNP per sequence, (9) individual DNA contamination based on high proportion of heterozygous SNPs, (10) duplicated samples, and (11) departure from Hardy-Weinberg equilibrium. We further filtered out RAD-tags that were too short (<60 bp), with low complexity (>6 nucleotide repeats), or with high cluster size (>10) to allow efficient capture with the RNA baits. Of the markers satisfying the filtering criteria, 3200 population genetic markers were randomly selected for DNA-capture bait development. Secondly, the unfiltered data were run through the 'sexy_markers' function from radiator (Devloo-Delva et al., 2022) to identify polymorphic sex-linked markers (putatively located on sex chromosomes), which were subsequently included in the DNA capture panel. One biotinylated RNA MYbait (Arbor Bioscience, USA) was synthesized per RAD-tag. DArTcap hybridization and washing followed the MYbaits standard protocol (https://arorbiosci.com/wp-content/uploads/2020/08/myBaits_v5.0_Manual.pdf). The DArTcap-enriched libraries were sequenced on a HiSeq 2500 (Illumina, San Diego, California, USA) as described by Feutry et al. (2020).

## 2.3 | DArTcap SNP and individual filtering

Overall, 1014 sample libraries, including 92 technical replicates, were genotyped with the DArTcap protocol. Additional filtering was performed to remove any unspecific enrichment or low-quality loci. This filtering was performed with the 'filter-rad' function as described previously (see Supporting Information sections 7-12). Sex-linked markers were removed prior to population genetic analyses.

The filtering steps were applied with different thresholds to six different subsets of the data (Table 1). A hierarchical approach was used for the clustering methods (e.g., Vähä et al., 2007). The first DArTcap data set was filtered under less stringent thresholds (e.g., allow more missing data and higher levels of heterozygosity) to permit the identification of species that were not Bull Shark (DATA1: 922 sharks). The second dataset only contained confirmed Bull Sharks (DATA2: 868 sharks). Since unequal sample sizes can introduce bias in clustering algorithms (Foster et al., 2018; Puechmaille, 2016), in the third dataset, all sampling locations were included but Australia was represented by a random subset of 60 individuals (DATA3: 430 sharks). The fourth and fifth datasets are comprised of sampling locations in the W-ATL and IWP, respectively (DATA4: 117 sharks; DATA5: 732 sharks). Lastly, within the IWP, samples from Japan and Fiji were removed and 60 Australian samples were again randomly selected to equalize the sample sizes and investigate unbiased signals of more subtle structure (DATA6: 221 sharks).

## 2.4 | Mitogenome amplification and sequencing

Samples from the E-ATL, IWP, and E-PAC were used to investigate their matrilineal evolutionary history (sample sizes indicated

in Figure 1). The full mitochondrial genomes of 384 putative Bull Sharks were amplified with two primer pairs (A and B fragments; Supporting Information section 13.1). Polymerase chain reactions (PCR) were performed in 30µL volumes, following the standard proofreading Takara LA Taq protocol (Takara, Otsu, Shiga, Japan). PCR conditions were set to 1 min at 94°C for initial denaturation, then 40 cycles of denaturation (94°C, 30s), annealing (55°C, 30s), and extension (68°C, 10 min); concluded with a 10 min extension at 72°C. PCR products were cleaned following the Agencourt AMPure XP magnetic bead protocol (Beckman Coulter Inc., Indianapolis, Indiana, USA). Amplicons were quantified with a NanoDrop 8000 Spectrophotometer and the purified A and B fragments were pooled at equimolar concentrations for each individual. Subsequently, these amplicons were simultaneously fragmented and barcoded with the Nextera XT DNA Sample Preparation kits and 96 sample Nextera Index kit (Illumina). The libraries were quantified with the Qubit dsDNA BR assay kit (Life Technologies, Carlsbad, California, USA) and normalized. Libraries were then pooled and sequenced on a MiSeq desktop sequencer using the 2x250 bp paired-end reads MiSeq reagent kit v2 (Illumina).

## 2.5 | Mitogenome assembly and alignment

Demultiplexed fastq files were imported into Geneious prime software v2020.1 (Biomatters Ltd., Auckland, New Zealand) and the reads were paired. The Nextera adapters were removed, and the reads were quality trimmed at a phred score < 20 for a k-mer of 20 using the BBDuk tool as implemented in Geneious. Reads shorter than 60 bp after trimming were discarded from subsequent analyses. Merged reads for each individual were then mapped onto a previously published Bull Shark reference sequence (Chen et al., 2015) using "Map to Reference" tool in Geneious with the "high sensitivity" parameters and 10 iterations. The majority rule consensus (>50% of mapped reads for any single SNP, insertion, or deletion) for each shark was exported. All mitogenome sequences were aligned with the 'multiple align' tool and the MUSCLE algorithm (Edgar, 2004).

## 2.6 | Species identification

The mitogenomes were blasted (megablast) in Geneious against the GenBank nr/nt database. All sequences with a match of identity ≥98% to a database entry were assigned as that species. Samples without mitogenome information (only DArTcap genotypes, i.e., DATA1) were assigned to species using a principal component analysis (PCA) as implemented in *adegenet* v2.1.1 (Supporting Information section 7; Jombart & Ahmed, 2011). Here, the species identified based on mitogenomes served as a baseline to detect species clusters from PCA. To provide certainty, individuals needed to cluster on multiple PC axes with a mitochondrial-verified species. Where this was not possible, we marked the species as 'unknown'. Only samples identified as Bull Shark were included in further analyses.
<table><caption>TABLE 1 Description of the different *Carcharhinus leucas* datasets used in this study.</caption><thead><tr><th>Data set</th><th>Genotyping method</th><th>Data description</th><th>Usage</th><th>Before filtering</th><th>After filtering</th></tr></thead><tbody><tr><td>DARTseq1</td><td>DARTseq</td><td>Full global dataset including replicates and other species</td><td>For species identification (allow more missing data and higher levels of heterozygosity)</td><td>218 sharks; 250,945 SNPs</td><td>218 sharks; 28,418 SNPs</td></tr><tr><td>DARTseq2</td><td>DARTseq</td><td>Full global dataset with only Bull Shark samples</td><td>For DARTcap marker selection</td><td>182 sharks; 250,945 SNPs</td><td>149 sharks; 17,029 SNPs</td></tr><tr><td>DATA1</td><td>DARTcap</td><td>Full global dataset including replicates and other species</td><td>For species identification (allow more missing data and higher levels of heterozygosity)</td><td>1014 sharks; 37,537 SNPs</td><td>1014 sharks; 5053 SNPs</td></tr><tr><td>DATA2</td><td>DARTcap</td><td>Global dataset with confirmed Bull Sharks</td><td>For genetic diversity, F<sub>ST</sub>, and provenance assignment</td><td>868 sharks; 37,537 SNPs</td><td>769 sharks; 3409 SNPs</td></tr><tr><td>DATA3</td><td>DARTcap</td><td>Global Bull Shark dataset with a subset of 60 samples from Australia</td><td>Equalized sample size to allow unbiased clustering analyses</td><td>430 sharks; 37,537 SNPs</td><td>382 sharks; 1849 SNPs</td></tr><tr><td>DATA4</td><td>DARTcap</td><td>All western Atlantic Bull Shark samples</td><td>For genetic diversity estimates without bias from population structure</td><td>117 sharks; 37,537 SNPs</td><td>91 sharks; 931 SNPs</td></tr><tr><td>DATA5</td><td>DARTcap</td><td>All Indo-West Pacific Bull Shark samples</td><td>For genetic diversity estimates without bias from population structure</td><td>732 sharks; 37,537 SNPs</td><td>635 sharks; 3416 SNPs</td></tr><tr><td>DATA6</td><td>DARTcap</td><td>IWP data with a subset of 60 samples from Australia</td><td>Equalized sample size to allow unbiased clustering analyses</td><td>221 sharks; 37,537 SNPs</td><td>189 sharks; 1785 SNPs</td></tr></tbody></table>

Note: The data generated by DARTseq and DARTcap genotyping were subdivided to account for unequal sample sizes and hierarchical clustering. The number of sharks and SNPs before and after quality filtering are presented.
## 2.7 | Genetic diversity

Allelic richness ($A_r$), observed heterozygosity ($H_O$), and unbiased expected heterozygosity ($H_E$) were calculated for each sampling location with $n > 1$ using the R package *diveRsity* v1.9.90 (Keenan et al., 2013). Heterozygosity was calculated across the global dataset (including monomorphic and polymorphic SNPs) to yield an unbiased estimate of $H_O$ (see Schmidt et al., 2021). Inbreeding coefficients ($F_{IS}$) were calculated with *hierfstat* v0.04-22 (Goudet, 2005), using 1000 bootstraps to determine a 95% confidence interval. Mitochondrial nucleotide ($\pi_{mt}$) and haplotype (h) diversities per location, and parsimony haplotype networks were calculated with the *pegas* v0.14 package (Paradis, 2010). All analyses were performed in R 4.0.2 (R Core Team, 2020) and provided in the Supporting Information.

## 2.8 | Population structure

Fixation indices (mtDNA: sequence-based $\Phi_{ST}$, SNPs: $F_{ST}$) were calculated between all sampling locations ($n \ge 1$) and between/ across ocean basins with the 'popStructTest' function in the *strataG* v2.4.905 package (Archer et al., 2017), and their significance assessed by 1000 permutations. Next, dimensionality-reduction clustering analyses were conducted with *adegenet* (PCA and Discriminant Analysis of Principal Components, DAPC; Jombart & Ahmed, 2011; Jombart et al., 2010). Individuals were grouped using the successive K-means algorithm implemented in the 'find.clusters' function. The goodness of fit, determined by the Bayesian information criterion (BIC), was employed to find the best number of clusters (K). To avoid overfitting, the optimal number of principal components was selected through cross-validation with a 10% hold-out set and 1000 replicates for all DAPC analyses, where individuals were grouped according to their sample location.

## 2.9 | Provenance assignment

The provenance assignment success of the DArTcap markers was tested with *asssignPOP* v. 1.2.2 (Chen, Marschall, et al., 2018) and *rubias* v.0.3.2 (Anderson et al., 2008; Moran & Anderson, 2019). Assignment accuracy was tested with *assignPOP*, using both the Monte-Carlo and K-fold cross-validation procedures to test the assignment of a hold-out data set with 1000 iterations. We tested power of the markers by selecting a subset of loci with the highest $F_{ST}$ values (5%, 10%, 50%, and 100% of all loci) to train the assignment model. Similarly, the assignment accuracy of simulated mixed groups, based on a reference leave-one-out dataset, was evaluated with *rubias* (Anderson et al., 2008). Known simulated proportions for each reporting unit were compared with the numbers estimated by *rubias*. Populations with a sample size of one (i.e., Sierra Leone, eastern Atlantic) were excluded from these

analyses. We also examined the minimum number of informative markers needed to assign provenance by subsampling 5-500 markers based on loading contributions of each principal component from the DAPC analysis and testing the assignment accuracy with *rubias*.

## 2.10 | Kinship assignment

To investigate fine-scale contemporary connectivity, close-kin relationships were examined in each identified genetic cluster with >100 individuals to allow accurate allele frequency estimation, namely W-ATL and IWP. Kinship was tested using a log-likelihood-ratio (LLR) approach developed by Bravington et al. (2016) and applied by Hillary et al. (2018) and Feutry et al. (2020). A statistical threshold was set to reduce the number of false positive detections (i.e., more distantly-related kin) due to the large number of pairwise comparisons. Replicated or recaptured individuals were already visually identified with the 'filter-rad' function based on the number of loci with the same genotype for each pair of individuals (<10% genotypic difference).

# 3 | RESULTS

## 3.1 | SNP genotyping, baiting success, and data filtering

One sample from Papua New Guinea (PNG) failed DArTseq library construction, resulting in 187 samples for DArTseq sequencing and genotyping. An average of 2,182,162 reads (of 69 bp length) per sample were obtained from the sequencing and a total of 250,945 SNPs, located on 168,810 unique RAD-tags, were called by the DArTsoft14 program. Data filtering discarded 233,916 SNPs and 33 sharks, leaving a total of 17,029 high-quality SNPs. Of those, we randomly selected 3200 loci with at least one SNP (*Supporting Information* sections 4–6). We also identified 469 sex-linked markers (three Y-linked and 466 X-linked markers; *Supporting Information* section 4.2). Two Y- and 208 polymorphic X-linked sequences were included for DArTcap bait design.

Five samples failed DArTcap library construction. On average, 583,809 reads per sample were obtained from the DArTcap sequencing. After sequence clustering and SNP calling, we obtained 37,537 SNPs found on 26,335 RAD-tags. Quality-filtering was applied to six subsets of the DArTcap data (Table 1): the full data set (DATA1: 1014 sharks including replicates and other species; 5053 SNPs), the global dataset with confirmed Bull Sharks (DATA2: 769 sharks; 3409 SNPs), the global dataset with equalized samples sizes (DATA3: 382 sharks; 1849 SNPs), and regional Bull Shark data sets: W-ATL data (DATA4: 91 sharks; 931 SNPs), IWP data (DATA5: 635 sharks; 3416 SNPs), and IWP data with equalized samples size (DATA6: 189 sharks; 1785 SNPs).
## 3.2 | Mitogenome sequencing and assembly

The mitogenome was sequenced for 384 individuals with an average of 49,766 reads. Six individuals had low sequence coverage (<5000 mapped reads and mtDNA regions with no sequence coverage relative to the reference genome) and were subsequently removed from analyses (from South Africa, Sri Lanka, PNG, and Australia: Clarence River and Wenlock River). All reads that mapped to the reference genome were checked for ambiguous base calls at an 85% threshold (where a base needed to be present in >85% of mapped reads to be called unambiguous) to detect DNA cross-contamination, barcode slippage, or heteroplasmy. Seven individuals that had ambiguous base calls were removed due to signs of DNA contamination (originating from South Africa, PNG, Seychelles, Réunion, and Australia). Overall, the mitogenome length was 16,707–16,708 bp (Supporting Information section 13).

## 3.3 | Species identification

Mitogenome sequencing revealed 11 individuals from the Arabian Peninsula, Réunion, Sri Lanka, PNG, Costa Rica, and Fiji that were not Bull Shark: Pigeye Shark (*Carcharhinus amboinensis*, *n* = 5), Spinner Shark (*Carcharhinus brevipinna*, *n* = 1), Graceful Shark (*Carcharhinus amblyrhynchoides*, *n* = 1), Gray Reef Shark (*Carcharhinus amblyrhynchos*, *n* = 1), Pacific Smalltail Shark (*Carcharhinus cerdale*, *n* = 1), Dusky Shark (*Carcharhinus obscurus*, *n* = 1), and Speartooth Shark (*Glyphis glyphis*, *n* = 1). Furthermore, in the PCA analysis carried out on DATA1 an additional 54 individuals grouped with one of these species and another four individuals were neither these nor Bull Shark (Supporting Information section 7.8). Consequently, these individuals were omitted, and 769 confirmed Bull Sharks were retained in the filtered DArTcap data (DATA2) and 361 in the filtered mitogenome data.

## 3.4 | Global genetic diversity

More than half the DArTseq and DArTcap (DATA2) markers were monomorphic at almost all locations (Table 2 and Supporting Information section 5.6). In both genotyping protocols, the $H_O$ for the E-PAC ($H_O$ = 0.036 in Costa Rica) and the W-ATL ($H_O$ range = 0.048–0.052) were lower than those for the E-ATL and IWP ($H_O$ range = 0.056–0.069; Table 2). Most sampling locations exhibited significantly positive inbreeding coefficients (Table 2), but this was most pronounced in the E-PAC (Costa Rica) and IWP (South Africa, Mozambique, and Fiji; $F_{IS}$ = 0.041–0.068).

A total of 165 polymorphic sites were identified across all 361 mitogenome sequences, with an average nucleotide diversity of 0.001 and a haplotype diversity of 0.890 (Supporting Information section 13.3.1). Mitochondrial nucleotide diversities per sampling location ranged from 0.0001 in northern Australia (e.g., Roper and Towns rivers) to 0.003 in Sri Lanka (Table 2). Of the 165 polymorphic

sites, most were only present within the IWP group (Supporting Information section 13.3.4). Haplotype diversities were high in almost all sampling locations (h range = 0.80–0.99), with the lowest values in the Seychelles, Trinity Inlet (Australia), and Japan (h = 0.50–0.70).

## 3.5 | Population structure

Based on the nuclear DNA data ($F_{ST}$, DAPC, and PCA), we identified four major genetic clusters and additional hierarchical structure within these clusters (Figure 2, Table 3). The DAPC indicated that the optimal number of clusters for DATA3 (i.e., equal sample sizes) was K = 3–4 (Supporting Information section 9.9). These clusters corresponded to the E-PAC, W-ATL, E-ATL, and IWP/Japan/Fiji ($F_{ST}$ = 0.10–0.51, *p* < .001; Table 3). Further investigation within the latter group (i.e., DATA5; Supporting Information section 11.9) revealed further clustering into three groups corresponding to IWP, Japan, and Fiji (Figure 2). Here, the $F_{ST}$ values ranged from 0.06 between Japan and Fiji to <0.02 among IWP locations (*p* < .001). We saw high divergence between the two Japan locations (Okinawa and Urauchi River; Figure 2c; $F_{ST}$ = 0.03, *p* < .001), where individuals from Okinawa clustered closer to the IWP group and individuals from the Urauchi River were split into two groups. While the DAPC grouped the E-ATL (one individual from Sierra Leone) with the IWP cluster, $F_{ST}$ showed this location as distinctly different (Table 3). Nuclear differentiation between locations in the W-ATL were small ($F_{ST}$ < 0.007), although the DAPC and PCA analyses showed signs of differentiation between the southern (Brazil) and the northern locations (Gulf of Mexico and Western North Atlantic; Supporting Information section 10). Overall, these results indicated significant genetic differentiation of the E-PAC, southern W-ATL, northern W-ATL, E-ATL, IWP, Japan, and Fiji.

The mitogenome distance-based $\Phi_{ST}$ values showed high differentiation between and across ocean basins ($\Phi_{ST}$ = 0.67–0.99, *p* < .001; Table 3) and the haplotype network demonstrated four major clusters (Figure 3): E-PAC, E-ATL, western Indian Ocean (W-IO), and the E-IO and western Pacific combined (E-IO/W-PAC). All haplotypes, except four, were unique to a single sampling location (Figure 3). Within the Indian Ocean, differentiation was strong between the W-IO and E-IO locations ($\Phi_{ST}$ = 0.75–0.96, *p* < .001), while the northern Indian Ocean locations (N-IO; i.e., Arabian Peninsula, Sri Lanka, and Thailand) showed intermediate divergence ($\Phi_{ST}$ = 0.46–0.52, *p* < .001; Supporting Information section 13.5.3). The latter is demonstrated in the haplotype network with individuals from the N-IO split between the W-IO and E-IO/W-PAC haplogroups. Between the W-IO and E-IO/W-PAC groups, the haplotypes of several individuals indicated matrilineal movement at an evolutionary-recent time scale: two haplotypes, sampled from the Fitzroy River (Australia) and Fiji, clustered with the W-IO group and one haplotype from South Africa grouped with the E-IO/W-PAC group (indicated in Figure 3). Haplotypes sampled from Japan (*n* = 4) and Fiji (*n* = 8) grouped with the E-IO/W-PAC
<table><caption>TABLE 2 Description of the genetic variation for the full *Carcharhinus leucas* data set (DATA2: 769 sharks; 3409 SNPs) and the full mitochondrial genomes (361 sharks; 16,708 bp).</caption><thead><tr><th></th><th>GOC</th><th>COR</th><th>W-ATL</th><th>BRZ</th><th>CAR</th><th>GOM</th><th>WNA</th><th>E-ATL/S/IL</th><th>IWP</th><th>SAF</th><th>MOZ</th><th>RUN</th><th>SEY</th><th>ARP</th><th>SRL</th><th>TAI</th><th>IND</th><th>PNG</th></tr></thead><tbody><tr><td>N<sub>nu</sub></td><td>16</td><td>0</td><td>16</td><td>97</td><td>50</td><td>2</td><td>37</td><td>8</td><td>1</td><td>586</td><td>22</td><td>11</td><td>28</td><td>36</td><td>15</td><td>12</td><td>5</td><td>4</td><td>9</td></tr><tr><td>A<sub>r</sub></td><td>1.1540</td><td>NA</td><td>1.1540</td><td>1.2580</td><td>1.0950</td><td>1.0650</td><td>1.0870</td><td>1.0840</td><td>NA</td><td>1.3450</td><td>1.1120</td><td>1.0860</td><td>1.1080</td><td>1.1120</td><td>1.1140</td><td>1.1100</td><td>1.1090</td><td>1.0970</td><td>1.1050</td></tr><tr><td>H<sub>o</sub></td><td>0.0360</td><td>NA</td><td>0.0360</td><td>0.0500</td><td>0.0520</td><td>0.0490</td><td>0.0480</td><td>0.0500</td><td>NA</td><td>0.0610</td><td>0.0600</td><td>0.0560</td><td>0.0590</td><td>0.0600</td><td>0.0620</td><td>0.0610</td><td>0.0660</td><td>0.0610</td><td>0.0590</td></tr><tr><td>uH<sub>E</sub></td><td>0.0410</td><td>NA</td><td>0.0410</td><td>0.0550</td><td>0.0570</td><td>0.0510</td><td>0.0530</td><td>0.0540</td><td>NA</td><td>0.0640</td><td>0.0630</td><td>0.0640</td><td>0.0630</td><td>0.0630</td><td>0.0650</td><td>0.0640</td><td>0.0660</td><td>0.0610</td><td>0.0620</td></tr><tr><td>F<sub>is</sub></td><td>0.0750</td><td>NA</td><td>0.0750</td><td>0.0536</td><td>0.0576</td><td>-0.0481</td><td>0.0638</td><td>0.0369</td><td>NA</td><td>0.0472</td><td>0.0465</td><td>0.0959</td><td>0.0572</td><td>0.0398</td><td>0.0370</td><td>0.0356</td><td>-0.0069</td><td>-0.0098</td><td>0.0285</td></tr><tr><td>F<sub>is2.5</sub></td><td>0.0781</td><td>NA</td><td>0.0781</td><td>0.0748</td><td>0.0694</td><td>-0.0306</td><td>0.0724</td><td>0.0327</td><td>NA</td><td>0.0421</td><td>0.0467</td><td>0.0985</td><td>0.0528</td><td>0.0312</td><td>0.0273</td><td>0.0126</td><td>-0.0320</td><td>-0.0293</td><td>0.0192</td></tr><tr><td>F<sub>is7.5</sub></td><td>0.1401</td><td>NA</td><td>0.1401</td><td>0.1082</td><td>0.1060</td><td>0.1257</td><td>0.1165</td><td>0.0984</td><td>NA</td><td>0.0579</td><td>0.0811</td><td>0.1568</td><td>0.0893</td><td>0.0606</td><td>0.0668</td><td>0.0614</td><td>0.0378</td><td>0.0551</td><td>0.0703</td></tr><tr><td>M<sub>o</sub></td><td>1370</td><td>NA</td><td>2562</td><td>3340</td><td>1771</td><td>2811</td><td>2075</td><td>2512</td><td>3340</td><td>301</td><td>2081</td><td>2443</td><td>2010</td><td>1914</td><td>2195</td><td>2353</td><td>2699</td><td>2835</td><td>2512</td></tr><tr><td>N<sub>int</sub></td><td>16</td><td>1</td><td>15</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>302</td><td>7</td><td>0</td><td>16</td><td>18</td><td>19</td><td>12</td><td>6</td><td>6</td><td>15</td></tr><tr><td>π<sub>int</sub></td><td>0.0003</td><td>NA</td><td>0.0003</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.0021</td><td>0.0019</td><td>NA</td><td>0.0003</td><td>0.0003</td><td>0.0018</td><td>0.0030</td><td>0.0029</td><td>0.0004</td><td>0.0006</td></tr><tr><td>H</td><td>10</td><td>1</td><td>9</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>1</td><td>142</td><td>6</td><td>NA</td><td>6</td><td>3</td><td>13</td><td>11</td><td>4</td><td>6</td><td>13</td></tr><tr><td>h</td><td>0.9250</td><td>NA</td><td>0.9143</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.9873</td><td>0.9524</td><td>NA</td><td>0.8083</td><td>0.6601</td><td>0.9415</td><td>0.9848</td><td>0.8000</td><td>1.0000</td><td>0.9810</td></tr><tr><td>AUS</td><td>444</td><td>3</td><td>0</td><td>49</td><td>18</td><td>15</td><td>88</td><td>47</td><td>18</td><td>9</td><td>10</td><td>24</td><td>5</td><td>83</td><td>65</td><td>10</td><td>44</td><td>9</td><td>35</td><td>25</td></tr><tr><td>N<sub>nu</sub></td><td>1.1160</td><td>1.0830</td><td>NA</td><td>1.1140</td><td>1.1120</td><td>1.1160</td><td>1.1130</td><td>1.1110</td><td>1.1160</td><td>1.1000</td><td>1.1060</td><td>1.1120</td><td>1.1110</td><td>1.1200</td><td>1.1130</td><td>1.1080</td><td>1.2800</td><td>1.1020</td><td>1.1090</td><td>1.2820</td></tr><tr><td>A<sub>r</sub></td><td>0.0610</td><td>0.0600</td><td>NA</td><td>0.0600</td><td>0.0600</td><td>0.0640</td><td>0.0600</td><td>0.0590</td><td>0.0650</td><td>0.0590</td><td>0.0590</td><td>0.0600</td><td>0.0690</td><td>0.0650</td><td>0.0600</td><td>0.0610</td><td>0.0610</td><td>0.0590</td><td>0.0620</td><td>0.0570</td></tr><tr><td>H<sub>o</sub></td><td>0.0640</td><td>0.0620</td><td>NA</td><td>0.0630</td><td>0.0630</td><td>0.0670</td><td>0.0630</td><td>0.0620</td><td>0.0660</td><td>0.0650</td><td>0.0630</td><td>0.0640</td><td>0.0680</td><td>0.0660</td><td>0.0630</td><td>0.0640</td><td>0.0610</td><td>0.0600</td><td>0.0600</td><td>0.0640</td></tr><tr><td>F<sub>is</sub></td><td>0.0440</td><td>-0.0006</td><td>NA</td><td>0.0470</td><td>0.0392</td><td>0.0275</td><td>0.0409</td><td>0.0383</td><td>0.0212</td><td>0.0648</td><td>0.0434</td><td>0.0456</td><td>-</td><td>0.0301</td><td>0.0448</td><td>0.0342</td><td>0.0155</td><td>-0.0026</td><td>0.0006</td><td>0.0809</td></tr><tr><td>F<sub>is2.5</sub></td><td>0.0387</td><td>-0.0016</td><td>NA</td><td>0.0448</td><td>0.0338</td><td>0.0190</td><td>0.0370</td><td>0.0333</td><td>0.0034</td><td>0.0614</td><td>0.0449</td><td>0.0441</td><td>0.0430</td><td>0.0062</td><td>0.0431</td><td>0.0270</td><td>-0.0151</td><td>-0.0273</td><td>-0.0440</td><td>0.0955</td></tr><tr><td>F<sub>is7.5</sub></td><td>0.0544</td><td>0.0940</td><td>NA</td><td>0.0704</td><td>0.0695</td><td>0.0588</td><td>0.0596</td><td>0.0591</td><td>0.0410</td><td>0.1174</td><td>0.0948</td><td>0.0744</td><td>0.0249</td><td>0.0296</td><td>0.0676</td><td>0.0768</td><td>0.0170</td><td>0.0315</td><td>-0.0057</td><td>0.1358</td></tr><tr><td>M<sub>o</sub></td><td>395</td><td>2914</td><td>NA</td><td>1727</td><td>2178</td><td>2204</td><td>1360</td><td>1779</td><td>2120</td><td>2517</td><td>2480</td><td>2033</td><td>2681</td><td>1320</td><td>1550</td><td>2452</td><td>1969</td><td>2653</td><td>2133</td><td>2114</td></tr><tr><td>N<sub>int</sub></td><td>203</td><td>5</td><td>21</td><td>25</td><td>18</td><td>0</td><td>38</td><td>17</td><td>13</td><td>4</td><td>5</td><td>19</td><td>4</td><td>21</td><td>13</td><td>0</td><td>31</td><td>10</td><td>21</td><td>9</td></tr><tr><td>π<sub>int</sub></td><td>0.0006</td><td>0.0031</td><td>0.0004</td><td>0.0006</td><td>0.0007</td><td>NA</td><td>0.0002</td><td>0.0005</td><td>0.0001</td><td>0.0001</td><td>0.0001</td><td>0.0001</td><td>0.0001</td><td>0.0001</td><td>0.0003</td><td>0.0003</td><td>0.0003</td><td>NA</td><td>0.0004</td><td>0.0006</td><td>0.0001</td></tr><tr><td>H</td><td>85</td><td>4</td><td>14</td><td>16</td><td>12</td><td>NA</td><td>18</td><td>12</td><td>7</td><td>3</td><td>3</td><td>14</td><td>2</td><td>12</td><td>11</td><td>NA</td><td>4</td><td>3</td><td>3</td><td>0.0014</td></tr><tr><td>h</td><td>0.9777</td><td>0.9000</td><td>0.9333</td><td>0.9367</td><td>0.9216</td><td>NA</td><td>0.9346</td><td>0.9412</td><td>0.8718</td><td>0.8333</td><td>0.7000</td><td>0.9649</td><td>0.5000</td><td>0.9333</td><td>0.9744</td><td>NA</td><td>0.6903</td><td>0.7111</td><td>0.5286</td><td>0.9722</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>FIGURE 2 Population clustering analysis for the global *Carcharhinus leucas* dataset with a subsample of Australian sharks (DATA3: 382 sharks; 1849 SNPs). Panel a shows the Discriminant Analysis of Principal Components (DAPC) assignment barplot for K = 6 and 51 principal components (PC). Each bar represents an individual and is colored according to the posterior membership probabilities. Panels b–d represent the principal component analysis (PCA) scatterplot, where each point represents an individual shark. triangles indicate the mean PCA score per sampling location, and colors represent the sampling country or oceanographic location. (b) PCA scatterplot with PC1 on the x-axis and PC2 on the y-axis. (c) PCA scatterplot with PC3 on the x-axis and PC4 on the y-axis. (d) PCA scatterplot with PC3 on the x-axis and PC5 on the y-axis.</figcaption></figure>
cluster, where the Japanese haplotypes did not correspond to the three Japanese clusters identified with the nuclear SNPs. Within Australia, the mitochondrial $\Phi_{ST}$ showed structure between western (Fitzroy River), northern (Victoria River, Daly River, Adelaide River, Alligator Rivers—encompassing the South Alligator River and East Alligator River—Blue Mud Bay, Roper River, Towns River, Wenlock River), and eastern Australian sites (Trinity Inlet, Clarence River, and Sydney Harbor; Table 3). The haplotype network of all Australian individuals indicated three major haplotype clusters with no obvious geographic pattern, yet most haplotypes were unique to a sampling site (Supporting Information section 13.4.2.3).

## 3.6 | Provenance assignment

After population structure was identified with the nuclear DArTcap markers, we tested the accuracy to assign provenance to a hold-out data set (Figure 4). Using the marker contributions of the DAPC analysis, we found that at least 250–500 highly differentiating markers were required for 100% provenance assignment accuracy (Supporting Information section 8.10.9). Individuals from Japan and Fiji were unlikely to be assigned to the inferred genetic cluster with less than 500 markers. A minimum of 50–100 markers were needed to obtain a reliable assignment (>80%) to each inferred population and only 5–50 markers were required to differentiate individuals from the E-PAC, W-PAC, and E-ATL, but sites within the IWP required up to 100 markers (Figure 4).

## 3.7 | Individual and kinship assignment

In total, 102 duplicate individuals were identified in the dataset. Most represented the technical replicates that were purposely included, but we also detected six recaptured sharks from Sydney Harbor, Japan, and Fiji (time-at-liberty = 38–741 days). These recaptures were always assigned to the same age-cohort where length data were available, thus providing confidence in the length-at-age function from Tillett et al. (2011). Within the W-ATL, we performed 4095 pairwise comparisons and found two full-sibling pairs (FSPs) in Brazil. The IWP revealed more kin pairs due to the large number of pairwise comparisons (201,295). Overall, we found 18 full-sibling pairs (FSPs), all within the same sampling locations (Brazil, Réunion, Indonesia, Fiji, and within Australia: Daly, South Alligator, East Alligator, Towns, and Clarence rivers), and 40 half-sibling pairs (HSPs). Of the 40 HSPs, 32 were found within the same locations (Seychelles, Réunion, PNG, and within Australia: Sydney Harbor, South Alligator, East Alligator, Wenlock, and Clarence rivers), where 10 and 8 HSPs had the same and different haplotypes, respectively; 14 pairs had missing haplotype information. Eight HSPs were distributed between rivers within northern Australia (Daly River/South Alligator River, *n* = 4; Adelaide River/South Alligator River, *n* = 1; South Alligator River/East Alligator River, *n* = 2; and Towns River/Wenlock River, *n* = 1). Seven of those eight ‘cross-river’ HSPs were

juveniles (<150 cm TL) from different age cohorts, and six pairs had different mitochondrial haplotypes.

## 4 | DISCUSSION

This study provides the first global assessment of genetic population structure of the Bull Shark, a cosmopolitan coastal predator, using thousands of SNP markers and full mitogenome data. We identify distinct genetic divergence driven by significant biogeographic barriers and philopatric behavior. These results further facilitate the development of a SNP panel for species, sex, provenance, and kinship identification.

### 4.1 | Global population structure across ocean barriers

Population structure in large coastal sharks tends to be driven by environmental barriers, movement ecology, and habitat preference (Dudgeon et al., 2012; Hirschfeld et al., 2021). The DArTcap results revealed four global genetic clusters (E-PAC, W-ATL, E-ATL, and IWP), with Japan and Fiji each supporting additional isolated Bull Shark populations. Our results concur with previous studies in the Western Atlantic and Indo-Pacific regions, which showed genetic differentiation between the W-ATL, IWP, and Fiji using microsatellite data (Pirog et al., 2019; Testerman, 2014). We show that the E-PAC is strongly differentiated from the W-ATL, likely coinciding with the closure of the Isthmus of Panama (~3 Myr; Knowlton et al., 1993; O'Dea et al., 2016). Previously, Bull Shark connectivity has been suggested between the E-PAC and W-ATL, based on microsatellite data, through the Panama Canal (Pirog et al., 2019; Testerman, 2014), but our results indicate these populations are demographically isolated and that gene flow through the Panama Canal is negligible, if existent at all. Rather, we believe the differences between our results are likely due to higher power delivered by thousands of SNPs compared to fewer microsatellites (e.g., Green et al., 2019; Layton et al., 2020). In addition, the latter can suffer from homoplasy caused by high mutation rates (i.e., identical by state, not identical by descent; Estoup et al., 2002).

The greater divergence between the W-ATL and E-ATL, compared to the E-ATL and IWP, demonstrates that the Mid Atlantic Barrier (an open ocean distance barrier) forms a stronger and more enduring barrier than the Benguela Upwelling System (a thermal barrier). The permeability of the latter barrier may fluctuate over time due to oscillations in climate (Hirschfeld et al., 2021). Fiji and New Caledonia have previously been identified as differentiated island populations, although New Caledonia showed less differentiation than Fiji from other sites (Glaus et al., 2020; Pirog et al., 2019). This was also attributed to the long distances across deep water trenches and is consistent with our results for Fiji. Similarly, individuals from Japan could be isolated by deep water trenches (>500 m deep, e.g., Okinawa
<table><caption>TABLE 3 Pairwise comparison of fixation indices for *Carcharhinus leucas*.</caption><thead><tr><th rowspan="2"></th><th colspan="2">E-PAC</th><th colspan="4">W-ATL</th><th colspan="2">E-ATL</th><th colspan="2">IWP</th><th rowspan="2">RUN</th><th rowspan="2">SEY</th><th rowspan="2">ARP</th><th rowspan="2">SRL</th><th rowspan="2">TAI</th><th rowspan="2">IND</th><th rowspan="2">PNG</th></tr><tr><th>GOC</th><th>COR</th><th>BRZ</th><th>CAR</th><th>GOM</th><th>WNA</th><th>SIL</th><th>SAF</th><th>MOZ</th><th>SEY</th></tr></thead><tbody><tr><td>GOC</td><td></td><td>0.124</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.788</td><td>NA</td><td>0.968</td><td>*0.967</td><td>0.804</td><td>0.664</td><td>0.677</td><td>0.951</td><td>0.931</td></tr><tr><td>COR</td><td>NA</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.958</td><td>*0.913</td><td>NA</td><td>*0.969</td><td>*0.968</td><td>*0.877</td><td>*0.832</td><td>*0.887</td><td>*0.965</td><td>*0.951</td></tr><tr><td>BRZ</td><td>*0.220</td><td>*0.237</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>CAR</td><td>*0.229</td><td>*0.248</td><td>*0.007</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>GOM</td><td>*0.225</td><td>*0.243</td><td>*0.004</td><td>*0.006</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>WNA</td><td>*0.229</td><td>*0.251</td><td>*0.005</td><td>*0.003</td><td>*0.004</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>SIL</td><td>NA</td><td>NA</td><td>0.467</td><td>0.473</td><td>0.473</td><td>0.476</td><td></td><td>0.788</td><td>NA</td><td>0.968</td><td>0.967</td><td>*0.802</td><td>0.662</td><td>0.676</td><td>0.950</td><td>0.930</td></tr><tr><td>SAF</td><td>*0.478</td><td>*0.497</td><td>*0.472</td><td>*0.473</td><td>*0.476</td><td>*0.479</td><td>0.098</td><td></td><td>NA</td><td>0.211</td><td>*0.159</td><td>*0.590</td><td>0.113</td><td>0.048</td><td>0.747</td><td>*0.782</td></tr><tr><td>MOZ</td><td>*0.489</td><td>*0.508</td><td>*0.481</td><td>*0.482</td><td>*0.485</td><td>*0.488</td><td>0.107</td><td>*0.004</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>RUN</td><td>*0.485</td><td>*0.502</td><td>*0.476</td><td>*0.478</td><td>*0.481</td><td>*0.483</td><td>0.104</td><td>*0.002</td><td>*0.006</td><td>*0.606</td><td>*0.776</td><td>*0.456</td><td>*0.520</td><td>*0.941</td><td>*0.919</td></tr><tr><td>SEY</td><td>*0.486</td><td>*0.502</td><td>*0.476</td><td>*0.478</td><td>*0.481</td><td>*0.483</td><td>0.107</td><td>*0.002</td><td>*0.006</td><td>*0.003</td><td>*0.776</td><td>*0.395</td><td>*0.399</td><td>*0.937</td><td>*0.916</td></tr><tr><td>ARP</td><td>*0.479</td><td>*0.496</td><td>*0.470</td><td>*0.472</td><td>*0.475</td><td>*0.477</td><td>0.102</td><td>*0.002</td><td>*0.005</td><td>*0.002</td><td>*0.002</td><td></td><td>*0.292</td><td>*0.387</td><td>0.135</td><td>*0.286</td></tr><tr><td>SRL</td><td>*0.492</td><td>*0.507</td><td>*0.480</td><td>*0.481</td><td>*0.484</td><td>*0.487</td><td>0.107</td><td>*0.003</td><td>*0.006</td><td>*0.003</td><td>*0.005</td><td>0.001</td><td></td><td>0.000</td><td>*0.391</td><td>*0.504</td></tr><tr><td>TAI</td><td>0.492</td><td>*0.509</td><td>*0.480</td><td>*0.482</td><td>*0.484</td><td>*0.487</td><td>0.122</td><td>*0.011</td><td>*0.013</td><td>*0.010</td><td>*0.012</td><td>*0.009</td><td>*0.011</td><td></td><td>*0.550</td><td>*0.648</td></tr><tr><td>IND</td><td>0.478</td><td>*0.497</td><td>*0.468</td><td>*0.470</td><td>*0.473</td><td>*0.475</td><td>0.102</td><td>*0.009</td><td>*0.013</td><td>*0.012</td><td>*0.014</td><td>*0.012</td><td>*0.014</td><td>*0.020</td><td></td><td>*0.424</td></tr><tr><td>PNG</td><td>*0.483</td><td>*0.497</td><td>*0.471</td><td>*0.472</td><td>*0.475</td><td>*0.478</td><td>0.108</td><td>*0.004</td><td>*0.012</td><td>*0.005</td><td>*0.006</td><td>*0.006</td><td>*0.006</td><td>*0.013</td><td>*0.015</td><td></td></tr><tr><td>FZR</td><td>NA</td><td>*0.512</td><td>0.483</td><td>0.483</td><td>0.487</td><td>0.489</td><td>NA</td><td>0.002</td><td>0.006</td><td>0.003</td><td>*0.007</td><td>0.005</td><td>*0.009</td><td>0.016</td><td>0.017</td><td>0.005</td></tr><tr><td>VIR</td><td>*0.483</td><td>*0.497</td><td>*0.472</td><td>*0.473</td><td>*0.476</td><td>*0.479</td><td>0.105</td><td>0.001</td><td>*0.005</td><td>*0.004</td><td>*0.004</td><td>*0.003</td><td>*0.004</td><td>*0.011</td><td>*0.014</td><td>*0.006</td></tr><tr><td>DAR</td><td>*0.485</td><td>*0.500</td><td>*0.476</td><td>*0.477</td><td>*0.480</td><td>*0.483</td><td>0.103</td><td>*0.001</td><td>*0.006</td><td>*0.002</td><td>*0.002</td><td>*0.002</td><td>*0.003</td><td>*0.011</td><td>*0.012</td><td>*0.003</td></tr><tr><td>ADR</td><td>*0.485</td><td>*0.501</td><td>*0.475</td><td>*0.477</td><td>*0.479</td><td>*0.482</td><td>0.108</td><td>*0.002</td><td>*0.005</td><td>*0.002</td><td>*0.002</td><td>0.001</td><td>*0.003</td><td>*0.011</td><td>*0.009</td><td>*0.004</td></tr><tr><td>DWC</td><td>*0.475</td><td>*0.493</td><td>*0.468</td><td>*0.469</td><td>*0.472</td><td>*0.475</td><td>0.099</td><td>0.001</td><td>*0.006</td><td>*0.004</td><td>*0.003</td><td>*0.003</td><td>*0.005</td><td>*0.012</td><td>*0.011</td><td>*0.005</td></tr><tr><td>SAR</td><td>*0.485</td><td>*0.500</td><td>*0.475</td><td>*0.476</td><td>*0.479</td><td>*0.482</td><td>0.105</td><td>0.001</td><td>*0.005</td><td>*0.003</td><td>*0.002</td><td>*0.002</td><td>*0.003</td><td>*0.009</td><td>*0.012</td><td>*0.004</td></tr><tr><td>EAR</td><td>*0.483</td><td>*0.499</td><td>*0.475</td><td>*0.476</td><td>*0.479</td><td>*0.482</td><td>0.104</td><td>*0.001</td><td>*0.006</td><td>*0.003</td><td>*0.002</td><td>*0.002</td><td>*0.002</td><td>*0.011</td><td>*0.013</td><td>*0.004</td></tr><tr><td>BMB</td><td>*0.486</td><td>*0.502</td><td>*0.476</td><td>*0.477</td><td>*0.480</td><td>*0.483</td><td>0.109</td><td>*0.005</td><td>*0.008</td><td>*0.005</td><td>*0.006</td><td>*0.006</td><td>*0.007</td><td>*0.013</td><td>*0.016</td><td>*0.008</td></tr><tr><td>ROR</td><td>*0.481</td><td>*0.502</td><td>*0.475</td><td>*0.476</td><td>*0.479</td><td>*0.482</td><td>0.098</td><td>*0.002</td><td>0.006</td><td>*0.004</td><td>*0.004</td><td>*0.003</td><td>0.003</td><td>*0.010</td><td>*0.013</td><td>*0.006</td></tr><tr><td>TOR</td><td>*0.495</td><td>*0.508</td><td>*0.481</td><td>*0.483</td><td>*0.485</td><td>*0.489</td><td>0.114</td><td>*0.006</td><td>*0.011</td><td>*0.007</td><td>*0.006</td><td>*0.004</td><td>*0.008</td><td>*0.015</td><td>*0.019</td><td>*0.010</td></tr><tr><td>WER</td><td>*0.488</td><td>*0.503</td><td>*0.477</td><td>*0.478</td><td>*0.481</td><td>*0.484</td><td>0.111</td><td>*0.003</td><td>*0.008</td><td>*0.002</td><td>*0.003</td><td>*0.002</td><td>*0.003</td><td>*0.012</td><td>*0.014</td><td>*0.005</td></tr><tr><td>TRI</td><td>0.478</td><td>*0.498</td><td>*0.470</td><td>*0.472</td><td>*0.475</td><td>*0.477</td><td>0.103</td><td>*0.004</td><td>0.010</td><td>*0.007</td><td>*0.006</td><td>*0.005</td><td>*0.011</td><td>0.013</td><td>0.013</td><td>0.008</td></tr><tr><td>CLR</td><td>*0.473</td><td>*0.488</td><td>*0.465</td><td>*0.466</td><td>*0.469</td><td>*0.472</td><td>0.103</td><td>*0.002</td><td>*0.005</td><td>*0.003</td><td>*0.003</td><td>*0.002</td><td>*0.004</td><td>*0.010</td><td>*0.012</td><td>*0.005</td></tr><tr><td>SYH</td><td>*0.485</td><td>*0.500</td><td>*0.475</td><td>*0.476</td><td>*0.479</td><td>*0.482</td><td>0.104</td><td>*0.002</td><td>*0.005</td><td>*0.003</td><td>*0.003</td><td>*0.002</td><td>*0.003</td><td>*0.011</td><td>*0.013</td><td>*0.004</td></tr><tr><td>UNK</td><td>*0.484</td><td>*0.499</td><td>*0.473</td><td>*0.474</td><td>*0.477</td><td>*0.480</td><td>0.104</td><td>*0.003</td><td>*0.008</td><td>*0.005</td><td>*0.004</td><td>*0.005</td><td>*0.006</td><td>*0.015</td><td>*0.013</td><td>*0.009</td></tr><tr><td>OKI</td><td>*0.507</td><td>*0.516</td><td>*0.486</td><td>*0.488</td><td>*0.490</td><td>*0.494</td><td>0.125</td><td>*0.017</td><td>*0.024</td><td>*0.017</td><td>*0.018</td><td>*0.018</td><td>*0.023</td><td>*0.028</td><td>*0.023</td><td>*0.022</td></tr><tr><td>URR</td><td>*0.509</td><td>*0.519</td><td>*0.492</td><td>*0.493</td><td>*0.496</td><td>*0.499</td><td>0.142</td><td>*0.021</td><td>*0.026</td><td>*0.023</td><td>*0.020</td><td>*0.021</td><td>*0.025</td><td>*0.033</td><td>*0.034</td><td>*0.024</td></tr><tr><td>FIJ</td><td>*0.499</td><td>*0.514</td><td>*0.486</td><td>*0.487</td><td>*0.491</td><td>*0.493</td><td>0.116</td><td>*0.018</td><td>*0.023</td><td>*0.017</td><td>*0.017</td><td>*0.016</td><td>*0.019</td><td>*0.027</td><td>*0.026</td><td>*0.019</td></tr></tbody></table>

Note: $F_{ST}$ values are presented in the bottom diagonal and sequence-based $\Phi_{ST}$ are in the top diagonal. 'NA' values indicate that no or too few samples were available to estimate genetic differentiation. Asterisks (*) indicate statistical significance of $p$ < .05. No pairwise comparisons were significant after Bonferroni correction. Fixation indices are colored from high (red) to medium (yellow) to low (blue) values. GOC, Gulf of California; COR, Costa Rica; BRZ, Brazil; CAR, Caribbean Sea; GOM, Gulf of Mexico; WNA, Western North Atlantic; SIL, Sierra Leone; SAF, South Africa; MOZ, Mozambique; RUN, Réunion Island; SEY, Seychelles; ARP, Arabian Peninsula; SRL, Sri Lanka; TAI, Thailand; IND, Indonesia; PNG, Papua New Guinea; FZR, Fitzroy River; VIR, Victoria River; DAR, Daly River; ADR, Adelaide River; DWC, Darwin Coastal; SAR, South Alligator River; EAR, East Alligator River; BMB, Blue Mud Bay; ROR, Roper River; TOR, Towns River; WER, Wenlock River; TRI, Trinity Inlet; CLR, Clarence River; SYH, Sydney Harbor; and UNK, Australian fisheries samples from unknown origin; OKI, Churaumi Aquarium, Okinawa; URR, Urauchi River; and FIJ, Fiji.

Trough), since the Bull Shark generally prefers shallow waters (<160 m; Rigby et al., 2021). Yet, sea levels between Japan and continental Asia were low until 17 Kyr ago (Voris, 2000), although the high latitudinal location of Japan may have formed a thermal barrier at that time. Additionally, historical land-bridges such as the Taiwan Strait, Tokara Strait, or Tsushima/Korea Strait could have restricted dispersal between Japan and the IWP (Chen, Wang, et al., 2018; Yin et al., 2009).

Historical land-bridges also play a potential role for the genetic structuring of Bull Sharks across northern Australia. Despite a lack

of structure inferred from the nuclear markers, the mitogenome $\Phi_{ST}$ differentiation showed three separate clusters (western, northern, and eastern Australia). Here, the Torres Strait land-bridge between Australia and Papua New Guinea has been implicated for genetic structuring in a variety of marine taxa across this region (Mirams et al., 2011). This signal from the Bull Shark haplotype network is less structured, and may indicate either recent separation or historical isolation followed by gene flow (e.g., Puckridge et al., 2013).
<table><thead><tr><th colspan="15">Australia</th><th colspan="3">Japan</th><th colspan="3">Fiji</th></tr><tr><th>FZR</th><th>VIR</th><th>DAR</th><th>ADR</th><th>DWC</th><th>SAR</th><th>EAR</th><th>BMB</th><th>ROR</th><th>TOR</th><th>WER</th><th>TRI</th><th>CLR</th><th>SYH</th><th>UNK</th><th>OKI</th><th>URR</th><th>FIJ</th></tr></thead><tbody><tr><td>0.649</td><td>*0.951</td><td>*0.934</td><td>0.923</td><td>NA</td><td>*0.977</td><td>0.944</td><td>*0.988</td><td>0.988</td><td>0.985</td><td>*0.932</td><td>1.000</td><td>*0.966</td><td>0.936</td><td>NA</td><td>0.937</td><td>0.989</td><td>0.848</td></tr><tr><td>*0.891</td><td>*0.959</td><td>*0.947</td><td>*0.944</td><td>NA</td><td>*0.975</td><td>*0.956</td><td>*0.979</td><td>*0.973</td><td>*0.974</td><td>*0.949</td><td>*0.975</td><td>*0.968</td><td>*0.954</td><td>NA</td><td>*0.957</td><td>0.981</td><td>*0.925</td></tr><tr><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>0.646</td><td>*0.950</td><td>*0.933</td><td>0.922</td><td>NA</td><td>*0.977</td><td>0.943</td><td>0.988</td><td>0.988</td><td>0.985</td><td>0.931</td><td>1.000</td><td>0.966</td><td>0.935</td><td>NA</td><td>0.937</td><td>*0.989</td><td>0.846</td></tr><tr><td>0.276</td><td>*0.828</td><td>*0.813</td><td>*0.784</td><td>NA</td><td>*0.897</td><td>*0.808</td><td>*0.846</td><td>*0.736</td><td>*0.757</td><td>*0.800</td><td>*0.741</td><td>*0.849</td><td>*0.781</td><td>NA</td><td>*0.777</td><td>*0.884</td><td>*0.643</td></tr><tr><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>*0.713</td><td>*0.933</td><td>*0.915</td><td>*0.911</td><td>NA</td><td>*0.958</td><td>*0.929</td><td>*0.963</td><td>*0.953</td><td>*0.954</td><td>*0.918</td><td>*0.956</td><td>*0.946</td><td>*0.925</td><td>NA</td><td>*0.931</td><td>*0.967</td><td>*0.864</td></tr><tr><td>*0.678</td><td>*0.929</td><td>*0.912</td><td>*0.907</td><td>NA</td><td>*0.955</td><td>*0.925</td><td>*0.959</td><td>*0.949</td><td>*0.950</td><td>*0.915</td><td>*0.952</td><td>*0.943</td><td>*0.922</td><td>NA</td><td>*0.928</td><td>*0.964</td><td>*0.862</td></tr><tr><td>*0.203</td><td>*0.322</td><td>*0.252</td><td>*0.224</td><td>NA</td><td>*0.483</td><td>*0.256</td><td>*0.392</td><td>0.256</td><td>*0.287</td><td>*0.227</td><td>0.228</td><td>*0.362</td><td>*0.272</td><td>NA</td><td>*0.256</td><td>*0.326</td><td>*0.219</td></tr><tr><td>0.000</td><td>*0.563</td><td>*0.552</td><td>*0.500</td><td>NA</td><td>*0.692</td><td>*0.524</td><td>*0.561</td><td>*0.404</td><td>*0.434</td><td>*0.520</td><td>*0.393</td><td>*0.585</td><td>*0.490</td><td>NA</td><td>*0.476</td><td>*0.610</td><td>*0.334</td></tr><tr><td>0.030</td><td>*0.713</td><td>*0.691</td><td>*0.646</td><td>NA</td><td>*0.825</td><td>*0.677</td><td>*0.729</td><td>*0.544</td><td>*0.580</td><td>*0.666</td><td>*0.539</td><td>*0.741</td><td>*0.638</td><td>NA</td><td>*0.621</td><td>*0.780</td><td>*0.451</td></tr><tr><td>*0.345</td><td>*0.463</td><td>*0.320</td><td>*0.255</td><td>NA</td><td>*0.705</td><td>*0.389</td><td>*0.762</td><td>*0.651</td><td>*0.673</td><td>0.308</td><td>*0.629</td><td>*0.571</td><td>*0.409</td><td>NA</td><td>*0.425</td><td>*0.660</td><td>*0.321</td></tr><tr><td>*0.405</td><td>*0.290</td><td>*0.245</td><td>*0.153</td><td>NA</td><td>*0.486</td><td>*0.277</td><td>*0.475</td><td>*0.343</td><td>*0.380</td><td>*0.257</td><td>*0.429</td><td>*0.427</td><td>*0.116</td><td>NA</td><td>*0.411</td><td>*0.639</td><td>*0.339</td></tr><tr><td>*0.439</td><td>*0.433</td><td>*0.375</td><td>NA</td><td>NA</td><td>*0.611</td><td>*0.404</td><td>*0.469</td><td>0.197</td><td>*0.256</td><td>0.405</td><td>0.288</td><td>*0.543</td><td>*0.399</td><td>NA</td><td>*0.431</td><td>*0.666</td><td>*0.210</td></tr><tr><td>0.004</td><td>*0.063</td><td>0.039</td><td>NA</td><td>NA</td><td>*0.042</td><td>0.030</td><td>*0.133</td><td>0.009</td><td>0.073</td><td>*0.069</td><td>*0.502</td><td>*0.440</td><td>*0.317</td><td>NA</td><td>*0.501</td><td>*0.667</td><td>*0.410</td></tr><tr><td>*0.006</td><td>*0.002</td><td>0.025</td><td>NA</td><td>NA</td><td>*0.206</td><td>0.005</td><td>*0.268</td><td>0.145</td><td>*0.194</td><td>*0.029</td><td>*0.377</td><td>*0.366</td><td>*0.259</td><td>NA</td><td>*0.386</td><td>*0.514</td><td>*0.358</td></tr><tr><td>0.001</td><td>*0.003</td><td>*0.001</td><td>NA</td><td>NA</td><td>*0.223</td><td>0.028</td><td>*0.234</td><td>0.110</td><td>*0.156</td><td>*0.018</td><td>*0.335</td><td>*0.338</td><td>*0.184</td><td>NA</td><td>*0.356</td><td>*0.514</td><td>*0.312</td></tr><tr><td>0.003</td><td>0.001</td><td>*0.002</td><td>*0.002</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>*0.004</td><td>*0.002</td><td>*0.001</td><td>*0.001</td><td>*0.002</td><td>*0.148</td><td>*0.147</td><td>0.062</td><td>*0.148</td><td>*0.233</td><td>*0.720</td><td>*0.603</td><td>*0.527</td><td>NA</td><td>*0.687</td><td>*0.802</td><td>*0.594</td></tr><tr><td>*0.004</td><td>*0.002</td><td>*0.001</td><td>0.001</td><td>*0.002</td><td>0.000</td><td>*0.245</td><td>0.103</td><td>0.163</td><td>*0.044</td><td>*0.449</td><td>*0.411</td><td>*0.294</td><td>NA</td><td>*0.435</td><td>*0.615</td><td>*0.363</td></tr><tr><td>*0.009</td><td>*0.006</td><td>*0.005</td><td>*0.004</td><td>*0.006</td><td>*0.005</td><td>*0.005</td><td>*0.296</td><td>*0.331</td><td>*0.236</td><td>*0.873</td><td>*0.646</td><td>*0.528</td><td>NA</td><td>*0.691</td><td>*0.883</td><td>*0.512</td></tr><tr><td>*0.006</td><td>*0.004</td><td>*0.002</td><td>0.003</td><td>0.003</td><td>*0.002</td><td>0.002</td><td>*0.006</td><td>0.000</td><td>*0.132</td><td>*0.911</td><td>*0.578</td><td>*0.394</td><td>NA</td><td>*0.574</td><td>*0.880</td><td>*0.326</td></tr><tr><td>0.006</td><td>*0.007</td><td>*0.006</td><td>*0.005</td><td>*0.005</td><td>*0.005</td><td>*0.005</td><td>*0.008</td><td>0.005</td><td>*0.175</td><td>*0.885</td><td>*0.597</td><td>*0.428</td><td>NA</td><td>*0.598</td><td>*0.879</td><td>*0.367</td></tr><tr><td>0.004</td><td>*0.004</td><td>*0.002</td><td>*0.002</td><td>0.002</td><td>*0.001</td><td>*0.001</td><td>*0.006</td><td>*0.004</td><td>*0.006</td><td>*0.357</td><td>*0.348</td><td>*0.255</td><td>NA</td><td>*0.378</td><td>*0.529</td><td>*0.328</td></tr><tr><td>0.002</td><td>0.006</td><td>0.002</td><td>*0.006</td><td>0.003</td><td>*0.004</td><td>0.002</td><td>0.007</td><td>0.002</td><td>0.010</td><td>0.004</td><td>0.000</td><td>0.186</td><td>NA</td><td>*0.551</td><td>*0.882</td><td>0.000</td></tr><tr><td>*0.004</td><td>*0.003</td><td>*0.002</td><td>*0.001</td><td>*0.002</td><td>*0.001</td><td>*0.001</td><td>*0.005</td><td>*0.003</td><td>*0.005</td><td>*0.002</td><td>*0.004</td><td>*0.188</td><td>NA</td><td>*0.570</td><td>*0.750</td><td>*0.178</td></tr><tr><td>0.003</td><td>*0.002</td><td>*0.001</td><td>0.000</td><td>0.002</td><td>*0.001</td><td>*0.001</td><td>*0.005</td><td>*0.003</td><td>*0.005</td><td>*0.002</td><td>*0.004</td><td>*0.001</td><td>NA</td><td>*0.446</td><td>*0.666</td><td>*0.178</td></tr><tr><td>0.005</td><td>*0.005</td><td>*0.003</td><td>0.003</td><td>0.003</td><td>*0.002</td><td>*0.002</td><td>*0.007</td><td>*0.004</td><td>0.005</td><td>*0.004</td><td>0.003</td><td>*0.002</td><td>*0.002</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>*0.018</td><td>*0.019</td><td>*0.016</td><td>*0.017</td><td>*0.018</td><td>*0.017</td><td>*0.017</td><td>*0.021</td><td>*0.019</td><td>*0.023</td><td>*0.018</td><td>*0.019</td><td>*0.016</td><td>*0.017</td><td>*0.016</td><td>*0.555</td><td>*0.611</td><td></td></tr><tr><td>*0.023</td><td>*0.024</td><td>*0.022</td><td>*0.021</td><td>*0.021</td><td>*0.021</td><td>*0.021</td><td>*0.027</td><td>*0.024</td><td>*0.028</td><td>*0.021</td><td>*0.024</td><td>*0.021</td><td>*0.022</td><td>*0.021</td><td>*0.032</td><td>*0.389</td><td></td></tr><tr><td>*0.017</td><td>*0.018</td><td>*0.018</td><td>*0.017</td><td>*0.019</td><td>*0.017</td><td>*0.016</td><td>*0.020</td><td>*0.019</td><td>*0.021</td><td>*0.018</td><td>*0.022</td><td>*0.016</td><td>*0.017</td><td>*0.018</td><td>*0.032</td><td>*0.037</td><td></td></tr></tbody></table>

Our study shows strong influences of land-bridges, long ocean distances, cold upwelling or high latitudes, and deep-water trenches on the genetic population structure of the Bull Shark. We hypothesize that these barriers will equally affect elasmobranch species with similar dispersal capacity, distribution, and habitat use. To date, the global population structure of only a few large-bodied shark species has been investigated (e.g., Galapagos Shark, *C. galapagensis*; Sandbar Shark, *C. plumbeus*; and Scalloped Hammerhead, *S. lewini*; Daly-Engel et al., 2012; Pazmiño et al., 2018; Portnoy et al., 2010).

Known biogeographic barriers have similarly affected the population structure of these species (reviewed in Hirschfeld et al., 2021). For example, the Galapagos Shark, which occupies insular habitats, shows limited gene flow across oceanic barriers, such as the East Pacific Barrier (Pazmiño et al., 2018). Similarly, both the Scalloped Hammerhead and Sandbar Shark exhibited increased connectivity along continuous coastlines within ocean basins but showed significant divergence due to the East Pacific Barrier and closure of the Isthmus of Panama (Daly-Engel et al., 2012; Portnoy et al., 2010).
<figure><img src="image_6.png" /><figcaption>FIGURE 3 *Carcharhinus leucas* mtDNA haplotype network, based on the full mitogenome (16,707-16,708 bp). Panel a presents the network for all 361 sharks from the eastern Pacific, eastern Atlantic, Indo-West Pacific, Japan, and Fiji. The distance between haplotypes reflects the number of mutations between them. Panel b shows the 'eastern Indian Ocean/western Pacific/Japan/Fiji' cluster in detail (285 sharks). Here, the number of mutations between haplotypes are represented by small black dots. In panels a and b, the size of the shape is equivalent to the square root of the number of individuals that share this haplotype. The color and shape of each haplotype corresponds to the sampling location where they were found. The three black arrows indicate haplotypes that represent recent maternal movement between haplotype clusters.</figcaption></figure>

Consequently, future studies could use our results as a foundation to construct hypotheses for other exploited and threatened species with similar vagility, ecology, and life-history characteristics (such as the Spinner Shark, the Graceful Shark, or the Dusky Shark).

## 4.2 | Fine-scale structure and sex-biased dispersal

Minor nuclear differentiation was detected between the southern and northern W-ATL locations. This concurs with other Bull Shark studies showing strong mitochondrial and weak nuclear population structure in the W-ATL (Karl et al., 2011; Sandoval Laurrabaquio-Alvarado et al., 2021) and that long-distance movement is rare in the Gulf of Mexico (Carlson et al., 2010). This would suggest that a weak

or recent barrier to dispersal may be present in this region. Such a barrier can be attributed to oceanographic features (e.g., Caribbean Current) or reproductive asynchrony caused by seasonal differences between the hemispheres (Carrillo et al., 2015, 2017; Castro, 2011; Jensen, 1976). Similarly, genetic differentiation may be present in the IWP (as suggested by Pirog et al., 2019), but our ability to detect significant nuclear differences is constrained by high gene flow and/or recent separation (Bailleul et al., 2018; Waples, 1998). A different set of genetic markers (e.g., whole genome analysis) combined with alternative approaches, like telemetry or parasitology, may provide a new perspective on population structure within the IWP (reviewed in Green, Simpfendorfer, & Devloo-Delva, 2022).

Within the IWP, the mitogenome results demonstrated maternal structure at a finer spatial scale than the SNP results, supporting the
<figure><img src="image_7.png" /><figcaption>FIGURE 4 Reassignment success of *Carcharinus leucas* respective to the number of informative markers to assign simulated mixtures to their population of origin (E-PAC = eastern Pacific; W-ATL = western Atlantic; E-ATL = eastern Atlantic; IWP = Indo-West Pacific; Japan = Urauchi River, Japan; and Fiji). The simulated mixtures are based on a leave-one-out resampling method in *rubias*. The assignment success was tested based on 5–500 informative markers that were selected according to their high DAPC loadings contribution from the *adegenet* package.</figcaption></figure>
growing evidence of long-term female philopatry and male-biased dispersal in Bull Sharks (Karl et al., 2011; Sandoval Laurrabaquio-Alvarado et al., 2019; Tillett et al., 2012) and other elasmobranchs (reviewed in Chapman et al., 2015; Phillips et al., 2021). Specifically, the mtDNA shows a separation between the W-IO and E-IO/W-PAC, not evident from the nuclear DNA, and greater differentiation of the E-ATL. The presence of haplotypes from multiple distinct mitochondrial lineages in South Africa, Sri Lanka, Thailand, Fiji, and the Fitzroy River (Australia) represents signals of secondary female contact after historical isolation or low ongoing female gene flow. However, we see that almost every haplotype is unique to a sampling location (except for locations within Australia). This could indicate that each location has experienced sufficient female isolation for the mitogenome to mutate (i.e., complete lineage sorting), and/or is a consequence of high haplotype diversity as a result of sequencing complete mitogenomes. In addition, demographic events affect the non-recombining haploid mtDNA and recombining diploid nuclear DNA differently (Heist, 2012; Lawson Handley & Perrin, 2007), and, therefore, the mito-nuclear discordance may not necessarily be driven by sex-biased dispersal (Toews & Brelsford, 2012).

The spatial distribution of close-kin, and the cross-cohort HSPs in particular, is valuable to identify contemporary dispersal and sex bias therein (e.g., Feutry et al., 2017, 2020). In the IWP, we found that all FSPs and most HSPs were sampled from the same sampling location. The 32 same-river, cross-cohort HSPs had an equal likelihood of being maternally (same haplotype) or paternally (different haplotype) related, indicating contemporary biparental philopatry to a finer scale than detected by conventional population genetic methods (Tillett et al., 2012). However, the seven cross-river, cross-cohort HSPs in northern Australia were most likely paternally related, suggesting a bias toward male dispersal between rivers. This unconfounded result supports the mito-nuclear findings and suggests that male-biased dispersal may be occurring at a broader scale, yet this needs to be confirmed with larger samples sizes.

Within Japan, the SNP data grouped the individuals into three separate clusters, with individuals from Okinawa clustering more closely with the IWP and those from the Urauchi River forming two separate groups. The sharks from the Churaumi Aquarium, Okinawa were originally sourced from Japanese waters near Okinawa Island and several individuals were known to be related. The close similarity of the Okinawa group to the IWP cluster may represent ongoing genetic connectivity or recent colonization, where the strong nuclear differentiation may be caused by an accumulation of mutations at the edge of a range expansion due to the more pronounced effect of drift on a small and recently established population (Peischl et al., 2013). This divergence may be reinforced through high relatedness within a small population and seasonal differences in partition, which would also explain why we observed multiple genetic clusters within the Urauchi River. The sample size within Japan was too small to accurately estimate kinship; thus, family structure due to extremely small population size cannot be excluded (e.g., Devloo-Delva et al., 2019; Feutry et al., 2017). Similarly, the Fiji population

contains many related individuals (see Glaus et al., 2020), yet these were not removed from our analyses as they could signify an artifact of a small population size (Waples & Anderson, 2017). Mitochondrial haplotypes from individuals in Japan and Fiji are part of the E-IO/W-PAC haplogroup, supporting that these are recently established populations.

## 4.3 | Species, sex, provenance, and kinship identification

Increasingly, there is need to monitor wild populations and the global trade of wildlife products through DNA forensics (Cardeñosa et al., 2022). Our DArTcap panel of 3400 SNPs can identify species, sex, provenance, and kinship, which allows future monitoring of multiple demographic aspects at low cost (~AU$15 per sample; see Feutry et al., 2020). We were able to identify most of the individuals to species-level using the DArTseq and DArTcap data (i.e., Bull Shark, Pigeye Shark, Spinner Shark, Graceful Shark, Gray Reef Shark, Smalltail Shark, Dusky Shark, and Speartooth Shark), validated alongside mitochondrial species verification and a complete mtDNA reference database. Given the global scale of shark product trade (Clarke et al., 2006) and the morphological similarity between juveniles of many carcharhinid species, a selection of these SNP markers could be developed into a rapid tool for species identification of various shark products (Johri et al., 2019; Liu et al., 2017). While the species composition of shark fin trade has been studied using mtDNA (Cardeñosa et al., 2022; Fields et al., 2018), some species, like the Galapagos Shark and Dusky Shark, show mitochondrial introgression (Corrigan et al., 2017; Naylor et al., 2012). Here, the use of nuclear SNPs can resolve the ability to separate species with such introgression or incomplete lineage sorting (e.g., Kyne & Feutry, 2017; Liu et al., 2017). Moreover, the presence of sex-linked markers, specifically Y-linked markers, in the DArTseq data allowed us to identify the genetic sex of our samples when the visual sex information was missing. Unfortunately, the RNA baits failed to capture these markers appropriately with DArTcap sequencing, although future studies could assess sex ratios from fisheries or trade samples by redesigning the RNA baits for these Y-linked markers (e.g., Stovall et al., 2018).

The DArTcap panel also shows great promise and power to assign individuals to their respective populations (ocean-basin scale), and we estimated that a minimum of 100 markers are needed to achieve an accuracy of >80%. Few studies have specifically investigated the power of genetic markers to assign samples from traded shark products to their population of origin, and, to date, most forensic studies of shark trade have employed mtDNA to identify broad geographic origins (e.g., Cardeñosa et al., 2021; Chapman et al., 2009; Fields et al., 2020). Here, by genetically characterizing the source populations with 769 sharks and 3409 SNPs, the regional contributions of Bull Shark products to the global trade market can be estimated in future studies and highlight key regions of harvest to inform management actions. Furthermore, our results suggest
a diagnostic SNP panel of 100–500 markers could be designed to identify species, sex, and provenance from tissue samples. New genotyping technologies, such as DNA microarrays, may provide a more cost-efficient method for monitoring of Bull Shark exploitation (e.g., Arenas et al., 2017; Wenne et al., 2016). However, this panel needs further testing prior to use as an enforcement tool. This will involve obtaining additional samples from undersampled locations (e.g., E-PAC and E-ATL) and testing the assignment accuracy on these newly acquired samples.

Finally, we were able to identify close-kin relationships with the selected DArTcap markers. These relationships are important to estimate contemporary dispersal patterns (Feutry et al., 2017, 2020) or total adult abundance of a population (Bravington et al., 2016). This application requires a large sample size (relative to the total population size) and sufficient biological information (such as sex, age-at-length, and age-at-maturity), which can pose a challenge in studies of rare or threatened species. While the number of kin was limited in this dataset, the identified kin pairs provided insights into local demographic and genetic connectivity in Australia.

## 4.4 | Conservation and management considerations

Coastal shark species, such as the Bull Shark, are particularly susceptible to anthropogenic pressures such as fishing, shark control programs, and habitat modification, due to their close proximity to human populations (Knip et al., 2010). These species fulfill an important role in maintaining the ecosystem dynamics and functions (Ferretti et al., 2010), where a reduction in the abundance of top predators may lead to cascading effects in the food web and ecosystem health (e.g., Dudley & Simpfendorfer, 2006). Understanding the patterns of reproductive isolation of coastal sharks is essential to assess the impact of these threats at the population level (Frankham et al., 2017). Given our results of large-scale gene flow, we advocate for international collaboration to monitor the impacts on dispersing individuals (predominantly males) as they maintain the gene flow and genetic diversity between populations, and likely experience threats across many different countries. Similarly, coastal development near nursery areas may affect philopatric females disproportionately, with potential consequences for a wider region if other countries depend on recruitment from those nurseries. In addition, genetically isolated and small island populations, such as Japan and Fiji, currently face a number of threats, such as targets or bycatch in fisheries (Glaus et al., 2015) and require close monitoring as even low levels of catch may cause a population reduction. Any decline in abundance is unlikely to be replenished by migration from neighboring populations. Overall, this study suggests that coastal species with comparable characteristics exhibit similar dispersal patterns, and consequently, risk and vulnerability assessments and management actions should be considered at the smallest spatial scales, in accordance with the philopatric behaviors of these species, in order to prevent local depletion.

## AUTHOR CONTRIBUTIONS

**Floriaan Devloo-Delva:** Conceptualization (equal); data curation (lead); formal analysis (lead); funding acquisition (equal); investigation (lead); methodology (lead); project administration (supporting); resources (supporting); visualization (lead); writing – original draft (lead); writing – review and editing (equal). **Christopher Burridge:** Conceptualization (equal); funding acquisition (equal); investigation (supporting); project administration (equal); supervision (equal); writing – review and editing (equal). **Peter Kyne:** Conceptualization (equal); funding acquisition (equal); project administration (equal); resources (equal); supervision (equal); writing – review and editing (equal). **Juerg Brunnschweiler:** Resources (equal); writing – review and editing (equal). **Demian Chapman:** Resources (equal); writing – review and editing (equal). **Patricia Charvet:** Resources (equal); writing – review and editing (equal). **Xiao Chen:** Resources (equal); writing – review and editing (equal). **Geremy Cliff:** Resources (equal); writing – review and editing (equal). **Ryan Daly:** Resources (equal); writing – review and editing (equal). **J. Marcus Drymon:** Resources (equal); writing – review and editing (equal). **Mario Espinoza:** Resources (equal); writing – review and editing (equal). **Daniel Fernando:** Resources (equal); writing – review and editing (equal). **Laura Garcia Barcia:** Resources (equal); writing – review and editing (equal). **Kerstin Glaus:** Resources (equal); writing – review and editing (equal). **Blanca I. Gonzalez-Garza:** Resources (equal); writing – review and editing (equal). **Michael Grant:** Resources (equal); writing – review and editing (equal). **Rasanthi M. Gunasekera:** Investigation (equal); project administration (supporting); writing – review and editing (equal). **Sebastian Hernandez:** Resources (equal); writing – review and editing (equal). **Susumu Hyodo:** Resources (equal); writing – review and editing (equal). **Rima W. Jabado:** Resources (equal); writing – review and editing (equal). **Sébastien Jaquemet:** Resources (equal); writing – review and editing (equal). **Grant Johnson:** Resources (equal); writing – review and editing (equal). **James T. Ketchum:** Resources (equal); writing – review and editing (equal). **Hélène Magalon:** Funding acquisition (equal); resources (equal); writing – review and editing (equal). **James R. Marthick:** Investigation (equal); writing – review and editing (equal). **Frederik Mollen:** Resources (equal); writing – review and editing (equal). **Stefano Mona:** Conceptualization (equal); funding acquisition (equal); supervision (supporting); writing – review and editing (equal). **Gavin J.P. Naylor:** Conceptualization (supporting); resources (equal); writing – review and editing (equal). **John E.G. Nevill:** Resources (equal); writing – review and editing (equal). **Nicole M. Phillips:** Resources (equal); writing – review and editing (equal). **Richard D. Pillans:** Funding acquisition (equal); resources (equal); writing – review and editing (equal). **Bautisse D. Postaire:** Conceptualization (equal); data curation (supporting); funding acquisition (supporting); writing – review and editing (equal). **Amy F. Smoothey:** Resources (equal); writing – review and editing (equal). **Katsunori Tachihara:** Resources (equal); writing – review and editing (equal). **Bree J. Tillett:** Resources (equal); writing – review and editing (equal). **Jorge A. Valerio-Vargas:** Resources (equal); writing – review and editing (equal). **Pierre Feutry:** Conceptualization (equal); funding
acquisition (equal); investigation (supporting); project administration (lead); supervision (lead); writing - review and editing (equal).

## ACKNOWLEDGMENTS

This work was supported by the Sea World Research and Rescue Foundation Inc, the Ord River Research Offset grant through CSIRO, and the Marine Biodiversity Hub, a collaborative partnership supported through funding from the Australian Government's National Environmental Science Program. FDD was supported by a joint UTAS/CSIRO scholarship and the Quantitative Marine Science program. Samples from Réunion Island were collected as part of projects EcoRecoRun and Eurraica funded by the DEAL/SEB. We would like to acknowledge Lei Yang, William White, Jennifer Ovenden, Jessica Morgan, Thomas Poirout, Arie de Jong, Jeffrey de Pauw, A. Bodycott, J. Young, D. Young, A. Murray, O. Masens, M. Miller, D. Johnson, Betty Laglbauer, Sierra Fisheries, Pelagios Kakunja, Blue Resources Trust, NSW DPI Fisheries Bull Shark research, Beqa Adventure Divers, the Operations staff of the KZN Sharks Board (SAF), and the many Traditional Owners, assistants, and volunteers for their help with obtaining samples. In respect to this, we also thank all fishers and landing sites workers for their help in providing access to the sampled specimens. We appreciate the work from Marc Dando for providing the Bull Shark illustration. FDD thanks Jessica Farley and Peter Grewe for the aid in importing the tissue samples.

## CONFLICT OF INTEREST STATEMENT

We note no conflict or competing interests among the authors in relation to the information provided in this manuscript.

## OPEN RESEARCH BADGES

<figure><img src="image_8.png" /></figure>

This article has earned an Open Data badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at [https://doi.org/10.5061/dryad.4qrfj6qdg](https://doi.org/10.5061/dryad.4qrfj6qdg).

## DATA AVAILABILITY STATEMENT

Genetic data: Mitochondrial genomes are uploaded to GenBank (accession numbers: OP006762-OP007122). Mitochondrial DNA alignment and SNP genotypes are available on DataDryad (https://doi.org/10.5061/dryad.4qrfj6qdg).

Sample metadata: Sample metadata and Rmarkdown document, which included all performed analyses discussed in this study, are available on DataDryad (https://doi.org/10.5061/dryad.4qrfj6qdg).

Benefit-Sharing Statement: Benefits Generated: A research collaboration was developed with scientists from the countries providing genetic samples, all collaborators are included as co-authors, the results of research have been shared with the provider communities and the broader scientific community (see above), and the research addresses a priority concern, in this case the conservation of the Bull Shark. In Australia and Papua New Guinea, we consulted with

the Indigenous communities providing the biodiversity resources and worked with members of the communities to collect samples based on local knowledge. More broadly, our group is committed to international scientific partnerships, as well as institutional capacity building.

## ORCID

Floríaan Devloo-Delva https://orcid.org/0000-0002-6757-2949

Christopher P. Burridge https://orcid.org/0000-0002-8185-6091

Peter M. Kyne https://orcid.org/0000-0003-4494-2625

Juerg M. Brunnschweiler https://orcid.org/0000-0002-9901-3279

Demian D. Chapman https://orcid.org/0000-0002-7976-0947

Patricia Charvet https://orcid.org/0000-0002-8801-433X

Geremy Cliff https://orcid.org/0000-0003-1790-6055

Ryan Daly https://orcid.org/0000-0002-4409-6951

J. Marcus Drymon https://orcid.org/0000-0002-2104-004X

Mario Espinoza https://orcid.org/0000-0002-8355-2411

Daniel Fernando https://orcid.org/0000-0003-2163-1832

Laura Garcia Barcia https://orcid.org/0000-0002-0402-7564

Kerstin Glaus https://orcid.org/0000-0001-9985-2243

Michael I. Grant https://orcid.org/0000-0002-6127-8968

Rasanthi M. Gunasekera https://orcid.org/0000-0003-1990-2752

Sebastian Hernandez https://orcid.org/0000-0002-2908-6050

Susumu Hyodo https://orcid.org/0000-0001-8146-2386

Rima W. Jabado https://orcid.org/0000-0001-6239-6723

Sébastien Jaquemet https://orcid.org/0000-0003-4199-4657

James T. Ketchum https://orcid.org/0000-0002-7340-4749

Hélène Magalon https://orcid.org/0000-0002-7061-955X

Frederik H. Mollen https://orcid.org/0000-0002-9934-1029

Stefano Mona https://orcid.org/0000-0001-9087-0656

Gavin J. P. Naylor https://orcid.org/0000-0002-8731-2626

Nicole M. Phillips https://orcid.org/0000-0002-4138-4966

Richard D. Pillans https://orcid.org/0000-0001-6234-8764

Bautisse D. Postaire https://orcid.org/0000-0001-8566-1040

Amy F. Smoothey https://orcid.org/0000-0001-8271-686X

Bree J. Tillet https://orcid.org/0000-0002-0159-6649

Jorge A. Valerio-Vargas https://orcid.org/0000-0002-7036-8586

Pierre Feuty https://orcid.org/0000-0001-8492-5382

## REFERENCES

Ali, O. A., O'Rourke, S. M., Amish, S. J., Meek, M. H., Luikart, G., Jeffres, C., & Miller, M. R. (2016). RAD capture (rapture): Flexible and efficient sequence-based genotyping. *Genetics*, *202*(2), 389-400.
https://doi.org/10.1534/genetics.115.183665

Anderson, E. C., Waples, R. S., & Kalinowski, S. T. (2008). An improved method for predicting the accuracy of genetic stock identification. *Canadian Journal of Fisheries and Aquatic Sciences*, *65*(7), 1475-1486.
https://doi.org/10.1139/F08-049

Archer, F. I., Adams, P. E., & Schneiders, B. B. (2017). strataG: An R package for manipulating, summarizing and analysing population genetic data. *Molecular Ecology Resources*, *17*(1), 5-11.
https://doi.org/10.1111/1755-0998.12559

Arenas, M., Pereira, F., Oliveira, M., Pinto, N., Lopes, A. M., Gomes, V., Carracedo, A., & Amorim, A. (2017). Forensic genetics and
<figure><img src="image_9.png" /></figure>

genomics: Much more than just a human affair. *PLoS Genetics*, 13(9),
e1006960. https://doi.org/10.1371/journal.pgen.1006960

Avise, J. C., Bowen, B. W., Lamb, T., Meylan, A. B., & Bermingham, E. (1992). Mitochondrial DNA evolution at a turtle's pace: Evidence for low genetic variability and reduced microevolutionary rate in the Testudines. *Molecular Biology and Evolution*, *9*(3), 457–473.
https://doi.org/10.1093/oxfordjournals.molbev.a040735

Bailleul, D., Mackenzie, A., Sacchi, O., Poisson, F., Bierne, N., & Arnaud-Haond, S. (2018). Large-scale genetic panmixia in the blue shark (*Prionace glauca*): A single worldwide population, or a genetic lag-time effect of the “grey zone” of differentiation? *Evolutionary Applications*, *11*(5), 614–630. https://doi.org/10.1111/eva.12591

Begg, G. A., Friedland, K. D., & Pearce, J. B. (1999). Stock identification and its role in stock assessment and fisheries management: An overview. *Fisheries Research*, *43*(1–3), 1–8. https://doi.org/10.1016/S0165-7836(99)00062-4

Blaison, A., Jaquemet, S., Guyomard, D., Vangrevelynghie, G., Gazzo, T., Cliff, G., Cotel, P., & Soria, M. (2015). Seasonal variability of bull and tiger shark presence on the west coast of Reunion Island, western Indian Ocean. *African Journal of Marine Science*, *37*(2), 199–208.
https://doi.org/10.2989/1814232X.2015.1050453

Bowen, B. W., Gaither, M. R., DiBattista, J. D., Iacchei, M., Andrews, K. R., Grant, W. S., Toonen, R. J., & Briggs, J. C. (2016). Comparative phylogeography of the ocean planet. *Proceedings of the National Academy of Sciences*, *113*(29), 7962–7969. https://doi.org/10.1073/pnas.1602404113

Bravington, M. V., Skaug, H. J., & Anderson, E. C. (2016). Close-kin mark-recapture. *Statistical Science*, *31*(2), 259–274. https://doi.org/10.1214/16-Sts552

Brunnschweiler, J. M., Queiroz, N., & Sims, D. W. (2010). Oceans apart? Short-term movements and behaviour of adult bull sharks *Carcharhinus leucas* in Atlantic and Pacific oceans determined from pop-off satellite archival tagging. *Journal of Fish Biology*, 77(6), 1343–1458. https://doi.org/10.1111/j.1095-8649.2010.02757.x

Cardeñosa, D., Fields, A. T., Shea, S. K. H., Feldheim, K. A., & Chapman, D. D. (2021). Relative contribution to the shark fin trade of Indo-Pacific and eastern Pacific pelagic thresher sharks. *Animal Conservation*, *24*(3), 367–372. https://doi.org/10.1111/acv.12644

Cardeñosa, D., Shea, S. K., Zhang, H., Fischer, G. A., Simpfendorfer, C. A., & Chapman, D. D. (2022). Two thirds of species in a global shark fin trade hub are threatened with extinction: Conservation potential of international trade regulations for coastal sharks. *Conservation Letters*, e12910. https://doi.org/10.1111/conl.12910

Carlson, J. K., Ribera, M. M., Conrath, C. L., Heupel, M. R., & Burgess, G. H. (2010). Habitat use and movement patterns of bull sharks *Carcharhinus leucas* determined using pop-up satellite archi- val tags. *Journal of Fish Biology*, 77(3), 661–675. https://doi.org/10.1111/j.1095-8649.2010.02707.x

Carrillo, L., Johns, E. M., Smith, R. H., Lamkin, J. T., & Lagier, J. L. (2015). Pathways and hydrography in the mesoamerican barrier reef sys- tem part 1: Circulation. *Continental Shelf Research*, *109*, 164–176.
https://doi.org/10.1016/j.csr.2015.09.014

Carrillo, L., Lamkin, J. T., Johns, E. M., Vásquez-Yeomans, L., Sosa-Cordero, F., Malca, E., Smith, R. H., & Gerard, T. (2017). Linking oceano- graphic processes and marine resources in the western Caribbean Sea large marine ecosystem subarea. *Environment and Development*, *22*, 84–96. https://doi.org/10.1016/j.envdev.2017.01.004

Castro, J. I. (2011). *The sharks of North America* (Vol. 1). Oxford University Press.

Chapman, D. D., Feldheim, K. A., Papastamatiou, Y. P., & Hueter, R. E. (2015). There and back again: A review of residency and return mi- grations in sharks, with implications for population structure and management. *Annual Review of Marine Science*, *7*, 547–570. https:// doi.org/10.1146/annurev-marine-010814-015730

Chapman, D. D., Pinhal, D., & Shivji, M. S. (2009). Tracking the fin trade: Genetic stock identification in western Atlantic scalped

hammerhead sharks *Sphyrna lewini*. *Endangered Species Research*, *9*(3), 221–228. https://doi.org/10.3354/esr00241

Chen, K. Y., Marschall, E. A., Sovic, M. G., Fries, A. C., Gibbs, H. L., & Ludsin, S. A. (2018). assignPOP: An R package for population assignment using genetic, non-genetic, or integrated data in a machine-learning framework. *Methods in Ecology and Evolution*, *9*(2), 439–446. https://doi.org/10.1111/2041-210X.12897

Chen, X., Liu, M., Peng, Z. Q., & Shi, X. F. (2015). Mitochondrial ge- nome of the bull shark *Carcharhinus leucas* (*Carcharhiniformes*: Carcharhinidae). *Mitochondrial DNA*, *26*(6), 813–814. https://doi. org/10.3109/19401736.2013.855906

Chen, X., Wang, J.-J., Ai, W.-M., Chen, H., & Lin, H.-D. (2018). Phylogeography and genetic population structure of the spad- enose shark (*Scoliodon macrorrhynchos*) from the Chinese coast. *Mitochondrial DNA Part A*, *29*(7), 1100–1107. https://doi. org/10.1080/24701394.2017.1413363

Clarke, S. C., McAllister, M. K., Milner-Gulland, E. J., Kirkwood, G. P., Michielsens, C. G. J., Agnew, D. J., Pikitch, E. K., Nakano, H., & Shivji, M. S. (2006). Global estimates of shark catches using trade records from commercial markets. *Ecology Letters*, *9*(10), 1115–1126. https://doi.org/10.1111/j.1523-1739.2005.00247.x

Corrigan, S., Delser, P. M., Eddy, C., Duffy, C., Yang, L., Li, C., Bazinet, A. L., Mona, S., & Naylor, G. J. P. (2017). Historical introgression drives pervasive mitochondrial admixture between two species of pelagic sharks. *Molecular Phylogenetics and Evolution*, *110*, 122–126.
https://doi.org/10.1016/j.ympev.2017.03.011

Cowman, P. F., & Bellwood, D. R. (2013). Vicariance across major ma- rine biogeographic barriers: Temporal concordance and the rela- tive intensity of hard versus soft barriers. *Proceedings of the Royal Society B: Biological Sciences*, *280*(1768), 20131541. https://doi. org/10.1098/rspb.2013.1541

Daly, R., Smale, M. J., Cowley, P. D., & Froneman, P. W. (2014). Residency patterns and migration dynamics of adult bull sharks (*Carcharhinus leucas*) on the east coast of southern Africa. *PLoS One*, *9*(10), e109357. https://doi.org/10.1371/journal.pone.0109357

Daly-Engel, T. S., Seraphin, K. D., Holland, K. N., Coffey, J. P., Nance, H. A., Toonen, R. J., & Bowen, B. W. (2012). Global phyloge- raphy with mixed-marker analysis reveals male-mediated dis- persal in the endangered scalloped hammerhead shark (*Sphyrna lewini*). *PLoS One*, *7*(1), e29986. https://doi.org/10.1371/journal.pone.0029986

Deng, Z., Chen, J., Song, N., Li, Y., & Han, Z. (2019). Genetic homoge- neity among bull sharks *Carcharhinus leucas* in the South China Sea. *Pakistan Journal of Zoology*, *51*(4), 1281–1288. https://doi.org/10.17582/journal.pjz.2019.51.4.12811288

Devloo-Delva, F., Gosselin, T., Butcher, P. A., Grewe, P. M., Huverneers, C., Thomson, R. B., Werry, J. M., & Feutry, P. (2022). An R-based tool for identifying sex-linked markers from restriction site-associated DNA sequencing with applications to elasmobranch conservation. *Research Square*, 1–8.

Devloo-Delva, F., Maes, G. E., Hernández, S. I., McAllister, J. D., Gunasekera, R. M., Grewe, P. M., Thomson, R. B., & Feutry, P. (2019). Accounting for kin sampling reveals genetic connectivity in Tasmanian and New Zealand school sharks, *Galeorhinus galeus*. *Ecology and Evolution*, *9*, 4465–4472. https://doi.org/10.1002/ece3.5012

Dudgeon, C. L., Blower, D. C., Broderick, D., Giles, J. L., Holmes, B. J., Kashiwagi, T., Kruck, N. C., Morgan, J. A., Tillett, B. J., & Ovenden, J. R. (2012). A review of the application of molecular genetics for fisheries management and conservation of sharks and rays. *Journal of Fish Biology*, *80*(5), 1789–1843. https://doi.org/10.1111/j.1095-8649.2012.03265.x

Dudley, S. F., & Simpfendorfer, C. A. (2006). Population status of 14 shark species caught in the protective gillnets off KwaZulu-Natal beaches, South Africa, 1978–2003. *Marine and Freshwater Research*, *57*(2), 225–240. https://doi.org/10.1071/MF05156
Dutton, P. H., Bowen, B. W., Owens, D. W., Barragan, A., & Davis, S. K. (1999). Global phylogeography of the leatherback turtle (*Dermochelys coriacea*). *Journal of Zoology, 248*(3), 397–409. https://doi.org/10.1111/j.1469-7998.1999.tb01038.x

Edgar, R. C. (2004). MUSC: Multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Research, 32*(5), 1792–1797. https://doi.org/10.1093/nar/gkh340

Espinoza, M., Heupel, M. R., Tobin, A. J., & Simpfendorfer, C. A. (2016). Evidence of partial migration in a large coastal predator: Opportunistic foraging and reproduction as key drivers? *PLoS One*, 11(2), e0147608. https://doi.org/10.1371/journal.pone.0147608

Espinoza, M., Lédée, E. J. I., Smoothey, A. F., Heupel, M. R., Peddemors, V. M., Tobin, A. J., & Simpfendorfer, C. A. (2021). Intra-specific variation in movement and habitat connectivity of a mobile predator revealed by acoustic telemetry and network analyses. *Marine Biology*, 168(6), 80. https://doi.org/10.1007/s00227-021-03886-z

Estoup, A., Jarne, P., & Cornuet, J.-M. (2002). Homoplasny and mutation model at microsatellite loci and their consequences for population genetics analysis. *Molecular Ecology*, 11(9), 1591–1604. https://doi.org/10.1046/j.1365-294X.2002.01576.x

Ferretti, F., Worm, B., Britten, G. L., Heithaus, M. R., & Lotze, H. K. (2010). Patterns and ecosystem consequences of shark declines in the ocean. *Ecology Letters*, 13(8), 1055–1071. https://doi.org/10.1111/j.1461-0248.2010.01489.x

Feutry, P., Berry, O., Kyne, P. M., Pillans, R. D., Hillary, R. M., Grewe, P. M., Marthick, J. R., Johnson, G., Gunasekera, R. M., Bax, N. J., & Bravington, M. (2017). Inferring contemporary and historical genetic connectivity from juveniles. *Molecular Ecology*, 26(2), 444–456. https://doi.org/10.1111/mec.13929

Feutry, P., Devloo-Delva, F., Tran Lu, Y. A., Mona, S., Gunasekera, R. M., Johnson, G., Pillans, R. D., Jaccoud, D., Kilian, A., Morgan, D. L., Saunders, T., Bax, N. J., & Kyne, P. M. (2020). One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark. *Molecular Ecology Resources*, 20(6), 1470–1485. https://doi.org/10.1111/1755-0998.13204

Feutry, P., Kyne, P. M., Pillans, R. D., Chen, X., Naylor, G. J., & Grewe, P. M. (2014). Mitogenomics of the Speartooth shark challenges ten years of control region sequencing. *BMC Evolutionary Biology*, 14, 232. https://doi.org/10.1186/s12862-014-0232-x

Fields, A. T., Fischer, G. A., Shea, S. K. H., Zhang, H., Abercrombie, D. L., Feldheim, K. A., Babcock, E. A., & Chapman, D. D. (2018). Species composition of the international shark fin trade assessed through a retail-market survey in Hong Kong. *Conservation Biology*, 32(2), 376–389. https://doi.org/10.1111/cobi.13043

Fields, A. T., Fischer, G. A., Shea, S. K. H., Zhang, H., Feldheim, K. A., & Chapman, D. D. (2020). DNA zip-coding: Identifying the source populations supplying the international trade of a critically endangered coastal shark. *Animal Conservation*, 23(6), 670–678. https://doi.org/10.1111/acv.12585

Fontaine, M. C., Baird, S. J., Piry, S., Ray, N., Tolley, K. A., Duke, S., Birkun, A., Ferreira, M., Jauniaux, T., Llavona, A., Öztürk, B., Öztürk, A. A., Ridoux, V., Rogan, E., Sequeira, M., Siebert, U., Vikingsson, G. A., Bouquegneau, J.-M., & Michaux, J. R. (2007). Rise of oceanographic barriers in continuous populations of a cetacean: The genetic structure of harbour porpoises in Old World waters. *BMC Biology*, 5(30), 1–16. https://doi.org/10.1186/1741-7007-5-30

Foster, S. D., Feutry, P., Grewe, P. M., Berry, O., Hui, F. K., & Davies, C. R. (2018). Reliably discriminating stock structure with genetic markers: Mixture models with robust and fast computation. *Molecular Ecology Resources*, 18(6), 1310–1325. https://doi.org/10.1111/1755-0998.12920

Frankham, R., Ballou, J. D., Ralls, K., Eldridge, M., Dudash, M. R., Fenster, C. B., Lacy, R. C., & Sunnucks, P. (2017). Genetic management of fragmented animal and plant populations (1st ed.). Oxford University Press. https://doi.org/10.1093/oso/9780198783398.001.0001

Glaus, K. B., Adrian-Kalchhauser, I., Burkhardt-Holm, P., White, W. T., & Brunnschweiler, J. M. (2015). Characteristics of the shark fisheries of Fiji. *Scientific Reports*, 5(117556), 1–11. https://doi.org/10.1038/srep17556

Glaus, K. B., Appleyard, S. A., Stockwell, B., Brunnschweiler, J. M., Shivji, M., Clua, E., Marie, A. D., & Rico, C. (2020). Insights into insular isolation of the bull shark, *Carcharhinus leucas* (Müller and Henle, 1839), in Fijian Waters. *Frontiers in Marine Science*, 7(586015), 1–10. https://doi.org/10.3389/fmars.2020.586015

Gosselin, T., Lamothe, M., Devloo-Delva, F., & Grewe, P. (2020). radiator: RADseq data exploration, manipulation and visualization using R. R package version 1.1.9. https://thierrygosselin.github.io/radiator/

Goudet, J. (2005). Hierfstat, a package for R to compute and test hierarchical F-statistics. *Molecular Ecology Notes*, 5(1), 184–186. https://doi.org/10.1111/j.1471-8286.2004.00828.x

Green, M. E., Appleyard, S. A., White, W. T., Tracey, S. R., Devloo-Delva, F., & Ovenden, J. R. (2019). Novel multi-marker comparisons address the genetic population structure of silvertip sharks (*Carcharhinus albimarginatus*). *Marine and Freshwater Research*, 70(7), 1007–1019. https://doi.org/10.1071/MF18296

Green, M. E., Appleyard, S. A., White, W. T., Tracey, S. R., Heupel, M. R., & Ovenden, J. R. (2022). Updated connectivity assessment for the scalloped hammerhead (*Sphyrna lewini*) in Pacific and Indian oceans using a multi-marker genetic approach. *Fisheries Research*, 251, 106305. https://doi.org/10.1016/j.fishres.2022.106305

Green, M. E., Simpfendorfer, C. A., & Devloo-Delva, F. (2022). Population structure and connectivity of chondrichthyans. In J. C. Carrier, C. A. Simpfendorfer, M. R. Heithaus, & K. E. Yopak (Eds.), *Biology of sharks and their relatives* (3rd ed., pp. 523–544). CRC Press. https://doi.org/10.1201/9781003262190-16

Heist, E. J. (2012). Genetics of sharks, skates and rays. In J. C. Carrier, J. A. Musick, & M. R. Heithaus (Eds.), *Biology of sharks and their relatives* (2nd ed., pp. 487–504). CRC press. https://doi.org/10.1201/b11867

Heupel, M. R., Simpfendorfer, C. A., Espinoza, M., Smoothey, A. F., Tobin, A., & Peddemors, V. (2015). Conservation challenges of sharks with continental scale migrations. *Frontiers in Marine Science*, 2, Article 12. https://doi.org/10.3389/fmars.2015.00012

Hillary, R. M., Bravington, M. V., Patterson, T. A., Grewe, P., Bradford, R., Feutry, P., Gunasekera, R., Peddemors, V., Werry, J., Francis, M. P., Duffy, C. A. J., & Bruce, B. D. (2018). Genetic relatedness reveals total population size of white sharks in eastern Australia and New Zealand. *Scientific Reports*, 8, 2261. https://doi.org/10.1038/s41598-018-20593-w

Hirschfeld, M., Dudgeon, C., Sheaves, M., & Barnett, A. (2021). Barriers in a sea of elasmobranchs: From fishing for populations to testing hypotheses in population genetics. *Global Ecology and Biogeography*, 30, 2147–2163. https://doi.org/10.1111/geb.13379

Holmes, B. H., Steinke, D., & Ward, R. D. (2009). Identification of shark and ray fins using DNA barcoding. *Fisheries Research*, 95(2–3), 280–288. https://doi.org/10.1016/j.fishres.2008.09.036

Jaccoud, D., Peng, K., Feinstein, D., & Kilian, A. (2001). Diversity arrays: A solid state technology for sequence information independent genotyping. *Nucleic Acids Research*, 29(4), e25. https://doi.org/10.1093/nar/29.4.e25

Jensen, N. H. (1976). Reproduction of the bull shark, *Carcharhinus leucas*, in the Lake Nicaragua-Rio San Juan system. In T. B. Thorson (Ed.), *Investigations of the ichthyofauna of Nicaraguan Lakes* (pp. 539–559). University of Nebraska-Lincoln.

Johri, S., Solanki, J., Cantu, V. A., Fellows, S. R., Edwards, R. A., Moreno, I., Vyas, A., & Dinsdale, E. A. (2019). 'Genome skimming' with the MinION hand-held sequencer identifies CITES-listed shark species in India's exports market. *Scientific Reports*, 9(4476). https://doi.org/10.1038/s41598-019-40940-9

Jombart, T., & Ahmed, I. (2011). Adegenet 1.3-1: New tools for the analysis of genome-wide SNP data. *Bioinformatics*, 27(21), 3070–3071. https://doi.org/10.1093/bioinformatics/btr521
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. https://doi.org/10.1186/1471-2156-11-94

Karl, S. A., Castro, A. L. F., Lopez, J. A., Charvet, P., & Burgess, G. H. (2011). Phylogeography and conservation of the bull shark (*Carcharhinus leucas*) inferred from mitochondrial and microsatellite DNA. *Conservation Genetics*, *12*(2), 371-382. https://doi.org/10.1007/s10592-010-0145-1

Keenan, K., McGinnity, P., Cross, T. F., Crozier, W. W., & Prodöhl, P. A. (2013). *diveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors. Methods in Ecology and Evolution*, *4*(8), 782-788. https://doi.org/10.1111/2041-210X.12067

Kitamura, T., Takemura, A., Watabe, S., Taniuchi, T., & Shimizu, M. (1996). Mitochondrial DNA analysis for the cytochrome b gene and D-loop region from the bull shark *Carcharhinus leucas*. *Fisheries Science*, *62*(1), 21-27. https://doi.org/10.2331/fishsci.62.21

Knip, D. M., Heupel, M. R., & Simpfendorfer, C. A. (2010). Sharks in nearshore environments: Models, importance, and consequences. *Marine Ecology Progress Series*, *402*, 1-11. https://doi.org/10.3354/meps08498

Knowlton, N., Weigt, L. A., Solorzano, L. A., Mills, D. K., & Bermingham, E. (1993). Divergence in proteins, mitochondrial DNA, and reproductive compatibility across the isthmus of Panama. *Science*, *260*(5114), 1629-1632. https://doi.org/10.1126/science.8503007

Komoroske, L. M., Miller, M. R., O'Rourke, S. M., Stewart, K. R., Jensen, M. P., & Dutton, P. H. (2019). A versatile rapture (RAD-capture) platform for genotyping marine turtles. *Molecular Ecology Resources*, *19*(2), 497-511. https://doi.org/10.1111/1755-0998.12980

Kyne, P. M., & Feugry, P. (2017). Recreational fishing impacts on threatened river sharks: A potential conservation issue. *Ecological Management and Restoration*, *18*(3), 209-213. https://doi.org/10.1111/emr.12266

Lawson Handley, L. J., & Perrin, N. (2007). Advances in our understanding of mammalian sex-biased dispersal. *Molecular Ecology*, *16*(8), 1559-1578. https://doi.org/10.1111/j.1365-294X.2006.03152.x

Layton, K. K. S., Dempson, B., Snelgrove, P. V. R., Duffy, S. J., Messmer, A. M., Paterson, I. G., Jeffery, N. W., Kess, T., Horne, J. B., & Salisbury, S. J. (2020). Solving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish. *Evolutionary Applications*, *13*(5), 1055-1068. https://doi.org/10.1111/eva.12922

Lea, J. S., Humphries, N. E., Clarke, C. R., & Sims, D. W. (2015). To Madagascar and back: Long-distance, return migration across open ocean by a pregnant female bull shark *Carcharhinus leucas*. *Journal of Fish Biology*, *87*(6), 1313-1321. https://doi.org/10.1111/jfb.12805

Liu, J., Jiang, J., Song, S., Tornabene, L., Chabarria, R., Naylor, G. J. P., & Li, C. (2017). Multilocus DNA barcoding—species identification with multilocus data. *Scientific Reports*, *7*(16601). https://doi.org/10.1038/s41598-017-16920-2

Luikart, G., Kardos, M., Hand, B. K., Rajora, O. P., Aitken, S. N., & Hohenlohe, P. A. (2019). Population genomics: Advancing understanding of nature. In O. P. Rajora (Ed.), *Population genomics* (pp. 3-79). Springer. https://doi.org/10.1007/13836_2018_60

MacNeil, M. A., Chapman, D. D., Heupel, M., Simpfendorfer, C. A., Heithaus, M., Meekan, M., Harvey, E., Goetze, J., Kiszka, J., Bond, M. E., Currey-Randall, L. M., Speed, C. W., Sherman, C. S., Rees, M. J., Udyawer, V., Flowers, K. I., Clementi, G., Valentin-Albanese, J., Gorham, T., ... Cinner, J. E. (2020). Global status and conservation potential of reef sharks. *Nature*, *583*(7818), 801-806. https://doi.org/10.1038/s41586-020-2519-y

Martin, A. P., Naylor, G. J. P., & Palumbi, S. R. (1992). Rates of mitochondrial-DNA evolution in sharks are slow compared with mammals. *Nature*, *357*(6374), 153-155. https://doi.org/10.1038/357153a0

Matich, P., Heithaus, M. R., & Layman, C. A. (2011). Contrasting patterns of individual specialization and trophic coupling in two marine apex predators. *Journal of Animal Ecology*, *80*(1), 294-305. https://doi.org/10.1111/j.1365-2656.2010.01753.x

Mirams, A. G. K., Treml, E. A., Shields, J. L., Liggins, L., & Riginos, C. (2011). Vicariance and dispersal across an intermittent barrier: Population genetic structure of marine animals across the Torres Strait land bridge. *Coral Reefs*, *30*(4), 937-949. https://doi.org/10.1007/s00338-011-0767-x

Moran, B. M., & Anderson, E. C. (2019). Bayesian inference from the conditional genetic stock identification model. *Canadian Journal of Fisheries and Aquatic Sciences*, *76*(4), 551-560. https://doi.org/10.1139/cjfas-2018-0016

Moritz, C. (1994). Defining 'evolutionarily significant units' for conservation. *Trends in Ecology & Evolution*, *9*(10), 373-375. https://doi.org/10.1016/0169-5347(94)90057-4

Naylor, G. J. P., Caira, J. N., Jensen, K., Rosana, K. A. M., White, W. T., & Last, P. R. (2012). A DNA sequence-based approach to the identification of shark and ray species and its implications for global elasmobranch diversity and parasitology. *Bulletin of the American Museum of Natural History*, 2012(367), 1-262. https://doi.org/10.1206/754.1

Niella, Y., Peddemors, V. M., Green, M., Smoothey, A. F., & Harcourt, R. (2021). A "wicked problem" reconciling human-shark conflict, shark bite mitigation, and threatened species. *Frontiers in Conservation Science*, 2, 720741. https://doi.org/10.3389/fcosc.2021.720741

Nielsen, E. E., Cariani, A., Aoidh, E. M., Maes, G. E., Milano, I., Ogden, R., Taylor, M., Hemmer-Hansen, J., Babbucci, M., Bargelloni, L., Bekkevold, D., Diopere, E., Grenfell, L., Helyar, S., Limborg, M. T., Martinsohn, J. T., McEwing, R., Panitz, F., Patarnello, T., ... Carvalho, G. R. (2012). Gene-associated markers provide tools for tackling illegal fishing and false eco-certification. *Nature Communications*, *3*, 851. https://doi.org/10.1038/ncomms1845

O'Dea, A., Lessios, H. A., Coates, A. G., Eytan, R. I., Restrepo-Moreno, S. A., Cione, A. L., Collins, L. S., De Queiroz, A., Farris, D. W., & Norris, R. D. (2016). Formation of the isthmus of Panama. *Science Advances*, *2*(8), e1600883. https://doi.org/10.1126/sciadv.1600883

Oleksiak, M. F., & Rajora, O. P. (2020). *Population genomics: Marine organisms*. Springer. https://doi.org/10.1007/978-3-030-37936-0

Ovenden, J. R., Dudgeon, C., Feugry, P., Feldheim, K., & Maes, G. E. (2018). Genetics and genomics for fundamental and applied research on elasmobranchs. In J. C. Carrier, M. R. Heithaus, & C. A. Simpfendorfer (Eds.), *Shark research: Emerging technologies and applications for the Field and laboratory* (Vol. 1, pp. 235-254). CRC Press. https://doi.org/10.1201/b21842

Paradis, E. (2010). Pegas: An R package for population genetics with an integrated-modular approach. *Bioinformatics*, *26*(3), 419-420. https://doi.org/10.1093/bioinformatics/btp696

Pazmiño, D. A., Maes, G. E., Green, M. E., Simpfendorfer, C. A., Hoyos-Padilla, E. M., Duffy, C. J. A., Meyer, C. G., Kerwath, S. E., Salinas-de-León, P., & van Herwerden, L. (2018). Strong trans-Pacific break and local conservation units in the Galapagos shark (*Carcharhinus galapagensis*) revealed by genome-wide cytonuclear markers. *Heredity*, *120*(5), 407-421. https://doi.org/10.1038/s41437-017-0025-2

Peischl, S., Dupanloup, I., Kirkpatrick, M., & Excoffier, L. (2013). On the accumulation of deleterious mutations during range expansions. *Molecular Ecology*, *22*(24), 5972-5982. https://doi.org/10.1111/mec.12524

Phillips, N. M., Devloo-Delva, F., McCall, C., & Daly-Engel, T. (2021). Reviewing the genetic evidence for sex-biased dispersal in elasmobranchs. *Reviews in Fish Biology and Fisheries*, *31*(4), 821-841. https://doi.org/10.1007/s11160-021-09673-9

Pillans, R. D., & Franklin, C. E. (2004). Plasma osmolyte concentrations and rectal gland mass of bull sharks *Carcharhinus leucas*, captured along a salinity gradient. *Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology*, *138*(3), 363-371. https://doi.org/10.1016/j.cbpb.2004.05.006
Pirog, A., Ravigné, V., Fontaine, M. C., Rieux, A., Gilabert, A., Cliff, G., Clua, E., Daly, R., Heithaus, M. R., Kiszka, J. J., Matich, P., Nevill, J. E. G., Smoothey, A. F., Temple, A. J., Berggren, P., Jaquemet, S., & Magalon, H. (2019). Population structure, connectivity, and demographic history of an apex marine predator, the bull shark *Carcharhinus leucas*. *Ecology and Evolution*, *9*(23), 12980–13000.
https://doi.org/10.1002/ece3.5597

Portnoy, D. S., McDowell, J. R., Heist, E. J., Musick, J. A., & Graves, J. E. (2010). World phylogeography and male-mediated gene flow in the sandbar shark, *Carcharhinus plumbeus*. *Molecular Ecology*, *19*(10), 1994–2010.
https://doi.org/10.1111/j.1365-294X.2010.04626.x

Puckridge, M., Last, P. R., White, W. T., & Andreakis, N. (2013). Phylogeography of the indo-West Pacific maskrays (Dasyatidae, Neotrygon): A complex example of chondrichthyan radiation in the Cenozoic. *Ecology and Evolution*, *3*(2), 217–232.
https://doi.org/10.1002/ece3.448

Puechmaille, S. J. (2016). The program STRUCTURE does not reliably recover the correct population structure when sampling is uneven: Subsampling and new estimators alleviate the problem. *Molecular Ecology Resources*, *16*(3), 608–627.
https://doi.org/10.1111/1755-0998.12512

R Core Team. (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing Available online at
https://www.r-project.org/

Rigby, C. L., Espinoza, M., Derrick, D., Pacoureau, N., & Dicken, M. (2021). *Carcharhinus Leucas*. Downloaded on 15 May 2021. The IUCN Red List of Threatened Species 2021: e.T39372A2910670. https://doi.org/10.2305/IUCN.UK.2009-2.RLTS.T39372A10187195.en

Rocha, L. A., Craig, M. T., & Bowen, B. W. (2007). Phylogeography and the conservation of coral reef fishes. *Coral Reefs*, *26*(3), 501–512.
https://doi.org/10.1007/s00338-007-0261-7

Ronce, O. (2007). How does it feel to be like a rolling stone? Ten questions about dispersal evolution. *Annual Review of Ecology, Evolution, and Systematics*, *38*, 231–253.
https://doi.org/10.1146/annur.evolcosys.38.091206.095611

Sandoval Laurrabaquio-Alvarado, N., Díaz-Jaimes, P., Hinojosa-Álvarez, S., del Pilar Blanco-Parra, M., Adams, D. H., Pérez-Jiménez, J. C., & Castillo-Géniz, J. L. (2021). Mitochondrial DNA genome evidence for the existence of a third divergent lineage in the western Atlantic Ocean for the bull shark (*Carcharhinus leucas*). *Journal of Fish Biology*, *99*(1), 275–282.
https://doi.org/10.1111/jfb.14698

Sandoval Laurrabaquio-Alvarado, N., Islas-Villanueva, V., Adams, D. H., Uribe-Alcocer, M., Alvarado-Bremer, J. R., & Díaz-Jaimes, P. (2019). Genetic evidence for regional philopatry of the bull shark (*Carcharhinus leucas*), to nursery areas in estuaries of the Gulf of Mexico and western North Atlantic Ocean. *Fisheries Research*, *209*, 67–74.
https://doi.org/10.1016/j.fishres.2018.09.013

Schmidt, T. L., Jasper, M.-E., Weeks, A. R., & Hoffmann, A. A. (2021). Unbiased population heterozygosity estimates from genome-wide sequence data. *Methods in Ecology and Evolution*, *12*(10), 1888–1898.
https://doi.org/10.1111/2041-210X.13659

Simpfendorfer, C. A., & Milward, N. E. (1993). Utilisation of a tropical bay as a nursery area by sharks of the families Carcharinidae and Sphyrnidae. *Environmental Biology of Fishes*, *37*(4), 337–345.
https://doi.org/10.1007/BF00005200

Smoothey, A. F., Lee, K. A., & Peddemors, V. M. (2019). Long-term patterns of abundance, residency and movements of bull sharks (*Carcharhinus leucas*) in Sydney harbour, Australia. *Scientific Reports*, *9*(18864). https://doi.org/10.1038/s41598-019-54365-x

Stovall, W. R., Taylor, H. R., Black, M., Grosser, S., Rutherford, K., & Gemmell, N. J. (2018). Genetic sex assignment in wild populations using genotyping-by-sequencing data: A statistical threshold approach. *Molecular Ecology Resources*, *18*(2), 179–190.
https://doi.org/10.1111/1755-0998.12767

Testerman, C. B. (2014). Molecular ecology of globally distributed sharks [thesis, Nova southeastern university]. Fort Lauderdale, FL, USA.
https://nsuworks.nova.edu/occ_stuetd/6/

Tillett, B. J., Meekan, M. G., Field, I. C., Hua, Q., & Bradshaw, C. J. A. (2011). Similar life history traits in bull (*Carcharhinus leucas*) and pig-eye (*C. amboinensis*) sharks. *Marine and Freshwater Research*, *62*(7), 850–860.
https://doi.org/10.1071/Mf10271

Tillett, B. J., Meekan, M. G., Field, I. C., Thorburn, D. C., & Ovenden, J. R. (2012). Evidence for reproductive philopatry in the bull shark *Carcharhinus leucas*. *Journal of Fish Biology*, *80*(6), 2140–2158.
https://doi.org/10.1111/j.1095-8649.2012.03228.x

Toews, D. P., & Brelsford, A. (2012). The biogeography of mitochondrial and nuclear discordance in animals. *Molecular Ecology*, *21*(16), 3907–3930.
https://doi.org/10.1111/j.1365-294X.2012.05664.x

Trystram, C., Rogers, K. M., Soria, M., & Jaquemet, S. (2017). Feeding patterns of two sympatric shark predators in coastal ecosystems of an oceanic Island. *Canadian Journal of Fisheries and Aquatic Sciences*, *74*(2), 216–227.
https://doi.org/10.1139/cjfas-2016-0105

Vähä, J. P., Erkinaro, J., Niemelä, E., & Primmer, C. R. (2007). Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. *Molecular Ecology*, *16*(13), 2638–2654.
https://doi.org/10.1111/j.1365-294X.2007.03329.x

Voris, H. K. (2000). Maps of Pleistocene sea levels in Southeast Asia: Shorelines, river systems and time durations. *Journal of Biogeography*, *27*(5), 1153–1167.
https://doi.org/10.1046/j.1365-2699.2000.00489.x

Waples, R. S. (1998). Separating the wheat from the chaff: Patterns of genetic differentiation in high gene flow species. *Journal of Heredity*, *89*(5), 438–450.
https://doi.org/10.1093/jhered/89.5.438

Waples, R. S., & Anderson, E. C. (2017). Purging putative siblings from population genetic data sets: A cautionary view. *Molecular Ecology*, *26*(5), 1211–1224.
https://doi.org/10.1111/mec.14022

Waples, R. S., & 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. *Molecular Ecology*, *15*(6), 1419–1439.
https://doi.org/10.1111/j.1365-294X.2006.02890.x

Wenne, R., Drywa, A., Kent, M., Sundsaasen, K. K., & Lien, S. (2016). SNP arrays for species identification in salmonids. In J. M. Walker (Ed.), *Marine genomics* (pp. 97–111). Springer.
https://doi.org/10.1007/978-1-4939-3774-5_6

Yin, W., Fu, C., Guo, L., He, Q., Li, J., Jin, B., Wu, Q., & Li, B. (2009). Species delimitation and historical biogeography in the genus *Helice* (Brachyura: Varunidae) in the northwestern Pacific. *Zoological Science*, *26*(7), 467–475.
https://doi.org/10.2108/zzj.26.467

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Supporting Information section at the end of this article.

**How to cite this article:** Devloo-Delva, F., Burridge, C P., Kyne, P M., Brunnschweiler, J M., Chapman, D D., Charvet, P., Chen, X., Cliff, G., Daly, R., Drymon, J M., Espinoza, M., Fernando, D., Barcia, L G., Glaus, K., González-Garza, B I., Grant, M I., Gunasekera, R M., Hernandez, S., Hyodo, S. ... Feutry, P. (2023). From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator. *Ecology and Evolution*, *13*, e9837.
https://doi.org/10.1002/ece3.9837