Can biochemical tracers reveal ontogenetic trophic shift and individual prey selection in white sharks from Guadalupe Island, Northeast Pacific?

Refining the role of apex predators in marine food webs is a necessary step in predicting the consequences of their global decline under the footprint of fishing activities. White sharks (Carcharodon carcharias) are vulnerable predators, performing large migrations and able to forage on a variety of prey in different habitats. In the Northeast Pacific, juvenile and adult white sharks are found seasonally at the same aggregation sites, such as Guadalupe Island off Mexico. While adults are thought to target local pinniped colonies, very few predator-prey interactions have been documented and the diet of juveniles in this area remains poorly understood. Here we used carbon/nitrogen stable isotopes and fatty acids to characterize the trophic ecology of white sharks at Guadalupe Island. In contrast to the ontogenetic trophic shift paradigm, we detected no influence of size on muscle stable isotope and fatty acid composition, revealing no significant dietary variation between juvenile and adult sharks. Stable isotopes did not allow definitive conclusions to be drawn regarding the diet of white sharks at Guadalupe Island, due to significant variability in the contribution of different potential prey depending on the trophic discrimination factors used. However, most sharks were rich in polyunsaturated fatty acids (such as long-chain omega 3), suggesting a local diet of mainly pelagic prey (potentially large fish or cephalopods). A few individuals appeared to show recent consumption of pinnipeds, with higher proportions of saturated and monounsaturated fatty acids. These individual differences in fatty acid composition could reflect an ecological trade-off between consumption of prey rich in fat (marine mammals) versus prey rich in polyunsaturated fatty acids (pelagic prey), respectively meeting the energetic and physiological needs of white sharks. Although ontogenetic trophic changes were not able to be discerned, our results provide new insights into the physiological drivers of predator-prey interactions, which can benefit the definition of conservation strategies in a changing ocean.

Do marine protected areas influence mercury exposure? Insights from a shark community in the tropical Northeast Pacific.

Biomass depletion caused by overfishing is likely to alter the structure of food webs and impact mercury transfer to marine predators. Although marine protected areas (MPAs) are spared from fishing pressure, their influence on biota mercury levels is poorly understood. Here, we used carbon and nitrogen stable isotope compositions as well as mercury concentrations in fin clips to characterize foraging habitat and mercury exposure of a shark community composed of migratory and resident species of the Revillagigedo archipelago, an offshore MPA in the Northeast Pacific off Mexico. We found that the probability of finding migratory sharks in the isotopic niche of Revillagigedo-resident sharks was low, likely reflecting the use of habitats outside the archipelago by highly mobile species. Community-wide variations in mercury were primarily explained by shark length, revealing that bioaccumulation was the main driver of Hg concentrations. We failed to detect a clear effect of foraging habitat on shark mercury exposure, which may be related to migratory species using both exploited and protected areas when moving outside the Revillagigedo MPA. More similar studies on the potential mitigation of Hg contamination by MPAs are needed in the future if fishing pressure increases to satisfy the growing global human population.

Are tunas relevant bioindicators of mercury concentrations in the global ocean?

Humans are exposed to toxic methylmercury mainly by consuming marine fish. The Minamata Convention aims at reducing anthropogenic mercury releases to protect human and ecosystem health, employing monitoring programs to meet its objectives. Tunas are suspected to be sentinels of mercury exposure in the ocean, though not evidenced yet. Here, we conducted a literature review of mercury concentrations in tropical tunas (bigeye, yellowfin, and skipjack) and albacore, the four most exploited tunas worldwide. Strong spatial patterns of tuna mercury concentrations were shown, mainly explained by fish size, and methylmercury bioavailability in marine food web, suggesting that tunas reflect spatial trends of mercury exposure in their ecosystem. The few mercury long-term trends in tunas were contrasted and sometimes disconnected to estimated regional changes in atmospheric emissions and deposition, highlighting potential confounding effects of legacy mercury, and complex reactions governing the fate of mercury in the ocean. Inter-species differences of tuna mercury concentrations associated with their distinct ecology suggest that tropical tunas and albacore could be used complementarily to assess the vertical and horizontal variability of methylmercury in the ocean. Overall, this review elevates tunas as relevant bioindicators for the Minamata Convention, and calls for large-scale and continuous mercury measurements within the international community. We provide guidelines for tuna sample collection, preparation, analyses and data standardization with recommended transdisciplinary approaches to explore tuna mercury content in parallel with observation abiotic data, and biogeochemical model outputs. Such global and transdisciplinary biomonitoring is essential to explore the complex mechanisms of the marine methylmercury cycle.

