Evolution of antibiotic resistance impacts optimal temperature and growth rate in Escherichia coli and Staphylococcus epidermidis.

AIMS: Bacterial response to temperature changes can influence their pathogenicity to plants and humans. Changes in temperature can affect cellular and physiological responses in bacteria that can in turn affect the evolution and prevalence of antibiotic-resistance genes. Yet, how antibiotic-resistance genes influence microbial temperature response is poorly understood. METHODS AND RESULTS: We examined growth rates and physiological responses to temperature in two species-E. coli and Staph. epidermidis-after evolved resistance to 13 antibiotics. We found that evolved resistance results in species-, strain- and antibiotic-specific shifts in optimal temperature. When E. coli evolves resistance to nucleic acid and cell wall inhibitors, their optimal growth temperature decreases, and when Staph. epidermidis and E. coli evolve resistance to protein synthesis and their optimal temperature increases. Intriguingly, when Staph. epidermidis evolves resistance to Teicoplanin, fitness also increases in drug-free environments, independent of temperature response. CONCLUSION: Our results highlight how the complexity of antibiotic resistance is amplified when considering physiological responses to temperature. SIGNIFICANCE: Bacteria continuously respond to changing temperatures-whether through increased body temperature during fever, climate change or other factors. It is crucial to understand the interactions between antibiotic resistance and temperature.

Climate shapes and shifts functional biodiversity in forests worldwide.

Much ecological research aims to explain how climate impacts biodiversity and ecosystem-level processes through functional traits that link environment with individual performance. However, the specific climatic drivers of functional diversity across space and time remain unclear due largely to limitations in the availability of paired trait and climate data. We compile and analyze a global forest dataset using a method based on abundance-weighted trait moments to assess how climate influences the shapes of whole-community trait distributions. Our approach combines abundance-weighted metrics with diverse climate factors to produce a comprehensive catalog of trait-climate relationships that differ dramatically-27% of significant results change in sign and 71% disagree on sign, significance, or both-from traditional species-weighted methods. We find that (i) functional diversity generally declines with increasing latitude and elevation, (ii) temperature variability and vapor pressure are the strongest drivers of geographic shifts in functional composition and ecological strategies, and (iii) functional composition may currently be shifting over time due to rapid climate warming. Our analysis demonstrates that climate strongly governs functional diversity and provides essential information needed to predict how biodiversity and ecosystem function will respond to climate change.

Using a newly introduced framework to measure ecological stressor interactions.

Understanding how stressors combine to affect population abundances and trajectories is a fundamental ecological problem with increasingly important implications worldwide. Generalisations about interactions among stressors are challenging due to different categorisation methods and how stressors vary across species and systems. Here, we propose using a newly introduced framework to analyse data from the last 25 years on ecological stressor interactions, for example combined effects of temperature, salinity and nutrients on population survival and growth. We contrast our results with the most commonly used existing method - analysis of variance (ANOVA) - and show that ANOVA assumptions are often violated and have inherent limitations for detecting interactions. Moreover, we argue that rescaling - examining relative rather than absolute responses - is critical for ensuring that any interaction measure is independent of the strength of single-stressor effects. In contrast, non-rescaled measures - like ANOVA - find fewer interactions when single-stressor effects are weak. After re-examining 840 two-stressor combinations, we conclude that antagonism and additivity are the most frequent interaction types, in strong contrast to previous reports that synergy dominates yet supportive of more recent studies that find more antagonism. Consequently, measuring and re-assessing the frequency of stressor interaction types is imperative for a better understanding of how stressors affect populations.

Self-Similar Processes Follow a Power Law in Discrete Logarithmic Space.

