Reconstructing transoceanic migration patterns of Pacific bluefin tuna using a chemical tracer toolbox

Amino acid δ15N in eye lens laminae reveals life-time ontogenetic trophic shifts of a highly migratory species

Investigating the feeding ecology through the ontogenesis of highly migratory species such as the Pacific Bluefin tuna (PBFT; Thunnus orientalis) is difficult due to its extensive home range and cross-oceanic migration. Here, we show the potential of conducting nitrogen stable isotope (δ15N) analyses in bulk tissue and amino acids (AAs) in consecutive eye lens laminae of PBFT to reconstruct the trophic life history for an individual tuna. The δ15Nbulk profiles between individuals caught in the wild and pen-raised were compared. For all individuals, δ15Nbulk values increased with increasing eye lens diameter or fork length, and exhibited low variation among individual profiles despite tuna being captured in different months. Large δ15Nbulk shifts (6.8‰-8.5‰) were quantified between the first and last deposited laminae for each individual, suggesting major ontogenetic changes in either foraging areas or trophic position. AA δ15N values indicate that this highly migratory schooling predator switches feeding areas from lower to higher δ15N baseline values, reflecting feeding on both sides of the north Pacific, and tends to feed on prey of higher trophic position as it grows. Together, stable isotope analysis in bulk tissue and individual AAs in eye lens laminae could be a powerful approach to investigate changes in the foraging habitat and trophic status of highly migratory species.

Body size and isotopic profiles enable discrimination between long-term resident and highly migrant contingents of Atlantic bluefin tuna

The Atlantic bluefin tuna (ABFT) population ranges throughout the Atlantic Ocean and Mediterranean Sea, and consists of multiple contingents that use diverse habitats and show different movement patterns over the life cycle. Based on body size, elemental and isotopic data of C and N in muscle and liver, we analysed eastern-stock ABFT by comparing mid-sized individuals caught by hook-and-line gears with larger individuals harvested from traps in the Strait of Gibraltar (SoG). Our results show that trophic-related chemical markers have potential for separating temporarily sympatric contingents throughout the ABFT population range, reflecting size-dependent spatial distribution and differential patterns of residency and migration. We present evidence of long-term residency of ABFT in the SoG that persists until the estimated age of ∼5 years (size of ∼140 cm in straight fork length). This age apparently marks a turning point in the life history, where there occurs an ontogenetic switch in the migratory behaviour and distributional pattern. This study contributes new insight into our knowledge about size structure and residency-movement patterns in eastern ABFT. It shows distinct size-dependent migratory and spatial dynamics. The present results encourage further investigation on poorly studied ABFT contingents for a better understanding of the population dynamics towards more comprehensive management plans.

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.

Nursery origin of yellowfin tuna in the western Atlantic Ocean: significance of Caribbean Sea and trans-Atlantic migrants

Natural geochemical markers in the otolith of yellowfin tuna (Thunnus albacares) were used to establish nursery-specific signatures for investigating the origin of fish captured in the western Atlantic Ocean (WAO). Two classes of chemical markers (trace elements, stable isotopes) were used to first establish nursery-specific signatures of age-0 yellowfin tuna from four primary production zones in the Atlantic Ocean: Gulf of Mexico, Caribbean Sea, Cape Verde, and Gulf of Guinea. Next, mixture and individual assignment methods were applied to predict the origin of sub-adult and adult yellowfin tuna from two regions in the WAO (Gulf of Mexico, Mid Atlantic Bight) by relating otolith core signatures (corresponding to age-0 period) to baseline signatures of age-0 fish from each nursery. Significant numbers of migrants from Caribbean Sea and eastern Atlantic Ocean (EAO) production zones (Gulf of Guinea, Cape Verde) were detected in the WAO, suggesting that fisheries in this region were subsidized by outside spawning/nursery areas. Contributions from local production (Gulf of Mexico) were also evident in samples from both WAO fisheries, but highly variable from year to year. High levels of mixing by yellowfin tuna from the different production zones and pronounced interannual trends in nursery-specific contribution rates in the WAO emphasize the complex and dynamic nature of this species' stock structure and population connectivity. Given that geographic shifts in distribution across national or political boundaries leads to governance and management challenges, this study highlights the need for temporally resolved estimates of nursery origin to refine assessment models and promote the sustainable harvest of this species.

