Emergent encoding of dispersal network topologies in spatial metapopulation models.

We address a generalization of the concept of metapopulation capacity for trees and networks acting as the template for ecological interactions. The original measure had been derived from an insightful phenomenological model and is based on the leading eigenvalue of a suitable landscape matrix. It yields a versatile predictor of metapopulation persistence through a threshold value of the eigenvalue determined by ecological features of the focal species. Here, we present an analytical solution to a fundamental microscopic model that incorporates key ingredients of metapopulation dynamics and explicitly distinguishes between individuals comprising the "settled population" and "explorers" seeking colonization. Our approach accounts for general network characteristics (in particular graph-driven directional dispersal which is known to significantly constrain many ecological estimates) and yields a generalized version of the original model, to which it reduces for particular cases. Through examples, including real landscapes used as the template, we compare the predictions from our approach with those of the standard model. Results suggest that in several cases of practical interest, differences are significant. We also examine, with both models, how changes in habitat fragmentation, including removal, addition, or alteration in size, affect metapopulation persistence. The current approach demonstrates a high level of flexibility, enabling the incorporation of diverse "microscopic" elements and their impact on the resulting biodiversity landscape pattern.

Delay effects on the stability of large ecosystems.

The common intuition among the ecologists of the midtwentieth century was that large ecosystems should be more stable than those with a smaller number of species. This view was challenged by Robert May, who found a stability bound for randomly assembled ecosystems; they become unstable for a sufficiently large number of species. In the present work, we show that May's bound greatly changes when the past population densities of a species affect its own current density. This is a common feature in real systems, where the effects of species' interactions may appear after a time lag rather than instantaneously. The local stability of these models with self-interaction is described by bounds, which we characterize in the parameter space. We find a critical delay curve that separates the region of stability from that of instability, and correspondingly, we identify a critical frequency curve that provides the characteristic frequencies of a system at the instability threshold. Finally, we calculate analytically the distributions of eigenvalues that generalize Wigner's as well as Girko's laws. Interestingly, we find that, for sufficiently large delays, the eigenvalues of a randomly coupled system are complex even when the interactions are symmetric.

III. Geometrical framework for thinking about globular proteins: Turns in proteins.

We have shown recently that the notion of poking pairwise interactions along a chain provides a unifying framework for understanding the formation of both secondary and the tertiary protein structure based on symmetry and geometry. α-helices and ß-sheets are found to be special geometries that have systematic poking contacts in a repetitive manner with the contacts being local along the α-helix and non-local along a pair of adjacent strands within a ß-sheet. Pairwise poking interactions also govern tertiary structure formation, but they are weaker and there are no special geometrical constraints as in secondary structure formation. Here we demonstrate that protein turns, the most prevalent non-repetitive structural element in proteins, are instances of local (as in α-helices) and isolated (non-repetitive) poking pairwise contacts for which the geometrical constraints are partially relaxed. This simple and purely geometrical definition of protein turns (also sometimes known as reverse turns, ß-turns, ß-bends, hairpin bends, 310 bends, kinks, widgets, etc.) provides a simple framework for unifying them. We present the results of a systematic analysis and identify their structural classes as well as their respective amino acid preferences.

Scaling of joint mass and metabolism fluctuations in in silico cell-laden spheroids.

Variations and fluctuations are characteristic features of biological systems and are also manifested in cell cultures. Here, we describe a computational pipeline for identifying the range of three-dimensional (3D) cell-aggregate sizes in which nonisometric scaling emerges in the presence of joint mass and metabolic rate fluctuations. The 3D cell-laden spheroids with size and single-cell metabolic rates described by probability density functions were randomly generated in silico. The distributions of the resulting metabolic rates of the spheroids were computed by modeling oxygen diffusion and reaction. Then, a method for estimating scaling exponents of correlated variables through statistically significant data collapse of joint probability distributions was developed. The method was used to identify a physiologically relevant range of spheroid sizes, where both nonisometric scaling and a minimum oxygen concentration (0.04 mol⋅m-3) is maintained. The in silico pipeline described enables the prediction of the number of experiments needed for an acceptable collapse and, thus, a consistent estimate of scaling parameters. Using the pipeline, we also show that scaling exponents may be significantly different in the presence of joint mass and metabolic-rate variations typically found in cells. Our study highlights the importance of incorporating fluctuations and variability in size and metabolic rates when estimating scaling exponents. It also suggests the need for taking into account their covariations for better understanding and interpreting experimental observations both in vitro and in vivo and brings insights for the design of more predictive and physiologically relevant in vitro models.

