Population structure and cryptic speciation in bonnethead sharks Sphyrna tiburo in the south-eastern U.S.A. and Caribbean

Population structure and lineage diversification within a small, non-dispersive hammerhead shark species, the bonnethead shark Sphyrna tiburo, was assessed. Sphyrna tiburo is currently described as one continuously distributed species along the Atlantic continental margins of North, Central and South America, but recent genetic analysis of an insular population (Trinidad) suggests the possibility of cryptic speciation. To address this issue S. tiburo were sampled at six sites along c. 6200 km of continuous, continental coastline and from one island location (Grand Bahama) across a discontinuity in their distribution (the Straits of Florida), in order to test if they constitute a single lineage over this distribution. A total of 1030 bp of the mitochondrial control region (CR) was obtained for 239 S. tiburo, revealing 73 distinct haplotypes, high nucleotide diversity (0·01089) and a pair of highly divergent lineages estimated to have separated 3·61-5·62 million years ago. Mitochondrial cytochrome oxidase I and nuclear internal transcribed spacer loci show the same pattern. Divergence is similar within S. tiburo to that observed between established elasmobranch sister species, providing further evidence of cryptic speciation. A global AMOVA based on CR confirms that genetic diversity is primarily partitioned among populations (Φ_ST = 0·828, P < 0·001) because the divergent lineages are almost perfectly segregated between Belize and North America-The Bahamas. An AMOVA consisting solely of the North American and Bahamian samples is also significantly different from zero (Φ_ST = 0·088, P < 0·001) and pairwise F_ST is significantly different between all sites. These findings suggest that S. tiburo comprises a species complex and supports previous research indicating fine population structure, which has implications for fisheries management and biodiversity conservation.

Into the deep: the functionality of mesopelagic excursions by an oceanic apex predator

Comprehension of ecological processes in marine animals requires information regarding dynamic vertical habitat use. While many pelagic predators primarily associate with epipelagic waters, some species routinely dive beyond the deep scattering layer. Actuation for exploiting these aphotic habitats remains largely unknown. Recent telemetry data from oceanic whitetip sharks (Carcharhinus longimanus) in the Atlantic show a strong association with warm waters (>20°C) less than 200 m. Yet, individuals regularly exhibit excursions into the meso‐ and bathypelagic zone. In order to examine deep‐diving behavior in oceanic whitetip sharks, we physically recovered 16 pop‐up satellite archival tags and analyzed the high‐resolution depth and temperature data. Diving behavior was evaluated in the context of plausible functional behavior hypotheses including interactive behaviors, energy conservation, thermoregulation, navigation, and foraging. Mesopelagic excursions (n = 610) occurred throughout the entire migratory circuit in all individuals, with no indication of site specificity. Six depth‐versus‐time descent and ascent profiles were identified. Descent profile shapes showed little association with examined environmental variables. Contrastingly, ascent profile shapes were related to environmental factors and appear to represent unique behavioral responses to abiotic conditions present at the dive apex. However, environmental conditions may not be the sole factors influencing ascents, as ascent mode may be linked to intentional behaviors. While dive functionality remains unconfirmed, our study suggests that mesopelagic excursions relate to active foraging behavior or navigation. Dive timing, prey constituents, and dive shape support foraging as the most viable hypothesis for mesopelagic excursions, indicating that the oceanic whitetip shark may regularly survey extreme environments (deep depths, low temperatures) as a foraging strategy. At the apex of these deep‐water excursions, sharks exhibit a variable behavioral response, perhaps, indicating the presence or absence of prey.

Using genetic inference to re-evaluate the minimum longevity of the lemon shark Negaprion brevirostris

A combination of mark-recapture and genetic sampling was used to extend the minimum longevity of an elasmobranch species and the life span estimate of the lemon shark Negaprion brevirostris was increased conservatively from 20·2 to 37 years. This increase in longevity means higher vulnerability and a longer recovery time from exploitation.

Facultative parthenogenesis in a critically endangered wild vertebrate

