Stochastic Turing patterns of trichomes in Arabidopsis leaves.

Biological patterns that emerge during the morphogenesis of multicellular organisms can display high precision at large scales, while at cellular scales, cells exhibit large fluctuations stemming from cell-cell differences in molecular copy numbers also called demographic noise. We study the conflicting interplay between high precision and demographic noise in trichome patterns on the epidermis of wild-type Arabidopsis thaliana leaves, as a two-dimensional model system. We carry out a statistical characterization of these patterns and show that their power spectra display fat tails-a signature compatible with noise-driven stochastic Turing patterns-which are absent in power spectra of patterns driven by deterministic instabilities. We then present a theoretical model that includes demographic noise stemming from birth-death processes of genetic regulators which we study analytically and by stochastic simulations. The model captures the observed experimental features of trichome patterns.

Deep learning and alignment of spatially resolved single-cell transcriptomes with Tangram.

Charting an organs' biological atlas requires us to spatially resolve the entire single-cell transcriptome, and to relate such cellular features to the anatomical scale. Single-cell and single-nucleus RNA-seq (sc/snRNA-seq) can profile cells comprehensively, but lose spatial information. Spatial transcriptomics allows for spatial measurements, but at lower resolution and with limited sensitivity. Targeted in situ technologies solve both issues, but are limited in gene throughput. To overcome these limitations we present Tangram, a method that aligns sc/snRNA-seq data to various forms of spatial data collected from the same region, including MERFISH, STARmap, smFISH, Spatial Transcriptomics (Visium) and histological images. Tangram can map any type of sc/snRNA-seq data, including multimodal data such as those from SHARE-seq, which we used to reveal spatial patterns of chromatin accessibility. We demonstrate Tangram on healthy mouse brain tissue, by reconstructing a genome-wide anatomically integrated spatial map at single-cell resolution of the visual and somatomotor areas.

Application limits of the scaling relations for Monte Carlo simulations in diffuse optics. Part 1: theory.

Monte Carlo (MC) is a powerful tool to study photon migration in scattering media, yet quite time-consuming to solve inverse problems. To speed up MC-simulations, scaling relations can be applied to an existing initial MC-simulation to generate a new data-set with different optical properties. We named this approach trajectory-based since it uses the knowledge of the detected photon trajectories of the initial MC-simulation, in opposition to the slower photon-based approach, where a novel MC-simulation is rerun with new optical properties. We investigated the convergence and applicability limits of the scaling relations, both related to the likelihood that the sample of trajectories considered is representative also for the new optical properties. For absorption, the scaling relation contains smoothly converging Lambert-Beer factors, whereas for scattering it is the product of two quickly diverging factors, whose ratio, for NIRS cases, can easily reach ten orders of magnitude. We investigated such instability by studying the probability-distribution for the number of scattering events in trajectories of given length. We propose a convergence test of the scattering scaling relation based on the minimum-maximum number of scattering events in recorded trajectories. We also studied the dependence of MC-simulations on optical properties, most critical in inverse problems, finding that scattering derivatives are ascribed to small deviations in the distribution of scattering events from a Poisson distribution. This paper, which can also serve as a tutorial, helps to understand the physics of the scaling relations with the causes of their limitations and devise new strategies to deal with them.

Cortical propagation tracks functional recovery after stroke.

Stroke is a debilitating condition affecting millions of people worldwide. The development of improved rehabilitation therapies rests on finding biomarkers suitable for tracking functional damage and recovery. To achieve this goal, we perform a spatiotemporal analysis of cortical activity obtained by wide-field calcium images in mice before and after stroke. We compare spontaneous recovery with three different post-stroke rehabilitation paradigms, motor training alone, pharmacological contralesional inactivation and both combined. We identify three novel indicators that are able to track how movement-evoked global activation patterns are impaired by stroke and evolve during rehabilitation: the duration, the smoothness, and the angle of individual propagation events. Results show that, compared to pre-stroke conditions, propagation of cortical activity in the subacute phase right after stroke is slowed down and more irregular. When comparing rehabilitation paradigms, we find that mice treated with both motor training and pharmacological intervention, the only group associated with generalized recovery, manifest new propagation patterns, that are even faster and smoother than before the stroke. In conclusion, our new spatiotemporal propagation indicators could represent promising biomarkers that are able to uncover neural correlates not only of motor deficits caused by stroke but also of functional recovery during rehabilitation. In turn, these insights could pave the way towards more targeted post-stroke therapies.

Reconstruction scheme for excitatory and inhibitory dynamics with quenched disorder: application to zebrafish imaging.

