Substrate evaporation drives collective construction in termites

Termites build complex nests which are an impressive example of self-organization. We know that the coordinated actions involved in the construction of these nests by multiple individuals are primarily mediated by signals and cues embedded in the structure of the nest itself. However, to date there is still no scientific consensus about the nature of the stimuli that guide termite construction, and how they are sensed by termites. In order to address these questions, we studied the early building behavior of Coptotermes gestroi termites in artificial arenas, decorated with topographic cues to stimulate construction. Pellet collections were evenly distributed across the experimental setup, compatible with a collection mechanism that is not affected by local topography, but only by the distribution of termite occupancy (termites pick pellets at the positions where they are). Conversely, pellet depositions were concentrated at locations of high surface curvature and at the boundaries between different types of substrate. The single feature shared by all pellet deposition regions was that they correspond to local maxima in the evaporation flux. We can show analytically and we confirm experimentally that evaporation flux is directly proportional to the local curvature of nest surfaces. Taken together, our results indicate that surface curvature is sufficient to organize termite building activity and that termites likely sense curvature indirectly through substrate evaporation. Our findings reconcile the apparently discordant results of previous studies.

Yeast cell responses and survival during periodic osmotic stress are controlled by glucose availability

Natural environments of living organisms are often dynamic and multifactorial, with multiple parameters fluctuating over time. To better understand how cells respond to dynamically interacting factors, we quantified the effects of dual fluctuations of osmotic stress and glucose deprivation on yeast cells using microfluidics and time-lapse microscopy. Strikingly, we observed that cell proliferation, survival, and signaling depend on the phasing of the two periodic stresses. Cells divided faster, survived longer, and showed decreased transcriptional response when fluctuations of hyperosmotic stress and glucose deprivation occurred in phase than when the two stresses occurred alternatively. Therefore, glucose availability regulates yeast responses to dynamic osmotic stress, showcasing the key role of metabolic fluctuations in cellular responses to dynamic stress. We also found that mutants with impaired osmotic stress response were better adapted to alternating stresses than wild-type cells, showing that genetic mechanisms of adaptation to a persistent stress factor can be detrimental under dynamically interacting conditions.

Unifying fluidic excretion across life from cicadas to elephants

Can insects weighing mere grams challenge our current understanding of fluid dynamics in urination, jetting fluids like their larger mammalian counterparts? Current fluid urination models, predominantly formulated for mammals, suggest that jetting is confined to animals over 3 kg, owing to viscous and surface tension constraints at microscales. Our findings defy this paradigm by demonstrating that cicadas-weighing just 2 g-possess the capability for jetting fluids through remarkably small orifices. Using dimensional analysis, we introduce a unifying fluid dynamics scaling framework that accommodates a broad range of taxa, from surface-tension-dominated insects to inertia and gravity-reliant mammals. This study not only refines our understanding of fluid excretion across various species but also highlights its potential relevance in diverse fields such as ecology, evolutionary biology, and biofluid dynamics.

Caveolin-1 protects endothelial cells from extensive expansion of transcellular tunnel by stiffening the plasma membrane

Large transcellular pores elicited by bacterial mono-ADP-ribosyltransferase (mART) exotoxins inhibiting the small RhoA GTPase compromise the endothelial barrier. Recent advances in biophysical modeling point toward membrane tension and bending rigidity as the minimal set of mechanical parameters determining the nucleation and maximal size of transendothelial cell macroaperture (TEM) tunnels induced by bacterial RhoA-targeting mART exotoxins. We report that cellular depletion of caveolin-1, the membrane-embedded building block of caveolae, and depletion of cavin-1, the master regulator of caveolae invaginations, increase the number of TEMs per cell. The enhanced occurrence of TEM nucleation events correlates with a reduction in cell height due to the increase in cell spreading and decrease in cell volume, which, together with the disruption of RhoA-driven F-actin meshwork, favor membrane apposition for TEM nucleation. Strikingly, caveolin-1 specifically controls the opening speed of TEMs, leading to their dramatic 5.4-fold larger widening. Consistent with the increase in TEM density and width in siCAV1 cells, we record a higher lethality in CAV1 KO mice subjected to a catalytically active mART exotoxin targeting RhoA during staphylococcal bloodstream infection. Combined theoretical modeling with independent biophysical measurements of plasma membrane bending rigidity points toward a specific contribution of caveolin-1 to membrane stiffening in addition to the role of cavin-1/caveolin-1-dependent caveolae in the control of membrane tension homeostasis.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

C9orf72 polyPR directly binds to various nuclear transport components

The disruption of nucleocytoplasmic transport (NCT) is an important mechanism in neurodegenerative diseases. In the case of C9orf72-ALS, trafficking of macromolecules through the nuclear pore complex (NPC) might get frustrated by the binding of C9orf72-translated arginine-containing dipeptide repeat proteins (R-DPRs) to the Kapβ family of nuclear transport receptors. Besides Kapβs, several other types of transport components have been linked to NCT impairments in R-DPR-expressed cells, but the molecular origin of these observations has not been clarified. Here, we adopt a coarse-grained molecular dynamics model at amino acid resolution to study the direct interaction between polyPR, the most toxic DPR, and various nuclear transport components to elucidate the binding mechanisms and provide a complete picture of potential polyPR-mediated NCT defects. We found polyPR to directly bind to several isoforms of the Impα family, CAS (the specific exporter of Impα) and RanGAP. We observe no binding between polyPR and Ran. Longer polyPRs at lower salt concentrations also make contact with RanGEF and NTF2. Analyzing the polyPR contact sites on the transport components reveals that polyPR potentially interferes with RanGTP/RanGDP binding, with nuclear localization signal (NLS)-containing cargoes (cargo-NLS) binding to Impα, with cargo-NLS release from Impα, and with Impα export from the nucleus. The abundance of polyPR-binding sites on multiple transport components combined with the inherent polyPR length dependence makes direct polyPR interference of NCT a potential mechanistic pathway of C9orf72 toxicity.