True scale-free networks hidden by finite size effects.

We analyze about 200 naturally occurring networks with distinct dynamical origins to formally test whether the commonly assumed hypothesis of an underlying scale-free structure is generally viable. This has recently been questioned on the basis of statistical testing of the validity of power law distributions of network degrees. Specifically, we analyze by finite size scaling analysis the datasets of real networks to check whether the purported departures from power law behavior are due to the finiteness of sample size. We find that a large number of the networks follows a finite size scaling hypothesis without any self-tuning. This is the case of biological protein interaction networks, technological computer and hyperlink networks, and informational networks in general. Marked deviations appear in other cases, especially involving infrastructure and transportation but also in social networks. We conclude that underlying scale invariance properties of many naturally occurring networks are extant features often clouded by finite size effects due to the nature of the sample data.

The Widened Pipe Model of plant hydraulic evolution.

Shaping global water and carbon cycles, plants lift water from roots to leaves through xylem conduits. The importance of xylem water conduction makes it crucial to understand how natural selection deploys conduit diameters within and across plants. Wider conduits transport more water but are likely more vulnerable to conduction-blocking gas embolisms and cost more for a plant to build, a tension necessarily shaping xylem conduit diameters along plant stems. We build on this expectation to present the Widened Pipe Model (WPM) of plant hydraulic evolution, testing it against a global dataset. The WPM predicts that xylem conduits should be narrowest at the stem tips, widening quickly before plateauing toward the stem base. This universal profile emerges from Pareto modeling of a trade-off between just two competing vectors of natural selection: one favoring rapid widening of conduits tip to base, minimizing hydraulic resistance, and another favoring slow widening of conduits, minimizing carbon cost and embolism risk. Our data spanning terrestrial plant orders, life forms, habitats, and sizes conform closely to WPM predictions. The WPM highlights carbon economy as a powerful vector of natural selection shaping plant function. It further implies that factors that cause resistance in plant conductive systems, such as conduit pit membrane resistance, should scale in exact harmony with tip-to-base conduit widening. Furthermore, the WPM implies that alterations in the environments of individual plants should lead to changes in plant height, for example, shedding terminal branches and resprouting at lower height under drier climates, thus achieving narrower and potentially more embolism-resistant conduits.

II. Geometrical framework for thinking about globular proteins: The power of poking.

Recently, we presented a framework for understanding protein structure based on the idea that simple constructs of holding hands or touching of objects can be used to rationalize the common characteristics of globular proteins. We developed a consistent approach for understanding the formation of the two key common building blocks of helices and sheets as well as the compatible assembly of secondary structures into the tertiary structure through the notion of poking pairwise interactions. Here we benchmark our predictions with a detailed analysis of structural data of over 4000 proteins from the Protein Data Bank. We also present the results of detailed computer simulations of a simplified model demonstrating a pre-sculpted free energy landscape, determined by geometry and symmetry, comprising numerous minima corresponding to putative native state structures. We explore the consequences of our model. Our results suggest that symmetry and geometry are a powerful guide to capture the simplicity underlying protein complexity.

A geometrical framework for thinking about proteins.

We present a model, based on symmetry and geometry, for proteins. Using elementary ideas from mathematics and physics, we derive the geometries of discrete helices and sheets. We postulate a compatible solvent-mediated emergent pairwise attraction that assembles these building blocks, while respecting their individual symmetries. Instead of seeking to mimic the complexity of proteins, we look for a simple abstraction of reality that yet captures the essence of proteins. We employ analytic calculations and detailed Monte Carlo simulations to explore some consequences of our theory. The predictions of our approach are in accord with experimental data. Our framework provides a rationalization for understanding the common characteristics of proteins. Our results show that the free energy landscape of a globular protein is pre-sculpted at the backbone level, sequences and functionalities evolve in the fixed backdrop of the folds determined by geometry and symmetry, and that protein structures are unique in being simultaneously characterized by stability, diversity, and sensitivity.