An inverse procedure is developed and tested to recover functional and structural information from global signals of brains activity. The method assumes a leaky-integrate and fire model with excitatory and inhibitory neurons, coupled via a directed network. Neurons are endowed with a heterogenous current value, which sets their associated dynamical regime. By making use of a heterogenous mean-field approximation, the method seeks to reconstructing from global activity patterns the distribution of in-coming degrees, for both excitatory and inhibitory neurons, as well as the distribution of the assigned currents. The proposed inverse scheme is first validated against synthetic data. Then, time-lapse acquisitions of a zebrafish larva recorded with a two-photon light sheet microscope are used as an input to the reconstruction algorithm. A power law distribution of the in-coming connectivity of the excitatory neurons is found. Local degree distributions are also computed by segmenting the whole brain in sub-regions traced from annotated atlas.

Robust stochastic Turing patterns in the development of a one-dimensional cyanobacterial organism.

Under nitrogen deprivation, the one-dimensional cyanobacterial organism Anabaena sp. PCC 7120 develops patterns of single, nitrogen-fixing cells separated by nearly regular intervals of photosynthetic vegetative cells. We study a minimal, stochastic model of developmental patterns in Anabaena that includes a nondiffusing activator, two diffusing inhibitor morphogens, demographic fluctuations in the number of morphogen molecules, and filament growth. By tracking developing filaments, we provide experimental evidence for different spatiotemporal roles of the two inhibitors during pattern maintenance and for small molecular copy numbers, justifying a stochastic approach. In the deterministic limit, the model yields Turing patterns within a region of parameter space that shrinks markedly as the inhibitor diffusivities become equal. Transient, noise-driven, stochastic Turing patterns are produced outside this region, which can then be fixed by downstream genetic commitment pathways, dramatically enhancing the robustness of pattern formation, also in the biologically relevant situation in which the inhibitors' diffusivities may be comparable.

Analysis and forecast of COVID-19 spreading in China, Italy and France.

In this note we analyze the temporal dynamics of the coronavirus disease 2019 outbreak in China, Italy and France in the time window 22 / 01 - 15 / 03 / 2020 . A first analysis of simple day-lag maps points to some universality in the epidemic spreading, suggesting that simple mean-field models can be meaningfully used to gather a quantitative picture of the epidemic spreading, and notably the height and time of the peak of confirmed infected individuals. The analysis of the same data within a simple susceptible-infected-recovered-deaths model indicates that the kinetic parameter that describes the rate of recovery seems to be the same, irrespective of the country, while the infection and death rates appear to be more variable. The model places the peak in Italy around March 21st 2020, with a peak number of infected individuals of about 26000 (not including recovered and dead) and a number of deaths at the end of the epidemics of about 18,000. Since the confirmed cases are believed to be between 10 and 20% of the real number of individuals who eventually get infected, the apparent mortality rate of COVID-19 falls between 4% and 8% in Italy, while it appears substantially lower, between 1% and 3% in China. Based on our calculations, we estimate that 2500 ventilation units should represent a fair figure for the peak requirement to be considered by health authorities in Italy for their strategic planning. Finally, a simulation of the effects of drastic containment measures on the outbreak in Italy indicates that a reduction of the infection rate indeed causes a quench of the epidemic peak. However, it is also seen that the infection rate needs to be cut down drastically and quickly to observe an appreciable decrease of the epidemic peak and mortality rate. This appears only possible through a concerted and disciplined, albeit painful, effort of the population as a whole.

Patterns of non-normality in networked systems.

Several mechanisms have been proposed to explain the spontaneous generation of self-organised patterns, hypothesised to play a role in the formation of many of the magnificent patterns observed in Nature. In several cases of interest, the system under scrutiny displays a homogeneous equilibrium, which is destabilised via a symmetry breaking instability which reflects the specificity of the problem being inspected. The Turing instability is among the most celebrated paradigms for pattern formation. In its original form, the diffusion constants of the two mobile species need to be quite different from each other for the instability to develop. Unfortunately, this condition limits the applicability of the theory. To overcome this impediment, and with the ambitious long term goal to eventually reconcile theory and experiments, we here propose an alternative mechanism for promoting the onset of pattern. To this end a multi-species reactive model is studied, assuming a generalized transport on a discrete and directed network-like support: the instability is triggered by the non-normality of the embedding network. The non-normal character of the dynamics instigates a short time amplification of the imposed perturbation, thus making the system unstable for a choice of parameters that would yield stability under the conventional scenario. In other words, non-normality promotes the emergence of patterns in cases where a classical linear analysis would not predict them. The importance of our result relies also on the fact that non-normal networks are pervasively found, motivating the general interest of the mechanism discussed here.