Facultative parthenogenesis - the ability of sexually reproducing species to sometimes produce offspring asexually - is known from a wide range of ordinarily sexually reproducing vertebrates in captivity, including some birds, reptiles and sharks [1-3]. Despite this, free-living parthenogens have never been observed in any of these taxa in the wild, although two free-living snakes were recently discovered each gestating a single parthenogen - one copperhead (Agkistrodon contortrix) and one cottonmouth (Agkistrodon piscivorus) [1]. Vertebrate parthenogens are characterized as being of the homogametic sex (e.g., females in sharks, males in birds) and by having elevated homozygosity compared to their mother [1-3], which may reduce their viability [4]. Although it is unknown if either of the parthenogenetic snakes would have been carried to term or survived in the wild, facultative parthenogenesis might have adaptive significance [1]. If this is true, it is reasonable to hypothesize that parthenogenesis would be found most often at low population density, when females risk reproductive failure because finding mates is difficult [5]. Here, we document the first examples of viable parthenogens living in a normally sexually reproducing wild vertebrate, the smalltooth sawfish (Pristis pectinata). We also provide a simple approach to screen any microsatellite DNA database for parthenogens, which will enable hypothesis-driven research on the significance of vertebrate parthenogenesis in the wild.

Genetic structure, population demography and seasonal occurrence of blacktip shark Carcharhinus limbatus in Bimini, the Bahamas

A longline survey was conducted from 2004 to 2014 to investigate the demographic population structure and seasonal abundance of the blacktip shark Carcharhinus limbatus in the Bimini Islands, the Bahamas. All individuals sampled (n = 242) were sub-adult or adults [70·1-145·1 cm pre-caudal length (LPC) range] with no neonates or YOY recorded in Bimini. Carcharhinus limbatus abundance peaked in September, coincident with the largest ratio of female to male sharks and a peak in fresh mating wounds on females. Mitochondrial control region (mtCR) DNA sequences were obtained from C. limbatus at Bimini to test whether Bimini C. limbatus are most closely related to geographically proximate populations sampled on the south-eastern coast of the U.S.A., the closest known nursery areas for this species. Nine mtCR haplotypes were observed in 32 individuals sampled at Bimini [haplotype diversity (h) = 0·821, nucleotide diversity (π) = 0·0015]. Four haplotypes observed from Bimini matched those previously found in the northern Yucatan (Mexico)-Belize and two matched a haplotype previously found in the U.S.A. Four haplotypes were novel but were closely related to the northern Yucatan-Belizean haplotypes. Pair-wise ΦST analysis showed that Bimini was significantly differentiated from all of the populations previously sampled (U.S.A. Atlantic, U.S.A. Gulf of Mexico, northern Yucatan, Belize and Brazil). This indicates that C. limbatus sampled from Bimini are unlikely from the described, proximate U.S.A. nurseries.

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.

A quantitative genetic approach to assess the evolutionary potential of a coastal marine fish to ocean acidification

Assessing the potential of marine organisms to adapt genetically to increasing oceanic CO₂ levels requires proxies such as heritability of fitness-related traits under ocean acidification (OA). We applied a quantitative genetic method to derive the first heritability estimate of survival under elevated CO₂ conditions in a metazoan. Specifically, we reared offspring, selected from a wild coastal fish population (Atlantic silverside, Menidia menidia), at high CO₂ conditions (∼2300 μatm) from fertilization to 15 days posthatch, which significantly reduced survival compared to controls. Perished and surviving offspring were quantitatively sampled and genotyped along with their parents, using eight polymorphic microsatellite loci, to reconstruct a parent–offspring pedigree and estimate variance components. Genetically related individuals were phenotypically more similar (i.e., survived similarly long at elevated CO₂ conditions) than unrelated individuals, which translated into a significantly nonzero heritability (0.20 ± 0.07). The contribution of maternal effects was surprisingly small (0.05 ± 0.04) and nonsignificant. Survival among replicates was positively correlated with genetic diversity, particularly with observed heterozygosity. We conclude that early life survival of M. menidia under high CO₂ levels has a significant additive genetic component that could elicit an evolutionary response to OA, depending on the strength and direction of future selection.

A novel mini-DNA barcoding assay to identify processed fins from internationally protected shark species