Neural criticality from effective latent variables

Observations of power laws in neural activity data have raised the intriguing notion that brains may operate in a critical state. One example of this critical state is 'avalanche criticality', which has been observed in various systems, including cultured neurons, zebrafish, rodent cortex, and human EEG. More recently, power laws were also observed in neural populations in the mouse under an activity coarse-graining procedure, and they were explained as a consequence of the neural activity being coupled to multiple latent dynamical variables. An intriguing possibility is that avalanche criticality emerges due to a similar mechanism. Here, we determine the conditions under which latent dynamical variables give rise to avalanche criticality. We find that populations coupled to multiple latent variables produce critical behavior across a broader parameter range than those coupled to a single, quasi-static latent variable, but in both cases, avalanche criticality is observed without fine-tuning of model parameters. We identify two regimes of avalanches, both critical but differing in the amount of information carried about the latent variable. Our results suggest that avalanche criticality arises in neural systems in which activity is effectively modeled as a population driven by a few dynamical variables and these variables can be inferred from the population activity.

The positioning mechanics of microtubule asters in Drosophila embryo explants

Microtubule asters are essential in localizing the action of microtubules in processes including mitosis and organelle positioning. In large cells, such as the one-cell sea urchin embryo, aster dynamics are dominated by hydrodynamic pulling forces. However, in systems with more densely positioned nuclei such as the early Drosophila embryo, which packs around 6000 nuclei within the syncytium in a crystalline-like order, it is unclear what processes dominate aster dynamics. Here, we take advantage of a cell cycle regulation Drosophila mutant to generate embryos with multiple asters, independent from nuclei. We use an ex vivo assay to further simplify this biological system to explore the forces generated by and between asters. Through live imaging, drug and optical perturbations, and theoretical modeling, we demonstrate that these asters likely generate an effective pushing force over short distances.

Energetics of the microsporidian polar tube invasion machinery

Microsporidia are eukaryotic, obligate intracellular parasites that infect a wide range of hosts, leading to health and economic burdens worldwide. Microsporidia use an unusual invasion organelle called the polar tube (PT), which is ejected from a dormant spore at ultra-fast speeds, to infect host cells. The mechanics of PT ejection are impressive. Anncaliia algerae microsporidia spores (3-4 μm in size) shoot out a 100-nm-wide PT at a speed of 300 μm/s, creating a shear rate of 3000 s-1. The infectious cargo, which contains two nuclei, is shot through this narrow tube for a distance of ∼60-140 μm (Jaroenlak et al, 2020) and into the host cell. Considering the large hydraulic resistance in an extremely thin tube and the low-Reynolds-number nature of the process, it is not known how microsporidia can achieve this ultrafast event. In this study, we use Serial Block-Face Scanning Electron Microscopy to capture 3-dimensional snapshots of A. algerae spores in different states of the PT ejection process. Grounded in these data, we propose a theoretical framework starting with a systematic exploration of possible topological connectivity amongst organelles, and assess the energy requirements of the resulting models. We perform PT firing experiments in media of varying viscosity, and use the results to rank our proposed hypotheses based on their predicted energy requirement. We also present a possible mechanism for cargo translocation, and quantitatively compare our predictions to experimental observations. Our study provides a comprehensive biophysical analysis of the energy dissipation of microsporidian infection process and demonstrates the extreme limits of cellular hydraulics.

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population's mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn.

Diameter dependence of transport through nuclear pore complex mimics studied using optical nanopores

The nuclear pore complex (NPC) regulates the selective transport of large biomolecules through the nuclear envelope. As a model system for nuclear transport, we construct NPC mimics by functionalizing the pore walls of freestanding palladium zero-mode waveguides with the FG-nucleoporin Nsp1. This approach enables the measurement of single-molecule translocations through individual pores using optical detection. We probe the selectivity of Nsp1-coated pores by quantitatively comparing the translocation rates of the nuclear transport receptor Kap95 to the inert probe BSA over a wide range of pore sizes from 35 nm to 160 nm. Pores below 55 ± 5 nm show significant selectivity that gradually decreases for larger pores. This finding is corroborated by coarse-grained molecular dynamics simulations of the Nsp1 mesh within the pore, which suggest that leakage of BSA occurs by diffusion through transient openings within the dynamic mesh. Furthermore, we experimentally observe a modulation of the BSA permeation when varying the concentration of Kap95. The results demonstrate the potential of single-molecule fluorescence measurements on biomimetic NPCs to elucidate the principles of nuclear transport.

Energy conservation by collective movement in schooling fish

Many animals moving through fluids exhibit highly coordinated group movement that is thought to reduce the cost of locomotion. However, direct energetic measurements demonstrating the energy-saving benefits of fluid-mediated collective movements remain elusive. By characterizing both aerobic and anaerobic metabolic energy contributions in schools of giant danio (Devario aequipinnatus), we discovered that fish schools have a concave upward shaped metabolism-speed curve, with a minimum metabolic cost at ~1 body length s-1. We demonstrate that fish schools reduce total energy expenditure (TEE) per tail beat by up to 56% compared to solitary fish. When reaching their maximum sustained swimming speed, fish swimming in schools had a 44% higher maximum aerobic performance and used 65% less non-aerobic energy compared to solitary individuals, which lowered the TEE and total cost of transport by up to 53%, near the lowest recorded for any aquatic organism. Fish in schools also recovered from exercise 43% faster than solitary fish. The non-aerobic energetic savings that occur when fish in schools actively swim at high speed can considerably improve both peak and repeated performance which is likely to be beneficial for evading predators. These energetic savings may underlie the prevalence of coordinated group locomotion in fishes.

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information-theoretic framework to quantify the distribution of sensing abilities from single-cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an 'average cell'. We verify these predictions using live-cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information-theoretic framework will significantly improve our understanding of how cells sense in their environment.

Everything, everywhere, almost at once

