An annotated checklist of the chondrichthyans of South Africa

An annotated checklist of chondrichthyan fishes (sharks, batoids, and chimaeras) occurring in South African waters is presented. The checklist is the result of decades of research and on-going systematic revisions of the regional fauna. The chondrichthyan fauna of South Africa is one of the richest in the world with 191 species, comprising 50 families and 103 genera. It consists of 30 families, 64 genera, and 111 species of sharks; 17 families, 36 genera, and 72 species of batoids; and, 3 families, 5 genera, and 8 species of chimaeras. The most species-rich shark families are the whaler sharks Carcharhinidae with 20 species followed by the deepwater catsharks Pentanchidae with 13 species. The most species-rich batoid families are the hardnose stakes Rajidae with at least 21 species followed by the stingrays Dasyatidae with 13 species. This monograph represents the first detailed annotated checklist of chondrichthyans from South Africa in over 30 years.

Determining the appropriate pretreatment procedures and the utility of liver tissue for bulk stable isotope (δ13 C and δ15 N) studies in sharks

Stable-isotope analysis (SIA) provides a valuable tool to address complex questions pertaining to elasmobranch ecology. Liver, a metabolically active, high turnover tissue (~166 days for 95% turnover), has the potential to reveal novel insights into recent feeding/movement behaviours of this diverse group. To date, limited work has used this tissue, but ecological application of SIA in liver requires consideration of tissue preparation techniques given the potential for high concentrations of urea and lipid that could bias δ¹³ C and δ¹⁵ N values (i.e., result in artificially lower δ¹³ C and δ¹⁵ N values). Here we investigated the effectiveness of (a) deionized water washing (WW) for urea removal from liver tissue and (b) chloroform-methanol for extraction of lipids from this lipid rich tissue. We then (a) established C:N thresholds for deriving ecologically relevant liver isotopic values given complications of removing all lipid and (b) undertook a preliminary comparison of δ¹³ C values between tissue pairs (muscle and liver) to test if observed isotopic differences correlated with known movement behaviour. Tests were conducted on four large shark species: the dusky (DUS, Carcharhinus obscurus), sand tiger (RAG, Carcharias taurus), scalloped hammerhead (SCA, Sphyrna lewini) and white shark (GRE, Carcharodon carcharias). There was no significant difference in δ¹⁵ N values between lipid-extracted (LE) liver and lipid-extracted/water washed (WW) treatments, however, WW resulted in significant increases in %N, δ¹³ C and %C. Following lipid extraction (repeated three times), some samples were still biased by lipids. Our species-specific "C:N thresholds" provide a method to derive ecologically viable isotope data given the complexities of this lipid rich tissue (C:N thresholds of 4.0, 3.6, 4.7 and 3.9 for DUS, RAG, SCA and GRE liver_LEWW tissue, respectively). The preliminary comparison of C:N threshold corrected liver and muscle δ¹³ C values corresponded with movement/habitat behaviours for each shark; minor differences in δ¹³ C values were observed for known regional movements of DUS and RAG (δ¹³ C_Diffs = 0.24 ± 0.99‰ and 0.57 ± 0.38‰, respectively), while SCA and GRE showed greater differences (1.24 ± 0.63‰ and 1.08 ± 0.71‰, respectively) correlated to large-scale movements between temperate/tropical and pelagic/coastal environments. These data provide an approach for the successful application of liver δ¹³ C and δ¹⁵ N values to examine elasmobranch ecology.

Spatio-temporal genetic tagging of a cosmopolitan planktivorous shark provides insight to gene flow, temporal variation and site-specific re-encounters

Migratory movements in response to seasonal resources often influence population structure and dynamics. Yet in mobile marine predators, population genetic consequences of such repetitious behaviour remain inaccessible without comprehensive sampling strategies. Temporal genetic sampling of seasonally recurring aggregations of planktivorous basking sharks, Cetorhinus maximus, in the Northeast Atlantic (NEA) affords an opportunity to resolve individual re-encounters at key sites with population connectivity and patterns of relatedness. Genetic tagging (19 microsatellites) revealed 18% of re-sampled individuals in the NEA demonstrated inter/multi-annual site-specific re-encounters. High genetic connectivity and migration between aggregation sites indicate the Irish Sea as an important movement corridor, with a contemporary effective population estimate (Ne) of 382 (CI = 241-830). We contrast the prevailing view of high gene flow across oceanic regions with evidence of population structure within the NEA, with early-season sharks off southwest Ireland possibly representing genetically distinct migrants. Finally, we found basking sharks surfacing together in the NEA are on average more related than expected by chance, suggesting a genetic consequence of, or a potential mechanism maintaining, site-specific re-encounters. Long-term temporal genetic monitoring is paramount in determining future viability of cosmopolitan marine species, identifying genetic units for conservation management, and for understanding aggregation structure and dynamics.