DNA sequence symmetries from randomness: the origin of the Chargaff's second parity rule.

Most living organisms rely on double-stranded DNA (dsDNA) to store their genetic information and perpetuate themselves. This biological information has been considered as the main target of evolution. However, here we show that symmetries and patterns in the dsDNA sequence can emerge from the physical peculiarities of the dsDNA molecule itself and the maximum entropy principle alone, rather than from biological or environmental evolutionary pressure. The randomness justifies the human codon biases and context-dependent mutation patterns in human populations. Thus, the DNA 'exceptional symmetries,' emerged from the randomness, have to be taken into account when looking for the DNA encoded information. Our results suggest that the double helix energy constraints and, more generally, the physical properties of the dsDNA are the hard drivers of the overall DNA sequence architecture, whereas the selective biological processes act as soft drivers, which only under extraordinary circumstances overtake the overall entropy content of the genome.

Pairing statistics and melting of random DNA oligomers: Finding your partner in superdiverse environments.

Understanding of the pairing statistics in solutions populated by a large number of distinct solute species with mutual interactions is a challenging topic, relevant in modeling the complexity of real biological systems. Here we describe, both experimentally and theoretically, the formation of duplexes in a solution of random-sequence DNA (rsDNA) oligomers of length L = 8, 12, 20 nucleotides. rsDNA solutions are formed by 4L distinct molecular species, leading to a variety of pairing motifs that depend on sequence complementarity and range from strongly bound, fully paired defectless helices to weakly interacting mismatched duplexes. Experiments and theory coherently combine revealing a hybridization statistics characterized by a prevalence of partially defected duplexes, with a distribution of type and number of pairing errors that depends on temperature. We find that despite the enormous multitude of inter-strand interactions, defectless duplexes are formed, involving a fraction up to 15% of the rsDNA chains at the lowest temperatures. Experiments and theory are limited here to equilibrium conditions.

Phase behaviour of semiflexible lattice polymers in poor-solvent solution: Mean-field theory and Monte Carlo simulations.

We study a solution of interacting semiflexible polymers with curvature energy in poor-solvent conditions on the d-dimensional cubic lattice using mean-field theory and Monte Carlo computer simulations. Building upon past studies on a single chain, we construct a field-theory representation of the system and solve it within a mean-field approximation supported by Monte Carlo simulations in d = 3. A gas-liquid transition is found in the temperature-density plane that is then interpreted in terms of real systems. Interestingly, we find this transition to be independent of the bending rigidity. Past classical Flory-Huggins and Flory mean-field results are shown to be particular cases of this more general framework. Perspectives in terms of guiding experimental results towards optimal conditions are also proposed.

Effect of delay on the emergent stability patterns in generalized Lotka-Volterra ecological dynamics.

Understanding the conditions of feasibility and stability in ecological systems is a major challenge in theoretical ecology. The seminal work of May in 1972 and recent developments based on the theory of random matrices have shown the existence of emergent universal patterns of both stability and feasibility in ecological dynamics. However, only a few studies have investigated the role of delay coupled with population dynamics in the emergence of feasible and stable states. In this work, we study the effects of delay on generalized Loka-Volterra population dynamics of several interacting species in closed ecological environments. First, we investigate the relation between feasibility and stability of the modelled ecological community in the absence of delay and find a simple analytical relation when intra-species interactions are dominant. We then show how, by increasing the time delay, there is a transition in the stability phases of the population dynamics: from an equilibrium state to a stable non-point attractor phase. We calculate analytically the critical delay of that transition and show that it is in excellent agreement with numerical simulations. Finally, following a similar approach to characterizing stability in empirical studies, we investigate the coefficient of variation, which quantifies the magnitude of population fluctuations. We show that in the oscillatory regime induced by the delay, the variability at community level decreases for increasing diversity. This article is part of the theme issue 'Emergent phenomena in complex physical and socio-technical systems: from cells to societies'.

Generalized size scaling of metabolic rates based on single-cell measurements with freshwater phytoplankton.