Mercury stable isotopes suggest reduced foraging depth in oxygen minimum zones for blue sharks.

Oxygen minimum zones (OMZs) are currently expanding across the global ocean due to climate change, leading to a compression of usable habitat for several marine species. Mercury stable isotope compositions provide a spatially and temporally integrated view of marine predator foraging habitat and its variability with environmental conditions. Here, we analyzed mercury isotopes in blue sharks Prionace glauca from normoxic waters in the northeastern Atlantic and from the world's largest and shallowest OMZ, located in the northeastern Pacific (NEP). Blue sharks from the NEP OMZ area showed higher Δ199Hg values compared to sharks from the northeastern Atlantic, indicating a reduction in foraging depth of approximately 200 m. Our study suggests for the first time that blue shark feeding depth is altered by expanding OMZs and illustrates the use of mercury isotopes to assess the impacts of ocean deoxygenation on the vertical foraging habitat of pelagic predators.

Mercury concentrations in tuna blood and muscle mirror seawater methylmercury in the Western and Central Pacific Ocean.

Understanding the relationship between mercury in seafood and the distribution of oceanic methylmercury is key to understand human mercury exposure. Here, we determined mercury concentrations in muscle and blood of bigeye and yellowfin tunas from the Western and Central Pacific. Results showed similar latitudinal patterns in tuna blood and muscle, indicating that both tissues are good candidates for mercury monitoring. Complementary tuna species analyses indicated species- and tissue- specific mercury patterns, highlighting differences in physiologic processes of mercury uptake and accumulation associated with tuna vertical habitat. Tuna mercury content was correlated to ambient seawater methylmercury concentrations, with blood being enriched at a higher rate than muscle with increasing habitat depth. The consideration of a significant uptake of dissolved methylmercury from seawater in tuna, in addition to assimilation from food, might be interesting to test in models to represent the spatiotemporal evolutions of mercury in tuna under different mercury emission scenarios.

Mercury concentrations, biomagnification and isotopic discrimination factors in two seabird species from the Humboldt Current ecosystem.

Assessing mercury (Hg) biomagnification requires the description of prey-predator relationships, for each species and ecosystem, usually based on carbon and nitrogen isotope analyses. Here, we analyzed two seabirds from the Humboldt Current ecosystem, the Guanay cormorant (Phalacrocorax bougainvillii) and the Peruvian booby (Sula variegata), as well as their main prey, the Peruvian anchovy (Engraulis ringens). We reported Hg concentrations, Hg biomagnification (BMF) and isotopic discrimination factors (Δ13C and Δ15N) in seabird whole blood. BMFs and Δ13C in our study (on wild birds where diet was not controlled) were similar to other piscivorous seabirds previously studied in captive settings, but Δ15N were lower than most captive experiments. We observed lower Hg concentrations in Humboldt seabirds compared to other oligotrophic ecosystems, possibly due to Hg biodilution in the high biomass of the first trophic levels. This work calls for a better characterization of Hg trophic dynamics in productive upwelling ecosystems.

Evidence that Pacific tuna mercury levels are driven by marine methylmercury production and anthropogenic inputs.

Pacific Ocean tuna is among the most-consumed seafood products but contains relatively high levels of the neurotoxin methylmercury. Limited observations suggest tuna mercury levels vary in space and time, yet the drivers are not well understood. Here, we map mercury concentrations in skipjack tuna across the Pacific Ocean and build generalized additive models to quantify the anthropogenic, ecological, and biogeochemical drivers. Skipjack mercury levels display a fivefold spatial gradient, with maximum concentrations in the northwest near Asia, intermediate values in the east, and the lowest levels in the west, southwest, and central Pacific. Large spatial differences can be explained by the depth of the seawater methylmercury peak near low-oxygen zones, leading to enhanced tuna mercury concentrations in regions where oxygen depletion is shallow. Despite this natural biogeochemical control, the mercury hotspot in tuna caught near Asia is explained by elevated atmospheric mercury concentrations and/or mercury river inputs to the coastal shelf. While we cannot ignore the legacy mercury contribution from other regions to the Pacific Ocean (e.g., North America and Europe), our results suggest that recent anthropogenic mercury release, which is currently largest in Asia, contributes directly to present-day human mercury exposure.

