Ecological roles and importance of sharks in the Anthropocene Ocean.

In ecosystems, sharks can be predators, competitors, facilitators, nutrient transporters, and food. However, overfishing and other threats have greatly reduced shark populations, altering their roles and effects on ecosystems. We review these changes and implications for ecosystem function and management. Macropredatory sharks are often disproportionately affected by humans but can influence prey and coastal ecosystems, including facilitating carbon sequestration. Like terrestrial predators, sharks may be crucial to ecosystem functioning under climate change. However, large ecosystem effects of sharks are not ubiquitous. Increasing human uses of oceans are changing shark roles, necessitating management consideration. Rebuilding key populations and incorporating shark ecological roles, including less obvious ones, into management efforts are critical for retaining sharks' functional value. Coupled social-ecological frameworks can facilitate these efforts.

Intraguild processes drive space-use patterns in a large-bodied marine predator community.

Interspecific interactions, including predator-prey, intraguild predation (IGP) and competition, may drive distribution and habitat use of predator communities. However, elucidating the relative importance of these interactions in shaping predator distributions is challenging, especially in marine communities comprising highly mobile species. We used individual-based models (IBMs) to predict the habitat distributions of apex predators, intraguild (IG) prey and prey. We then used passive acoustic telemetry to test these predictions in a subtropical marine predator community consisting of eight elasmobranch (i.e. shark and ray) species in Bimini, The Bahamas. IBMs predicted that prey and IG prey will preferentially select habitats based on safety over resources (food), with stronger selection for safe habitat by smaller prey. Elasmobranch space-use patterns matched these predictions. Species with predator-prey and asymmetrical IGP (between apex and small mesopredators) interactions showed the clearest spatial separation, followed by asymmetrical IGP among apex and large mesopredators. Competitors showed greater spatial overlap although with finer-scale differences in microhabitat use. Our study suggests space-use patterns in elasmobranchs are at least partially driven by interspecific interactions, with stronger spatial separation occurring where interactions include predator-prey relationships or IGP.

Directed conservation of the world's reef sharks and rays.

Many shark populations are in decline around the world, with severe ecological and economic consequences. Fisheries management and marine protected areas (MPAs) have both been heralded as solutions. However, the effectiveness of MPAs alone is questionable, particularly for globally threatened sharks and rays ('elasmobranchs'), with little known about how fisheries management and MPAs interact to conserve these species. Here we use a dedicated global survey of coral reef elasmobranchs to assess 66 fully protected areas embedded within a range of fisheries management regimes across 36 countries. We show that conservation benefits were primarily for reef-associated sharks, which were twice as abundant in fully protected areas compared with areas open to fishing. Conservation benefits were greatest in large protected areas that incorporate distinct reefs. However, the same benefits were not evident for rays or wide-ranging sharks that are both economically and ecologically important while also threatened with extinction. We show that conservation benefits from fully protected areas are close to doubled when embedded within areas of effective fisheries management, highlighting the importance of a mixed management approach of both effective fisheries management and well-designed fully protected areas to conserve tropical elasmobranch assemblages globally.

Widespread diversity deficits of coral reef sharks and rays.

A global survey of coral reefs reveals that overfishing is driving resident shark species toward extinction, causing diversity deficits in reef elasmobranch (shark and ray) assemblages. Our species-level analysis revealed global declines of 60 to 73% for five common resident reef shark species and that individual shark species were not detected at 34 to 47% of surveyed reefs. As reefs become more shark-depleted, rays begin to dominate assemblages. Shark-dominated assemblages persist in wealthy nations with strong governance and in highly protected areas, whereas poverty, weak governance, and a lack of shark management are associated with depauperate assemblages mainly composed of rays. Without action to address these diversity deficits, loss of ecological function and ecosystem services will increasingly affect human communities.

An apex predator engineers wetland food-web heterogeneity through nutrient enrichment and habitat modification.