A new platform that can follow the movement of individual proteins inside millions of cells in a single day will help contribute to existing knowledge of cell biology and identify new therapeutics.

Synchronizing beta cells in the pancreas

The secretion of insulin from the pancreas relies on both gap junctions and subpopulations of beta cells with specific intrinsic properties.

Investigating macroecological patterns in coarse-grained microbial communities using the stochastic logistic model of growth

The structure and diversity of microbial communities are intrinsically hierarchical due to the shared evolutionary history of their constituents. This history is typically captured through taxonomic assignment and phylogenetic reconstruction, sources of information that are frequently used to group microbes into higher levels of organization in experimental and natural communities. Connecting community diversity to the joint ecological dynamics of the abundances of these groups is a central problem of community ecology. However, how microbial diversity depends on the scale of observation at which groups are defined has never been systematically examined. Here, we used a macroecological approach to quantitatively characterize the structure and diversity of microbial communities among disparate environments across taxonomic and phylogenetic scales. We found that measures of biodiversity at a given scale can be consistently predicted using a minimal model of ecology, the Stochastic Logistic Model of growth (SLM). This result suggests that the SLM is a more appropriate null-model for microbial biodiversity than alternatives such as the Unified Neutral Theory of Biodiversity. Extending these within-scale results, we examined the relationship between measures of biodiversity calculated at different scales (e.g. genus vs. family), an empirical pattern previously evaluated in the context of the Diversity Begets Diversity (DBD) hypothesis (Madi et al., 2020). We found that the relationship between richness estimates at different scales can be quantitatively predicted assuming independence among community members, demonstrating that the DBD can be sufficiently explained using the SLM as a null model of ecology. Contrastingly, only by including correlations between the abundances of community members (e.g. as the consequence of interactions) can we predict the relationship between estimates of diversity at different scales. The results of this study characterize novel microbial patterns across scales of organization and establish a sharp demarcation between recently proposed macroecological patterns that are not and are affected by ecological interactions.

Recursive self-embedded vocal motifs in wild orangutans

Recursive procedures that allow placing a vocal signal inside another of a similar kind provide a neuro-computational blueprint for syntax and phonology in spoken language and human song. There are, however, no known vocal sequences among nonhuman primates arranged in self-embedded patterns that evince vocal recursion or potential incipient or evolutionary transitional forms thereof, suggesting a neuro-cognitive transformation exclusive to humans. Here, we uncover that wild flanged male orangutan long calls feature rhythmically isochronous call sequences nested within isochronous call sequences, consistent with two hierarchical strata. Remarkably, three temporally and acoustically distinct call rhythms in the lower stratum were not related to the overarching rhythm at the higher stratum by any low multiples, which suggests that these recursive structures were neither the result of parallel non-hierarchical procedures nor anatomical artifacts of bodily constraints or resonances. Findings represent a case of temporally recursive hominid vocal combinatorics in the absence of syntax, semantics, phonology, or music. Second-order combinatorics, 'sequences within sequences', involving hierarchically organized and cyclically structured vocal sounds in ancient hominids may have preluded the evolution of recursion in modern language-able humans.

Free volume theory explains the unusual behavior of viscosity in a non-confluent tissue during morphogenesis

A recent experiment on zebrafish blastoderm morphogenesis showed that the viscosity (η) of a non-confluent embryonic tissue grows sharply until a critical cell packing fraction (ϕS). The increase in η up to ϕS is similar to the behavior observed in several glass-forming materials, which suggests that the cell dynamics is sluggish or glass-like. Surprisingly, η is a constant above ϕS. To determine the mechanism of this unusual dependence of η on ϕ, we performed extensive simulations using an agent-based model of a dense non-confluent two-dimensional tissue. We show that polydispersity in the cell size, and the propensity of the cells to deform, results in the saturation of the available free area per cell beyond a critical packing fraction. Saturation in the free space not only explains the viscosity plateau above ϕS but also provides a relationship between equilibrium geometrical packing to the dramatic increase in the relaxation dynamics.

Probe-free optical chromatin deformation and measurement of differential mechanical properties in the nucleus

The nucleus is highly organized to facilitate coordinated gene transcription. Measuring the rheological properties of the nucleus and its sub-compartments will be crucial to understand the principles underlying nuclear organization. Here, we show that strongly localized temperature gradients (approaching 1°C/µm) can lead to substantial intra-nuclear chromatin displacements (>1 µm), while nuclear area and lamina shape remain unaffected. Using particle image velocimetry (PIV), intra-nuclear displacement fields can be calculated and converted into spatio-temporally resolved maps of various strain components. Using this approach, we show that chromatin displacements are highly reversible, indicating that elastic contributions are dominant in maintaining nuclear organization on the time scale of seconds. In genetically inverted nuclei, centrally compacted heterochromatin displays high resistance to deformation, giving a rigid, solid-like appearance. Correlating spatially resolved strain maps with fluorescent reporters in conventional interphase nuclei reveals that various nuclear compartments possess distinct mechanical identities. Surprisingly, both densely and loosely packed chromatin showed high resistance to deformation, compared to medium dense chromatin. Equally, nucleoli display particularly high resistance and strong local anchoring to heterochromatin. Our results establish how localized temperature gradients can be used to drive nuclear compartments out of mechanical equilibrium to obtain spatial maps of their material responses.

Hydrodynamics and multiscale order in confluent epithelia

We formulate a hydrodynamic theory of confluent epithelia: i.e. monolayers of epithelial cells adhering to each other without gaps. Taking advantage of recent progresses toward establishing a general hydrodynamic theory of p-atic liquid crystals, we demonstrate that collectively migrating epithelia feature both nematic (i.e. p = 2) and hexatic (i.e. p = 6) orders, with the former being dominant at large and the latter at small length scales. Such a remarkable multiscale liquid crystal order leaves a distinct signature in the system's structure factor, which exhibits two different power-law scaling regimes, reflecting both the hexagonal geometry of small cells clusters and the uniaxial structure of the global cellular flow. We support these analytical predictions with two different cell-resolved models of epithelia - i.e. the self-propelled Voronoi model and the multiphase field model - and highlight how momentum dissipation and noise influence the range of fluctuations at small length scales, thereby affecting the degree of cooperativity between cells. Our construction provides a theoretical framework to conceptualize the recent observation of multiscale order in layers of Madin-Darby canine kidney cells and pave the way for further theoretical developments.