Cities, wealth, and earthquakes follow continuous power-law probability distributions such as the Pareto distribution, which are canonically associated with scale-free behavior and self-similarity. However, many self-similar processes manifest as discrete steps that do not produce a continuous scale-free distribution. We construct a discrete power-law distribution that arises naturally from a simple model of hierarchical self-similar processes such as turbulence and vasculature, and we derive the maximum-likelihood estimate (MLE) for its exponent. Our distribution is self-similar, in contrast to previously studied discrete power laws such as the Zipf distribution. We show that the widely used MLE derived from the Pareto distribution leads to inaccurate estimates in systems that lack continuous scale invariance such as branching networks and data subject to logarithmic binning. We apply our MLE to data from bronchial tubes, blood vessels, and earthquakes to produce new estimates of scaling exponents and resolve contradictions among previous studies.

Social tipping points in animal societies.

Animal social groups are complex systems that are likely to exhibit tipping points-which are defined as drastic shifts in the dynamics of systems that arise from small changes in environmental conditions-yet this concept has not been carefully applied to these systems. Here, we summarize the concepts behind tipping points and describe instances in which they are likely to occur in animal societies. We also offer ways in which the study of social tipping points can open up new lines of inquiry in behavioural ecology and generate novel questions, methods, and approaches in animal behaviour and other fields, including community and ecosystem ecology. While some behaviours of living systems are hard to predict, we argue that probing tipping points across animal societies and across tiers of biological organization-populations, communities, ecosystems-may help to reveal principles that transcend traditional disciplinary boundaries.

A general model for metabolic scaling in self-similar asymmetric networks.

How a particular attribute of an organism changes or scales with its body size is known as an allometry. Biological allometries, such as metabolic scaling, have been hypothesized to result from selection to maximize how vascular networks fill space yet minimize internal transport distances and resistances. The West, Brown, Enquist (WBE) model argues that these two principles (space-filling and energy minimization) are (i) general principles underlying the evolution of the diversity of biological networks across plants and animals and (ii) can be used to predict how the resulting geometry of biological networks then governs their allometric scaling. Perhaps the most central biological allometry is how metabolic rate scales with body size. A core assumption of the WBE model is that networks are symmetric with respect to their geometric properties. That is, any two given branches within the same generation in the network are assumed to have identical lengths and radii. However, biological networks are rarely if ever symmetric. An open question is: Does incorporating asymmetric branching change or influence the predictions of the WBE model? We derive a general network model that relaxes the symmetric assumption and define two classes of asymmetrically bifurcating networks. We show that asymmetric branching can be incorporated into the WBE model. This asymmetric version of the WBE model results in several theoretical predictions for the structure, physiology, and metabolism of organisms, specifically in the case for the cardiovascular system. We show how network asymmetry can now be incorporated in the many allometric scaling relationships via total network volume. Most importantly, we show that the 3/4 metabolic scaling exponent from Kleiber's Law can still be attained within many asymmetric networks.

Optimal occlusion uniformly partitions red blood cells fluxes within a microvascular network.

In animals, gas exchange between blood and tissues occurs in narrow vessels, whose diameter is comparable to that of a red blood cell. Red blood cells must deform to squeeze through these narrow vessels, transiently blocking or occluding the vessels they pass through. Although the dynamics of vessel occlusion have been studied extensively, it remains an open question why microvessels need to be so narrow. We study occlusive dynamics within a model microvascular network: the embryonic zebrafish trunk. We show that pressure feedbacks created when red blood cells enter the finest vessels of the trunk act together to uniformly partition red blood cells through the microvasculature. Using mathematical models as well as direct observation, we show that these occlusive feedbacks are tuned throughout the trunk network to prevent the vessels closest to the heart from short-circuiting the network. Thus occlusion is linked with another open question of microvascular function: how are red blood cells delivered at the same rate to each micro-vessel? Our analysis shows that tuning of occlusive feedbacks increase the total dissipation within the network by a factor of 11, showing that uniformity of flows rather than minimization of transport costs may be prioritized by the microvascular network.

Dimensionality of consumer search space drives trophic interaction strengths.