Toxicity and Starvation Induce Major Trophic Isotope Variation in Daphnia Individuals: A Diet Switch Experiment Using Eight Phytoplankton Species of Differing Nutritional Quality

Stable isotope values can express resource usage by organisms, but their precise interpretation is predicated using a controlled experiment-based validation process. Here, we develop a stable isotope tracking approach towards exploring resource shifts in a key primary consumer species Daphnia magna. We used a diet switch experiment and model fitting to quantify the stable carbon (δ¹³C) and nitrogen (δ¹⁵N) isotope turnover rates and discrimination factors for eight dietary sources of the plankton species that differ in their cellular organization (unicellular or filamentous), pigment and nutrient compositions (sterols and polyunsaturated fatty acids), and secondary metabolite production rates. We also conduct a starvation experiment. We evaluate nine tissue turnover models using Akaike's information criterion and estimate the repetitive trophic discrimination factors. Using the parameter estimates, we calculate the hourly stable isotope turnover rates. We report an exceedingly faster turnover value following dietary switching (72 to 96 h) and a measurable variation in trophic discrimination factors. The results show that toxic stress and the dietary quantity and quality induce trophic isotope variation in Daphnia individuals. This study provides insight into the physiological processes that underpin stable isotope patterns. We explicitly test multiple alternative dietary sources and fasting and discuss the parameters that are fundamental for field- and laboratory-based stable isotope studies.

Philopatry as a Tool to Define Tentative Closed Migration Cycles and Conservation Areas for Large Pelagic Fishes in the Pacific

Migrations of large pelagic fishes across the Pacific are usually inferred from tagging or genetic studies. Even though these techniques have improved over time, they still fail to demonstrate large transoceanic migrations, usually proposing ‘routes’ that do not cycle seasonally. The current study uses the concept of ‘philopatry’ in 11 large pelagic fish species, i.e., the tendency for animals to return to their natal site to reproduce. Tentative migration routes and maps emerge by applying this concept to the movements extracted through a comprehensive review of the literature on satellite and conventional tagging, and population and subpopulation linkages inferred from genetic and/or genomic studies. Moreover, when comparing these proposed migration routes and the mapped reconstructed catch (1950–2016, Sea Around Us) of each species in the Pacific, similarities emerge, reinforcing the accuracy of these migration cycles informed by philopatry. Finally, by superposing the migration routes of our 11 species, we identified areas of the Pacific that are part of the inferred migration routes of multiple species, leading to a discussion of possible ‘blue corridors’ that would protect the studied species’ key migration routes and stocks, which are important for the fisheries, culture and nutrition of Pacific islanders.

Flexible use of a dynamic energy landscape buffers a marine predator against extreme climate variability

Animal migrations track predictable seasonal patterns of resource availability and suitable thermal habitat. As climate change alters this ‘energy landscape’, some migratory species may struggle to adapt. We examined how climate variability influences movements, thermal habitat selection and energy intake by juvenile Pacific bluefin tuna (Thunnus orientalis) during seasonal foraging migrations in the California Current. We tracked 242 tuna across 15 years (2002–2016) with high-resolution archival tags, estimating their daily energy intake via abdominal warming associated with digestion (the ‘heat increment of feeding’). The poleward extent of foraging migrations was flexible in response to climate variability, allowing tuna to track poleward displacements of thermal habitat where their standard metabolic rates were minimized. During a marine heatwave that saw temperature anomalies of up to +2.5°C in the California Current, spatially explicit energy intake by tuna was approximately 15% lower than average. However, by shifting their mean seasonal migration approximately 900 km poleward, tuna remained in waters within their optimal temperature range and increased their energy intake. Our findings illustrate how tradeoffs between physiology and prey availability structure migration in a highly mobile vertebrate, and suggest that flexible migration strategies can buffer animals against energetic costs associated with climate variability and change.

Integrating multiple chemical tracers to elucidate the diet and habitat of Cookiecutter Sharks