Foraging plasticity diversifies mercury exposure sources and bioaccumulation patterns in the world's largest predatory fish.

Large marine predators exhibit high concentrations of mercury (Hg) as neurotoxic methylmercury, and the potential impacts of global change on Hg contamination in these species remain highly debated. Current contaminant model predictions do not account for intraspecific variability in Hg exposure and may fail to reflect the diversity of future Hg levels among conspecific populations or individuals, especially for top predators displaying a wide range of ecological traits. Here, we used Hg isotopic compositions to show that Hg exposure sources varied significantly between and within three populations of white sharks (Carcharodon carcharias) with contrasting ecology: the north-eastern Pacific, eastern Australasian, and south-western Australasian populations. Through Δ200Hg signatures in shark tissues, we found that atmospheric Hg deposition pathways to the marine environment differed between coastal and offshore habitats. Discrepancies in δ202Hg and Δ199Hg signatures among white sharks provided evidence for intraspecific exposure to distinct sources of marine methylmercury, attributed to population and ontogenetic shifts in foraging habitat and prey composition. We finally observed a strong divergence in Hg accumulation rates between populations, leading to three times higher Hg concentrations in large Australasian sharks compared to north-eastern Pacific sharks, and likely due to different trophic strategies adopted by adult sharks across populations. This study illustrates the variety of Hg exposure sources and bioaccumulation patterns that can be found within a single species and suggests that intraspecific variability needs to be considered when assessing future trajectories of Hg levels in marine predators.

ENSO Climate Forcing of the Marine Mercury Cycle in the Peruvian Upwelling Zone Does Not Affect Methylmercury Levels of Marine Avian Top Predators.

Climate change is expected to affect marine mercury (Hg) biogeochemistry and biomagnification. Recent modeling work suggested that ocean warming increases methylmercury (MeHg) levels in fish. Here, we studied the influence of El Niño Southern Oscillations (ENSO) on Hg concentrations and stable isotopes in time series of seabird blood from the Peruvian upwelling and oxygen minimum zone. Between 2009 and 2016, La Niña (2011) and El Niño conditions (2015-2016) were accompanied by sea surface temperature anomalies up to 3 °C, oxycline depth change (20-100 m), and strong primary production gradients. Seabird Hg levels were stable and did not co-vary significantly with oceanographic parameters, nor with anchovy biomass, the primary dietary source to seabirds (90%). In contrast, seabird Δ199Hg, proxy for marine photochemical MeHg breakdown, and δ15N showed strong interannual variability (up to 0.8 and 3‰, respectively) and sharply decreased during El Niño. We suggest that lower Δ199Hg during El Niño represents reduced MeHg photodegradation due to the deepening of the oxycline. This process was balanced by equally reduced Hg methylation due to reduced productivity, carbon export, and remineralization. The non-dependence of seabird MeHg levels on strong ENSO variability suggests that marine predator MeHg levels may not be as sensitive to climate change as is currently thought.

Comment on Trophic strategy and bleaching resistance in reef-building corals.

In an era of major environmental changes, understanding corals' resistance to bleaching is as crucial as it is challenging. A promising framework for inferring corals' trophic strategies from Stable Isotope Bayesian Ellipses has been recently proposed to this end. As a contribution to this framework, we quantify a risk of bias inherent in its application and propose three alternative adjustments.

Lipid-free tuna muscle samples are suitable for total mercury analysis.

Tropical tunas are largely consumed worldwide, providing major nutritional benefits to humans, but also representing the main exposure to methylmercury, a potent neurotoxin that biomagnifies along food webs. The combination of ecological tracers (nitrogen and carbon stable isotopes, δ15N and δ13C) to mercury concentrations in tunas is scarce yet crucial to better characterize the influence of tuna foraging ecology on mercury exposure and bioaccumulation. Given the difficulties to get modern and historical tuna samples, analyses have to be done on available and unique samples. However, δ13C values are often analysed on lipid-free samples to avoid bias related to lipid content. While lipid extraction with non-polar solvents is known to have no effect on δ15N values, its impact on mercury concentrations is still unclear. We used white muscle tissues of three tropical tuna species to evaluate the efficiency and repeatability of different lipid extraction protocols commonly used in δ13C and δ15N analysis. Dichloromethane was more efficient than cyclohexane in extracting lipids in tuna muscle, while the automated method appeared more efficient but as repeatable as the manual method. Lipid extraction with dichloromethane had no effect on mercury concentrations. This may result from i) the affinity of methylmercury to proteins in tuna flesh, ii) the low lipid content in tropical tuna muscle samples, and iii) the non-polar nature of dichloromethane. Our study suggests that lipid-free samples, usually prepared for tropical tuna foraging ecology research, can be used equivalently to bulk samples to document in parallel mercury concentrations at a global scale.

