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.

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.

Neuronal branching is increasingly asymmetric near synapses, potentially enabling plasticity while minimizing energy dissipation and conduction time.

Neurons' primary function is to encode and transmit information in the brain and body. The branching architecture of axons and dendrites must compute, respond and make decisions while obeying the rules of the substrate in which they are enmeshed. Thus, it is important to delineate and understand the principles that govern these branching patterns. Here, we present evidence that asymmetric branching is a key factor in understanding the functional properties of neurons. First, we derive novel predictions for asymmetric scaling exponents that encapsulate branching architecture associated with crucial principles such as conduction time, power minimization and material costs. We compare our predictions with extensive data extracted from images to associate specific principles with specific biophysical functions and cell types. Notably, we find that asymmetric branching models lead to predictions and empirical findings that correspond to different weightings of the importance of maximum, minimum or total path lengths from the soma to the synapses. These different path lengths quantitatively and qualitatively affect energy, time and materials. Moreover, we generally observe that higher degrees of asymmetric branching-potentially arising from extrinsic environmental cues and synaptic plasticity in response to activity-occur closer to the tips than the soma (cell body).

Neuronal Branching is Increasingly Asymmetric Near Synapses, Potentially Enabling Plasticity While Minimizing Energy Dissipation and Conduction Time.

Neurons' primary function is to encode and transmit information in the brain and body. The branching architecture of axons and dendrites must compute, respond, and make decisions while obeying the rules of the substrate in which they are enmeshed. Thus, it is important to delineate and understand the principles that govern these branching patterns. Here, we present evidence that asymmetric branching is a key factor in understanding the functional properties of neurons. First, we derive novel predictions for asymmetric scaling exponents that encapsulate branching architecture associated with crucial principles such as conduction time, power minimization, and material costs. We compare our predictions with extensive data extracted from images to associate specific principles with specific biophysical functions and cell types. Notably, we find that asymmetric branching models lead to predictions and empirical findings that correspond to different weightings of the importance of maximum, minimum, or total path lengths from the soma to the synapses. These different path lengths quantitatively and qualitatively affect energy, time, and materials. Moreover, we generally observe that higher degrees of asymmetric branching- potentially arising from extrinsic environmental cues and synaptic plasticity in response to activity- occur closer to the tips than the soma (cell body).

Stimuli-sensitive nano-drug delivery with programmable size changes to enhance accumulation of therapeutic agents in tumors.

Nano-based drug delivery systems hold significant promise for cancer therapies. Presently, the poor accumulation of drug-carrying nanoparticles in tumors has limited their success. In this study, based on a combination of the paradigms of intravascular and extravascular drug release, an efficient nanosized drug delivery system with programmable size changes is introduced. Drug-loaded smaller nanoparticles (secondary nanoparticles), which are loaded inside larger nanoparticles (primary nanoparticles), are released within the microvascular network due to temperature field resulting from focused ultrasound. This leads to the scale of the drug delivery system decreasing by 7.5 to 150 times. Subsequently, smaller nanoparticles enter the tissue at high transvascular rates and achieve higher accumulation, leading to higher penetration depths. In response to the acidic pH of tumor microenvironment (according to the distribution of oxygen), they begin to release the drug doxorubicin at very slow rates (i.e., sustained release). To predict the performance and distribution of therapeutic agents, a semi-realistic microvascular network is first generated based on a sprouting angiogenesis model and the transport of therapeutic agents is then investigated based on a developed multi-compartment model. The results show that reducing the size of the primary and secondary nanoparticles can lead to higher cell death rate. In addition, tumor growth can be inhibited for a longer time by enhancing the bioavailability of the drug in the extracellular space. The proposed drug delivery system can be very promising in clinical applications. Furthermore, the proposed mathematical model is applicable to broader applications to predict the performance of drug delivery systems.

Meta-analysis of three-stressor combinations on population-level fitness reveal substantial higher-order interactions.