The Cookiecutter shark (Isistius brasiliensis) is an ectoparasitic, mesopelagic shark that is known for removing plugs of tissue from larger prey, including teleosts, chondrichthyans, cephalopods, and marine mammals. Although this species is widely distributed throughout the world’s tropical and subtropical oceanic waters, like many deep-water species, it remains very poorly understood due to its mesopelagic distribution. We used a suite of biochemical tracers, including stable isotope analysis (SIA), fatty acid analysis (FAA), and environmental DNA (eDNA), to investigate the trophic ecology of this species in the Central Pacific around Hawaii. We found that large epipelagic prey constituted a relatively minor part of the overall diet. Surprisingly, small micronektonic and forage species (meso- and epipelagic) are the most important prey group for Cookiecutter sharks across the studied size range (17–43 cm total length), with larger mesopelagic species or species that exhibit diel vertical migration also being important prey. These results were consistent across all the tracer techniques employed. Our results indicate that Cookiecutter sharks play a unique role in pelagic food webs, feeding on prey ranging from the largest apex predators to small, low trophic level species, in particular those that overlap with the depth distribution of the sharks throughout the diel cycle. We also found evidence of a potential shift in diet and/or habitat with size and season. Environmental DNA metabarcoding revealed new prey items for Cookiecutter sharks while also demonstrating that eDNA can be used to identify recent prey in stomachs frozen for extended periods. Integrating across chemical tracers is a powerful tool for investigating the ecology of elusive and difficult to study species, such as meso- and bathypelagic chondrichthyans, and can increase the amount of information gained from small sample sizes. Better resolving the foraging ecology of these mesopelagic predators is critical for effective conservation and management of these taxa and ecosystems, which are intrinsically vulnerable to overfishing and exploitation.

Hydrodynamic Analysis for the Morphing Median Fins of Tuna during Yaw Motions

Tuna can change the area and shape of the median fins, including the first dorsal, second dorsal, and anal fins. The morphing median fins have the ability of adjusting the hydrodynamic forces, thereby affecting the yaw mobility of tuna to a certain extent. In this paper, the hydrodynamic analysis of the median fins under different morphing states is carried out by the numerical method, so as to clarify the influence of the erected median fins on the yaw maneuvers. By comparing the two morphing states of erected and depressed, it can be concluded that the erected median fins can increase their own hydrodynamic forces during the yaw movement. However, the second dorsal and anal fins have limited influence on the yaw maneuverability, and they tend to maintain the stability of tuna. The first dorsal fin has more lift increment in the erection state, which can obviously affect the hydrodynamic performance of tuna. Moreover, as the median fins are erected, the hydrodynamic forces of the tuna's body increase synchronously due to the interaction between the body and the median fins, which is also very beneficial to the yaw motion. This study indicates that tuna can use the morphing median fins to adjust its mobility and stability, which provides a new idea for the design of robotic fish.

The potential of fatty acid isotopes to trace trophic transfer in aquatic food-webs

Compound-specific isotope analyses (CSIA) of fatty acids (FA) constitute a promising tool for tracing energy flows in food-webs. However, past applications of FA-specific carbon isotope analyses have been restricted to a relatively coarse food-source separation and mainly quantified dietary contributions from different habitats. Our aim was to evaluate the potential of FA-CSIA to provide high-resolution data on within-system energy flows using algae and zooplankton as model organisms. First, we investigated the power of FA-CSIA to distinguish among four different algae groups, namely cyanobacteria, chlorophytes, haptophytes and diatoms. We found substantial within-group variation but also demonstrated that δ¹³C of several FA (e.g. 18:3ω3 or 18:4ω3) differed among taxa, resulting in group-specific isotopic fingerprints. Second, we assessed changes in FA isotope ratios with trophic transfer. Isotope fractionation was highly variable in daphnids and rotifers exposed to different food sources. Only δ¹³C of nutritionally valuable poly-unsaturated FA remained relatively constant, highlighting their potential as dietary tracers. The variability in fractionation was partly driven by the identity of food sources. Such systematic effects likely reflect the impact of dietary quality on consumers' metabolism and suggest that FA isotopes could be useful nutritional indicators in the field. Overall, our results reveal that the variability of FA isotope ratios provides a substantial challenge, but that FA-CSIA nevertheless have several promising applications in food-web ecology. This article is part of the theme issue 'The next horizons for lipids as 'trophic biomarkers': evidence and significance of consumer modification of dietary fatty acids'.