New polymorphic microsatellite loci revealed for the dusky shark Carcharhinus obscurus through Ion Proton double-digest RAD sequencing

The non-model shark species, dusky shark Carcharhinus obscurus, is a bio-economically and recreationally important shark in many areas of its range. Despite of the fishery importance of C. obscurus few genetic resources are currently available for the species. Here, we report on the isolation of eight novel microsatellite loci from C. obscurus using a double-digest restriction site associated DNA (RAD) sequencing approach on the Ion Proton semiconductor platform (ddRADseq-ion). We characterised the loci in 26 individuals and all loci were polymorphic, exhibiting 5-10 alleles (average 6.6), and observed and expected heterozygosities of 0.385-0.962 and 0.479-0.847, respectively. We found that all pairs of loci were in linkage equilibrium and conformed to Hardy-Weinberg expectations. The loci reported in this study are only the second set of microsatellite loci ever characterized for C. obscurus and will be valuable for molecular ecology studies for this vulnerable species.

The tooth, the whole tooth and nothing but the tooth: tooth shape and ontogenetic shift dynamics in the white shark Carcharodon carcharias

Results from this study of the white shark Carcharodon carcharias include measurements obtained using a novel photographic method that reveal significant differences between the sexes in the relationship between tooth cuspidity and shark total length, and a novel ontogenetic change in male tooth shape. Males exhibit broader upper first teeth and increased distal inclination of upper third teeth with increasing length, while females do not present a consistent morphological change. Substantial individual variation, with implications for pace of life syndrome, was present in males and tooth polymorphism was suggested in females. Sexual differences and individual variation may play major roles in ontogenetic changes in tooth morphology in C. carcharias, with potential implications for their foraging biology. Such individual and sexual differences should be included in studies of ontogenetic shift dynamics in other species and systems.

Diet and trophic ecology of the tiger shark (Galeocerdo cuvier) from South African waters

Knowledge of the diet and trophic ecology of apex predators is key for the implementation of effective ecosystem as well as species-based management initiatives. Using a combination of stomach content data and stable isotope analysis (δ15N and δ13C) the current study provides information on size-based and sex-specific variations in diet, trophic position (TP) and foraging habitat of tiger sharks (Galeocerdo cuvier) caught in the KwaZulu-Natal Sharks Board bather protection program. This study presents the longest time-series and most detailed analysis of stomach content data for G. cuvier worldwide. Prey identified from 628 non-empty stomachs revealed a size-based shift in diet. Reptiles, birds, mysticetes, and large shark species increased in dietary importance with G. cuvier size, concomitant with a decrease in smaller prey such as batoids and teleosts. Seasonal and decadal shifts in diet driven primarily by changes in the importance of elasmobranchs and mammal (cetacean) prey were recorded for medium sized (150-220 cm) G. cuvier. Both stomach content and stable isotope analysis indicated that G. cuvier is a generalist feeder at the population level. Size-based δ13C profiles indicated a movement to offshore foraging habitats by larger G. cuvier. Calculated TP varied by method ranging from 4.0 to 5.0 (TPSCA for stomach contents) and from 3.6 to 4.5 (TPscaled and TPadditive for δ15N). Large (> 220 cm) G. cuvier did not feed at discrete trophic levels, but rather throughout the food web. These data provide key information on the ecological role of G. cuvier to improve the accuracy of regional food web modelling. This will enable a better understanding of the ecological impacts related to changes in the abundance of this predator.

Assessing multiple paternity in three commercially exploited shark species: Mustelus mustelus, Carcharhinus obscurus and Sphyrna lewini

In this study, multiple paternity (MP) was investigated in three commercially important shark species, common smoothhound Mustelus mustelus, dusky shark Carcharhinus obscurus and scalloped hammerhead Sphyrna lewini occurring in southern Africa. Reduced marker panels of between five and six microsatellite loci were constructed for each species and used to genotype and assess the presence of MP in a total of 60 M. mustelus individuals from six litters, 90 C. obscurus individuals from 14 litters and 54 S. lewini individuals from 13 litters. Analysis in GERUD and COLONY revealed the presence of MP in all three species. Multiple paternities were observed in 67, 35 and 46% of the litters of M. mustelus, C. obscurus and S. lewini, with corresponding average sire size of 1·6, 1·4 and 2·0, respectively. The variation in the rate of MP among the three species is in accordance with previous studies whilst the comparatively high frequency of MP observed for M. mustelus, matches what has previously been reported for shark species demonstrating aggregation behaviour.

