Search

Methane (CH 4 ) is a potent greenhouse gas but also an important carbon and energy substrate for some lake food webs. Understanding how CH 4 incorporates into food webs is, therefore, crucial for unraveling CH 4 cycling and its impacts on climate and ecosystems. However, CH 4 -fueled lake food webs from pre-Holocene intervals, particularly during greenhouse climates in Earth history, have received relatively little attention. Here, we present a long-term record of CH 4 -fueled pelagic food webs across the Cretaceous Oceanic Anoxic Event 1a (~120 Mya) that serves as a geological analog to future warming. We show an exceptionally strong expansion of both methanogens and CH 4 -oxidizing bacteria (up to 87% of hopanoid-producing bacteria) during this Event. Grazing on CH 4 -oxidizing bacteria by zooplankton (up to 47% of ciliate diets) within the chemocline transferred substantial CH 4 -derived carbon to the higher trophic levels, representing an important CH 4 sink in the water column. Our findings suggest that as Earth warms, microbial CH 4 cycling could restructure food webs and fundamentally alter carbon and energy flows and trophic pathways in lake ecosystems., Competing Interests: Competing interests statement:The authors declare no competing interest.

Indeterminacy of ecological networks-the unpredictability of ecosystem responses to persistent perturbations-is an emergent property of indirect effects a species has on another through interaction chains. Thus, numerous indirect pathways in large, complex ecological communities could make forecasting the long-term outcomes of environmental changes challenging. However, a comprehensive understanding of ecological structures causing indeterminacy has not yet been reached. Here, using random matrix theory (RMT), we provide mathematical criteria determining whether network indeterminacy emerges across various ecological communities. Our analytical and simulation results show that indeterminacy intricately depends on the characteristics of species interaction. Specifically, contrary to conventional wisdom, network indeterminacy is unlikely to emerge in large competitive and mutualistic communities, while it is a common feature in top-down regulated food webs. Furthermore, we found that predictable and unpredictable perturbations can coexist in the same community and that indeterminate responses to environmental changes arise more frequently in networks where predator-prey relationships predominate than competitive and mutualistic ones. These findings highlight the importance of elucidating direct species relationships and analyzing them with an RMT perspective on two fronts: It aids in 1) determining whether the network's responses to environmental changes are ultimately indeterminate and 2) identifying the types of perturbations causing less predictable outcomes in a complex ecosystem. In addition, our framework should apply to the inverse problem of network identification, i.e., determining whether observed responses to sustained perturbations can reconstruct their proximate causalities, potentially impacting other fields such as microbial and medical sciences., Competing Interests: Competing interests statement:The author declares no competing interest.

Mercury (Hg) is a heterogeneously distributed toxicant affecting wildlife and human health. Yet, the spatial distribution of Hg remains poorly documented, especially in food webs, even though this knowledge is essential to assess large-scale risk of toxicity for the biota and human populations. Here, we used seabirds to assess, at an unprecedented population and geographic magnitude and high resolution, the spatial distribution of Hg in North Atlantic marine food webs. To this end, we combined tracking data of 837 seabirds from seven different species and 27 breeding colonies located across the North Atlantic and Atlantic Arctic together with Hg analyses in feathers representing individual seabird contamination based on their winter distribution. Our results highlight an east-west gradient in Hg concentrations with hot spots around southern Greenland and the east coast of Canada and a cold spot in the Barents and Kara Seas. We hypothesize that those gradients are influenced by eastern (Norwegian Atlantic Current and West Spitsbergen Current) and western (East Greenland Current) oceanic currents and melting of the Greenland Ice Sheet. By tracking spatial Hg contamination in marine ecosystems and through the identification of areas at risk of Hg toxicity, this study provides essential knowledge for international decisions about where the regulation of pollutants should be prioritized., Competing Interests: Competing interests statement:The authors declare no competing interest.

Consumers range from specialists that feed on few resources to generalists that feed on many. Generalism has the clear advantage of having more resources to exploit, but the costs that limit generalism are less clear. We explore two understudied costs of generalism in a generalist amoeba predator, Dictyostelium discoideum , feeding on naturally co-occurring bacterial prey. Both involve costs of combining prey that are suitable on their own. First, amoebas exhibit a reduction in growth rate when they switched to one species of prey bacteria from another compared to controls that experience only the second prey. The effect was consistent across all six tested species of bacteria. These switching costs typically disappear within a day, indicating adjustment to new prey bacteria. This suggests that these costs are physiological. Second, amoebas usually grow more slowly on mixtures of prey bacteria compared to the expectation based on their growth on single prey. There were clear mixing costs in three of the six tested prey mixtures, and none showed significant mixing benefits. These results support the idea that, although amoebas can consume a variety of prey, they must use partially different methods and thus must pay costs to handle multiple prey, either sequentially or simultaneously., Competing Interests: Competing interests statement:The authors declare no competing interest.

Cross-ecosystem subsidies are critical to ecosystem structure and function, especially in recipient ecosystems where they are the primary source of organic matter to the food web. Subsidies are indicative of processes connecting ecosystems and can couple ecological dynamics across system boundaries. However, the degree to which such flows can induce cross-ecosystem cascades of spatial synchrony, the tendency for system fluctuations to be correlated across locations, is not well understood. Synchrony has destabilizing effects on ecosystems, adding to the importance of understanding spatiotemporal patterns of synchrony transmission. In order to understand whether and how spatial synchrony cascades across the marine-terrestrial boundary via resource subsidies, we studied the relationship between giant kelp forests on rocky nearshore reefs and sandy beach ecosystems that receive resource subsidies in the form of kelp wrack (detritus). We found that synchrony cascades from rocky reefs to sandy beaches, with spatiotemporal patterns mediated by fluctuations in live kelp biomass, wave action, and beach width. Moreover, wrack deposition synchronized local abundances of shorebirds that move among beaches seeking to forage on wrack-associated invertebrates, demonstrating that synchrony due to subsidies propagates across trophic levels in the recipient ecosystem. Synchronizing resource subsidies likely play an underappreciated role in the spatiotemporal structure, functioning, and stability of ecosystems., Competing Interests: Competing interests statement:The authors declare no competing interest.