Stable mercury concentrations of tropical tuna in the south western Pacific ocean: An 18-year monitoring study.

Global anthropogenic mercury (Hg) emissions to the atmosphere since industrialization are widely considered to be responsible for a significant increase in surface ocean Hg concentrations. Still unclear is how those inputs are converted into toxic methylmercury (MeHg) then transferred and biomagnified in oceanic food webs. We used a unique long-term and continuous dataset to explore the temporal Hg trend and variability of three tropical tuna species (yellowfin, bigeye, and skipjack) from the southwestern Pacific Ocean between 2001 and 2018 (n = 590). Temporal trends of muscle nitrogen (δ15N) and carbon (δ13C) stable isotope ratios, amino acid (AA) δ15N values and oceanographic variables were also investigated to examine the potential influence of trophic, biogeochemical and physical processes on the temporal variability of tuna Hg concentrations. For the three species, we detected significant inter-annual variability but no significant long-term trend for Hg concentrations. Inter-annual variability was related to the variability in tuna sampled lengths among years and to tuna muscle δ15N and δ13C values. Complementary AA- and model-estimated phytoplankton δ15N values suggested the influence of baseline processes with enhanced tuna Hg concentrations observed when dinitrogen fixers prevail, possibly fuelling baseline Hg methylation and/or MeHg bioavailability at the base of the food web. Our results show that MeHg trends in top predators do not necessary capture the increasing Hg concentrations in surface waters suspected at the global oceanic scale due to the complex and variable processes governing Hg deposition, methylation, bioavailability and biomagnification. This illustrates the need for long-term standardized monitoring programs of marine biota worldwide.

Global data set for nitrogen and carbon stable isotopes of tunas.

Nitrogen and carbon stable isotope data sets are commonly used to assess complex population to ecosystem responses to natural or anthropogenic changes at regional to global spatial scales, and monthly to decadal timescales. Measured in the tissues of consumers, nitrogen isotopes (δ15 N) are primarily used to estimate trophic position while carbon isotopes (δ13 C) describe habitat associations and feeding pathways. Models of both δ15 N and δ13 C values and their associated variance can be used to estimate likely dietary contributions and niche width and provide inferences about consumer movement and migration. Stable isotope data have added utility when used in combination with other empirical data sets (e.g., stomach content, movement tracking, bioregionalization, contaminant, or fisheries data) and are increasingly relied upon in food web and ecosystem models. While numerous regional studies publish tables of mean δ15 N and δ13 C values, limited individual records have been made available for wider use. Such a deficiency has impeded full utility of the data, which otherwise would facilitate identification of macroscale patterns. The data provided here consist of 4,498 records of individuals of three tuna species, Thunnus alalunga, T. obesus, and T. albacares sampled from all major ocean basins from 2000 to 2015. For each individual tuna, we provide a record of the following: species name, sampling date, sampling location, tuna length, muscle bulk and baseline corrected δ15 N values, and muscle bulk and, where available, lipid corrected δ13 C values. We provide these individual records to support comparative studies and more robust modeling projects seeking to improve understanding of complex marine ecosystem dynamics and their responses to a changing environment. There are no copyright restrictions for research and/or teaching purposes. Users are requested to acknowledge their use of the data in publications, research proposals, websites, and other outlets following the citation instructions in Class III, Section B.

The Twilight Zone as a Major Foraging Habitat and Mercury Source for the Great White Shark.

The twilight zone contains the largest biomass of the world's ocean. Identifying its role in the trophic supply and contaminant exposure of marine megafauna constitutes a critical challenge in the context of global change. The white shark (Carcharodon carcharias) is a threatened species with some of the highest concentrations of neurotoxin methylmercury (MeHg) among marine top predators. Large white sharks migrate seasonally from coastal habitats, where they primarily forage on pinnipeds, to oceanic offshore habitats. Tagging studies suggest that while offshore, white sharks may forage at depth on mesopelagic species, yet no biochemical evidence exists. Here, we used mercury isotopic composition to assess the dietary origin of MeHg contamination in white sharks from the Northeast Pacific Ocean. We estimated that a minimum of 72% of the MeHg accumulated by white sharks originates from the consumption of mesopelagic prey, while a maximum of 25% derives from pinnipeds. In addition to highlighting the potential of mercury isotopes to decipher the complex ecological cycle of marine predators, our study provides evidence that the twilight zone constitutes a crucial foraging habitat for these large predators, which had been suspected for over a decade. Climate change is predicted to expand the production of mesopelagic MeHg and modify the mesopelagic biomass globally. Considering the pivotal role of the twilight zone is therefore essential to better predict both MeHg exposure and trophic supply to white sharks, and effectively protect these key vulnerable predators.