Trophic interactions govern biomass fluxes in ecosystems, and stability in food webs. Knowledge of how trophic interaction strengths are affected by differences among habitats is crucial for understanding variation in ecological systems. Here we show how substantial variation in consumption-rate data, and hence trophic interaction strengths, arises because consumers tend to encounter resources more frequently in three dimensions (3D) (for example, arboreal and pelagic zones) than two dimensions (2D) (for example, terrestrial and benthic zones). By combining new theory with extensive data (376 species, with body masses ranging from 5.24 × 10(-14) kg to 800 kg), we find that consumption rates scale sublinearly with consumer body mass (exponent of approximately 0.85) for 2D interactions, but superlinearly (exponent of approximately 1.06) for 3D interactions. These results contradict the currently widespread assumption of a single exponent (of approximately 0.75) in consumer-resource and food-web research. Further analysis of 2,929 consumer-resource interactions shows that dimensionality of consumer search space is probably a major driver of species coexistence, and the stability and abundance of populations.

When more is less: Emergent suppressive interactions in three-drug combinations.

BACKGROUND: In drug-drug interactions, there are surprising cases in which the growth inhibition of bacteria by a single antibiotic decreases when a second antibiotic is added. These interactions are termed suppressive and have been argued to have the potential to limit the evolution of resistance. Nevertheless, little attention has been given to suppressive interactions because clinical studies typically search for increases in killing efficiency and because suppressive interactions are believed to be rare based on pairwise studies. RESULTS: Here, we quantify the effects of single-, double-, and triple-drug combinations from a set of 14 antibiotics and 3 bacteria strains, totaling 364 unique three-drug combinations per bacteria strain. We find that increasing the number of drugs can increase the prevalence of suppressive interactions: 17% of three-drug combinations are suppressive compared to 5% of two-drug combinations in this study. Most cases of suppression we find (97%) are "hidden" cases for which the triple-drug bacterial growth is less than the single-drug treatments but exceeds that of a pairwise combination. CONCLUSIONS: We find a surprising number of suppressive interactions in higher-order drug combinations. Without examining lower-order (pairwise) bacterial growth, emergent suppressive effects would be missed, potentially affecting our understanding of evolution of resistance and treatment strategies for resistant pathogens. These findings suggest that careful examination of the full factorial of drug combinations is needed to uncover suppressive interactions in higher-order combinations.

Do Vascular Networks Branch Optimally or Randomly across Spatial Scales?

Modern models that derive allometric relationships between metabolic rate and body mass are based on the architectural design of the cardiovascular system and presume sibling vessels are symmetric in terms of radius, length, flow rate, and pressure. Here, we study the cardiovascular structure of the human head and torso and of a mouse lung based on three-dimensional images processed via our software Angicart. In contrast to modern allometric theories, we find systematic patterns of asymmetry in vascular branching, potentially explaining previously documented mismatches between predictions (power-law or concave curvature) and observed empirical data (convex curvature) for the allometric scaling of metabolic rate. To examine why these systematic asymmetries in vascular branching might arise, we construct a mathematical framework to derive predictions based on local, junction-level optimality principles that have been proposed to be favored in the course of natural selection and development. The two most commonly used principles are material-cost optimizations (construction materials or blood volume) and optimization of efficient flow via minimization of power loss. We show that material-cost optimization solutions match with distributions for asymmetric branching across the whole network but do not match well for individual junctions. Consequently, we also explore random branching that is constrained at scales that range from local (junction-level) to global (whole network). We find that material-cost optimizations are the strongest predictor of vascular branching in the human head and torso, whereas locally or intermediately constrained random branching is comparable to material-cost optimizations for the mouse lung. These differences could be attributable to developmentally-programmed local branching for larger vessels and constrained random branching for smaller vessels.

Real versus Artificial Variation in the Thermal Sensitivity of Biological Traits.