Competing Interests: Competing interests statement:The authors declare no competing interest.

Many predator species make regular excursions from near-surface waters to the twilight (200 to 1,000 m) and midnight (1,000 to 3,000 m) zones of the deep pelagic ocean. While the occurrence of significant vertical movements into the deep ocean has evolved independently across taxonomic groups, the functional role(s) and ecological significance of these movements remain poorly understood. Here, we integrate results from satellite tagging efforts with model predictions of deep prey layers in the North Atlantic Ocean to determine whether prey distributions are correlated with vertical habitat use across 12 species of predators. Using 3D movement data for 344 individuals who traversed nearly 1.5 million km of pelagic ocean in [Formula: see text]42,000 d, we found that nearly every tagged predator frequented the twilight zone and many made regular trips to the midnight zone. Using a predictive model, we found clear alignment of predator depth use with the expected location of deep pelagic prey for at least half of the predator species. We compared high-resolution predator data with shipboard acoustics and selected representative matches that highlight the opportunities and challenges in the analysis and synthesis of these data. While not all observed behavior was consistent with estimated prey availability at depth, our results suggest that deep pelagic biomass likely has high ecological value for a suite of commercially important predators in the open ocean. Careful consideration of the disruption to ecosystem services provided by pelagic food webs is needed before the potential costs and benefits of proceeding with extractive activities in the deep ocean can be evaluated.

Picophytoplankton populations [ Prochlorococcus , Synechococcus (SYN), and picoeukaryotes] are dominant primary producers in the open ocean and projected to become more important with climate change. Their fates can vary, however, with microbial food web complexities. In the California Current Ecosystem, picophytoplankton biomass and abundance peak in waters of intermediate productivity and decrease at higher production. Using experimental data from eight cruises crossing the pronounced CCE trophic gradient, we tested the hypothesis that these declines are driven by intensified grazing on heterotrophic bacteria (HBAC) passed to similarly sized picophytoplankton via shared predators. Results confirm previously observed distributions as well as significant increases in bacterial abundance, cell growth, and grazing mortality with primary production. Mortalities of picophytoplankton, however, diverge from the bacterial mortality trend such that relative grazing rates on SYN compared to HBAC decline by 12-fold between low and high productivity waters. The large shifts in mortality rate ratios for coexisting populations are not explained by size variability but rather suggest high selectivity of grazer assemblages or tightly coupled tradeoffs in microbial growth advantages and grazing vulnerabilities. These findings challenge the long-held view that protistan grazing mainly determines overall biomass of microbial communities while viruses uniquely regulate diversity by "killing the winners".

Ecologists have long sought to understand how diversity and structure mediate the stability of whole ecosystems. For high-diversity food webs, the interactions between species are typically represented using matrices with randomly chosen interaction strengths. Unfortunately, this procedure tends to produce ecological systems with no underlying equilibrium solution, and so ecological inferences from this approach may be biased by nonbiological outcomes. Using recent computationally efficient methodological advances from metabolic networks, we employ for the first time an inverse approach to diversity-stability research. We compare classical random interaction matrices of realistic food web topology (hereafter the classical model) to feasible, biologically constrained, webs produced using the inverse approach. We show that an energetically constrained feasible model yields a far higher proportion of stable high-diversity webs than the classical random matrix approach. When we examine the energetically constrained interaction strength distributions of these matrix models, we find that although these diverse webs have consistent negative self-regulation, they do not require strong self-regulation to persist. These energetically constrained diverse webs instead show an increasing preponderance of weak interactions that are known to increase local stability. Further examination shows that some of these weak interactions naturally appear to arise in the model food webs from a constraint-generated realistic generalist-specialist trade-off, whereby generalist predators have weaker interactions than more specialized species. Additionally, the inverse technique we present here has enormous promise for understanding the role of the biological structure behind stable high-diversity webs and for linking empirical data to the theory.

The idea that changing environmental conditions drive adaptive evolution is a pillar of evolutionary ecology. But, the opposite-that adaptive evolution alters ecological processes-has received far less attention yet is critical for eco-evolutionary dynamics. We assessed the ecological impact of divergent values in a key adaptive trait using 16 populations of the brown anole lizard ( Anolis sagrei ). Mirroring natural variation, we established islands with short- or long-limbed lizards at both low and high densities. We then monitored changes in lower trophic levels, finding that on islands with a high density of short-limbed lizards, web-spider densities decreased and plants grew more via an indirect positive effect, likely through an herbivore-mediated trophic cascade. Our experiment provides strong support for evolution-to-ecology connections in nature, likely closing an otherwise well-characterized eco-evolutionary feedback loop.