Tunas as a high-performance fish platform for inspiring the next generation of autonomous underwater vehicles

Tunas of the genus Thunnus are a group of high-performance pelagic fishes with many locomotor traits that are convergently shared with other high-performance fish groups. Because of their swimming abilities, tunas continue to be an inspiration for both comparative biomechanics and the design of biomimetic autonomous underwater vehicles (AUVs). Despite the strong history of studies in tuna physiology and current interest in tuna biomechanics and bio-inspired design, we lack quantitative data on the function of many features of tunas. Here we present data on the morphology, behavior, and function of tunas, focusing especially on experimentally examining the function of tuna lateral keels, finlets, and pectoral fins by using simple physical models. We find that both triangular lateral keels and flexible finlets decrease power requirements during swimming, likely by reducing lateral forces and yaw torques (compared to models either without keels or with rectangular keels, and models with stiff finlets or strip fins of equal area, respectively). However, both triangular keels and flexible finlets generate less thrust than other models either without these features or with modified keels or finlets, leading to a tradeoff between power consumption and thrust. In addition, we use micro computed tomography (µCT) to show that the flexible lateral keels possess a lateral line canal, suggesting these keels have a sensory function. The curved and fully-attached base of tuna pectoral fins provides high lift-to-drag ratio at low angles of attack, and generates the highest torques across speeds and angles of attack. Therefore, curved, fully-attached pectoral fins grant both better gliding and maneuvering performance compared to flat or curved, partially-attached designs. We provide both 3D models of tuna morphology derived from µCT scans and conclusions about the performance effects of tuna-like features as a resource for future biological and engineering work for next-generation tuna-inspired AUV designs.

Long-term trends in the foraging ecology and habitat use of an endangered species: an isotopic perspective

Evaluating long-term drivers of foraging ecology and population productivity is crucial for providing ecological baselines and forecasting species responses to future environmental conditions. Here, we examine the trophic ecology and habitat use of North Atlantic leatherback turtles (St. Croix nesting population) and investigate the effects of large-scale oceanographic conditions on leatherback foraging dynamics. We used bulk and compound-specific nitrogen isotope analysis of amino acids (CSIA-AA) to estimate leatherback trophic position (TP) over an 18-year period, compare these estimates with TP estimates from a Pacific leatherback population, and elucidate the pre-nesting habitat use patterns of leatherbacks. Our secondary objective was to use oceanographic indices and nesting information from St. Croix leatherbacks to evaluate relationships between trophic ecology, nesting parameters, and regional environmental conditions measured by the North Atlantic Oscillation (NAO) and Atlantic Multidecadal Oscillation. We found no change in leatherback TP over time and no difference in TP between Atlantic and Pacific leatherbacks, indicating that differences in trophic ecology between populations are an unlikely driver of the population dichotomy between Pacific and Atlantic leatherbacks. Isotope data suggested that St. Croix leatherbacks inhabit multiple oceanic regions prior to nesting, although, like their conspecifics in the Pacific, individuals exhibit fidelity to specific foraging regions. Leatherback nesting parameters were weakly related to the NAO, which may suggest that positive NAO phases benefit St. Croix leatherbacks, potentially through increases in resource availability in their foraging areas. Our data contribute to the understanding of leatherback turtle ecology and potential mechanistic drivers of the dichotomy between populations of this protected species.

Assessing seasonal changes in animal diets with stable-isotope analysis of amino acids: a migratory boreal songbird switches diet over its annual cycle