The potential for animals to modify spatial patterns of nutrient limitation for autotrophs and habitat availability for other members of their communities is increasingly recognized. However, net trophic effects of consumers acting as ecosystem engineers remain poorly known. The American Alligator Alligator mississippiensis is an abundant predator capable of dramatic modifications of physical habitat through the creation and maintenance of pond-like basins, but its role in influencing community structure and nutrient dynamics is less appreciated. We investigated if alligators engineer differences in nutrient availability and changes to community structure by their creation of 'alligator ponds' compared to the surrounding phosphorus (P)-limited oligotrophic marsh. We used a halo sampling design of three distinct habitats extending outward from 10 active alligator ponds across a hydrological gradient in the Everglades, USA. We performed nutrient analysis on basal food-web resources and quantitative community analyses, and stoichiometric analyses on plants and animals. Our findings demonstrate that alligators act as ecosystem engineers and enhance food-web heterogeneity by increasing nutrient availability, manipulating physical structure and altering algal, plant and animal communities. Flocculent detritus, an unconsolidated layer of particulate organic matter and soil, showed strong patterns of P enrichment in ponds. Higher P availability in alligator ponds also resulted in bottom-up trophic transfer of nutrients as evidenced by higher growth rates (lower N:P) for plants and aquatic consumers. Edge habitats surrounding alligator ponds contained the most diverse communities of invertebrates and plants, but low total abundance of fishes, likely driven by high densities of emergent macrophytes. Pond communities exhibited higher abundance of fish compared to edge habitat and were dominated by compositions of small invertebrates that track high nutrient availability in the water column. Marshes contained high numbers of animals that are closely tied to periphyton mats, which were absent from other habitats. Alligator-engineered habitats are ecologically important by providing nutrient-enriched 'hotspots' in an oligotrophic system, habitat heterogeneity to marshes, and refuges for other fauna during seasonal disturbances. This work adds to growing evidence that efforts to model community dynamics should routinely consider animal-mediated bottom-up processes like ecosystem engineering.

Critically endangered Rice's whales (Balaenoptera ricei) selectively feed on high-quality prey in the Gulf of Mexico.

Determining the drivers of prey selection in marine predators is critical when investigating ecosystem structure and function. The newly recognized Rice's whale (Balaenoptera ricei) is one of the most critically endangered large whales in the world and endemic to the industrialized Gulf of Mexico. Here, we investigated the drivers of resource selection by Rice's whales in relation to prey availability and energy density. Bayesian stable isotope (δ13C, δ15N) mixing models suggest that Rice's whales feed primarily on a schooling fish, Ariomma bondi (66.8% relative contribution). Prey selection using the Chesson's index revealed that active prey selection was found to be positive for three out of the four potential prey identified in the mixing model. A low degree of overlap between prey availability and diet inferred from the mixing model (Pianka Index: 0.333) suggests that prey abundance is not the primary driver of prey selection. Energy density data suggest that prey selection may be primarily driven by the energy content. Results from this study indicate that Rice's whales are selective predators consuming schooling prey with the highest energy content. Environmental changes in the region have the potential to influence prey species that would make them less available to Rice's whales.

Elucidating the role of competition in driving spatial and trophic niche patterns in sympatric juvenile sharks.

The coexistence of ecologically and morphologically similar species is often facilitated by the partitioning of ecological niches. While subordinate species can reduce competition with dominant competitors through spatial and/or trophic segregation, empirical support from wild settings, particularly those involving large-bodied taxa in marine ecosystems, are rare. Shark nursery areas provide an opportunity to investigate the mechanisms of coexistence. We used experimental and field studies of sympatric juvenile sharks (blacktip reef shark, Carcharhinus melanopterus; sicklefin lemon shark, Negaprion acutidens) to investigate how competitive ability influenced realized niches at St. Joseph Atoll, Seychelles. Captive trials revealed that sicklefin lemon sharks were dominant over blacktip reef sharks, consistently taking food rewards. In the field, blacktip reef sharks were captured over a broader area than sicklefin lemon sharks, but daily space use of actively tracked sharks showed a high degree of overlap across microhabitats. While stomach contents analysis revealed that blacktip reef shark diets included a broader range of prey items, stable isotope analysis demonstrated significantly higher mean δ13C values for sicklefin lemon sharks, suggesting diverging dietary preferences. Overall, our results matched theoretical predictions of subordinate competitors using a greater range of habitats and displaying broader feeding niches than competitively dominant species. While separating the realized and fundamental niche of marine predators is complicated, we provide evidence that resource partitioning is at least partially driven by interspecific competition.

A general pattern of trade-offs between ecosystem resistance and resilience to tropical cyclones.