In most ecosystems, microbes are the dominant consumers, commandeering much of the heterotrophic biomass circulating through food webs. Characterizing functional diversity within the microbiome, therefore, is critical to understanding ecosystem functioning, particularly in an era of global biodiversity loss. Using isotopic fingerprinting, we investigated the trophic positions of a broad diversity of heterotrophic organisms. Specifically, we examined the naturally occurring stable isotopes of nitrogen (15N:14N) within amino acids extracted from proteobacteria, actinomycetes, ascomycetes, and basidiomycetes, as well as from vertebrate and invertebrate macrofauna (crustaceans, fish, insects, and mammals). Here, we report that patterns of intertrophic 15N-discrimination were remarkably similar among bacteria, fungi, and animals, which permitted unambiguous measurement of consumer trophic position, independent of phylogeny or ecosystem type. The observed similarities among bacterial, fungal, and animal consumers suggest that within a trophic hierarchy, microbiota are equivalent to, and can be interdigitated with, macrobiota. To further test the universality of this finding, we examined Neotropical fungus gardens, communities in which bacteria, fungi, and animals are entwined in an ancient, quadripartite symbiosis. We reveal that this symbiosis is a discrete fourlevel food chain, wherein bacteria function as the apex carnivores, animals and fungi are meso-consumers, and the sole herbivores are fungi. Together, our findings demonstrate that bacteria, fungi, and animals can be integrated within a food chain, effectively uniting the macro- and microbiome in food web ecology and facilitating greater inclusion of the microbiome in studies of functional diversity. [ABSTRACT FROM AUTHOR]

Sea otters ( Enhydra lutris ) and wolves ( Canis lupus ) are two apex predators with strong and cascading effects on ecosystem structure and function. After decades of recovery from near extirpation, their ranges now overlap, allowing sea otters and wolves to interact for the first time in the scientific record. We intensively studied wolves during 2015 to 2021 in an island system colonized by sea otters in the 2000s and by wolves in 2013. After wolf colonization, we quantified shifts in foraging behavior with DNA metabarcoding of 689 wolf scats and stable isotope analyses, both revealing a dietary switch from Sitka black-tailed deer ( Odocoileus hemionus ), the terrestrial in situ primary prey, to sea otters. Here we show an unexpected result of the reintroduction and restoration of sea otters, which became an abundant marine subsidy for wolves following population recovery. The availability of sea otters allowed wolves to persist and continue to reproduce, subsequently nearly eliminating deer. Genotypes from 390 wolf scats and telemetry data from 13 wolves confirmed island fidelity constituting one of the highest known wolf population densities and upending standardly accepted wolf density predictions based on ungulate abundance. Whereas marine subsidies in other systems are generally derived from lower trophic levels, here an apex nearshore predator became a key prey species and linked nearshore and terrestrial food webs in a recently deglaciated and rapidly changing ecosystem. These results underscore that species restoration may serve as an unanticipated nutrient pathway for recipient ecosystems even resulting in cross-boundary subsidy cascades.

Human exposure to monomethylmercury (CH 3 Hg), a potent neurotoxin, is principally through the consumption of seafood. The formation of CH 3 Hg and its bioaccumulation in marine food webs experience ongoing impacts of global climate warming and ocean biogeochemistry alterations. Employing a series of sensitivity experiments, here we explicitly consider the effects of climate change on marine mercury (Hg) cycling within a global ocean model in the hypothesized twenty-first century under the business-as-usual scenario. Even though the overall prediction is subjected to significant uncertainty, we identify several important climate change impact pathways. Elevated seawater temperature exacerbates elemental Hg (Hg 0 ) evasion, while decreased surface wind speed reduces air-sea exchange rates. The reduced export of particulate organic carbon shrinks the pool of potentially bioavailable divalent Hg (Hg II ) that can be methylated in the subsurface ocean, where shallower remineralization depth associated with lower productivity causes impairment of methylation activity. We also simulate an increase in CH 3 Hg photodemethylation potential caused by increased incident shortwave radiation and less attenuation by decreased sea ice and chlorophyll. The model suggests that these impacts can also be propagated to the CH 3 Hg concentration in the base of the marine food web. Our results offer insight into synergisms/antagonisms in the marine Hg cycling among different climate change stressors.

Viruses impact host cells and have indirect effects on ecosystem processes. Plankton such as ciliates can reduce the abundance of virions in water, but whether virus consumption translates into demographic consequences for the grazers is unknown. Here, we show that small protists not only can consume viruses they also can grow and divide given only viruses to eat. Moreover, the ciliate Halteria sp. foraging on chloroviruses displays dynamics and interaction parameters that are similar to other microbial trophic interactions. These results suggest that the effect of viruses on ecosystems extends beyond (and in contrast to) the viral shunt by redirecting energy up food chains.

The marine pelagic compartment spans numerous trophic levels and consists of numerous reticulate connections between species from primary producers to iconic apex predators, while the benthic compartment is perceived to be simpler in structure and comprised of only low trophic level species. Here, we challenge this paradigm by illustrating that the benthic compartment is home to a subweb of similar structure and complexity to that of the pelagic realm, including the benthic equivalent to iconic polar bears: megafaunal-predatory sea stars.

Phosphorus (P) is a key nutrient limiting bacterial growth and primary production in the oceans. Unsurprisingly, marine microbes have evolved sophisticated strategies to adapt to P limitation, one of which involves the remodeling of membrane lipids by replacing phospholipids with non-P-containing surrogate lipids. This strategy is adopted by both cosmopolitan marine phytoplankton and heterotrophic bacteria and serves to reduce the cellular P quota. However, little, if anything, is known of the biological consequences of lipid remodeling. Here, using the marine bacterium Phaeobacter sp. MED193 and the ciliate Uronema marinum as a model, we sought to assess the effect of remodeling on bacteria-protist interactions. We discovered an important trade-off between either escape from ingestion or resistance to digestion. Thus, Phaeobacter grown under P-replete conditions was readily ingested by Uronema , but not easily digested, supporting only limited predator growth. In contrast, following membrane lipid remodeling in response to P depletion, Phaeobacter was less likely to be captured by Uronema , thanks to the reduced expression of mannosylated glycoconjugates. However, once ingested, membrane-remodeled cells were unable to prevent phagosome acidification, became more susceptible to digestion, and, as such, allowed rapid growth of the ciliate predator. This trade-off between adapting to a P-limited environment and susceptibility to protist grazing suggests the more efficient removal of low-P prey that potentially has important implications for the functioning of the marine microbial food web in terms of trophic energy transfer and nutrient export efficiency.