There is a growing need to identify shark products in trade, in part due to the recent listing of five commercially important species on the Appendices of the Convention on International Trade in Endangered Species (CITES; porbeagle, Lamna nasus, oceanic whitetip, Carcharhinus longimanus scalloped hammerhead, Sphyrna lewini, smooth hammerhead, S. zygaena and great hammerhead S. mokarran) in addition to three species listed in the early part of this century (whale, Rhincodon typus, basking, Cetorhinus maximus, and white, Carcharodon carcharias). Shark fins are traded internationally to supply the Asian dried seafood market, in which they are used to make the luxury dish shark fin soup. Shark fins usually enter international trade with their skin still intact and can be identified using morphological characters or standard DNA-barcoding approaches. Once they reach Asia and are traded in this region the skin is removed and they are treated with chemicals that eliminate many key diagnostic characters and degrade their DNA ("processed fins"). Here, we present a validated mini-barcode assay based on partial sequences of the cytochrome oxidase I gene that can reliably identify the processed fins of seven of the eight CITES listed shark species. We also demonstrate that the assay can even frequently identify the species or genus of origin of shark fin soup (31 out of 50 samples).

There and back again: a review of residency and return migrations in sharks, with implications for population structure and management

The overexploitation of sharks has become a global environmental issue in need of a comprehensive and multifaceted management response. Tracking studies are beginning to elucidate how shark movements shape the internal dynamics and structure of populations, which determine the most appropriate scale of these management efforts. Tracked sharks frequently either remain in a restricted geographic area for an extended period of time (residency) or return to a previously resided-in area after making long-distance movements (site fidelity). Genetic studies have shown that some individuals of certain species preferentially return to their exact birthplaces (natal philopatry) or birth regions (regional philopatry) for either parturition or mating, even though they make long-distance movements that would allow them to breed elsewhere. More than 80 peer-reviewed articles, constituting the majority of published shark tracking and population genetic studies, provide evidence of at least one of these behaviors in a combined 31 shark species from six of the eight extant orders. Residency, site fidelity, and philopatry can alone or in combination structure many coastal shark populations on finer geographic scales than expected based on their potential for dispersal. This information should therefore be used to scale and inform assessment, management, and conservation activities intended to restore depleted shark populations.

Methylmercury in dried shark fins and shark fin soup from American restaurants

Consumption of meat from large predatory sharks exposes human consumers to high levels of toxic monomethylmercury (MMHg). There also have been claims that shark fins, and hence the Asian delicacy shark fin soup, contain harmful levels of neurotoxic chemicals in combination with MMHg, although concentrations of MMHg in shark fins are unknown. We measured MMHg in dried, unprocessed fins (n=50) of 13 shark species that occur in the international trade of dried shark fins as well as 50 samples of shark fin soup prepared by restaurants from around the United States. Concentrations of MMHg in fins ranged from 9 to 1720 ng/g dry wt. MMHg in shark fin soup ranged from <0.01 to 34 ng/mL, with MMHg averaging 62 ± 7% of total Hg. The highest concentrations of MMHg and total Hg were observed in both fins and soup from large, high trophic level sharks such as hammerheads (Sphyrna spp.). Consumption of a 240 mL bowl of shark fin soup containing the average concentration of MMHg (4.6 ng/mL) would result in a dose of 1.1 μg MMHg, which is 16% of the U.S. EPA's reference dose (0.1 μg MMHg per 1 kg per day in adults) of 7.4 μg per day for a 74 kg person. If consumed, the soup containing the highest measured MMHg concentration would exceed the reference dose by 17%. While shark fin soup represents a potentially important source of MMHg to human consumers, other seafood products, particularly the flesh of apex marine predators, contain much higher MMHg concentrations and can result in substantially greater exposures of this contaminant for people.

Regional-scale migrations and habitat use of juvenile lemon sharks (Negaprion brevirostris) in the US South Atlantic

Resolving the geographic extent and timing of coastal shark migrations, as well as their environmental cues, is essential for refining shark management strategies in anticipation of increasing anthropogenic stressors to coastal ecosystems. We employed a regional-scale passive acoustic telemetry array encompassing 300 km of the east Florida coast to assess what factors influence site fidelity of juvenile lemon sharks (Negaprion brevirostris) to an exposed coastal nursery at Cape Canaveral, and to document the timing and rate of their seasonal migrations. Movements of 54 juvenile lemon sharks were monitored for three years with individuals tracked for up to 751 days. While most sharks demonstrated site fidelity to the Cape Canaveral region December through February under typical winter water temperatures, historically extreme declines in ocean temperature were accompanied by rapid and often temporary, southward displacements of up to 190 km along the Florida east coast. From late February through April each year, most sharks initiated a northward migration at speeds of up to 64 km day⁻¹ with several individuals then detected in compatible estuarine telemetry arrays in Georgia and South Carolina up to 472 km from release locations. Nineteen sharks returned for a second or even third consecutive winter, thus demonstrating strong seasonal philopatry to the Cape Canaveral region. The long distance movements and habitat associations of immature lemon sharks along the US southeast coast contrast sharply with the natal site fidelity observed in this species at other sites in the western Atlantic Ocean. These findings validate the existing multi-state management strategies now in place. Results also affirm the value of collaborative passive arrays for resolving seasonal movements and habitat preferences of migratory coastal shark species not easily studied with other tagging techniques.

