The influence of entropic crowding in cell monolayers

Cell-cell interaction dictates cell morphology and organization, which play a crucial role in the micro-architecture of tissues that guides their biological and mechanical functioning. Here, we investigate the effect of cell density on the responses of cells seeded on flat substrates using a novel statistical thermodynamics framework. The framework recognizes the existence of nonthermal fluctuations in cellular response and thereby naturally captures entropic interactions between cells in monolayers. In line with observations, the model predicts that cell area and elongation decrease with increasing cell seeding density—both are a direct outcome of the fluctuating nature of the cellular response that gives rise to enhanced cell-cell interactions with increasing cell crowding. The modeling framework also predicts the increase in cell alignment with increasing cell density: this cellular ordering is also due to enhanced entropic interactions and is akin to nematic ordering in liquid crystals. Our simulations provide physical insights that suggest that entropic cell-cell interactions play a crucial role in governing the responses of cell monolayers.     

Cholesterol Depletion by MβCD Enhances Cell Membrane Tension and Its Variations-Reducing Integrity

Cholesterol depletion by methyl-β-cyclodextrin (MβCD) remodels the plasma membrane’s mechanics in cells and its interactions with the underlying cytoskeleton, whereas in red blood cells, it is also known to cause lysis. Currently it’s unclear if MβCD alters membrane tension or only enhances membrane-cytoskeleton interactions—and how this relates to cell lysis. We map membrane height fluctuations in single cells and observe that MβCD reduces temporal fluctuations robustly but flattens spatial membrane undulations only slightly. Utilizing models explicitly incorporating membrane confinement besides other viscoelastic factors, we estimate membrane mechanical parameters from the fluctuations’ frequency spectrum. This helps us conclude that MβCD enhances membrane tension and does so even on ATP-depleted cell membranes where this occurs despite reduction in confinement. Additionally, on cholesterol depletion, cell membranes display higher intracellular heterogeneity in the amplitude of spatial undulations and membrane tension. MβCD also has a strong impact on the cell membrane’s tenacity to mechanical stress, making cells strongly prone to rupture on hypo-osmotic shock with larger rupture diameters—an effect not hindered by actomyosin perturbations. Our study thus demonstrates that cholesterol depletion increases membrane tension and its variability, making cells prone to rupture independent of the cytoskeletal state of the cell.     

ECM Cross-Linking Regulates Invadopodia Dynamics

Invadopodia are membrane protrusions dynamically assembled by invasive cancer cells in contact with the extracellular matrix (ECM). Invadopodia are enriched by the structural proteins actin and cortactin as well as metalloproteases such as MT1-MMP, whose function is to degrade the surrounding ECM. During metastasis, invadopodia are necessary for cancer cell intravasation and extravasation. Although signaling pathways involved in the assembly and function of invadopodia are well studied, few studies address invadopodia dynamics and how the cell-ECM interactions contribute to cell invasion. Using iterative analysis based on time-lapse microscopy and mathematical modeling of invasive cancer cells, we found that cells oscillate between invadopodia presence and cell stasis—termed the “invadopodia state”—and invadopodia absence during cell translocation—termed the “migration state.” Our data suggest that β1-integrin-ECM binding and ECM cross-linking control the duration of each of the two states. By changing the concentration of cross-linkers in two-dimensional and three-dimensional cultures, we generate an ECM in which 0–0.92 of total lysine residues are cross-linked. Using an ECM with a range of cross-linking degrees, we demonstrate that the dynamics of invadopodia-related functions have a biphasic relationship to ECM cross-linking. At intermediate levels of ECM cross-linking (0.39), cells exhibit rapid invadopodia protrusion-retraction cycles and rapid calcium spikes, which lead to more frequent MT1-MMP delivery, causing maximal invadopodia-mediated ECM degradation. In contrast, both extremely high or low levels of cross-linking lead to slower invadopodia-related dynamics and lower ECM degradation. Additionally, β1-integrin inhibition modifies the dynamics of invadopodia-related functions as well as the length of time cells spend in either of the states. Collectively, these data suggest that β1-integrin-ECM binding nonlinearly translates small physical differences in the extracellular environment to differences in the dynamics of cancer cell behaviors. Understanding the conditions under which invadopodia can be reduced by subtle environment-targeting treatments may lead to combination therapies for preventing metastatic spread.     