Distribution of Earth's biomes is structured by the match between climate and plant traits, which in turn shape associated communities and ecosystem processes and services. However, that climate-trait match can be disrupted by historical events, with lasting ecosystem impacts. As Earth's environment changes faster than at any time in human history, critical questions are whether and how organismal traits and ecosystems can adjust to altered conditions. We quantified the relative importance of current environmental forcing versus evolutionary history in shaping the growth form (stature and biomass) and associated community of eelgrass ( Zostera marina ), a widespread foundation plant of marine ecosystems along Northern Hemisphere coastlines, which experienced major shifts in distribution and genetic composition during the Pleistocene. We found that eelgrass stature and biomass retain a legacy of the Pleistocene colonization of the Atlantic from the ancestral Pacific range and of more recent within-basin bottlenecks and genetic differentiation. This evolutionary legacy in turn influences the biomass of associated algae and invertebrates that fuel coastal food webs, with effects comparable to or stronger than effects of current environmental forcing. Such historical lags in phenotypic acclimatization may constrain ecosystem adjustments to rapid anthropogenic climate change, thus altering predictions about the future functioning of ecosystems.

Many studies have aimed to understand food webs by investigating components such as trophic links (one consumer taxon eats one resource taxon), tritrophic interactions (one consumer eats an intermediate taxon, which eats a resource), or longer chains of links. We show here that none of these components (links, tritrophic interactions, and longer chains), individually or as an ensemble, accounts fully for the properties of the next higher level of organization. As a cell is more than its molecules, as an organ is more than its cells, and as an organism is more than its organs, in a food web, new structure emerges at every organizational level up to and including the whole web. We demonstrate the emergence of properties at progressively higher levels of structure by using all of the directly observed, appropriately organized, publicly available food web datasets with relatively complete trophic link data and with average body mass and population density data for each taxon. There are only three such webs, those of Tuesday Lake, Michigan, in 1984 and 1986, and Ythan Estuary, Scotland. We make the data freely available online with this report. Differences in web patterns between Tuesday Lake and Ythan Estuary, and similarities of Tuesday Lake in 1984 and 1986 despite 50% turnover of species, suggest that the patterns we describe respond to major differences between ecosystem types. [ABSTRACT FROM AUTHOR]

Parasitic infections are common, but how they shape ecosystem-level processes is understudied. Using a mathematical model and meta-analysis, we explored the potential for helminth parasites to trigger trophic cascades through lethal and sublethal effects imposed on herbivorous ruminant hosts after infection. First, using the model, we linked negative effects of parasitic infection on host survival, fecundity, and feeding rate to host and producer biomass. Our model, parameterized with data from a well-documented producer–caribou–helminth system, reveals that even moderate impacts of parasites on host survival, fecundity, or feeding rate can have cascading effects on ruminant host and producer biomass. Second, using meta-analysis, we investigated the links between helminth infections and traits of free-living ruminant hosts in nature. We found that helminth infections tend to exert negative but sublethal effects on ruminant hosts. Specifically, infection significantly reduces host feeding rates, body mass, and body condition but has weak and highly variable effects on survival and fecundity. Together, these findings suggest that while helminth parasites can trigger trophic cascades through multiple mechanisms, overlooked sublethal effects on nonreproductive traits likely dominate their impacts on ecosystems. In particular, by reducing ruminant herbivory, pervasive helminth infections may contribute to a greener world.

Ecological communities are constantly exposed to multiple natural and anthropogenic disturbances. Multivariate composition (if recovered) has been found to need significantly more time to be regained after pulsed disturbance compared to univariate diversity metrics and functional endpoints. However, the mechanisms driving the different recovery times of communities to single and multiple disturbances remain unexplored. Here, we apply quantitative ecological network analyses to try to elucidate the mechanisms driving long-term community-composition dissimilarity and late-stage disturbance interactions at the community level. For this, we evaluate the effects of two pesticides, nutrient enrichment, and their interactions in outdoor mesocosms containing a complex freshwater community. We found changes in interactions strength to be strongly related to compositional changes and identified postdisturbance interaction-strength rewiring to be responsible for most of the observed compositional changes. Additionally, we found pesticide interactions to be significant in the long term only when both interaction strength and food-web architecture are reshaped by the disturbances. We suggest that quantitative network analysis has the potential to unveil ecological processes that prevent long-term community recovery.

The factors that determine why ecosystems exhibit abrupt shifts in state are of paramount importance for management, conservation, and restoration efforts. Kelp forests are emblematic of such abruptly shifting ecosystems, transitioning from kelp-dominated to urchin-dominated states around the world with increasing frequency, yet the underlying processes and mechanisms that control their dynamics remain unclear. Here, we analyze four decades of data from biannual monitoring around San Nicolas Island, CA, to show that substrate complexity controls both the number of possible (alternative) states and the velocity with which shifts between states occur. The superposition of community dynamics with reconstructions of system stability landscapes reveals that shifts between alternative states at low-complexity sites reflect abrupt, high-velocity events initiated by pulse perturbations that rapidly propel species across dynamically unstable state-space. In contrast, high-complexity sites exhibit a single state of resilient kelp-urchin coexistence. Our analyses suggest that substrate complexity influences both top-down and bottom-up regulatory processes in kelp forests, highlight its influence on kelp-forest stability at both large (island-wide) and small (<10 m) spatial scales, and could be valuable for holistic kelp-forest management., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