Global versus local causes and health implications of high mercury concentrations in sharks from the east coast of South Africa

Conservation concern regarding the overharvest of global shark populations for meat and fin consumption largely surrounds documented deleterious ecosystem effects, but may be further supported by improved knowledge of possibly high levels in their edible tissues (particularly meat) of the neurotoxin, methylmercury (CH3Hg). For many regions, however, little data exist on shark tissue Hg concentrations, and reasons for Hg variation within and among species or across regions are poorly understood. We quantified total Hg (THg) in 17 shark species (total n=283) from the east coast of South Africa, a top Hg emitter globally. Concentrations varied from means of around 0.1 mg kg(-1) dry weight (dw) THg in hardnose smoothhound (Mustelus mosis) and whale (Rhincodon typus) sharks to means of over 10 mg kg(-1) dw in shortfin mako (Isurus oxyrinchus), scalloped hammerhead (Sphyrna lewini), white (Carcharodon carcharias) and ragged-tooth (Carcharias taurus) sharks. These sharks had higher THg levels than conspecifics sampled from coastal waters of the North Atlantic and North, mid-, and South Pacific, and although sampling year and shark size may play a confounding role, this result suggests the potential importance of elevated local emissions. Values of THg showed strong, species-specific correlations with length, and nearly half the remaining variation was explained by trophic position (using nitrogen stable isotopes, δ(15)N), whereas measures of foraging habitat (using carbon stable isotopes, δ(13)C) were not significant. Mercury concentrations were above the regulatory guidelines for fish health effects and safe human consumption for 88% and 70% of species, respectively, suggesting on-going cause for concern for shark health, and human consumers of shark meat.

Comparative organochlorine accumulation in two ecologically similar shark species (Carcharodon carcharias and Carcharhinus obscurus) with divergent uptake based on different life history

