Search

Vagal sensory neurons monitor mechanical and chemical stimuli in the gastrointestinal tract. Major efforts are underway to assign physiological functions to the many distinct subtypes of vagal sensory neurons. Here, we use genetically guided anatomical tracing, optogenetics, and electrophysiology to identify and characterize vagal sensory neuron subtypes expressing Prox2 and Runx3 in mice. We show that three of these neuronal subtypes innervate the esophagus and stomach in regionalized patterns, where they form intraganglionic laminar endings. Electrophysiological analysis revealed that they are low-threshold mechanoreceptors but possess different adaptation properties. Lastly, genetic ablation of Prox2 and Runx3 neurons demonstrated their essential roles for esophageal peristalsis in freely behaving mice. Our work defines the identity and function of the vagal neurons that provide mechanosensory feedback from the esophagus to the brain and could lead to better understanding and treatment of esophageal motility disorders.

Minerals are solid state nuclear track detectors - nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred degrees C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials at room temperature). The damage features from the $O(50)$ MeV fission fragments left by spontaneous fission of $^{238}$U and other heavy unstable isotopes have long been used for fission track dating of geological samples. Laboratory studies have demonstrated the readout of defects caused by nuclear recoils with energies as small as $O(1)$ keV. This whitepaper discusses a wide range of possible applications of minerals as detectors for $E_R \gtrsim O(1)$ keV nuclear recoils: Using natural minerals, one could use the damage features accumulated over $O(10)$ Myr$-O(1)$ Gyr to measure astrophysical neutrino fluxes (from the Sun, supernovae, or cosmic rays interacting with the atmosphere) as well as search for Dark Matter. Using signals accumulated over months to few-years timescales in laboratory-manufactured minerals, one could measure reactor neutrinos or use them as Dark Matter detectors, potentially with directional sensitivity. Research groups in Europe, Asia, and America have started developing microscopy techniques to read out the $O(1) - O(100)$ nm damage features in crystals left by $O(0.1) - O(100)$ keV nuclear recoils. We report on the status and plans of these programs. The research program towards the realization of such detectors is highly interdisciplinary, combining geoscience, material science, applied and fundamental physics with techniques from quantum information and Artificial Intelligence., Comment: 115 pages, many pictures of tracks. Please see the source file for higher resolution versions of some plots. v2: matches the published version

SummarySensory neurons of the vagus nerve monitor distention and stretch in the gastrointestinal tract. Major efforts are underway to assign physiological functions to the many distinct subtypes of vagal sensory neurons. Here, we used genetically guided anatomical tracing, optogenetics and electrophysiology to identify and characterize three vagal sensory neuronal subtypes expressing Prox2 and Runx3. We show that these neuronal subtypes innervate the esophagus and stomach where they display regionalized innervation patterns. Their electrophysiological analysis showed that they are all low threshold mechanoreceptors, but possess different adaptation properties. Lastly, genetic ablation of Prox2+ and Runx3+ neurons demonstrated their essential roles for esophageal peristalsis and swallowing in freely behaving animals. Our work reveals the identity and function of the vagal neurons that provide mechanosensory feedback from the esophagus to the brain, and could lead to better understanding and treatment of esophageal motility disorders.

Here we discuss barriers to reproducibility in regard to microscopes and related hardware, along with best practices for sharing novel designs created using computer-aided design (CAD). We hope to start a fruitful community discussion on how instrument development, especially in microscopy, can become more open and reproducible, ultimately leading to better, more trustworthy science.

Immunotherapy constitutes a paradigm shift in cancer treatment. Its FDA approval for several indications has yielded improved prognosis for cases where traditional therapy has shown limited efficiency. However, many patients still fail to benefit from this treatment modality, and the exact mechanisms responsible for tumor response are unknown. Noninvasive treatment monitoring is crucial for longitudinal tumor characterization and the early detection of non-responders. While various medical imaging techniques can provide a morphological picture of the lesion and its surrounding tissue, a molecular-oriented imaging approach holds the key to unraveling biological effects that occur much earlier in the immunotherapy timeline. Magnetic resonance imaging (MRI) is a highly versatile imaging modality, where the image contrast can be tailored to emphasize a particular biophysical property of interest using advanced engineering of the imaging pipeline. In this review, recent advances in molecular-MRI based cancer immunotherapy monitoring are described. Next, the presentation of the underlying physics, computational, and biological features are complemented by a critical analysis of the results obtained in preclinical and clinical studies. Finally, emerging artificial intelligence (AI)-based strategies to further distill, quantify, and interpret the image-based molecular MRI information are discussed in terms of perspectives for the future.