Intricate networks of single-celled eukaryotes (protists) dominate carbon flow in the ocean. Their growth, demise, and interactions with other microorganisms drive the fluxes of biogeochemical elements through marine ecosystems. Mixotrophic protists are capable of both photosynthesis and ingestion of prey and are dominant components of open-ocean planktonic communities. Yet the role of mixotrophs in elemental cycling is obscured by their capacity to act as primary producers or heterotrophic consumers depending on factors that remain largely uncharacterized. Here, we develop and apply a machine learning model that predicts the in situ trophic mode of aquatic protists based on their patterns of gene expression. This approach leverages a public collection of protist transcriptomes as a training set to identify a subset of gene families whose transcriptional profiles predict trophic mode. We applied our model to nearly 100 metatranscriptomes obtained during two oceanographic cruises in the North Pacific and found community-level and population-specific evidence that abundant open-ocean mixotrophic populations shift their predominant mode of nutrient and carbon acquisition in response to natural gradients in nutrient supply and sea surface temperature. Metatranscriptomic data from ship-board incubation experiments revealed that abundant mixotrophic prymnesiophytes from the oligotrophic North Pacific subtropical gyre rapidly remodeled their transcriptome to enhance photosynthesis when supplied with limiting nutrients. Coupling this approach with experiments designed to reveal the mechanisms driving mixotroph physiology provides an avenue toward understanding the ecology of mixotrophy in the natural environment., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

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

Global climate warming is causing the loss of freshwater ice around the Northern Hemisphere. Although the timing and duration of ice covers are known to regulate ecological processes in seasonally ice-covered ecosystems, the consequences of shortening winters for freshwater biota are poorly understood owing to the scarcity of under-ice research. Here, we present one of the first in-lake experiments to postpone ice-cover onset (by ≤21 d), thereby extending light availability (by ≤40 d) in early winter, and explicitly demonstrate cascading effects on pelagic food web processes and phenologies. Delaying ice-on elicited a sequence of events from winter to spring: 1) relatively greater densities of algal resources and primary consumers in early winter; 2) an enhanced prevalence of winter-active (overwintering) consumers throughout the ice-covered period, associated with augmented storage of high-quality fats likely due to a longer access to algal resources in early winter; and 3) an altered trophic structure after ice-off, with greater initial springtime densities of overwintering consumers driving stronger, earlier top-down regulation, effectively reducing the spring algal bloom. Increasingly later ice onset may thus promote consumer overwintering, which can confer a competitive advantage on taxa capable of surviving winters upon ice-off; a process that may diminish spring food availability for other consumers, potentially disrupting trophic linkages and energy flow pathways over the subsequent open-water season. In considering a future with warmer winters, these results provide empirical evidence that may help anticipate phenological responses to freshwater ice loss and, more broadly, constitute a case of climate-induced cross-seasonal cascade on realized food web processes., Competing Interests: The authors declare no competing interest.

Species invasions can have substantial impacts on native species and ecosystems, with important consequences for biodiversity. How these disturbances drive changes in the trophic structure of native food webs through time is poorly understood. Here, we quantify trophic disruption in freshwater food webs to invasion by an apex fish predator, lake trout, using an extensive stable isotope dataset across a natural gradient of uninvaded and invaded lakes in the northern Rocky Mountains, USA. Lake trout invasion increased fish diet variability (trophic dispersion), displaced native fishes from their reference diets (trophic displacement), and reorganized macroinvertebrate communities, indicating strong food web disruption. Trophic dispersion was greatest 25 to 50 y after colonization and dissipated as food webs stabilized in later stages of invasion (>50 y). For the native apex predator, bull trout, trophic dispersion preceded trophic displacement, leading to their functional loss in late-invasion food webs. Our results demonstrate how invasive species progressively disrupt native food webs via trophic dispersion and displacement, ultimately yielding biological communities strongly divergent from those in uninvaded ecosystems., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Ecological interactions uphold ecosystem structure and functioning. However, as species richness increases, the number of possible interactions rises exponentially. More than 6,000 species of coral reef fishes exist across the world's tropical oceans, resulting in an almost innumerable array of possible trophic interactions. Distilling general patterns in these interactions across different bioregions stands to improve our understanding of the processes that govern coral reef functioning. Here, we show that across bioregions, tropical coral reef food webs exhibit a remarkable congruence in their trophic interactions. Specifically, by compiling and investigating the structure of six coral reef food webs across distinct bioregions, we show that when accounting for consumer size and resource availability, these food webs share more trophic interactions than expected by chance. In addition, coral reef food webs are dominated by dietary specialists, which makes trophic pathways vulnerable to biodiversity loss. Prey partitioning among these specialists is geographically consistent, and this pattern intensifies when weak interactions are disregarded. Our results suggest that energy flows through coral reef communities along broadly comparable trophic pathways. Yet, these critical pathways are maintained by species with narrow, specialized diets, which threatens the existence of coral reef functioning in the face of biodiversity loss., Competing Interests: The authors declare no competing interest.

Bluefin tuna (BFT), highly prized among consumers, accumulate high levels of mercury (Hg) as neurotoxic methylmercury (MeHg). However, how Hg bioaccumulation varies among globally distributed BFT populations is not understood. Here, we show mercury accumulation rates (MARs) in BFT are highest in the Mediterranean Sea and decrease as North Pacific Ocean > Indian Ocean > North Atlantic Ocean. Moreover, MARs increase in proportion to the concentrations of MeHg in regional seawater and zooplankton, linking MeHg accumulation in BFT to MeHg bioavailability at the base of each subbasin's food web. Observed global patterns correspond to levels of Hg in each ocean subbasin; the Mediterranean, North Pacific, and Indian Oceans are subject to geogenic enrichment and anthropogenic contamination, while the North Atlantic Ocean is less so. MAR in BFT as a global pollution index reflects natural and human sources and global thermohaline circulation., Competing Interests: The authors declare no competing interest.