Criticality supports cross-frequency cortical-thalamic information transfer during conscious states

Consciousness is thought to be regulated by bidirectional information transfer between the cortex and thalamus, but the nature of this bidirectional communication - and its possible disruption in unconsciousness - remains poorly understood. Here, we present two main findings elucidating mechanisms of corticothalamic information transfer during conscious states. First, we identify a highly preserved spectral channel of cortical-thalamic communication that is present during conscious states, but which is diminished during the loss of consciousness and enhanced during psychedelic states. Specifically, we show that in humans, mice, and rats, information sent from either the cortex or thalamus via δ/θ/α waves (∼1-13 Hz) is consistently encoded by the other brain region by high γ waves (52-104 Hz); moreover, unconsciousness induced by propofol anesthesia or generalized spike-and-wave seizures diminishes this cross-frequency communication, whereas the psychedelic 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT) enhances this low-to-high frequency interregional communication. Second, we leverage numerical simulations and neural electrophysiology recordings from the thalamus and cortex of human patients, rats, and mice to show that these changes in cross-frequency cortical-thalamic information transfer may be mediated by excursions of low-frequency thalamocortical electrodynamics toward/away from edge-of-chaos criticality, or the phase transition from stability to chaos. Overall, our findings link thalamic-cortical communication to consciousness, and further offer a novel, mathematically well-defined framework to explain the disruption to thalamic-cortical information transfer during unconscious states.

Mechanical coupling coordinates microtubule growth

During mitosis, kinetochore-attached microtubules form bundles (k-fibers) in which many filaments grow and shorten in near-perfect unison to align and segregate each chromosome. However, individual microtubules grow at intrinsically variable rates, which must be tightly regulated for a k-fiber to behave as a single unit. This exquisite coordination might be achieved biochemically, via selective binding of polymerases and depolymerases, or mechanically, because k-fiber microtubules are coupled through a shared load that influences their growth. Here, we use a novel dual laser trap assay to show that microtubule pairs growing in vitro are coordinated by mechanical coupling. Kinetic analyses show that microtubule growth is interrupted by stochastic, force-dependent pauses and indicate persistent heterogeneity in growth speed during non-pauses. A simple model incorporating both force-dependent pausing and persistent growth speed heterogeneity explains the measured coordination of microtubule pairs without any free fit parameters. Our findings illustrate how microtubule growth may be synchronized during mitosis and provide a basis for modeling k-fiber bundles with three or more microtubules, as found in many eukaryotes.

Core PCP mutations affect short-time mechanical properties but not tissue morphogenesis in the Drosophila pupal wing

How morphogenetic movements are robustly coordinated in space and time is a fundamental open question in biology. We study this question using the wing of Drosophila melanogaster, an epithelial tissue that undergoes large-scale tissue flows during pupal stages. Previously, we showed that pupal wing morphogenesis involves both cellular behaviors that allow relaxation of mechanical tissue stress, as well as cellular behaviors that appear to be actively patterned (Etournay et al., 2015). Here, we show that these active cellular behaviors are not guided by the core planar cell polarity (PCP) pathway, a conserved signaling system that guides tissue development in many other contexts. We find no significant phenotype on the cellular dynamics underlying pupal morphogenesis in mutants of core PCP. Furthermore, using laser ablation experiments, coupled with a rheological model to describe the dynamics of the response to laser ablation, we conclude that while core PCP mutations affect the fast timescale response to laser ablation they do not significantly affect overall tissue mechanics. In conclusion, our work shows that cellular dynamics and tissue shape changes during Drosophila pupal wing morphogenesis do not require core PCP as an orientational guiding cue.

Packing up the genome

Nucleotide and force-dependent mechanisms control how the viral genome of lambda bacteriophage is inserted into capsids.

An in silico testbed for fast and accurate MR labeling of orthopedic implants

One limitation on the ability to monitor health in older adults using magnetic resonance (MR) imaging is the presence of implants, where the prevalence of implantable devices (orthopedic, cardiac, neuromodulation) increases in the population, as does the pervasiveness of conditions requiring MRI studies for diagnosis (musculoskeletal diseases, infections, or cancer). The present study describes a novel multiphysics implant modeling testbed using the following approaches with two examples: (1) an in silico human model based on the widely available Visible Human Project (VHP) cryo-section dataset; (2) a finite element method (FEM) modeling software workbench from Ansys (Electronics Desktop/Mechanical) to model MR radio frequency (RF) coils and the temperature rise modeling in heterogeneous media. The in silico VHP-Female model (250 parts with an additional 40 components specifically characterizing embedded implants and resultant surrounding tissues) corresponds to a 60-year-old female with a body mass index of 36. The testbed includes the FEM-compatible in silico human model, an implant embedding procedure, a generic parameterizable MRI RF birdcage two-port coil model, a workflow for computing heat sources on the implant surface and in adjacent tissues, and a thermal FEM solver directly linked to the MR coil simulator to determine implant heating based on an MR imaging study protocol. The primary target is MR labeling of large orthopedic implants. The testbed has very recently been approved by the US Food and Drug Administration (FDA) as a medical device development tool for 1.5 T orthopedic implant examinations.

A unified approach to dissecting biphasic responses in cell signaling