Funding: This work was supported by the Howard Hughes Medical Institute and by the Simons Foundation (Simons Collaboration on the Global Brain #542943SPI). To track and control self-location, animals integrate their movements through space. Representations of self-location are observed in the mammalian hippocampal formation, but it is unknown if positional representations exist in more ancient brain regions, how they arise from integrated self-motion, and by what pathways they control locomotion. Here, in a head-fixed, fictive-swimming, virtual-reality preparation, we exposed larval zebrafish to a variety of involuntary displacements. They tracked these displacements and, many seconds later, moved toward their earlier location through corrective swimming (“positional homeostasis”). Whole-brain functional imaging revealed a network in the medulla that stores a memory of location and induces an error signal in the inferior olive to drive future corrective swimming. Optogenetically manipulating medullary integrator cells evoked displacement-memory behavior. Ablating them, or downstream olivary neurons, abolished displacement corrections. These results reveal a multiregional hindbrain circuit in vertebrates that integrates self-motion and stores self-location to control locomotor behavior. Publisher PDF

To accurately track self-location, animals need to integrate their movements through space. In amniotes, representations of self-location have been found in regions such as the hippocampus. It is unknown whether more ancient brain regions contain such representations and by which pathways they may drive locomotion. Fish displaced by water currents must prevent uncontrolled drift to potentially dangerous areas. We found that larval zebrafish track such movements and can later swim back to their earlier location. Whole-brain functional imaging revealed the circuit enabling this process of positional homeostasis. Position-encoding brainstem neurons integrate optic flow, then bias future swimming to correct for past displacements by modulating inferior olive and cerebellar activity. Manipulation of position-encoding or olivary neurons abolished positional homeostasis or evoked behavior as if animals had experienced positional shifts. These results reveal a multiregional hindbrain circuit in vertebrates for optic flow integration, memory of self-location, and its neural pathway to behavior.

Light-sheet microscopy has become one of the primary tools for imaging live developing organisms because of its high speed, low phototoxicity, and optical sectioning capabilities. Detection from multiple sides (multi-view imaging) additionally allows nearly isotropic resolution via computational merging of the views. However, conventional light-sheet microscopes require that the sample is suspended in a gel to allow optical access from two or more sides. At the same time, the use of microfluidic devices is highly desirable for many experiments, but geometric constrains and strong optical aberrations caused by the coverslip titled relative to objectives make the use of multi-view lightsheet challenging for microfluidics.In this paper we describe the use of adaptive optics (AO) to enable multi-view light-sheet microscopy in such microfluidic setup by correcting optical aberrations introduced by the tilted coverslip. The optimal shape of deformable mirror is computed by an iterative stochastic gradient-descent algorithm that optimizes PSF in two orthogonal planes simultaneously. Simultaneous AO correction in two optical arms is achieved via a knife-edge mirror that splits excitation path and combines the detection path.We characterize the performance of this novel microscope setup and, by dual-view light-sheet imaging of C.elegans inside a microfluidic channel, demonstrate a drastic improvement of image quality due to AO and dual-view reconstruction. Our microscope design allows multi-view light-sheet microscopy with microfluidic devices for precisely controlled experimental conditions and high-content screening.

Whole-brain imaging allows for comprehensive functional mapping of distributed neural pathways, but neuronal perturbation experiments are usually limited to targeting predefined regions or genetically identifiable cell types. To complement whole-brain measures of activity with brain-wide manipulations for testing causal interactions, we introduce a system that uses measured activity patterns to guide optical perturbations of any subset of neurons in the same fictively behaving larval zebrafish. First, a light-sheet microscope collects whole-brain data that are rapidly analyzed by a distributed computing system to generate functional brain maps. On the basis of these maps, the experimenter can then optically ablate neurons and image activity changes across the brain. We applied this method to characterize contributions of behaviorally tuned populations to the optomotor response. We extended the system to optogenetically stimulate arbitrary subsets of neurons during whole-brain imaging. These open-source methods enable delineating the contributions of neurons to brain-wide circuit dynamics and behavior in individual animals.

Fast ripples (FRs) are network oscillations, defined variously as having frequencies of > 150 to > 250 Hz, with a controversial mechanism. FRs appear to indicate a propensity of cortical tissue to originate seizures. Here, we demonstrate field oscillations, at up to 400 Hz, in spontaneously epileptic human cortical tissue in vitro, and present a network model that could explain FRs themselves, and their relation to ‘ordinary’ (slower) ripples. We performed network simulations with model pyramidal neurons, having axons electrically coupled. Ripples (< 250 Hz) were favored when conduction of action potentials, axon to axon, was reliable. Whereas ripple population activity was periodic, firing of individual axons varied in relative phase. A switch from ripples to FRs took place when an ectopic spike occurred in a cell coupled to another cell, itself multiply coupled to others. Propagation could then start in one direction only, a condition suitable for re-entry. The resulting oscillations were > 250 Hz, were sustained or interrupted, and had little jitter in the firing of individual axons. The form of model FR was similar to spontaneously occurring FRs in excised human epileptic tissue. In vitro, FRs were suppressed by a gap junction blocker. Our data suggest that a given network can produce ripples, FRs, or both, via gap junctions, and that FRs are favored by clusters of axonal gap junctions. If axonal gap junctions indeed occur in epileptic tissue, and are mediated by connexin 26 (recently shown to mediate coupling between immature neocortical pyramidal cells), then this prediction is testable.

Mechanisms of place cell replay occurring during sharp-wave ripples (SPW-Rs) remain obscure due to the fact that ripples in vitro depend on non-synaptic mechanisms, presumably via axo-axonal gap junctions between pyramidal cells. We suggest a model of in vivo SPW-Rs in which synaptic excitatory post-synaptic potentials (EPSPs) control the axonal spiking of cells in SPW-Rs: ripple activity remains hidden in the network of axonal collaterals (connected by gap junctions) due to conduction failures, unless there is a sufficient dendritic EPSP. The EPSP brings the axonal branching point to threshold, and action potentials from the collateral start to propagate to the soma and to the distal axon. The model coherently explains multiple experimental data on SPW-Rs, both in vitro and in vivo. The mechanism of synaptic gating leads to the following implication: a sequence of pyramidal cells can be replayed at ripple frequency by the superposition of subthreshold dendritic EPSPs and ripple activity in the axonal plexus. Replay is demonstrated in both forward and reverse directions. We discuss several testable predictions. In general, the mechanism of synaptic gating suggests that pyramidal cells under certain conditions can act like a transistor.

Bacterial chemotaxis represents one of the simplest and best studied examples of unicellular behavior. Chemotaxis allows swimming bacterial cells to follow chemical gradients in the environment by performing temporal comparisons of ligand concentrations. The process of chemotaxis in the model bacterium Escherichia coli has been studied in great molecular detail over the past 40 years, using a large range of experimental tools to investigate physiology, genetics and biochemistry of the system. The abundance of quantitative experimental data enabled detailed computational modeling of the pathway and theoretical analyses of such properties as robustness and signal amplification. Because of the temporal mode of gradient sensing in bacterial chemotaxis, molecular memory is an essential component of the chemotaxis pathway. Recent studies suggest that the memory time scale has been evolutionary optimized to perform optimal comparisons of stimuli while swimming in the gradient. Moreover, noise in the adaptation system, which results from variations of the adaptation rate both over time and among cells, might be beneficial for the overall chemotactic performance of the population.

Gene expression is known to correlate with the degree of codon bias in many unicellular organisms. However, such a correlation is not observed in some organisms. It was demonstrated that inverted complementary repeats within coding DNA sequences (ORFs) should be considered for proper estimation of the translation efficiency because they can form secondary structures that obstruct ribosome movement. A program was developed for estimating the potential expression of ORFs in unicellular organisms on the basis of their genome sequences. The program computes the elongation efficiency index (EEI) and takes into account three key factors: codon bias, the average number of inverted complementary repeats, and the free energies of potential stem-loop structures formed by these repeats. The influence of these factors on translation was numerically estimated. Their optimal ratio was computed for each organism. EEIs of 384 unicellular organisms (351 bacteria, 28 archaea, and 5 eukaryotes) were computed using the annotated genomes available from GenBank. Five potential evolutionary strategies of translational optimization were determined in the organisms studied. A considerable difference in preferential translational strategies was observed between bacteria and archaea. Significant correlations between EEIs and gene expression levels were shown for two species (Saccharomyces cerevisiae and Helicobacter pylori), using the available microarray data. The method allows the numerical estimation of the translation efficiency of an ORF and optimization of the nucleotide composition of heterologous genes in specified unicellular organisms. The program is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/eei-calculator.

Understanding brain function requires monitoring and interpreting the activity of large networks of neurons during behavior. Advances in recording technology are greatly increasing the size and complexity of neural data. Analyzing such data will pose a fundamental bottleneck for neuroscience. We present a library of analytical tools called Thunder built on the open-source Apache Spark platform for large-scale distributed computing. The library implements a variety of univariate and multivariate analyses with a modular, extendable structure well-suited to interactive exploration and analysis development. We demonstrate how these analyses find structure in large-scale neural data, including whole-brain light-sheet imaging data from fictively behaving larval zebrafish, and two-photon imaging data from behaving mouse. The analyses relate neuronal responses to sensory input and behavior, run in minutes or less and can be used on a private cluster or in the cloud. Our open-source framework thus holds promise for turning brain activity mapping efforts into biological insights.

Increasing the volumetric imaging speed of light-sheet microscopy will improve its ability to detect fast changes in neural activity. Here, a system is introduced for brain-wide imaging of neural activity in the larval zebrafish by coupling structured illumination with cubic phase extended depth-of-field (EDoF) pupil encoding. This microscope enables faster light-sheet imaging and facilitates arbitrary plane scanning—removing constraints on acquisition speed, alignment tolerances, and physical motion near the sample. The usefulness of this method is demonstrated by performing multi-plane calcium imaging in the fish brain with a 416 × 832 × 160 µm field of view at 33 Hz. The optomotor response behavior of the zebrafish is monitored at high speeds, and time-locked correlations of neuronal activity are resolved across its brain.

High-frequency oscillations (HFOs) are an important part of brain activity in health and disease. However, their origins remain obscure and controversial. One possible mechanism depends on the presence of sparsely distributed gap junctions that electrically couple the axons of principal cells. A plexus of electrically coupled axons is modeled as a random network with bidirectional connections between its nodes. Under certain conditions the network can demonstrate one of two types of oscillatory activity. Type I oscillations (100-200 Hz) are predicted to be caused by spontaneously spiking axons in a network with strong (high-conductance) gap junctions. Type II oscillations (200-300 Hz) require no spontaneous spiking and relatively weak (low-conductance) gap junctions, across which spike propagation failures occur. The type II oscillations are reentrant and self-sustained. Here we examine what determines the frequency of type II oscillations. Using simulations we show that the distribution of loop lengths is the key factor for determining frequency in type II network oscillations. We first analyze spike failure between two electrically coupled cells using a model of anatomically reconstructed CA1 pyramidal neuron. Then network oscillations are studied by a cellular automaton model with random network connectivity, in which we control loop statistics. We show that oscillation periods can be predicted from the network's loop statistics. The shortest loop, around which a spike can travel, is the most likely pacemaker candidate.The principle of one loop as a pacemaker is remarkable, because random networks contain a large number of loops juxtaposed and superimposed, and their number rapidly grows with network size. This principle allows us to predict the frequency of oscillations from network connectivity and visa versa. We finally propose that type I oscillations may correspond to ripples, while type II oscillations correspond to so-called fast ripples.

Very fast oscillations (VFO) in neocortex are widely observed before epileptic seizures, and there is growing evidence that they are caused by networks of pyramidal neurons connected by gap junctions between their axons. We are motivated by the spatio-temporal waves of activity recorded using electrocorticography (ECoG), and study the speed of activity propagation through a network of neurons axonally coupled by gap junctions. We simulate wave propagation by excitable cellular automata (CA) on random (Erdös-Rényi) networks of special type, with spatially constrained connections. From the cellular automaton model, we derive a mean field theory to predict wave propagation. The governing equation resolved by the Fisher-Kolmogorov PDE fails to describe wave speed. A new (hyperbolic) PDE is suggested, which provides adequate wave speed v() that saturates with network degree , in agreement with intuitive expectations and CA simulations. We further show that the maximum length of connection is a much better predictor of the wave speed than the mean length. When tested in networks with various degree distributions, wave speeds are found to strongly depend on the ratio of network moments / rather than on mean degree , which is explained by general network theory. The wave speeds are strikingly similar in a diverse set of networks, including regular, Poisson, exponential and power law distributions, supporting our theory for various network topologies. Our results suggest practical predictions for networks of electrically coupled neurons, and our mean field method can be readily applied for a wide class of similar problems, such as spread of epidemics through spatial networks.

Chemotactic movement of Escherichia coli is one of the most thoroughly studied paradigms of simple behavior. Due to significant competitive advantage conferred by chemotaxis and to high evolution rates in bacteria, the chemotaxis system is expected to be strongly optimized. Bacteria follow gradients by performing temporal comparisons of chemoeffector concentrations along their runs, a strategy which is most efficient given their size and swimming speed. Concentration differences are detected by a sensory system and transmitted to modulate rotation of flagellar motors, decreasing the probability of a tumble and reorientation if the perceived concentration change during a run is positive. Such regulation of tumble probability is of itself sufficient to explain chemotactic drift of a population up the gradient, and is commonly assumed to be the only navigation mechanism of chemotactic E. coli. Here we use computer simulations to predict existence of an additional mechanism of gradient navigation in E. coli. Based on the experimentally observed dependence of cell tumbling angle on the number of switching motors, we suggest that not only the tumbling probability but also the degree of reorientation during a tumble depend on the swimming direction along the gradient. Although the difference in mean tumbling angles up and down the gradient predicted by our model is small, it results in a dramatic enhancement of the cellular drift velocity along the gradient. We thus demonstrate a new level of optimization in E. coli chemotaxis, which arises from the switching of several flagellar motors and a resulting fine tuning of tumbling angle. Similar strategy is likely to be used by other peritrichously flagellated bacteria, and indicates yet another level of evolutionary development of bacterial chemotaxis., Author Summary Chemotaxis of bacteria plays an important role in their life, providing them with the ability to actively search for an optimal growth environment. The chemotaxis system is supposed to be highly optimized, because on the evolutionary time scale even a modest enhancement of its efficiency can give cells a large competitive advantage. For a long time it was believed that the only navigation mechanism of bacteria is increasing the run length toward the preferred direction. The tumble was assumed to be a purely random change of direction between runs. We analysed recently published experimental data that demonstrate a dependence of tumbling angle on the number of CW-switched motors. We introduced such a dependence into our model of chemotactic E. coli, and simulated it under different conditions. Our simulations show that this dependence is an important additional mechanism of bacterial navigation, which was previously unrecognized because it lays below the experimental errors of conventional single-cell tracking. We show that such a fine tuning of tumbling significantly improves efficiency of chemotaxis, and represents a new level of evolutionary optimization of bacteria.

Adaptability is an essential property of many sensory systems, enabling maintenance of a sensitive response over a range of background stimulus levels. In bacterial chemotaxis, adaptation to the preset level of pathway activity is achieved through an integral feedback mechanism based on activity-dependent methylation of chemoreceptors. It has been argued that this architecture ensures precise and robust adaptation regardless of the ambient ligand concentration, making perfect adaptation a celebrated property of the chemotaxis system. However, possible deviations from such ideal adaptive behavior and its consequences for chemotaxis have not been explored in detail. Here we show that the chemotaxis pathway in Escherichia coli shows increasingly imprecise adaptation to higher concentrations of attractants, with a clear correlation between the time of adaptation to a step-like stimulus and the extent of imprecision. Our analysis suggests that this imprecision results from a gradual saturation of receptor methylation sites at high levels of stimulation, which prevents full recovery of the pathway activity by violating the conditions required for precise adaptation. We further use computer simulations to show that limited imprecision of adaptation has little effect on the rate of chemotactic drift of a bacterial population in gradients, but hinders precise accumulation at the peak of the gradient. Finally, we show that for two major chemoeffectors, serine and cysteine, failure of adaptation at concentrations above 1 mM might prevent bacteria from accumulating at toxic concentrations of these amino acids.

Evolutionary selection for robustness of signaling output in the face of stochastic variations in protein expression may explain the organization of bacterial chemotaxis genes., Chemotaxis allows bacteria to colonize their environment more efficiently and to find optimal growth conditions, and is consequently under strong evolutionary selection. Theoretical and experimental analyses of bacterial chemotaxis suggested that the pathway has been evolutionarily optimized to produce robust output under conditions of such physiological perturbations as stochastic intercellular variations in protein levels while at the same time minimizing complexity and cost of protein expression. Pathway topology in Escherichia coli apparently evolved to produce an invariant output under concerted variations in protein levels, consistent with experimentally observed transcriptional coupling of chemotaxis genes. Here, we show that the pathway robustness is further enhanced through the pairwise translational coupling of adjacent genes. Computer simulations predicted that the robustness of the pathway against the uncorrelated variations in protein levels can be enhanced by a selective pairwise coupling of individual chemotaxis genes on one mRNA, with the order of genes in E. coli ranking among the best in terms of noise compensation. Translational coupling between chemotaxis genes was experimentally confirmed, and coupled expression of these genes was shown to improve chemotaxis. Bioinformatics analysis further revealed that E. coli gene order corresponds to consensus in sequenced bacterial genomes, confirming evolutionary selection for noise reduction. Since polycistronic gene organization is common in bacteria, translational coupling between adjacent genes may provide a general mechanism to enhance robustness of their signaling and metabolic networks. Moreover, coupling between expression of neighboring genes is also present in eukaryotes, and similar principles of noise reduction might thus apply to all cellular networks., Author Summary All cellular networks are subject to fluctuations in the levels of their components. Robustness of the network output in the face of stochastic gene expression, or gene expression noise, is therefore essential to ensure proper function. Selection for robustness might thus have shaped much of the cellular evolution. We have used Escherichia coli chemotaxis, one of the most thoroughly studied model systems for signal transduction, to analyze the role of gene organization in robustness. Our mathematical modeling predicted that coupling the expression of chemotaxis proteins with opposing functions should buffer the output of the signaling pathway against stochastic variations in protein production. Consistent with this model, protein coexpression was indeed observed to improve chemotaxis and to be under selection during chemotaxis-driven spreading of a cell population. We show that tight coexpression is ensured by both transcriptional and translational gene coupling. We conclude that evolutionary selection for pathway robustness in the presence of gene expression noise can explain, not only the polycistronic organization of chemotaxis genes, but also the gene order within chemotaxis operons. Selection on the gene order was further confirmed by the observation of a strong bias towards specific pairwise occurrences of chemotaxis genes in sequenced prokaryotic genomes.

Simulation of cellular behavior on multiple scales requires models that are sufficiently detailed to capture central intracellular processes but at the same time enable the simulation of entire cell populations in a computationally cheap way. In this paper we present RapidCell, a hybrid model of chemotactic Escherichia coli that combines the Monod-Wyman-Changeux signal processing by mixed chemoreceptor clusters, the adaptation dynamics described by ordinary differential equations, and a detailed model of cell tumbling. Our model dramatically reduces computational costs and allows the highly efficient simulation of E. coli chemotaxis. We use the model to investigate chemotaxis in different gradients, and suggest a new, constant-activity type of gradient to systematically study chemotactic behavior of virtual bacteria. Using the unique properties of this gradient, we show that optimal chemotaxis is observed in a narrow range of CheA kinase activity, where concentration of the response regulator CheY-P falls into the operating range of flagellar motors. Our simulations also confirm that the CheB phosphorylation feedback improves chemotactic efficiency by shifting the average CheY-P concentration to fit the motor operating range. Our results suggest that in liquid media the variability in adaptation times among cells may be evolutionary favorable to ensure coexistence of subpopulations that will be optimally tactic in different gradients. However, in a porous medium (agar) such variability appears to be less important, because agar structure poses mainly negative selection against subpopulations with low levels of adaptation enzymes. RapidCell is available from the authors upon request., Author Summary Chemotaxis plays an important role in bacterial lifestyle, providing bacteria with the ability to actively search for an optimal growth environment. The chemotaxis system is likely to be highly optimized, because on the evolutionary time scale even a modest enhancement of its efficiency can give cells a large competitive advantage. In this study, we use up-to-date experimental and modeling information to construct a new computational model of chemotactic E. coli and implement it in a computationally efficient way to simulate large bacterial populations. Our simulations are performed in a new type of attractant gradient that ensures a constant level of chemotactic excitation at any position. We show that optimal chemotactic movement in a gradient results from a fine balance between excitation and adaptation. As a consequence, steeper gradients require higher optimal rates of adaptation. Simulations demonstrate that the observed intercellular variability of adaptation times, which is caused by gene expression noise, may play a positive role for the bacterial population as a whole, by allowing its subpopulations to be optimally tactic in gradients of different strengths. We further show that optimal chemotactic properties in a porous medium (agar) are different from those in a liquid.

We study wave propagation in excitable Greenberg-Hastings cellular automaton on random (Erdos-Renyi) graph, where connections are spatially constrained within radius r_c. We show that governing equation resolved by parabolic Fisher-Kolmogorov equation fails to describe wave speed, but a hyperbolic reaction-diffusion equation provides adequate wave speed for arbitrary mean network degree k. In more general case, wave speed depends rather on ratio of network moments (second moment/first moment). The wave speed is strikingly similar in different network topologies, including constant-degree, exponential and power-law degree distributions. Inspired by problem of firing propagation in networks of electrically coupled pyramidal neurons, our results may suggest spatial mean-field solution for a wide class of similar problems, such as spread of epidemics or information through real-world networks of arbitrary topology.