Comparative population genomics of manta rays has global implications for management.

Understanding population connectivity and genetic diversity is of fundamental importance to conservation. However, in globally threatened marine megafauna, challenges remain due to their elusive nature and wide-ranging distributions. As overexploitation continues to threaten biodiversity across the globe, such knowledge gaps compromise both the suitability and effectiveness of management actions. Here, we use a comparative framework to investigate genetic differentiation and diversity of manta rays, one of the most iconic yet vulnerable groups of elasmobranchs on the planet. Despite their recent divergence, we show how oceanic manta rays (Mobula birostris) display significantly higher heterozygosity than reef manta rays (Mobula alfredi) and that M. birostris populations display higher connectivity worldwide. Through inferring modes of colonization, we reveal how both contemporary and historical forces have likely influenced these patterns, with important implications for population management. Our findings highlight the potential for fisheries to disrupt population dynamics at both local and global scales and therefore have direct relevance for international conservation of marine species.

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 δ13 C and δ15 N values (i.e., result in artificially lower δ13 C and δ15 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 δ13 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 δ15 N values between lipid-extracted (LE) liver and lipid-extracted/water washed (WW) treatments, however, WW resulted in significant increases in %N, δ13 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 liverLEWW tissue, respectively). The preliminary comparison of C:N threshold corrected liver and muscle δ13 C values corresponded with movement/habitat behaviours for each shark; minor differences in δ13 C values were observed for known regional movements of DUS and RAG (δ13 CDiffs = 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 δ13 C and δ15 N values to examine elasmobranch ecology.

Phylogenomics and species delimitation for effective conservation of manta and devil rays.

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Are you really what you eat? Stomach content analysis and stable isotope ratios do not uniformly estimate dietary niche characteristics in three marine predators.

Calculation of dietary niche characteristics using stable isotopes has become a popular approach to understand the functional role of taxa across food webs. An underlying assumption of this approach is that stable isotopes accurately reflect the dietary breadth of a species over a temporal duration defined by tissue-specific isotopic turnover rates. In theory, dietary niche estimates derived from fast turnover rate tissues (e.g., blood plasma and liver) may augment stomach content-derived estimates more agreeably than slower turnover rate tissues (e.g., muscle or fin). We tested this hypothesis by comparing commonly used dietary niche estimates derived from stomach contents (nicheSCA: Levins', Shannon-Wiener's, and Smith's), with those estimated using stable isotopes [nicheSIA: standard ellipse area (SEA), convex hull total area (TA), theta (θ), and ellipse eccentricity (E)] of liver and muscle tissue. Model species were three large-bodied sharks: white (Carcharodon carcharias), dusky (Carcharhinus obscurus), and scalloped hammerhead (Sphyrna lewini). Within-technique comparisons for nicheSCA and nicheSIA metrics (i.e., SEA vs. TA) were often correlated; however, we did not observe any statistically significant correlations between nicheSCA and liver/muscle tissue nicheSIA (i.e., Levins' vs. SEA). We conclude that nicheSCA and nicheSIA do not provide comparable estimates of dietary niche, at least for the three predator species examined. This fundamental discrepancy highlights technique-specific limitations to estimating organismal dietary niche and identifies a need for the use of clearly defined niche metrics, i.e., the standardized use and reporting of the term isotopic niche as proposed by Newsome et al. (Front Ecol Environ 5:429-436, 2007). Finally, further investigation into the factors underpinning nicheSIA is required to better contextualize this popular ecological metric when compared to nicheSCA.

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.

Global population genetic dynamics of a highly migratory, apex predator shark.

Knowledge of genetic connectivity dynamics in the world's large-bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top-down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo-Pacific ocean basins (microsatellite FST  > 0.129; CR ΦST  > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex-biased gene flow. Coalescent- and genetic distance-based estimates of divergence from CR sequences were largely congruent (dcorr  = 0.0015-0.0050), indicating a separation of Indo-Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller-than-anticipated spatial scales.

A dated molecular phylogeny of manta and devil rays (Mobulidae) based on mitogenome and nuclear sequences.

Manta and devil rays are an iconic group of globally distributed pelagic filter feeders, yet their evolutionary history remains enigmatic. We employed next generation sequencing of mitogenomes for nine of the 11 recognized species and two outgroups; as well as additional Sanger sequencing of two mitochondrial and two nuclear genes in an extended taxon sampling set. Analysis of the mitogenome coding regions in a Maximum Likelihood and Bayesian framework provided a well-resolved phylogeny. The deepest divergences distinguished three clades with high support, one containing Manta birostris, Manta alfredi, Mobula tarapacana, Mobula japanica and Mobula mobular; one containing Mobula kuhlii, Mobula eregoodootenkee and Mobula thurstoni; and one containing Mobula munkiana, Mobula hypostoma and Mobula rochebrunei. Mobula remains paraphyletic with the inclusion of Manta, a result that is in agreement with previous studies based on molecular and morphological data. A fossil-calibrated Bayesian random local clock analysis suggests that mobulids diverged from Rhinoptera around 30 Mya. Subsequent divergences are characterized by long internodes followed by short bursts of speciation extending from an initial episode of divergence in the Early and Middle Miocene (19-17 Mya) to a second episode during the Pliocene and Pleistocene (3.6 Mya - recent). Estimates of divergence dates overlap significantly with periods of global warming, during which upwelling intensity - and related high primary productivity in upwelling regions - decreased markedly. These periods are hypothesized to have led to fragmentation and isolation of feeding regions leading to possible regional extinctions, as well as the promotion of allopatric speciation. The closely shared evolutionary history of mobulids in combination with ongoing threats from fisheries and climate change effects on upwelling and food supply, reinforces the case for greater protection of this charismatic family of pelagic filter feeders.

Genetic Diversity of White Sharks, Carcharodon carcharias, in the Northwest Atlantic and Southern Africa.

The white shark, Carcharodon carcharias, is both one of the largest apex predators in the world and among the most heavily protected marine fish. Population genetic diversity is in part shaped by recent demographic history and can thus provide information complementary to more traditional population assessments, which are difficult to obtain for white sharks and have at times been controversial. Here, we use the mitochondrial control region and 14 nuclear-encoded microsatellite loci to assess white shark genetic diversity in 2 regions: the Northwest Atlantic (NWA, N = 35) and southern Africa (SA, N = 131). We find that these 2 regions harbor genetically distinct white shark populations (Φ ST = 0.10, P < 0.00001; microsatellite F ST = 0.1057, P < 0.021). M-ratios were low and indicative of a genetic bottleneck in the NWA (M-ratio = 0.71, P < 0.004) but not SA (M-ratio = 0.85, P = 0.39). This is consistent with other evidence showing a steep population decline occurring in the mid to late 20th century in the NWA, whereas the SA population appears to have been relatively stable. Estimates of effective population size ranged from 22.6 to 66.3 (NWA) and 188 to 1998.3 (SA) and evidence of inbreeding was found (primarily in NWA). Overall, our findings indicate that white population dynamics within NWA and SA are determined more by intrinsic reproduction than immigration and there is genetic evidence of a population decline in the NWA, further justifying the strong domestic protective measures that have been taken for this species in this region. Our study also highlights how assessment of genetic diversity can complement other sources of information to better understand the status of threatened marine fish populations.

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.

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.

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.

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.

Maternal investment and size-specific reproductive output in carcharhinid sharks.

1. Life-history theory predicts that organisms will provide an optimal level of parental investment for offspring survival balanced against the effects on their own survival and future reproductive potential. 2. Optimal resource allocation models also predict an increase in reproductive output with age as expected future reproductive effort decreases. To date, maternal investment in sharks has received limited attention. 3. We found that neonatal dusky sharks (Carcharhinus obscurus) are not independent from maternal resource allocation at the point of parturition but instead are provisioned with energy reserves in the form of an enlarged liver that constitutes approximately 20% of total body mass. 4. Analysis of long-term archived data sets showed that a large proportion of this enlarged liver is utilized during the first weeks or months of life suggesting that the reported weight loss of newborn sharks signifies a natural orientation process and is not necessarily related to prey abundance and/or indicative of high mortality rates. 5. Interrogation of near-term pup mass in two carcharhinids, the dusky and spinner shark (Carcharhinus brevipinna), further revealed an increase in reproductive output with maternal size, with evidence for a moderate decline in the largest mothers. 6. For the dusky shark, there was a trade-off between increasing litter size and near-term pup mass in support of optimal offspring size theory. 7. For both the dusky and spinner shark, there was a linear increase in near-term pup mass with month, which may indicate variable parturition strategies and/or that carcharhinids are able to adjust the length of the gestation period. 8. The identification of optimal size-specific reproductive output has direct implications for improving the reproductive potential of exploited shark populations and for structuring future management strategies.

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.

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.

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.

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 δ(15)N diet-tissue discrimination factors (∆(15)N). As ∆(15)N values have been experimentally shown to vary amongst consumers, tissues and diet composition, resolving appropriate species-specific ∆(15)N values can be complex. Given the logistical and ethical challenges of controlled feeding experiments for determining ∆(15)N values for large and/or endangered species, our objective was to conduct an assessment of a range of reported ∆(15)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 δ(15)N dietary values). Overall, the most suitable species-specific ∆(15)N values decreased with increasing dietary-δ(15)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 ∆(15)N value was not supported for this speciose group of marine predatory fishes. For example, the ∆(15)N value of 3.7‰ provided the highest percent overlap between prey and predator isotope ellipses for the bonnethead shark (mean diet δ(15)N = 9‰) whereas a ∆(15)N value < 2.3‰ provided the highest percent overlap between prey and predator isotope ellipses for the white shark (mean diet δ(15)N = 15‰). These data corroborate the previously reported inverse ∆(15)N-dietary δ(15)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 ∆(15)N values that reflect the predators' δ(15)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.