Biphasic responses are encountered at all levels in biological systems. At the cellular level, biphasic dose-responses are widely encountered in cell signaling and post-translational modification systems and represent safeguards against overactivation or overexpression of species. In this paper, we provide a unified theoretical synthesis of biphasic responses in cell signaling systems, by assessing signaling systems ranging from basic biochemical building blocks to canonical network structures to well-characterized exemplars on one hand, and examining different types of doses on the other. By using analytical and computational approaches applied to a range of systems across levels (described by broadly employed models), we reveal (i) design principles enabling the presence of biphasic responses, including in almost all instances, an explicit characterization of the parameter space (ii) structural factors which preclude the possibility of biphasic responses (iii) different combinations of the presence or absence of enzyme-biphasic and substrate-biphasic responses, representing safeguards against overactivation and overexpression, respectively (iv) the possibility of broadly robust biphasic responses (v) the complete alteration of signaling behavior in a network due to biphasic interactions between species (biphasic regulation) (vi) the propensity of different co-existing biphasic responses in the Erk signaling network. These results both individually and in totality have a number of important consequences for systems and synthetic biology.

β-cell intrinsic dynamics rather than gap junction structure dictates subpopulations in the islet functional network

Diabetes is caused by the inability of electrically coupled, functionally heterogeneous β-cells within the pancreatic islet to provide adequate insulin secretion. Functional networks have been used to represent synchronized oscillatory [Ca2+] dynamics and to study β-cell subpopulations, which play an important role in driving islet function. The mechanism by which highly synchronized β-cell subpopulations drive islet function is unclear. We used experimental and computational techniques to investigate the relationship between functional networks, structural (gap junction) networks, and intrinsic β-cell dynamics in slow and fast oscillating islets. Highly synchronized subpopulations in the functional network were differentiated by intrinsic dynamics, including metabolic activity and KATP channel conductance, more than structural coupling. Consistent with this, intrinsic dynamics were more predictive of high synchronization in the islet functional network as compared to high levels of structural coupling. Finally, dysfunction of gap junctions, which can occur in diabetes, caused decreases in the efficiency and clustering of the functional network. These results indicate that intrinsic dynamics rather than structure drive connections in the functional network and highly synchronized subpopulations, but gap junctions are still essential for overall network efficiency. These findings deepen our interpretation of functional networks and the formation of functional subpopulations in dynamic tissues such as the islet.

Mechanical vibration patterns elicit behavioral transitions and habituation in crawling Drosophila larvae

How animals respond to repeatedly applied stimuli, and how animals respond to mechanical stimuli in particular, are important questions in behavioral neuroscience. We study adaptation to repeated mechanical agitation using the Drosophila larva. Vertical vibration stimuli elicit a discrete set of responses in crawling larvae: continuation, pause, turn, and reversal. Through high-throughput larva tracking, we characterize how the likelihood of each response depends on vibration intensity and on the timing of repeated vibration pulses. By examining transitions between behavioral states at the population and individual levels, we investigate how the animals habituate to the stimulus patterns. We identify time constants associated with desensitization to prolonged vibration, with re-sensitization during removal of a stimulus, and additional layers of habituation that operate in the overall response. Known memory-deficient mutants exhibit distinct behavior profiles and habituation time constants. An analogous simple electrical circuit suggests possible neural and molecular processes behind adaptive behavior.

The Reissner fiber under tension in vivo shows dynamic interaction with ciliated cells contacting the cerebrospinal fluid

The Reissner fiber (RF) is an acellular thread positioned in the midline of the central canal that aggregates thanks to the beating of numerous cilia from ependymal radial glial cells (ERGs) generating flow in the central canal of the spinal cord. RF together with cerebrospinal fluid (CSF)-contacting neurons (CSF-cNs) form an axial sensory system detecting curvature. How RF, CSF-cNs and the multitude of motile cilia from ERGs interact in vivo appears critical for maintenance of RF and sensory functions of CSF-cNs to keep a straight body axis, but is not well-understood. Using in vivo imaging in larval zebrafish, we show that RF is under tension and resonates dorsoventrally. Focal RF ablations trigger retraction and relaxation of the fiber's cut ends, with larger retraction speeds for rostral ablations. We built a mechanical model that estimates RF stress diffusion coefficient D at 5 mm2/s and reveals that tension builds up rostrally along the fiber. After RF ablation, spontaneous CSF-cN activity decreased and ciliary motility changed, suggesting physical interactions between RF and cilia projecting into the central canal. We observed that motile cilia were caudally-tilted and frequently interacted with RF. We propose that the numerous ependymal motile monocilia contribute to RF's heterogenous tension via weak interactions. Our work demonstrates that under tension, the Reissner fiber dynamically interacts with motile cilia generating CSF flow and spinal sensory neurons.

A mechanosensing mechanism controls plasma membrane shape homeostasis at the nanoscale

As cells migrate and experience forces from their surroundings, they constantly undergo mechanical deformations which reshape their plasma membrane (PM). To maintain homeostasis, cells need to detect and restore such changes, not only in terms of overall PM area and tension as previously described, but also in terms of local, nanoscale topography. Here, we describe a novel phenomenon, by which cells sense and restore mechanically induced PM nanoscale deformations. We show that cell stretch and subsequent compression reshape the PM in a way that generates local membrane evaginations in the 100 nm scale. These evaginations are recognized by I-BAR proteins, which triggers a burst of actin polymerization mediated by Rac1 and Arp2/3. The actin polymerization burst subsequently re-flattens the evagination, completing the mechanochemical feedback loop. Our results demonstrate a new mechanosensing mechanism for PM shape homeostasis, with potential applicability in different physiological scenarios.

Force propagation between epithelial cells depends on active coupling and mechano-structural polarization