Hippo pathway regulation by cell morphology and stress fibers

The Hippo signaling pathway plays an important role in regulation of cell proliferation. Cell density regulates the Hippo pathway in cultured cells; however, the mechanism by which cells detect density remains unclear. In this study, we demonstrated that changes in cell morphology are a key factor. Morphological manipulation of single cells without cell-cell contact resulted in flat spread or round compact cells with nuclear or cytoplasmic Yap, respectively. Stress fibers increased in response to expanded cell areas, and F-actin regulated Yap downstream of cell morphology. Cell morphology- and F-actin-regulated phosphorylation of Yap, and the effects of F-actin were suppressed by modulation of Lats. Our results suggest that cell morphology is an important factor in the regulation of the Hippo pathway, which is mediated by stress fibers consisting of F-actin acting upstream of, or on Lats, and that cells can detect density through their resulting morphology. This cell morphology (stress-fiber)-mediated mechanism probably cooperates with a cell-cell contact (adhesion)-mediated mechanism involving the Hippo pathway to achieve density-dependent control of cell proliferation.     

A thermodynamically motivated model for stress-fiber reorganization

We present a model for stress-fiber reorganization and the associated contractility that includes both the kinetics of stress-fiber formation and dissociation as well as the kinetics of stress-fiber remodeling. These kinetics are motivated by considering the enthalpies of the actin/myosin functional units that constitute the stress fibers. The stress, strain and strain rate dependence of the stress-fiber dynamics are natural outcomes of the approach. The model is presented in a general 3D framework and includes the transport of the unbound stress-fiber proteins. Predictions of the model for a range of cyclic loadings are illustrated to rationalize hitherto apparently contrasting observations. These observations include: (1) For strain amplitudes around 10 % and cyclic frequencies of about 1 Hz, stress fibers align perpendicular to the straining direction in cells subjected to cyclic straining on a 2D substrate while the stress fibers align parallel with the straining direction in cells constrained in a 3D tissue. (2) At lower applied cyclic frequencies, stress fibers in cells on 2D substrates display no sensitivity to symmetric applied strain versus time waveforms but realign in response to applied loadings with a fast lengthening rate and slow shortening. (3) At very low applied cyclic frequencies (on the order of mHz) with symmetric strain versus time waveforms, cells on 2D substrates orient perpendicular to the direction of cyclic straining above a critical strain amplitude.     

Distribution and Dynamics of Rat Basophilic Leukemia Immunoglobulin E Receptors (Fc?RI) on Planar Ligand-Presenting Surfaces

There is considerable interest in the signaling mechanisms of immunoreceptors, especially when triggered with membrane-bound ligands. We have quantified the spatiotemporal dynamics of the redistribution of immunoglobulin E-loaded receptors (IgE-FcɛRI) on rat basophilic leukemia-2H3 mast cells in contact with fluid and gel-phase membranes displaying ligands for immunoglobulin E, using total internal reflection fluorescence microscopy. To clearly separate the kinetics of receptor redistribution from cell spreading, and to precisely define the initial contact time (±50 ms), micropipette cell manipulation was used to bring individual cells into contact with surfaces. On ligand-free surfaces, there are micron-scale heterogeneities in fluorescence that likely reflect regions of the cell that are more closely apposed to the substrate. When ligands are present, receptor clusters form with this same size scale. The initial rate of accumulation of receptors into the clusters is consistent with diffusion-limited trapping with D ∼10⁻¹μm²/s. These results support the hypothesis that clusters form by diffusion to cell-surface contact regions. Over longer timescales (>10 s), individual clusters moved with both diffusive and directed motion components. The dynamics of the cluster motion is similar to the dynamics of membrane fluctuations of cells on ligand-free fluid membranes. Thus, the same cellular machinery may be responsible for both processes.     

Analysis of cell locomotion

