- Research article
- Open Access
The extent of population genetic subdivision differs among four co-distributed shark species in the Indo-Australian archipelago
BMC Evolutionary Biology volume 9, Article number: 40 (2009)
The territorial fishing zones of Australia and Indonesia are contiguous to the north of Australia in the Timor and Arafura Seas and in the Indian Ocean to the north of Christmas Island. The area surrounding the shared boundary consists of a variety of bio-diverse marine habitats including shallow continental shelf waters, oceanic trenches and numerous offshore islands. Both countries exploit a variety of fisheries species, including whaler (Carcharhinus spp.) and hammerhead sharks (Sphyrna spp.). Despite their differences in social and financial arrangements, the two countries are motivated to develop complementary co-management practices to achieve resource sustainability. An essential starting point is knowledge of the degree of population subdivision, and hence fisheries stock status, in exploited species.
Populations of four commercially harvested shark species (Carcharhinus obscurus, Carcharhinus sorrah, Prionace glauca, Sphyrna lewini) were sampled from northern Australia and central Indonesia. Neutral genetic markers (mitochondrial DNA control region sequence and allelic variation at co-dominant microsatellite loci) revealed genetic subdivision between Australian and Indonesian populations of C. sorrah. Further research is needed to address the possibility of genetic subdivision among C. obscurus populations. There was no evidence of genetic subdivision for P. glauca and S. lewini populations, but the sampling represented a relatively small part of their distributional range. For these species, more detailed analyses of population genetic structure is recommended in the future.
Cooperative management between Australia and Indonesia is the best option at present for P. glauca and S. lewini, while C. sorrah and C. obscurus should be managed independently. On-going research on these and other exploited shark and ray species is strongly recommended. Biological and ecological similarity between species may not be a predictor of population genetic structure, so species-specific studies are recommended to provide new data to assist with sustainable fisheries management.
The Indo-Australian archipelago is a biogeograhically complex region encompassing a series of continental shelves, volcanic mountainous islands and deep-sea trenches. The Indonesian section straddles the equator and extends about 5000 km from east to west. It has the world's highest marine endemism; including most the diverse sea-grass meadows, greatest expanses of mangroves and most extensive coral reef communities . The Australian section is dominated by an extensive continental shelf extending offshore from three Australian state jurisdictions; Queensland, Northern Territory and Western Australia. As the world's fourth most populous nation, Indonesia has a high demand for food from the sea. The landed tonnage of Indonesian fisheries is orders of magnitude larger than Australian fisheries. Indonesian and Australian fishing zones are contiguous along much of the northern Australian Exclusive Economic Zone that extends 200 km from shore. There is considerable goodwill between the countries for cooperative management of shared fisheries stocks .
The Indo-Australian archipelago contains about 30% of the one thousand shark and ray species in the world . Many of the species are endemic and new species of sharks and rays continue to be described [eg. ]. Indonesia has the world's largest elasmobranch catch [118,000 t in 2003, 1], but data on catch and effort in the commercial and artisanal fisheries is severely lacking. In Australia, the elasmobranch catch was approximately 3000 tonnes in 2001 consisting of about 180 species as target and incidental (bycatch) catch . Only seven of these species are regarded as adequately managed and monitored, with the remainder being data-deficient for the purposes of assessing the sustainability of current catch rates . Faced with minimal species-specific information, fisheries managers have to generalise across biologically similar groups when making predictions about the effects of fisheries on shark and ray populations.
Population genetic analysis at and below the species level is an important starting point for species conservation and sustainable exploitation. The biological unit described by fisheries mathematical population models (sensu stock assessment) is self-maintained by reproduction and subsequent recruitment and experiences natural and human-mediated mortality. Individuals within the unit share the same gene-pool and are genetically isolated from other such groups . Population genetic analyses can identify the number of such management units within the geographical distribution of a target species . Population genetic analyses can also identify groups of populations that represent a significant part of the evolutionary legacy of the species, which may have a higher conservation status than other populations in the range of the species . Other methods, such as otolith microchemical analyses  and analyses of parasite distribution and abundance  can detect population groupings below the species level that may exist in response to, or in addition to, genetic differentiation between populations.
The aim of this study was to assess intra-specific population genetic subdivision for four shark species that are co-distributed in the Indo-Pacific and that support fisheries in northern Australia and Indonesia; Dusky Shark Carcharhinus obscurus, Spot-tail Shark Carcharhinus sorrah, Blue Shark Prionace glauca (Fa. Carcharhinidae) and Scalloped Hammerhead Sphyrna lewini (Fa. Sphyrnidae). Carcharhinids and Sphyrnids are placental livebearers with low intrinsic rates of increase. Prionace glauca and S. lewini, however, have notably higher fecundities than C. sorrah and C. obscurus . All are predators with naturally low abundances. Except C. sorrah, the study species have worldwide distributions; C. obscurus is found in coastal and offshore, but not oceanic, waters; P. glauca is oceanic but may be found close inshore when the continental shelf is narrow and S. lewini is semi-oceanic occurring over continental and insular shelves and adjacent deep water . Carcharhinus sorrah is restricted to the tropical Indo-west Pacific and is found on continental and insular shelves, primarily near reefs. The IUCN Red List of Threatened Species (Global) lists C. obscurus, P. glauca and S. lewini as 'near threatened' and C. sorrah as data deficient.
We expected that the extent of genetic subdivision between shark populations would be low because of their capacity to move extensively within the geographic constraints of the study area. Three of the four species (C. obscurus, P. glauca, S. lewini) attain large body sizes [greater than 3.0 m; C. sorrah are smaller, 1.5–3.0 m; 12]. All species are regarded as strong and active swimmers, but this may not lead to dispersal. The average distance moved was 50 km for 23 species of sharks tagged and recaptured off the northern Australia coastline between 1983 and 1997 , although the maximum distance moved was greater than 1100 km. Dispersal occurs during juvenile to adult phases, as there is no larval stage. Genetic subdivision has been demonstrated for S. lewini, but it was reported on a larger spatial scale than our study area [eg. between Pacific and Atlantic Ocean basins, ]. Genetic subdivision between populations was observed for another Carcharhinus species (C. limbatus) linked to the occurrence of female philopatry to nursery areas in the Atlantic coast of the US . There are few studies on elasmobranch population structure in the Indo-west Pacific. Dudgeon et al  reports population subdivision in the leopard shark (Stegastoma fasciatum) in Australia and south-east Asia, however this species is demersal and its capacity for dispersal is unknown. Barriers to dispersal in the Indo-Pacific could include deep-sea trenches and strong currents. Various phylogenetic 'breaks' for marine species from sea-horses to mackerel have been proposed [18–21]. The results are discussed as the type of the information needed by management agencies for sustainable harvesting plans of elasmobranchs in the bio-diverse Indo-Australian region.
Collection locations and DNA extraction
Tissue samples from C. obscurus, C. sorrah, P. glauca and S. lewini were collected from two areas in Australia and one area in central Indonesia (Fig. 1). An additional Australian location was included from the Gulf of Carpentaria for C. sorrah to provide another test of gene flow restriction between northern Australia and Indonesia. A more distant location (mid-north Pacific) was included for P. glauca to test the spatial extent of genetic homogeneity (Fig. 1). The mtDNA sequence of one sample of S. lewini from the Atlantic Ocean was included as a link between this study and that of Duncan et al .
In Australia, sharks were sampled from commercial catches by on-board observers and fisheries biologists. In Indonesia, samples were taken from landed catch at local markets (Muara Angke, western Java and Tanjung Luar, Lombok) by fisheries biologists. The fishing grounds represented were assumed to be within a 100 – 300 km radius of the markets. White muscle (~200 mg) was sampled from each specimen and preserved in either DMSO solution (20% dimethylsulphoxide in saturated NaCl solution) or 70% ethanol solution.
Total genomic DNA was extracted from ten to 50 mg of preserved tissue using the Chelex method . After heat denaturation and brief centrifugation (1200 g for 5 min), the supernatant containing genomic DNA was removed to a fresh tube for subsequent manipulation and storage.
Taxonomic field identifications were based on morphological features, verified with photographs where necessary. While identification of the four species in this study is relatively unproblematic for Carcharhinid and Sphyrnid sharks, samples lacking photo verification were confirmed by comparing mtDNA sequence to a reference data-base for Carcharhiniformes species maintained by JO's research group and by principal coordinates analysis of genetic distance using GenAlex 6.1  between all pairs of samples within a species based on their microsatellite genotype. Outliers and misidentified samples were removed from the study.
Mitochondrial DNA methods – laboratory
The 5' end of the control region (CR) of the mitochondrial DNA (mtDNA) was amplified and sequenced. Approximately 1145 bp (base pairs) were amplified using the forward primers ProL2 (5' CTGCCCTTGGCTCCCAAAGC 3') and the reverse primer PheCacaH2 (5' CTTAGCATCTTCAGTGCCAT 3') designed by Pardini et al. . Each 50 μl PCR reaction contained 1× PCR buffer, 200 μM of each dNTP, 1 μM of each primer, 2.5 mM of MgCl2, 4 units of Taq DNA polymerase (Qiagen, Doncaster, Vic, Australia) and 5 μl of DNA template. The PCR cycling conditions were 90 s at 94°C followed by 35 cycles of 5 s at 94°C, 30 s at 50°C and 30 s at 72°C with a final extension of 72°C for 5 min. Amplified DNA was purified using a commercial purification kit. The cycle sequencing reaction used ABI Big Dye Terminators v3.1®. Approximately 400 bp of the control region was sequenced in one direction using the forward primer. Sequences were obtained using an automated sequencer (Applied Biosystems 3130xl).
Microsatellite methods – laboratory
Shark samples were genotyped with three to five microsatellite loci (C. obscurus, four loci; C. sorrah, five loci; P. glauca, five loci and S. lewini, three loci). Microsatellite loci were sourced from Ovenden et al , Feldheim et al  and Keeney and Heist  and were generally applied to non-target species (ie. as cross-species amplifications). Microsatellite PCRs were performed in 96-well plates using a Perkin Elmer 9700 thermocycler. Reactions (10 μl) contained 1 μl of PCR buffer ® (Qiagen P/L, Doncaster, Vic, Australia) containing Tris-HCl (pH 8.7), KCl and (NH4)2 SO4; 4 mM MgCl2; 0.02 μM forward primer with an M13 extension [28, 29]; 0.2 μM reverse primer; 1.9 μM fluoro-labeled M13 primer; 0.3 units Taq DNA polymerase (Qiagen); 200 μM dNTP (Pharmacia Biotech, GE Healthcare Bio-Sciences P/L, Rydalmere, Australia); 1% bovine serum albumin and approximately 25 ng genomic DNA template. The DNA template and enzyme were denatured at 94°C for 1 min 30 s, followed by 35 cycles consisting of 94°C for 5 s, 60°C for 20 s and 72°C for 30 s and a final extension at 72°C for 30 min. Loci were amplified in separate reactions and then combined for fragment separation according to label colour and fragment size. Microsatellite fragment separation and scoring was performed using capillary electrophoresis (ABI3130xl). The size in base pairs of microsatellite amplicons was calculated to two decimal places and amplicons were allocated to a 'bin' that represented the mean allele size.
Mitochondrial DNA methods – data analyses
MtDNA control region sequence data was edited, aligned and checked by eye using the software Sequencher v 4.1 (Genecodes, Ann Arbor, MI, US). Haplotypes and numbers of variable sites for each species were determined using MacClade v3.08. Haplotype and nucleotide diversities were calculated using Arlequin version 3.11 . Intraspecific phylogenies among haplotypes were estimated using Bayesian [MrBayes v 3.1, 31], maximum likelihood (ML) and parsimony methods [PAUP* v 4.0b10, 32]. Candidate species for use as outgroups were evaluated on the availability of CR sequence, ease of sequence alignment, interspecies genetic similarity based on cytochrome oxidase I sequence  and robustness of topologies to alternate combinations of outgroups. For Bayesian and ML analyses, the most appropriate model of DNA evolution was selected using Akaike Information Criteria using ModelTest v 3.7 . The General Time Reversible (GTR) model with gamma distributed among-site rate variation was preferred for P. glauca (gamma 0.608). The GTR was also preferred for C. sorrah with equal rates of variation among sites. The Hasegawa Kishino Yano 85 model (HYK85) was preferred for C. obscurus and S. lewini. Gamma was 0.479 for S. lewini and zero for C. obscurus. ML phylogenies were estimated using PAUP* v 4.0b10  and nodal support was assessed using non-parametric bootstrapping involving 2000 pseudo-replicates of the original sequence alignment. Bayesian phylogenetic analysis was implemented using MrBayes v 3.1 . Four MCMC chains were run, starting from different random trees that were sampled every 100 generations. Each MCMC run consisted of one million generations after which the average standard deviation of split frequencies was less than 0.01 and the potential scale reduction factor was reasonably close to 1.0 for all parameters, as recommended by the MrBayes manual. The first 25% percent of generations were discarded as burn in. The posterior probability values, after burn in and across the four MCMC runs, were used to construct a 50% majority rule consensus phylogeny. Intraspecific hapotype networks were constructed using the statistical parsimony method of Templeton et al  with a 95% connection limit using the TCS v 1.21 software .
The proportion of total mtDNA CR sequence variation that was due to genetic differentiation between populations (pairwise ΦST) was measured for the one Indonesian and two Australian populations. Where available, additional populations were included in the pairwise analyses. Tamura-Nei distances  were used to describe sequence variation between haplotypes, with gamma values determined by ModelTest v 3.7 (above). The significance of ΦST values was tested by comparing the real ΦST with the ΦST values produced from 1023 random permutations of the data using Arlequin version 3.11 .
Microsatellite methods – data analyses
The null hypothesis of Hardy-Weinberg equilibrium was tested using Genepop v 4.0.7 . The program Micro-checker  was used to investigate likely causes for possible deviation from Hardy-Weinberg equilibrium. Microsatellite genetic diversity was characterised by the number of alleles per locus, expected (HE) and unbiased (UHE) heterozygosity, observed heterozygosity (HO) and fixation index using GenAlex 6.1 . The probability of rejecting the null hypothesis of genotypic disequilibrium between pairs of loci across populations was tested using Genepop v 4.0.7 . Bonferroni corrections to p-values for performing multiple tests was determined by the method of Sankoh et al . Microsatellite allelic diversity was used to investigate the degree of genetic subdivision between Australia and Indonesia for the four shark species. This was done using a standard population-pairwise FST approach  implemented in GenAlex v 6.1  and Genepop v 4.0.7 . Non-parametric bootstrapping was implemented in GenAlex v6.1 to estimate p-values over 999 randomizations of the data set. Missing data for population pair-wise comparisons was handled by interpolation in GenAlex v6.1.
Seven C. obscurus mtDNA haplotypes were found from 28 individuals sequenced. Haplotype frequencies varied among collection locations. For instance, haplotype CO7 was found in 15% of the West Australia samples and 28% of the samples from East Australia, but was not found in Indonesia (Table 1). Nucleotide and haplotype diversity also varied among locations. Both measures were highest in Indonesia and lowest in the two Australian locations (Table 1). Two haplotypes, which were only found in Indonesia (CO01 and CO02), were placed in a well-supported clade (Bayesian posterior probability, 100; ML bootstrap, 86), but other haplotypes only found in Indonesia (CO4, CO5) were not part of this clade (Fig. 2). There was evidence that haplotypes CO01 and CO02 formed a distinct lineage within a parsimony network of haplotypes (Fig. 2).
Carcharhinus obscurus microsatellite genotypes were obtained from four loci. Low numbers of microsatellite genotypes were obtained from East Australia, which had downstream analysis implications. For that population, small sample sizes increased the confidence limits around allele frequencies  and most likely caused disequilibrium in genotype proportions at two of the four loci (Table 2). One of these loci (CS02) was out of HW equilibrium among samples from other locations. The mean unbiased expected heterozygosity across populations and loci was 0.762. The mean number of alleles across populations was four for locus Cli108, 10.67 ± 3.33 (mean ± SE) for locus CS02, 13.33 ± 2.67 for locus CS06 and 3.33 ± 0.33 for locus LS24. Genotypic disequilibrium was not detected for pairwise comparisons of the three loci (ie. excluding CS02) in the three populations.
Among C. obscurus collected from Indonesia and West Australia, genetic population subdivision was found using mtDNA control region sequence (ΦST 0.191), but not microsatellite allele frequencies (Table 3). This did not change when the analysis was repeated without the locus that was not in Hardy-Weinberg proportions (CS02). Comparisons involving East Australia across the two classes of genetic markers were not significant.
Twelve mtDNA haplotypes were found for C. sorrah, and there were pronounced differences in their frequencies between the three collection locations. Two haplotypes (CS01 and CS10) were only found in Indonesia, nine haplotypes were only found in Australian samples and only haplotype (CS02) was shared between Indonesia and Australia (Table 4). Phylogenetic analysis separated the Australian and Indonesian sequences into two well-supported clades (Fig. 3). The haplotype that was shared between Australia and Indonesia (CS02) showed close similarity to the Indonesian group of haplotypes. A single shark with haplotype CS02 was collected from Australian waters (Gulf of Carpentaria) suggesting that it, or its maternal ancestor, may have dispersed from Indonesia. Haplotypes found in Indonesia (CS01, CS02 and CS10) formed a cluster on the parsimony network and their connection to the remainder was equally parsimonious via either CS03 or CS04.
Carcharhinus sorrah samples were assayed in sufficient numbers (43.3 ± 2.6, mean ± SE) for statistical analyses from four populations using five microsatellite loci (Table 5). There were three instances of genotypic proportions not following Hardy-Weinberg expectations, but as they were not associated with particular loci or collection locations, their effect on subsequent analyses were judged to be slight. Unbiassed heterozygosity was 0.55 ± 0.07 and the number of alleles was 3.50 ± 0.96 for locus LS15, 12.50 ± 1.66 for locus CT05, 6.25 ± 0.18 for locus CS12, 5.50 ± 0.96 for locus Cli100 and 26.25 ± 1.44 for locus CS08 across loci and populations. Genotypic disequilibrium was not detected for pairwise comparisons of the five loci in the four populations.
Population subdivision for C. sorrah was detected between Indonesia and the three Australian collection locations, but not within Australian waters, and the pattern of genetic population subdivision was similar for mtDNA control region sequences and microsatellite allele frequencies. F-statistics based on mtDNA ranged from 0.751 to 0.903 and for microsatellite data from 0.038 to 0.047 (Table 6).
Sixteen mtDNA control region haplotypes were identified (Table 7) with little or no evidence of partitioning among populations. The most common haplotypes (PG03 and PG06) were evenly distributed across the four collection locations. Only the rare haplotypes; for example PG10, PG15 and PG16, were found in one geographic location; this apparent range restriction was because they were only sampled once each (Table 7). Haplotype and nucleotide diversity was similar across collection locations. There are no well-supported clades on the phylogenetic tree of haplotypes that show distinctiveness, and the haplotype network has no distinct structure (Fig. 4).
Prionace glauca samples from four populations were assayed with five microsatellite loci. Sample sizes were adequate, except for the number of genotypes produced using locus Cli100 on samples from the Mid North Pacific (Table 8). Two populations (West Australia and East Australia) showed evidence of lack of equilibrium for genotype proportions for two loci, CS02 and CT06. Widespread disequilibrium in genotype proportions for locus CS02, which was reported here for C. obscurus samples, was not a feature of the results for P. glauca, except for samples from West Australia where null alleles were implicated. Detailed analyses of locus CT06 raw data suggested that allele stuttering might have resulted in scoring errors, as there was a significant shortage of heterozygote genotypes with alleles of one repeat unit difference. Instances of departures from Hardy-Weinberg proportions for P. glauca data were judged to be minor and unlikely to bias results. Unbiassed heterozygosity for P. glauca loci was 0.58 ± 0.05 and the number of alleles was 15.00 ± 1.08 for locus CS02, 5.00 for locus CT04, 5.75 ± 0.25 for locus CT06, 3.00 ± 0.41 for locus Cli107 and 4.25 ± 0.75 for locus Cli100 across loci and populations. Genotypic disequilibrium was not detected for pairwise comparisons of the five loci in the four populations.
There were no instances of genetic subdivision between collection locations for P. glauca from Australia, Indonesia or the northern Pacific Ocean with either class of genetic marker (Table 9).
Eight mtDNA control region haplotypes were identified amongst Indonesian and Australian samples (Table 10). Haplotype SL02 was the most common and was found in all collection locations at similar frequencies. As such it provided little to no evidence of population genetic structure on this geographic scale. Likewise, the distribution of remaining haplotypes provided no consistent information about gene flow between populations. An Atlantic Ocean haplotype [number 16, 15] was included here as SL09 to highlight the divergence between the Indo-Pacific and Atlantic Ocean areas. Haplotypes SL03 and SL07 formed a well-supported clade in the phylogeny and were outliers within the haplotype network (Fig. 5). Suprisingly, they were as distinct from the remainder as the Atlantic haplotype (SL09). They were sampled between one (SL07) and four times (SL03) from West Australia and Indonesia populations.
Sample sizes of genotypes assayed with microsatellite loci were adequate for S. lewini, with the mean (± SE) number being 52.78 (6.17) over the three collection locations and three microsatellite loci (Table 11). However, the degree of polymorphism of one locus (CT06) was poor. It was monomorphic in Indonesia and East Australia, had only two alleles in West Australia and heterozygotes were not detected. Lack of Hardy-Weinberg equilibrium was detected for locus Cli100 in Indonesia and locus CT07 in West Australia, which detailed analyses suggested was due to cryptic allelic stuttering or presence of null alleles. Genotypic disequilibrium was not detected for pairwise comparisons of the three loci in the three populations.
As for P. glauca, genetic subdivision among S. lewini populations from Indonesia, East Australia and West Australia was not detected with mtDNA control region sequences or microsatellite allele frequencies (Table 12).
We have shown here that fisheries stocks of two species (P. glauca and S. lewini) most likely extend across the economic zones of Australia and Indonesia, suggesting that both countries are likely to be exploiting the same resource. Conversely, there is evidence of genetic population subdivision between the fisheries stocks of the spot-tail shark (C. sorrah) of the two countries, and there is some genetic evidence that a fourth species (dusky shark, C. obscurus) may also be subdivided along these lines. Prionace glauca and S. lewini are relatively large, ocean-going shark species, which would be expected to move freely across the Arafura and Timor Seas that are bisected by the fishing zones of the two countries.
While neither C. sorrah or C. obscurus are off-shore or truly pelagic species the finding of population subdivision was not expected. However, the robust design of this study adds weight to this finding. Two classes of independent genetic markers were used. Mitochondrial DNA is cytoplasmic, whereas microsatellite loci are derived from DNA in the cell nucleus. MtDNA is maternally inherited, while microsatellite loci are biparentally inherited . Nucleotide sequencing was used for mtDNA, while variation in allelic frequencies was assessed for microsatellite loci.
Concordance between markers is taken as robust evidence of population structure [eg. [44, 45]] and is considered best-practise when defining fisheries stocks . The finding of genetic subdivsion between C. sorrah populations is emphasised by the apparent lack of subdivision among samples of other species (eg. P. glauca and S. lewini) that were collected from the same locations and analysed with the same suite of genetic methods. The genetic population subdivision between northern Australia and Indonesia reported here for C. sorrah needs to be tested with other methods, such as parasite distribution and abundance [eg. ] or stable isotope analysis of hard parts [eg. ].
The application of a range of non-genetic methods in further testing fisheries stock boundaries in all four shark species studied here would be valuable in view of the challenges faced in this study with the application of microsatellite markers. Cross-species, rather than species-specific, loci were used due to the limited availability of microsatellite loci for sharks and rays. This may have been responsible for comparatively inconsistent amplification of loci across samples within a species, leading to lower than ideal sample sizes. Other problems that may be potentially clarified in future by using species-specific loci may include Hardy-Weinberg disequilibrium and null alleles.
There are 31 species in the genus Carcharhinus . The majority are found close to shore and are susceptible to exploitation. Nine globally distributed species (C. albimarginatus, C. altimus, C. brevipinna, C. falciformis, C. leucas, C limbatus, C. longimanus, C. obscurus, and C. plumbeus) occur in the Indo-Pacific, and the distribution of a further 14 species is centred on this region; C. amblyrhynchoides (also found in Indian Ocean), C. amboinensis (with Indian), C. amblyrhynchos (with west Pacific)C. borneensis, C. cautus, C. dussumieri (with Indian), C. fitzroyensis, C. hemiodon (with Indian), C. macloti (with Indian), C. melanopterus (with Indian and west Pacific), C. sealei (with Indian), C. sorrah (with Indian and west Pacific), and C. tilstoni and Carcharhinus sp. A. Population genetic subdivision has not been investigated for any of these species in the Indo-Pacific or elsewhere, except for C. sorrah (this study) and C. limbatus . Keeney and Heist  found two major clades between ocean basins (western Atlantic, Gulf of Mexico, and Caribbean Sea clades, and eastern Atlantic, Indian, and Pacific Ocean clades) based on mtDNA control region sequence, and shallow population structure within ocean basins. We strongly recommend the assessment of population structure of other Carcharhinus species in the Indo-Pacific to assist with their conservation. The group of species has high regional biodiversity, economic importance and high levels of exploitation, which predisposes them to non-sustainable harvesting.
This is underlined by our finding that at least one Indo-Pacific species (C. sorrah) is potentially exposed to local population depletion. There is concern for the status of C. sorrah populations in Indonesia as fishing methods target this species, but the relative catch rate is consistently below that of the fishery in northern Australia . Here we report that Indonesian and northern Australian populations are genetically subdivided. Consequently, the exploited Indonesian populations may not be replenished by migration, which justifies concerns for their susceptibility to overexploitation. As C. sorrah is commonly reported from inshore waters less than 200 m, the deep water of the Timor Trench (2 – 3000 m) is implicated as the isolating mechanism preventing gene flow between northern Australia and central Indonesia. The Timor Trench divides the northwestern edge of the Sahul continental shelf from the islands of central Indonesia. Deep water was also implicated in the stock structure of commercial snappers in this region [44, 50]. The degree of population genetic subdivision in C. sorrah in northern Australia in this study was similar to the overall low levels of genetic subdivision reported in a previous allozyme study [data not shown, ]. We agree with Lavery and Shaklee  that the species should be managed as a single unit in Australia.
Our mtDNA sequence data suggested that genetic population subdivision might exist between C. obscurus Australia and Indonesia. We strongly recommend that this hypothesis is tested with further studies as our sampling strategy, particularly on the Australian west coast was sparse. However, Indonesian samples possessed two haplotypes (CO1 and CO2) that were not found in Australia, and there was evidence that those haplotypes had a common evolutionary history. The same conclusion would have been reached if samples from West Australia were grouped into two sampling locations rather than one. Unfortunately, the power of microsatellite data to resolve population structure in this species may have been compromised for the comparison between East Australia and Indonesia, as sample sizes were low, but sample sizes were adequate for the comparison between West Australia and Indonesia.
In contrast to the genetic differentiation between northern Australian and Indonesian populations for C. sorrah and possibly C. obscurus, no differentiation is reported here on this scale for two other shark species; P. glauca and S. lewini. Furthermore, the inclusion of a northern Pacific Ocean sample of P. glauca suggests that genetic homogeneity in this species extends over large geographic scales. Both P. glauca and S. lewini are distributed worldwide, have large body sizes and hence are likely to have high dispersal capacity. Duncan et al  confirmed this for S. lewini. They found little mtDNA control region sequence variaton among samples from in-shore nursery grounds within ocean basins. However, they did discover an upper limit to gene flow in females, as genetic subdivision was pronounced between Pacific, Indian and Atlantic Ocean basins. Duncan et al  also found that control region haplotype diversity was pronounced within the Indo-Pacific and our study confirms this. Indonesian and west Australian haplotypes SL03 and SL07 were as distinct from the majority of S. lewini haplotypes as SL09, which was sampled from the North Atlantic Ocean and which has the same sequence as haplotype 16 from Duncan et al . Prionace glauca populations have not been previously assessed for degree of genetic subdivision. This study and the work by Duncan et al  on a similarly oceanic species (S. lewini) provides an hypothesis that can be tested in future studies: genetic homogeneity in P. glauca extends within, but not between, ocean basins. Further research that contributes to the sustainable management of this species is urgently needed. Clarke et al  used DNA-based and morphological species identification methods to show that P. glauca dominated (17%, by weight) dried shark fin auctioned at a Hong Kong market.
This study on shark population subdivision parallels previous research on red and gold snapper in northern Australia and Indonesia, which led to cooperative fisheries management actions directed at ensuring sustainability. Four Lutjanidae snapper species with similar biology were compared from the same collection locations. The species had the same habitat preferences (marine, reefs to 180 m) except for L. argentimaculatus, whose juveniles inhabit freshwater before becoming fully marine on maturation. Unlike shark species, snapper species are highly fecund (106 eggs/female/spawning season) with pelagic eggs and larvae. But like sharks, there was no expectation of population subdivision in the Indo-Pacific. One species (Pristipomoides multidens) had pronounced genetic subdivision [50, 53]. The two red snapper species had moderate [Lutjanus erythropterus, ] and low [Lutjanus malabaricus, ] degrees of subdivision. The fourth species [Lutjanus argentimaculatus, ] was not genetically subdivided. The presence of genetic subdivision was explained by previously unrecognised site fidelity across life history stages. New knowledge about the geographic extent of L. malabaricus fisheries stocks in the region allowed the development of a biomass dynamic stock assessment model, which showed that current levels of harvesting were unsustainable. Australia and Indonesia are working together to co-manage shared stocks and fishing practises for snapper species have been modified . It is expected that the work reported here will contribute to co-management actions on shared shark stocks.
Genetic analyses of population structure are an important component of fisheries management. When populations are genetically subdivided this information can be used to define the geographic boundaries of a fisheries stocks. Boundaries provide the confidence to use estimates of recruitment, growth and mortality (natural and fishing) to describe fisheries populations in mathematical models. However, when genetic analyses of population structure suggest the absence of genetic subdivision, this implies fisheries stocks are widely distributed. If this occurs, two courses of action are appropriate. Firstly, further research should be sponsored to test the hypothesis of no subdivision using a range of genetic and non-genetic methods. Researchers need to keep management agencies up-to-date with latest results, in case stock subdivision is found and management can occur independently. Secondly, management arrangements need to be applied and coordinated among authorities responsible. It is not sufficient for exploitation in only one part of the stock to be regulated, as unregulated exploitation elsewhere in the stock could cause uniform depletion across the entire stock. Thus agencies across Australian states (eg. Queensland, Northern Territory and Western Australia) should cooperate to manage C. sorrah. Likewise, Australian and Indonesian management agencies should cooperate to manage Indo-Pacific P. glauca and S. lewini populations.
Using mtDNA and microsatellite loci, this study has contributed to the definition of stock boundaries for one (C. sorrah) and maybe two (C. obscurus) exploited shark species, and has shown that the nations of Australia and Indonesia are most likely exploiting the same stocks of P. glauca and S. lewini. Australian authorities are aware of the need for catch limits on C. sorrah and Indonesia would be wise to take similar action. Indonesia is drafting a national plan of action for shark exploitation. Generally, we have also demonstrated the difficulties associated with predicting the presence of genetic subdivision even for co-distributed and ecologically similar species. For sustainable management and conservation purposes, we advocate a species-by-species assessment of population structure, at least until clear patterns emerge in the bio-diverse Indo-Pacific region.
White WT, Last PR, JD S, Yearsley GK, Fahmi , Dharmadi : Economically important sharks and rays of Indonesia. 2006, Canberra, Australia, Australian Centre for International Agricultural Research
Blaber SJM, Dichmont CM, Buckworth RC, Badrudin , Sumiono B, Nurhakim S, Iskandar B, Fegan B, Ramm DC, Salini JP: Shared stocks of snappers (Lutjanidae) in Australia and Indonesia: Integrating biology, population dynamics and socio-economics to examine management scenarios. Reviews in Fish Biology and Fisheries. 2005, 15: 111-127. 10.1007/s11160-005-3887-y.
Last PR, Stevens JD: Sharks and Rays of Australia. CSIRO Australia. 1994
White WT, Last PR, Compagno LJV: Description of a new species of weasel shark, Hemigaleus australiensis n sp (Carcharhiniformes: Hemigaleidae) from Australian waters. Zootaxa. 2005, 37-49.
Stevens J, Simpfendorfer C, Francis M: Regional overview – southwest Pacific, Australasia and Oceania. Sharks, Rays and Chimeras: The Status of the Chondrichthyan Fishes. Edited by: Fowler S, Cavanagh R, Camhi M, Burgess G, Cailliet G, Fordham S, Simpfendorfer C, Musick J. 2005, Cambridge, UK: International Union for the Conservation of Nature and Natural Resources, 161-171.
Shark Advisory Group, Lack M: National Plan of Action for the Conservation and Management of Sharks (Shark-plan). 2004, Canberra, Australian Government Department of Agriculture, Fisheries and Forestry
King M: Fisheries Biology, Assessment and Management. 2007, Oxford, Blackwell, Second
Ovenden JR: Mitochondrial DNA and marine stock assessment: a review. Australian Journal of Marine and Freshwater Research. 1990, 41: 835-853. 10.1071/MF9900835.
Moritz C: Defining "Evolutionary Significant Units" for conservation. Trends in Ecology and Evolution. 1994
Cappo M, Eden P, Newman SJ, Robertson S: A new approach to validation of periodicity and timing of opaque zone formation in the otoliths of eleven species of Lutjanus from the central Great Barrier Reef. Fishery Bulletin. 2000, 98: 474-488.
Moore BR, Buckworth RC, Moss H, Lester RJG: Stock discrimination and movements of narrow-barred Spanish mackerel across northern Australia as indicated by parasites. Journal of Fish Biology. 2003, 63: 765-779. 10.1046/j.1095-8649.2003.00190.x.
Cortes E: Life history patterns, demography, and population dynamics. Biology of sharks and their relatives. Edited by: Carrier JC, Musick JA, Heithaus MR. 2004, CRC Press, 449-469.
Compagno LJV, Dando M, Fowler S: Sharks of the World. 2005, Princeton and Oxford, Princeton University Press
Stevens JD, West GJ, McLoughlin KJ: Movements, recapture patterns, and factors affecting the return rate of carcharhinid and other sharks tagged off northern Australia. Marine and Freshwater Research. 2000, 51: 127-141. 10.1071/MF98158.
Duncan KM, Martin AP, Bowen BW, De Couet HG: Global phylogeography of the scalloped hammerhead shark (Sphyrna lewini). Molecular Ecology. 2006, 15: 2239-2251. 10.1111/j.1365-294X.2006.02933.x.
Keeney DB, Heupel M, Hueter RE, Heist EJ: Genetic heterogeneity among blacktip shark, Carcharhinus limbatus, continental nurseries along the US Atlantic and Gulf of Mexico. Marine Biology. 2003, 143: 1039-1046. 10.1007/s00227-003-1166-9.
Dudgeon CL, Broderick D, Ovenden JR: IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegastoma fasciatum in the Indo-West Pacific. Mol Ecol. 18 (2): 248-261. 10.1111/j.1365-294X.2008.04025.x.
Benzie JAH: Genetic structure of marine organisms and SE Asian biogeography. Biogeography and Geological Evolution. Edited by: Hall R, Holloway JD. 1998, Leiden, The Netherlands: Backhuys Publishers, 197-209.
Barber PH, Palumbi SR, Erdman MV, Moosa MK: A marine Wallace's line?. Nature. 2000, 406: 692-693. 10.1038/35021135.
Lourie SA, Vincent ACJ: A marine fish follows Wallace's Line: the phylogeography of the three-spot seahorse (Hippocampus trimaculatus, Syngnathidae, Teleostei) in Southeast Asia. Journal of Biogeography. 2004, 31: 1975-1985. 10.1111/j.1365-2699.2004.01153.x.
Sulaiman ZH, Ovenden JR: Population genetic evidence for the east-west division of the Spanish mackerel (Scomberomorus commerson, Perciformes, Teleostei) along Wallace's Line. Biodiversity and Conservation.
Walsh P, Metzger D, Higuchi R: Chelex 100 as a medium for simple extraction of DNA for PCR-based typing forensic material. Biotechniques. 1991, 10 (4): 506-513.
Peakall R, Smouse PE: GenAlEx 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.
Pardini AT, Jones CS, Noble LR, Kreiser B, Malcolm H, Bruce BD, Stevens JD, Cliff G, Scholl MC, Francis M, et al: Sex-biased dispersal of great white sharks – In some respects, these sharks behave more like whales and dolphins than other fish. Nature. 2001, 412: 139-140. 10.1038/35084125.
Ovenden JR, Broderick D, Street R: Microsatellite primers for two carcharinid sharks (Carcharinus tilstoni and C. sorrah) and their usefulness across a wide range of shark species. Molecular Ecology Notes. 2006, 6: 415-418. 10.1111/j.1471-8286.2005.01254.x.
Feldheim KA, Gruber SH, Ashley MV: Population genetic structure of the lemon shark (Negaprion brevirostris) in the western Atlantic: DNA microsatellite variation. Molecular Ecology. 2001, 10: 295-303. 10.1046/j.1365-294x.2001.01182.x.
Keeney DB, Heist EJ: Characterization of microsatellite loci isolated from the blacktip shark and their utility in requiem and hammerhead sharks. Mol Ecol Notes. 2003, 3: 501-504. 10.1046/j.1471-8286.2003.00492.x.
Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nature Biotechnology. 2000, 18: 233-234. 10.1038/72708.
Broderick D, Idaghdour Y, Korrida A, Hellmich J: Gene flow in great bustard populations across the Strait of Gibraltar as elucidated from excremental PCR and mtDNA sequencing. Conservation Genetics. 2003, 4: 793-800. 10.1023/B:COGE.0000006111.65204.c9.
Schneider S, Kueffer JM, Roessli D, Excoffier L: Arlequin: A software for population genetic data analysis. Version 2.0. 2000, Geneva: Genetics and Biometry Lab, Department of Anthropology, University of Geneva
Ronquist F, Huelsebeck J: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony. 1999, Champaign, Illinios., Sinauer Associates, 4.0
Ward RD, Holmes BH, White WT, Last PR: DNA barcoding in Australasian chondrichthyans: results and potential uses in conservation. Marine and Freshwater Research. 2008, 59: 57-71. 10.1071/MF07148.
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132 (2): 619-633.
Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Molecular Ecology. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.
Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10 (3): 512-526.
Rousset F: Genepop'007: a complete re-implementation of the genepop software for Windows and Linux. Molecular Ecology Resources. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICROCHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.
Sankoh A, Huque M, Dubey S: Some comments on frequently used multiple endpoint adjustments methods in clinical trials. Statistics in Medicine. 1997, 16: 2529-2542. 10.1002/(SICI)1097-0258(19971130)16:22<2529::AID-SIM692>3.0.CO;2-J.
Weir BS, Cockerham CC: Estimating F-stastistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Bartley DM, Kent DB, Drawbridge MA: Conservation of genetic diversity in a white seabass hatchery enhancement program in southern California. Uses and effects of cultured fishes in aquatic ecosystems. Edited by: Schramm HL, Piper RG. 1995, Betheseda, MD, USA: American Fisheries Society, 15: 249-258.
Frankham R, Ballou JD, Briscoe DA: Introduction to Conservation Genetics. 2002, Cambridge, UK, Cambridge University Press
Salini JP, Ovenden JR, Street R, Pendrey R, Haryanti , Ngurah : Genetic population structure of Lutjanus malabaricus and L. erythropterus between eastern Indonesia and Australia: reconciling allozyme and mtDNA evidence. J Fish Biol. 2006, 68: 217-234. 10.1111/j.0022-1112.2006.001060.x.
Hoarau G, Piquet AMT, Veer van der HW, Rijnsdorp AD, Stam W, Olsen JL: Population structure of plaice (Pleuronectes platessa L.) in northern Europe: a comparison of resolving power between microsatellites and mitochondrial DNA data. Journal of Sea Research. 2004, 51: 183-190. 10.1016/j.seares.2003.12.002.
Shaklee JB, Currens KP: Genetic stock identification and risk assessment. Population genetics: principles and applications for fisheries scientists. Edited by: Hallerman EM. 2003, Bethesda, Maryland: American Fisheries Society, 291-328.
Newman S, Steckis R, Edmonds J, Lloyd J: Stock structure of the goldband snapper, Pristipomoides multidens (Pisces:Lutjanidae) from the waters of northern and western Australia by stable isotope ratio analysis of sagittal otolith carbonate. Mar Ecol Prog Ser. 2000, 198: 239-247. 10.3354/meps198239.
Keeney DB, Heist EJ: Worldwide phylogeography of the blacktip shark (Carcharhinus limbatus) inferred from mitochondrial DNA reveals isolation of western Atlantic populations coupled with recent Pacific dispersal. Molecular Ecology. 2006, 15: 3669-3679. 10.1111/j.1365-294X.2006.03036.x.
Cavanagh R, Kyne P, Fowler S, Musick J, Bennett M, eds: The Conservation Status of Australasian Chondrichthyans. Report of the IUCN Shark Specialist Group Australia and Oceania Region Red List Workshop. 2003, Brisbane, Australia: University of Queensland, School of Biomedical Sciences
Ovenden J, Salini JP, Street R, O'Connor S: Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle. Molecular Ecology. 2004, 13: 1991-1999. 10.1111/j.1365-294X.2004.02210.x.
Lavery S, Shaklee JB: Population genetics of two tropical sharks, Carcharhinus tilstoni and C. sorrah, in northern Australia. Australian Journal of Marine and Freshwater Research. 1989, 40: 541-557. 10.1071/MF9890541.
Clarke SC, Magnussen JE, Abercrombie DL, McAllister MK, Shivji MS: Identification of shark species composition and proportion in the Hong Kong shark fin market based on molecular genetics and trade records. Conservation Biology. 2006, 20: 201-211. 10.1111/j.1523-1739.2005.00247.x.
Ovenden JR, Lloyd J, Newman SJ, Keenan CP, Slater LS: Spatial genetic subdivision between northern Australian and southeast Asian populations of Pristipomoides multidens: a tropical marine reef fish species. Fisheries Research. 2002, 59: 57-69. 10.1016/S0165-7836(01)00415-5.
Ovenden JR, Street R: Genetic population structure of mangrove jack, Lutjanus argentimaculatus (Forsskäl, 1775). Marine and Freshwater Research. 2003, 54: 127-137. 10.1071/MF02142.
We sincerely thank tissue sample collectors including Fahmi and Dharmadi (Agency of Marine and Fisheries Research, Republic of Indonesia), William White (CSIRO Marine Laboratories, Hobart), Stirling Peverell (Queensland Department of Primary Industries and Fisheries), Rik Buckworth and Chris Tarca (Northern Territory Department of Business, Industry and Resource Development), Rory McAuley (Western Australia Department of Fisheries), Richard Pillans (CSIRO Marine Laboratories) and Mahmood Shivji (Guy Harvey Research Institute, US). Raewyn Street performed the laboratory analyses. We are grateful to Wayne Sumpton and Jess Morgan for comments on earlier versions of this manuscript. The Australian Fisheries Research and Development Corporation and the Australian Centre for Agricultural Research made financial contributions to this study.
JO and JS developed, designed and coordinated the study. JO wrote and revised the manuscript and performed statistical analyses. TK carried out mtDNA sequence alignments and statistical analyses. DB and JG participated in data analyses and interpretation. All authors edited the manuscript drafts and approved the final version.
About this article
Cite this article
Ovenden, J.R., Kashiwagi, T., Broderick, D. et al. The extent of population genetic subdivision differs among four co-distributed shark species in the Indo-Australian archipelago. BMC Evol Biol 9, 40 (2009). https://doi.org/10.1186/1471-2148-9-40
- Microsatellite Locus
- Continental Shelf
- Collection Location
- Shark Species
- Genetic Subdivision