Cell-generated forces play a major role in coordinating the large-scale behavior of cell assemblies, in particular during development, wound healing, and cancer. Mechanical signals propagate faster than biochemical signals, but can have similar effects, especially in epithelial tissues with strong cell-cell adhesion. However, a quantitative description of the transmission chain from force generation in a sender cell, force propagation across cell-cell boundaries, and the concomitant response of receiver cells is missing. For a quantitative analysis of this important situation, here we propose a minimal model system of two epithelial cells on an H-pattern ('cell doublet'). After optogenetically activating RhoA, a major regulator of cell contractility, in the sender cell, we measure the mechanical response of the receiver cell by traction force and monolayer stress microscopies. In general, we find that the receiver cells show an active response so that the cell doublet forms a coherent unit. However, force propagation and response of the receiver cell also strongly depend on the mechano-structural polarization in the cell assembly, which is controlled by cell-matrix adhesion to the adhesive micropattern. We find that the response of the receiver cell is stronger when the mechano-structural polarization axis is oriented perpendicular to the direction of force propagation, reminiscent of the Poisson effect in passive materials. We finally show that the same effects are at work in small tissues. Our work demonstrates that cellular organization and active mechanical response of a tissue are key to maintain signal strength and lead to the emergence of elasticity, which means that signals are not dissipated like in a viscous system, but can propagate over large distances.

Continuous, long-term crawling behavior characterized by a robotic transport system

Detailed descriptions of behavior provide critical insight into the structure and function of nervous systems. In Drosophila larvae and many other systems, short behavioral experiments have been successful in characterizing rapid responses to a range of stimuli at the population level. However, the lack of long-term continuous observation makes it difficult to dissect comprehensive behavioral dynamics of individual animals and how behavior (and therefore the nervous system) develops over time. To allow for long-term continuous observations in individual fly larvae, we have engineered a robotic instrument that automatically tracks and transports larvae throughout an arena. The flexibility and reliability of its design enables controlled stimulus delivery and continuous measurement over developmental time scales, yielding an unprecedented level of detailed locomotion data. We utilize the new system's capabilities to perform continuous observation of exploratory search behavior over a duration of 6 hr with and without a thermal gradient present, and in a single larva for over 30 hr. Long-term free-roaming behavior and analogous short-term experiments show similar dynamics that take place at the beginning of each experiment. Finally, characterization of larval thermotaxis in individuals reveals a bimodal distribution in navigation efficiency, identifying distinct phenotypes that are obfuscated when only analyzing population averages.

Continuous odor profile monitoring to study olfactory navigation in small animals

Olfactory navigation is observed across species and plays a crucial role in locating resources for survival. In the laboratory, understanding the behavioral strategies and neural circuits underlying odor-taxis requires a detailed understanding of the animal's sensory environment. For small model organisms like Caenorhabditis elegans and larval Drosophila melanogaster, controlling and measuring the odor environment experienced by the animal can be challenging, especially for airborne odors, which are subject to subtle effects from airflow, temperature variation, and from the odor's adhesion, adsorption, or reemission. Here, we present a method to control and measure airborne odor concentration in an arena compatible with an agar substrate. Our method allows continuous controlling and monitoring of the odor profile while imaging animal behavior. We construct stationary chemical landscapes in an odor flow chamber through spatially patterned odorized air. The odor concentration is measured with a spatially distributed array of digital gas sensors. Careful placement of the sensors allows the odor concentration across the arena to be continuously inferred in space and monitored through time. We use this approach to measure the odor concentration that each animal experiences as it undergoes chemotaxis behavior and report chemotaxis strategies for C. elegans and D. melanogaster larvae populations as they navigate spatial odor landscapes.

Multiplexed microfluidic screening of bacterial chemotaxis

Microorganism sensing of and responding to ambient chemical gradients regulates a myriad of microbial processes that are fundamental to ecosystem function and human health and disease. The development of efficient, high-throughput screening tools for microbial chemotaxis is essential to disentangling the roles of diverse chemical compounds and concentrations that control cell nutrient uptake, chemorepulsion from toxins, and microbial pathogenesis. Here, we present a novel microfluidic multiplexed chemotaxis device (MCD) which uses serial dilution to simultaneously perform six parallel bacterial chemotaxis assays that span five orders of magnitude in chemostimulant concentration on a single chip. We first validated the dilution and gradient generation performance of the MCD, and then compared the measured chemotactic response of an established bacterial chemotaxis system (Vibrio alginolyticus) to a standard microfluidic assay. Next, the MCD's versatility was assessed by quantifying the chemotactic responses of different bacteria (Psuedoalteromonas haloplanktis, Escherichia coli) to different chemoattractants and chemorepellents. The MCD vastly accelerates the chemotactic screening process, which is critical to deciphering the complex sea of chemical stimuli underlying microbial responses.

Fluid extraction from the left-right organizer uncovers mechanical properties needed for symmetry breaking

Humans and other vertebrates define body axis left-right asymmetry in the early stages of embryo development. The mechanism behind left-right establishment is not fully understood. Symmetry breaking occurs in a dedicated organ called the left-right organizer (LRO) and involves motile cilia generating fluid-flow therein. However, it has been a matter of debate whether the process of symmetry breaking relies on a chemosensory or a mechanosensory mechanism (Shinohara et al., 2012). Novel tailored manipulations for LRO fluid extraction in living zebrafish embryos allowed us to pinpoint a physiological developmental period for breaking left-right symmetry during development. The shortest critical time-window was narrowed to one hour and characterized by a mild counterclockwise flow. The experimental challenge consisted in emptying the LRO of its fluid, abrogating simultaneously flow force and chemical determinants. Our findings revealed an unprecedented recovery capacity of the embryo to re-fil and re-circulate new LRO fluid. The embryos that later developed laterality problems were found to be those that had lower anterior angular velocity and thus less anterior-posterior heterogeneity. Next, aiming to test the presence of any secreted determinant, we replaced the extracted LRO fluid by a physiological buffer. Despite some transitory flow homogenization, laterality defects were absent unless viscosity was altered, demonstrating that symmetry breaking does not depend on the nature of the fluid content but is rather sensitive to fluid mechanics. Altogether, we conclude that the zebrafish LRO is more sensitive to fluid dynamics for symmetry breaking.

Protein compactness and interaction valency define the architecture of a biomolecular condensate across scales