Tropical cyclones drive coastal ecosystem dynamics, and their frequency, intensity, and spatial distribution are predicted to shift with climate change. Patterns of resistance and resilience were synthesized for 4138 ecosystem time series from n = 26 storms occurring between 1985 and 2018 in the Northern Hemisphere to predict how coastal ecosystems will respond to future disturbance regimes. Data were grouped by ecosystems (fresh water, salt water, terrestrial, and wetland) and response categories (biogeochemistry, hydrography, mobile biota, sedentary fauna, and vascular plants). We observed a repeated pattern of trade-offs between resistance and resilience across analyses. These patterns are likely the outcomes of evolutionary adaptation, they conform to disturbance theories, and they indicate that consistent rules may govern ecosystem susceptibility to tropical cyclones.

Conservation implications of forage base requirements of a marine predator population at carrying capacity.

Prey depletion may contribute to marine predator declines, yet the forage base required to sustain an unfished population of predatory fish at carrying capacity is unknown. We integrated demographic and physiological data within a Bayesian bioenergetic model to estimate annual consumption of a gray reef shark (Carcharhinus amblyrhynchos) population at a remote Pacific atoll (Palmyra Atoll) that are at carrying capacity. Furthermore, we estimated the proportion of the atoll's reef fish biomass production consumed by the gray reef sharks, assuming sharks either partially foraged pelagically (mean 7%), or solely within the reef environment (mean 52%). We then predicted the gray reef shark population potential of other, less remote Pacific Ocean coral reef islands, illustrating that current populations are substantially smaller than could be supported by their forage base. Our research highlights the utility of modeling how far predator population sizes are from their expected carrying capacity in informing marine conservation.

Moray eels are more common on coral reefs subject to higher human pressure in the greater Caribbean.

Proximity and size of the nearest market ('market gravity') have been shown to have strong negative effects on coral reef fish communities that can be mitigated by the establishment of closed areas. However, moray eels are functionally unique predators that are generally not subject to targeted fishing and should therefore not directly be affected by these factors. We used baited remote underwater video systems to investigate associations between morays and anthropogenic, habitat, and ecological factors in the Caribbean region. Market gravity had a positive effect on morays, while the opposite pattern was observed in a predator group subject to exploitation (sharks). Environmental DNA analyses corroborated the positive effect of market gravity on morays. We hypothesize that the observed pattern could be the indirect result of the depletion of moray competitors and predators near humans.

Loss of predation risk from apex predators can exacerbate marine tropicalization caused by extreme climatic events.

Extreme climatic events (ECEs) and predator removal represent some of the most widespread stressors to ecosystems. Though species interactions can alter ecological effects of climate change (and vice versa), it is less understood whether, when and how predator removal can interact with ECEs to exacerbate their effects. Understanding the circumstances under which such interactions might occur is critical because predator loss is widespread and ECEs can generate rapid phase shifts in ecosystems which can ultimately lead to tropicalization. Our goal was to determine whether loss of predation risk may be an important mechanism governing ecosystem responses to extreme events, and whether the effects of such events, such as tropicalization, can occur even when species range shifts do not. Specifically, our goal was to experimentally simulate the loss of an apex predator, the tiger shark Galeocerdo cuvier effects on a recently damaged seagrass ecosystem of Shark Bay, Australia by applying documented changes to risk-sensitive grazing of dugong Dugong dugon herbivores. Using a 16-month-field experiment established in recently disturbed seagrass meadows, we used previous estimates of risk-sensitive dugong foraging behaviour to simulate altered risk-sensitive foraging densities and strategies of dugongs consistent with apex predator loss, and tracked seagrass responses to the simulated grazing. Grazing treatments targeted and removed tropical seagrasses, which declined. However, like in other mixed-bed habitats where dugongs forage, treatments also incidentally accelerated temperate seagrass losses, revealing that herbivore behavioural changes in response to predator loss can exacerbate ECE and promote tropicalization, even without range expansions or introductions of novel species. Our results suggest that changes to herbivore behaviours triggered by loss of predation risk can undermine ecological resilience to ECEs, particularly where long-lived herbivores are abundant. By implication, ongoing losses of apex predators may combine with increasingly frequent ECEs to amplify climate change impacts across diverse ecosystems and large spatial scales.

Elucidating shark diets with DNA metabarcoding from cloacal swabs.