Resilience for stochastic systems interacting via a quasi-degenerate network.

A stochastic reaction-diffusion model is studied on a networked support. In each patch of the network, two species are assumed to interact following a non-normal reaction scheme. When the interaction unit is replicated on a directed linear lattice, noise gets amplified via a self-consistent process, which we trace back to the degenerate spectrum of the embedding support. The same phenomenon holds when the system is bound to explore a quasidegenerate network. In this case, the eigenvalues of the Laplacian operator, which governs species diffusion, accumulate over a limited portion of the complex plane. The larger the network, the more pronounced the amplification. Beyond a critical network size, a system deemed deterministically stable, hence resilient, can develop seemingly regular patterns in the concentration amount. Non-normality and quasidegenerate networks may, therefore, amplify the inherent stochasticity and so contribute to altering the perception of resilience, as quantified via conventional deterministic methods.

Hopping in the Crowd to Unveil Network Topology.

We introduce a nonlinear operator to model diffusion on a complex undirected network under crowded conditions. We show that the asymptotic distribution of diffusing agents is a nonlinear function of the nodes' degree and saturates to a constant value for sufficiently large connectivities, at variance with standard diffusion in the absence of excluded-volume effects. Building on this observation, we define and solve an inverse problem, aimed at reconstructing the a priori unknown connectivity distribution. The method gathers all the necessary information by repeating a limited number of independent measurements of the asymptotic density at a single node, which can be chosen randomly. The technique is successfully tested against both synthetic and real data and is also shown to estimate with great accuracy the total number of nodes.

Theory of Turing Patterns on Time Varying Networks.

The process of pattern formation for a multispecies model anchored on a time varying network is studied. A nonhomogeneous perturbation superposed to an homogeneous stable fixed point can be amplified following the Turing mechanism of instability, solely instigated by the network dynamics. By properly tuning the frequency of the imposed network evolution, one can make the examined system behave as its averaged counterpart, over a finite time window. This is the key observation to derive a closed analytical prediction for the onset of the instability in the time dependent framework. Continuously and piecewise constant periodic time varying networks are analyzed, setting the framework for the proposed approach. The extension to nonperiodic settings is also discussed.

Correction: Theory of diffusion-influenced reactions in complex geometries.

Correction for 'Theory of diffusion-influenced reactions in complex geometries' by Marta Galanti et al., Phys. Chem. Chem. Phys., 2016, DOI: .

Theory of diffusion-influenced reactions in complex geometries.

Chemical transformations involving the diffusion of reactants and subsequent chemical fixation steps are generally termed "diffusion-influenced reactions" (DIR). Virtually all biochemical processes in living media can be counted among them, together with those occurring in an ever-growing number of emerging nano-technologies. The role of the environment's geometry (obstacles, compartmentalization) and distributed reactivity (competitive reactants, traps) is key in modulating the rate constants of DIRs, and is therefore a prime design parameter. Yet, it is a formidable challenge to build a comprehensive theory that is able to describe the environment's "reactive geometry". Here we show that such a theory can be built by unfolding this many-body problem through addition theorems for special functions. Our method is powerful and general and allows one to study a given DIR reaction occurring in arbitrary "reactive landscapes", made of multiple spherical boundaries of given size and reactivity. Importantly, ready-to-use analytical formulas can be derived easily in most cases.

Reaction rate of a composite core-shell nanoreactor with multiple nanocatalysts.

We present a detailed theory for the total reaction rate constant of a composite core-shell nanoreactor, consisting of a central solid core surrounded by a hydrogel layer of variable thickness, where a given number of small catalytic nanoparticles are embedded at prescribed positions and are endowed with a prescribed surface reaction rate constant. Besides the precise geometry of the assembly, our theory accounts explicitly for the diffusion coefficients of the reactants in the hydrogel and in the bulk as well as for their transfer free energy jump upon entering the hydrogel shell. Moreover, we work out an approximate analytical formula for the overall rate constant, which is valid in the physically relevant range of geometrical and chemical parameters. We discuss in depth how the diffusion-controlled part of the rate depends on the essential variables, including the size of the central core. In particular, we derive some simple rules for estimating the number of nanocatalysts per nanoreactor for an efficient catalytic performance in the case of small to intermediate core sizes. Our theoretical treatment promises to provide a very useful and flexible tool for the design of superior performing nanoreactor geometries with optimized nanoparticle load.

Linear noise approximation for stochastic oscillations of intracellular calcium.