Severe inbreeding and small effective number of breeders in a formerly abundant marine fish

In contrast to freshwater fish it is presumed that marine fish are unlikely to spawn with close relatives due to the dilution effect of large breeding populations and their propensity for movement and reproductive mixing. Inbreeding is therefore not typically a focal concern of marine fish management. We measured the effective number of breeders in 6 New York estuaries for winter flounder (Pseudopleuronectes americanus), a formerly abundant fish, using 11 microsatellite markers (6–56 alleles per locus). The effective number of breeders for 1–2 years was remarkably small, with point estimates ranging from 65–289 individuals. Excess homozygosity was detected at 10 loci in all bays (F_IS = 0.169–0.283) and individuals exhibited high average internal relatedness (IR; mean = 0.226). These both indicate that inbreeding is very common in all bays, after testing for and ruling out alternative explanations such as technical and sampling artifacts. This study demonstrates that even historically common marine fish can be prone to inbreeding, a factor that should be considered in fisheries management and conservation plans.

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.

Supplemental feeding for ecotourism reverses diel activity and alters movement patterns and spatial distribution of the southern stingray, Dasyatis americana

Southern stingrays, Dasyatis americana, have been provided supplemental food in ecotourism operations at Stingray City Sandbar (SCS), Grand Cayman since 1986, with this site becoming one of the world’s most famous and heavily visited marine wildlife interaction venues. Given expansion of marine wildlife interactive tourism worldwide, there are questions about the effects of such activities on the focal species and their ecosystems. We used a combination of acoustic telemetry and tag-recapture efforts to test the hypothesis that human-sourced supplemental feeding has altered stingray activity patterns and habitat use at SCS relative to wild animals at control sites. Secondarily, we also qualitatively estimated the population size of stingrays supporting this major ecotourism venue. Tag-recapture data indicated that a population of at least 164 stingrays, over 80% female, utilized the small area at SCS for prolonged periods of time. Examination of comparative movements of mature female stingrays at SCS and control sites revealed strong differences between the two groups: The fed animals demonstrated a notable inversion of diel activity, being constantly active during the day with little movement at night compared to the nocturnally active wild stingrays; The fed stingrays utilized significantly (p<0.05) smaller 24 hour activity spaces compared to wild conspecifics, staying in close proximity to the ecotourism site; Fed stingrays showed a high degree of overlap in their core activity spaces compared to wild stingrays which were largely solitary in the spaces utilized (72% vs. 3% overlap respectively). Supplemental feeding has strikingly altered movement behavior and spatial distribution of the stingrays, and generated an atypically high density of animals at SCS which could have downstream fitness costs for individuals and potentially broader ecosystem effects. These findings should help environmental managers plan mitigating measures for existing operations, and develop precautionary policies regarding proposed feeding sites.

A streamlined DNA tool for global identification of heavily exploited coastal shark species (genus Rhizoprionodon)

Obtaining accurate species-specific landings data is an essential step toward achieving sustainable shark fisheries. Globally distributed sharpnose sharks (genus Rhizoprionodon) exhibit life-history characteristics (rapid growth, early maturity, annual reproduction) that suggests that they could be fished in a sustainable manner assuming an investment in monitoring, assessment and careful management. However, obtaining species-specific landings data for sharpnose sharks is problematic because they are morphologically very similar to one another. Moreover, sharpnose sharks may also be confused with other small sharks (either small species or juveniles of large species) once they are processed (i.e., the head and fins are removed). Here we present a highly streamlined molecular genetics approach based on seven species-specific PCR primers in a multiplex format that can simultaneously discriminate body parts from the seven described sharpnose shark species commonly occurring in coastal fisheries worldwide. The species-specific primers are based on nucleotide sequence differences among species in the nuclear ribosomal internal transcribed spacer 2 locus (ITS2). This approach also distinguishes sharpnose sharks from a wide range of other sharks (52 species) and can therefore assist in the regulation of coastal shark fisheries around the world.