Whether the thermal sensitivity of an organism's traits follows the simple Boltzmann-Arrhenius model remains a contentious issue that centers around consideration of its operational temperature range and whether the sensitivity corresponds to one or a few underlying rate-limiting enzymes. Resolving this issue is crucial, because mechanistic models for temperature dependence of traits are required to predict the biological effects of climate change. Here, by combining theory with data on 1,085 thermal responses from a wide range of traits and organisms, we show that substantial variation in thermal sensitivity (activation energy) estimates can arise simply because of variation in the range of measured temperatures. Furthermore, when thermal responses deviate systematically from the Boltzmann-Arrhenius model, variation in measured temperature ranges across studies can bias estimated activation energy distributions toward higher mean, median, variance, and skewness. Remarkably, this bias alone can yield activation energies that encompass the range expected from biochemical reactions (from ~0.2 to 1.2 eV), making it difficult to establish whether a single activation energy appropriately captures thermal sensitivity. We provide guidelines and a simple equation for partially correcting for such artifacts. Our results have important implications for understanding the mechanistic basis of thermal responses of biological traits and for accurately modeling effects of variation in thermal sensitivity on responses of individuals, populations, and ecological communities to changing climatic temperatures.

Testing Foundations of Biological Scaling Theory Using Automated Measurements of Vascular Networks.

Scientists have long sought to understand how vascular networks supply blood and oxygen to cells throughout the body. Recent work focuses on principles that constrain how vessel size changes through branching generations from the aorta to capillaries and uses scaling exponents to quantify these changes. Prominent scaling theories predict that combinations of these exponents explain how metabolic, growth, and other biological rates vary with body size. Nevertheless, direct measurements of individual vessel segments have been limited because existing techniques for measuring vasculature are invasive, time consuming, and technically difficult. We developed software that extracts the length, radius, and connectivity of in vivo vessels from contrast-enhanced 3D Magnetic Resonance Angiography. Using data from 20 human subjects, we calculated scaling exponents by four methods-two derived from local properties of branching junctions and two from whole-network properties. Although these methods are often used interchangeably in the literature, we do not find general agreement between these methods, particularly for vessel lengths. Measurements for length of vessels also diverge from theoretical values, but those for radius show stronger agreement. Our results demonstrate that vascular network models cannot ignore certain complexities of real vascular systems and indicate the need to discover new principles regarding vessel lengths.

The body size dependence of trophic cascades.

Trophic cascades are indirect positive effects of predators on resources via control of intermediate consumers. Larger-bodied predators appear to induce stronger trophic cascades (a greater rebound of resource density toward carrying capacity), but how this happens is unknown because we lack a clear depiction of how the strength of trophic cascades is determined. Using consumer resource models, we first show that the strength of a trophic cascade has an upper limit set by the interaction strength between the basal trophic group and its consumer and that this limit is approached as the interaction strength between the consumer and its predator increases. We then express the strength of a trophic cascade explicitly in terms of predator body size and use two independent parameter sets to calculate how the strength of a trophic cascade depends on predator size. Both parameter sets predict a positive effect of predator size on the strength of a trophic cascade, driven mostly by the body size dependence of the interaction strength between the first two trophic levels. Our results support previous empirical findings and suggest that the loss of larger predators will have greater consequences on trophic control and biomass structure in food webs than the loss of smaller predators.

Deviation from symmetrically self-similar branching in trees predicts altered hydraulics, mechanics, light interception and metabolic scaling.