Ecologists are still puzzled by the diverse population dynamics of herbivorous small mammals that range from high-amplitude, multiannual cycles to stable dynamics. Theory predicts that this diversity results from combinations of climatic seasonality, weather stochasticity, and density-dependent food web interactions. The almost ubiquitous 3- to 5-y cycles in boreal and arctic climates may theoretically result from bottom-up (plant-herbivore) and top-down (predator-prey) interactions. Assessing, empirically, the roles of such interactions and how they are influenced by environmental stochasticity has been hampered by food web complexity. Here, we take advantage of a uniquely simple High Arctic food web, which allowed us to analyze the dynamics of a graminivorous vole population not subjected to top-down regulation. This population exhibited high-amplitude, noncyclic fluctuations-partly driven by weather stochasticity. However, the predominant driver of the dynamics was overcompensatory density dependence in winter that caused the population to frequently crash. Model simulations showed that the seasonal pattern of density dependence would yield regular 2-y cycles in the absence of stochasticity. While such short cycles have not yet been observed in mammals, they are theoretically plausible if graminivorous vole populations are deterministically bottom-up regulated. When incorporating weather stochasticity in the model simulations, cyclicity became disrupted and the amplitude was increased-akin to the observed dynamics. Our findings contrast with the 3- to 5-y population cycles that are typical of graminivorous small mammals in more complex food webs, suggesting that top-down regulation is normally an important component of such dynamics., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Mesopredator release theory suggests that dominant predators suppress subordinate carnivores and ultimately shape community dynamics, but the assumption that subordinate species are only negatively affected ignores the possibility of facilitation through scavenging. We examined the interplay within a carnivore community consisting of cougars, coyotes, black bears, and bobcats using contemporaneous Global Positioning System telemetry data from 51 individuals; diet analysis from 972 DNA-metabarcoded scats; and data from 128 physical investigations of cougar kill sites, 28 of which were monitored with remote cameras. Resource provisioning from competitively dominant cougars to coyotes through scavenging was so prolific as to be an overwhelming determinant of coyote behavior, space use, and resource acquisition. This was evident via the strong attraction of coyotes to cougar kill sites, frequent scavenging of cougar-killed prey, and coyote diets that nearly matched cougars in the magnitude of ungulate consumption. Yet coyotes were often killed by cougars and used space to minimize encounters, complicating the fitness benefits gained from scavenging. We estimated that 23% (95% CI: 8 to 55%) of the coyote population in our study area was killed by cougars annually, suggesting that coyote interactions with cougars are a complex behavioral game of risk and reward. In contrast, we found no indication that bobcat space use or diet was influenced by cougars. Black bears avoided cougars, but there was no evidence of attraction to cougar kill sites and much lower levels of ungulate consumption and carcass visitation than for coyotes. Interspecific interactions among carnivores are multifaceted, encompassing both suppression and facilitation., Competing Interests: The authors declare no competing interest.

Metapopulation capacity provides an analytic tool to quantify the impact of landscape configuration on metapopulation persistence, which has proven powerful in biological conservation. Yet surprisingly few efforts have been made to apply this approach to multispecies systems. Here, we extend metapopulation capacity theory to predict the persistence of trophically interacting species. Our results demonstrate that metapopulation capacity could be used to predict the persistence of trophic systems such as prey-predator pairs and food chains in fragmented landscapes. In particular, we derive explicit predictions for food chain length as a function of metapopulation capacity, top-down control, and population dynamical parameters. Under certain assumptions, we show that the fraction of empty patches for the basal species provides a useful indicator to predict the length of food chains that a fragmented landscape can support and confirm this prediction for a host-parasitoid interaction. We further show that the impact of habitat changes on biodiversity can be predicted from changes in metapopulation capacity or approximately by changes in the fraction of empty patches. Our study provides an important step toward a spatially explicit theory of trophic metacommunities and a useful tool for predicting their responses to habitat changes., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

While the important role of animal-mediated interactions in the top-down restructuring of plant communities is well documented, less is known of their ensuing repercussions at higher trophic levels. We demonstrate how typically decoupled ecological interactions may become intertwined such that the impact of an insect pest on forest structure and composition alters predator-prey interactions among large mammals. Specifically, we show how irruptions in a common, cyclic insect pest of the boreal forest, the spruce budworm ( Choristoneura fumiferana ), modulated an indirect trophic interaction by initiating a flush in deciduous vegetation that benefited moose ( Alces alces ), in turn strengthening apparent competition between moose and threatened boreal caribou ( Rangifer tarandus caribou ) via wolf ( Canis lupus ) predation. Critically, predation on caribou postoutbreak was exacerbated by human activity (salvage logging). We believe our observations of significant, large-scale reverberating consumer-producer-consumer interactions are likely to be common in nature., Competing Interests: The authors declare no competing interest.