Non-membrane-bound biomolecular condensates have been proposed to represent an important mode of subcellular organization in diverse biological settings. However, the fundamental principles governing the spatial organization and dynamics of condensates at the atomistic level remain unclear. The Saccharomyces cerevisiae Lge1 protein is required for histone H2B ubiquitination and its N-terminal intrinsically disordered fragment (Lge11-80) undergoes robust phase separation. This study connects single- and multi-chain all-atom molecular dynamics simulations of Lge11-80 with the in vitro behavior of Lge11-80 condensates. Analysis of modeled protein-protein interactions elucidates the key determinants of Lge11-80 condensate formation and links configurational entropy, valency, and compactness of proteins inside the condensates. A newly derived analytical formalism, related to colloid fractal cluster formation, describes condensate architecture across length scales as a function of protein valency and compactness. In particular, the formalism provides an atomistically resolved model of Lge11-80 condensates on the scale of hundreds of nanometers starting from individual protein conformers captured in simulations. The simulation-derived fractal dimensions of condensates of Lge11-80 and its mutants agree with their in vitro morphologies. The presented framework enables a multiscale description of biomolecular condensates and embeds their study in a wider context of colloid self-organization.

Regulation of chromatin microphase separation by binding of protein complexes

We show evidence of the association of RNA polymerase II (RNAP) with chromatin in a core-shell organization, reminiscent of microphase separation where the cores comprise dense chromatin and the shell, RNAP and chromatin with low density. These observations motivate our physical model for the regulation of core-shell chromatin organization. Here, we model chromatin as a multiblock copolymer, comprising active and inactive regions (blocks) that are both in poor solvent and tend to be condensed in the absence of binding proteins. However, we show that the solvent quality for the active regions of chromatin can be regulated by the binding of protein complexes (e.g., RNAP and transcription factors). Using the theory of polymer brushes, we find that such binding leads to swelling of the active chromatin regions which in turn modifies the spatial organization of the inactive regions. In addition, we use simulations to study spherical chromatin micelles, whose cores comprise inactive regions and shells comprise active regions and bound protein complexes. In spherical micelles the swelling increases the number of inactive cores and controls their size. Thus, genetic modifications affecting the binding strength of chromatin-binding protein complexes may modulate the solvent quality experienced by chromatin and regulate the physical organization of the genome.

Kinetochore-fiber lengths are maintained locally but coordinated globally by polesin the mammalian spindle

At each cell division, nanometer-scale components self-organize to build a micron-scale spindle. In mammalian spindles, microtubule bundles called kinetochore-fibers attach to chromosomes and focus into spindle poles. Despite evidence suggesting that poles can set spindle length, their role remains poorly understood. In fact, many species do not have spindle poles. Here, we probe the pole's contribution to mammalian spindle length, dynamics, and function by inhibiting dynein to generate spindles whose kinetochore-fibers do not focus into poles, yet maintain a metaphase steady-state length. We find that unfocused kinetochore-fibers have a mean length indistinguishable from control, but a broader length distribution, and reduced length coordination between sisters and neighbors. Further, we show that unfocused kinetochore-fibers, like control, can grow back to their steady-state length if acutely shortened by drug treatment or laser ablation: they recover their length by tuning their end dynamics, albeit slower due to their reduced baseline dynamics. Thus, kinetochore-fiber dynamics are regulated by their length, not just pole-focusing forces. Finally, we show that spindles with unfocused kinetochore-fibers can segregate chromosomes but fail to correctly do so. We propose that mammalian spindle length emerges locally from individual k-fibers while spindle poles globally coordinate k-fibers across space and time.

Spontaneous body wall contractions stabilize the fluid microenvironment that shapes host-microbe associations

The freshwater polyp Hydra is a popular biological model system; however, we still do not understand one of its most salient behaviours, the generation of spontaneous body wall contractions. Here, by applying experimental fluid dynamics analysis and mathematical modelling, we provide functional evidence that spontaneous contractions of body walls enhance the transport of chemical compounds from and to the tissue surface where symbiotic bacteria reside. Experimentally, a reduction in the frequency of spontaneous body wall contractions is associated with a changed composition of the colonizing microbiota. Together, our findings suggest that spontaneous body wall contractions create an important fluid transport mechanism that (1) may shape and stabilize specific host-microbe associations and (2) create fluid microhabitats that may modulate the spatial distribution of the colonizing microbes. This mechanism may be more broadly applicable to animal-microbe interactions since research has shown that rhythmic spontaneous contractions in the gastrointestinal tracts are essential for maintaining normal microbiota.

Putting scaling laws on a physical foundation

As a cell changes size during the cell cycle, why does its density remain constant?

Exploring the dynamics of vascular adaptation

A combination of in toto imaging and theory suggests a new mechanism for the remodeling of veins in vascular networks.

Substrate stiffness impacts early biofilm formation by modulating Pseudomonas aeruginosa twitching motility

Surface-associated lifestyles dominate in the bacterial world. Large multicellular assemblies, called biofilms, are essential to the survival of bacteria in harsh environments, and are closely linked to antibiotic resistance in pathogenic strains. Biofilms stem from the surface colonization of a wide variety of substrates encountered by bacteria, from living tissues to inert materials. Here, we demonstrate experimentally that the promiscuous opportunistic pathogen Pseudomonas aeruginosa explores substrates differently based on their rigidity, leading to striking variations in biofilm structure, exopolysaccharides (EPS) distribution, strain mixing during co-colonization and phenotypic expression. Using simple kinetic models, we show that these phenotypes arise through a mechanical interaction between the elasticity of the substrate and the type IV pilus (T4P) machinery, that mediates the surface-based motility called twitching. Together, our findings reveal a new role for substrate softness in the spatial organization of bacteria in complex microenvironments, with far-reaching consequences on efficient biofilm formation.

Physical basis of the cell size scaling laws