Animal dietary information provides the foundation for understanding trophic relationships, which is essential for ecosystem management. Yet, in marine systems, high-resolution diet reconstruction tools are currently under-developed. This is particularly pertinent for large marine vertebrates, for which direct foraging behaviour is difficult or impossible to observe and, due to their conservation status, the collection of stomach contents at adequate sample sizes is frequently impossible. Consequently, the diets of many groups, such as sharks, have largely remained unresolved. To address this knowledge gap, we applied metabarcoding to prey DNA in faecal residues (fDNA) collected on cotton swabs from the inside of a shark's cloaca. We used a previously published primer set targeting a small section of the 12S rRNA mitochondrial gene to amplify teleost prey species DNA. We tested the utility of this method in a controlled feeding experiment with captive juvenile lemon sharks (Negaprion brevirostris) and on free-ranging juvenile bull sharks (Carcharhinus leucas). In the captive trial, we successfully isolated and correctly identified teleost prey DNA without incurring environmental DNA contamination from the surrounding seawater. In the field, we were able to reconstruct high-resolution teleost dietary information from juvenile C. leucas fDNA that was generally consistent with expectations based on published diet studies of this species. While further investigation is needed to validate the method for larger sharks and other species, it is expected to be broadly applicable to aquatic vertebrates and provides an opportunity to advance our understanding of trophic interactions in marine and freshwater systems.

The context dependence of non-consumptive predator effects.

Non-consumptive predator effects (NCEs) are now widely recognised for their capacity to shape ecosystem structure and function. Yet, forecasting the propagation of these predator-induced trait changes through particular communities remains a challenge. Accordingly, focusing on plasticity in prey anti-predator behaviours, we conceptualise the multi-stage process by which predators trigger direct and indirect NCEs, review and distil potential drivers of contingencies into three key categories (properties of the prey, predator and setting), and then provide a general framework for predicting both the nature and strength of direct NCEs. Our review underscores the myriad factors that can generate NCE contingencies while guiding how research might better anticipate and account for them. Moreover, our synthesis highlights the value of mapping both habitat domains and prey-specific patterns of evasion success ('evasion landscapes') as the basis for predicting how direct NCEs are likely to manifest in any particular community. Looking ahead, we highlight two key knowledge gaps that continue to impede a comprehensive understanding of non-consumptive predator-prey interactions and their ecosystem consequences; namely, insufficient empirical exploration of (1) context-dependent indirect NCEs and (2) the ways in which direct and indirect NCEs are shaped interactively by multiple drivers of context dependence.

Population structure, connectivity, and demographic history of an apex marine predator, the bull shark Carcharhinus leucas.

Knowledge of population structure, connectivity, and effective population size remains limited for many marine apex predators, including the bull shark Carcharhinus leucas. This large-bodied coastal shark is distributed worldwide in warm temperate and tropical waters, and uses estuaries and rivers as nurseries. As an apex predator, the bull shark likely plays a vital ecological role within marine food webs, but is at risk due to inshore habitat degradation and various fishing pressures. We investigated the bull shark's global population structure and demographic history by analyzing the genetic diversity of 370 individuals from 11 different locations using 25 microsatellite loci and three mitochondrial genes (CR, nd4, and cytb). Both types of markers revealed clustering between sharks from the Western Atlantic and those from the Western Pacific and the Western Indian Ocean, with no contemporary gene flow. Microsatellite data suggested low differentiation between the Western Indian Ocean and the Western Pacific, but substantial differentiation was found using mitochondrial DNA. Integrating information from both types of markers and using Bayesian computation with a random forest procedure (ABC-RF), this discordance was found to be due to a complete lack of contemporary gene flow. High genetic connectivity was found both within the Western Indian Ocean and within the Western Pacific. In conclusion, these results suggest important structuring of bull shark populations globally with important gene flow occurring along coastlines, highlighting the need for management and conservation plans on regional scales rather than oceanic basin scale.

Effects of anticoagulants on stable-isotope values (δ13 C and δ15 N) of shark blood components.

The effects of anticoagulant EDTA and sodium heparin (SH) on stable carbon δ13 C and nitrogen δ15 N isotopic values of red blood cells (RBC) and blood plasma in juvenile blacktip reef sharks Carcharhinus melanopterus were analysed. Plasma preserved with anticoagulants was not isotopically distinct from plasma stored in no-additive control tubes but RBC δ15 N values exhibited small enrichments when preserved with EDTA and SH. Results suggest EDTA and SH are viable anticoagulants for stable isotopic analyses of blood fractions but further studies are advised to validate results.