Stable isotopes and elasmobranchs: tissue types, methods, applications and assumptions

Stable-isotope analysis (SIA) can act as a powerful ecological tracer with which to examine diet, trophic position and movement, as well as more complex questions pertaining to community dynamics and feeding strategies or behaviour among aquatic organisms. With major advances in the understanding of the methodological approaches and assumptions of SIA through dedicated experimental work in the broader literature coupled with the inherent difficulty of studying typically large, highly mobile marine predators, SIA is increasingly being used to investigate the ecology of elasmobranchs (sharks, skates and rays). Here, the current state of SIA in elasmobranchs is reviewed, focusing on available tissues for analysis, methodological issues relating to the effects of lipid extraction and urea, the experimental dynamics of isotopic incorporation, diet-tissue discrimination factors, estimating trophic position, diet and mixing models and individual specialization and niche-width analyses. These areas are discussed in terms of assumptions made when applying SIA to the study of elasmobranch ecology and the requirement that investigators standardize analytical approaches. Recommendations are made for future SIA experimental work that would improve understanding of stable-isotope dynamics and advance their application in the study of sharks, skates and rays.

Reef sharks exhibit site-fidelity and higher relative abundance in marine reserves on the Mesoamerican Barrier Reef

Carcharhinid sharks can make up a large fraction of the top predators inhabiting tropical marine ecosystems and have declined in many regions due to intense fishing pressure. There is some support for the hypothesis that carcharhinid species that complete their life-cycle within coral reef ecosystems, hereafter referred to as “reef sharks”, are more abundant inside no-take marine reserves due to a reduction in fishing pressure (i.e., they benefit from marine reserves). Key predictions of this hypothesis are that (a) individual reef sharks exhibit high site-fidelity to these protected areas and (b) their relative abundance will generally be higher in these areas compared to fished reefs. To test this hypothesis for the first time in Caribbean coral reef ecosystems we combined acoustic monitoring and baited remote underwater video (BRUV) surveys to measure reef shark site-fidelity and relative abundance, respectively. We focused on the Caribbean reef shark (Carcharhinus perezi), the most common reef shark in the Western Atlantic, at Glover's Reef Marine Reserve (GRMR), Belize. Acoustically tagged sharks (N = 34) were detected throughout the year at this location and exhibited strong site-fidelity. Shark presence or absence on 200 BRUVs deployed at GRMR and three other sites (another reserve site and two fished reefs) showed that the factor “marine reserve” had a significant positive effect on reef shark presence. We rejected environmental factors or site-environment interactions as predominant drivers of this pattern. These results are consistent with the hypothesis that marine reserves can benefit reef shark populations and we suggest new hypotheses to determine the underlying mechanism(s) involved: reduced fishing mortality or enhanced prey availability.

Genetic diversity despite population collapse in a critically endangered marine fish: the smalltooth sawfish (Pristis pectinata)

Sawfish (family Pristidae) are among the most critically endangered marine fish in the world, yet very little is known about how genetic bottlenecks, genetic drift, and inbreeding depression may be affecting these elasmobranchs. In the US Atlantic, the smalltooth sawfish (Pristis pectinata) has declined to 1-5% of its abundance in the 1900s, and its core distribution has contracted to southwest Florida. We used 8 polymorphic microsatellite markers to show that this remnant population still exhibits high genetic diversity in terms of average allelic richness (18.23), average alleles per locus (18.75, standard deviation [SD] 6.6) and observed heterozygosity (0.43-0.98). Inbreeding is rare (mean individual internal relatedness = -0.02, SD 0.14; F(IS) = -0.011, 95% confidence interval [CI] = -0.039 to 0.011), even though the estimated effective population size (N(e)) is modest (250-350, 95% CI = 142-955). Simulations suggest that the remnant smalltooth sawfish population will probably retain >90% of its current genetic diversity over the next century even at the lower estimate of N(e). There is no evidence of a genetic bottleneck accompanying last century's demographic bottleneck, and we discuss hypotheses that could explain this. We also discuss features of elasmobranch life history and population biology that could make them less vulnerable than other large marine vertebrates to genetic change associated with reduced population size.