Microbial interactions in aquatic environments profoundly affect global biogeochemical cycles, but the role of microparasites has been largely overlooked. Using a model pathosystem, we studied hitherto cryptic interactions between microparasitic fungi (chytrid Rhizophydiales ), their diatom host Asterionella, and cell-associated and free-living bacteria. We analyzed the effect of fungal infections on microbial abundances, bacterial taxonomy, cell-to-cell carbon transfer, and cell-specific nitrate-based growth using microscopy (e.g., fluorescence in situ hybridization), 16S rRNA gene amplicon sequencing, and secondary ion mass spectrometry. Bacterial abundances were 2 to 4 times higher on individual fungal-infected diatoms compared to healthy diatoms, particularly involving Burkholderiales. Furthermore, taxonomic compositions of both diatom-associated and free-living bacteria were significantly different between noninfected and fungal-infected cocultures. The fungal microparasite, including diatom-associated sporangia and free-swimming zoospores, derived ∼100% of their carbon content from the diatom. By comparison, transfer efficiencies of photosynthetic carbon were lower to diatom-associated bacteria (67 to 98%), with a high cell-to-cell variability, and even lower to free-living bacteria (32%). Likewise, nitrate-based growth for the diatom and fungi was synchronized and faster than for diatom-associated and free-living bacteria. In a natural lacustrine system, where infection prevalence reached 54%, we calculated that 20% of the total diatom-derived photosynthetic carbon was shunted to the parasitic fungi, which can be grazed by zooplankton, thereby accelerating carbon transfer to higher trophic levels and bypassing the microbial loop. The herein termed "fungal shunt" can thus significantly modify the fate of photosynthetic carbon and the nature of phytoplankton-bacteria interactions, with implications for diverse pelagic food webs and global biogeochemical cycles., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Natural ecological communities are diverse, complex, and often surprisingly stable, but the mechanisms underlying their stability remain a theoretical enigma. Interactions such as competition and predation presumably structure communities, yet theory predicts that complex communities are stable only when species growth rates are mostly limited by intraspecific self-regulation rather than by interactions with resources, competitors, and predators. Current theory, however, considers only the network topology of population-level interactions between species and ignores within-population differences, such as between juvenile and adult individuals. Here, using model simulations and analysis, I show that including commonly observed differences in vulnerability to predation and foraging efficiency between juvenile and adult individuals results in up to 10 times larger, more complex communities than observed in simulations without population stage structure. These diverse communities are stable or fluctuate with limited amplitude, although in the model only a single basal species is self-regulated, and the population-level interaction network is highly connected. Analysis of the species interaction matrix predicts the simulated communities to be unstable but for the interaction with the population-structure subsystem, which completely cancels out these instabilities through dynamic changes in population stage structure. Common differences between juveniles and adults and fluctuations in their relative abundance may hence have a decisive influence on the stability of complex natural communities and their vulnerability when environmental conditions change. To explain community persistence, it may not be sufficient to consider only the network of interactions between the constituting species., Competing Interests: The author declares no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Population-level scaling in ecological systems arises from individual growth and death with competitive constraints. We build on a minimal dynamical model of metabolic growth where the tension between individual growth and mortality determines population size distribution. We then separately include resource competition based on shared capture area. By varying rates of growth, death, and competitive attrition, we connect regular and random spatial patterns across sessile organisms from forests to ants, termites, and fairy circles. Then, we consider transient temporal dynamics in the context of asymmetric competition, such as canopy shading or large colony dominance, whose effects primarily weaken the smaller of two competitors. When such competition couples slow timescales of growth to fast competitive death, it generates population shocks and demographic oscillations similar to those observed in forest data. Our minimal quantitative theory unifies spatiotemporal patterns across sessile organisms through local competition mediated by the laws of metabolic growth, which in turn, are the result of long-term evolutionary dynamics., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Consumer and predator foraging behavior can impart profound trait-mediated constraints on community regulation that scale up to influence the structure and stability of ecosystems. Here, we demonstrate how the behavioral response of an apex predator to changes in prey behavior and condition can dramatically alter the role and relative contribution of top-down forcing, depending on the spatial organization of ecosystem states. In 2014, a rapid and dramatic decline in the abundance of a mesopredator ( Pycnopodia helianthoides ) and primary producer ( Macrocystis pyrifera ) coincided with a fundamental change in purple sea urchin ( Strongylocentrotus purpuratus ) foraging behavior and condition, resulting in a spatial mosaic of kelp forests interspersed with patches of sea urchin barrens. We show that this mosaic of adjacent alternative ecosystem states led to an increase in the number of sea otters ( Enhydra lutris nereis ) specializing on urchin prey, a population-level increase in urchin consumption, and an increase in sea otter survivorship. We further show that the spatial distribution of sea otter foraging efforts for urchin prey was not directly linked to high prey density but rather was predicted by the distribution of energetically profitable prey. Therefore, we infer that spatially explicit sea otter foraging enhances the resistance of remnant forests to overgrazing but does not directly contribute to the resilience (recovery) of forests. These results highlight the role of consumer and predator trait-mediated responses to resource mosaics that are common throughout natural ecosystems and enhance understanding of reciprocal feedbacks between top-down and bottom-up forcing on the regional stability of ecosystems., Competing Interests: The authors declare no competing interest.

The first Caribbean settlers were Amerindians from South America. Great Abaco and Grand Bahama, the final islands colonized in the northernmost Bahamas, were inhabited by the Lucayans when Europeans arrived. The timing of Lucayan arrival in the northern Bahamas has been uncertain because direct archaeological evidence is limited. We document Lucayan arrival on Great Abaco Island through a detailed record of vegetation, fire, and landscape dynamics based on proxy data from Blackwood Sinkhole. From about 3,000 to 1,000 y ago, forests dominated by hardwoods and palms were resilient to the effects of hurricanes and cooling sea surface temperatures. The arrival of Lucayans by about 830 CE (2σ range: 720 to 920 CE) is demarcated by increased burning and followed by landscape disturbance and a time-transgressive shift from hardwoods and palms to the modern pine forest. Considering that Lucayan settlements in the southern Bahamian archipelago are dated to about 750 CE (2σ range: 600 to 900 CE), these results demonstrate that Lucayans spread rapidly through the archipelago in less than 100 y. Although precontact landscapes would have been influenced by storms and climatic trends, the most pronounced changes follow more directly from landscape burning and ecosystem shifts after Lucayan arrival. The pine forests of Abaco declined substantially between 1500 and 1670 CE, a period of increased regional hurricane activity, coupled with fires on an already human-impacted landscape. Any future intensification of hurricane activity in the tropical North Atlantic Ocean threatens the sustainability of modern pine forests in the northern Bahamas., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