The methods of statistical physics have been applied to the analysis of cell movement. Human polymorphonuclear leukocytes were exposed to different surfaces possessing parallel oriented physical stuctures (scratched glass surface, machine drilled aluminium surface, optical grid and stretched polyethylene foil) and cell migration was observed using time-lapse photography.We demonstrate that in cell migration along physical structures, referred to as contact guidance, two subgroups can be distinguished: 1) The nematic type where the cell size is large in relation to the grid distance of the undulate surface. 2) The smectic type where the cell size is small in relation to the grid distance of the substrate.Nematic contact guidance is characterized by an anisotropic random walk. In all substrates investigated the diffusion process parallel to the lines was faster than the diffusion process perpendicular to them. The angular dependent diffusion coefficient was described by an ellipse. Deviation from a circle defined an apolar order parameter, whose value was about 0.3. The amount of information which the cells collected from, the undulate surface was very low, between 0.1 and 0.2 bits. We demonstrate that cells do not recognize all the details of their surroundings and that their migration can be compared to the groping around of a short sighted man. The blurred environment can be described by a mean field whose strength is proportional to the aploar order parameter. It is argued that the anisotropic surface tension is the basic source for nematic contact guidance.Smectic contact guidance is characterized by an anisotropic random walk and is quantified by a density order paramter which is 0.28 in the case of the scratched glass surface of a Neubauer counting chamber. The information which the cells collect from their environment is very low (0.03 bits). The lines seen by the cell can be described by a mean field whose strength is proportional to the density oder parameter.Finally, we demonstrate that the locomotion of granulocytes is governed by an internal clock and internal programs. After migrating for a certain time (32 s) in a particular direction, a new direction of locomotion is determined by an internal program. The cell decides basically between left or right, thereby preferring a turn angle such that the cell migrates either parallel or perpendicular to the lines. the angles are nearly equally probable but the cell moves, in the case of nematic guidance, with different velocities in the +or-direction. The cell also has directional memories with characteristic times of 32 s and greater than 100 s.     

How Entorhinal Grid Cells May Learn Multiple Spatial Scales from a Dorsoventral Gradient of Cell Response Rates in a Self-organizing Map

Place cells in the hippocampus of higher mammals are critical for spatial navigation. Recent modeling clarifies how this may be achieved by how grid cells in the medial entorhinal cortex (MEC) input to place cells. Grid cells exhibit hexagonal grid firing patterns across space in multiple spatial scales along the MEC dorsoventral axis. Signals from grid cells of multiple scales combine adaptively to activate place cells that represent much larger spaces than grid cells. But how do grid cells learn to fire at multiple positions that form a hexagonal grid, and with spatial scales that increase along the dorsoventral axis? In vitro recordings of medial entorhinal layer II stellate cells have revealed subthreshold membrane potential oscillations (MPOs) whose temporal periods, and time constants of excitatory postsynaptic potentials (EPSPs), both increase along this axis. Slower (faster) subthreshold MPOs and slower (faster) EPSPs correlate with larger (smaller) grid spacings and field widths. A self-organizing map neural model explains how the anatomical gradient of grid spatial scales can be learned by cells that respond more slowly along the gradient to their inputs from stripe cells of multiple scales, which perform linear velocity path integration. The model cells also exhibit MPO frequencies that covary with their response rates. The gradient in intrinsic rhythmicity is thus not compelling evidence for oscillatory interference as a mechanism of grid cell firing. A response rate gradient combined with input stripe cells that have normalized receptive fields can reproduce all known spatial and temporal properties of grid cells along the MEC dorsoventral axis. This spatial gradient mechanism is homologous to a gradient mechanism for temporal learning in the lateral entorhinal cortex and its hippocampal projections. Spatial and temporal representations may hereby arise from homologous mechanisms, thereby embodying a mechanistic “neural relativity” that may clarify how episodic memories are learned.     

Data-driven and model-guided systematic framework for media development in CHO cell culture