Tools to study seasonal changes in animal diets are needed to address a wide range of ecological questions. This is especially true of migratory animals that experience distinct environments where diets may be substantially different. However, tracking diets of individuals that move vast distances has proven difficult. Compound-specific isotope analysis has emerged as a valuable tool to study diets but has been little used to study dietary changes of migratory animals. Using this technique, we quantify seasonal variation in the annual diet of a migratory songbird (gray-cheeked thrush, Catharus minimus) and test the hypothesis that migrants change their diet in response to the energetic requirements of different periods of the annual cycle. By measuring δ¹³C and δ¹⁵N values of amino acids from feathers grown on the breeding grounds, blood formed during migration and claw grown on the wintering grounds, we found that migration is associated with greater consumption of fruit, compared to the breeding or wintering periods. This was confirmed by the lower trophic position of blood compared to feather and claw, by a decrease in the δ¹⁵N value of the source amino acid phenylalanine in blood as a function of days of stopover, and by the positive correlation between δ¹⁵N and δ¹³C values of phenylalanine in blood, and not in feather or claw. This study illustrates how isotopic analysis of amino acids can contribute to understand food webs, seasonal dietary changes and metabolic routing of nutrients in migratory animals.

Scale diversity in bigeye tuna (Thunnus obesus): Fat-filled trabecular scales made of cellular bone

Tunas of the genus Thunnus possess many morphological and physiological adaptations for their high-performance epipelagic ecology. Although Thunnus anatomy has been studied, there are no quantitative studies on the structure of their scales. We investigated the scales of bigeye tuna (Thunnus obesus) from ten regions of the body using micro computed tomography (µCT)-scanning and histology to quantitatively and qualitatively compare regional scale morphology. We found a diversity of scale sizes and shapes across the body of bigeye tuna and discriminant function analysis on variables derived from µCT-data showed that scales across the body differ quantitatively in shape and size. We also report the discovery of a novel scale type in corselet, tail, and cheek regions. These modified scales are ossified shells supported by internal trabeculae, filled with fat, and possessing an internal blood supply. Histological analysis showed that the outer lamellar layers of these thickened scales are composed of cellular bone, unexpected for a perciform fish in which bone is typically acellular. In the fairing region of the anterior body, these fat-filled scales are stacked in layers up to five scales deep, forming a thickened bony casing. Cheek scales also possess a fat-filled internal trabecular structure, while most posterior body scales are more plate-like and similar to typical teleost scales. While the function of these novel fat-filled scales is unknown, we explore several possible hypotheses for their function.

Assessing Fukushima-Derived Radiocesium in Migratory Pacific Predators

The 2011 release of Fukushima-derived radionuclides into the Pacific Ocean made migratory sharks, teleosts, and marine mammals a source of speculation and anxiety regarding radiocesium (^134+137Cs) contamination, despite a lack of actual radiocesium measurements for these taxa. We measured radiocesium in a diverse suite of large predators from the North Pacific Ocean and report no detectable (i.e., ≥ 0.1 Bq kg^-1 dry wt) Fukushima-derived ¹³⁴Cs in all samples, except in one olive ridley sea turtle (Lepidochelys olivacea) with trace levels (0.1 Bq kg^-1). Levels of ¹³⁷Cs varied within and across taxa, but were generally consistent with pre-Fukushima levels and were lower than naturally occurring ⁴⁰K by one to one to two orders of magnitude. Predator size had a weaker effect on ¹³⁷Cs and ⁴⁰K levels than tissue lipid content. Predator stable isotope values (δ¹³C and δ¹⁵N) were used to infer recent migration patterns, and showed that predators in the central, eastern, and western Pacific should not be assumed to accumulate detectable levels of radiocesium a priori. Nondetection of ¹³⁴Cs and low levels of ¹³⁷Cs in diverse marine megafauna far from Fukushima confirms negligible increases in radiocesium, with levels comparable to those prior to the release from Fukushima. Reported levels can inform recently developed models of cesium transport and bioaccumulation in marine species.

Stable isotope analyses of feather amino acids identify penguin migration strategies at ocean basin scales