Although natural populations are typically subjected to multiple stressors, most past research has focused on single-stressor and two-stressor interactions, with little attention paid to higher-order interactions among three or more stressors. However, higher-order interactions increasingly appear to be widespread. Consequently, we used a recently introduced and improved framework to re-analyze higher-order ecological interactions. We conducted a literature review of the last 100 years (1920-2020) and reanalyzed 142 ecological three-stressor interactions on species' populations from 38 published papers; the vast majority of these studies were from the past 10 years. We found that 95.8 % (n = 136) of the three-stressor combinations had either not been categorized before or resulted in different interactions than previously reported. We also found substantial levels of emergent properties-interactions that are not due to strong pairwise interactions within the combination but rather uniquely due to all three stressors being combined. Calculating net interactions-the overall accounting for all possible interactions within a combination including the emergent and all pairwise interactions-we found that the most prevalent interaction type is antagonism, corresponding to a smaller than expected effect based on single stressor effects. In contrast, for emergent interactions, the most prevalent interaction type is synergistic, resulting in a larger than expected effect based on single stressor effects. Additionally, we found that hidden suppressive interactions-where a pairwise interaction is suppressed by a third stressor-are found in the majority of combinations (74 %). Collectively, understanding multiple stressor interactions through applying an appropriate framework is crucial for answering fundamental questions in ecology and has implications for conservation biology and population management. Crucially, identifying emergent properties can reveal hidden suppressive interactions that could be particularly important for the ecological management of at-risk populations.

How axon and dendrite branching are guided by time, energy, and spatial constraints.

Neurons are connected by complex branching processes-axons and dendrites-that process information for organisms to respond to their environment. Classifying neurons according to differences in structure or function is a fundamental part of neuroscience. Here, by constructing biophysical theory and testing against empirical measures of branching structure, we develop a general model that establishes a correspondence between neuron structure and function as mediated by principles such as time or power minimization for information processing as well as spatial constraints for forming connections. We test our predictions for radius scale factors against those extracted from neuronal images, measured for species that range from insects to whales, including data from light and electron microscopy studies. Notably, our findings reveal that the branching of axons and peripheral nervous system neurons is mainly determined by time minimization, while dendritic branching is determined by power minimization. Our model also predicts a quarter-power scaling relationship between conduction time delay and body size.

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.

Antibiotics Shift the Temperature Response Curve of Escherichia coli Growth.

Temperature variation-through time and across climatic gradients-affects individuals, populations, and communities. Yet how the thermal response of biological systems is altered by environmental stressors is poorly understood. Here, we quantify two key features-optimal temperature and temperature breadth-to investigate how temperature responses vary in the presence of antibiotics. We use high-throughput screening to measure growth of Escherichia coli under single and pairwise combinations of 12 antibiotics across seven temperatures that range from 22°C to 46°C. We find that antibiotic stress often results in considerable changes in the optimal temperature for growth and a narrower temperature breadth. The direction of the optimal temperature shifts can be explained by the similarities between antibiotic-induced and temperature-induced damage to the physiology of the bacterium. We also find that the effects of pairs of stressors in the temperature response can often be explained by just one antibiotic out of the pair. Our study has implications for a general understanding of how ecological systems adapt and evolve to environmental changes. IMPORTANCE The growth of living organisms varies with temperature. This dependence is described by a temperature response curve that is described by an optimal temperature where growth is maximized and a temperature range (termed breadth) across which the organism can grow. Because an organism's temperature response evolves or acclimates to its environment, it is often assumed to change over only evolutionary or developmental timescales. Counter to this, we show here that antibiotics can quickly (over hours) change the optimal growth temperature and temperature breadth for the bacterium Escherichia coli. Moreover, our results suggest a shared-damage hypothesis: when an antibiotic damages similar cellular components as hot (or cold) temperatures do, this shared damage will combine and compound to more greatly reduce growth when that antibiotic is administered at hot (or cold) temperatures. This hypothesis could potentially also explain how temperature responses are modified by stressors other than antibiotics.

The allometry of locomotion.