Behavioral and trophic segregations help the Tahiti petrel to cope with the abundance of wedge-tailed shearwater when foraging in oligotrophic tropical waters.

Two species breeding in sympatry are more likely to coexist if their ecological niches are segregated either in time, space or in trophic habits. Here, we combined GPS-tracking, stable isotope analysis and DNA metabarcoding analysis to understand how the rare Tahiti petrel Pseudobulweria rostrata (TP) copes with the very abundant (i.e. 500,000 breeding pairs) wedge-tailed shearwater Ardenna pacifica (WTS) when breeding in sympatry in a tropical area. WTS foraged in restricted areas along their path, while TP predominantly foraged using extensive search behavior, suggesting a more opportunistic foraging strategy. Interspecific overlap of foraging areas was higher than intraspecific overlap. Breeding seasons largely overlap between species during the study, but TP seems to be asynchronous breeders. TP fed upon prey with higher δ15N values than WTS, and their diet was mainly composed of deep-sea organisms. TP could feed upon dead prey floating at the surface while WTS preyed mainly upon fish species that generally move in schools. Our study highlights several segregating mechanisms (temporal, behavioral and trophic) that could facilitate the coexistence of the two species despite the predominant number of WTS, and provides the very first information on the foraging and trophic ecology of the poorly-known TP.

Mercury isotopes as tracers of ecology and metabolism in two sympatric shark species.

In coastal ecosystems, top predators are exposed to a wide variety of nutrient and contaminant sources due to the diversity of trophic webs within inshore marine habitats. Mercury contamination could represent an additional threat to shark populations that are declining worldwide. Here we measured total mercury, carbon and nitrogen isotopes, as well as mercury isotopes, in two co-occurring shark species (the bull shark Carcharhinus leucas and the tiger shark Galeocerdo cuvier) and their potential prey from a coastal ecosystem of the western Indian Ocean (La Réunion Island). Our primary goals were to (i) determine the main trophic Hg sources for sharks and (ii) better characterize their diet composition and foraging habitat. Hg isotope signatures (Δ199Hg and δ202Hg) of shark prey suggested that bull sharks were exposed to methylmercury (MeHg) produced in offshore epipelagic waters, while tiger sharks were exposed to offshore mesopelagic MeHg with additional microbial transformation in slope sediments. Δ199Hg values efficiently traced the ecology of the two predators, demonstrating that bull sharks targeted coastal prey in shallow waters while tiger sharks were mainly foraging on mesopelagic species in the deeper waters of the island slope. Unexpectedly, we found a positive shift in δ202Hg (>1‰) between sharks and their prey, leading to high δ202Hg values in the two shark species (e.g. 1.91 ± 0.52‰ in bull sharks). This large shift in δ202Hg indicates that sharks may display strong MeHg demethylation abilities, possibly reflecting evolutionary pathways for mitigating their MeHg contamination.

Trophic resources and mercury exposure of two silvertip shark populations in the Northeast Pacific Ocean.

Worldwide shark populations have experienced rapid declines over the last decades, mainly due to overfishing. Marine protected areas (MPAs) have thus become an indispensable tool for the protection of these marine predators. Two recently-created MPAs in the Northeast Pacific Ocean, the Revillagigedo National Park and Clipperton Atoll, are characterized by different trophic structures potentially influencing the trophic niche and contaminant exposure of resident sharks in these two sites. In this context, we used carbon (δ13C) and nitrogen (δ15N) stable isotope analyzes as well as total mercury concentrations ([THg]) to assess the effect of foraging site on the trophic niche and Hg levels of juvenile silvertip (ST) sharks Carcharhinus albimarginatus. Analyzing fin clip samples from Revillagigedo and Clipperton, we found that shark δ15N varied spatially in relation to δ15N baselines, suggesting similar trophic position in both MPAs. Moreover, δ13C values indicated that ST sharks from Revillagigedo would feed on different food webs (i.e. both benthic and pelagic) while individuals from Clipperton would only rely on benthic food webs. These differences between MPAs led to a weak overlap of isotopic niches between the two populations, highlighting the site residency of juvenile ST sharks. Within each population, [THg] was not correlated with trophic tracers (δ15N and δ13C) and was also similar between populations. This study revealed no influence of site or food web in [THg] and raises the question of the origin of Hg exposure for reef shark populations in the Northeast Pacific Ocean.