One of the most prominent features of life on Earth is the uneven number of species across large spatial scales. Despite being inherently linked to energetic constraints, these gradients in species richness distribution have rarely been examined from a trophic perspective. Here we dissect the global diversity of over 3,600 coral reef fishes to reveal patterns across major trophic groups. By analyzing multiple nested spatial scales, we show that planktivores contribute disproportionally to the formation of the Indo-Australian Archipelago (IAA) marine biodiversity hotspot. Besides being "hotter" at the hotspot, planktivorous fishes display the steepest decline in species numbers with distance from the IAA when compared to other trophic groups. Surprisingly, we did not detect differences in diversification, transition, and dispersal rates in extant species phylogenies that would explain this remarkable gradient in planktivorous fish richness. Thus, we identify two potential complementary drivers for this pattern. First, exceptional levels of partitioning among planktivorous coral reef fishes were driven by temporally stable oceanographic conditions and abundant planktonic resources in the IAA. Second, extinctions of planktivores outside the IAA have been particularly pronounced during Quaternary climate fluctuations. Overall, our results highlight trophic ecology as an important component of global species richness gradients., Competing Interests: The authors declare no competing interest.

Predator specialization has often been considered an evolutionary "dead end" due to the constraints associated with the evolution of morphological and functional optimizations throughout the organism. However, in some predators, these changes are localized in separate structures dedicated to prey capture. One of the most extreme cases of this modularity can be observed in siphonophores, a clade of pelagic colonial cnidarians that use tentilla (tentacle side branches armed with nematocysts) exclusively for prey capture. Here we study how siphonophore specialists and generalists evolve, and what morphological changes are associated with these transitions. To answer these questions, we: a) Measured 29 morphological characters of tentacles from 45 siphonophore species, b) mapped these data to a phylogenetic tree, and c) analyzed the evolutionary associations between morphological characters and prey-type data from the literature. Instead of a dead end, we found that siphonophore specialists can evolve into generalists, and that specialists on one prey type have directly evolved into specialists on other prey types. Our results show that siphonophore tentillum morphology has strong evolutionary associations with prey type, and suggest that shifts between prey types are linked to shifts in the morphology, mode of evolution, and evolutionary correlations of tentilla and their nematocysts. The evolutionary history of siphonophore specialization helps build a broader perspective on predatory niche diversification via morphological innovation and evolution. These findings contribute to understanding how specialization and morphological evolution have shaped present-day food webs., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

The evolution of taste perception is usually associated with the ecology and dietary changes of organisms. However, the association between feeding ecology and taste receptor evolution is unclear in some lineages of vertebrate animals. One example is the sweet taste receptor gene Tas1r2 Previous analysis of partial sequences has revealed that Tas1r2 has undergone equally strong purifying selection between insectivorous and frugivorous bats. To test whether the sweet taste function is also important in bats with contrasting diets, we examined the complete coding sequences of both sweet taste receptor genes ( Tas1r2 and Tas1r3 ) in 34 representative bat species. Although these two genes are highly conserved between frugivorous and insectivorous bats at the sequence level, our behavioral experiments revealed that an insectivorous bat ( Myotis ricketti ) showed no preference for natural sugars, whereas the frugivorous species ( Rousettus leschenaultii ) showed strong preferences for sucrose and fructose. Furthermore, while both sweet taste receptor genes are expressed in the taste tissue of insectivorous and frugivorous bats, our cell-based assays revealed striking functional divergence: the sweet taste receptors of frugivorous bats are able to respond to natural sugars whereas those of insectivorous bats are not, which is consistent with the behavioral preference tests, suggesting that functional evolution of sweet taste receptors is closely related to diet. This comprehensive study suggests that using sequence conservation alone could be misleading in inferring protein and physiological function and highlights the power of combining behavioral experiments, expression analysis, and functional assays in molecular evolutionary studies., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

An evolutionarily stable strategy (ESS) is an evolutionary strategy that, if adapted by a population, cannot be invaded by any deviating (mutant) strategy. The concept of ESS has been extensively studied and widely applied in ecology and evolutionary biology [M. Smith, On Evolution (1972)] but typically on the assumption that the system is ecologically stable. With reference to a Rosenzweig-MacArthur predator-prey model [M. Rosenzweig, R. MacArthur, Am. Nat. 97, 209-223 (1963)], we derive the mathematical conditions for the existence of an ESS when the ecological dynamics have asymptotically stable limit points as well as limit cycles. By extending the framework of Reed and Stenseth [J. Reed, N. C. Stenseth, J. Theoret. Biol. 108, 491-508 (1984)], we find that ESSs occur at values of the evolutionary strategies that are local optima of certain functions of the model parameters. These functions are identified and shown to have a similar form for both stable and fluctuating populations. We illustrate these results with a concrete example., Competing Interests: Competing interest statement: R.D.H. and N.C.S. are coauthors on a 2019 review article., (Copyright © 2021 the Author(s). Published by PNAS.)