Towards better management of Australia’s shark fishery: genetic analyses reveal unexpected ratios of cryptic blacktip species Carcharhinus tilstoni and C. limbatus

The common blacktip shark (Carcharhinus limbatus) and the Australian blacktip shark (C. tilstoni) are morphologically similar species that co-occur in subtropical and tropical Australia. In striking contrast to what has been previously reported, we demonstrate that the common blacktip shark is not rare in northern Australia but occurs in approximately equal frequencies with the Australian blacktip shark. Management of shark resources in northern Australia needs to take account of this new information. Species identification was performed using nucleotide sequences of the control, NADH dehydrogenase subunit 4 (ND4) and cytochrome oxidase I (COI) regions in the mitochondrial genome. The proportion of overall genetic variation (FST) between the two species was small (0.042, P < 0.01) based on allele frequencies at five microsatellite loci. We confirm that a third blacktip species (C. amblyrhynchoides, graceful shark) is closely related to C. tilstoni and C. limbatus and can be distinguished from them on the basis of mtDNA sequences from two gene regions. The Australian blacktip shark (C. tilstoni) was not encountered among 20 samples from central Indonesia that were later confirmed to be common blacktip and graceful sharks. Fisheries regulators urgently need new information on life history, population structure and morphological characters for species identification of blacktip shark species in Australia.

Evidence for extensive population structure in the white-spotted eagle ray within the Indo-Pacific inferred from mitochondrial gene sequences

The white-spotted eagle ray Aetobatus narinari is a species complex that occurs circumglobally throughout warm-temperate waters. Aetobatus narinari is semi-pelagic and large (up to 300 cm disc width), suggesting high dispersal capabilities and gene flow on a wide spatial scale. Sequence data from two mitochondrial genes, cytochrome b (cytb) and NADH dehydrogenase subunit 4 (ND4), were used to determine the genetic variability within and among 18 sampling locations in the central Indo-Pacific biogeographical region. Populations in the Indo-Pacific were highly genetically structured with c. 70% of the total genetic variation found among three geographical regions (East China Sea, Southeast Asia and Australia). FST was 0.64 for cytb and 0.53 for ND4, with φST values being even larger, that is, 0.78 for cytb and 0.65 for ND4. This high-level genetic partitioning provides strong evidence against extensive gene flow in A. narinari. The degree of genetic population structuring in the Indo-Pacific was similar to that found on a global scale. Global FST was 0.63 for cytb and 0.57 for ND4, and global φST values were 0.94 for cytb and 0.82 for ND4. This suggests that the A. narinari complex may be more speciose than the two or three species proposed to date. Further sampling and genetic analyses are likely to uncover the ‘evolutionarily significant’ and ‘management’ units that are critical to determine the susceptibilities of individual populations to regional fishing pressures and to provide advice on management options. Network analyses showed a close genetic relationship between haplotypes from the central Indo-Pacific and South Africa, providing support for a proposed dispersal pathway from the possible centre of origin of the A. narinari species complex in the Indo-Pacific into the Atlantic Ocean.

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum, is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: ‘Least Concern’ and ‘Vulnerable’. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric FST analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.

A mitochondrial species identification assay for Australian blacktip sharks (Carcharhinus tilstoni, C. limbatus and C. amblyrhynchoides) using real-time PCR and high-resolution melt analysis

Tropical Australian shark fisheries target two morphologically indistinguishable blacktip sharks, the Australian blacktip (Carcharhinus tilstoni) and the common blacktip (C. limbatus). Their relative contributions to northern and eastern Australian coastal fisheries are unclear because of species identification difficulties. The two species differ in their number of precaudal vertebrae, which is difficult and time consuming to obtain in the field. But, the two species can be distinguished genetically with diagnostic mutations in their mitochondrial DNA ND4 gene. A third closely related sister species, the graceful shark C. amblyrhynchoides, can also be distinguished by species-specific mutations in this gene. DNA sequencing is an effective diagnostic tool, but is relatively expensive and time consuming. In contrast, real-time high-resolution melt (HRM) PCR assays are rapid and relatively inexpensive. These assays amplify regions of DNA with species-specific genetic mutations that result in PCR products with unique melt profiles. A real-time HRM PCR species-diagnostic assay (RT-HRM-PCR) has been developed based on the mtDNA ND4 gene for rapid typing of C. tilstoni, C. limbatus and C. amblyrhynchoides. The assay was developed using ND4 sequences from 66 C. tilstoni, 33. C. limbatus and five C. amblyrhynchoides collected from Indonesia and Australian states and territories; Western Australia, the Northern Territory, Queensland and New South Wales. The assay was shown to be 100% accurate on 160 unknown blacktip shark tissue samples by full mtDNA ND4 sequencing.