Organismal locomotion mediates ecological interactions and shapes community dynamics. Locomotion is constrained by intrinsic and environmental factors and integrating these factors should clarify how locomotion affects ecology across scales. We extended general theory based on metabolic scaling and biomechanics to predict the scaling of five locomotor performance traits: routine speed, maximum speed, maximum acceleration, minimum powered turn radius, and angular speed. To test these predictions, we used phylogenetically informed analyses of a new database with 884 species and found support for our quantitative predictions. Larger organisms were faster but less maneuverable than smaller organisms. Routine and maximum speeds scaled with body mass to 0.20 and 0.17 powers, respectively, and plateaued at higher body masses, especially for maximum speed. Acceleration was unaffected by body mass. Minimum turn radius scaled to a 0.19 power, and the 95% CI included our theoretical prediction, as we predicted. Maximum angular speed scaled higher than predicted but in the same direction. We observed universal scaling among locomotor modes for routine and maximum speeds but the intercepts varied; flying organisms were faster than those that swam or ran. Acceleration was independent of size in flying and aquatic taxa but decreased with body mass in land animals, possibly due to the risk of injury large, terrestrial organisms face at high speeds and accelerations. Terrestrial mammals inhabiting structurally simple habitats tended to be faster than those in complex habitats. Despite effects of body size, locomotor mode, and habitat complexity, universal scaling of locomotory performance reveals the general ways organisms move across Earth's complex environments.

Hidden suppressive interactions are common in higher-order drug combinations.

The rapid increase of multi-drug resistant bacteria has led to a greater emphasis on multi-drug combination treatments. However, some combinations can be suppressive-that is, bacteria grow faster in some drug combinations than when treated with a single drug. Typically, when studying interactions, the overall effect of the combination is only compared with the single-drug effects. However, doing so could miss "hidden" cases of suppression, which occur when the highest order is suppressive compared with a lower-order combination but not to a single drug. We examined an extensive dataset of 5-drug combinations and all lower-order-single, 2-, 3-, and 4-drug-combinations. We found that a majority of all combinations-54%-contain hidden suppression. Examining hidden interactions is critical to understanding the architecture of higher-order interactions and can substantially affect our understanding and predictions of the evolution of antibiotic resistance under multi-drug treatments.

Branching principles of animal and plant networks identified by combining extensive data, machine learning and modelling.

Branching in vascular networks and in overall organismic form is one of the most common and ancient features of multicellular plants, fungi and animals. By combining machine-learning techniques with new theory that relates vascular form to metabolic function, we enable novel classification of diverse branching networks-mouse lung, human head and torso, angiosperm and gymnosperm plants. We find that ratios of limb radii-which dictate essential biologic functions related to resource transport and supply-are best at distinguishing branching networks. We also show how variation in vascular and branching geometry persists despite observing a convergent relationship across organisms for how metabolic rate depends on body mass.

Unraveling why we sleep: Quantitative analysis reveals abrupt transition from neural reorganization to repair in early development.

Sleep serves disparate functions, most notably neural repair, metabolite clearance and circuit reorganization. Yet the relative importance remains hotly debated. Here, we create a novel mechanistic framework for understanding and predicting how sleep changes during ontogeny and across phylogeny. We use this theory to quantitatively distinguish between sleep used for neural reorganization versus repair. Our findings reveal an abrupt transition, between 2 and 3 years of age in humans. Specifically, our results show that differences in sleep across phylogeny and during late ontogeny (after 2 or 3 years in humans) are primarily due to sleep functioning for repair or clearance, while changes in sleep during early ontogeny (before 2 or 3 years) primarily support neural reorganization and learning. Moreover, our analysis shows that neuroplastic reorganization occurs primarily in REM sleep but not in NREM. This developmental transition suggests a complex interplay between developmental and evolutionary constraints on sleep.

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.

Predicting collapse of complex ecological systems: quantifying the stability-complexity continuum.

