Aquatic biodiversity enhances multiple nutritional benefits to humans
See allHide authors and affiliations
Edited by James J. Elser, University of Montana, Polson, MT, and approved January 8, 2021 (received for review October 13, 2019)

Significance
Food security is not simply about maintaining yields, but it is also about the need for a stable supply of nutritionally diverse foods. Obtaining nutritious food is a major challenge facing humanity, and diverse aquatic ecosystems can help meet this goal. To test how aquatic biodiversity affects human health, we assembled a dataset of nutrients, contaminants, and ecological traits of 801 aquatic species. We used ecological models to quantify the role of species richness and ecological functional diversity and found that these biodiversity dimensions enhanced seafood micronutrient and fatty acid provisioning by the same mechanisms that link biodiversity to productivity in grasslands, forests, and other systems. Our results underscore the need to minimize aquatic biodiversity loss to sustain and improve human well-being.
Abstract
Humanity depends on biodiversity for health, well-being, and a stable environment. As biodiversity change accelerates, we are still discovering the full range of consequences for human health and well-being. Here, we test the hypothesis—derived from biodiversity–ecosystem functioning theory—that species richness and ecological functional diversity allow seafood diets to fulfill multiple nutritional requirements, a condition necessary for human health. We analyzed a newly synthesized dataset of 7,245 observations of nutrient and contaminant concentrations in 801 aquatic animal taxa and found that species with different ecological traits have distinct and complementary micronutrient profiles but little difference in protein content. The same complementarity mechanisms that generate positive biodiversity effects on ecosystem functioning in terrestrial ecosystems also operate in seafood assemblages, allowing more diverse diets to yield increased nutritional benefits independent of total biomass consumed. Notably, nutritional metrics that capture multiple micronutrients and fatty acids essential for human well-being depend more strongly on biodiversity than common ecological measures of function such as productivity, typically reported for grasslands and forests. Furthermore, we found that increasing species richness did not increase the amount of protein in seafood diets and also increased concentrations of toxic metal contaminants in the diet. Seafood-derived micronutrients and fatty acids are important for human health and are a pillar of global food and nutrition security. By drawing upon biodiversity–ecosystem functioning theory, we demonstrate that ecological concepts of biodiversity can deepen our understanding of nature’s benefits to people and unite sustainability goals for biodiversity and human well-being.
Species losses and range shifts because of climate change, harvesting, and other human activities are altering aquatic biodiversity locally and globally (1⇓⇓⇓–5). In aquatic ecosystems, not only are some species severely depleted because of overfishing or habitat loss (3, 6⇓–8), the ecosystem-level dimensions of biodiversity such as the total number of species and their functional diversity have also changed (9). Beyond the loss of particular species, changes in ecosystem-level dimensions of biodiversity threaten numerous ecosystem services to humans, which include the cultural, economic, or health benefits people derive from nature (10⇓⇓–13). In many regions, such as tropical coastal systems, the cumulative impacts of human activities are severe and associated with strong declines in taxonomic and ecological functional diversity (6) and coincide with regions with a high dependence of people upon wild-caught seafood for food and nutrition (14). In temperate regions, where some coastal communities depend on local wild seafood harvests to meet their nutritional needs (15, 16), species richness may be increasing as species recover from exploitation and warmer oceans allow species to expand their ranges into new territory (1, 2, 17).
There is growing concern that biodiversity change leads to changes in human health and well-being (10, 13, 18). Specific and quantitative links between aquatic biodiversity and human health that distinguish contributions of species diversity from those of biomass, as predicted by biodiversity–ecosystem functioning theory, have not been established. At a time of unprecedented global change and increasing reliance on seafood to meet nutritional demands (19), there is an urgent need to understand how changing aquatic ecosystem structure may alter the provisioning of seafood-derived human nutrition.
Seafood, consisting of wild-caught marine and freshwater finfish and invertebrates, provides an important source of protein and calories to humans. Additionally, unlike staple foods such as rice or other grains, seafood can address multiple dimensions of food and nutritional security simultaneously by providing essential micronutrients, such as vitamins, minerals, and polyunsaturated essential fatty acids critical to human health (19⇓⇓–22). Given the multiple attributes of seafood that are valuable to human health, it is possible that the diversity of an aquatic assemblage, distinct from the inclusion of any particularly nutritious species, could support human well-being consistent with a large body of evidence for biodiversity’s major contributions to ecological functions (11, 23⇓⇓–26). Dietary diversity is a basic tenet of a nutritious diet (27) and it is widely appreciated that diets composed of more food groups and more species are more nutritious (28⇓⇓–31). Ecological measures of dietary diversity (diet diversity, species richness, functional diversity, and Simpson’s index of evenness) have been associated with the nutritional value of diets in a range of contexts (27, 29, 32⇓⇓⇓⇓⇓–38). These studies rely on relationships between species included in the diet (or other food intake measures) and nutritional adequacy of reported diets. However, a simple correlation between dietary diversity and a measure of dietary benefits provides only partial support for a claim that biodiversity benefits human well-being, consistent with the same ecological processes by which biodiversity supports numerous ecosystem functions and services (23, 26). We build upon this foundation of empirical relationships between diet diversity and diet quality by placing this question in the quantitative ecological theoretical framework that relates biodiversity to function (24, 25), thereby laying the groundwork for additional development of links between biodiversity science and our understanding of human well-being.
Ecological theory predicts that biodiversity can be ecologically and economically important, apart from the importance of total biomass or the presence of particular species (23, 39). According to theory and over 500 explicit experimental tests (23, 40, 41), diversity in ecological communities and agricultural systems enhances ecosystem functioning by two mechanisms: 1) more diverse assemblages may outperform less diverse assemblages of the same density or biomass of individuals because more diverse assemblages will include more of the possible species and are therefore more likely to include high-performing species, assuming random processes of including species from the species pool (a selection effect), or 2) more diverse assemblages of a given density (or biomass) contain species with complementary functional traits, allowing them to function more efficiently (a complementarity effect) (25, 39). For aquatic animals, increased diversity enhances productivity of fish biomass (42) and also enhances temporal stability of biomass production and total yields (43, 44), providing economic and nutritional benefits to humans related to increased stability of harvests and production of biomass for consumption (43). However, when considering aquatic species from the perspective of human nutrition, functions other than biomass production become relevant because total seafood biomass consumption is not predictive of micronutrient benefits from seafood (45, 46).
Here, we test a hypothesis central to ecological theory in the 21st century: whether biodiversity per se (species richness and ecological functional diversity), distinct from the identities and abundance of species, enhances human well-being (Fig. 1). We chose a measure of human well-being distinct from provision of protein, calories, or total yields—the micronutrient and essential fatty acid benefits of seafood. For increasing biodiversity per se (as opposed to increasing total seafood consumption) to enhance nutritional benefits as predicted by biodiversity–ecosystem functioning theory (25, 47), the amounts of various nutrients within edible tissues must differ among species, and furthermore, nutrient concentrations must trade off among species, such that species that have relatively high concentrations of some nutrients also have relatively low concentrations of others (25). Specifically, a “biodiversity effect” (sensu ref. 25) on nutritional benefits requires that concentrations of multiple nutrients are negatively correlated with each other, or uncorrelated, when compared among species, creating a complementary distribution of nutrients across species. In contrast, if nutrient concentrations in edible tissue are positively correlated for multiple nutrients across species such that, for example, a species containing high amounts of iron also has a high essential fatty acid concentration, thereby containing multiple nutrients in high concentrations simultaneously, seafood species or ecological functional diversity in the diet would not be important. In the case of positive correlations among nutrient concentrations, the ecosystem service of nutritional benefits would be enhanced by consuming more fish biomass or by selecting a few highly nutritious species, without considering species richness or ecological functional diversity.
Aquatic biodiversity increases human well-being because edible species have distinct and complementary multinutrient profiles (A) and differ in mean micro- and macronutrient content (shown here relative to 10 and 25% thresholds of recommended dietary allowance, RDA, guidelines) for representative finfish (Abramis brama, Mullus surmuletus), mollusc (Mytilus galloprovincialis), and crustacean species (Nephrops norvegicus). Biodiversity–ecosystem functioning theory predicts that nutritional benefits, including the number of nutrient RDA targets met per 100 g portion (NT; i, iii) and minimum portion size (Pmin; ii, iv) (B and E), are enhanced with increasing seafood species richness. Orange dots in B and E correspond to potential diets of high and low biodiversity levels. Seafood consumers with limited access to seafood each day may not reach RDA targets if diets are low in diversity (D–F versus A–C; gray shading indicates proportion of population that meets nutrient requirements). DHA: docosahexaenoic acid, EPA: eicosapentaenoic acid.
We aimed to bridge two distinct theoretical frameworks—the biodiversity–ecosystem functioning theory and human nutrition science—by quantitatively testing for effects of aquatic species richness and ecological functional diversity (48, 49) in seafood diets on nutritional benefits via complementarity or selection effects. We used the public health measure of recommended dietary allowance (RDA) index to quantify nutritional benefits. RDAs are nutrient-based reference values that indicate the average daily dietary intake level that is sufficient to meet the nutrient requirement of nearly all (97 to 98%) healthy individuals in a particular life stage and gender group (50). Here, we used the RDA for females aged 19 to 50 y (SI Appendix, Tables S1 and S2; see SI Appendix, Table S1 for definitions of key terms). We measured nutritional value in terms of concentrations relative to RDAs, and we refer to these recommended amounts (or portions thereof) as “RDA targets” (SI Appendix, Tables S1 and S2 and Metrics). We quantified nutritional value in two ways: 1) the minimum amount of seafood tissue (in grams) required to meet given RDA targets (either for a single nutrient or the five micronutrients and fatty acids simultaneously; referred to as “minimum portion size required,” Pmin [SI Appendix, Table S1, Eq. A1, and Metrics]) and 2) the number of nutrients that meet an RDA target in a single 100 g seafood portion (NT, SI Appendix, Table S1, Eq. A2). By considering nutritional value per unit biomass in both metrics, we avoided confounding diversity of seafood consumed with the total amount consumed (Metrics). We first tested two hypotheses: 1) seafood species richness increases NT because of complementarity in nutrient concentrations among species, and 2) seafood species richness increases the nutritional value of a 100 g edible portion of seafood, thereby lowering the minimum portion size, Pmin, and improving the efficiency with which seafood consumers reach nutritional targets (Fig. 1). Following biodiversity–ecosystem functioning theory, we predicted that increased species richness is correlated with ecological functional diversity (51) in potential seafood diets and that ecological functional diversity is related to diversity in the concentration of essential elements and fatty acids that have nutritional value to human consumers, such that species and ecological functional diversity yields increased nutritional benefits. We also tested the hypothesis that seafood diversity increases total intake of heavy metal contaminants because some aquatic animals are known to bioaccumulate toxic metals in their tissues. For this reason, variation in bioaccumulation among species could lead to a biodiversity effect on contaminant intake that is detrimental to human health.
In a global analysis of over 5,040 observations of nutrient concentrations in 547 aquatic species (SI Appendix, Fig. S1), we considered the provision of nutritional benefits to human consumers. To assess whether the relationships between biodiversity and human nutrition benefits depend on the geographic extent (global or local) over which seafood are harvested or accessed (11), we tested whether seafood species richness is associated with higher nutritional value at local scales (versus global scale) in traditional Indigenous seafood diets in North America (SI Appendix, Methods 1.4). Seafood is critical for Indigenous groups, who on average consume seafood at a rate that is 15 times higher than the global average per capita consumption rate (16). To test our hypotheses at the geographic scale of local consumer communities, we complemented our global analysis with additional analyses of 25 to 57 species in 14 geographically constrained groups of species consumed together as part of traditional Indigenous diets (SI Appendix, Methods 1.4).
Results
Diversity of Seafood Nutrient Concentrations.
Biodiversity effects via complementarity or selection require that species differ in their functional traits. We found that the global species pool was highly diverse with regard to concentrations of the micronutrients iron, zinc, calcium, and two fatty acids, docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA), in edible fish tissue relative to RDAs for those micronutrients and fatty acids (Fig. 2; micronutrient and fatty acid geometric coefficients of variation [CV]: ln[iron] = 3.28, ln[calcium] = 3.56, ln[EPA] = 2.62, ln[zinc] = 3.02, and ln[DHA] = 2.17; note the log scale). We observed limited variation in protein concentrations (ln[protein] CV = 0.04). The frequency distribution of trait values such as nutrient concentrations across species may indicate the potential strength of biodiversity effects, with lognormal distributions (such as we observed for micronutrients and fatty acids) more likely to confer strong effects of biodiversity than normal distributions with low dispersion (as we observed for protein). Most species did not meet a single micronutrient or fatty acid RDA in a 100 g portion: fewer than half of the 547 species we examined reached an RDA target of 10% RDA for calcium, iron, and the essential fatty acid EPA in a standard 100 g edible portion of a single species (SI Appendix, Table S3).
Variation in nutrient concentrations differs among taxonomic groups. (Upper) Frequency of reported protein, fat, micronutrient and fatty acid content in 100 g of the edible portion of 547 seafood species. Note the x-axis is plotted on a log scale. (Lower) Proportion of species, and number shown on each bar, with available data that reach 10% of RDA targets for any one, two, or up to five of the micronutrients and fatty acids examined here.
Biodiversity Increased the Nutritional Content of an Edible Portion of Seafood.
We found that seafood species richness not only enhanced nutritional value for consumers selecting seafood from our global species dataset but that seafood species richness per se was essential to meeting nutritional targets for seafood diets with limited biomass consumption. Increasing seafood species richness allowed simulated diets to reach more RDA targets per 100 g of tissue so that nutritional value increased with species richness even as total biomass consumption (e.g., total seafood portion size) remained constant. We quantified the minimum amount of seafood, in grams, that would be required to reach an RDA target (SI Appendix, Table S2) at each of 10 levels of species richness (referred to as minimum portion size, Pmin, for which lower values signify higher nutrition benefits to consumers per gram seafood consumed, Metrics). We then estimated the biodiversity effect using Eq. 3 (Statistical Analyses and Hypothesis Testing), in which bPmin is the scaling coefficient that describes how function (here, Pmin) varies with species richness; higher absolute values of bf (where f is an ecosystem function) indicate a steeper relationship between biodiversity and function and can be used to compare “biodiversity effects” among studies and systems (12, 52). As species richness increased in potential diets, Pmin declined, and RDA targets for each micronutrient or fatty acid were achieved with less total seafood intake (Fig. 3A, bPmin < 0 for every micronutrient and fatty acid: calcium −0.32 [95% CI −0.35, −0.28], iron −0.24 [95% CI −0.27, −0.22], zinc −0.26 [95% CI −0.28, −0.23], EPA −0.25 [95% CI −0.27, −0.23], and DHA −0.22 [95% CI −0.23, −0.21]). Increasing species richness reduced the minimum portion size required, Pmin, in our sample diets, independent of systematic changes in the identity of species included (Stastical Analysis and Hypothesis Testing), because the diets were assembled using random samples of the species pool. The restricted variation and symmetrical distribution in protein concentrations (Fig. 2 Upper), combined with high levels of protein in all edible tissues, lead to no benefit and a minimal detrimental effect of seafood species richness on protein provisioning (Fig. 3A, bPmin = 0.0071 95% CI 0.0062, 0.0080). In other words, the ecosystem service of protein provisioning was adequately provided by total seafood edible biomass and not improved by species richness or even species identity. The findings for the micronutrients and fatty acids are consistent with demonstrations that variety and diversity in diets is important for nutrition (28, 29), but we extend these findings to show seafood species richness allows consumers to gain more nutritional benefit without consuming more total seafood biomass, and explicitly relate this pattern with general effects of biodiversity in ecological systems.
Aquatic biodiversity enhances nutritional benefits. (A) Seafood species richness improves the efficiency with which human diets can meet RDA targets by reducing the minimum portion size required, Pmin, to meet RDA targets (measured in grams of seafood). Pmin is shown for micronutrients, fatty acids, and protein separately (points are median values for calcium, iron, zinc, EPA, DHA, and protein, lines show the fit of Eq. 3 to the data, and shading refers to 95% CI) as well as for the five micronutrients and fatty acids simultaneously (top purple line). Colors corresponding to each nutrient in A are shown in B. (B) Estimates (± 95% CI) for the scaling parameter that relates species richness to Pmin (bPmin) (Eq. 3). (C) Species richness increases the number of distinct nutrient RDA targets met in a 100 g seafood portion (NT); black line and 95% CI correspond to the fit of Eq. 2 to mean NT derived from resampled diets from the global seafood species pool. Flower plots in C summarize the micronutrient and fatty acid concentrations relative to 10% RDA (gray circle) in two representative diets at low and high species richness levels. Data shown in A and C are derived from n = 1,000 resampled diets.
We then considered the effects of seafood species richness on the provisioning of multiple nutrients simultaneously (Metrics). This is referred to as a multifunctional benefit of biodiversity (53, 54) and takes into account possible trade-offs or correlations among functions; in this case, concentrations of micronutrients and fatty acids. For some ecosystem services (e.g., water quality or ecotourism), benefits of biodiversity accumulate when multiple ecosystem functions are considered simultaneously (54⇓–56). Consistent with biodiversity–ecosystem functioning theory, we found that in the case of a multifunctional metric of an ecosystem service defined from the human beneficiary’s perspective (i.e., multiple micronutrient and fatty acid targets reached simultaneously), biodiversity benefits for the multifunctional service are greater than for individual functions (bPmin for all five micronutrients and fatty acids simultaneously = −0.42 [95% CI −0.47, −0.38] versus single nutrients bPmin range from −0.32 [95% CI −0.35, −0.28] for calcium to −0.22 for DHA [95% CI −0.23, −0.21]) that comprise the ecosystem service (Fig. 3 A and B). Increasing seafood species richness from one to 10 species in 1,000 simulated, resampled diets drawn from our global species pool allowed diets to meet RDA targets for five essential microelements and fatty acids simultaneously more than twice as efficiently (i.e., a median of 494.19 g of tissue required with one species versus median of 213.34 g of tissue required with 10 species) (Fig. 3 A and B). Then, we assessed the effects of biodiversity when the total biomass of seafood consumed was held constant, at 100 g, by counting the number of nutrients for which RDA targets were reached in a 100 g portion (NT). We found positive effects of biodiversity on the number of nutrients that met RDA targets, NT, in a single 100 g portion (Fig. 3C): diets with higher species richness reached more nutritional targets (higher NT) per 100 g serving than diets of the same fish biomass comprising fewer species (bNT = 0.20 [95% CI: 0.20 to 0.21]; Fig. 3C). These findings were robust to different RDA target levels (for Pmin, they were independent of RDA target level, and for NT, they were positive over the range from 1 to 40% RDA per 100 g portion; SI Appendix, Fig. S2). These results reveal biodiversity effects of seafood quantitatively comparable with the widely recognized relationship between biodiversity and productivity (26, 47, 52) and demonstrate a benefit of biodiversity for human nutritional well-being over and above the benefits of consuming a particular amount (biomass) or identity of aquatic species.
Increasing Seafood Diversity Increased Contaminant Exposure.
We considered a range of trace elements, which, at high concentrations, are known to be harmful to human health (57). We focused on four heavy metals considered as contaminants (methylmercury, cadmium, arsenic, and lead) for which there exist public health guidelines for upper tolerable limits and for which seafood is potentially a major source of dietary intake (SI Appendix, Fig. S1 and Table S2). We examined the concentrations of these elements in the muscle tissues of 353 seafood species (thereby excluding parts such as viscera, liver, and bones). We found that the same mechanisms that lead to a positive relationship between seafood species richness and nutritional benefits also contributed to exposure to a wider range of contaminants with higher species richness. We observed high levels of variation in contaminant concentrations across species (Fig. 4 A–D), and as we observed for the micronutrients and fatty acids, these distributions were often right skewed such that most species contained low contaminant concentrations, and few contained high concentrations. On average, muscle tissue concentrations of contaminants were weakly positively correlated across species, such that species that contained high concentrations of one contaminant also contained high concentrations of another (average pairwise correlation = 0.09) (SI Appendix, Fig. S3). When considering multiple contaminants simultaneously, increasing seafood species richness increased the number of contaminants that exceeded their upper tolerable limits (PTDI) in a 100 g portion, referred to as NC (SI Appendix, Table S1, Eq. A4), (bNC = 0.10, 95% CI 0.084, 0.12) (Fig. 4E). When we considered the effects of biodiversity on exposure to each contaminant separately, we found that increasing species richness generally increased contaminant content per 100 g, but the strength of this effect varied among contaminants. For example, increasing species richness from one to 10 species was associated with doubling methylmercury concentrations on average, thereby reducing the maximum portion size before exceeding PTDI, referred to as Pmax (SI Appendix, Table S1, Eq. A3) (bPmax = −0.25 95% CI −0.27, −0.22). For lead, however, the biodiversity effect was weaker (an order of magnitude smaller than for methylmercury). Increasing species richness from one to 10 species was associated with only a 10% increase in average lead concentrations and a small reduction in maximum portion size, Pmax (bPmax = −0.039, 95% CI −0.049, −0.034). For cadmium, bPmax = −0.12 (95% CI −0.13, −0.12), and for arsenic, bPmax = −0.17 (95% CI −0.19, −0.14). These differences in the strength of biodiversity effects among contaminants are linked to the shape of the distribution of their concentrations among fish species (i.e., normal versus lognormal) (Fig. 4 A–D), and the more skewed the distribution, the stronger the biodiversity effect.
Frequency histograms of concentration of arsenic (A), cadmium (B), lead (C), and methylmercury (D) in edible muscle tissues of North American aquatic species relative to PTDI. (E) Increasing seafood species richness increases the number of contaminants which exceed the upper tolerable limit (PTDI) in a 100 g portion (NC). Black line indicates the mean NC from 1,000 bootstrapped samples of seafood diets sampled from the North American seafood contaminant dataset, and gray shading refers to 95% CIs. Slope, bNC = 0.10, 95% CI 0.084, 0.12.
Biodiversity Benefits Were Consistent at Local and Global Scales.
Consistent with the positive biodiversity effects we observed when sampling diets from a global seafood species pool, we also found benefits of seafood diversity in a local context. We analyzed the effects of biodiversity on nutrient content in 14 traditional Indigenous North American diets of seafood harvested and consumed locally. We found a consistent, beneficial effect of biodiversity on NT and Pmin, although the magnitude of the biodiversity effect was generally lower at the local scale than the global scale (SI Appendix, Fig. S4) (global bPmin = −0.42 [95% CI −0.47, −0.38] versus mean local bPmin = −0.25 ±0.0091 SE and global bNT = 0.20 [95% CI 0.20, 0.21] versus mean local bNT = 0.097 ±0.0082 SE). This finding is consistent with lower nutritional functional diversity (NFD, SI Appendix, Table S1) sensu ref. 34 (mean local NFD = 2.77 ±0.17 SE versus global NFD = 3.87) and higher nutritional functional evenness in local diets (mean local NFEve = 0.82 ±0.0037 SE versus global NFEve = 0.57) (SI Appendix, Methods 3 and Fig. S5), suggesting that functional consequences of changes to diversity in local seafood diets may be buffered by higher nutritional redundancy among species.
Nutritional Traits Covary with Ecological Traits.
We found that nutrient concentrations varied substantially among species in ways that differed for different nutrients (Fig. 2)—a condition necessary for biodiversity per se to increase nutritional benefits. The diversity that we observed in the nutrient content of edible portions (Fig. 2) was partly explained by ecological attributes and functional traits: habitat, trophic position, body size, diet source, and feeding mode (SI Appendix, Tables S4–S8). When considering all five micronutrients and fatty acids (calcium, iron, zinc, EPA, and DHA) together, finfish, crustaceans, and molluscs differed significantly in their multinutrient profiles (SI Appendix, Table S1, permutational multivariate ANOVA [PERMANOVA], F 2,103 = 3.429, P = 0.006). Among finfish, nutrient concentrations depended on which tissues were included in the edible portions (significant “body part” effect shown in SI Appendix, Tables S4 and S6–S11). Finfish species whose edible portions included organs such as liver or bones had higher nutrient concentrations in the edible portion than those whose edible portions were restricted to muscle tissue (ANOVA P < 0.01 for calcium, iron, and zinc concentrations; SI Appendix, Fig. S6). Principal components analysis of multiple nutrient concentrations in species showed that essential element (calcium, iron, and zinc) concentrations were typically negatively correlated with essential fatty acid concentrations (EPA and DHA) (SI Appendix, Fig. S7), allowing complementarity among species to increase nutritional benefits. Specifically, high EPA and DHA concentrations traded off against low calcium and zinc concentrations, and vice-versa (negative pairwise Pearson correlation coefficients; SI Appendix, Fig. S7). When considering muscle tissues or muscle and skin tissues of finfish only (thereby eliminating the influence of body parts such as bones on nutrient concentration), concentrations of calcium, iron, zinc, EPA, and DHA were associated with ecological traits across species, including habitat and diet source (e.g., demersal versus pelagic), body size, and trophic position (Fig. 5). Relationships between species’ nutrient tissue concentrations and their habitats and trophic positions have been predicted by ecological stoichiometry theory (58). Relationships between ecological traits and nutrient concentrations differed for different nutrients. For example, tissue concentrations decreased with body size for calcium but not for the other microelements, EPA or DHA (Fig. 5 and SI Appendix, Figs. S8–S10 and Tables S4–S8). Species at lower trophic positions had higher zinc concentrations in their tissues than species at higher trophic positions, but we did not observe this relationship for other nutrients (SI Appendix, Tables S4–S8). These examples illustrate that trade-offs and variation in nutrient concentrations across species were associated with variation in different ecological traits and roles that species play in ecosystems.
Nutrient concentrations in finfish muscle tissue vary with ecological traits in ways that differ among the essential trace elements (calcium, iron, zinc) and the essential fatty acids (EPA and DHA). Model-averaged standardized regression coefficients and 95% CIs from phylogenetic least squares regression (see SI Appendix, Methods 4 for full model description) are shown for samples including muscle or muscle and skin tissues only. Traits for which there is no symbol did not appear in any of the models in the top 95% set (with cumulative sum of Akaike weights ≤ 0.95). Open symbols indicate coefficient estimates for which 95% CIs do not encompass zero. Note that x-axes differ across panels for clarity of presentation. Number of species: n = 155 for EPA, n = 159 for DHA, n = 104 for calcium, n = 99 for iron, and n = 90 for zinc. For model results including tissues other than muscle tissue, see SI Appendix, Tables S4–S11.
Nutritional Value Was Linked to Diversity of Ecological Functions.
Ecological functional diversity of a species assemblage captures variation in traits and ecological roles of species (48, 49) and is understood to play an important role in the relationship between biodiversity and ecosystem function (Fig. 6A). We assessed the relationship between ecological functional diversity of seafood species diets (independent from species richness) and nutritional value. Consistent with observations for ecosystem functions such as productivity and biomass, nutritional benefits and ecological functional diversity were positively related, such that seafood diets with higher ecological functional diversity also provided higher nutritional value (i.e., were more likely to reach five micronutrient and fatty acid RDA targets simultaneously, NT = 5) (Fig. 6B). Ecological functional diversity increased with species richness (SI Appendix, Fig. S11), and higher levels of ecological functional diversity were also associated with lower minimum portion size required (SI Appendix, Fig. S12). Because aquatic assemblages with higher ecological functional diversity have been shown to exploit more diverse resources, transform and transport energy and materials more efficiently, produce higher yields, and be more productive and resilient over time (42, 59⇓–61), it is possible that the provisioning of multiple micronutrients and fatty acids occurs in tandem with a range of other ecological functions.
(A) Ecological functional diversity, EFD, captures variation in the ecological traits and ecological roles of species and can be quantified as the average of the distances between species’ trait values (circles) and the center of functional trait space (square, indicated by z, adapted from ref. 93). (B) Seafood diets with higher levels of EFD are associated with higher levels of nutritional benefits (i.e., higher number of RDA targets per 100 g portion). Probability of reaching five micronutrient and fatty acid RDA targets simultaneously in a single 100 g portion (NT = 5, light orange line) increases as the EFD of the diet increases. Probabilities predicted from ordinal logistic regression (log odds of functional diversity = 2.08, 95% CI [0.52, 3.64], n = 1,000 resampled diets).
Discussion
Biodiversity of seafood provides high levels of nutritional benefits to humans because nutrient concentrations vary substantially among species in ways that differ for different nutrients (Figs. 2 and 5 and SI Appendix, Fig. S7). The effects of species richness observed for nutritional benefits equal or exceed mean observed diversity effects documented for plant and forest species richness and productivity (12, 40, 52, 62). Our findings build on evidence showing that biodiversity in marine systems enhances ecological functions by extending this paradigm to human nutrition. In agro-ecosystems, specific combinations of species such as corn, pumpkin, and beans have been planted to exploit complementary traits (shade provisioning, nitrogen fixation, and biomass production) to attain higher yields, more resilient crops, and enhanced nutritional benefits (35). By demonstrating that nutritional benefits of seafood can be understood as a consequence of seafood species richness and ecological functional diversity, we have shown that diverse seafood assemblages also provide nutritional and other ecological benefits simultaneously.
Our analysis provides robust evidence that biodiversity is critical to multifunctionality of ecosystem services when function thresholds are grounded a priori in multivariate metrics meaningful for human well-being such as RDA. Our approach overcomes the critique that multifunctionality is not enhanced by biodiversity (63) but rather a statistical artifact of how multifunctionality was commonly estimated. Our findings are robust to a range of RDA target levels (SI Appendix, Fig. S2). In the case of Pmin, we found that the benefit of biodiversity was consistent across all RDA target levels considered, and in the case of NT, we found beneficial effects of biodiversity across a range of levels of nutritional value that are significant for human nutrition (i.e., 1 to 40% RDA per 100 g portion), highlighting the importance of species richness in seafood diets. More generally, ecosystem service benefits, as defined in metrics of human well-being rather than the traits of the species pool under consideration (e.g., biomass or stability of the food web), typically are produced by several underlying ecosystem functions (54). The strong effects of diversity on multifunctional benefits observed here may also apply to relationships between diversity and other services, for example, desired filtration rates of pollutants in wetlands (64), or desired pest consumption rates in agricultural systems (65).
Physiological and Ecological Mechanisms Underlying Nutritional Benefits and Risks.
When comparing the magnitude of the biodiversity–ecosystem function scaling parameter b among the nutrients and contaminants, we found that the magnitude of b (Eq. 1) was higher for nutrients with more skewed distributions of tissue concentrations (Figs. 2 and 4 and SI Appendix, Fig. S7). For all of these highly skewed distributions, we observed a strong biodiversity effect (i.e., high values of b). However, for protein and lead in muscle tissues, we observed symmetrical tissue concentration distributions. In these two cases, the biodiversity effect was either nonexistent (in the case of protein) or very weak (in the case of lead). This finding suggests that biodiversity effects may increase when the trait distribution is skewed and includes some species with extreme trait values, known as “functionally specialized” and “functionally unique” species (9, 66, 67). Understanding the drivers of the tissue concentration distribution among species may lend insights into how biodiversity effects may change as the nutrient under consideration changes or as the environment changes.
Increasing seafood species richness increases nutritional benefits as well as contaminant exposure. Increasing biodiversity can have negative consequences for some ecosystem functions, despite positive consequences for others. On balance, the benefits of diversity may outweigh the negative effects on function when considering multiple functions together (68). Here, we found that the same mechanisms (lognormal and complementary tissue concentration distributions among species) that contributed to the positive relationship between biodiversity and human nutritional benefits also applied to contaminant exposure. As a result, increasing seafood species richness comes with both benefits and risks (69⇓–71). However, interpreting our results in the context of public health outcomes is complicated by the fact that epidemiological evidence on health outcomes of contaminant exposure is mixed and likely dependent on complex social and health risks. The health risks of contaminant exposure depend on other health factors such as smoking or nutrients in the diet (72, 73) and disease status. Complex interactions among multiple diet and health risk factors (74) were beyond the scope of this study but would be necessary to understand exposure risks from seafood consumption in any particular community. Nonetheless, our results suggest that while increasing biodiversity increases contaminant exposure, it also increases nutritional benefit and reduces the portion sizes required to meet nutritional demands. Finding a balance between seafood biodiversity, seafood biomass consumption, and the resulting risks and benefits will be critical for both human and ecosystem health.
Aquatic Biodiversity and Food Security in a Changing World.
Seafood-derived nutrition plays an important role in food security. The link that we have demonstrated between seafood biodiversity (species richness and ecological functional diversity) and nutrition in an ecological framework unites three of the United Nations Sustainable Development Goals focused on biodiversity, hunger, and well-being (75). More than two billion people suffer from micronutrient deficiencies (76, 77), and many of the most nutritionally vulnerable populations—those that are deficient in essential micronutrients and fatty acids during particularly sensitive stages of life (i.e., pregnancy, breastfeeding, and childhood)—may rely heavily on local aquatic ecosystems to meet their nutritional demands (15, 19, 22, 78). These populations may have access to a limited amount of locally available seafood tissue each day, suggesting that nutritional efficiency (i.e., lower Pmin) provided by biodiversity in wild-caught seafood may be particularly important for these populations. Regions of high nutritional vulnerability continue to experience major changes in biodiversity and ecosystem structure (6, 17, 79), and climate change is compounding threats to the sustainability of capture fisheries (2, 80). Our results suggest that declines in these aspects of diversity in wild aquatic ecosystems could make achieving sustainability goals for food security via seafood even more difficult. Seafood diets composed of more species, and groups of species with higher levels of ecological functional diversity, are more likely to provide more nutrients per unit biomass than less diverse seafood diets while also maintaining high levels of ecosystem function (26). This finding bridges the growing understanding of hidden hunger and food security with the large theoretical and empirical understanding of relationships between biodiversity, ecosystem function, and benefits to people. Biodiversity in natural aquatic systems can be maintained by reducing pollution and overharvest and by allowing ecosystems to adapt to climate change, and these measures could also benefit humanity directly through seafood provisioning.
Conclusions
Nutritional value appears to be derived from ecological diversity, suggesting links between the complexity of aquatic ecosystems and their capacity to produce nutritional benefits. While the role of seafood is well recognized as an important source of protein in the human diet, the role of seafood biodiversity as an important aspect of the provision of essential micronutrients has been overlooked. Our results reveal that aspects of ecological structure including species and ecological functional diversity enhance nutritional benefits while also increasing contaminant exposure, thereby linking the processes that structure ecosystems with their potential benefits and risks to human nutrition and health.
Methods
To test our hypotheses about aquatic species diversity and potential nutritional benefits for human well-being, we assembled a database of nutritional values by synthesizing observations from existing data (SI Appendix, Methods and Fig. S1). To build the database, we identified quantitative and comparable measures of nutritional content for the edible portions of aquatic species (Quantifying Nutritional Value for Single Nutrients). We then identified metrics for relating nutrient content to human health (Defining Nutritional Benefits and Risks for Multiple Nutrients or Contaminants Across Diverse Species Groups). We repeated this exercise for contaminants (Quantifying Contaminant Exposure in Terms of Human Health Risks). Next, we adapted an approach from biodiversity–ecosystem functioning theory for quantitatively assessing the potential nutritional value of seafood diets varying in their species composition and diversity for multiple nutrients separately and simultaneously (Modeling the Biodiversity Effect for Sample Diets). To test predictions of biodiversity–ecosystem functioning theory, we simulated potential diets at global and local scales and fit models to test for effects of species richness and ecological functional diversity (Simulating Seafood Diets and Estimating Nutritional Benefits at Different Seafood Species Richness Levels; and Testing Hypotheses That Biodiversity Enhances Nutritional Benefits). SI Appendix, Fig. S1 provides a graphical overview of our analyses.
Metrics.
Quantifying nutritional value for single nutrients.
We characterized an aquatic species’ nutritional value by drawing on two well-established nutritional metrics: nutrient concentration (mg/100 g edible portion) and RDAs. RDAs are developed following health guidelines to quantify the recommended amount of a particular nutrient required to maintain health (50). We used RDAs established by the Food and Nutrition Board of the US Institute of Medicine to quantify the daily intake level of a nutrient needed to meet the requirements of 97 to 98% of healthy adults (females aged 19 to 50 y) (50) (SI Appendix, Table S2). We refer to an RDA-defined threshold as an RDA target (SI Appendix, Tables S1 and S2). We calculated a ratio of the nutrient content in a 100 g edible portion relative to the RDA (or fraction thereof) for that nutrient. For many species, nutrient concentrations in edible tissues provided only small fractions of the RDA (SI Appendix, Table S2). Following the Institute of Medicine, we chose an RDA target as 10% of RDA for a given nutrient in a single portion (SI Appendix, Table S2) because this is a minimum threshold for a food to be considered of nutritional benefit (50). We defined the minimum portion size required, Pmin, as the minimum amount of edible seafood tissue (g) required to reach the RDA target for a given nutrient (or a set of nutrients, Defining Nutritional Benefits and Risks for Multiple Nutrients or Contaminants Across Diverse Species Groups) (SI Appendix, Table S1, Eq. A1). We quantified the sensitivity of our biodiversity findings to different threshold levels of RDA for both single nutrient analyses and multinutrient analyses and found that for Pmin, they did not vary with threshold level, and for NT, they varied with threshold but were robust to RDA threshold levels between 1 and 40% RDA (SI Appendix, Fig. S2).
We examined data on concentrations of macronutrients including protein and fat as well as five micronutrients and essential fatty acids (n = 5,041 observations of nutrient concentrations, SI Appendix, Fig. S1A): metals beneficial at low concentrations but toxic at high concentrations (zinc and iron), one beneficial mineral (calcium), and the polyunsaturated fatty acids EPA and DHA. We chose these five micronutrients and fatty acids because we required that RDA standards exist for each one (50) and that each nutrient has known functions in organismal physiology and is considered “biologically essential” because it is required by organisms to grow or reproduce. The concentration of biologically essential nutrients in organisms’ tissues are controlled homeostatically by organisms (81, 82) and therefore might be biologically related to ecological trait values we considered in our trait analysis.
Quantifying contaminant exposure in terms of human health risks.
We characterized an aquatic species’ contaminant content relative to established public health guidelines for exposure. We used the provisional tolerable weekly intake (PTWI) developed by the Food and Agriculture Organization of the United Nations (FAO)/World Health Organization Expert Committee on Food Additives (57), which estimates the amount of a substance in air, food, soil, or drinking water that can be assimilated weekly per unit body weight over a lifetime without appreciable health risk (57). To make this metric comparable to RDA targets defined above, we quantified contaminant exposure per day by dividing the PTWI by seven and using the body weight of a 70-kg person to calculate a daily tolerable limit for use in our analyses, which we refer to as provisional tolerable daily intake (PTDI; SI Appendix, Table S1). In our analyses, we consider a tolerable upper limit to be exceeded if an edible portion contains 100% or more of the PTDI. We examined the sensitivity of this choice of 100% threshold by considering lower values (e.g., 50%) and found that for all values between 50 and 100% of PTDI, our findings remain consistent.
Defining nutritional benefits and risks for multiple nutrients or contaminants across diverse species groups.
To test our hypothesis that nutritional value may depend on seafood diversity in a way that considers multiple nutrients at the same time, we used two metrics that considered multiple micronutrients and fatty acids simultaneously (Fig. 1): 1) the minimum portion size introduced above, Pmin, (SI Appendix, Table S1, Eq. A1), which quantifies the amount of tissue, in grams, required to reach the RDA target for five nutrients simultaneously and 2) the number of nutrients for which RDA targets are met in a standard 100 g edible portion, NT (SI Appendix, Table S1, Eq. A2). Note that our NT differs from another measure, NFD (34), which is a measure of nutritional profile in multivariate trait space and does not explicitly quantify the number of nutrient RDA targets met in a portion. Our NT is a measure of potential nutritional value and is derived from the biodiversity–ecosystem function perspective to allow us to compare nutritional value with other functions that depend on biodiversity (e.g., productivity, resource cycling). This measure does not consider the potential physiological interactions or co-benefits of consuming multiple nutrients that would be considered in metabolic models for human health. To test our hypothesis that contaminant exposure also varies with seafood diversity, we defined Pmax as the amount of tissue beyond which a tolerable upper limit (PTDI) for a given contaminant would be reached (SI Appendix, Table S1, Eq. A3 and Methods 2). Analogous to the NT metric defined above for the nutrients, we defined the number of distinct PTDIs in a 100 g portion as NC (SI Appendix, Table S1, Eq. A4 and Methods 2). We chose to standardize NT and NC relative to 100 g of seafood per day because 100 g edible portion is a standard metric used in the food and nutrition literature (37, 83), and consumption rates of 100 g per day is within the range of daily consumption rates for many communities that rely heavily on seafood to meet nutritional needs and for whom seafood may also pose substantial contaminant exposure risks (16, 45, 84, 85).
Statistical Analyses and Hypothesis Testing.
Modeling the biodiversity effect for sample diets.
To quantitatively compare effects of biodiversity among possible seafood diets varying in species richness, we modeled species richness effects as a power function:
where the parameter b is referred to as the “biodiversity effect” and describes the relationship between a change in species richness, S, and a measure of function, Y, such as Pmin or NT, where a is a constant (12, 52).
Simulating seafood diets and estimating nutritional benefits at different seafood species richness levels.
We tested the effect of species richness, S, on nutritional value by randomly assembling diets from the global seafood species pool at varying levels of species richness (SI Appendix, Fig. S1B). In our analyses of NT, we kept the total biomass constant at each level of species richness (100 g). In analyses using the global dataset (547 species) (SI Appendix, Methods 1.1), we assembled diets from the entire global species pool, choosing species at random without replacement. This way of simulating diets certainly ignores economic, social, and cultural factors that affect which species people consume but allows us to consider the potential effect of biodiversity on diets before diets are filtered by these other processes. To assess potential effects of biodiversity on nutritional value for populations that consume seafood locally from a restricted species pool, we sampled diets from species contained within traditional diets in 14 Indigenous cultures in North America (SI Appendix, Methods 1.4).
We created sample diets by sampling 10 species at random from the global species pool and then assembling seafood diets from all possible combinations of these 10 randomly chosen species at 10 levels of species richness (1 to 10 species) to generate 1,023 simulated diets (SI Appendix, Fig. S1 B, i). At each level of species richness, we assembled diets following the typical experimental design employed to test the hypothesis that biodiversity affects ecosystem functioning, analogous to a biodiversity–ecosystem function experiment with a replacement design (86), where species’ abundances in the diet decline proportionally as species richness increases such that each species contributed an equal proportion of biomass. For each diet at each level of species richness (i.e., n = 1,023), we calculated Pmin (SI Appendix, Table S1, Eq. A1 and Fig. S1 C, ii; for either, one of six possible nutrients targets individually [protein, calcium, iron, zinc, EPA, and DHA] or five micronutrient and fatty acid [calcium, iron, zinc, EPA, and DHA] targets simultaneously). To estimate NT (SI Appendix, Fig. S1 C, ii), we quantified the number of distinct nutrient RDA targets by assigning each diet (n = 1,023) a set of zeros or ones according to whether that combination met the RDA target for each nutrient (SI Appendix, Table S1, Eq. A2). This approach allowed us to explore how likely it would be for potential human diets containing different numbers of seafood species to reach RDA targets for a given number of micronutrients and fatty acids (NT ranges between 0 and 5), assuming that seafood species were included in the human diet at different levels of richness at random. At each level of species richness, we averaged Pmin and NT. We then repeated this process of random sampling 10 species from the species pool, assembling diets at each level of richness, estimating metrics of nutritional value, and averaging at each richness level 1,000 times, yielding 1,000 estimates of each metric at each richness level (SI Appendix, Fig. S1 C, iii).
Testing hypotheses that biodiversity enhances nutritional benefits.
We tested the hypothesis that complementarity in nutrient concentrations among species increases nutritional benefits by increasing NT (SI Appendix, Fig. S1D). We quantified the effect of seafood species richness, S, on NT at each richness level estimated in Simulating Seafood Diets and Estimating Nutritional Benefits at Different Seafood Species Richness Levels (n = 1,000 estimates of NT per richness level) by fitting a power function of the form shown in Eq. 1,
where the parameter bNT describes the relationship between a change in species richness, S, and a change in NT, and a is a constant.
We tested the hypothesis that complementarity in nutrient concentrations among species reduces the minimum portion size required, Pmin, by estimating the effect of species richness, S, on Pmin (n = 1,000 estimates of Pmin per richness level) at each richness level using
where the parameter bPmin describes the relationship between a change in species richness, S, and a change in Pmin, and a is a constant (SI Appendix, Fig. S1D).
We estimated a and b in Eqs. 2 and 3 using nonlinear regression using the nls.lm function in the minpack.lm package in R (87). We conducted all analyses in R version 3.3.2 (88). To quantify uncertainty in parameter estimates associated with the fit of Eqs. 2 and 3 to our estimates of function (Simulating Seafood Diets and Estimating Nutritional Benefits at Different Seafood Species Richness Levels), we calculated bootstrapped CIs (n = 1,000 bootstraps) using nonparametric bootstrapping of mean centered residuals using the nlsBoot function in the R package nlstools (89). For both Pmin and NT, we tested the hypothesis that biodiversity enhances nutritional benefits by assessing whether the estimate of the scaling exponent, b, had CIs not overlapping zero.
We tested the hypothesis that nutrient concentrations are related to species’ ecological traits in two ways: 1) testing whether multinutrient profiles (i.e., concentrations of all five micronutrients and fatty acids) differ among major phylogenetic groups using PERMANOVA (SI Appendix, Methods 4) and 2) whether differences in single nutrient concentrations differ with species’ ecological traits using phylogenetic least squares regression (SI Appendix, Methods 4). We quantified the relationship between nutritional benefits and ecological functional diversity (EFD, SI Appendix, Table S1) by estimating NT, Pmin, and EFD, estimated as functional dispersion, of diets simulated from the global species pool (SI Appendix, Methods 5).
Uncertainties.
There are several sources of uncertainty in our analyses. First, there are substantial sources of uncertainty in food composition estimates. The data in our dataset meet international standards for data quality and standardization, meaning that we followed guidelines for checking food composition data and converting units, denominators, and expressions (90). Still, tissue concentrations may vary depending on analytical techniques, laboratories, season, diet of the animal, life stage, etc. Some of these sources of uncertainty (e.g., differences in analytical techniques) are unavoidable consequences of synthesizing previously published data collected across many laboratories. We assumed that these uncertainties in the data were randomly distributed over our geographically and taxonomically diverse dataset. Further uncertainty is associated with how well our set of 547 species represents the global pool of seafood consumed. We do not know whether our sample is random or biased, though we can say that our dataset includes 41 of the 67 most consumed species worldwide [as determined by FAO production volumes (91), species with capture production of 150,000 tons or more, after removing species for which the majority of production volume is diverted to fish meal and oil (92), SI Appendix, Table S13]. A remaining source of variation among samples is likely due to natural sources of variation associated with seasonal and other sources of temporal variability, which we consider to be an important component of biodiversity.
Data Availability
Data and code are available at GitHub (https://github.com/JoeyBernhardt/Nutrient_analysis) and are archived using Zenodo with DOI: 10.5281/zenodo.4474988. Data are available on Dryad (https://doi.org/10.5061/dryad.rn8pk0p8t) (94).
Change History
April 07, 2021: Due to a production error, the wrong Fig. 1 was published. It has been updated to the correct version.
Acknowledgments
We thank P. Thompson, A. Gehman, W. Cheung, M. Ruckelshaus, C. Harley, M. Deith, L. Dee, and F. Isbell for helpful comments on previous versions of this manuscript. Portions of this manuscript were adapted from and based on work contained in the doctoral thesis of J.R.B. from the University of British Columbia. Funding was provided by a Vanier Canada Graduate Scholarship and a Four Year Fellowship from the University of British Columbia to J.R.B.
Footnotes
- ↵1To whom correspondence may be addressed. Email: joey.bernhardt{at}biodiversity.ubc.ca.
Author contributions: J.R.B. conceived and designed the study, developed software, collected and curated the data; J.R.B. analyzed and visualized the data with input from M.I.O.; J.R.B. wrote the paper with help from M.I.O.; and J.R.B. and M.I.O. supervised the study.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1917487118/-/DCSupplemental.
- Copyright © 2021 the Author(s). Published by PNAS.
This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).
References
- ↵
- ↵
- M. L. Pinsky,
- B. Worm,
- M. J. Fogarty,
- J. L. Sarmiento,
- S. A. Levin
- ↵
- D. J. McCauley et al
- ↵
- J. Bélanger,
- D. Pilling
- FAO
- ↵
- E. S. Brondizio,
- J. Settele,
- S. Díaz,
- H. T. Ngo
- IPBES
- ↵
- G. P. Jones,
- M. I. McCormick,
- M. Srinivasan,
- J. V. Eagle
- ↵
- ↵
- J. A. Hutchings,
- R. A. Myers
- ↵
- C. Pimiento et al
- ↵
- S. S. Myers et al
- ↵
- ↵
- M. I. O’Connor et al
- ↵
- E. Chivian
- ↵
- FAO
- ↵
- H. V. Kuhnlein,
- O. Receveur
- ↵
- A. M. Cisneros-Montemayor,
- D. Pauly,
- L. V. Weatherdon,
- Y. Ota
- ↵
- J. G. Molinos et al
- ↵
- WHO
- ↵
- ↵
- B. Belton,
- S. H. Thilsted
- ↵
- FAO
- ↵
- N. Kawarazuka,
- C. Béné
- ↵
- ↵
- D. Tilman et al
- ↵
- ↵
- ↵
- M. T. Ruel
- ↵
- E. A. Frison,
- J. Cherfas,
- T. Hodgkin
- ↵
- J. Fanzo,
- D. Hunter,
- T. Borelli,
- F. Mattei
- ↵
- ↵
- ↵
- A. D. Jones
- ↵
- C. Lachat et al
- ↵
- ↵
- ↵
- ↵
- R. Remans,
- S. A. Wood,
- N. Saha,
- T. L. Anderman,
- R. S. DeFries
- ↵
- S. A. Wood
- ↵
- ↵
- P. B. Reich et al
- ↵
- ↵
- ↵
- R. Hilborn,
- T. P. Quinn,
- D. E. Schindler,
- D. E. Rogers
- ↵
- ↵
- J. R. Bogard et al
- ↵
- ↵
- ↵
- D. Tilman
- ↵
- ↵
- IOM
- ↵
- ↵
- B. J. Cardinale et al
- ↵
- ↵
- P. Manning et al
- ↵
- ↵
- ↵
- WHO
- ↵
- R. W. Sterner,
- J. J. Elser
- ↵
- J. S. Lefcheck et al
- ↵
- ↵
- J. S. Lefcheck,
- J. E. Duffy
- ↵
- J. Liang et al
- ↵
- L. Gamfeldt,
- F. Roger
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- M. O. Gribble et al
- ↵
- UN
- ↵
- ↵
- J. Fanzo et al.
- ↵
- J. R. Bogard et al
- ↵
- W. W. L. Cheung et al
- ↵
- L. A. Rogers et al
- ↵
- ↵
- M. T. Arts,
- M. T. Brett,
- M. Kainz
- ↵
- R. DeFries et al
- ↵
- N. Roos,
- M. A. Wahab,
- C. Chamnan,
- S. H. Thilsted
- ↵
- EPA
- ↵
- ↵
- T. V. Elzhov,
- K. M. Mullen,
- A.-N. Spiess,
- B. Bolker
- ↵
- R Core Team
- ↵
- ↵
- Food and Agriculture Organization/International Network of Food Data Systems (INFOODS)
- ↵
- Food and Agriculture Organization of the United Nations
- ↵
- T. Cashion,
- F. Le Manach,
- D. Zeller,
- D. Pauly
- ↵
- ↵
- J. R. Bernhardt,
- M. I. O'Connor
Citation Manager Formats
Article Classifications
- Biological Sciences
- Sustainability Science