Kleiber's law describes the scaling of metabolic rate with body size across several orders of magnitude in size and across taxa and is widely regarded as a fundamental law in biology. The physiological origins of Kleiber's law are still debated and generalizations of the law accounting for deviations from the scaling behavior have been proposed. Most theoretical and experimental studies of Kleiber's law, however, have focused on the relationship between the average body size of a species and its mean metabolic rate, neglecting intraspecific variation of these 2 traits. Here, we propose a theoretical characterization of such variation and report on proof-of-concept experiments with freshwater phytoplankton supporting such framework. We performed joint measurements at the single-cell level of cell volume and nitrogen/carbon uptake rates, as proxies of metabolic rates, of 3 phytoplankton species using nanoscale secondary ion mass spectrometry (NanoSIMS) and stable isotope labeling. Common scaling features of the distribution of nutrient uptake rates and cell volume are found to hold across 3 orders of magnitude in cell size. Once individual measurements of cell volume and nutrient uptake rate within a species are appropriately rescaled by a function of the average cell volume within each species, we find that intraspecific distributions of cell volume and metabolic rates collapse onto a universal curve. Based on the experimental results, this work provides the building blocks for a generalized form of Kleiber's law incorporating intraspecific, correlated variations of nutrient-uptake rates and body sizes.

Effective Resource Competition Model for Species Coexistence.

Local coexistence of species in large ecosystems is traditionally explained within the broad framework of niche theory. However, its rationale hardly justifies rich biodiversity observed in nearly homogeneous environments. Here we consider a consumer-resource model in which a coarse-graining procedure accounts for a variety of ecological mechanisms and leads to effective spatial effects which favor species coexistence. Herein, we provide conditions for several species to live in an environment with very few resources. In fact, the model displays two different phases depending on whether the number of surviving species is larger or smaller than the number of resources. We obtain conditions whereby a species can successfully colonize a pool of coexisting species. Finally, we analytically compute the distribution of the population sizes of coexisting species. Numerical simulations as well as empirical distributions of population sizes support our analytical findings.

Dynamic metabolic adaptation can promote species coexistence in competitive microbial communities.

Microbes are capable of physiologically adapting to diverse environmental conditions by differentially varying the rates at which they uptake different nutrients. In particular, microbes can switch hierarchically between different energy sources, consuming first those that ensure the highest growth rate. Experimentally, this can result in biphasic growth curves called "diauxic shifts" that typically arise when microbes are grown in media containing several nutrients. Despite these observations are well known in microbiology and molecular biology, the mathematical models generally used to describe the population dynamics of microbial communities do not account for dynamic metabolic adaptation, thus implicitly assuming that microbes cannot switch dynamically from one resource to another. Here, we introduce dynamic metabolic adaptation in the framework of consumer-resource models, which are commonly used to describe competitive microbial communities, allowing each species to temporally change its preferred energy source to maximize its own relative fitness. We show that dynamic metabolic adaptation enables the community to self-organize, allowing several species to coexist even in the presence of few resources, and to respond optimally to a time-dependent environment, thus showing that dynamic metabolic adaptation could be an important mechanism for maintaining high levels of diversity even in environments with few energy sources. We show that introducing dynamic metabolic strategies in consumer-resource models is necessary for reproducing experimental growth curves of the baker's yeast Saccharomyces cerevisiae growing in the presence of two carbon sources. Even though diauxic shifts emerge naturally from the model when two resources are qualitatively very different, the model predicts that the existence of such shifts is not a prerequisite for species coexistence in competitive communities.

River landscapes and optimal channel networks.

We study tree structures termed optimal channel networks (OCNs) that minimize the total gravitational energy loss in the system, an exact property of steady-state landscape configurations that prove dynamically accessible and strikingly similar to natural forms. Here, we show that every OCN is a so-called natural river tree, in the sense that there exists a height function such that the flow directions are always directed along steepest descent. We also study the natural river trees in an arbitrary graph in terms of forbidden substructures, which we call k-path obstacles, and OCNs on a d-dimensional lattice, improving earlier results by determining the minimum energy up to a constant factor for every [Formula: see text] Results extend our capabilities in environmental statistical mechanics.

Covariations in ecological scaling laws fostered by community dynamics.