Trophic position and body mass are traits commonly used to predict organochlorine burdens. Sharks, however, have a variety of feeding and life history strategies and metabolize lipid uniquely. Because of this diversity, and the lipid-association of organochlorines, the dynamics of organochlorine accumulation in sharks may be predicted ineffectively by stable isotope-derived trophic position and body mass, as is typical for other taxa. The present study compared ontogenetic organochlorine profiles in the dusky shark (Carcharhinus obscurus) and white shark (Carcharodon carcharias), which differ in metabolic thermoregulation and trophic position throughout their ontogeny. Although greater organochlorine concentrations were observed in the larger bodied and higher trophic position white shark (e.g., p,p'-dichlorodiphenyldichloroethylene: 20.2 ± 2.7 ng/g vs 9.3 ± 2.2 ng/g in the dusky shark), slopes of growth-dilution corrected concentrations with age were equal to those of the dusky shark. Similar ontogenetic trophic position increases in both species, less frequent white shark seal predation than previously assumed, or inaccurate species-specific growth parameters are possible explanations. Inshore habitat use (indicated by δ(13)C values) and mass were important predictors in white and dusky sharks, respectively, of both overall compound profiles and select organochlorine concentrations. The present study clarified understanding of trophic position and body mass as reliable predictors of interspecific organochlorine accumulation in sharks, whereas regional endothermy and diet shifting were shown to have less impact on overall rates of accumulation.

Laser photogrammetry improves size and demographic estimates for whale sharks

Whale sharks Rhincodon typus are globally threatened, but a lack of biological and demographic information hampers an accurate assessment of their vulnerability to further decline or capacity to recover. We used laser photogrammetry at two aggregation sites to obtain more accurate size estimates of free-swimming whale sharks compared to visual estimates, allowing improved estimates of biological parameters. Individual whale sharks ranged from 432-917 cm total length (TL) (mean ± SD = 673 ± 118.8 cm, N = 122) in southern Mozambique and from 420-990 cm TL (mean ± SD = 641 ± 133 cm, N = 46) in Tanzania. By combining measurements of stranded individuals with photogrammetry measurements of free-swimming sharks, we calculated length at 50% maturity for males in Mozambique at 916 cm TL. Repeat measurements of individual whale sharks measured over periods from 347-1,068 days yielded implausible growth rates, suggesting that the growth increment over this period was not large enough to be detected using laser photogrammetry, and that the method is best applied to estimating growth rates over longer (decadal) time periods. The sex ratio of both populations was biased towards males (74% in Mozambique, 89% in Tanzania), the majority of which were immature (98% in Mozambique, 94% in Tanzania). The population structure for these two aggregations was similar to most other documented whale shark aggregations around the world. Information on small (<400 cm) whale sharks, mature individuals, and females in this region is lacking, but necessary to inform conservation initiatives for this globally threatened species.

The last frontier: catch records of white sharks (Carcharodon carcharias) in the Northwest Pacific Ocean

White sharks are highly migratory apex predators, globally distributed in temperate, sub-tropical, and tropical waters. Knowledge of white shark biology and ecology has increased recently based on research at known aggregation sites in the Indian, Atlantic, and Northeast Pacific Oceans; however, few data are available for the Northwest Pacific Ocean. This study provides a meta-analysis of 240 observations of white sharks from the Northwest Pacific Ocean between 1951 and 2012. Records comprise reports of bycatch in commercial fisheries, media accounts, personal communications, and documentation of shark-human interactions from Russia (n = 8), Republic of Korea (22), Japan (129), China (32), Taiwan (45), Philippines (1) and Vietnam (3). Observations occurred in all months, excluding October-January in the north (Russia and Republic of Korea) and July-August in the south (China, Taiwan, Philippines, and Vietnam). Population trend analysis indicated that the relative abundance of white sharks in the region has remained relatively stable, but parameterization of a 75% increase in observer effort found evidence of a minor decline since 2002. Reliably measured sharks ranged from 126–602 cm total length (TL) and 16–2530 kg total weight. The largest shark in this study (602 cm TL) represents the largest measured shark on record worldwide. For all countries combined the sex ratio was non-significantly biased towards females (1∶1.1; n = 113). Of 60 females examined, 11 were confirmed pregnant ranging from the beginning stages of pregnancy (egg cases) to near term (140 cm TL embryos). On average, 6.0±2.2 embryos were found per litter (maximum of 10) and gestation period was estimated to be 20 months. These observations confirm that white sharks are present in the Northwest Pacific Ocean year-round. While acknowledging the difficulties of studying little known populations of a naturally low abundance species, these results highlight the need for dedicated research to inform regional conservation and management planning.

Effect of sample preparation techniques for carbon and nitrogen stable isotope analysis of hydroxyapatite structures in the form of elasmobranch vertebral centra

RATIONALE

Bulk stable isotope analysis (SIA) provides an important tool for the study of animal ecology. Elasmobranch vertebral centra can be serially sampled to obtain an isotopic history of an individual over ontogeny. The measured total δ(13)C value, however, may be misinterpreted due to the inclusion of the (13)C-rich inorganic portion. Hydrochloric acid (HCl) is commonly used to remove the inorganic portion of hydroxyapatite structures before undertaking SIA, but more recently ethylenediaminetetraacetic acid (EDTA) has been recommended for elasmobranch vertebrae. These acid treatments may introduce uncertainty on measured δ(13)C and δ(15)N values above instrument precision and the effect of small sample size remains untested for elasmobranch vertebrae.

METHODS

Using a non-dilution program on an isotope ratio mass spectrometer the minimum sample weight of vertebrae required to obtain accurate isotopic values was determined for three shark species: white (Carcharodon carcharias), tiger (Galeocerdo cuvier), and sand tiger (Carcharias taurus). To examine if acid treatment completely removes the inorganic component of the vertebrae or whether the technique introduces its own uncertainty on measured δ(13)C and δ(15)N values, vertebrae samples were analyzed untreated and following EDTA treatment.

RESULTS

The minimum sample weight required for accurate stable isotope values and the percentage sample yield following EDTA treatment varied within and among species. After EDTA treatment, white shark vertebrae were all enriched in (13)C and depleted in (15) N, tiger shark vertebrae showed both enrichment and depletion of (13)C and (15)N, and sand tiger shark vertebrae were all depleted in (13)C and (15)N.

CONCLUSIONS

EDTA treatment of elasmobranch vertebrae produces unpredictable effects (i.e. non-linear and non-correctable) among species in both the percentage sample yield and the measured δ(13)C and δ(15)N values. Prior to initiating a large-scale study, we strongly recommend investigating (i) the minimum weight of vertebral material required to obtain consistent isotopic values and (ii) the effects of EDTA treatment, specific to the study species and the isotope ratio mass spectrometer employed.

Rescaling the trophic structure of marine food webs

Measures of trophic position (TP) are critical for understanding food web interactions and human-mediated ecosystem disturbance. Nitrogen stable isotopes (δ(15) N) provide a powerful tool to estimate TP but are limited by a pragmatic assumption that isotope discrimination is constant (change in δ(15) N between predator and prey, Δ(15) N = 3.4‰), resulting in an additive framework that omits known Δ(15) N variation. Through meta-analysis, we determine narrowing discrimination from an empirical linear relationship between experimental Δ(15) N and δ(15) N values of prey consumed. The resulting scaled Δ(15) N framework estimated reliable TPs of zooplanktivores to tertiary piscivores congruent with known feeding relationships that radically alters the conventional structure of marine food webs. Apex predator TP estimates were markedly higher than currently assumed by whole-ecosystem models, indicating perceived food webs have been truncated and species-interactions over simplified. The scaled Δ(15) N framework will greatly improve the accuracy of trophic estimates widely used in ecosystem-based management.

Variable δ(15)N diet-tissue discrimination factors among sharks: implications for trophic position, diet and food web models

The application of stable isotopes to characterize the complexities of a species foraging behavior and trophic relationships is dependent on assumptions of δ¹⁵N diet-tissue discrimination factors (∆¹⁵N). As ∆¹⁵N values have been experimentally shown to vary amongst consumers, tissues and diet composition, resolving appropriate species-specific ∆¹⁵N values can be complex. Given the logistical and ethical challenges of controlled feeding experiments for determining ∆¹⁵N values for large and/or endangered species, our objective was to conduct an assessment of a range of reported ∆¹⁵N values that can hypothetically serve as surrogates for describing the predator-prey relationships of four shark species that feed on prey from different trophic levels (i.e., different mean δ¹⁵N dietary values). Overall, the most suitable species-specific ∆¹⁵N values decreased with increasing dietary-δ¹⁵N values based on stable isotope Bayesian ellipse overlap estimates of shark and the principal prey functional groups contributing to the diet determined from stomach content analyses. Thus, a single ∆¹⁵N value was not supported for this speciose group of marine predatory fishes. For example, the ∆¹⁵N value of 3.7‰ provided the highest percent overlap between prey and predator isotope ellipses for the bonnethead shark (mean diet δ¹⁵N = 9‰) whereas a ∆¹⁵N value < 2.3‰ provided the highest percent overlap between prey and predator isotope ellipses for the white shark (mean diet δ¹⁵N = 15‰). These data corroborate the previously reported inverse ∆¹⁵N-dietary δ¹⁵N relationship when both isotope ellipses of principal prey functional groups and the broader identified diet of each species were considered supporting the adoption of different ∆¹⁵N values that reflect the predators’ δ¹⁵N-dietary value. These findings are critical for refining the application of stable isotope modeling approaches as inferences regarding a species’ ecological role in their community will be influenced with consequences for conservation and management actions.

Population expansion and genetic structure in Carcharhinus brevipinna in the southern Indo-Pacific

BACKGROUND

Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa.

METHODS AND FINDINGS

Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and F ST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717-0.03508, p values ≤ 0.0013; pairwise F ST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia's east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).

CONCLUSIONS

The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

The behavioural and genetic mating system of the sand tiger shark, Carcharias taurus, an intrauterine cannibal

Sand tiger sharks (Carcharias taurus) have an unusual mode of reproduction, whereby the first embryos in each of the paired uteri to reach a certain size ('hatchlings') consume all of their smaller siblings during gestation ('embryonic cannibalism' or EC). If females commonly mate with multiple males ('behavioural polyandry') then litters could initially have multiple sires. It is possible, however, that EC could exclude of all but one of these sires from producing offspring thus influencing the species genetic mating system ('genetic monogamy'). Here, we use microsatellite DNA profiling of mothers and their litters (n = 15, from two to nine embryos per litter) to quantify the frequency of behavioural and genetic polyandry in this system. We conservatively estimate that nine of the females we examined (60%) were behaviourally polyandrous. The genetic mating system was characterized by assessing sibling relationships between hatchlings and revealed only 40 per cent genetic polyandry (i.e. hatchlings were full siblings in 60% of litters). The discrepancy stemmed from three females that were initially fertilized by multiple males but only produced hatchlings with one of them. This reveals that males can be excluded even after fertilizing ova and that some instances of genetic monogamy in this population arise from the reduction in litter size by EC. More research is needed on how cryptic post-copulatory and post-zygotic processes contribute to determining paternity and bridging the behavioural and genetic mating systems of viviparous species.