Trends in tuna carbon isotopes suggest global changes in pelagic phytoplankton communities.

Considerable uncertainty remains over how increasing atmospheric CO2 and anthropogenic climate changes are affecting open-ocean marine ecosystems from phytoplankton to top predators. Biological time series data are thus urgently needed for the world's oceans. Here, we use the carbon stable isotope composition of tuna to provide a first insight into the existence of global trends in complex ecosystem dynamics and changes in the oceanic carbon cycle. From 2000 to 2015, considerable declines in δ13 C values of 0.8‰-2.5‰ were observed across three tuna species sampled globally, with more substantial changes in the Pacific Ocean compared to the Atlantic and Indian Oceans. Tuna recorded not only the Suess effect, that is, fossil fuel-derived and isotopically light carbon being incorporated into marine ecosystems, but also recorded profound changes at the base of marine food webs. We suggest a global shift in phytoplankton community structure, for example, a reduction in 13 C-rich phytoplankton such as diatoms, and/or a change in phytoplankton physiology during this period, although this does not rule out other concomitant changes at higher levels in the food webs. Our study establishes tuna δ13 C values as a candidate essential ocean variable to assess complex ecosystem responses to climate change at regional to global scales and over decadal timescales. Finally, this time series will be invaluable in calibrating and validating global earth system models to project changes in marine biota.

Bleaching forces coral's heterotrophy on diazotrophs and Synechococcus.

Coral reefs are threatened by global warming, which disrupts the symbiosis between corals and their photosynthetic symbionts (Symbiodiniaceae), leading to mass coral bleaching. Planktonic diazotrophs or dinitrogen (N2)-fixing prokaryotes are abundant in coral lagoon waters and could be an alternative nutrient source for corals. Here we incubated untreated and bleached coral colonies of Stylophora pistillata with a 15N2-pre-labelled natural plankton assemblage containing diazotrophs. 15N2 assimilation rates in Symbiodiniaceae cells and tissues of bleached corals were 5- and 30-fold higher, respectively, than those measured in untreated corals, demonstrating that corals incorporate more nitrogen derived from planktonic diazotrophs under bleaching conditions. Bleached corals also preferentially fed on Synechococcus, nitrogen-rich picophytoplanktonic cells, instead of Prochlorococcus and picoeukaryotes, which have a lower cellular nitrogen content. By providing an alternative source of bioavailable nitrogen, both the incorporation of nitrogen derived from planktonic diazotrophs and the ingestion of Synechococcus may have profound consequences for coral bleaching recovery, especially for the many coral reef ecosystems characterized by high abundance and activity of planktonic diazotrophs.

A Model of Mercury Distribution in Tuna from the Western and Central Pacific Ocean: Influence of Physiology, Ecology and Environmental Factors.

Information on ocean scale drivers of methylmercury levels and variability in tuna is scarce, yet crucial in the context of anthropogenic mercury (Hg) inputs and potential threats to human health. Here we assess Hg concentrations in three commercial tuna species (bigeye, yellowfin, and albacore, n = 1000) from the Western and Central Pacific Ocean (WCPO). Models were developed to map regional Hg variance and understand the main drivers. Mercury concentrations are enriched in southern latitudes (10°S-20°S) relative to the equator (0°-10°S) for each species, with bigeye exhibiting the strongest spatial gradients. Fish size is the primary factor explaining Hg variance but physical oceanography also contributes, with higher Hg concentrations in regions exhibiting deeper thermoclines. Tuna trophic position and oceanic primary productivity were of weaker importance. Predictive models perform well in the Central Equatorial Pacific and Hawaii, but underestimate Hg concentrations in the Eastern Pacific. A literature review from the global ocean indicates that size tends to govern tuna Hg concentrations, however regional information on vertical habitats, methylmercury production, and/or Hg inputs are needed to understand Hg distribution at a broader scale. Finally, this study establishes a geographical context of Hg levels to weigh the risks and benefits of tuna consumption in the WCPO.