Scaling laws in ecology, intended both as functional relationships among ecologically relevant quantities and the probability distributions that characterize their occurrence, have long attracted the interest of empiricists and theoreticians. Empirical evidence exists of power laws associated with the number of species inhabiting an ecosystem, their abundances, and traits. Although their functional form appears to be ubiquitous, empirical scaling exponents vary with ecosystem type and resource supply rate. The idea that ecological scaling laws are linked has been entertained before, but the full extent of macroecological pattern covariations, the role of the constraints imposed by finite resource supply, and a comprehensive empirical verification are still unexplored. Here, we propose a theoretical scaling framework that predicts the linkages of several macroecological patterns related to species' abundances and body sizes. We show that such a framework is consistent with the stationary-state statistics of a broad class of resource-limited community dynamics models, regardless of parameterization and model assumptions. We verify predicted theoretical covariations by contrasting empirical data and provide testable hypotheses for yet unexplored patterns. We thus place the observed variability of ecological scaling exponents into a coherent statistical framework where patterns in ecology embed constrained fluctuations.

Emergence of Network Motifs in Deep Neural Networks.

Network science can offer fundamental insights into the structural and functional properties of complex systems. For example, it is widely known that neuronal circuits tend to organize into basic functional topological modules, called network motifs. In this article, we show that network science tools can be successfully applied also to the study of artificial neural networks operating according to self-organizing (learning) principles. In particular, we study the emergence of network motifs in multi-layer perceptrons, whose initial connectivity is defined as a stack of fully-connected, bipartite graphs. Simulations show that the final network topology is shaped by learning dynamics, but can be strongly biased by choosing appropriate weight initialization schemes. Overall, our results suggest that non-trivial initialization strategies can make learning more effective by promoting the development of useful network motifs, which are often surprisingly consistent with those observed in general transduction networks.

The elixir phase of chain molecules.

A phase of matter is a familiar notion for inanimate physical matter. The nature of a phase of matter transcends the microscopic material properties. For example, materials in the liquid phase have certain common properties independent of the chemistry of the constituents: liquids take the shape of the container; they flow; and they can be poured-alcohol, oil, and water as well as a Lennard-Jones computer model exhibit similar behavior when poised in the liquid phase. Here, we identify a hitherto unstudied "phase" of matter, the elixir phase, in a simple model of a polymeric chain whose backbone has the correct local cylindrical symmetry induced by the tangent to the chain. The elixir phase appears on breaking the cylindrical symmetry by adding side spheres along the negative normal direction, as in proteins. This phase, nestled between other phases, has multiple ground states made up of building blocks of helices and almost planar sheets akin to protein native folds. We discuss the similarities of this "phase" of a finite size system to the liquid crystal and spin glass phases. Our findings are relevant for understanding proteins; the creation of novel bioinspired nanomachines; and also may have implications for life elsewhere in the cosmos.

Archetypes of human cognition defined by time preference for reward and their brain correlates: An evolutionary trade-off approach.

Biological systems carry out multiple tasks in their lifetime, which, in the course of evolution, may lead to trade-offs. In fact phenotypes (different species, individuals within a species, circuits, bacteria, proteins, etc.) cannot be optimal at all tasks, and, according to Pareto optimality theory, lay into a well-defined geometrical distribution (polygons and/or polyhedrons) in the space of traits. The vertices of this distribution contain archetypes, namely phenotypes that are specialists at one of the tasks, whereas phenotypes toward the center of the geometrical distribution show average performance across tasks. We applied this theory to the variability of cognitive and behavioral scores measured in 1206 individuals from the Human Connectome Project. Among all possible combinations of pairs of traits, we found the best fit to Pareto optimality when individuals were plotted in the trait-space of time preferences for reward, evaluated with the Delay Discounting Task (DDT). The DDT measures subjects' preference in choosing either immediate smaller rewards or delayed larger rewards. Time preference for reward was described by a triangular distribution in which each of the three vertices included individuals who used a particular strategy to discount reward. These archetypes accounted for variability on many cognitive, personality, and socioeconomic status variables, as well as differences in brain structure and functional connectivity, with only a weak influence of genetics. In summary, time preference for reward reflects a core variable that biases human phenotypes via natural and cultural selection.