Proposed herein is a systematic media design framework that combines multivariate statistical approaches with in silico analysis of a genome-scale metabolic model of Chinese hamster ovary cell. The framework comprises sequential modules including cell culture and metabolite data collection, multivariate data analysis, in silico modeling and flux prediction, and knowledge-based identification of target media components. Two monoclonal antibody-producing cell lines under two different media conditions were used to demonstrate the applicability of the framework. First, the cell culture and metabolite profiles from all conditions were generated, and then statistically and mechanistically analyzed to explore combinatorial effects of cell line and media on intracellular metabolism. As a result, we found a metabolic bottleneck via a redox imbalance in the TCA cycle in the poorest growth condition, plausibly due to inefficient coenzyme q10-q10h2 recycling. Subsequent in silico simulation allowed us to suggest q10 supplementation to debottleneck the imbalance for the enhanced cellular energy state and TCA cycle activity. Finally, experimental validation was successfully conducted by adding q10 in the media, resulting in increased cell growth. Taken together, the proposed framework rationally identified target nutrients for cell line-specific media design and reformulation, which could greatly improve cell culture performance.     

Non-linear Membrane Properties in Entorhinal Cortical Stellate Cells Reduce Modulation of Input-Output Responses by Voltage Fluctuations

The presence of voltage fluctuations arising from synaptic activity is a critical component in models of gain control, neuronal output gating, and spike rate coding. The degree to which individual neuronal input-output functions are modulated by voltage fluctuations, however, is not well established across different cortical areas. Additionally, the extent and mechanisms of input-output modulation through fluctuations have been explored largely in simplified models of spike generation, and with limited consideration for the role of non-linear and voltage-dependent membrane properties. To address these issues, we studied fluctuation-based modulation of input-output responses in medial entorhinal cortical (MEC) stellate cells of rats, which express strong sub-threshold non-linear membrane properties. Using in vitro recordings, dynamic clamp and modeling, we show that the modulation of input-output responses by random voltage fluctuations in stellate cells is significantly limited. In stellate cells, a voltage-dependent increase in membrane resistance at sub-threshold voltages mediated by Na⁺ conductance activation limits the ability of fluctuations to elicit spikes. Similarly, in exponential leaky integrate-and-fire models using a shallow voltage-dependence for the exponential term that matches stellate cell membrane properties, a low degree of fluctuation-based modulation of input-output responses can be attained. These results demonstrate that fluctuation-based modulation of input-output responses is not a universal feature of neurons and can be significantly limited by subthreshold voltage-gated conductances.     

Size and stripes: how fish clients recognize cleaners

Little is known of how individuals find each other in interspecific mutualisms involving free-living partners. We tested the importance of two factors, namely body size and the presence of a lateral body stripe, in the recognition of cleanerfish by their fish clients. Clients on an Indonesian reef flat readily approached wooden models of the bluestreak cleaner wrasse, Labroides dimidiatus, which varied in size and stripe characteristics. The composition of the clientele of models was not significantly different from that of natural cleaning stations, suggesting that fish visiting the models were seeking to be cleaned. Normal-sized models of cleaner wrasses attracted significantly more clients, which showed more intense interest and stayed with the models for significantly longer, than super-sized models. For normal-sized models, the number of clients increased as the length of the cleaner's lateral stripe increased (from 0, to 44, 67 and 100% of body length). However, there was no effect of stripe length on client numbers for super-sized models. Client interest also did not vary with stripe length for models of either size. Small body size and the presence of a lateral stripe therefore appear to be long-distance signals that their bearer may be a cleaner, but after initial attraction, client interest is maintained by other cues. Alternative short-distance signals may include colour, other visual signals such as cleanerfish dances, or physical contact between cleaner and client.     

Sample size calculation for phylogenetic case linkage

Sample size calculations are an essential component of the design and evaluation of scientific studies. However, there is a lack of clear guidance for determining the sample size needed for phylogenetic studies, which are becoming an essential part of studying pathogen transmission. We introduce a statistical framework for determining the number of true infector-infectee transmission pairs identified by a phylogenetic study, given the size and population coverage of that study. We then show how characteristics of the criteria used to determine linkage and aspects of the study design can influence our ability to correctly identify transmission links, in sometimes counterintuitive ways. We test the overall approach using outbreak simulations and provide guidance for calculating the sensitivity and specificity of the linkage criteria, the key inputs to our approach. The framework is freely available as the R package phylosamp, and is broadly applicable to designing and evaluating a wide array of pathogen phylogenetic studies.     

Clasp-mediated microtubule bundling regulates persistent motility and contact repulsion in Drosophila macrophages in vivo