Inter-individual differences in ontogenetic trophic shifts among three marine predators.

Ontogenetic niche shifts are widespread. However, individual differences in size at birth, morphology, sex, and personalities can cause variability in behavior. As such, inherent inter-individual differences within populations may lead to context-dependent changes in behavior with animal body size, which is of concern for understanding population dynamics and optimizing ecological monitoring. Using stable carbon and nitrogen isotope values from concurrently sampled tissues, we quantified the direction and magnitude of intraspecific variation in trophic shifts among three shark species, and how these changed with body size: spurdogs (Squalus spp.) in deep-sea habitats off La Réunion, bull sharks (Carcharhinus leucas) in estuarine habitats of the Florida Everglades, and blacktip reef sharks (Carcharhinus melanopterus) in coral reef ecosystems of Moorea, French Polynesia. Intraspecific variation in trophic shifts was limited among spurdogs, and decreased with body size, while bull sharks exhibited greater individual differences in trophic shifts, but also decreased in variability through ontogeny. In contrast, blacktip reef sharks exhibited increased intraspecific variation in trophic interactions with body size. Variability in trophic interactions and ontogenetic shifts are known to be associated with changes in energetic requirements, but can vary with ecological context. Our results suggest that environmental stability may affect variability within populations, and ecosystems with greater spatial and/or temporal variability in environmental conditions, and those with more diverse food webs may facilitate greater individual differences in trophic interactions, and thus ontogenetic trophic shifts. In light of concerns over environmental disturbance, elucidating the contexts that promote or dampen phenotypic variability is invaluable for predicting population- and community-level responses to environmental changes.

Effect of body length, trophic position and habitat use on mercury concentrations of sharks from contrasted ecosystems in the southwestern Indian Ocean.

The non-essential metal mercury (Hg) can have deleterious effects on health of organisms, and tends to bioaccumulate with age in long-lived organisms and to biomagnify along food webs. Because elasmobranchs are fished for human consumption and their Hg levels are frequently above the maximum Hg concentration recommended for fish consumption, understanding the drivers of Hg concentration is of considerable interest. Total Hg concentrations were analysed in muscle tissues of 14 shark and 2 batoid species (n = 339 individuals) sampled across multiple habitats (coastal, open ocean and bathyal) in the southwestern Indian Ocean. Stable isotope ratios of carbon (δ13C) and nitrogen (δ15N) were analysed to assess whether relative trophic position and foraging habitats affected Hg concentrations. Hg concentrations increased with δ15N and body length, highlighting the mechanisms of bioaccumulation and biomagnification in relation with the trophic position and size of the individuals. Habitats where elasmobranchs were collected also affected their Hg concentrations. Bathyal sharks had high Hg concentrations that were almost similar to those of oceanic species, despite their lower relative trophic position. Higher bioavailability of Hg due to its enhanced methylation in deeper waters was considered as the most likely explanation for this result. These results highlight that multiple factors contribute to mercury accumulation in elasmobranchs.

Habitat use of sympatric prey suggests divergent anti-predator responses to recolonizing gray wolves.

The non-consumptive effects of predators on prey are now widely recognized, but the need remains for studies identifying the factors that determine how particular prey species respond behaviorally when threatened with predation. We took advantage of ongoing gray wolf (Canis lupus) recolonization in eastern Washington, USA, to contrast habitat use of two sympatric prey species-mule (Odocoileus hemionus) and white-tailed (O. virginianus) deer-at sites with and without established wolf packs. Under the hypothesis that the nature and scale of responses by these ungulates to wolf predation risk depend on their divergent flight tactics (i.e., modes of fleeing from an approaching predator), we predicted that (1) mule deer would respond to wolves with coarse-scale spatial shifts to rugged terrain favoring their stotting tactic; (2) white-tailed deer would manage wolf risk with fine-scale shifts toward gentle terrain facilitating their galloping tactic within their current home range. Resource selection functions based on 61 mule deer and 59 white-tailed deer equipped with GPS radio-collars from 2013 to 2016 revealed that habitat use for each species was altered by wolf presence, but in divergent ways that supported our predictions. Our findings add to a growing literature highlighting flight behavior as a viable predictor of prey responses to predation risk across multiple ecosystem types. Consequently, they suggest that predators could initiate multiple indirect non-consumptive effects in the same ecosystem that are transmitted by divergent responses of sympatric prey with different flight tactics.