The West, Brown, Enquist (WBE) model derives symmetrically self-similar branching to predict metabolic scaling from hydraulic conductance, K, (a metabolism proxy) and tree mass (or volume, V). The original prediction was Kα V(0.75). We ask whether trees differ from WBE symmetry and if it matters for plant function and scaling. We measure tree branching and model how architecture influences K, V, mechanical stability, light interception and metabolic scaling. We quantified branching architecture by measuring the path fraction, Pf : mean/maximum trunk-to-twig pathlength. WBE symmetry produces the maximum, Pf = 1.0. We explored tree morphospace using a probability-based numerical model constrained only by biomechanical principles. Real tree Pf ranged from 0.930 (nearly symmetric) to 0.357 (very asymmetric). At each modeled tree size, a reduction in Pf led to: increased K; decreased V; increased mechanical stability; and decreased light absorption. When Pf was ontogenetically constant, strong asymmetry only slightly steepened metabolic scaling. The Pf ontogeny of real trees, however, was 'U' shaped, resulting in size-dependent metabolic scaling that exceeded 0.75 in small trees before falling below 0.65. Architectural diversity appears to matter considerably for whole-tree hydraulics, mechanics, photosynthesis and potentially metabolic scaling. Optimal architectures likely exist that maximize carbon gain per structural investment.

Temperature dependence of trophic interactions are driven by asymmetry of species responses and foraging strategy.

Environmental temperature has systematic effects on rates of species interactions, primarily through its influence on organismal physiology. We present a mechanistic model for the thermal response of consumer-resource interactions. We focus on how temperature affects species interactions via key traits - body velocity, detection distance, search rate and handling time - that underlie per capita consumption rate. The model is general because it applies to all foraging strategies: active-capture (both consumer and resource body velocity are important), sit-and-wait (resource velocity dominates) and grazing (consumer velocity dominates). The model predicts that temperature influences consumer-resource interactions primarily through its effects on body velocity (either of the consumer, resource or both), which determines how often consumers and resources encounter each other, and that asymmetries in the thermal responses of interacting species can introduce qualitative, not just quantitative, changes in consumer-resource dynamics. We illustrate this by showing how asymmetries in thermal responses determine equilibrium population densities in interacting consumer-resource pairs. We test for the existence of asymmetries in consumer-resource thermal responses by analysing an extensive database on thermal response curves of ecological traits for 309 species spanning 15 orders of magnitude in body size from terrestrial, marine and freshwater habitats. We find that asymmetries in consumer-resource thermal responses are likely to be a common occurrence. Overall, our study reveals the importance of asymmetric thermal responses in consumer-resource dynamics. In particular, we identify three general types of asymmetries: (i) different levels of performance of the response, (ii) different rates of response (e.g. activation energies) and (iii) different peak or optimal temperatures. Such asymmetries should occur more frequently as the climate changes and species' geographical distributions and phenologies are altered, such that previously noninteracting species come into contact. 6. By using characteristics of trophic interactions that are often well known, such as body size, foraging strategy, thermy and environmental temperature, our framework should allow more accurate predictions about the thermal dependence of consumer-resource interactions. Ultimately, integration of our theory into models of food web and ecosystem dynamics should be useful in understanding how natural systems will respond to current and future temperature change.

Systematic variation in the temperature dependence of physiological and ecological traits.

To understand the effects of temperature on biological systems, we compile, organize, and analyze a database of 1,072 thermal responses for microbes, plants, and animals. The unprecedented diversity of traits (n = 112), species (n = 309), body sizes (15 orders of magnitude), and habitats (all major biomes) in our database allows us to quantify novel features of the temperature response of biological traits. In particular, analysis of the rising component of within-species (intraspecific) responses reveals that 87% are fit well by the Boltzmann-Arrhenius model. The mean activation energy for these rises is 0.66 ± 0.05 eV, similar to the reported across-species (interspecific) value of 0.65 eV. However, systematic variation in the distribution of rise activation energies is evident, including previously unrecognized right skewness around a median of 0.55 eV. This skewness exists across levels of organization, taxa, trophic groups, and habitats, and it is partially explained by prey having increased trait performance at lower temperatures relative to predators, suggesting a thermal version of the life-dinner principle-stronger selection on running for your life than running for your dinner. For unimodal responses, habitat (marine, freshwater, and terrestrial) largely explains the mean temperature at which trait values are optimal but not variation around the mean. The distribution of activation energies for trait falls has a mean of 1.15 ± 0.39 eV (significantly higher than rises) and is also right-skewed. Our results highlight generalities and deviations in the thermal response of biological traits and help to provide a basis to predict better how biological systems, from cells to communities, respond to temperature change.