Drosophila melanogaster macrophages are highly migratory cells that lend themselves beautifully to high resolution in vivo imaging experiments. By expressing fluorescent probes to reveal actin and microtubules, we can observe the dynamic interplay of these two cytoskeletal networks as macrophages migrate and interact with one another within a living organism. We show that before an episode of persistent motility, whether responding to developmental guidance or wound cues, macrophages assemble a polarized array of microtubules that bundle into a compass-like arm that appears to anticipate the direction of migration. Whenever cells collide with one another, their microtubule arms transiently align just before cell–cell repulsion, and we show that forcing depolymerization of microtubules by expression of Spastin leads to their defective polarity and failure to contact inhibit from one another. The same is true in orbit/clasp mutants, indicating a pivotal role for this microtubule-binding protein in the assembly and/or functioning of the microtubule arm during polarized migration and contact repulsion.     

Assessment of cell alignment by fibronectin multi-fibre cables capable of large scale production

The repair of damaged human tissue will be enhanced greatly by a capacity to organise the arrangement of the cells. We have demonstrated an approach to quantifying the capacity of fibronectin multi-fibre cables to align human cells. It is based on staining and subsequent image analysis of cells in the neighbourhood of the cables, and the application of an Orientation Index (S), to give a quantitative measure of alignment of the cells. Alignment of human dermal fibroblasts, parallel to the fibronectin cable axis, was observed by light microscopy after 3 days in culture. At 7 days, alignment was determined by image analysis after hematoxylin and eosin staining. This showed that the lateral extent of cell alignment i.e. the number of cell layers lined up beside a cable, termed the cell docking band width, was independent of cable diameter. Scanning electron microscopy confirmed cell alignment and the deposition of fine collagen fibrils by the cells. Orientation of cells to the cable, as measured by an orientation index (S), was S = 0.97 - 0.05 in cell docking bands denoting almost perfect alignment. Cell seeding levels ranging from 14 K cells/cm² up to 54 K cells/cm² resulted in cell docking band widths of aligned cells from 335 wm to 890 wm respectively, increasing as a function of cell seeding levels. This method represents a quantitative measure of quality for potential contact guidance materials. The guidance system is capable of large-scale manufacture.     

Four-dimensional imaging with virtual reality to quantitatively explore jigsaw puzzle-like morphogenesis of Arabidopsis cotyledon pavement cells

In most dicotyledonous plants, leaf pavement cells exhibit complex jigsaw puzzle-like cell morphogenesis during leaf expansion. Although detailed molecular biological information and mathematical modeling of this jigsaw puzzle-like cell morphogenesis are now available, a full understanding of this process remains elusive. Recent reports have highlighted the importance of three-dimensional (3D) structures (i.e., anticlinal and periclinal cell wall) in understanding the mechanical models that describe this morphogenetic process. We believe that it is important to acquire 3D shapes of pavement cells over time, i.e., acquire and analyze four-dimensional (4D) information when studying the relationship between mechanical modeling and simulations and the actual cell shape. In this report, we have developed a framework to capture and analyze 4D morphological information of Arabidopsis thaliana cotyledon pavement cells by using both direct water immersion observations and computational image analyses, including segmentation, surface modeling, virtual reality and morphometry. The 4D cell models allowed us to perform time-lapse 3D morphometrical analysis, providing detailed quantitative information about changes in cell growth rate and shape, with cellular complexity observed to increase during cell growth. The framework should enable analysis of various phenotypes (e.g., mutants) in greater detail, especially in the 3D deformation of the cotyledon surface, and evaluation of theoretical models that describe pavement cell morphogenesis using computational simulations. Additionally, our accurate and high-throughput acquisition of growing cell structures should be suitable for use in generating in silico model cell structures.     

The CC-NB-LRR-Type Rdg2a Resistance Gene Confers Immunity to the Seed-Borne Barley Leaf Stripe Pathogen in the Absence of Hypersensitive Cell Death

Background