Dynamical shifts between the extremes of stability and collapse are hallmarks of ecological systems. These shifts are limited by and change with biodiversity, complexity, and the topology and hierarchy of interactions. Most ecological research has focused on identifying conditions for a system to shift from stability to any degree of instability-species abundances do not return to exact same values after perturbation. Real ecosystems likely have a continuum of shifting between stability and collapse that depends on the specifics of how the interactions are structured, as well as the type and degree of disturbance due to environmental change. Here we map boundaries for the extremes of strict stability and collapse. In between these boundaries, we find an intermediate regime that consists of single-species extinctions, which we call the extinction continuum. We also develop a metric that locates the position of the system within the extinction continuum-thus quantifying proximity to stability or collapse-in terms of ecologically measurable quantities such as growth rates and interaction strengths. Furthermore, we provide analytical and numerical techniques for estimating our new metric. We show that our metric does an excellent job of capturing the system's behaviour in comparison with other existing methods-such as May's stability criteria or critical slowdown. Our metric should thus enable deeper insights about how to classify real systems in terms of their overall dynamics and their limits of stability and collapse.

Stability of ecosystems enhanced by species-interaction constraints.

Ecosystem stability is a central question both in theoretical and applied biology. Dynamical systems theory can be used to analyze how growth rates, carrying capacities, and patterns of species interactions affect the stability of an ecosystem. The response to increasing complexity has been extensively studied and the general conclusion is that there is a limit. While there is a complexity limit to stability at which global destabilisation occurs, the collapse rarely happens suddenly if a system is fully viable (no species is extinct). In fact, when complexity is successively increased, we find that the generic response is to go through multiple single-species extinctions before a global collapse. In this paper we demonstrate this finding via both numerical simulations and elaborations of theoretical predictions. We explore more biological interaction patterns, and, perhaps most importantly, we show that constrained interaction structures-a constant row sum in the interaction matrix-prevent extinctions from occurring. This makes an ecosystem more robust in terms of allowed complexity, but it also means singles-species extinctions do not precede or signal collapse-a drastically different behavior compared to the generic and commonly assumed case. We further argue that this constrained interaction structure-limiting the total interactions for each species-is biologically plausible.

Informing trait-based ecology by assessing remotely sensed functional diversity across a broad tropical temperature gradient.

Spatially continuous data on functional diversity will improve our ability to predict global change impacts on ecosystem properties. We applied methods that combine imaging spectroscopy and foliar traits to estimate remotely sensed functional diversity in tropical forests across an Amazon-to-Andes elevation gradient (215 to 3537 m). We evaluated the scale dependency of community assembly processes and examined whether tropical forest productivity could be predicted by remotely sensed functional diversity. Functional richness of the community decreased with increasing elevation. Scale-dependent signals of trait convergence, consistent with environmental filtering, play an important role in explaining the range of trait variation within each site and along elevation. Single- and multitrait remotely sensed measures of functional diversity were important predictors of variation in rates of net and gross primary productivity. Our findings highlight the potential of remotely sensed functional diversity to inform trait-based ecology and trait diversity-ecosystem function linkages in hyperdiverse tropical forests.

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.

Stressor interaction networks suggest antibiotic resistance co-opted from stress responses to temperature.

Environmental factors like temperature, pressure, and pH partly shaped the evolution of life. As life progressed, new stressors (e.g., poisons and antibiotics) arose as part of an arms race among organisms. Here we ask if cells co-opted existing mechanisms to respond to new stressors, or whether new responses evolved de novo. We use a network-clustering approach based purely on phenotypic growth measurements and interactions among the effects of stressors on population growth. We apply this method to two types of stressors-temperature and antibiotics-to discover the extent to which their cellular responses overlap in Escherichia coli. Our clustering reveals that responses to low and high temperatures are clearly separated, and each is grouped with responses to antibiotics that have similar effects to cold or heat, respectively. As further support, we use a library of transcriptional fluorescent reporters to confirm heat-shock and cold-shock genes are induced by antibiotics. We also show strains evolved at high temperatures are more sensitive to antibiotics that mimic the effects of cold. Taken together, our results strongly suggest that temperature stress responses have been co-opted to deal with antibiotic stress.

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.