Multilevel relations among plankton stitched together with an eco-evolutionary needle.

Power-law relationships between population abundances, energy use, and other factors are often referred to as macroecological scaling. A recent study convincingly shows that these relationships emerge from individual physiology but only after the population distribution is shaped by trophic interactions that are subject to both ecological and evolutionary pressures.

A flexible model for thermal performance curves.

Temperature responses of many biological traits-including population growth, survival, and development-are described by thermal performance curves (TPCs) with phenomenological models like the Briere function or mechanistic models related to chemical kinetics. Existing TPC models are either simple but inflexible in shape, or flexible yet difficult to interpret in biological terms. Here we present flexTPC: a model that is parameterized exclusively in terms of biologically interpretable quantities, including the thermal minimum, optimum, and maximum, and the maximum trait value. FlexTPC can describe unimodal temperature responses of any skewness and thermal breadth, enabling direct comparisons across populations, traits, or taxa with a single model. We apply flexTPC to various microbial and entomological datasets, compare results with the Briere model, and find that flexTPC often has better predictive performance. The interpretability of flexTPC makes it ideal for modeling how thermal responses change with ecological stressors or evolve over time.

Mean daily temperatures can predict the thermal limits of malaria transmission better than rate summation.

Temperature shapes the distribution, seasonality, and magnitude of mosquito-borne disease outbreaks. Mechanistic models predicting transmission often use mosquito and pathogen thermal responses from constant temperature experiments. However, mosquitoes live in fluctuating environments. Rate summation (nonlinear averaging) is a common approach to infer performance in fluctuating environments, but its accuracy is rarely validated. We measured three mosquito traits that impact transmission (bite rate, survival, fecundity) in a malaria mosquito (Anopheles stephensi) across temperature gradients with three diurnal temperature ranges (0, 9 and 12°C). We compared thermal suitability models with temperature-trait relationships observed under constant temperatures, fluctuating temperatures, and those predicted by rate summation. We mapped results across An. stephenesi's native Asian and invasive African ranges. We found: 1) daily temperature fluctuation significantly altered trait thermal responses; 2) rate summation partially captured decreases in performance near thermal optima, but also incorrectly predicted increases near thermal limits; and 3) while thermal suitability characterized across constant temperatures did not perfectly capture suitability in fluctuating environments, it was more accurate for estimating and mapping thermal limits than predictions from rate summation. Our study provides insight into methods for predicting mosquito-borne disease risk and emphasizes the need to improve understanding of organismal performance under fluctuating conditions.

Evolution Under Thermal Stress Affects Escherichia coli's Resistance to Antibiotics.

Exposure to both antibiotics and temperature changes can induce similar physiological responses in bacteria. Thus, changes in growth temperature may affect antibiotic resistance. Previous studies have found that evolution under antibiotic stress causes shifts in the optimal growth temperature of bacteria. However, little is known about how evolution under thermal stress affects antibiotic resistance. We examined 100+ heat-evolved strains of Escherichia coli that evolved under thermal stress. We asked whether evolution under thermal stress affects optimal growth temperature, if there are any correlations between evolving in high temperatures and antibiotic resistance, and if these strains' antibiotic efficacy changes depending on the local environment's temperature. We found that: (1) surprisingly, most of the heat-evolved strains displayed a decrease in optimal growth temperature and overall growth relative to the ancestor strain, (2) there were complex patterns of changes in antibiotic resistance when comparing the heat-evolved strains to the ancestor strain, and (3) there were few significant correlations among changes in antibiotic resistance, optimal growth temperature, and overall growth.