Cellular growth is the result of passive physical constraints and active biological processes. Their interplay leads to the appearance of robust and ubiquitous scaling laws relating linearly cell size, dry mass, and nuclear size. Despite accumulating experimental evidence, their origin is still unclear. Here, we show that these laws can be explained quantitatively by a single model of size regulation based on three simple, yet generic, physical constraints defining altogether the Pump-Leak model. Based on quantitative estimates, we clearly map the Pump-Leak model coarse-grained parameters with the dominant cellular components. We propose that dry mass density homeostasis arises from the scaling between proteins and small osmolytes, mainly amino-acids and ions. Our model predicts this scaling to naturally fail, both at senescence when DNA and RNAs are saturated by RNA polymerases and ribosomes respectively, and at mitotic entry due to the counterion release following histone tail modifications. Based on the same physical laws, we further show that nuclear scaling results from a osmotic balance at the nuclear envelope and a large pool of metabolites, which dilutes chromatin counterions that do not scale during growth.

Making it hard to replicate

Understanding how to harden liquid condensates produced by influenza A virus could accelerate the development of novel antiviral drugs.

Comparable in vivo joint kinematics between self-reported stable and unstable knees after TKA can be explained by muscular adaptation strategies: a retrospective observational study

Background: Postoperative knee instability is one of the major reasons accounting for unsatisfactory outcomes, as well as a major failure mechanism leading to total knee arthroplasty (TKA) revision. Nevertheless, subjective knee instability is not well defined clinically, plausibly because the relationships between instability and implant kinematics during functional activities of daily living remain unclear. Although muscles play a critical role in supporting the dynamic stability of the knee joint, the influence of joint instability on muscle synergy patterns is poorly understood. Therefore, this study aimed to understand the impact of self-reported joint instability on tibiofemoral kinematics and muscle synergy patterns after TKA during functional gait activities of daily living.

Methods: Tibiofemoral kinematics and muscle synergy patterns were examined during level walking, downhill walking, and stair descent in eight self-reported unstable knees after TKA (3M:5F, 68.9±8.3 years, BMI 26.1±3.2 kg/m², 31.9±20.4 months postoperatively), and compared against ten stable TKA knees (7M:3F, 62.6±6.8 years, 33.9±8.5 months postoperatively, BMI 29.4±4.8 kg/m²). For each knee joint, clinical assessments of postoperative outcome were performed, while joint kinematics were evaluated using moving video-fluoroscopy, and muscle synergy patterns were recorded using electromyography.

Results: Our results reveal that average condylar A-P translations, rotations, as well as their ranges of motion were comparable between stable and unstable groups. However, the unstable group exhibited more heterogeneous muscle synergy patterns and prolonged activation of knee flexors compared to the stable group. In addition, subjects who reported instability events during measurement showed distinct, subject-specific tibiofemoral kinematic patterns in the early/mid-swing phase of gait.

Conclusions: Our findings suggest that accurate movement analysis is sensitive for detecting acute instability events, but might be less robust in identifying general joint instability. Conversely, muscle synergy patterns seem to be able to identify muscular adaptation associated with underlying chronic knee instability.

Funding: This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Mechanical basis and topological routes to cell elimination

Cell layers eliminate unwanted cells through the extrusion process, which underlines healthy versus flawed tissue behaviors. Although several biochemical pathways have been identified, the underlying mechanical basis including the forces involved in cellular extrusion remains largely unexplored. Utilizing a phase-field model of a three-dimensional cell layer, we study the interplay of cell extrusion with cell-cell and cell-substrate interactions in a flat monolayer. Independent tuning of cell-cell versus cell-substrate adhesion forces reveals that extrusion events can be distinctly linked to defects in nematic and hexatic orders associated with cellular arrangements. Specifically, we show that by increasing relative cell-cell adhesion forces the cell monolayer can switch between the collective tendency towards fivefold, hexatic, disclinations relative to half-integer, nematic, defects for extruding a cell. We unify our findings by accessing three-dimensional mechanical stress fields to show that an extrusion event acts as a mechanism to relieve localized stress concentration.

Generating active T1 transitions through mechanochemical feedback

Convergence-extension in embryos is controlled by chemical and mechanical signalling. A key cellular process is the exchange of neighbours via T1 transitions. We propose and analyse a model with positive feedback between recruitment of myosin motors and mechanical tension in cell junctions. The model produces active T1 events, which act to elongate the tissue perpendicular to the main direction of tissue stress. Using an idealized tissue patch comprising several active cells embedded in a matrix of passive hexagonal cells we identified an optimal range of mechanical stresses to trigger an active T1 event. We show that directed stresses also generate tension chains in a realistic patch made entirely of active cells of random shapes, and leads to convergence-extension over a range of parameters. Our findings show that active intercalations can generate stress that activates T1 events in neighbouring cells resulting in tension dependent tissue reorganisation, in qualitative agreement with experiments on gastrulation in chick embryos.

Defining basic rules for hardening influenza A virus liquid condensates

In biological systems, liquid and solid-like biomolecular condensates may contain the same molecules but their behaviour, including movement, elasticity and viscosity, is different on account of distinct physicochemical properties. As such, it is known that phase transitions affect the function of biological condensates and that material properties can be tuned by several factors including temperature, concentration and valency. It is, however, unclear if some factors are more efficient than others at regulating their behaviour. Viral infections are good systems to address this question as they form condensates de novo as part of their replication programmes. Here, we used influenza A virus liquid cytosolic condensates, A.K.A viral inclusions, to provide a proof of concept that liquid condensate hardening via changes in the valency of its components is more efficient than altering their concentration or the temperature of the cell. Liquid IAV inclusions may be hardened by targeting vRNP interactions via the known NP oligomerizing molecule, nucleozin, both in vitro and in vivo without affecting host proteome abundance nor solubility. This study is a starting point for understanding how to pharmacologically modulate the material properties of IAV inclusions and may offer opportunities for alternative antiviral strategies.

The benefits of swimming together

When a fish beats its tail, it produces vortices in the water that other fish could take advantage of to save energy while swimming.