More stories

  • in

    Redox-informed models of global biogeochemical cycles

    Microorganisms drive biogeochemical cycling in the earth system1 (Fig. 1). Photoautotrophic microorganisms are responsible for about half of CO2 fixation and O2 production on earth, and heterotrophic microorganisms are responsible for much of the return reaction: the oxidation of organic matter back into CO2. The temporal and spatial separation of photoautotrophy and heterotrophy in the global environment drives the biological sequestration of carbon, the reduction of atmospheric CO2, and the maintenance of elevated atmospheric and oceanic O22,3,4,5. Chemoautotrophic microorganisms also fix CO2 and, together with anaerobic heterotrophic metabolisms, carry out diverse chemical transformations including the fluxes of nitrogen to and from biologically available states and the formation of the potent greenhouse gas nitrous oxide (N2O)6,7. Since these transformations respond to, and feedback on, changes in climate (Fig. 1), estimating microbial activity accurately at global scales is important for climate science.
    Fig. 1: Key microbially driven redox transformations that mediate the atmospheric fluxes of climatically relevant gases.

    Radiatively active gases are notated with red type. The processes in black type are represented in some way (though not necessarily with electron balancing) in both the marine and terrestrial biospheres in earth system models within the Coupled Model Intercomparison Project (land: NCAR Community Earth System Model103; ocean: GFDL COBALTv2104), which are used for projections of climate change in reports by the Intergovernmental Panel on Climate Change. Processes in green type are represented in only the terrestrial model. Current models do not yet include other relevant reactions, some of which are represented in gray type, such as anaerobic ammonia oxidation (anammox), the marine production and consumption of methane, the redox cycling of iron, manganese, and other metals, and the methane-relevant redox chemistry of phosphorus105. COBALTv2 does account for sulfate reduction in marine sediments, but sulfate is not represented. Image courtesy of NASA.

    Full size image

    However, understanding and projecting the impacts of microbial processes are limited in part due to oversimplified representation in earth system models. For example, in marine biogeochemical models, much attention is given to the complex impacts of phytoplankton—the photoautotrophic microorganisms responsible for primary production—and their small zooplankton predators8,9,10. The bacterial and archaeal activities responsible for other critical aspects of biogeochemical cycling in the land and ocean—remineralization, denitrification, nitrogen fixation, methanogenesis, etc.—are often crudely parameterized10. Such models have limited prognostic capability. For example, models typically prescribe the ecological niche of a given metabolism with imposed, empirically determined parameters that are site- or organism-specific. These parameterizations may or may not apply to other environments, including past and future ecosystems.
    These simplistic approaches have been largely necessary due to the difficulties of characterizing the taxonomy and metabolic capabilities of natural microbial communities. However, the rapid expansion of genetic sequencing capabilities has enabled a clearer view of microbial biogeography and activity in the environment. In consequence, computational biogeochemistry is opening up the black box of remineralization and other microbially mediated processes in marine and terrestrial environments11,12,13,14,15,16,17,18.
    As we expand models to include the full metabolic potential of microorganisms, how can we organize and reduce the complexity of the descriptions of metabolic diversity? Non-photosynthetic organisms oxidize chemical species for energy, and thus their respiration is biogeochemically significant19. Here, we explain how the key reduction-oxidation (redox) reactions that supply energy for metabolisms can provide an additional organizing principle for explicit descriptions of microbial populations in ecosystem models. This redox basis can be exploited to quantitatively resolve chemical transformations in terms of assimilatory and respiratory fluxes. While not yet incorporated into earth system models, this view has been advocated for such applications20, and has been embraced and employed in the field of environmental biotechnology, such as in the interpretation and modeling of wastewater bioreactors21. Just as models of ocean and atmospheric circulation are constrained by conservation of energy and potential vorticity, complementing mass balance with powerful redox and energetic constraints enables self-consistent descriptions of diverse microbial metabolisms.
    This approach aims to advance ecological modeling beyond species-specific descriptions to those that matter for biogeochemical function, in line with trait-based modeling approaches9. In analogy to the use of redox chemistry, trait-based functional type models of phytoplankton have used cell size as an organizing principle for understanding phytoplankton biogeography, biodiversity, and impact on biogeochemistry9,22,23. These types of theoretical constraints allow for the inclusion of more functional types without introducing as many degrees of freedom as would be necessary if each were empirically described. The guiding perspective is that organizing complex biological behavior by its underlying chemical and physical constraints gives more universally applicable descriptions of large-scale biogeochemistry.
    When incorporating a redox-balanced approach into ecosystem models, microbial function emerges from underlying chemistry as a consequence of interactions between populations modeled as metabolic functional types and their environment. Resulting theoretically grounded ecosystem models independently simulate microbial growth, respiration, and abundances in ways that we can compare with observations such as sequencing datasets. Thus, sequencing datasets are used as critical tests for the models, as external constraints rather than as input to the models, allowing for an iterative relationship between theory, observations, and models.
    In contrast with empirically informed models, this approach involves constructing a model of microbial activity theoretically, and then comparing the results with the observations in order to gain an understanding of the system. The goal is to understand why biology functions as it does, in addition to anticipating global impacts. From a first-principles biogeochemical perspective with respect to physical and chemical forcing, genes are an intermediate step between forcing and function, with the detailed complexity of biological reality following the underlying chemical and physical constraints (analogous to the form follows function principle of architect Louis Sullivan). This does not equate to thinking that biology (or genetic information) does not matter or can be replaced entirely in models by physics and chemistry. Rather, we want to fundamentally understand biological activity as an integrated part of an ecosystem, and physics and chemistry become tools for doing so.
    Here, we outline the basis for using redox chemistry as an organizing principle and its translation into quantitative descriptions of microbial activity that are simple enough for global earth system models. We then discuss the benefits of this approach in the context of their implications for improved understanding and projections of global change impacts. Finally, we discuss limitations and possible future developments.
    Predicting microbial activity
    From one perspective, microbial communities are characterized by interactions at the micro-scale: gene expression, enzymatic capabilities, metabolites, species-specific interdependencies, etc., as well as the physical and chemical environment surrounding small cells24,25,26,27,28. The information from sequencing in particular has allowed for a huge expansion of insight into the detailed in situ activity of uncultivated species. When investigating global-scale impacts, how do we decide which of these details may be bypassed for simplicity? Or, if this simplification is impossible, must we incrementally construct a microbial ecosystem model that incorporates all known micro-scale detail?
    Another way forward arises from a macro-scale perspective, which examines how ecosystem function relates to the chemical potential utilized by organisms for energy20,29,30. For example, it is well known that microbial communities in sediments and anoxic zones organize according to the redox tower – the ranking of half-reactions by electrochemical potential14,31,32. Furthermore, respiration by living organisms increases the entropy of the environment by dissipating concentrated sources of chemical energy in accordance with the Second Law of Thermodynamics30,33,34.
    This perspective suggests that chemical potential can be used to predict the activity of microbial communities and their biogeochemical impact. However, given the notorious complexity of microbial cells and systems, which is many steps away from governing chemical or physical equations, how can we be sure that this activity is indeed predictable? Frentz et al.35 demonstrated that external conditions cause the seemingly random fluctuations observed in microbial growth, rather than stochastic variation in gene expression. This provides direct evidence of deterministic behavior, and so the authors conclude that microbial systems can in principle be determined by macroscopic laws.
    How is this determinism manifested? If microbial communities can respond relatively quickly to changes to their local environment, they may predictably optimize the exploitation of locally available resources. In the ocean, dispersal in microbes is thought to be a highly efficient process such that microbial communities can in effect draw from an extensive seed bank36,37, as captured in the phrase “everything is everywhere, the environment selects”38. Furthermore, recent evidence also shows that gene acquisitions and deletions happen quickly enough to allow for horizontal gene transfer to dominate bacterial adaptation39,40,41,42,43, implying that evolution can occur within few generations and thus on timescales similar to ecological interactions. Perhaps consequentially, similar geochemical environments have been demonstrated to have high-microbial functional redundancy despite different taxonomic compositions17,44. This may be interpreted with the hypothesis that physics and chemistry selects for metabolic traits, and that these traits can be housed in different organisms with taxonomic composition shaped by micro-scale or biotic interactions17,44,45.
    The prediction of microbial activity from environmental chemical potential has a long history in microbiology20,21,46,47,48,49,50,51,52, and is conceptually similar to other redox-balanced approaches to understanding microbial activity in sediments, soils, subsurfaces, and aquatic systems13,31,53,54,55,56,57,58. Illustrating the power of these approaches, anticipating metabolism from chemical potential resulted in a prediction that anaerobic ammonia oxidation (anammox) should exist decades before it was observed59,60. Quantitatively understanding microbially mediated rates of conversion of substrates has practical implications for wastewater treatment, and thus the field of biotechnology has established methodologies for an approach in textbook form21. Flux balance analysis (FBA) models can be considered as much more highly detailed analogs of this approach that resolve the mass and electron balances among a multitude of chemical reactions within a cell61,62.
    Redox-balanced metabolic functional types
    We can resolve microbial activity in global ecosystem models using the underlying redox chemistry of diverse metabolisms as a constraint. One specific way forward is to model distinct metabolisms as populations of metabolic functional types. This systematically quantifies relative rates of substrate consumption, biomass synthesis, and excretions of transformed products associated with each metabolism. Coupled with estimates of substrate uptake, this replaces implicit parameterizations of processes such as organic matter consumption, oxygen depletion, and denitrification with electron-balanced respiratory fluxes of dynamic microbial populations. Box 1 provides a detailed description of this methodology for multi-dimensional models.
    A particular set of redox reactions may distinguish a functional type, such as the oxidation of organic matter using oxygen (aerobic heterotrophy), or the oxidation of ammonia or nitrite using oxygen (chemoautotrophic nitrification) as exemplified in Table 1. For each metabolism, an electron-balanced description consists of multiple half-reactions: biomass synthesis, oxidation of an electron donor, and reduction of an electron acceptor21,48. The ratio of anabolism and catabolism can then be represented by the fraction f of electrons fueling cell synthesis vs. respiration for energy, following ref. 21. This provides a yield y (moles biomass synthesized per mole substrate utilized) of each required substrate that reflects two inputs: electron fraction f and the coefficients of the half-reactions (Fig. 2). The interlinked yields reflect the energy supplied by the redox reaction, the energy required for synthesis and other cellular demands, and the inefficiencies of energy conversion. Either f or y for any one of the substrates may be estimated theoretically with Gibbs free energies of reaction21 or with a combination of theoretical and empirical strategies63.
    Table 1 Simplified equations describing two exemplary metabolic functional types.
    Full size table

    Fig. 2: Schematic of a single cell represented as a metabolic functional type carrying out the aerobic oxidation of ammonia.

    The redox balance informs the elemental ratios of substrates utilized, biomass synthesized, and waste products excreted (Table 1).

    Full size image

    The result is a stoichiometric budget of the metabolism of the whole organism (Table 1). These descriptions quantify the elemental ratios of utilized substrates, biomass, and the excretion of waste products. For example, the descriptions account for the CO2 produced by heterotrophic metabolisms as well as the CO2 fixed by chemoautotrophic metabolisms (Table 1, Fig. 2, and Supplementary Fig. 2), linking microbial activity directly to global carbon cycling.
    To estimate the growth rate of each functional type, the yields from the metabolic budgets are combined with the uptake rates of the required substrates (Box 1). Limiting uptake rates may rely on empirically derived uptake kinetic parameters, or they can be estimated theoretically from diffusive supply, cell size, membrane physiology, and other physical constraints64,65,66. If theoretical models of uptake are used, the physical constraints on substrate acquisition and the redox chemical constraints on energy acquisition can provide an entirely theoretical estimate of the growth of each metabolic functional type.
    One strategy is to represent the populations carrying out each of these discrete metabolisms as one functional type population, which aggregates the diverse community of many species that are fueled by the same (or a similar) redox reaction (Fig. 3). Such aggregation has been deemed a useful strategy for representing the biogeochemical impacts of microbial communities for certain research questions67,68. However, for other questions this wipes out critical diversity among the aggregated populations. For example, diverse aerobic heterotrophic populations consume organic matter over a wide range of rates, and these rates dictate the amount of biologically sequestered carbon in the ocean. Redox chemistry and physical limitations alone may not inform the heterogeneity among similar metabolisms. One additional constraint is the limited capacity of the cell and thus its allocation of proteome towards different functions69. While the electrons supplied to the cell must be conserved following the redox balance, the electrons may be partitioned differently into machinery for substrate uptake vs. biomass synthesis, for instance, for different phenotypes. This partitioning can be quantitatively related to ecological fitness and biogeochemical impact via uptake kinetics, effective yields, and other traits9 (Supplementary Note 1 and Supplementary Fig. 1).
    Fig. 3: Solutions from a global simulation resolving multiple metabolic functional types.

    Net primary productivity (NPP), the biomasses of the metabolic functional types, and the sinking particulate organic carbon (POC) flux are resolved along a transect of a global microbial ecosystem model coupled with an estimate of the ocean circulation (Darwin-MITgcm18).

    Full size image

    Box 1 Incorporating metabolic functional types into ecosystem models

    A metabolic functional type can be represented as a population with a growth rate that is limited or co-limited by multiple required substrates. If Liebig’s Law of the Minimum is employed, the limiting growth rate μ is described as

    $$mu = min (V_{mathrm{i}}y_{mathrm{i}})$$
    (1)

    where Vi is the specific uptake rate of substrate i, and yield yi is the biomass yield with respect to that substrate. Yields for the different substrates and elements are interlinked in the metabolic budget derived from the underlying redox chemistry. Yields reflect Gibbs free energies of reaction among other factors. In the simplest model, non-limiting substrates are consumed in proportion to the limiting resource according to the metabolic budget, although in reality they may accumulate in the form of storage molecules.
    Each metabolic functional type population can be incorporated into a multi-dimensional environmental model (e.g., an ocean simulation) with physical transport as

    $$frac{{dB}}{{dt}} = mu B – Lleft( B right)B – underbrace {nabla cdot left( {bf{u}}{B} right)}_{{mathrm{advection}}} + underbrace {nabla cdot left( {{mathbf{kappa }}nabla B} right)}_{{mathrm{diffusion}}}$$
    (2)

    for biomass concentration B, loss rate L, velocity u, and diffusion coefficient κ. The loss rate function varies with biomass and represents a combination of processes, including predation, viral lysis, maintenance, and senescence. These processes remain largely unconstrained, although efforts have been made to relate losses to ecological dynamics107,108.
    The yield partitions the amount of substrate taken up by the population into that used for growth, Viyi, versus that exiting the cell in modified form as a waste product, Vi(1 − yi) (Fig. 2 and Table 1). Equation 1 suggests a correlation between μ and y, but yields may be further modified by other factors. For example, accounting for maintenance energy decreases the ratio of growth to respiration, contributing to a decoupling between growth rate and yield particularly at low growth rates109. Furthermore, a trade-off between uptake rate and yield at the cellular level reflects the allocation of enzyme towards machinery for substrate uptake vs. biomass synthesis, among other factors. Considering a proteome constraint can incorporate this trade-off (Supplementary Note 1).

    Benefits and implications for anticipating global change
    Redox chemistry aids in reducing the number of degrees of freedom in descriptions of diverse microbial metabolisms. We next discuss the benefits of this electron-balanced approach, each contextualized by specific projected impacts of global change due to microbial activity and broad challenges in the fields of microbial ecology and biogeochemistry.
    Flexible and broadly applicable metabolic thresholds: A key question for microbial biogeochemical studies, for which biogeochemical models are primed to answer, is how the biogeographies of diverse, active metabolisms vary with changes in the physical and chemical environment. What threshold determines the viability of a given metabolism?
    Redox-balanced metabolic budgets obviate the need to impose critical concentrations or other thresholds that determine the presence of any given metabolism. Rather than being imposed following empirical relationships, metabolic biogeography emerges dynamically from ecological interactions and reflects environmental chemical potential. This flexibility aids in understanding metabolic thresholds more fundamentally, and it expands model applicability to diverse and unobserved environments. This is of particular importance for understanding global change, since past and future worlds may include very different ecosystems that do not reflect current empirical trends.
    For example, the oceans are currently losing oxygen due to global warming5,70. If anoxic zones continue to expand, this will increase the habitat of anaerobic microorganisms, whose respiration results in emissions of N2 and N2O to the atmosphere7. Many biogeochemical models prescribe O2 concentrations that inhibit anaerobic activity in accordance with observations of specific organisms or communities in experimental conditions. This assumes that the same O2 concentrations limit metabolism similarly in all environments, and often trades mechanistic understanding of oxygen limitation for empirical correlations that may reflect a variety of natural and introduced biases, such as micro-scale heterogeneity, physical mixing in the ocean, and experimental bottle effects.
    In contrast, a metabolic functional type model does not require imposed oxygen threshold concentrations (Supplementary Fig. 3). When oxygen supply is abundant, anaerobic types are competitively excluded because growth using alternative electron acceptors is lower than with oxygen. When oxygen supply is low, aerobic populations may persist and continue to deplete any available oxygen even as their growth is limited by oxygen, allowing for a steady state stable coexistence of aerobic and anaerobic metabolisms, which is consistent with a variety of observations71.
    Descriptions of microbial growth that reflect underlying chemical potential can enable predictions of many other metabolic transitions, such as nitrogen fixation, nitrification, and the transition to sulfur oxidation and reduction13,18,72,73,74. As another example, this approach predicts the restriction of nitrification from the sunlit surface ocean as a consequence of competitive exclusion by phytoplankton in many environments (Fig. 4 and Supplementary Fig. 2), as well as active nitrification in some surface locations where phytoplankton are limited by another factor not affecting the chemoautotrophs, such as at high latitudes where phytoplankton are limited by light18. The emergent exclusion from most of the surface ocean anticipates that many clades of nitrifying microorganisms have adapted to long-term exclusion from the surface and consequentially lost (or did not develop) photoprotective cellular machinery.
    Fig. 4: Model simulation and observations of the marine nitrification system.

    Biogeochemistry is driven by microbial metabolic functional types in a vertical water column model18. Lines are model solutions, and marked points are observations from two stations in the Pacific Ocean75,106 (see Supplementary Fig. 2 for more detail) (a). Chlorophyll a concentrations and abundances of ammonia-oxidizing organisms (AOO) and nitrite-oxidizing organisms (NOO). Observed abundances are of the 16S rRNA abundances of archaeal Marine Group I and Nitrospina-like bacteria75,106. Model abundances are converted from biomass with 0.1 fmol N cell-1 for AOO, 0.2 fmol N cell-1 for NOO76, and one gene copy per cell. (b). Light (solar irradiance) and bulk nitrification rates.

    Full size image

    Replacing implicit descriptions of organic matter remineralization: The fate of organic matter dictates the amount of carbon sequestered in the marine and terrestrial biospheres. Microbial consumption mediates the carbon stored in soils, the carbon stored in the ocean as dissolved organic matter (DOM), and the sinking flux of organic carbon that constitutes the marine biological carbon pump4, without which atmospheric CO2 would be 100–200 ppm higher than current levels. We want to understand how these carbon reservoirs respond to changes in climate, such as increased temperatures and changes in precipitation patterns. However, in biogeochemical models, simple rate constants often dictate the remineralization of elements from organic back into inorganic constituents.
    Replacing simplistic parameterizations with dynamic metabolic functional types means that electron-balanced descriptions of growth and respiration instead drive the fate of organic matter in earth system models (Fig. 3). In addition to a more sophisticated and responsive description of carbon sequestration, non-living organic matter is fully integrated into ecosystem frameworks, enabling theoretical studies of phytoplankton-bacteria interactions to complement observational and experimental approaches.
    Much work remains in the development of these descriptions. As we discuss below, accurate estimates of organic matter turnover rates require more accurate descriptions of the complex processes governing microbial uptake rates of organic matter. However, the redox-informed yields are still useful for quantifying the relative amount of CO2 excreted and the absolute amount of biomass sustained on a given substrate, independent of uptake kinetics (Supplementary Note 2).
    Relationships between abundances, rates, nutrient concentrations, and elemental ratios: An overarching puzzle challenging microbial ecology is to understand how chemical transformations in the environment are set by the ecological interactions at the organism level, among individual microscopic cells. It is clear that abundances of populations are not simply and directly correlated with biogeochemical impact (i.e., higher abundance does not necessarily imply an associated higher rate of chemical transformation). Untangling the relationship between abundances and biogeochemical function is also necessary for interpretation of genetic evidence that provide insight into this complex ecosystem structure.
    Redox-balanced metabolic functional type modeling links rates of biomass synthesis associated with a particular metabolism to its rate of respiration as well as the standing stock of limiting nutrients. As functional type modeling is coupled with estimates of population loss rates due to grazing, viral lysis, or other mortality, simulations also resolve the standing stocks of functional biomass. This quantifies the relationship between biomass concentrations and volumetric rates of chemical transformations, emphasizing how relatively low biomass may be associated with relatively high bulk rates71.
    For example, the approach has revealed a clear example of the signature of chemical potential in the ecology of marine nitrification18 (Fig. 4). In this model, the two steps of nitrification are represented by two functional type populations. This predicts about a three-fold difference in the abundances of the organisms responsible for each of the two steps of nitrification, despite the fact that the two populations carry out the same rate of subsurface N-cycling at steady state6,18. A three-fold or greater difference in abundance and associated ammonium (NH4+) and nitrite (NO2−) concentrations is consistent with observed differences18,75, and it reflects that the oxidation of one mole of NH4+ generates three times more electrons than the oxidation of one mole of NO2−, with differences in cell size further contributing to differences in abundances (Fig. 4 and Supplementary Fig. 2). Recent observations confirm the redox-based difference in NH4+ and NO2− biomass yield76,77, although measured rates from a nonsteady environment suggest that NO2−-oxidizing bacteria can partition electrons more efficiently than NH4+-oxidizing archaea76 (i.e., higher fraction f despite lower yield y; see Supplementary Note 3).
    As redox-based descriptions resolve the stoichiometry of whole organism metabolism, they also link together elemental cycles. Explicit description of relative elemental flow through the ecosystem, and specifically their variation from average values, is critical for understanding climate-biogeochemical feedbacks78,79,80. For example, the nitrification model also estimates the CO2 fixation rates associated with nitrification rates (Supplementary Fig. 2), enabling global-scale, electron-balanced projections of the amount of carbon converted to organic form by chemoautotrophic nitrifying microorganisms.
    Connections with sequencing datasets: How do we relate metabolic functional type models to sequencing datasets measuring genetic, transcriptomic, and proteomic diversity? Connecting biogeochemical models with sequencing data is critical because this data provides an enormous amount of information about ecosystem structure and function. Genes (or transcripts) themselves are not necessarily the most concise or useful currency given functional redundancies, unattributed function, and variation in gene dosage from horizontal gene transfer as well as growth rate40. Recent gene-centric models aim to resolve the abundances of key genes as proxies for a predetermined set of metabolic pathways13,14,17. However, the parameters used to describe metabolic pathways in these models are estimated similarly to the redox-balanced yields and efficiencies described here.
    The innovation of gene-centric models is the sophisticated conversion of estimates of biogeochemical activity and biomass to genes. For example, the model of Coles et al.17 resolves biomass and nutrient concentrations prognostically, and then uses a three-part formula—representing constitutive, regulated, and steady state transcription—to diagnostically calculate transcription rates from modeled biomass and growth rates17. Thus, the two types of modeling are complimentary, with redox chemistry providing estimates of metabolic activity from fundamental principles, and the careful calibrations between activity and sequencing providing a comparative metric.
    The examples here externalize the conversion between modeled activity and sequencing information as a transparent process. In Fig. 4, the predicted functional biomass of ammonia-oxidizing population is related to archaeal Marine Group I (MGI) and Nitrospina-like 16S rRNA genes with two conversion factors: the cell elemental quota (fmol N cell−1) and the number of cellular gene copies. Conversion error arises since cell mass and size vary with growth rate81,82. Maintaining transparency of the conversion from predicted microbial activity to genes and transcripts allows interdisciplinary audiences to understand and critique the models.
    Limitations and possible extensions
    Using chemical potential as a theoretically grounding organizing principle for the resolution of diverse metabolisms can greatly improve microbial descriptions in global biogeochemical models. However, the approach does have its limitations, which generally increase in significance with increased temporal or spatial resolution.
    Modeling metabolic diversity with functional type populations requires choosing how metabolisms are distributed among the populations. This has consequences when interpreting time-varying states: model solutions become dependent on the partitioning of metabolism among the functional types as the timescales of physical change approach the timescales of microbial growth (see Supplementary Note 4, Supplementary Fig. 3, and Supplementary Fig. 4 for a detailed example). Other species-specific time-varying phenomena such as the lag response of organisms to substrate availability also become relevant83. On one hand, this is beneficial for resolution of microbial processes in fine-grained ocean circulation models where flow can vary on the order of days. However, incorporating another constraint, such as proteome allocation69, is necessary to inform these choices. For example, considering enzymatic allocation in combination with energetics allowed for the prediction of both the division of nitrification into a two-step process in mixed environments and the combined, complete pathway in one organism (comammox) in biofilms, which preceded observations of the latter84,85,86.
    Uncertainty in distributions of metabolism lies not only in the length of a metabolic pathway, but also in the degree of metabolic versatility (metabolic mixotrophy). Such versatility characterizes key players in large-scale biogeochemistry, such as nitrite-oxidizing bacteria and photoheterotrophs87,88,89. Mixotrophic lifestyles can increase the fitness of populations in their environments, impacting overall ecosystem function90. In one sense, the approach here provides a prediction of where we might expect such mixotrophy by resolving stable coexistences of diverse metabolisms. In Fig. 3, for example, syntrophic coexistence occurs at depth among heterotrophs, ammonia oxidizers, and nitrite oxidizers, and future work could investigate what determines which combinations of these coexistences remain as passive interactions, which develop into mutualistic dependencies as active interactions91, and which evolve into mixotrophic phenotypes or endosymbionts. Additionally, by considering the potential to carry out a metabolism as a trait, we can use the current framework along with an additional constraint to investigate implications of metabolic mixotrophy. For example, Coles et al.17 impose a trade-off between the degree of metabolic diversity of a single functional type and growth rate, enabling the exploration the consequences of distribution of metabolism on the biogeochemical state.
    Also, the metabolic functional type approach resolves only active functional biomass, while evidence suggests that less than 10% to more than 75% of the microbial community may be inactive92. Some seemingly inactive populations may slowly metabolize over long timescales, requiring longer model integration times and careful attention to their loss rates for resolution, while some populations are periodically active as revealed by high-resolution observations in time37.
    The proposed modeling approach relies on estimates of the limiting uptake rates of required substrates. In lieu of suitable theoretical descriptions, the use of empirically derived uptake kinetic parameters still employs the benefits of the redox-informed yields (Supplementary Note 2). However, underlying physical constraints to substrate acquisition can in principle be exploited to develop more universally applicable descriptions for a variety of substrates and contexts. Uptake kinetics are complex, but for many limiting resources, encounter effectively controls the uptake, and the physics of encounter has been relatively well described. For example, uptake rates estimated from diffusive supply of substrate, cellular geometry, and membrane physiology64,65,66 have been empirically supported93. For organic matter, future work is needed to develop suitable descriptions of consumption rates, whether empirical or theoretical. For example, descriptions require attention to the hydrolysis of organic compounds by extracellular enzymes and the ecology of sinking marine particles—the diffusion of monomer away from the particle, within-particle transport, and dynamic ecological interactions on particle surfaces, among other processes11,74,94,95.
    In Fig. 3, the electron-balanced description consists of an average stoichiometry and electron fraction for one sinking pool of organic matter in the ocean, which, as mentioned above, is not sufficient to accurately resolve the carbon storage that is shaped by a distribution of turnover rates. As one of the many factors impacting the rates, an energetics-based perspective can serve as a tool for further deciphering organic matter complexity. For example, organic matter may be partially organized by the nominal oxidation state of its carbon atoms, which relates to a measure of free energy and accessibility96,97. This could be used to improve the phenomenological description of organic matter in models as labile vs. non-labile, for example, with a more mechanistic underpinning.
    Descriptions of phytoplankton are currently much more sophisticated than of bacteria and archaea in models, reflecting a longer history of comprehensive sets of observations. However, further work could develop simple descriptions of photoautotrophic metabolisms from underlying energetics by connecting the supply of photons to available energy for biosynthesis within the cell. Many biogeochemical models account for an inefficiency of phytoplankton metabolism with a parameter that dictates their excretions of dissolved organic matter98. Incorporating this excretion into an energetic framework would enhance studies of phytoplankton ecology, such as studies of photoautotrophic-heterotrophic interactions in the ocean surface or photoautotrophic-chemoautotrophic interactions at the base of the euphotic zone where some phytoplankton excrete nitrite due to incomplete reduction of nitrate99.
    As a more radical extension, can we progress past population modeling and model microbial consortia as one aggregate community biomass34,100? This may improve resolution of time-varying metabolic versatility. However, if both steps of nitrification were a part of such a consortium, would the characteristic accumulation of nitrite be predicted (Supplementary Fig. 2)? We leave these questions for future research and conclude that the best choice for the degree of resolution of metabolism will depend on the specific research question and the available observations.
    We have described a useful approach for understanding and anticipating microbial control of biogeochemical cycling that is suitable for global applications. The approach aims to represent microbial growth and respiration explicitly and consistently from knowledge of chemical gradients in the environment, towards a goal of building an independently constructed theoretical ecosystem model that can then be compared to observations. Describing microbial communities with underlying energetic constraints connects metabolisms dynamically with global geochemical distributions, such as those of carbon dioxide, oxygen, and biologically available nitrogen. This deepens our understanding of microbial ecosystems and enables the incorporation of the feedbacks of microbial activity to changes in global biogeochemistry and the climate system. More

  • in

    Experimental tests of bivalve shell shape reveal potential tradeoffs between mechanical and behavioral defenses

    Experimental limitations and implications
    Many factors are known to contribute to mollusk shell strength including thickness, microstructure, previous shell damage, and ornamentation19,21,29,30,31. Thus, it is challenging to separate the influence of one feature from another when testing the strength of real shells. Fortunately, there are alternative methods by which to model shells22,32,33,34. 3D printing is an effective tool to create model shells that exhibit brittle behavior in compression, serving as an accurate first-order approximation for shell breakage (Fig. 2). 3D printed models are not created to replicate shell microstructure; rather they serve to normalize confounding factors encountered with real shells, like size, variations in thickness, and taphonomy/degradation of microstructure, to isolate variables of interest while maintaining the predominantly brittle behavior seen in real shells. As a result, the failure of different 3D printed shapes under bulk mechanical compression can be used for relative comparisons between morphologies without needing to reproduce the exact magnitude of load to failure of real shells22,33 (see “Methods”, Fig. 2). Natural shells ultimately break after cracking through microstructural layers. Analogously, in this study, 3D printed shells cracked predominantly through layers of printed materials, not along printing boundaries, demonstrating an appropriate proxy for shell failure (Figs. 1, 2).
    Resistance to compression by a predator can be tested directly using a variety of loading tests, the effectiveness of which has been demonstrated repeatedly by investigators studying topics ranging from predation to climate change22,33,34. Compression experiments of bivalve shells, specifically, have been used extensively for this purpose35,36,37,38,39,40,41,42,43,44,45,46. In this study, the mechanical strength of model bivalve shell shapes was analyzed to understand potential defensive value against durophagous (crushing) predators using hypothetical flat teeth and jaws to crush prey (see Crofts and Summers (2014)22 for examples of flat crushing morphologies) (Figs. 1, 2). These experiments are most analogous to small shell-crushers with flat dentitions, like fishes (e.g., guitar fishes, stingrays, etc.), rather than large predators like the modern walrus or extinct marine reptiles (e.g., placodonts, mosasaurs) which are so much larger than their prey that the differences in shell strength resulting from shape are likely insignificant. Bivalves are also preyed upon by many other predators including asteroids, gastropods, birds, and mammals47. The ability to escape is likely dependent upon predator capabilities, where not all escape mechanisms are equally effective against different predators. The results of these experiments are not intended to be general proxies of predation resistance—they are only applicable as a proxy for predators that use flat crushing dentitions. Additionally, it is important to note that compression of a shell by vertebrate predators is a different mechanical process from compression by an arthropod, as claws localize forces differently than teeth and jaws48. Therefore, the experiments used in this study are not meant to model predation by invertebrate durophages. This experimental setup represents an idealized case of shell compression assuming consistent shell thickness and a common predator using a quasistatic loading regime to isolate the influence of shell shape. Therefore, this study makes no conclusions as to the effects of shell shape on strength under conditions of point loading by claws or impact. Furthermore, while bivalve shape is undoubtedly influenced by many factors including fabrication (shell growth and construction)49, phylogenetics, location of soft tissue, etc., these experiments were designed to study an idealized case in which defense against vertebrate shell crushing predators is the most important functional constraint on shell shape.
    Advantageously, the use of mathematically generated theoretical bivalve shells, rather than real shells, enables testing a range of shapes—some of which can be or have been found in nature—and others that have yet to exist. For example, these methods enable physical testing of morphologies which are only found in the fossil record. Fossils cannot be used for mechanical experimentation due to changes in the integrity of the shell structure resulting from the fossilization process. Furthermore, testing shapes which do not exist due to biological constraints like fabrication, can provide valuable insights when combined with the study of shapes that have evolved naturally. For example, testing theoretical shapes can reveal morphologies that perform better than natural morphologies in specific functional settings. Thus, testing shapes which perform well, but do not exist, is one way to identify the potential influence of evolutionary constraints such as fabrication (shell growth patterns), phylogenetics, or other necessary biological functions. In this study, the extremely perpendicular-elongate shells represent theoretical shapes that are not seen in nature (likely due to space limitations for soft tissue). Developing theoretical physical models also allows for the elimination of many confounding variables that cannot be avoided when testing real shells. The effects of variations in shell thickness and ornamentation, which also play a role in shell strength, can thus be separated from the strength imparted by gross shell shape when theoretical models are generated. However, because overall shell strength is derived from a combination of these many factors (microstructure, thickness, etc.31.) the conclusions from these experiments represent a first-order approximation of the effects of shape on shell strength alone. Because the models for this study were generated mathematically, and are therefore not exact replicas of specific taxa, the discussion below aims to address potential tradeoffs based on general shell shapes that could be further studied for specific bivalve taxa in future experiments.
    Shell shape and strength
    Three parameters of shell shape were modified for this study. (1) Generating curve shape is both a modeling parameter and biological feature of mollusk shells27. In mathematical models, generating curve shape describes the shape that is rotated around an axis to create a surface that represents the overall shell shape. Biologically, the generating curve refers to the portion of the shell upon which new material is secreted by the mantle. In bivalves, this is the shape of the commissure. Thus, by changing the shape of the generating curve in a model bivalve, i.e., changing the shape of the ellipse that generates the surface of the shell, models can be perpendicular-elongated (a  > b) or parallel-elongated (a  More

  • in

    Spatial patterns of pathogen prevalence in questing Ixodes ricinus nymphs in southern Scandinavia, 2016

    1.
    Pfäffle, M., Littwin, N., Muders, S. V. & Petney, T. N. The ecology of tick-borne diseases. Int. J. Parasitol. 43, 1059–1077 (2013).
    Article  Google Scholar 
    2.
    Han, B. A. & Yang, L. Predicting novel tick vectors of zoonotic disease. in ICML Workshop on #Data4Good: Machine Learning in Social Good Applications 71–75 (2016).

    3.
    de la Fuente, J., Estrada-Pena, A., Venzal, J. M., Kocan, K. M. & Sonenshine, D. E. Overview: ticks as vectors of pathogens that cause disease in humans and animals. Front. Biosci. 13, 6938–6946 (2008).
    Article  Google Scholar 

    4.
    Michelet, L. et al. High-throughput screening of tick-borne pathogens in Europe. Front. Cell. Infect. Microbiol. 4, 103 (2014).
    Article  Google Scholar 

    5.
    Estrada-Peña, A. & de la Fuente, J. The ecology of ticks and epidemiology of tick-borne viral diseases. Antiviral Res. 108, 104–128 (2014).
    Article  CAS  Google Scholar 

    6.
    Paul, R. E. L. et al. Environmental factors influencing tick densities over seven years in a French suburban forest. Parasit. Vectors 9, 309 (2016).
    Article  CAS  Google Scholar 

    7.
    Randolph, S. E. Tick-borne disease systems emerge from the shadows: the beauty lies in molecular detail, the message in epidemiology. Parasitology 136, 1403 (2009).
    CAS  Article  Google Scholar 

    8.
    Jore, S. et al. Multi-source analysis reveals latitudinal and altitudinal shifts in range of Ixodes ricinus at its northern distribution limit. Parasit. Vectors 4, 1–11 (2011).
    Article  Google Scholar 

    9.
    Bernstein, L. et al. Climate Change 2007: Synthesis Report. Contribution of working groups I, II and III to the fourth assessment report of the Intergovernmental Panel on Climate Change. Intergovernmental Panel on Climate Change [Core Writing Team, Pachauri, R.K and Reisinger, A. (eds.)] https://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr_full_report.pdf (2007).

    10.
    Kovats, R. S., Campbell-Lendrum, D. H., McMichael, A. J., Woodward, A. & Cox, J. S. Early effects of climate change: do they include changes in vector-borne disease?. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 356, 1057–1068 (2001).
    CAS  PubMed  PubMed Central  Google Scholar 

    11.
    Gage, K. L., Burkot, T. R., Eisen, R. J. & Hayes, E. B. Climate and vectorborne diseases. Am. J. Prev. Med. 35, 436–450 (2008).
    Article  Google Scholar 

    12.
    Medlock, J. M. et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit. Vectors 6, 1–11 (2013).
    Article  Google Scholar 

    13.
    Andreassen, A. et al. Prevalence of tick borne encephalitis virus in tick nymphs in relation to climatic factors on the southern coast of Norway. Parasit. Vectors 5, 1–12 (2012).
    Article  Google Scholar 

    14.
    Soleng, A. et al. Distribution of Ixodes ricinus ticks and prevalence of tick-borne encephalitis virus among questing ticks in the Arctic Circle region of northern Norway. Ticks Tick. Borne. Dis. 9, 97–103 (2018).
    CAS  Article  Google Scholar 

    15.
    Kjelland, V. et al. Tick-borne encephalitis virus, Borrelia burgdorferi sensu lato, Borrelia miyamotoi, Anaplasma phagocytophilum and Candidatus Neoehrlichia mikurensis in Ixodes ricinus ticks collected from recreational islands in southern Norway. Ticks Tick. Borne. Dis. 9, 1098–1102 (2018).
    Article  Google Scholar 

    16.
    Paulsen, K. M. et al. Prevalence of tick-borne encephalitis virus in Ixodes ricinus ticks from three islands in north-western Norway. APMIS 123, 759–764 (2015).
    Article  Google Scholar 

    17.
    Kjær, L. J. et al. A large-scale screening for the taiga tick, Ixodes persulcatus, and the meadow tick, Dermacentor reticulatus, in southern Scandinavia, 2016. Parasit. Vectors 12, 338 (2019).
    Article  Google Scholar 

    18.
    Oechslin, C. P. et al. Prevalence of tick-borne pathogens in questing Ixodes ricinus ticks in urban and suburban areas of Switzerland. Parasit. Vectors 10, 558 (2017).
    Article  CAS  Google Scholar 

    19.
    Becker, N. S. et al. Recurrent evolution of host and vector association in bacteria of the Borrelia burgdorferi sensu lato species complex. BMC Genom. 17, 734 (2016).
    Article  CAS  Google Scholar 

    20.
    Bowman, A. S. & Nuttall, P. A. Ticks: Biology, Disease and Control (Cambridge University Press, Cambridge, 2004).
    Google Scholar 

    21.
    Hasle, G. et al. Transport of ticks by migratory passerine birds to Norway. J. Parasitol. 95, 1342–1351 (2009).
    PubMed  Google Scholar 

    22.
    Klitgaard, K. et al. Screening for multiple tick-borne pathogens in Ixodes ricinus ticks from birds in Denmark during spring and autumn migration seasons. Ticks Tick. Borne. Dis. 10, 546–552 (2019).
    PubMed  Google Scholar 

    23.
    Skarphédinsson, S. et al. Detection and identification of Anaplasma phagocytophilum, Borrelia burgdorferi, and Rickettsia helvetica in Danish Ixodes ricinus ticks. APMIS 115, 225–230 (2007).
    PubMed  Google Scholar 

    24.
    Fraenkel, C.-J., Garpmo, U. & Berglund, J. Determination of novel Borrelia genospecies in Swedish Ixodes ricinus ticks. J. Clin. Microbiol. 40, 3308–3312 (2002).
    CAS  PubMed  PubMed Central  Google Scholar 

    25.
    Wilhelmsson, P. et al. Prevalence, diversity, and load of Borrelia species in ticks that have fed on humans in regions of Sweden and Åland Islands, Finland with different Lyme borreliosis incidences. PLoS ONE 8, e81433 (2013).
    ADS  PubMed  PubMed Central  Google Scholar 

    26.
    Vennestrøm, J., Egholm, H. & Jensen, P. M. Occurrence of multiple infections with different Borrelia burgdorferi genospecies in Danish Ixodes ricinus nymphs. Parasitol. Int. 57, 32–37 (2008).
    PubMed  Google Scholar 

    27.
    Kjelland, V., Stuen, S., Skarpaas, T. & Slettan, A. Prevalence and genotypes of Borrelia burgdorferi sensu lato infection in Ixodes ricinus ticks in southern Norway. Scand. J. Infect. Dis. 42, 579–585 (2010).
    CAS  PubMed  Google Scholar 

    28.
    Klitgaard, K., Kjær, L. J., Isbrand, A., Hansen, M. F. & Bødker, R. Multiple infections in questing nymphs and adult female Ixodes ricinus ticks collected in a recreational forest in Denmark. Ticks Tick. Borne. Dis. 10, 1060–1065 (2019).
    PubMed  Google Scholar 

    29.
    Maraspin, V., Ruzic-Sabljic, E. & Strle, F. Lyme borreliosis and Borrelia spielmanii. Emerg. Infect. Dis. 12, 1177–1177 (2006).
    PubMed  PubMed Central  Google Scholar 

    30.
    Rudenko, N., Golovchenko, M., Grubhoffer, L. & Oliver, J. H. Updates on Borrelia burgdorferi sensu lato complex with respect to public health. Ticks Tick. Borne. Dis. 2, 123–128 (2011).
    Article  Google Scholar 

    31.
    Fertner, M. E., Mølbak, L., Pihl, T. P. B., Fomsgaard, A. & Bødker, R. First detection of tick-borne “Candidatus Neoehrlichia mikurensis” in Denmark 2011. Eurosurveillance 17, 20096 (2012).
    Google Scholar 

    32.
    Quarsten, H. et al. Candidatus Neoehrlichia mikurensis and Borrelia burgdorferi sensu lato detected in the blood of Norwegian patients with erythema migrans. Ticks Tick. Borne. Dis. 8, 715–720 (2017).
    CAS  Article  Google Scholar 

    33.
    Stuen, S., Granquist, E. G. & Silaghi, C. Anaplasma phagocytophilum—a widespread multi-host pathogen with highly adaptive strategies. Front. Cell. Infect. Microbiol. 3, 31 (2013).
    Article  CAS  Google Scholar 

    34.
    Fomsgaard, A. et al. Tick-borne encephalitis virus, Zealand, Denmark, 2011. Emerg. Infect. Dis. 19, 1171–1173 (2013).
    Article  Google Scholar 

    35.
    Jensen, P. M. et al. Transmission differentials for multiple pathogens as inferred from their prevalence in larva, nymph and adult of Ixodes ricinus (Acari: Ixodidae). Exp. Appl. Acarol. 71, 171–182 (2017).
    Article  Google Scholar 

    36.
    Lundkvist, Å., Wallensten, A., Vene, S. & Hjertqvist, M. Tick-borne encephalitis increasing in Sweden, 2011. Eurosurveillance 16, 19981 (2011).
    Article  Google Scholar 

    37.
    Svensson, J., Hunfeld, K.-P. & Persson, K. E. M. High seroprevalence of Babesia antibodies among Borrelia burgdorferi-infected humans in Sweden. Ticks Tick. Borne. Dis. 10, 186–190 (2019).
    Article  Google Scholar 

    38.
    Mørch, K., Holmaas, G., Frolander, P. S. & Kristoffersen, E. K. Severe human Babesia divergens infection in Norway. Int. J. Infect. Dis. 33, 37–38 (2015).
    Article  Google Scholar 

    39.
    Uhnoo, I. et al. First documented case of human babesiosis in Sweden. Scand. J. Infect. Dis. 24, 541–547 (2009).
    Article  Google Scholar 

    40.
    Dumler, J. S., Barat, N. C., Barat, C. E. & Bakken, J. S. Human granulocytic anaplasmosis and macrophage activation. Clin. Infect. Dis. 45, 199–204 (2007).
    CAS  Article  Google Scholar 

    41.
    Nilsson, K., Elfving, K. & Påhlson, C. Rickettsia helvetica in patient with meningitis, Sweden, 2006. Emerg. Infect. Dis. 16, 490–492 (2010).
    CAS  Article  Google Scholar 

    42.
    Frivik, J. O., Noraas, S., Grankvist, A., Wennerås, C. & Quarsten, H. En mann i 60-årene fra Sørlandet med intermitterende feber (In Norwegian). Tidsskr. Den Nor. legeforening 137, (2017).

    43.
    Grankvist, A. et al. Infections with the tick-borne bacterium ‘Candidatus Neoehrlichia mikurensis’ mimic noninfectious conditions in patients with B cell malignancies or autoimmune diseases. Clin. Infect. Dis. 58, 1716–1722 (2014).
    CAS  Article  Google Scholar 

    44.
    Welinder-Olsson, C., Kjellin, E., Vaht, K., Jacobsson, S. & Wenneras, C. First case of human ‘Candidatus Neoehrlichia mikurensis’ infection in a febrile patient with chronic lymphocytic leukemia. J. Clin. Microbiol. 48, 1956–1959 (2010).
    Article  Google Scholar 

    45.
    Rizzoli, A. et al. Ixodes ricinus and its transmitted pathogens in urban and peri-urban areas in Europe: new hazards and relevance for public health. Front. Public Health 2, 251 (2014).
    Article  Google Scholar 

    46.
    Michelitsch, A., Wernike, K., Klaus, C., Dobler, G. & Beer, M. Exploring the reservoir hosts of tick-borne encephalitis virus. Viruses vol. 11 (2019).

    47.
    Keesing, F. et al. Reservoir competence of vertebrate hosts for Anaplasma phagocytophilum. Emerg. Infect. Dis. 18, 2013–2016 (2012).
    Article  Google Scholar 

    48.
    Zhan, L. et al. Anaplasma phagocytophilum in livestock and small rodents. Vet. Microbiol. 144, 405–408 (2010).
    ADS  Article  Google Scholar 

    49.
    Portillo, A., Santibáñez, P., Palomar, A. M., Santibáñez, S. & Oteo, J. A. Candidatus Neoehrlichia mikurensis, Europe. New Microbes New Infect. 22, 30–36 (2018).
    CAS  Article  Google Scholar 

    50.
    Jenkins, A. et al. Detection of Candidatus Neoehrlichia mikurensis in Norway up to the northern limit of Ixodes ricinus distribution using a novel real time PCR test targeting the groEL gene. BMC Microbiol. 19, 199 (2019).
    PubMed  PubMed Central  Google Scholar 

    51.
    Obiegala, A. & Silaghi, C. Candidatus Neoehrlichia mikurensis—recent insights and future perspectives on clinical cases, vectors, and reservoirs in Europe. Curr. Clin. Microbiol. Rep. 5, 1–9 (2018).
    Google Scholar 

    52.
    Yabsley, M. J. & Shock, B. C. Natural history of zoonotic Babesia: role of wildlife reservoirs. Int. J. Parasitol. Parasites Wildl. 2, 18–31 (2013).
    PubMed  Google Scholar 

    53.
    Sprong, H. et al. Ixodes ricinus ticks are reservoir hosts for Rickettsia helvetica and potentially carry flea-borne Rickettsia species. Parasit. Vectors 2, 41 (2009).
    PubMed  PubMed Central  Google Scholar 

    54.
    Jaenson, T. G. T. et al. Risk indicators for the tick Ixodes ricinus and Borrelia burgdorferi sensu lato in Sweden. Med. Vet. Entomol. 23, 226–237 (2009).
    CAS  PubMed  Google Scholar 

    55.
    Hudson, P. J. et al. Tick-borne encephalitis virus in northern Italy: molecular analysis, relationships with density and seasonal dynamics of Ixodes ricinus. Med. Vet. Entomol. 15, 304–313 (2001).
    MathSciNet  CAS  Article  Google Scholar 

    56.
    Nazzi, F. et al. Ticks and Lyme borreliosis in an alpine area in northeast Italy. Med. Vet. Entomol. 24, 220–226 (2010).
    CAS  PubMed  Google Scholar 

    57.
    Hubalek, Z., Halouzka, J. & Juricova, Z. Longitudinal surveillance of the tick Ixodes ricinus for Borreliae. Med. Vet. Entomol. 17, 46–51 (2003).
    CAS  Article  Google Scholar 

    58.
    Lindström, A. & Jaenson, T. G. T. Distribution of the common tick, Ixodes ricinus (Acari: Ixodidae), in different vegetation types in southern Sweden. J. Med. Entomol. 40, 375–378 (2003).
    Article  Google Scholar 

    59.
    Mejlon, H. A. & Jaenson, T. G. T. Jaenson (1993) Seasonal prevalence of Borrelia burgdorferi in Ixodes ricinus in different vegetation types in Sweden. Scand. J. Infect. Dis. 25, 449–456 (2009).
    Article  Google Scholar 

    60.
    Tack, W. et al. Local habitat and landscape affect Ixodes ricinus tick abundances in forests on poor, sandy soils. For. Ecol. Manag. 265, 30–36 (2012).
    Google Scholar 

    61.
    Walhström, L. K. & Kjellander, P. Ideal free distribution and natal dispersal in female roe deer. Oecologia 103, 302–308 (1995).
    ADS  PubMed  Google Scholar 

    62.
    Zeman, P. Objective assessment of risk maps of tick-borne encephalitis and Lyme borreliosis based on spatial patterns of located cases. Int. J. Epidemiol. 26, 1121–1129 (1997).
    CAS  PubMed  Google Scholar 

    63.
    Jat, M. K. & Mala, S. Application of GIS and space-time scan statistic for vector born disease clustering. In ICEGOV ’17 Proceedings of the 10th International Conference on Theory and Practice of Electronic Governance (2017) https://doi.org/10.1145/3047273.3047361.

    64.
    Hönig, V. et al. Model of risk of exposure to Lyme borreliosis and tick-borne encephalitis virus-infected ticks in the border area of the Czech Republic (South Bohemia) and Germany (Lower Bavaria and Upper Palatinate). Int. J. Environ. Res. Public Health 16, 1173 (2019).
    PubMed Central  Google Scholar 

    65.
    Randolph, S. E. & Rogers, D. J. Fragile transmission cycles of tick-borne encephalitis virus may be disrupted by predicted climate change. Proc. Biol. Sci. 267, 1741–1744 (2000).
    CAS  PubMed  PubMed Central  Google Scholar 

    66.
    Stefanoff, P. et al. A Predictive model has identified tick-borne encephalitis high-risk areas in regions where no cases were reported previously, Poland, 1999–2012. Int. J. Environ. Res. Public Health 15, 677 (2018).
    PubMed Central  Google Scholar 

    67.
    Kjær, L. J. et al. Predicting and mapping human risk of exposure to Ixodes ricinus nymphs using climatic and environmental data, Denmark, Norway and Sweden, 2016. Eurosurveillance 24, 1800101 (2019).
    PubMed Central  Google Scholar 

    68.
    Kjær, L. J. et al. Predicting the spatial abundance of Ixodes ricinus ticks in southern Scandinavia using environmental and climatic data. Sci. Rep. 9, 18144 (2019).
    ADS  Google Scholar 

    69.
    Kjær, L. J. et al. Spatial data of Ixodes ricinus instar abundance and nymph pathogen prevalence, Scandinavia, 2016–2017. Collection https://doi.org/10.6084/m9.figshare.c.4938270.v1 (2020).
    Article  Google Scholar 

    70.
    Kjær, L. J. et al. Spatial data of Ixodes ricinus instar abundance and nymph pathogen prevalence, Scandinavia, 2016–2017. Sci. Data 7, 1–7 (2020).
    Google Scholar 

    71.
    Scharlemann, J. P. W. et al. Global data for ecology and epidemiology: a novel algorithm for temporal Fourier processing MODIS data. PLoS ONE 3, e1408 (2008).
    ADS  Article  Google Scholar 

    72.
    Corine Land Cover 2006 raster data. European Environment Agency https://www.eea.europa.eu/data-and-maps/data/clc-2006-raster (2010).

    73.
    Klitgaard, K., Chriél, M., Isbrand, A., Jensen, T. K. & Bødker, R. Identification of Dermacentor reticulatus ticks carrying Rickettsia raoultii on migrating jackal, Denmark. Emerg. Infect. Dis. 23, 2072–2074 (2017).
    Article  Google Scholar 

    74.
    Moutailler, S. et al. Co-infection of ticks: the rule rather than the exception. PLoS Negl. Trop. Dis. 10, e0004539 (2016).
    Article  CAS  Google Scholar 

    75.
    Reye, A. L. et al. Prevalence of tick-borne pathogens in Ixodes ricinus and Dermacentor reticulatus ticks from different geographical locations in Belarus. PLoS ONE 8, e54476 (2013).
    ADS  CAS  Article  Google Scholar 

    76.
    R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing https://www.r-project.org (2018).

    77.
    Cowling, D. W., Gardner, I. A. & Johnson, W. O. Comparison of methods for estimation of individual-level prevalence based on pooled samples. Prev. Vet. Med. 39, 211–225 (1999).
    CAS  Article  Google Scholar 

    78.
    ESRI. ArcGIS Desktop: Release 10.6.1. Redlands, CA: Environmental Systems Research Institute. (2017).

    79.
    Kulldorff M. and Information Management Services, I. SaTScanTM v9.6: Software for the spatial and space-time scan statistics www.satscan.org, 2018.

    80.
    Kleinman, K. rsatscan: Tools, classes, and methods for interfacing with SaTScan stand-alone software. (2015).

    81.
    Kulldorff, M. A spatial scan statistic. Communications in Statistics – Theory and Methods vol. 26 https://www.tandfonline.com/doi/abs/10.1080/03610929708831995 (1997).

    82.
    Han, J. et al. Using Gini coefficient to determining optimal cluster reporting sizes for spatial scan statistics. Int. J. Health Geogr. 15, 27 (2016).
    Article  Google Scholar 

    83.
    Kuhn., M., Contributions from Jed Wing, Steve Weston, Andre Williams, Chris Keefer, Allan Engelhardt, T., Cooper, Zachary Mayer, Brenton Kenkel, the R Core Team, Michael Benesty, Reynald Lescarbeau, A. Z. & Luca Scrucca, Yuan Tang, C. C. and T. H. caret: Classification and regression training. R package version 6.0-81. https://CRAN.R-project.org/package=caret. (2018).

    84.
    Elith, J., Leathwick, J. R. & Hastie, T. A working guide to boosted regression trees. J. Anim. Ecol. 77, 802–813 (2008).
    CAS  Article  Google Scholar 

    85.
    Vapnik, V., Golowich, S. E. & Smola, A. Support vector method for function approximation, regression estimation, and signal processing. in Advances in Neural Information Processing Systems 9 (eds. Mozer, M., Jordan, M. & Petsche, T.) 281–287 (MIT Press., 1997).

    86.
    Sanz, H., Valim, C., Vegas, E., Oller, J. M. & Reverter, F. SVM-RFE: selection and visualization of the most relevant features through non-linear kernels. BMC Bioinform. 19, 432 (2018).
    Article  Google Scholar 

    87.
    Ghojogh, B., Ca, B., Crowley, M. & Ca, M. The theory behind overfitting, cross validation, rRegularization, bagging, and boosting: tTutorial. https://arxiv.org/abs/1905.12787 [stat.ML] 1–23 (2019).

    88.
    Skarphédinsson, S., Jensen, P. M. & Kristiansen, K. Survey of tickborne infections in Denmark. Emerg. Infect. Dis. 11, 1055–1061 (2005).
    Article  Google Scholar 

    89.
    Quarsten, H., Skarpaas, T., Fajs, L., Noraas, S. & Kjelland, V. Tick-borne bacteria in Ixodes ricinus collected in southern Norway evaluated by a commercial kit and established real-time PCR protocols. Ticks Tick. Borne. Dis. 6, 538–544 (2015).
    CAS  Article  Google Scholar 

    90.
    Wilhelmsson, P. et al. Prevalence and diversity of Borrelia species in ticks that have bitten humans in Sweden. J. Clin. Microbiol. 48, 4169–4176 (2010).
    CAS  PubMed  PubMed Central  Google Scholar 

    91.
    Strnad, M., Hönig, V., Růžek, D., Grubhoffer, L. & Rego, R. O. M. Europe-wWide meta-analysis of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks. Appl. Environ. Microbiol. 83, 3838 (2017).
    Google Scholar 

    92.
    Mysterud, A. et al. Tick abundance, pathogen prevalence, and disease incidence in two contrasting regions at the northern distribution range of Europe. Parasit. Vectors 11, 309 (2018).
    PubMed  PubMed Central  Google Scholar 

    93.
    Severinsson, K., Jaenson, T. G., Pettersson, J., Falk, K. & Nilsson, K. Detection and prevalence of Anaplasma phagocytophilum and Rickettsia helvetica in Ixodes ricinus ticks in seven study areas in Sweden. Parasit. Vectors 3, 66 (2010).
    PubMed  PubMed Central  Google Scholar 

    94.
    Karlsson, M. E. & Andersson, M. O. Babesia species in questing Ixodes ricinus, Sweden. Ticks Tick. Borne. Dis. 7, 10–12 (2016).
    PubMed  Google Scholar 

    95.
    Øines, Ø., Radzijevskaja, J., Paulauskas, A. & Rosef, O. Prevalence and diversity of Babesia spp. in questing Ixodes ricinus ticks from Norway. Parasit. Vectors 5, 156 (2012).
    PubMed  PubMed Central  Google Scholar 

    96.
    Andersson, M., Bartkova, S., Lindestad, O. & Råberg, L. Co-Infection with ‘Candidatus Neoehrlichia mikurensis’ and Borrelia afzelii in Ixodes ricinus Ticks in Southern Sweden. Vector-Borne Zoonotic Dis. 13, 438–442 (2013).
    PubMed  Google Scholar 

    97.
    Pedersen, B. N. et al. Distribution of Neoehrlichia mikurensis in Ixodes ricinus ticks along the coast of Norway: the western seaboard is a low-prevalence region. Zoonoses Public Health https://doi.org/10.1111/zph.12662 (2019).
    Article  PubMed  Google Scholar 

    98.
    Kantsø, B., Bo Svendsen, C., Moestrup Jensen, P., Vennestrøm, J. & Krogfelt, K. A. Seasonal and habitat variation in the prevalence of Rickettsia helvetica in Ixodes ricinus ticks from Denmark. Ticks Tick. Borne. Dis. 1, 101–103 (2010).
    PubMed  Google Scholar 

    99.
    Solano-Gallego, L., Sainz, Á., Roura, X., Estrada-Peña, A. & Miró, G. A review of canine babesiosis: the European perspective. Parasit. Vectors 9, 336 (2016).
    Article  CAS  Google Scholar 

    100.
    Randolph, S. E. The shifting landscape of tick-borne zoonoses: tick-borne encephalitis and Lyme borreliosis in Europe. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 356, 1045–1056 (2001).
    CAS  PubMed  PubMed Central  Google Scholar 

    101.
    Sumilo, D. et al. Tick-borne encephalitis in the Baltic States : Identifying risk factors in space and time. Int. J. Med. Microbiol. 296(Suppl), 76–79 (2006).
    PubMed  Google Scholar 

    102.
    Sumilo, D. et al. Socio-economic factors in the differential upsurge of tick-borne encephalitis in Central and Eastern Europe. Rev. Med. Virol. 18, 81–95 (2008).
    PubMed  Google Scholar 

    103.
    Randolph, S. E., Green, R. M., Peacey, M. F. & Rogers, D. J. Seasonal synchrony : the key to tick-borne encephalitis foci identified by satellite data. Parasitology 121, 15–23 (2000).
    PubMed  Google Scholar 

    104.
    Halos, L. et al. Ecological factors characterizing the prevalence of bacterial tick-borne pathogens in Ixodes ricinus ticks in pastures and woodlands. Appl. Environ. Microbiol. 76, 4413–4420 (2010).
    CAS  PubMed  PubMed Central  Google Scholar 

    105.
    Sjörs, H. Nordisk växtgeografi (in Swedish) (Bonniers, Scandinavian University Books, 1967).
    Google Scholar  More

  • in

    A first look at the metabolic rate of Greenland sharks (Somniosus microcephalus) in the Canadian Arctic

    1.
    Cavicchioli, R., Amils, R., Wagner, D. & McGenity, T. Life and applications of extremophiles. Environ. Microbiol. 13, 1903–1907 (2011).
    Article  Google Scholar 
    2.
    Riesch, R., Tobler, M. & Plath, M. Extremophile Fishes (Springer, New York, 2015).
    Google Scholar 

    3.
    Wharton, D. A. Life at the Limits: Organisms in Extreme Environments (Cambridge University Press, Cambridge, 2007).
    Google Scholar 

    4.
    Lear, K. O. et al. Divergent field metabolic rates highlight the challenges of increasing temperatures and energy limitation in aquatic ectotherms. Oecologia 193, 311–323 (2020).
    ADS  Article  Google Scholar 

    5.
    Elliott, K. H. et al. High flight costs, but low dive costs, in auks support the biomechanical hypothesis for flightlessness in penguins. Proc. Natl. Acad. Sci. 110, 9380–9384 (2013).
    ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

    6.
    Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Toward a metabolic theory of ecology. Ecology 85, 1771–1789 (2004).
    Article  Google Scholar 

    7.
    Clarke, A. & Johnston, N. M. Scaling of metabolic rate with body mass and temperature in teleost fish. J. Anim. Ecol. 68, 893–905 (1999).
    Article  Google Scholar 

    8.
    Schulte, P. M. The effects of temperature on aerobic metabolism: towards a mechanistic understanding of the responses of ectotherms to a changing environment. J. Exp. Biol. 218, 1856–1866 (2015).
    Article  Google Scholar 

    9.
    Kleiber, M. Body size and metabolism. ENE 1, 315–353 (1932).
    Google Scholar 

    10.
    Glazier, D. S. A unifying explanation for diverse metabolic scaling in animals and plants. Biol. Rev. 85, 111–138 (2010).
    Article  PubMed  Google Scholar 

    11.
    Jerde, C. L. et al. Strong evidence for an intraspecific metabolic scaling coefficient near 0.89 in fish. Front. Physiol. 10, 1166 (2019).
    Article  PubMed  PubMed Central  Google Scholar 

    12.
    van der Meer, J. Metabolic theories in ecology. Trends Ecol. Evol. 21, 136–140 (2006).
    Article  PubMed  Google Scholar 

    13.
    Luongo, S. M. & Lowe, C. G. Seasonally acclimated metabolic Q10 of the California horn shark, Heterodontus francisci. J. Exp. Mar. Bio. Ecol. 503, 129–135 (2018).
    Article  Google Scholar 

    14.
    White, C. R., Alton, L. A. & Frappell, P. B. Metabolic cold adaptation in fishes occurs at the level of whole animal, mitochondria and enzyme. Proc. R. Soc. B Biol. Sci. 279, 1740–1747 (2011).
    Article  CAS  Google Scholar 

    15.
    Krogh, A. The Quantitative Relation Between Temperature and Standard Metabolism in Animals (Internationale Zeitschrift fuÈr Physikalisch-Chemische Biologie, New York, 1914).
    Google Scholar 

    16.
    Messamah, B., Kellermann, V., Malte, H., Loeschcke, V. & Overgaard, J. Metabolic cold adaptation contributes little to the interspecific variation in metabolic rates of 65 species of Drosophilidae. J. Insect Physiol. 98, 309–316 (2017).
    CAS  Article  PubMed  Google Scholar 

    17.
    Holeton, G. F. Metabolic cold adaptation of polar fish: fact or artefact?. Physiol. Zool. 47, 137–152 (1974).
    Article  Google Scholar 

    18.
    Steffensen, J. F. Metabolic cold adaptation of polar fish based on measurements of aerobic oxygen consumption: fact or artefact? Artefact!. Comp. Biochem. Physiol. A. 132, 789–795 (2002).
    Article  Google Scholar 

    19.
    Peck, L. S. A cold limit to adaptation in the sea. Trends Ecol. Evol. 31, 13–26 (2016).
    Article  PubMed  Google Scholar 

    20.
    Chabot, D., Steffensen, J. F. & Farrell, A. P. The determination of standard metabolic rate in fishes. J. Fish Biol. 88, 81–121 (2016).
    CAS  Article  PubMed  Google Scholar 

    21.
    Lawson, C. L. et al. Powering ocean giants : the energetics of shark and ray megafauna. Trends Ecol. Evol. 34, 1–13 (2019).
    MathSciNet  Article  Google Scholar 

    22.
    Lowe, C. Metabolic rates of juvenile scalloped hammerhead sharks (Sphyrna lewini). Mar. Biol. 139, 447–453 (2001).
    Article  Google Scholar 

    23.
    Payne, N. L. et al. A new method for resolving uncertainty of energy requirements in large water breathers: the ‘mega-flume’ seagoing swim-tunnel respirometer. Methods Ecol. Evol. 6, 668–677 (2015).
    Article  Google Scholar 

    24.
    Byrnes, E. E., Lear, K. O., Morgan, D. L. & Gleiss, A. C. Respirometer in a box: development and use of a portable field respirometer for estimating oxygen consumption of large-bodied fishes. J. Fish Biol. 96, 1045–1050 (2020).
    CAS  Article  PubMed  Google Scholar 

    25.
    MacNeil, M. A. et al. Biology of the greenland shark Somniosus microcephalus. J. Fish Biol. 80, 991–1018 (2012).
    CAS  Article  PubMed  Google Scholar 

    26.
    Edwards, J. E. et al. Advancing research for the management of long-lived species: a case study on the Greenland shark. Front. Mar. Sci. 6, 12 (2019).
    Article  Google Scholar 

    27.
    Augustine, S., Lika, K. & Kooijman, S. A. L. M. Comment on the ecophysiology of the Greenland shark, Somniosus microcephalus. Polar Biol. 40, 2429–2433 (2017).
    Article  Google Scholar 

    28.
    Nielsen, J. et al. Eye lens radiocarbon reveals centuries of longevity in the Greenland shark (Somniosus microcephalus). Science 353, 702–704 (2016).
    ADS  CAS  Article  PubMed  Google Scholar 

    29.
    Watanabe, Y. Y., Lydersen, C., Fisk, A. T. & Kovacs, K. M. The slowest fish: Swim speed and tail-beat frequency of Greenland sharks. J. Exp. Mar. Biol. Ecol. 426–427, 5–11 (2012).
    Article  Google Scholar 

    30.
    Hussey, N. E. et al. Rescaling the trophic structure of marine food webs. Ecol. Lett. 17, 239–250 (2014).
    Article  Google Scholar 

    31.
    Devine, B. M., Wheeland, L. J. & Fisher, J. A. D. First estimates of Greenland shark (Somniosus microcephalus) local abundances in Arctic waters. Sci. Rep. 8, 1–10 (2018).
    CAS  Article  Google Scholar 

    32.
    Wilson, E. E. & Wolkovich, E. M. Scavenging: how carnivores and carrion structure communities. Trends Ecol. Evol. 26, 129–135 (2011).
    Article  Google Scholar 

    33.
    Lear, K. O. et al. Correlations of metabolic rate and body acceleration in three species of coastal sharks under contrasting temperature regimes. J. Exp. Biol. 220, 397–407 (2017).
    Article  Google Scholar 

    34.
    Killen, S. S., Atkinson, D. & Glazier, D. S. The intraspecific scaling of metabolic rate with body mass in fishes depends on lifestyle and temperature. Ecol. Lett. 13, 184–193 (2010).
    Article  Google Scholar 

    35.
    Lear, K. O., Whitney, N. M., Brewster, L. R. & Gleiss, A. C. Treading water: respirometer choice may hamper comparative studies of energetics in fishes. Mar. Freshw. Res. 70, 437–448 (2018).
    Article  Google Scholar 

    36.
    Whitney, N. M., Lear, K. O., Gaskins, L. C. & Gleiss, A. C. The effects of temperature and swimming speed on the metabolic rate of the nurse shark (Ginglymostoma cirratum, Bonaterre). J. Exp. Mar. Bio. Ecol. 477, 40–46 (2016).
    Article  Google Scholar 

    37.
    Sims, D. W. The effect of body size on the standard metabolic rate of the lesser spotted dogfish. J. Fish Biol. 48, 542–544 (1996).
    Article  Google Scholar 

    38.
    Semmens, J. M., Payne, N. L., Huveneers, C., Sims, D. W. & Bruce, B. D. Feeding requirements of white sharks may be higher than originally thought. Sci. Rep. 3, 10–13 (2013).
    Article  CAS  Google Scholar 

    39.
    Giacomin, M., Schulte, P. M. & Wood, C. M. Differential effects of temperature on oxygen consumption and branchial fluxes of urea, ammonia, and water in the dogfish shark (Squalus acanthias suckleyi). Physiol. Biochem. Zool. 90, 627–637 (2017).
    Article  PubMed  PubMed Central  Google Scholar 

    40.
    Lowe, C. G. Bioenergetics of free-ranging juvenile scalloped hammerhead sharks (Sphyrna lewini) in Kāne’ohe Bay, Ō’ahu, HI. J. Exp. Mar. Biol. Ecol. 278, 141–156 (2002).
    Article  Google Scholar 

    41.
    Ezcurra, J. M., Lowe, C. G., Mollet, H. F., Ferry, L. A. & O’Sullivan, J. B. Oxygen consumption rate of young-of-the-year white sharks, Carcharodon carcharias during transport to the Monterey Bay Aquarium. Glob. Perspect. Biol. Life Hist. 1, 17–26 (2012).
    Article  Google Scholar 

    42.
    Barnett, A. et al. The utility of bioenergetics modelling in quantifying predation rates of marine apex predators: ecological and fisheries implications. Sci. Rep. 7, 12982 (2017).
    ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

    43.
    Watanabe, Y. Y., Payne, N. L., Semmens, J. M., Fox, A. & Huveneers, C. Swimming strategies and energetics of endothermic white sharks during foraging. J. Exp. Biol. 222, 4 (2019).
    Article  Google Scholar 

    44.
    Secor, S. M. Specific dynamic action: a review of the postprandial metabolic response. J. Comp. Physiol. B 179, 1–56 (2009).
    ADS  Article  PubMed  Google Scholar 

    45.
    Auer, S. K., Dick, C. A., Metcalfe, N. B. & Reznick, D. N. Metabolic rate evolves rapidly and in parallel with the pace of life history. Nat. Commun. 9, 8–13 (2018).
    ADS  Article  CAS  Google Scholar 

    46.
    Drazen, J. C. & Seibel, B. A. Depth-related trends in metabolism of benthic and benthopelagic deep-sea fishes. Limnol. Oceanogr. 52, 2306–2316 (2007).
    ADS  CAS  Article  Google Scholar 

    47.
    Brett, J. R. & Groves, T. D. D. Physiological energetics. Fish Physiol. 8, 280–352 (1979).
    Google Scholar 

    48.
    Widdows, J. Application of calorimetric methods in ecological studies. Therm. Energy. Stud. Cell. Biol. Syst. 1, 182–215 (1987).
    Article  Google Scholar 

    49.
    Armstrong, J. B. & Schindler, D. E. Excess digestive capacity in predators reflects a life of feast and famine. Nature 476, 84–87 (2011).
    CAS  Article  PubMed  Google Scholar 

    50.
    Stirling, I. & McEwan, E. Caloric value of whole ringed seals (Phoca hispida) in relation to Polar Bear (Ursus maritimus) ecology and hunting behavior. Can. J. Zool. 53, 1021–1027 (1975).
    CAS  Article  PubMed  Google Scholar 

    51.
    Furey, N. B., Hinch, S. G., Mesa, M. G. & Beauchamp, D. A. Piscivorous fish exhibit temperature-influenced binge feeding during an annual prey pulse. J. Anim. Ecol. 85, 1307–1317 (2016).
    Article  PubMed  Google Scholar 

    52.
    Svendsen, M. B. S., Bushnell, P. G. & Steffensen, J. F. Design and setup of intermittent-flow respirometry system for aquatic organisms. J. Fish Biol. 88, 26–50 (2016).
    CAS  Article  PubMed  Google Scholar 

    53.
    Clark, T. D., Sandblom, E. & Jutfelt, F. Aerobic scope measurements of fishes in an era of climate change: respirometry, relevance and recommendations. J. Exp. Biol. 216, 2771–2782 (2013).
    Article  PubMed  Google Scholar 

    54.
    Leclerc, L.-M.E. et al. A missing piece in the Arctic food web puzzle? Stomach contents of Greenland sharks sampled in Svalbard, Norway. Polar Biol. 35, 1197–1208 (2012).
    Article  Google Scholar  More

  • in

    The resilience of weed seedbank regulation by carabid beetles, at continental scales, to alternative prey

    1.
    Begg, G. S. et al. A functional overview of conservation biological control. Crop Prot. 97, 145–158 (2017).
    Article  Google Scholar 
    2.
    Shields, M. W. et al. History, current situation and challenges for conservation biological control. Biol. Control 131, 25–35 (2019).
    Article  Google Scholar 

    3.
    Petit, S., Boursault, A. & Bohan, D. A. Weed seed choice by carabid beetles (Coleoptera: Carabidae): Linking field measurements with laboratory diet assessments. Eur. J. Entomol. 111, 1–6 (2014).
    Article  Google Scholar 

    4.
    Saska, P., Honěk, A. & Martinková, Z. Preferences of carabid beetles (Coleoptera: Carabidae) for herbaceous seeds. Acta Zool. Acad. Sci. Hung. 65, 57–76 (2019).
    Article  Google Scholar 

    5.
    Honěk, A., Martinkova, Z., Saska, P. & Pekar, S. Size and taxonomic constraints determine the seed preferences of Carabidae (Coleoptera). Basic Appl. Ecol. 8, 343–353 (2007).
    Article  Google Scholar 

    6.
    Honěk, A., Martinkova, Z. & Jarosik, V. Ground beetles (Carabidae) as seed predators. Eur. J. Entomol. 100, 531–544 (2003).
    Article  Google Scholar 

    7.
    Kulkarni, S. S., Dosdall, L. M. & Willenborg, C. J. The role of ground beetles (Coleoptera: Carabidae) in weed seed consumption: A review. Weed Sci. 63, 355–376 (2015).
    Article  Google Scholar 

    8.
    Petit, S., Trichard, A., Biju-Duval, L., McLaughlin, B. & Bohan, D. A. Interactions between conservation agricultural practice and landscape composition promote weed seed predation by invertebrates. Agric. Ecosyst. Environ. 240, 45–53 (2017).
    Article  Google Scholar 

    9.
    Kromp, B. Carabid beetles in sustainable agriculture: A review on pest control efficacy, cultivation impacts and enhancement. Agric. Ecosyst. Environ. 74, 187–228 (1999).
    Article  Google Scholar 

    10.
    Firbank, L. G. & Watkinson, A. R. On the analysis of competition within two-species mixtures of plants. J. Appl. Ecol. 22, 503–517 (1985).
    Article  Google Scholar 

    11.
    Westerman, P. R. et al. Are many little hammers effective? Velvetleaf (Abutilon theophrasti) population dynamics in two- and four-year crop rotation systems. Weed Sci. 53, 382–392 (2005).
    CAS  Article  Google Scholar 

    12.
    Petit, S. et al. Biodiversity-based options for arable weed management. A review. Agron. Sustain. Dev. 38, 48 (2018).
    Article  Google Scholar 

    13.
    Westerman, P. R., Dixon, P. M. & Liebman, M. Burial rates of surrogate seeds in arable fields. Weed Res. 49, 142–152 (2009).
    Article  Google Scholar 

    14.
    Trichard, A., Ricci, B., Ducourtieux, C. & Petit, S. The spatio-temporal distribution of weed seed predation differs between conservation agriculture and conventional tillage. Agric. Ecosyst. Environ. 188, 40–47 (2014).
    Article  Google Scholar 

    15.
    Carbonne, B., Bohan, D. A. & Petit, S. Key carabid species drive spring weed seed predation of Viola arvensis. Biol. Control 141, 104148 (2020).
    CAS  Article  Google Scholar 

    16.
    Westerman, P. R., Wes, J. S., Kropff, M. J. & Van Der Werf, W. Annual losses of weed seeds due to predation in organic cereal fields. J. Appl. Ecol. 40, 824–836 (2003).
    Article  Google Scholar 

    17.
    Blubaugh, C. K. & Kaplan, I. Invertebrate seed predators reduce weed emergence following seed rain. Weed Sci. 64, 80–86 (2016).
    Article  Google Scholar 

    18.
    Pannwitt, H., Westerman, P. R. & Gerowitt, B. Post-dispersal seed predation can limit the number of seedlings of Echinochloa crus-galli. Biol. Control 143, 95–98 (2019).
    Google Scholar 

    19.
    Bohan, D. A., Boursault, A., Brooks, D. R. & Petit, S. National-scale regulation of the weed seedbank by carabid predators. J. Appl. Ecol. 48, 888–898 (2011).
    Article  Google Scholar 

    20.
    Saska, P., Van Der Werf, W., De Vries, E. & Westerman, P. R. Spatial and temporal patterns of carabid activity-density in cereals do not explain levels of predation on weed seeds. Bull. Entomol. Res. 98, 169–181 (2008).
    CAS  Article  PubMed  Google Scholar 

    21.
    Mauchline, A. L., Watson, S. J., Brown, V. K. & Froud-Williams, R. J. Post-dispersal seed predation of non-target weeds in arable crops. Weed Res. 45, 157–164 (2005).
    Article  Google Scholar 

    22.
    Davis, A. S. & Raghu, S. Weighing abiotic and biotic influences on weed seed predation. Weed Res. 50, 402–412 (2010).
    Article  Google Scholar 

    23.
    Davis, A. S., Taylor, E. C., Haramoto, E. R. & Renner, K. A. Annual postdispersal weed seed predation in contrasting field environments. Weed Sci. 61, 296–302 (2013).
    CAS  Article  Google Scholar 

    24.
    Lövei, G. L. & Szentkiralyi, F. Carabids climbing maize plants. Z. Angew. Entomol. 97, 107–110 (1984).
    Article  Google Scholar 

    25.
    Frei, B., Guenay, Y., Bohan, D. A., Traugott, M. & Wallinger, C. Molecular analysis indicates high levels of carabid weed seed consumption in cereal fields across Central Europe. J. Pest Sci. 2004(92), 935–942 (2019).
    Article  Google Scholar 

    26.
    Roubinet, E. et al. High redundancy as well as complementary prey choice characterize generalist predator food webs in agroecosystems. Sci. Rep. 8, 8054 (2018).
    ADS  Article  CAS  PubMed  PubMed Central  Google Scholar 

    27.
    Staudacher, K. et al. Habitat heterogeneity induces rapid changes in the feeding behaviour of generalist arthropod predators. Funct. Ecol. 32, 809–819 (2018).
    Article  PubMed  PubMed Central  Google Scholar 

    28.
    Evans, E. W. Multitrophic interactions among plants, aphids, alternate prey and shared natural enemies – A review. Eur. J. Entomol. 105, 369–380 (2008).
    Article  Google Scholar 

    29.
    Snyder, W. E. Give predators a complement: Conserving natural enemy biodiversity to improve biocontrol. Biol. Control 135, 73–82 (2019).
    Article  Google Scholar 

    30.
    Harwood, J. D. et al. Invertebrate biodiversity affects predator fitness and hence potential to control pests in crops. Biol. Control 51, 499–506 (2009).
    Article  Google Scholar 

    31.
    Chailleux, A., Mohl, E. K., Teixeira Alves, M., Messelink, G. J. & Desneux, N. Natural enemy-mediated indirect interactions among prey species: Potential for enhancing biocontrol services in agroecosystems. Pest Manag. Sci. 70, 1769–1779 (2014).

    32.
    von Berg, K., Thies, C., Tscharntke, T. & Scheu, S. Cereal aphid control by generalist predators in presence of belowground alternative prey: Complementary predation as affected by prey density. Pedobiologia (Jena). 53, 41–48 (2009).
    Article  Google Scholar 

    33.
    Mair, J. & Port, G. R. Predation by the carabid beetles Pterostichus madidus and Nebria brevicollis is affected by size and condition of the prey slug Deroceras reticulatum. Agric. For. Entomol. 3, 99–106 (2001).
    Article  Google Scholar 

    34.
    Symondson, W. O. C. et al. Biodiversity vs. biocontrol: positive and negative effects of alternative prey on control of slugs by carabid beetles. Bull. Entomol. Res. 96, 637–645 (2006).

    35.
    Prasad, R. P. & Snyder, W. E. Polyphagy complicates conservation biological control that targets generalist predators. J. Appl. Ecol. 43, 343–352 (2006).
    Article  Google Scholar 

    36.
    Renkema, J. M., Lynch, D. H., Cutler, G. C., MacKenzie, K. & Walde, S. J. Predation by Pterostichus melanarius (Illiger) (Coleoptera: Carabidae) on immature Rhagoletis mendax Curran (Diptera: Tephritidae) in semi-field and field conditions. Biol. Control 60, 46–53 (2012).
    Article  Google Scholar 

    37.
    Roubinet, E. et al. Diet of generalist predators reflects effects of cropping period and farming system on extra- and intraguild prey. Ecol. Appl. 27, 1167–1177 (2017).
    Article  Google Scholar 

    38.
    Honěk, A., Saska, P. & Martinkova, Z. Seasonal variation in seed predation by adult carabid beetles. Entomol. Exp. Appl. 118, 157–162 (2006).
    Article  Google Scholar 

    39.
    Talarico, F., Giglio, A., Pizzolotto, R. & Brandmayr, P. A synthesis of feeding habits and reproduction rhythm in Italian seed-feeding ground beetles (Coleoptera: Carabidae). Eur. J. Entomol. 113, 325–336 (2016).
    Article  Google Scholar 

    40.
    Charalabidis, A., Dechaume-Moncharmont, F.-X., Carbonne, B., Bohan, D. A. & Petit, S. Diversity of foraging strategies and responses to predator interference in seed-eating carabid beetles. Basic Appl. Ecol. 36, 13–24 (2019).
    Article  Google Scholar 

    41.
    Pilipaviius, V. Weed seed rain dynamics and ecological control ability in agrophytocenosis. in Herbicides—Advances in Research (ed. Price, A.) 51–83 (InTech, 2013). https://doi.org/10.5772/55972.

    42.
    Saska, P., Koprdová, S., Martinková, Z. & Honěk, A. Comparing methods of weed seed exposure to predators. Ann. Appl. Biol. 164, 301–312 (2014).
    Article  Google Scholar 

    43.
    Johnson, N. E. & Cameron, R. S. Phytophagous ground beetles. Ann. Entomol. Soc. Am. 62, 909–914 (1969).
    Article  Google Scholar 

    44.
    Russell, M. C., Lambrinos, J., Records, E. & Ellen, G. Seasonal shifts in ground beetle (Coleoptera: Carabidae) species and functional composition maintain prey consumption in Western Oregon agricultural landscapes. Biol. Control 106, 54–63 (2017).
    Article  Google Scholar 

    45.
    Williams, C. L. et al. Over-winter predation of Abutilon theophrasti and Setaria faberi seeds in arable land. Weed Res. 49, 439–447 (2009).
    Article  Google Scholar 

    46.
    Westerman, P., Luijendijk, C. D., Wevers, J. D. A. & Van Der Werf, W. Weed seed predation in a phenologically late crop. Weed Res. 51, 157–164 (2011).
    Article  Google Scholar 

    47.
    Winder, L. et al. Predatory activity and spatial pattern: The response of generalist carabids to their aphid prey. J. Anim. Ecol. 74, 443–454 (2005).
    Article  Google Scholar 

    48.
    Bohan, D. A. et al. Spatial dynamics of predation by carabid beetles on slugs. J. Anim. Ecol. 69, 367–379 (2000).
    Article  Google Scholar 

    49.
    Frank, S. D., Shrewsbury, P. M. & Denno, R. F. Plant versus prey resources: Influence on omnivore behavior and herbivore suppression. Biol. Control 57, 229–235 (2011).
    Article  Google Scholar 

    50.
    Abrams, P. A. & Matsuda, H. Positive indirect effects between prey species that share predators. Ecology 77, 610–616 (1996).
    Article  Google Scholar 

    51.
    Boetzl, F. A., Konle, A. & Krauss, J. Aphid cards – Useful model for assessing predation rates or bias prone nonsense?. J. Appl. Entomol. 144, 74–80 (2020).
    Article  Google Scholar 

    52.
    Bilde, T. & Toft, S. Consumption by carabid beetles of three cereal aphid species relative to other prey types. Entomophaga 42, 21–32 (1997).
    Article  Google Scholar 

    53.
    Madsen, M., Terkildsen, S. & Toft, S. Microcosm studies on control of aphids by generalist arthropod predators: Effects of alternative prey. Biocontrol 49, 483–504 (2004).
    Article  Google Scholar 

    54.
    Fawki, S. & Toft, S. Food preferences and the value of animal food for the carabid beetle Amara similata (Gyll.) (Col., Carabidae). J. Appl. Entomol. 129, 551–556 (2005).

    55.
    Saska, P. Effect of diet on the fecundity of three carabid beetles. Physiol. Entomol. 33, 188–192 (2008).
    ADS  Article  Google Scholar 

    56.
    Haschek, C., Drapela, T., Schuller, N., Fiedler, K. & Frank, T. Carabid beetle condition, reproduction and density in winter oilseed rape affected by field and landscape parameters. J. Appl. Entomol. 136, 665–674 (2012).
    Article  Google Scholar 

    57.
    Symondson, W. O. C., Sunderland, K. D. & Greenstone, M. H. Can generalist predators be effective biocontrol agents?. Annu. Rev. Entomol. 47, 561–594 (2002).
    CAS  Article  PubMed  Google Scholar 

    58.
    Lundgren, J. G. Chapter 18: Biological control of weed seeds in agriculture using omnivorous insects. in Relationships of Natural Enemies and Non-Prey Foods 333–351 (Springer Netherlands, 2009).

    59.
    Löbl, I. & Smetana, A. Catalogue of Palaearctic Colcoptera. Vol. 1 (2003).

    60.
    Homburg, K., Homburg, N., Schäfer, F., Schuldt, A. & Assmann, T. Carabids.org—A dynamic online database of ground beetle species traits (Coleoptera, Carabidae). Insect Conserv. Divers. 7, 195–205 (2014).

    61.
    Penell, A., Raub, F. & Höfer, H. Estimating biomass from body size of European spiders based on regression models. J. Arachnol. 46, 413 (2018).
    Article  Google Scholar 

    62.
    Pey, B. et al. A thesaurus for soil invertebrate trait-based approaches. PLoS ONE 9, e108985 (2014).
    ADS  Article  CAS  PubMed  PubMed Central  Google Scholar 

    63.
    Nentwig, W., Blick, T., Gloor, D., Hänggi, A. & Kropf, C. Araneae: Spiders of Europe. https://araneae.nmbe.ch, https://www.araneae.nmbe.ch (2019).

    64.
    Caballero, M., Baquero, E., Ariño, A. H. & Jordana, R. Indirect biomass estimations in Collembola. Pedobiologia (Jena). 48, 551–557 (2004).
    Article  Google Scholar 

    65.
    Migui, S. M. & Lamb, R. J. Sources of variation in the interaction between three cereal aphids (Hemiptera: Aphididae) and wheat (Poaceae). Bull. Entomol. Res. 96, 235–241 (2006).
    CAS  Article  PubMed  Google Scholar 

    66.
    Brooks, D. R. et al. Invertebrate responses to the management of genetically modified herbicide-tolerant and conventional spring crops. I. Soil-surface-active invertebrates. Philos. Trans. R. Soc. B Biol. Sci. 358, 1847–1862 (2003).

    67.
    Bohan, D. A. et al. Effects on weed and invertebrate abundance and diversity of herbicide management in genetically modified herbicide-tolerant winter-sown oilseed rape. Proc. R. Soc. B Biol. Sci. 272, 463–474 (2005).
    Article  Google Scholar 

    68.
    John, F. & Weisberg, S. An R Companion to Applied Regression. (Sage, 2019).

    69.
    Long, J. jtools: Analysis and Presentation of Social Scientific Data. R package version 2.0.1. (2019).

    70.
    Lenth, R. emmeans: Estimated Marginal Means, aka Least-Squares Means. (2020).

    71.
    Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
    Article  Google Scholar 

    72.
    R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. https://www.R-project.org (2019). More

  • in

    Life history and population regulation shape demographic competence and influence the maintenance of endemic disease

    1.
    Daszak, P., Cunningham, A. A. & Hyatt, A. D. Emerging infectious diseases of wildlife—threats to biodiversity and human health. Science 287, 443–449 (2000).
    CAS  PubMed  Google Scholar 
    2.
    Wiethoelter, A. K., Beltrán-Alcrudo, D., Kock, R. & Mor, S. M. Global trends in infectious diseases at the wildlife–livestock interface. Proc. Natl Acad. Sci. USA 112, 9662–9667 (2015).
    CAS  PubMed  Google Scholar 

    3.
    Jones, K. E. et al. Global trends in emerging infectious diseases. Nature 451, 990–993 (2008).
    CAS  PubMed  PubMed Central  Google Scholar 

    4.
    Narrod, C., Zinsstag, J. & Tiongco, M. A One Health framework for estimating the economic costs of zoonotic diseases on society. Ecohealth 9, 150–162 (2012).
    PubMed  PubMed Central  Google Scholar 

    5.
    Cunningham, A. A., Daszak, P. & Wood, J. L. N. One Health, emerging infectious diseases and wildlife: two decades of progress? Phil. Trans. R. Soc. B 372, 20160167 (2017).
    PubMed  Google Scholar 

    6.
    Park, A. W. Phylogenetic aggregation increases zoonotic potential of mammalian viruses. Biol. Lett. 15, 20190668 (2019).
    PubMed  Google Scholar 

    7.
    Davies, T. J. & Pedersen, A. B. Phylogeny and geography predict pathogen community similarity in wild primates and humans. Proc. R. Soc. B 275, 1695–1701 (2008).
    PubMed  Google Scholar 

    8.
    Dallas, T., Park, A. W. & Drake, J. M. Predictability of helminth parasite host range using information on geography, host traits and parasite community structure. Parasitology 144, 200–205 (2017).
    PubMed  Google Scholar 

    9.
    Martin, L. B., Burgan, S. C., Adelman, J. S. & Gervasi, S. S. Host competence: an organismal trait to integrate immunology and epidemiology. Integr. Comp. Biol. 56, 1225–1237 (2016).
    PubMed  Google Scholar 

    10.
    Olival, K. J. et al. Host and viral traits predict zoonotic spillover from mammals. Nature 546, 646–650 (2017).
    PubMed  PubMed Central  Google Scholar 

    11.
    Faust, C. L. et al. Pathogen spillover during land conversion. Ecol. Lett. 21, 471–483 (2018).
    PubMed  Google Scholar 

    12.
    Martin, G. et al. Climate change could increase the geographic extent of Hendra virus spillover risk. Ecohealth 15, 509–525 (2018).
    PubMed  PubMed Central  Google Scholar 

    13.
    Plowright, R. K. et al. Pathways to zoonotic spillover. Nat. Rev. Microbiol. 15, 502–510 (2017).
    CAS  PubMed  PubMed Central  Google Scholar 

    14.
    Viana, M. et al. Assembling evidence for identifying reservoirs of infection. Trends Ecol. Evol. 29, 270–279 (2014).
    PubMed  PubMed Central  Google Scholar 

    15.
    Thrall, P. H., Antonovics, J. & Hall, D. W. Host and pathogen coexistence in sexually transmitted and vector-borne diseases characterized by frequency-dependent disease transmission. Am. Nat. 142, 543–552 (1993).
    Google Scholar 

    16.
    Anderson, R. M. & May, R. M. The population dynamics of microparasites and their invertebrate hosts. Phil. Trans. R. Soc. Lond. B 291, 451–524 (1981).
    Google Scholar 

    17.
    Stearns, S. C. Life-history tactics: a review of the ideas. Q. Rev. Biol. 51, 3–47 (1976).
    CAS  PubMed  Google Scholar 

    18.
    Solomon, M. E. The natural control of animal populations. J. Anim. Ecol. 18, 1–35 (1949).

    19.
    McDonald, J. L. et al. Demographic buffering and compensatory recruitment promotes the persistence of disease in a wildlife population. Ecol. Lett. 19, 443–449 (2016).

    20.
    Promislow, D. E. L. & Harvey, P. H. Living fast and dying young: a comparative analysis of life-history variation among mammals. J. Zool. 220, 417–437 (1990).
    Google Scholar 

    21.
    Pfister, C. A. Patterns of variance in stage-structured populations: evolutionary predictions and ecological implications. Proc. Natl Acad. Sci. USA 95, 213–218 (1998).
    CAS  PubMed  Google Scholar 

    22.
    Johnson, P. T. J. et al. Living fast and dying of infection: host life history drives interspecific variation in infection and disease risk. Ecol. Lett. 15, 235–242 (2012).
    PubMed  Google Scholar 

    23.
    Cronin, J. P., Rúa, M. A. & Mitchell, C. E. Why is living fast dangerous? Disentangling the roles of resistance and tolerance of disease. Am. Nat. 184, 172–187 (2014).
    PubMed  Google Scholar 

    24.
    Lachish, S., McCallum, H. & Jones, M. Demography, disease and the devil: life-history changes in a disease-affected population of Tasmanian devils (Sarcophilus harrisii). J. Anim. Ecol. 78, 427–436 (2009).
    PubMed  Google Scholar 

    25.
    Muths, E., Scherer, R. D. & Pilliod, D. S. Compensatory effects of recruitment and survival when amphibian populations are perturbed by disease. J. Appl. Ecol. 48, 873–879 (2011).
    Google Scholar 

    26.
    Arthur, A., Ramsey, D. & Efford, M. Impact of bovine tuberculosis on a population of brushtail possums (Trichosurus vulpecula Kerr) in the Orongorongo Valley, New Zealand. Wildl. Res. 31, 389–395 (2004).
    Google Scholar 

    27.
    Gaillard, J.-M. & Yoccoz, N. G. Temporal variation in survival of mammals: a case of environmental canalization? Ecology 84, 3294–3306 (2003).
    Google Scholar 

    28.
    Korpimäki, E., Brown, P. R., Jacob, J. & Pech, R. P. The puzzles of population cycles and outbreaks of small mammals solved? BioScience 54, 1071–1079 (2004).
    Google Scholar 

    29.
    Gaillard, J.-M., Festa-Bianchet, M. & Yoccoz, N. G. Population dynamics of large herbivores: variable recruitment with constant adult survival. Trends Ecol. Evol. 13, 58–63 (1998).
    CAS  PubMed  Google Scholar 

    30.
    McCallum, H., Barlow, N. & Hone, J. How should pathogen transmission be modelled? Trends Ecol. Evol. 16, 295–300 (2001).
    CAS  PubMed  Google Scholar 

    31.
    De Castro, F. & Bolker, B. Mechanisms of disease‐induced extinction. Ecol. Lett. 8, 117–126 (2005).
    Google Scholar 

    32.
    McCallum, H. Disease and the dynamics of extinction. Phil. Trans. R. Soc. B 367, 2828–2839 (2012).
    PubMed  Google Scholar 

    33.
    Anderson, R. M. & May, R. M. Coevolution of hosts and parasites. Parasitology 85, 411–426 (1982).
    PubMed  Google Scholar 

    34.
    Alizon, S., Hurford, A., Mideo, N. & Van Baalen, M. Virulence evolution and the trade-off hypothesis: history, current state of affairs and the future. J. Evol. Biol. 22, 245–259 (2009).
    CAS  PubMed  Google Scholar 

    35.
    De Roode, J. C., Yates, A. J. & Altizer, S. Virulence-transmission trade-offs and population divergence in virulence in a naturally occurring butterfly parasite. Proc. Natl Acad. Sci. USA 105, 7489–7494 (2008).
    PubMed  Google Scholar 

    36.
    Anderson, R. M. Parasite pathogenicity and the depression of host population equilibria. Nature 279, 150 (1979).
    Google Scholar 

    37.
    Anderson, R. M. & May, R. M. Infectious Diseases of Humans: Dynamics and Control (Oxford Univ. Press, 1991).

    38.
    Lloyd-Smith, J. O. et al. Should we expect population thresholds for wildlife disease? Trends Ecol. Evol. 20, 511–519 (2005).
    PubMed  Google Scholar 

    39.
    Calvete, C. Modeling the effect of population dynamics on the impact of rabbit hemorrhagic disease. Conserv. Biol. 20, 1232–1241 (2006).
    PubMed  Google Scholar 

    40.
    Gutiérrez, J. S., Piersma, T. & Thieltges, D. W. Micro- and macroparasite species richness in birds: the role of host life history and ecology. J. Anim. Ecol. 88, 1226–1239 (2019).
    PubMed  Google Scholar 

    41.
    Calisher, C. H. et al. Do unusual site-specific population dynamics of rodent reservoirs provide clues to the natural history of hantaviruses? J. Wildl. Dis. 37, 280–288 (2001).
    CAS  PubMed  Google Scholar 

    42.
    Rickman, S. J., Dulvy, N. K., Jennings, S. & Reynolds, J. D. Recruitment variation related to fecundity in marine fishes. Can. J. Fish. Aquat. Sci. 57, 116–124 (2000).
    Google Scholar 

    43.
    Sæther, B.-E. & Bakke, Ø. Avian life history variation and contribution of demographic traits to the population growth rate. Ecology 81, 642–653 (2000).
    Google Scholar 

    44.
    Bielby, J. et al. The fast–slow continuum in mammalian life history: an empirical reevaluation. Am. Nat. 169, 748–757 (2007).
    CAS  PubMed  Google Scholar 

    45.
    Godfray, H. C. J. et al. A restatement of the natural science evidence base relevant to the control of bovine tuberculosis in Great Britain. Proc. R. Soc. Lond. B 280, 20131634 (2013).
    Google Scholar 

    46.
    Gaillard, J.-M. et al. Generation time: a reliable metric to measure life-history variation among mammalian populations. Am. Nat. 166, 119–123 (2005).
    PubMed  Google Scholar 

    47.
    Earn, D. J. D., Rohani, P. & Grenfell, B. T. Persistence, chaos and synchrony in ecology and epidemiology. Proc. R. Soc. Lond. B 265, 7–10 (1998).
    CAS  Google Scholar 

    48.
    Rohani, P., Earn, D. J. D. & Grenfell, B. T. Opposite patterns of synchrony in sympatric disease metapopulations. Science 286, 968–971 (1999).
    CAS  PubMed  Google Scholar 

    49.
    Salathé, M. & Jones, J. H. Dynamics and control of diseases in networks with community structure. PLoS Comput. Biol. 6, e1000736 (2010).
    PubMed  PubMed Central  Google Scholar 

    50.
    Sah, P., Leu, S. T., Cross, P. C., Hudson, P. J. & Bansal, S. Unraveling the disease consequences and mechanisms of modular structure in animal social networks. Proc. Natl Acad. Sci. USA 114, 4165–4170 (2017).

    51.
    Silk, M. J. et al. Integrating social behaviour, demography and disease dynamics in network models: applications to disease management in declining wildlife populations. Phil. Trans. R. Soc. B 374, 20180211 (2019).
    PubMed  Google Scholar 

    52.
    Hopkins, S. R., Fleming‐Davies, A. E., Belden, L. K. & Wojdak, J. M. Systematic review of modelling assumptions and empirical evidence: does parasite transmission increase nonlinearly with host density? Methods Ecol. Evol. 11, 476–486 (2020).
    Google Scholar 

    53.
    Froissart, R., Doumayrou, J., Vuillaume, F., Alizon, S. & Michalakis, Y. The virulence–transmission trade-off in vector-borne plant viruses: a review of (non-) existing studies. Phil. Trans. R. Soc. B 365, 1907–1918 (2010).
    CAS  PubMed  Google Scholar 

    54.
    Wickham, M. E., Brown, N. F., Boyle, E. C., Coombes, B. K. & Finlay, B. B. Virulence is positively selected by transmission success between mammalian hosts. Curr. Biol. 17, 783–788 (2007).
    CAS  PubMed  Google Scholar 

    55.
    Paul, R. E. L. et al. Experimental evaluation of the relationship between lethal or non-lethal virulence and transmission success in malaria parasite infections. BMC Evol. Biol. 4, 30 (2004).
    CAS  PubMed  PubMed Central  Google Scholar 

    56.
    Agnew, P. & Koella, J. C. Virulence, parasite mode of transmission, and host fluctuating asymmetry. Proc. R. Soc. Lond. B 264, 9–15 (1997).
    CAS  Google Scholar 

    57.
    Medica, D. L. & Sukhdeo, M. V. K. Estimating transmission potential in gastrointestinal nematodes (Order: Strongylida). J. Parasitol. 87, 442–446 (2001).
    CAS  PubMed  Google Scholar 

    58.
    Jäkel, T. et al. Reduction of transmission stages concomitant with increased host immune responses to hypervirulent Sarcocystis singaporensis, and natural selection for intermediate virulence. Int. J. Parasitol. 31, 1639–1647 (2001).
    PubMed  Google Scholar 

    59.
    Stott, I., Hodgson, D. J. & Townley, S. popdemo: an R package for population demography using projection matrix analysis. Methods Ecol. Evol. 3, 797–802 (2012).
    Google Scholar 

    60.
    R Core Team R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019). More

  • in

    The politics of biodiversity offsetting across time and institutional scales

    1.
    Bull, J. W., Gordon, A., Watson, J. E. M. & Maron, M. Seeking convergence on the key concepts in ‘no net loss’ policy. J. Appl. Ecol. 53, 1686–1693 (2016).
    Article  Google Scholar 
    2.
    Maron, M. et al. The many meanings of no net loss in environmental policy. Nat. Sustain 1, 19–27 (2018).
    Article  Google Scholar 

    3.
    Bull, J. W. & Strange, N. The global extent of biodiversity offset implementation under no net loss policies. Nat. Sustain 1, 790–798 (2018).
    Article  Google Scholar 

    4.
    von Hase, A. & ten Kate, K. Correct framing of biodiversity offsets and conservation: a response to Apostolopoulou & Adams. Oryx 51, 32–34 (2017).
    Article  Google Scholar 

    5.
    Standard on Biodiversity Offsets (BBOP, Forest Trends, 2012).

    6.
    Zu Ermgassen, S. O. S. E. et al. The ecological outcomes of biodiversity offsets under “no net loss” policies: a global review. Conserv. Lett. https://doi.org/10.1111/conl.12664 (2019).

    7.
    Moreno-Mateos, D., Maris, V., Béchet, A. & Curran, M. The true loss caused by biodiversity offsets. Biol. Conserv. 192, 552–559 (2015).
    Article  Google Scholar 

    8.
    Bonneuil, C. Tell me where you come from, I will tell you who you are: a genealogy of biodiversity offsetting mechanisms in historical context. Biol. Conserv. 192, 485–491 (2015).
    Article  Google Scholar 

    9.
    Boon, P. I. & Prahalad, V. Ecologists, economics and politics: problems and contradictions in applying neoliberal ideology to nature conservation in Australia. Pac. Conserv. Biol. 23, 115–132 (2017).
    Article  Google Scholar 

    10.
    Penca, J. Marketing the market: the ideology of market mechanisms for biodiversity conservation. Transnatl Environ. Law 2, 235–257 (2013).
    Article  Google Scholar 

    11.
    Lapeyre, R., Froger, G. & Hrabanski, M. Biodiversity offsets as market-based instruments for ecosystem services? From discourses to practices. Ecosyst. Serv. 15, 125–133 (2015).
    Article  Google Scholar 

    12.
    Hackett, R. Market-based environmental governance and public resources in Alberta, Canada. Ecosyst. Serv. 15, 174–180 (2015).
    Article  Google Scholar 

    13.
    Zero Draft of the Post-2020 Global Biodiversity Framework (CBD, 2020).

    14.
    Feindt, P. H. & Oels, A. Does discourse matter? Discourse analysis in environmental policy making. J. Environ. Policy Plan. 7, 161–173 (2005).
    Article  Google Scholar 

    15.
    Leipold, S., Feindt, P. H., Winkel, G. & Keller, R. Discourse analysis of environmental policy revisited: traditions, trends, perspectives. J. Environ. Policy Plan. 21, 445–463 (2019).
    Article  Google Scholar 

    16.
    Hajer, M. & Versteeg, W. A decade of discourse analysis of environmental politics: achievements, challenges, perspectives. J. Environ. Policy Plan. 7, 175–184 (2005).
    Article  Google Scholar 

    17.
    Bacchi, C. & Goodwin, S. Poststructural Policy Analysis: A Guide to Practice (Springer, 2016).

    18.
    Dryzek, J. S. The Politics of the Earth: Environmental Discourses (Oxford Univ. Press, 2013).

    19.
    Foucault, M. The History of Sexuality: An Introduction Vol. 1(Penguin Group, 2008).

    20.
    Hajer, M. A. in Words Matter in Policy and Planning: Discource Theory and Method in the Social Sciences (eds Van den Brink, M. & Metze, T.) 65–76 (Netherlands Graduate School of Urban and Regional Research, 2006).

    21.
    Hajer, M. A. The Politics of Environmental Discourse: Ecological Modernization and the Policy Process (Clarendon Press, 1995).

    22.
    Hopwood, B., Mellor, M. & O’Brien, G. Sustainable development: mapping different approaches. Sustain. Dev. 13, 38–52 (2005).
    Article  Google Scholar 

    23.
    Carson, R. Silent Spring (Houghton Mifflin, 1962).

    24.
    Our Common Future (United Nations World Commission on Environment and Development, Oxford Univ. Press, 1987).

    25.
    Robertson, M. M. The neoliberalization of ecosystem services: wetland mitigation banking and problems in environmental governance. Geoforum 35, 361–373 (2004).
    Article  Google Scholar 

    26.
    Clapp, J. & Dauvergne, P. in Paths to a Green World: The Political Economy of the Global Environment (eds Clapp, J. & Dauvergne, P.) 161–191 (MIT Press, 2011).

    27.
    Werksman, J. The clean development mechanism: unwrapping the Kyoto surprise. Rev. Eur. Comp. Int. Environ. Law 7, 147–158 (1998).
    Article  Google Scholar 

    28.
    Christoff, P. Ecological modernisation, ecological modernities. Environ. Polit. 5, 476–500 (1996).
    Article  Google Scholar 

    29.
    Breaking New Ground: The Report of the Mining, Minerals and Sustainable Development Project (International Institute for Environment and Development, World Business Council for Sustainable Development, Earthscan Publications, 2002).

    30.
    Hrabanski, M. The biodiversity offsets as market-based instruments in global governance: origins, success and controversies. Ecosyst. Serv. 15, 143–151 (2015).
    Article  Google Scholar 

    31.
    Ecosystems and Human Well-Being: Our Human Planet: Summary for Decision Makers (Millennium Ecosystem Assessment, 2005); http://millenniumassessment.org/en/Global.html

    32.
    The Economics of Ecosystems and Biodiversity: Mainstreaming the Economics of Nature: A Synthesis of the Approach, Conclusions and Recommendations of TEEB (TEEB, 2010); http://teebweb.org/publications/teeb-for/synthesis/

    33.
    Bassey, N. et al. IUCN Withdrawal (Friends of the Earth International, 2009).

    34.
    WWC 10 Final Resolution 12: Building a Global Alliance to Assert ‘No-Go Areas’ for Mining and Other Extractive Industries and Destructive Activities Threatening World Heritage Sites, and Protected Areas, including Indigenous Peoples’ and Local Communities Conserved Areas and Territories (ICCAs) and Sacred Natural Sites and Territories (WWC, 2013).

    35.
    WCC-2012-Res-110-EN: Biodiversity Offsets and Related Compensatory Approaches (WCC, 2012).

    36.
    IUCN Resolutions, Recommendations and Other Decisions (IUCN, 2016).

    37.
    Permitted Clearing of Native Vegetation: Biodiversity Assessment Guidelines (The Victorian Government, Department of Environment and Primary Industries, 2013).

    38.
    NSW Biodiversity Offsets Policy for Major Projects (State of NSW, Office of Environment and Heritage, 2014).

    39.
    Our Evolving Approach to Biodiversity: The Next Chapter in Biodiversity Management (Rio Tinto, 2017); http://www.riotinto.com/ourcommitment/spotlight-18130_21621.aspx

    40.
    Konisky, D. M. & Woods, N. D. Environmental federalism and the Trump presidency: a preliminary assessment. Publius 48, 345–371 (2018).
    Article  Google Scholar 

    41.
    Working for Biodiversity Net Gain: An Overview of the Business and Biodiversity Offsets Programme (BBOP) 2004–2018 (BBOP, Forest Trends, 2018).

    42.
    Summary for Policymakers of the Global Assessment Report on Biodiversity and Ecosystem Services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES Secretariat, 2019).

    43.
    Leipold, S. & Winkel, G. Discursive agency: (re-)conceptualizing actors and practices in the analysis of discursive policymaking. Policy Stud. J. 45, 510–534 (2017).
    Article  Google Scholar 

    44.
    Foucault, M. The History of Sexuality: The Will to Knowledge Vol. I (Penguin Group, 2008).

    45.
    Walker, S., Brower, A. L., Stephens, R. T. T. & Lee, W. G. Why bartering biodiversity fails. Conserv. Lett. 2, 149–157 (2009).
    Article  Google Scholar 

    46.
    Bäckstrand, K. & Lövbrand, E. The road to Paris: contending climate governance discourses in the post-Copenhagen era. J. Environ. Policy Plan. 21, 519–532 (2016).
    Article  Google Scholar 

    47.
    Griggs, S. & Howarth, D. Discourse, policy and the environment: hegemony, statements and the analysis of UK airport expansion. J. Environ. Policy Plan. 21, 464–478 (2019).
    Article  Google Scholar 

    48.
    Reflections on the Zero Draft Post-2020 Global Biodiversity Framework (BirdLife International, 2020).

    49.
    IUCN Position: Zero Draft of the Post-2020 Global Biodiversity Framework (IUCN, 2020); https://go.nature.com/3kEA4rP

    50.
    Dingler, J. The discursive nature of nature: towards a post-modern concept of nature. J. Environ. Policy Plan. 7, 209–225 (2005).
    Article  Google Scholar 

    51.
    Sharp, L. & Richardson, T. Reflections on Foucauldian discourse analysis in planning and environmental policy research. J. Environ. Policy Plan. 3, 193–209 (2001).
    Article  Google Scholar 

    52.
    Moon, K. & Blackman, D. A guide to understanding social science research for natural scientists. Conserv. Biol. 28, 1167–1177 (2014).
    Article  Google Scholar 

    53.
    Fairclough, N. Critical discourse analysis. Marges Linguist. 9, 76–94 (2005).
    Google Scholar 

    54.
    NVivo qualitative data analysis software (QSR International, 2019); https://www.qsrinternational.com/nvivo/home

    55.
    Calvet, C., Ollivier, G. & Napoleone, C. Tracking the origins and development of biodiversity offsetting in academic research and its implications for conservation: a review. Biol. Conserv. 192, 492–503 (2015).
    Article  Google Scholar 

    56.
    Darier, É. (ed.) Discourses of the Environment (Blackwell, 1999).

    57.
    Bäckstrand, K. & Lövbrand, E. in The Social Construction of Climate Change: Power, Knowledge, Norms, Discourses (ed. Pettenger, M. E.) 123–147 (Taylor & Francis Group, 2007).

    58.
    Mol, A. P. J., Spaargaren, G. & Sonnenfeld, D. A. in The Ecological Modernisation Reader. Environmental Reform in Theory and Practice (eds Mol, A. P. J. et al.) 3–14 (Routledge, 2009).

    59.
    Nilsen, H. R. The joint discourse ‘reflexive sustainable development’—from weak towards strong sustainable development. Ecol. Econ. 69, 495–501 (2010).
    Article  Google Scholar 

    60.
    Jacobs, M. in The Handbook of Global Climate and Environment Policy (ed. Falkner, R.) 197–214 (John Wiley & Sons, 2013); https://doi.org/10.1002/9781118326213.ch12

    61.
    Ferguson, P. The green economy agenda: business as usual or transformational discourse? Environ. Polit. 24, 17–37 (2015).
    Article  Google Scholar 

    62.
    Coffey, B. Unpacking the politics of natural capital and economic metaphors in environmental policy discourse. Environ. Polit. 25, 203–222 (2016).
    Article  Google Scholar 

    63.
    Bakker, K. The limits of ‘neoliberal natures’: debating green neoliberalism. Prog. Hum. Geogr. 34, 715–735 (2010).
    Article  Google Scholar  More

  • in

    Cryopreservation of testicular tissue from Murray River Rainbowfish, Melanotaenia fluviatilis

    Animal husbandry and sample collection
    All animal handling and experimental procedures were approved by the Animal Ethics Committee B at Monash Medical Centre (MMCB/2017/39) and conducted in accordance with the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes. Melanotaenia fluviatilis (Aquarium Industries, Victoria, Australia) were held at 25 °C ± 1 °C on a 12:12 light–dark cycle. At the time of experimentation, fish 5.76 cm ± 1.00 cm in length and weighing 3.25 g ± 1.38 g, were humanely killed by anesthetic overdose using aquatic anaesthetic AQUI-S (Primo Aquaculture, Queensland, Australia) and death was confirmed by destruction of the brain. The gonads were removed and placed into handling medium composed of Eagles minimum essential media (EMEM, SigmaAldrich) supplemented with 5% FBS (ThermoFisher Scientific, Victoria Australia), and 25 mM HEPES (ThermoFisher Scientific; pH 7.8) and kept on ice.
    Histology and immunohistochemistry
    Whole testes were fixed in 10% neutral buffered formalin (Merck, Victoria, Australia) for 48 h and processed by the Monash Histology Platform which included standard hematoxylin and eosin staining. Unstained sections were stained for Vasa using a zebrafish-specific anti-Vasa antibody (Sapphire Bioscience Pty. Ltd, New South Wales, Australia) and counter-stained with Hoechst (ThermoFisher Scientific). De-paraffinised sections were rehydrated through changes of xylene and a standard series of decreasing ethanol dilutions before antigen retrieval in 10 mM citrate buffer (pH 6), microwaved to boiling point for 10 min. Sections were rested in citrate buffer for 30 min prior to blocking with CAS Block (Invitrogen) for one hour followed by incubation with anti-Vasa antibody (1:200) in 5% BSA in PBS at 4 °C overnight. Sections were washed in PBS and incubated with secondary antibody, Alexa Fluor 488-conjugated goat anti-rabbit IgG (1:500; Invitrogen), and Hoechst nuclear counterstain (1:1000) in 5% BSA and PBS for one hour at room temperature.
    Images were captured using the EVOS FL Auto 2 Imaging system (ThermoFisher Scientific) and an Olympus BX43 Upright Microscope with an X-Cite Series 120 Q laser (Lumen Dynamics). Approximate cell sizes were measured using cellSens Standard imaging software (Software version: 1.16, build 15,404, Olympus) and images were analysed in FIJI23 (Software version: 2.0.0-rc-69/1.52p, Image J).
    Validation of size-based cell sorting by flow cytometry
    Using cell measurements taken from histological analysis as a guide, a size-based cell sorting method was developed to isolate our target spermatogonial cells. A set of five size-specific beads (16.5 μm, 10.2 μm, 7.56 μm, 5.11 μm, 3.3 μm, Spherotech, Lake Forest, IL, USA) were analysed on a FACS Aria Fusion flow cytometer (BD Biosciences, New South Wales, Australia). These sizes cover the range of cell sizes seen in the testis, with sperm heads being approximately 2–3 μm and spermatogonia being over 10 μm in M.fluviatilis. Due to differences in the light scattering properties of plastic beads in comparison to live cells, these bead sizes can only be interpreted as a guide of scale and not as an exact size indication for cells in suspension. Using the scatter profile produced by these beads, two gates were set: the “A” gate surrounded events in the high forward scatter region on the scatter plot, approximately 9 μm and larger to capture larger cells such as spermatogonia; the “B” gate surrounded events in a low forward scatter region, between 2—5 μm, to capture smaller germ cells such as spermatids and spermatocytes. An unstained cell suspension was then sorted through these gates and sorted cells were pelleted by centrifugation (500 g for 15mins). Images were taken of live cells in suspension using the EVOS FL Auto 2 Imaging system (ThermoFisher Scientific) and cell sizes were measured in FIJI. Samples were then fixed in 2% PFA (Thermo Fisher Scientific) for 10 min and suspended in PBS.
    Aliquots of each sample (A gate, B gate and an unsorted control) were smeared onto Superfrost Plus slides (ThermoFisher Scientific), baked overnight at 37 °C and stained with anti-Vasa antibody to determine the number of Vasa-positive cells in each sample. Briefly, the slides were washed with MilliQ water to remove any salt that was present and irrigated with wash buffer (0.1% BSA in PBS) before blocking with 10% goat serum, 0.1% Triton X in PBS for 45 min. Sections were stained with anti-Vasa antibody (1:200) in PBS containing 5% BSA for 1 h at room temperature, washed with wash buffer, incubated with Alexa Fluor 488-conjugated goat anti-rabbit IgG (1:500), and counterstained with Hoechst (1:1000). Sections were imaged on the EVOS FL Auto 2 and analysed using FIJI.
    Cryopreservation protocol
    This cryopreservation method was adapted from research described by Lee et al.14,15. Whole gonads weighing 0.0124 g ± 0.0095 g were transferred into 1.2-ml CryoTubes with 500 μl of cryomedia containing a permeating cryoprotectant, dimethyl sulfoxide (DMSO), ethylene glycol (EG), methanol or glycerol (all purchased from Merck), at concentrations ranging between 1.0 M and 2.0 M, with 0.1 M trehalose (Merck), and 1.5% BSA (Bovogen Biologicals Pty. Ltd, Victoria, Australia) in a mixed salt solution (~ 296 mOsm, pH 7.8) previously described by Lee et al.14. Control samples contained all components except the permeating cryoprotectant. Samples were equilibrated on ice for one hour and then cooled at a rate of -1 °C/minute in a CoolCell (Merck) in a -80 °C freezer for at least 3 h before being plunged into liquid nitrogen. Samples were held in liquid nitrogen for at least 24 h before thawing.
    Thawing and cell suspension preparation
    Samples were thawed in a 30 °C water bath for 1 min. The gonad was removed and gently blotted on a Kim-wipe to remove excess cryoprotectant residue and then rehydrated in three changes of handling medium (as described under “Animal husbandry and sample collection”) for 20 min per change (60 min total). After rehydration, the testis was placed in a tissue grinder with 500 μl of PBS and crushed. The tissue grinder was washed with another 500 μl of PBS resulting in a final volume of 1 ml. The cell suspension was passed through a 40 μm nylon filter to remove any large particulates prior to flow cytometry.
    Viability assessment by flow cytometry
    Cell suspensions were stained with the LIVE/DEAD Sperm Viability Kit (ThermoFisher Scientific) which included a membrane-permeating SYBR14 nucleic acid dye for detecting live cells and membrane-impermeable Propidium Iodide (PI) nucleic acid dye to detect membrane-compromised, presumably dead cells. SYBR14 was added and incubated for 5 min in the dark, followed by PI for a further 5-min incubation.
    Prior to the assessment of experimental samples, the sized beads (Spherotech) were analysed on the FACS Aria Fusion flow cytometer. Using these beads as a guide, a gate was set for the approximate size of the spermatogonial cells based on our own histological analysis of this species and previous publications on fish in general24. An unstained control and two single stain controls (PI only or SYBR14 only) were included with the experimental samples in the analysis. The sample used for the PI-only control was flash frozen in liquid nitrogen three times to ensure a high percentage of dead cell to provide an adequate count for PI staining. Flow cytometry output was analysed in FlowJoTM25. Events captured by the gate were analysed for SYBR14 and PI spectra and divided into quartiles based on the absorbance of single stain controls (Fig. 1).
    Figure 1

    Flow cytometry scatter plots and gating method. (a) Analysis of size-specific beads shows five distinct clusters. (b) A gate is set to capture events from the 9 μm measurement and above. (c) Events detected in this region are replotted to determine SYB14 and PI absorbance. Events in the Q3 region are SYB14 positive and PI negative and therefore viable. In samples treated with a negative control (d), the majority of events falls in the Q1 region, with only propidium iodide detected (e).

    Full size image

    Statistical analysis
    Statistical analysis was performed using GraphPad Prism version 8.1.2 for MacOS, GraphPad Software, La Jolla California USA, www.graphpad.com. Data is presented as mean ± standard deviation, with a p-value less than 0.05 considered statistically significant.
    For cell gating data, the proportion of cell sizes in live cell suspensions in each treatment group was analysed using a chi-square. The percentage of Vasa-positive cells in the unsorted sample and the “A” gate was analysed using an un-paired t-test; data for the “B” gate was excluded as no Vasa-positive cells were detected.
    For percentage viability data assumptions for normality and variance were met using the Shapiro–Wilk test and the Brown-Forsythe test, respectively. Following this, treatment groups were compared by one-way ANOVA and Tukey’s post hoc test. More