A stochastic model of intracellular calcium oscillations is analytically studied. The governing master equation is expanded under the linear noise approximation and a closed prediction for the power spectrum of fluctuations analytically derived. A peak in the obtained power spectrum profile signals the presence of stochastic, noise induced oscillations which extend also outside the region where a deterministic limit cycle is predicted to occur.

Stochastic amplification of spatial modes in a system with one diffusing species.

The problem of pattern formation in a generic two species reaction-diffusion model is studied, under the hypothesis that only one species can diffuse. For such a system, the classical Turing instability cannot take place. At variance, by working in the generalized setting of a stochastic formulation to the inspected problem, spatially organized patterns can develop, seeded by finite size corrections. General conditions are given for the stochastic patterns to occur. The predictions of the theory are tested for a specific case study.

Force and kinetics of fast and slow muscle myosin determined with a synthetic sarcomere-like nanomachine.

Myosin II is the muscle molecular motor that works in two bipolar arrays in each thick filament of the striated (skeletal and cardiac) muscle, converting the chemical energy into steady force and shortening by cyclic ATP-driven interactions with the nearby actin filaments. Different isoforms of the myosin motor in the skeletal muscles account for the different functional requirements of the slow muscles (primarily responsible for the posture) and fast muscles (responsible for voluntary movements). To clarify the molecular basis of the differences, here the isoform-dependent mechanokinetic parameters underpinning the force of slow and fast muscles are defined with a unidimensional synthetic nanomachine powered by pure myosin isoforms from either slow or fast rabbit skeletal muscle. Data fitting with a stochastic model provides a self-consistent estimate of all the mechanokinetic properties of the motor ensemble including the motor force, the fraction of actin-attached motors and the rate of transition through the attachment-detachment cycle. The achievements in this paper set the stage for any future study on the emergent mechanokinetic properties of an ensemble of myosin molecules either engineered or purified from mutant animal models or human biopsies.

Integrins regulate hERG1 dynamics by girdin-dependent Gαi3: signaling and modeling in cancer cells.

The hERG1 potassium channel is aberrantly over expressed in tumors and regulates the cancer cell response to integrin-dependent adhesion. We unravel a novel signaling pathway by which integrin engagement by the ECM protein fibronectin promotes hERG1 translocation to the plasma membrane and its association with ß1 integrins, by activating girdin-dependent Gαi3 proteins and protein kinase B (Akt). By sequestering hERG1, ß1 integrins make it avoid Rab5-mediated endocytosis, where unbound channels are degraded. The cycle of hERG1 expression determines the resting potential (Vrest) oscillations and drives the cortical f-actin dynamics and thus cell motility. To interpret the slow biphasic kinetics of hERG1/ß1 integrin interplay, we developed a mathematical model based on a generic balanced inactivation-like module. Integrin-mediated cell adhesion triggers two contrary responses: a rapid stimulation of hERG1/ß1 complex formation, followed by a slow inhibition which restores the initial condition. The protracted hERG1/ß1 integrin cycle determines the slow time course and cyclic behavior of cell migration in cancer cells.

Brain-wide neuron quantification toolkit reveals strong sexual dimorphism in the evolution of fear memory.

Fear responses are functionally adaptive behaviors that are strengthened as memories. Indeed, detailed knowledge of the neural circuitry modulating fear memory could be the turning point for the comprehension of this emotion and its pathological states. A comprehensive understanding of the circuits mediating memory encoding, consolidation, and retrieval presents the fundamental technological challenge of analyzing activity in the entire brain with single-neuron resolution. In this context, we develop the brain-wide neuron quantification toolkit (BRANT) for mapping whole-brain neuronal activation at micron-scale resolution, combining tissue clearing, high-resolution light-sheet microscopy, and automated image analysis. The robustness and scalability of this method allow us to quantify the evolution of activity patterns across multiple phases of memory in mice. This approach highlights a strong sexual dimorphism in recruited circuits, which has no counterpart in the behavior. The methodology presented here paves the way for a comprehensive characterization of the evolution of fear memory.

Prestrain-induced contraction in one-dimensional random elastic chains.

Prestrained elastic networks arise in a number of biological and technological systems ranging from the cytoskeleton of cells to tensegrity structures. Motivated by this observation, we here consider a minimal model in one dimension to set the stage for understanding the response of such networks as a function of the prestrain. To this end we consider a chain [one-dimensional (1D) network] of elastic springs upon which a random, zero mean, finite variance prestrain is imposed. Numerical simulations and analytical predictions quantify the magnitude of the contraction as a function of the variance of the prestrain, and show that the chain always shrinks. To test these predictions, we vary the topology of the chain, consider more complex connectivity and show that our results are relatively robust to these changes.