Identifying the at-sea distribution of wide-ranging marine predators is critical to understanding their ecology. Advances in electronic tracking devices and intrinsic biogeochemical markers have greatly improved our ability to track animal movements on ocean-wide scales. Here, we show that, in combination with direct tracking, stable carbon isotope analysis of essential amino acids in tail feathers provides the ability to track the movement patterns of two, wide-ranging penguin species over ocean basin scales. In addition, we use this isotopic approach across multiple breeding colonies in the Scotia Arc to evaluate migration trends at a regional scale that would be logistically challenging using direct tracking alone.

Intrapopulation variability in the timing of ontogenetic habitat shifts in sea turtles revealed using δ15 N values from bone growth rings

Determining location and timing of ontogenetic shifts in the habitat use of highly migratory species, along with possible intrapopulation variation in these shifts, is essential for understanding mechanisms driving alternate life histories and assessing overall population trends. Measuring variations in multi-year habitat-use patterns is especially difficult for remote oceanic species. To investigate the potential for differential habitat use among migratory marine vertebrates, we measured the naturally occurring stable nitrogen isotope (δ¹⁵ N) patterns that differentiate distinct ocean regions to create a 'regional isotope characterization', analysed the δ¹⁵ N values from annual bone growth layer rings from dead-stranded animals, and then combined the bone and regional isotope data to track individual animal movement patterns over multiple years. We used humeri from juvenile North Pacific loggerhead turtles (Caretta caretta), animals that undergo long migrations across the North Pacific Ocean (NPO), using multiple discrete regions as they develop to adulthood. Typical of many migratory marine species, ontogenetic changes in habitat use throughout their decades-long juvenile stage is poorly understood, but each potential habitat has unique foraging opportunities and spatially explicit natural and anthropogenic threats that could affect key life-history parameters. We found a bimodal size/age distribution in the timing that juveniles underwent an ontogenetic habitat shift from the oceanic central North Pacific (CNP) to the neritic east Pacific region near the Baja California Peninsula (BCP) (42·7 ± 7·2 vs. 68·3 ± 3·4 cm carapace length, 7·5 ± 2·7 vs. 15·6 ± 1·7 years). Important to the survival of this population, these disparate habitats differ considerably in their food availability, energy requirements and threats, and these differences can influence life-history parameters such as growth, survival and future fecundity. This is the first evidence of alternative ontogenetic shifts and habitat-use patterns for juveniles foraging in the eastern NPO. We combine two techniques, skeletochronology and stable isotope analysis, to reconstruct multi-year habitat-use patterns of a remote migratory species, linked to estimated ages and body sizes of individuals, to reveal variable ontogeny during the juvenile life stage that could drive alternate life histories and that has the potential to illuminate the migration patterns for other species with accretionary tissues.

The Global Food System as a Transport Pathway for Hazardous Chemicals: The Missing Link between Emissions and Exposure

BACKGROUND

Food is a major pathway for human exposure to hazardous chemicals. The modern food system is becoming increasingly complex and globalized, but models for food-borne exposure typically assume locally derived diets or use concentrations directly measured in foods without accounting for food origin. Such approaches may not reflect actual chemical intakes because concentrations depend on food origin, and representative analysis is seldom available. Processing, packaging, storage, and transportation also impart different chemicals to food and are not yet adequately addressed. Thus, the link between environmental emissions and realistic human exposure is effectively broken.

OBJECTIVES

We discuss the need for a fully integrated treatment of the modern industrialized food system, and we propose strategies for using existing models and relevant supporting data sources to track chemicals during production, processing, packaging, storage, and transport.

DISCUSSION

Fate and bioaccumulation models describe how chemicals distribute in the environment and accumulate through local food webs. Human exposure models can use concentrations in food to determine body burdens based on individual or population characteristics. New models now include the impacts of processing and packaging but are far from comprehensive. We propose to close the gap between emissions and exposure by utilizing a wider variety of models and data sources, including global food trade data, processing, and packaging models.

CONCLUSIONS

A comprehensive approach that takes into account the complexity of the modern global food system is essential to enable better prediction of human exposure to chemicals in food, sound risk assessments, and more focused risk abatement strategies. Citation: Ng CA, von Goetz N. 2017. The global food system as a transport pathway for hazardous chemicals: the missing link between emissions and exposure. Environ Health Perspect 125:1-7; http://dx.doi.org/10.1289/EHP168.