Leaf stripe disease on barley (Hordeum vulgare) is caused by the seed-transmitted hemi-biotrophic fungus Pyrenophora graminea. Race-specific resistance to leaf stripe is controlled by two known Rdg (Resistance to Drechslera graminea) genes: the H. spontaneum-derived Rdg1a and Rdg2a, identified in H. vulgare. The aim of the present work was to isolate the Rdg2a leaf stripe resistance gene, to characterize the Rdg2a locus organization and evolution and to elucidate the histological bases of Rdg2a-based leaf stripe resistance.

Principal Findings

We describe here the positional cloning and functional characterization of the leaf stripe resistance gene Rdg2a. At the Rdg2a locus, three sequence-related coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NB-LRR) encoding genes were identified. Sequence comparisons suggested that paralogs of this resistance locus evolved through recent gene duplication, and were subjected to frequent sequence exchange. Transformation of the leaf stripe susceptible cv. Golden Promise with two Rdg2a-candidates under the control of their native 5′ regulatory sequences identified a member of the CC-NB-LRR gene family that conferred resistance against the Dg2 leaf stripe isolate, against which the Rdg2a-gene is effective. Histological analysis demonstrated that Rdg2a-mediated leaf stripe resistance involves autofluorescing cells and prevents pathogen colonization in the embryos without any detectable hypersensitive cell death response, supporting a cell wall reinforcement-based resistance mechanism.

Conclusions

This work reports about the cloning of a resistance gene effective against a seed borne disease. We observed that Rdg2a was subjected to diversifying selection which is consistent with a model in which the R gene co-evolves with a pathogen effector(s) gene. We propose that inducible responses giving rise to physical and chemical barriers to infection in the cell walls and intercellular spaces of the barley embryo tissues represent mechanisms by which the CC-NB-LRR-encoding Rdg2a gene mediates resistance to leaf stripe in the absence of hypersensitive cell death.     

Intensity Correlation-Based Calibration of FRET

Dual-laser flow cytometric resonance energy transfer (FCET) is a statistically efficient and accurate way of determining proximity relationships for molecules of cells even under living conditions. In the framework of this algorithm, absolute fluorescence resonance energy transfer (FRET) efficiency is determined by the simultaneous measurement of donor-quenching and sensitized emission. A crucial point is the determination of the scaling factor α responsible for balancing the different sensitivities of the donor and acceptor signal channels. The determination of α is not simple, requiring preparation of special samples that are generally different from a double-labeled FRET sample, or by the use of sophisticated statistical estimation (least-squares) procedures. We present an alternative, free-from-spectral-constants approach for the determination of α and the absolute FRET efficiency, by an extension of the presented framework of the FCET algorithm with an analysis of the second moments (variances and covariances) of the detected intensity distributions. A quadratic equation for α is formulated with the intensity fluctuations, which is proved sufficiently robust to give accurate α-values on a cell-by-cell basis in a wide system of conditions using the same double-labeled sample from which the FRET efficiency itself is determined. This seemingly new approach is illustrated by FRET measurements between epitopes of the MHCI receptor on the cell surface of two cell lines, FT and LS174T. The figures show that whereas the common way of α determination fails at large dye-per-protein labeling ratios of mAbs, this presented-as-new approach has sufficient ability to give accurate results. Although introduced in a flow cytometer, the new approach can also be straightforwardly used with fluorescence microscopes.     

A deterministic model of admixture and genetic introgression: the case of Neanderthal and Cro-Magnon

There is an ongoing debate in the field of human evolution about the possible contribution of Neanderthals to the modern human gene pool. To study how the Neanderthal private alleles may have spread over the genes of Homo sapiens, we propose a deterministic model based on recursive equations and ordinary differential equations. If the Neanderthal population was large compared to the Homo sapiens population at the beginning of the contact period, we show that genetic introgression should have been fast and complete meaning that most of the Neanderthal private alleles should be found in the modern human gene pool in case of ancient admixture. In order to test/reject ancient admixture from genome-wide data, we incorporate the model of genetic introgression into a statistical hypothesis-testing framework. We show that the power to reject ancient admixture increases as the ratio, at the time of putative admixture, of the population size of Homo sapiens over that of Neanderthal decreases. We find that the power to reject ancient admixture might be particularly low if the population size of Homo sapiens was comparable to the Neanderthal population size.